Optimizing stellarators for large flows

Plasma flow is damped in stellarators because they are not intrinsically ambipolar, unlike tokamaks, in which the flux-surface averaged radial electric current vanishes for any value of the radial electric field. Only quasisymmetric stellarators are intrinsically ambipolar, but exact quasisymmetry is impossible to achieve in non-axisymmetric toroidal configurations. By calculating the violation of intrinsic ambipolarity due to deviations from quasisymmetry, one can derive criteria to assess when a stellarator can be considered quasisymmetric in practice, i.e. when the flow damping is weak enough. Let us denote by $\alpha$ a small parameter that controls the size of a perturbation to an exactly quasisymmetric magnetic field. Recently, it has been shown that if the gradient of the perturbation is sufficiently small, the flux-surface averaged radial electric current scales as $\alpha^2$ for any value of the collisionality. It was also argued that when the gradient of the perturbation is large, the quadratic scaling is replaced by a more unfavorable one. In this paper, perturbations with large gradients are rigorously treated. In particular, it is proven that for low collisionality a perturbation with large gradient yields, at best, an $O(|\alpha|)$ deviation from quasisymmetry. Heuristic estimations in the literature incorrectly predicted an $O(|\alpha|^{3/2})$ deviation.Comment: 24 pages, 2 figures. To appear in Plasma Physics and Controlled Fusio

Flow damping in stellarators close to quasisymmetry

Quasisymmetric stellarators are a type of optimized stellarators for which flows are undamped to lowest order in an expansion in the normalized Larmor radius. However, perfect quasisymmetry is impossible. Since large flows may be desirable as a means to reduce turbulent transport, it is important to know when a stellarator can be considered to be sufficiently close to quasisymmetry. The answer to this question depends strongly on the size of the spatial gradients of the deviation from quasisymmetry and on the collisionality regime. Recently, formal criteria for closeness to quasisymmetry have been derived in a variety of situations. In particular, the case of deviations with large gradients was solved in the $1/\nu$ regime. Denoting by $\alpha$ a parameter that gives the size of the deviation from quasisymmetry, it was proven that particle fluxes do not scale with $\alpha^{3/2}$, as typically claimed, but with $\alpha$. It was also shown that ripple wells are not necessarily the main cause of transport. This paper reviews those works and presents a new result in another collisionality regime, in which particles trapped in ripple wells are collisional and the rest are collisionless.Comment: 14 pages, 2 figures. To appear in Plasma Physics and Controlled Fusio

The effect of tangential drifts on neoclassical transport in stellarators close to omnigeneity

In general, the orbit-averaged radial magnetic drift of trapped particles in stellarators is non-zero due to the three-dimensional nature of the magnetic field. Stellarators in which the orbit-averaged radial magnetic drift vanishes are called omnigeneous, and they exhibit neoclassical transport levels comparable to those of axisymmetric tokamaks. However, the effect of deviations from omnigeneity cannot be neglected in practice. For sufficiently low collision frequencies (below the values that define the $1/\nu$ regime), the components of the drifts tangential to the flux surface become relevant. This article focuses on the study of such collisionality regimes in stellarators close to omnigeneity when the gradient of the non-omnigeneous perturbation is small. First, it is proven that closeness to omnigeneity is required to preserve radial locality in the drift-kinetic equation for collisionalities below the $1/\nu$ regime. Then, it is shown that neoclassical transport is determined by two layers in phase space. One of the layers corresponds to the $\sqrt{\nu}$ regime and the other to the superbanana-plateau regime. The importance of the superbanana-plateau layer for the calculation of the tangential electric field is emphasized, as well as the relevance of the latter for neoclassical transport in the collisionality regimes considered in this paper. In particular, the tangential electric field is essential for the emergence of a new subregime of superbanana-plateau transport when the radial electric field is small. A formula for the ion energy flux that includes the $\sqrt{\nu}$ regime and the superbanana-plateau regime is given. The energy flux scales with the square of the size of the deviation from omnigeneity. Finally, it is explained why below a certain collisionality value the formulation presented in this article ceases to be valid.Comment: 36 pages. Version to be published in Plasma Physics and Controlled Fusio

When omnigeneity fails

A generic non-symmetric magnetic field does not confine magnetized charged particles for long times due to secular magnetic drifts. Stellarator magnetic fields should be omnigeneous (that is, designed such that the secular drifts vanish), but perfect omnigeneity is technically impossible. There always are small deviations from omnigeneity that necessarily have large gradients. The amplification of the energy flux caused by a deviation of size $\epsilon$ is calculated and it is shown that the scaling with $\epsilon$ of the amplification factor can be as large as linear. In opposition to common wisdom, most of the transport is not due to particles trapped in ripple wells, but to the perturbed motion of particles trapped in the omnigeneous magnetic wells around their bounce points.Comment: 6 pages, 2 figure

Electrostatic potential variations on stellarator magnetic surfaces in low collisionality regimes

The component of the neoclassical electrostatic potential that is non-constant on the magnetic surface, that we denote by $\tilde\varphi$, can affect radial transport of highly charged impurities, and this has motivated its inclusion in some modern neoclassical codes. The number of neoclassical simulations in which $\tilde\varphi$ is calculated is still scarce, partly because they are usually demanding in terms of computational resources, especially at low collisionality. In this paper the size, the scaling with collisionality and with aspect ratio, and the structure of $\tilde\varphi$ on the magnetic surface are analytically derived in the $1/\nu$, $\sqrt{\nu}$ and superbanana-plateau regimes of stellarators close to omnigeneity; i. e. stellarators that have been optimized for neoclassical transport. It is found that the largest $\tilde\varphi$ that the neoclassical equations admit scales linearly with the inverse aspect ratio and with the size of the deviation from omnigeneity. Using a model for a perturbed omnigeneous configuration, the analytical results are verified and illustrated with calculations by the code KNOSOS. The techniques, results and numerical tools employed in this paper can be applied to neoclassical transport problems in tokamaks with broken axisymmetry.Comment: 30 pages, 12 figures, 1 table. Published versio

Radiative decay Z_H-> \gamma A_H in the little Higgs model with T-parity

In the little Higgs model with T-parity (LHTM), the only tree-level kinematically allowed two-body decay of the Z_H boson is Z_H-> A_H H and thus one-loop induced two-body decays may have a significant rate. We study the Z_H-> \gamma A_H decay, which is induced at the one-loop level by a fermion triangle and is interesting as it depends on the mechanism of anomaly cancellation of the model. All the relevant two- and three-body decays of the Z_H gauge boson arising at the tree-level are also calculated. We consider a small region of the parameter space where the scale of the symmetry breaking f is still allowed to be as low as 500 GeV by electroweak precision data. We first analyze the scenario of a Higgs boson with a mass of 120 GeV. We found that the Z_H->\gamma A_H branching ratio can be of the order of a tree-level three-body decay and may be at the reach of detection at the LHC for f close to 500 GeV, but it may be difficult to detect for f=1 TeV. There is also an scenario where the Higgs boson has an intermediate mass such that the Z_H-> A_H H decay is closed, the Z_H-> \gamma A_H gets considerably enhanced and the chances of detection get a large boost.Comment: 19 pages, 9 figures, 2 table