Search

In omnigeneous magnetic fields, charged particles are perfectly confined in the absence of collisions and turbulence. For this reason, the magnetic configuration is optimized to be close to omnigenity in any candidate for a stellarator fusion reactor. However, approaching omnigenity imposes severe constraints on the spatial variation of the magnetic field. In particular, the topology of the contours of constant magnetic-field-strength on each magnetic surface must be such that there are no particles transitioning between different types of wells. This, in turn, usually leads to complicated plasma shapes and coils. This Letter presents a new family of optimized fields that display tokamak-like collisional energy transport while having transitioning particles. This result radically broadens the space of accessible reactor-relevant configurations.

CIEMAT-QI4 is a quasi-isodynamic stellarator configuration that simultaneously features very good fast-ion confinement in a broad range of $\beta$ values, low neoclassical transport and bootstrap current, and ideal magnetohydrodynamic stability up to $\beta=5\%$. In this paper it is shown that CIEMAT-QI4 also exhibits reduced turbulent transport. This is demonstrated through nonlinear electrostatic simulations with the gyrokinetic code stella, including kinetic ions and electrons. The relation between reduced turbulent transport and the fact that CIEMAT-QI4 very approximately satisfies the so-called maximum-J property is discussed., Comment: 12 pages, 9 figures

Efficient control of turbulent heat transport is crucial for magnetic confinement fusion reactors. This work discusses the complex interplay between density gradients and micro-instabilities, shedding light on their impact on turbulent heat transport in different fusion devices. In particular, the influence of density gradients on turbulent heat transport is investigated through an extensive inter-machine study, including various stellarators such as W7-X, LHD, TJ-II and NCSX, along with the Asdex-Upgrade tokamak and the tokamak geometry of the Cyclone Base Case (CBC). Linear and nonlinear simulations are performed employing the $\delta$f-gyrokinetic code stella across a wide range of parameters to explore the effects of density gradients, temperature gradients, and kinetic electrons. A strong reduction in ion heat flux with increasing density gradients is found in NCSX and W7-X due to the stabilization of temperature-gradient-driven modes without significantly destabilizing density-gradient-driven modes. In contrast, the tokamaks exhibit an increase in ion heat flux with density gradients. Notably, the behavior of ion heat fluxes in stellarators does not align with that of linear growth rates. Additionally, this study provides physical insights into the micro-instabilities, emphasizing the dominance of trapped-electron-modes in CBC, AUG, TJ-II, LHD and NCSX, while both the trapped-electron-mode and the passing-particle-driven universal instability contribute significantly in W7-X., Comment: 26 pages, 30 figures

MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES/

The first experimental campaigns of Wendelstein 7-X (W7-X) have shown that turbulence plays a decisive role in the performance of neoclassically optimized stellarators. This stresses the importance of understanding microturbulence from the theoretical and experimental points of view. To this end, this paper addresses a comprehensive characterization of the turbulent fluctuations by means of nonlinear gyrokinetic simulations performed with the code stella in two W7-X scenarios. In the first part of the paper, the amplitude of the density fluctuations is calculated and compared with measurements obtained by Doppler reflectometry (DR) in the OP1 experimental campaigns. It is found that the trend of the fluctuations along the radius is explained by the access of the DR system to different regions of the turbulence wavenumber spectrum. In the second part of the article, frequency spectra of the density fluctuations and the zonal component of the turbulent flow are numerically characterized for comparisons against future experimental analyses. Both quantities feature broad frequency spectra with dominant frequencies of O(1)-O(10) kHz., Comment: 27 pages, 16 figures

Stellarator magnetic configurations need to be optimized in order to meet all the required properties of a fusion reactor. In this work, it is shown that a flat-mirror quasi-isodynamic configuration (i.e. a quasi-isodynamic configuration with sufficiently small radial variation of the mirror term) can achieve small radial transport of energy and good confinement of bulk and fast ions even if it is not very close to perfect omnigeneity, and for a wide range of plasma scenarios, including low $\beta$ and small radial electric field. This opens the door to constructing better stellarator reactors. On the one hand, they would be easier to design, as they would be robust against error fields. On the other hand, they would be easier to operate since, both during startup and steady-state operation, they would require less auxiliary power, and the damage to plasma-facing components caused by fast ion losses would be reduced to acceptable levels.

A systematic study of the impact of impurities on the turbulent heat fluxes is presented for the stellarator Wendelstein 7-X (W7-X) and, for comparison, the Large Helical Device and ITER. By means of nonlinear multispecies gyrokinetic simulations, it is shown that impurities, depending on the sign of their density gradient, can significantly enhance or reduce turbulent ion heat losses. For the relevant scenario of turbulence reduction, an optimal impurity concentration that minimizes the ion heat diffusivity emerges as a universal feature. This result demonstrates the potential of impurities for controlling turbulence and accessing enhanced confinement regimes in fusion plasmas and, in particular, in W7-X., Comment: 8 pages, 8 figures

A new quasi-isodynamic stellarator configuration optimized for the confinement of energetic ions at low plasma $\beta$ is obtained. The numerical optimization is carried out using the STELLOPT suite of codes. New proxies to measure closeness to quasi-isodynamicity and quality of fast ion confinement have been included. The new configuration has poloidally closed contours of magnetic field strength, low magnetic shear and a rotational transform profile allowing an island divertor. It shows ideal and ballooning magnetohydrodynamic stability up to $\beta = 5%$, reduced effective ripple, with $\epsilon_{eff} < 0.5%$ in the plasma core. Even at low $\beta$, the configuration approximately satisfies the maximum-$J$ property, and the confinement of fast ions is good at $\beta \sim 1.5%$ and becomes excellent at reactor values, $\beta \sim 4%$. An evaluation of the $D_{31}$ neoclassical mono-energetic coefficient supports the expectation of a reduced bootstrap current for plasmas confined in quasi-isodynamic configurations. A set of filamentary coils that preserve the good confinement of fast ions in the core is presented., Comment: 21 pages, 16 figures

In reactor-relevant plasmas, neoclassical transport drives an outward particle flux in the core of large stellarators and predicts strongly hollow density profiles. However, this theoretical prediction is contradicted by experiments. In particular, in Wendelstein 7-X, the first large optimized stellarator, flat or weakly peaked density profiles are generally measured, indicating that neoclassical theory is not sufficient and that an inward contribution to the particle flux is missing in the core. In this Research Letter, it is shown that the turbulent contribution to the particle flux can explain the difference between experimental measurements and neoclassical predictions. The results of this Research Letter also prove that theoretical and numerical tools are approaching the level of maturity needed for the prediction of equilibrium density profiles in stellarator plasmas, which is a fundamental requirement for the design of operation scenarios of present devices and future reactors., Comment: 8 pages, 6 figures, updated the title and incorporated the appendix into the main letter

A recent characterization of core turbulence carried out with a Doppler reflectometer in the optimized stellarator Wendelstein 7-X (W7-X) found that discharges achieving high ion temperatures at the core featured an ITG-like suppression of density fluctuations driven by a reduction of the gradient ratio $\eta_i = L_n/L_{T_i}$ [D. Carralero et al., Nucl. Fusion, 2021]. In order to confirm the role of ITG turbulence in this process, we set out to establish experimentally the relation between core density fluctuations, turbulent heat flux and global confinement. With this aim, we consider the scenarios found in the previous work and carry out power balance analysis for a number of representative ones, including some featuring high ion temperature. As well, we evaluate the global energy confinement time and discuss it in the context of the ISS04 inter-stellarator scaling. We find that, when turbulence is suppressed as a result of a reduction of $\eta_i$, there is a reduction of ion turbulent transport, and global performance is improved as a result. This is consistent with ITG turbulence limiting the ion temperature at the core of W7-X. In contrast, when turbulence is reduced following a decrease in collisionality, no changes are observed in transport or confinement. This could be explained by ITG modes being combined with TEM turbulence when the later is destabilized at low collisionalities., Comment: Submitted to PPCF, 47th EPS special issue

The ongoing development of electromagnets based on High Temperature Superconductors has led to the conceptual exploration of high-magnetic-field fusion reactors of the tokamak type, operating at on-axis fields above 10 T. In this work we explore the consequences of the potential future availability of high-field three-dimensional electromagnets on the physics design point of a stellarator reactor. We find that, when an increase in the magnetic field strength $B$ is used to maximally reduce the device linear size $R\sim B^{-4/3}$ (with otherwise fixed magnetic geometry), the physics design point is largely independent of the chosen field strength/device size. A similar degree of optimization is to be imposed on the magnetohydrodynamic, transport and fast ion confinement properties of the magnetic configuration of that family of reactor design points. Additionally, we show that the family shares an invariant operation map of fusion power output as a function of the auxiliary power and relative density variation. The effects of magnetic field over-engineering and the $R(B)$ scaling of design points with constant neutron wall loading are also inspected. In this study we use geometric parameters characteristic of the helias reactor, but most results apply to other stellarator configurations

The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and for the preparation of the next campaigns. At the same time, new stellarator gyrokinetic codes have just been or are being developed. It is therefore desirable to have a sufficiently complete, documented and verified set of gyrokinetic simulations in W7-X geometry against which new codes or upgrades of existing codes can be tested and benchmarked. This paper attemps to provide such a set of simulations in the form of a comprehensive benchmark between the recently developed code stella and the well-established code GENE. The benchmark consists of electrostatic gyrokinetic simulations in W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode (TEM)} stability analyses, computation of linear zonal flow responses and calculation of ITG-driven heat fluxes.

A good understanding of the confinement of energetic ions in non-axisymmetric magnetic fields is key for the design of reactors based on the stellarator concept. In this work, we develop a model that, based on the radially-local bounce-averaged drift-kinetic equation, classifies orbits and succeeds in predicting configuration-dependent aspects of the prompt losses of energetic ions in stellarators. Such a model could in turn be employed in the optimization stage of the design of new devices.

A large diffusive turbulent contribution to the radial impurity transport in Wendelstein 7-X (W7-X) plasmas has been experimentally inferred during the first campaigns and numerically confirmed by means of gyrokinetic simulations with the code stella. In general, the absence of strong impurity accumulation during the initial W7-X campaigns is attributed to this diffusive term. In the present work the diffusive contribution is also calculated in other stellarator plasmas. In particular, the diffusion (D) and convection (V) coefficients of carbon and iron impurities produced by ion-temperature-gradient (ITG) turbulence are obtained for W7-X, LHD, TJ-II and NCSX. The results show that, although the size of D and V can differ across the four devices, inward convection is found for all of them. For W7-X, TJ-II and NCSX the two coefficients are comparable and the turbulent peaking factor is surprisingly similar. In LHD, appreciably weaker diffusive and convective impurity transport and significantly larger turbulent peaking factor are predicted. All this suggests that ITG turbulence, although not strongly, would lead to negative impurity density gradients in stellarators. Then, considering mixed ITG/Trapped Electron Mode (TEM) turbulence for the specific case of W7-X, it has been quantitatively assessed to what degree pellet fueled reduced turbulence scenarios feature reduced turbulent transport of impurities as well. The results for trace iron impurities show that, although their turbulent transport is not entirely suppressed, a significant reduction of V and a stronger decrease of D are found. Although the diffusion is still above neoclassical levels, the neoclassical convection would gain under such conditions a greater specific weight on the dynamics of impurities in comparison with standard ECRH scenarios., Comment: 16 pages, 8 figures

In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube, full-flux-surface, and radially global domains. Two problems are studied: the linear relaxation of zonal flows and the linear stability of ion temperature gradient (ITG) modes. Simulations are carried out with the codes EUTERPE, GENE, GENE-3D, and stella in magnetic configurations of LHD and W7-X using adiabatic electrons. The zonal flow relaxation properties obtained in different flux tubes are found to differ with each other and with the radially global result, except for sufficiently long flux tubes, in general. The flux tube length required for convergence is configuration-dependent. Similarly, for ITG instabilities, different flux tubes provide different results, but the discrepancy between them diminishes with increasing flux tube length. Full-flux-surface and flux tube simulations show good agreement in the calculation of the growth rate and frequency of the most unstable modes in LHD, while for W7-X differences in the growth rates are found between the flux tube and the full-flux-surface domains. Radially global simulations provide results close to the full-flux-surface ones. The radial scale of unstable ITG modes is studied in global and flux tube simulations finding that in W7-X, the radial scale of the most unstable modes depends on the binormal wavenumber, while in LHD no clear dependency is found., Comment: submitted to Nuclear Fusion

In this work, a new version of KNOSOS is presented. KNOSOS is a low-collisionality radially-local, bounce-averaged neoclassical code that is extremely fast, and at the same time, includes physical effects often neglected by more standard codes: the component of the magnetic drift that is tangent to the flux-surface and the variation of the electrostatic potential on the flux-surface. An earlier version of the code could only describe configurations that were sufficiently optimized with respect to neoclassical transport. KNOSOS can now be applied to any large aspect ratio stellarator, and its performance is demonstrated by means of detailed simulations in the configuration space of Wendelstein 7-X., Comment: IAEA paper, submitted to Nuclear Fusion

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, turbulent transport poses a serious limitation to the global performance of the machine. Several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of $T_{i} \simeq 1.7$ keV, except in the few scenarios in which turbulence can be suppressed. In order to understand turbulent mechanisms at play, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. In this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ($k_\perp\rho_i\simeq 1$) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients $a/L_{ne}$ and $a/L_{Ti}$. Furthermore, we discuss the importance of other potentially relevant parameters such as $T_e/T_i$, $E_r$ or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio $\eta_i=L_{ne}/L_{Ti}$, as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which $T_{i,core} > 1.7$ keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X., Comment: Manuscript submitted to Nuclear Fusion in May 2021

The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate $k_x$, but the oscillation and damping time scales are longer in local representations, particularly for small $k_x$ approaching the system size., Comment: The following article has been accepted by Physics of Plasmas. After it is published, it will be found at https://aip.scitation.org/journal/php. 33 pages, 18 figures

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ITG, TEM and ETG modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. It is also confirmed that thermo-diffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively., Comment: 19 pages, 10 figures, 2 tables

KNOSOS (KiNetic Orbit-averaging SOlver for Stellarators) is a freely available, open-source code (\href{https://github.com/joseluisvelasco/KNOSOS}{https://github.com/joseluisvelasco/KNOSOS}) that calculates neoclassical transport in low-collisionality plasmas of three-dimensional magnetic confinement devices by solving the radially local drift-kinetic and quasineutrality equations. The main feature of KNOSOS is that it relies on orbit-averaging to solve the drift-kinetic equation very fast. KNOSOS treats rigorously the effect of the component of the magnetic drift that is tangent to magnetic surfaces, and of the component of the electrostatic potential that varies on the flux surface, {\varphi}_1. Furthermore, the equation solved is linear in {\varphi}_1, which permits an efficient solution of the quasineutrality equation. As long as the radially local approach is valid, KNOSOS can be applied to the calculation of neoclassical transport in stellarators (helias, heliotrons, heliacs, etc.) and tokamaks with broken axisymmetry. In this paper, we show several calculations for the stellarators W7-X, LHD, NCSX and TJ-II that provide benchmark with standard local codes and demonstrate the advantages of this approach.

In the present work we report recent radial electric field measurements carried out with the Doppler reflectometry system in the TJ-II stellarator. The study focuses on the fact that, under some conditions, the radial electric field measured at different points over the same flux surface shows significantly different values. A numerical analysis is carried out considering the contribution arising from the radial dependence of $\Phi_1$ as a possible correction term to the total radial electric field. Here $\Phi_1$ is the neoclassical electrostatic potential variation over the surface. The comparison shows good agreement in some aspects, like the conditions under which this correction is large (electron-root conditions) or negligible (ion-root conditions). But it disagrees in others like the sign of the correction. The results are discussed together with the underlying reasons of this partial disagreement. In addition, motivated by the recent installation of the dual Doppler reflectometry system in Wendelstein 7-X (W7-X), $\Phi_1$ estimations for W7-X are revisited considering Core-Electron-Root-Plasma (CERC) plasmas from its first experimental campaign. The simulations show larger values of $\Phi_1$ under electron-root conditions than under ion root ones. The contribution from the kinetic electron response is shown to become important at some radii. All this results in a potential variation size noticeably larger than estimated in our previous work in W7-X \cite{Regana_nf_57_056004_2017} for other plasma parameters and another configuration., Comment: 23 pages, 10 figures

The low frequency oscillatory relaxation of zonal potential perturbations is studied numerically in the TJ-II stellarator (where it was experimentally detected for the first time). It is studied in full global gyrokinetic simulations of multi-species plasmas. The oscillation frequency obtained is compared with predictions based on single-species simulations using simplified analytical relations. It is shown that the frequency of this oscillation for a multi-species plasma can be accurately obtained from single-species calculations using extrapolation formulas. The damping of the oscillation and the influence of the different inter-species collisions is studied in detail. It is concluded that taking into account multiple kinetic ions and electrons with impurity concentrations realistic for TJ-II plasmas allows to account for the values of frequency and damping rate in zonal flows relaxations observed experimentally., Comment: 11 figures, 22 pages

Low-collisionality stellarator plasmas usually display a large negative radial electric field that has been expected to cause accumulation of impurities due to their high charge number. In this paper, two combined effects that can potentially modify this scenario are discussed. First, it is shown that, in low collisionality plasmas, the kinetic contribution of the electrons to the radial electric field can make it negative but small, bringing the plasma close to impurity temperature screening (i.e., to a situation in which the ion temperature gradient is the main drive of impurity transport and causes outward flux); in plasmas of very low collisionality, such as those of the Large Helical Device displaying impurity hole, screening may actually occur. Second, the component of the electric field that is tangent to the flux surface (in other words, the variation of the electrostatic potential on the flux surface), although smaller than the radial component, has recently been suggested to be an additional relevant drive for radial impurity transport. Here, it is explained that, especially when the radial electric field is small, the tangential magnetic drift has to be kept in order to correctly compute the tangential electric field, that can be larger than previously expected. This can have a strong impact on impurity transport, as we illustrate by means of simulations using the newly-developed code KNOSOS (KiNetic Orbit-averaging-SOlver for Stellarators).

Measurements of the relaxation of a zonal electrostatic potential perturbation in a non-axisymmetric magnetically confined plasma are presented. A sudden perturbation of the plasma equilibrium is induced by the injection of a cryogenic hydrogen pellet in the TJ-II stellarator, which is observed to be followed by a damped oscillation in the electrostatic potential. The waveform of the relaxation is consistent with theoretical calculations of zonal potential relaxation in a non-axisymmetric magnetic geometry. The turbulent transport properties of a magnetic confinement configuration are expected to depend on the features of the collisionless damping of zonal flows, of which the present letter is the first direct observation.

Achieving impurity and helium ash control is a crucial issue in the path towards fusion-grade magnetic confinement devices, and this is particularly the case of helical reactors, whose low-collisionality ion-root operation scenarios usually display a negative radial electric field which is expected to cause inwards impurity pinch. In these work we discuss, based on experimental measurements and standard predictions of neoclassical theory, how plasmas of very low ion collisionality, similar to those observed in the impurity hole of the Large Helical Device, can be an exception to this general rule, and how a negative radial electric field can coexist with an outward impurity flux. This interpretation is supported by comparison with documented discharges available in the International Stellarator-Heliotron Profile Database, and it can be extrapolated to show that achievement of high ion temperature in the core of helical devices is not fundamentally incompatible with low core impurity content.

We review in a tutorial fashion some of the causes of impurity density variations along field lines and radial impurity transport in the moment approach framework. An explicit and compact form of the parallel inertia force valid for arbitrary toroidal geometry and magnetic coordinates is derived and shown to be non-negligible for typical TJ-II plasma conditions. In the second part of the article, we apply the fluid model including main ion-impurity friction and inertia to observations of asymmetric emissivity patterns in neutral beam heated plasmas of the TJ-II stellarator. The model is able to explain qualitatively several features of the radiation asymmetry, both in stationary and transient conditions, based on the calculated in-surface variations of the impurity density.

We study radial particle transport in stellarator plasmas using cryogenic pellet injection. By means of perturbative experiments, we estimate the experimental particle flux and compare it with neoclassical simulations. Experimental evidence is obtained of the fact that core depletion in helical devices can be slowed-down even by pellets that do not reach the core region. This phenomenon is well captured by neoclassical predictions with DKES and FORTEC-3D.

In the linear collisionless limit, a zonal potential perturbation in a toroidal plasma relaxes, in general, to a non-zero residual value. Expressions for the residual value in tokamak and stellarator geometries, and for arbitrary wavelengths, are derived. These expressions involve averages over the lowest order particle trajectories, that typically cannot be evaluated analytically. In this work, an efficient numerical method for the evaluation of such expressions is reported. It is shown that this method is faster than direct gyrokinetic simulations performed with the GENE and EUTERPE codes. Calculations of the residual value in stellarators are provided for much shorter wavelengths than previously available in the literature. Electrons must be treated kinetically in stellarators because, unlike in tokamaks, kinetic electrons modify the residual value even at long wavelengths. This effect, that had already been predicted theoretically, is confirmed by gyrokinetic simulations., Comment: 31 pages, 15 figures

The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {\bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {\bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated explicitly., Comment: 6 pages

First observations of electrostatic potential variations within the flux surfaces of a toroidal magnetic confinement device are presented. Measurements are taken in the TJ-II stellarator with two distant Langmuir probe arrays. The edge floating potentials display differences of several tens of Volts in electron-root wave-heated plasmas. The differences are reduced for higher densities and lower electron temperatures after the ion-root electric field forms at the plasma edge. Neoclassical Monte Carlo simulations estimate the correct order of magnitude for the overall variation in potential and predict the trend observed with the radial electric field. However, for the specific location of the probes, the simulations give differences smaller than those observed experimentally.

The drift kinetic equation is solved for low density TJ-II plasmas employing slowly varying, time-dependent profiles. This allows to simulate density ramp-up experiments and to describe from first principles the formation and physics of the radial electric field shear, which is associated to the transition from electron to ion root. We show that the range of frequencies of plasma potential fluctuations in which zonal flows are experimentally observed is neoclassically undamped in a neighborhood of the transition. This makes the electron root regime of stellarators, close to the transition to ion root, a propitious regime for the study of zonal-flow evolution. We present simulations of collisionless relaxation of zonal flows, in the sense of the Rosenbluth and Hinton test, that show an oscillatory behaviour in qualitative agreement with the experiment close to the transition.

The drift kinetic equation is solved for low density TJ-II plasmas employing slowly varying, time-dependent profiles. This allows to simulate density ramp-up experiments and describe from first principles the formation and physics of the radial electric field shear layer. The main features of the transition are perfectly captured by the calculation, and good quantitative agreement is also found. The results presented here, that should be valid for other non-quasisymmetric stellarators, provide a fundamental explanation for a wealth of experimental observations connected to the shear layer emergence in TJ-II. The key quantity is the neoclassical viscosity, which is shown to go smoothly to zero when the critical density is approached from below. This makes it possible for turbulence-related phenomena, and particularly zonal flows, to arise in the neighborhood of the transition., Comment: 5 pages, 5 figures, accepted for publication in Physical Review Letters

We introduce a simple instance of the renormalization group transformation in the Banach space of probability densities. By changing the scaling of the renormalized variables we obtain, as fixed points of the transformation, the L\'evy strictly stable laws. We also investigate the behavior of the transformation around these fixed points and the domain of attraction for different values of the scaling parameter. The physical interest of a renormalization group approach to the generalized central limit theorem is discussed., Comment: 16 pages, to appear in J. Stat. Phys

The papers shows an algorithm to search for approximations of reals to rationals of the form a/b^2 that runs on \sqrt(b) polynomial time steps., Comment: 4 pages, 1 figure

Given a Dirac subbundle and an isotropic subbundle of a Courant algebroid, we provide a canonical method to obtain a new Dirac subbundle. When the original Dirac subbundle is involutive (i.e., a Dirac structure) this construction has interesting applications, for instance to Dirac's theory of constraints and to the Marsden-Ratiu reduction in Poisson geometry., Comment: 10 pages. Paper reformulated and reorganized. Added Section 5.3 on Poisson fibrations

In relation to Hall's conjecture, a new algorithm is presented to search for small k = |x^3 -y^2| values. Thirteen new values of k < x^(1/2) are reported., Comment: 9 pages, 1 figure

Plots of quadratic residues display some features that are analyzed mathematically., Comment: 3 figures