Abstract
Research on magnetic confinement of hightemperature plasmas has the ultimate goal of harnessing nuclear fusion for the production of electricity. Although the tokamak^{1} is the leading toroidal magneticconfinement concept, it is not without shortcomings and the fusion community has therefore also pursued alternative concepts such as the stellarator. Unlike axisymmetric tokamaks, stellarators possess a threedimensional (3D) magnetic field geometry. The availability of this additional dimension opens up an extensive configuration space for computational optimization of both the field geometry itself and the currentcarrying coils that produce it. Such an optimization was undertaken in designing Wendelstein 7X (W7X)^{2}, a large helicalaxis advanced stellarator (HELIAS), which began operation in 2015 at Greifswald, Germany. A major drawback of 3D magnetic field geometry, however, is that it introduces a strong temperature dependence into the stellarator’s nonturbulent ‘neoclassical’ energy transport. Indeed, such energy losses will become prohibitive in hightemperature reactor plasmas unless a strong reduction of the geometrical factor associated with this transport can be achieved; such a reduction was therefore a principal goal of the design of W7X. In spite of the modest heating power currently available, W7X has already been able to achieve hightemperature plasma conditions during its 2017 and 2018 experimental campaigns, producing record values of the fusion triple product for such stellarator plasmas^{3,4}. The triple product of plasma density, ion temperature and energy confinement time is used in fusion research as a figure of merit, as it must attain a certain threshold value before netenergyproducing operation of a reactor becomes possible^{1,5}. Here we demonstrate that such record values provide evidence for reduced neoclassical energy transport in W7X, as the plasma profiles that produced these results could not have been obtained in stellarators lacking a comparably high level of neoclassical optimization.
Main
In the quest for a viable fusion reactor, consideration of the plasma energy balance shows that—regardless of the confinement concept—a minimum value of the fusion triple product, nTτ_{E}, must be attained before netenergyproducing operation becomes possible^{1,5}. Here, n is the fuel density, T is its temperature and τ_{E} is the energy confinement time, defined by the ratio W/P, where W is the stored plasma energy and P is the heating power provided by fusion reactions. The temperature dependence of the fuel’s fusion reactivity provides an additional constraint; for deuterium–tritium fusion, this reactivity falls rapidly below a temperature of 10 keV (≈1.2 × 10^{8} K). High temperatures are thus mandatory in fusion plasmas but must be simultaneously consistent with a tolerable level of energy transport if the required τ_{E} is to be achieved.
Toroidal magnetic confinement of fully ionized fusion plasmas requires that field lines spiral around the minor axis of the torus poloidally as they encircle the major axis toroidally, tracing out magnetic flux surfaces in the course of numerous transits about the device. For a tokamak these toroidal and poloidal components of B are provided, respectively, by planar currentcarrying coils situated outside the plasma and by a toroidal plasma current induced with a central solenoid. The strength of the magnetic field in a tokamak plasma scales inversely with the distance from the major axis of the torus and B ≡ B thus varies along field lines, being largest on the inboard side of the torus and smallest on the outboard side. For strongly magnetized plasmas, both the total energy and the magnetic moment are constants of the particle motion, so that particles having only a small portion of their velocity vector aligned with the magnetic field will become trapped in this variation of B, and in axisymmetric tokamaks they perform socalled banana orbits, which are a consequence of the vertical drift caused by the B × ∇B terms in the particles’ guiding centre equation of motion^{6,7}. In the absence of collisions, these banana orbits experience no net radial displacement on average over the course of their periodic ‘bounce’ motion. However, as in the ‘classical’ case of plasma immersed in a homogeneous magnetic field, Coulomb collisions will cause such particle orbits to undergo a randomwalk diffusive process with a repetition rate linearly proportional to the collision frequency, ν. As the ‘width’ of banana orbits is larger than the particle’s gyroradius, the resulting transport will exceed the classical level and the adjective ‘neoclassical’ is used to signify that the inhomogeneity of B has been accounted for^{8,9}.
Hightemperature fusion plasmas are characterized by low collisionality (ν*, the ratio of collision frequency to banana bounce frequency), which scales as ν* ∝ nT^{−2}. For such plasmas, aside from the large geometrical factor due to tokamak banana orbits^{9}, the neoclassical energy flux through the magnetic surface with a minor radius of r will obey \(V{\prime} {Q}_{{\rm{neo}}}\propto {n}^{2}{T}^{1/2}{B}_{0}^{2}\), where the lefthand side of this expression is the product of magnetic surface area, V′, and fluxsurfaceaveraged neoclassical energyflux density, Q_{neo}, and where B_{0} is the average magnetic field strength at the major radius of the plasma axis (denoted by R_{0}). This is noteworthy as these scalings are identical to those of the classical case, and the temperature dependence is therefore benign.
In the majority of tokamak experiments, energy confinement is found to be worse than predicted by neoclassical theory, which is thought to be predominantly due to plasma turbulence. If the turbulence is of gyroBohm nature, the energy flux will scale with T^{5/2} (ref. ^{10}), and as neoclassical and turbulent losses are additive, the former are generally ignored when assessing the overall energy confinement to be expected in a tokamak reactor. The situation is very different in a hightemperature stellarator plasma, however, as will be discussed next.
The need to and means of reducing the stellarator Q _{neo}
Stellarators have the advantage over tokamaks of producing both the toroidal and poloidal components of their magnetic field with currentcarrying coils external to the plasma, thereby possessing an inherent steadystate capability that tokamaks lack. Numerous possibilities exist for the placement of coils so as to provide the spiraling of field lines needed for the formation of flux surfaces, the most intuitive of which is the use of continuous helical windings as employed in heliotrons such as the Large Helical Device (LHD)^{11}. Strong plasma shaping, however, is more readily produced with a set of modular coils^{12} appropriately twisted into nonplanar forms, and it is this concept that underlies W7X. In either case, the coil shapes combine both toroidal and poloidal excursions so that the magnetic field they produce cannot be axisymmetric and, indeed, an additional corrugation of B arises, which is commonly referred to as the stellarator’s ‘helical’ ripple. Particles that become trapped in these ripples are said to be localized, as they experience only a small variation of poloidal angle between successive reflection points and this leads to a nonzero time average of their vertical drift over this bounce period, unlike the case of tokamak banana orbits^{13}. As a consequence, the collisionless trajectories of localized particles may be able to leave the plasma after numerous reflections, and confinement of such particles will occur only if collisions are frequent enough to limit their radial excursions. Collisions are thus beneficial, and in such a case, typical for electrons in hightemperature stellarator plasmas, the randomwalk diffusive transport becomes inversely proportional to the collision frequency and electrons are therefore said to be in the ‘1/ν regime’. This regime is most remarkable for the very unfavourable temperature dependence of its neoclassical energy flux which scales as \(V{\prime} {Q}_{{\rm{neo}}}\propto {{\epsilon }}_{{\rm{eff}}}^{3/2}{T}^{9/2}{R}_{0}^{1}{B}_{0}^{2}\) (ref. ^{14}). The quantity ϵ_{eff} is referred to as the effective helical ripple for 1/ν transport^{15,16} and is a figure of merit devised to allow comparisons of devices that have different magnetic field geometries. As the name implies, ϵ_{eff} will have the same value as the helicalripple amplitude, ϵ_{h}, for the limiting case in which this amplitude is a constant over the entire flux surface.
The temperature dependence of this result has long been considered a serious hindrance to the prospects of stellarator reactors, and clearly favours operation at the highest collisionality tolerable. High magnetic field strength and large aspect ratio, R_{0}/a (r = a denotes the radius of the plasma’s last closed flux surface), are also beneficial but imply substantial capital cost for reactor construction as such expenditure scales with the stored magnetic energy of the device. The most ‘economical’ option offered by the 1/ν scalings is reduction of the effective helical ripple, which requires appropriate tailoring of the magnetic field so as to diminish the timeaveraged radial drifts experienced by localized particles, and this has become a common goal of stellarator optimization since the inception of W7X^{17}. This device has five field periods (the integer number of times the magnetic field geometry repeats itself when going around the device once in the toroidal direction) and these periods are oriented such that the coil system appears nearly pentagonal when viewed from above (see Extended Data Fig. 1). The largest values of B are located in the pentagon’s ‘corners’, where the field curvature is particularly strong, but the majority of localized particles are thereby trapped in regions that have small B × ∇B, which reduces the associated radial drifts. Making use of this, W7X was optimized for sufficiently reduced neoclassical energy transport that future HELIAS reactors with plasma volumes of 1,500 m^{3} become conceivable, requiring ϵ_{eff} values of a few per cent at most^{18,19,20}.
In addition to causing 1/ν transport, localized particles are also responsible for the appearance of a radial electric field, E_{r}, which must arise in stellarators to establish ambipolarity of the neoclassical particle fluxes (meaning that no net radial current flows in the plasma). This electric field introduces an E × B drift into the particles’ equation of motion, which causes localized particles to drift poloidally with a precession frequency Ω_{E} ≈ E_{r}/(rB_{0}), thereby placing an additional limit on the radial excursion of orbits. This poloidal precession is more important than collisions for particles that have Ω_{E} > ν_{eff}, where ν_{eff} is the frequency with which collisional removal from the ripple occurs. For fusion plasmas, the collision frequency of electrons will exceed that of the ions by roughly two orders of magnitude and the ions will not be subject to 1/ν transport but will instead have their orbits constrained by Ω_{E}, with a value of E_{r} such that the neoclassical ion particle flux is reduced to the value for electrons. For the interested reader, a more detailed presentation of neoclassical results is provided in Methods, where, in particular, it is demonstrated that strong reduction of electron 1/ν transport is also of direct benefit for reducing energy fluxes in the ion channel.
Given these ingredients, the recipe for reducing the sum of electron and ion neoclassical energy losses becomes readily apparent. Paramount is a minimization of electron 1/ν transport, which is best achieved in magnetic fields having small ϵ_{eff} and for plasmas at the highest tolerable collisionality. This strong reduction of the neoclassical electron energy transport is accompanied by a comparable decrease of the electron particle transport and ion neoclassical losses are then beneficially influenced by the ambipolar radial electric field. Reduction of ion transport is thus a knockon effect relying on the actual realization of this predicted value of E_{r}. Theoretically, such expectations are justified as particle transport is known to be intrinsically ambipolar for any turbulence satisfying the gyrokinetic orderings^{21,22}, but experimental verification in W7X is required for certainty.
In addition to reduced neoclassical transport, W7X was also optimized for improved magnetohydrodynamic equilibrium and stability, negligible bootstrap current and good confinement of collisionless fastparticle orbits in plasmas with reactorrelevant pressure profiles^{2}. The W7X coil set was designed with a great deal of flexibility so as to provide access to a large configuration space, having candidates for which each of these goals is weighted to varying degrees^{23}. Of particular interest for neoclassical transport studies is the portion of this space that has extremely small values of ϵ_{eff}(ρ = r/a), an excellent example of which is the ‘standard’ configuration, which has equal currents in all five different types of nonplanar coils. However, this configuration is expected to have larger bootstrap current and poorer fastparticle confinement than are desirable. Both deficiencies can be addressed by choosing a portion of the configuration space that has a larger variation of B along the magnetic axis. Such a ‘highmirror’ configuration is achieved by adjusting the current ratios in nonplanar coils of the same type so that these ratios are largest at the beginning of each field period (at the corners of the pentagon—see Extended Data Fig. 1) and smallest at its midpoint. The larger mirror term increases the fraction of trapped particles, however, leading to ϵ_{eff} values for the highmirror configuration that exceed those of the standard configuration by factors of between 2.5 and 3.5, as can be seen by comparison of the radial profiles plotted in Fig. 1. These values are nevertheless small enough to pose no obstacle to reactor operation and, indeed, all HELIAS reactor studies are based on variants of the highmirror configuration^{19,20}.
Also shown in Fig. 1 are ϵ_{eff} profiles for two reference configurations of LHD^{11}, a heliotron that has been in operation in Toki, Japan, since 1998. Such a comparison is a natural undertaking as W7X and LHD are the largest members of the stellarator family currently in existence. These two devices have nearly the same plasma volume, although the aspect ratio of LHD (R_{0}/a ≈ 6) is more compact than that of W7X (R_{0}/a ≳ 10). At LHD, vertical field coils can be used to shift the plasma column in the radial direction and configurations are typically differentiated according to the position of their major axis in vacuum; configurations with R_{0} = 3.6 m and R_{0} = 3.75 m are depicted here. As is typical for heliotrons, ϵ_{eff}(ρ) is reduced by shifting the plasma inwards, whereas the case with the larger major radius is a good approximation to a stellarator field that has ϵ_{eff} ≈ ϵ_{h}. Neoclassical transport coefficients for both LHD and W7X have been calculated and benchmarked within an international collaboration^{24} and provide a basis for the calculations of neoclassical energy transport presented in the next section.
To visualize the benefit of optimizing the magnetic field geometry for reduced neoclassical transport, collisionless singleparticle orbits of ions in W7X standard and LHD R_{0} = 3.75 m are also provided here in two videos (see Supplementary Videos 1, 2).
Q _{neo} for W7X discharge 20180918.045
In many stellarator plasmas, the electron and ion temperatures are too low to cause neoclassical energy losses of a magnitude relevant for the plasma energy balance. As in tokamaks, this is attributed to turbulent transport, and the global energy confinement exhibited by the two concepts is quite comparable for devices of the same volume (see figure 36 of ref. ^{25}). So far, plasma performance in W7X is commonly limited by turbulence as well, and an experimental assessment of the neoclassical energy confinement in this device therefore requires a plasma scenario for which the turbulent transport is reduced.
In this regard, a substantial transient improvement of the energy confinement has been observed in certain W7X experiments fuelled with hydrogen pellets^{26,27}; the discharge 20180918.045, performed in the W7X standard configuration, provides such an example. Time traces of interest for this experiment are provided in Extended Data Fig. 2, showing that a density ramp was initiated at t = 1.86 s by injecting a series of twentyeight frozen hydrogen pellets into the plasma at a frequency of 30 Hz. In the aftermath of pellet fuelling the central density exceeds 10^{20} m^{−3} and the stored energy measured by a diamagnetic loop reaches a maximum value of W_{dia} = 1.02 MJ with central electron and ion temperatures both in excess of 2.5 keV, although the plasma is heated with only 4.5 MW of electron cyclotron resonance heating (ECRH). Taking fully ionized carbon as the predominant plasma impurity and a core value of the effective charge state equal to the global experimental value, Z_{eff} = 1.4, yields a central triple product of n^{i}T^{i}τ_{E} > 5.3 × 10^{19} m^{−3} keV s (where n^{i} and T^{i} indicate the density and temperature of ions, respectively), which is in the typical range of results achieved in W7X during the highperformance phase of pelletfuelled discharges^{3,27}. The strong temperature dependence of stellarator neoclassical energy transport immediately raises the question of whether such record tripleproduct results are also experimental evidence for the reduction of neoclassical energy losses by optimization of the W7X magnetic field. This discharge is particularly attractive for such an investigation as the 1.02 MJ are maintained for a full energy confinement time of 230 ms, thereby simplifying the plasma energy balance as the ∂W/∂t term appearing in this equation will be of negligible importance by the end of the highenergy phase. For interested readers, further details of this experiment can be found in Methods.
Calculation of the neoclassical losses requires knowledge of the plasma profiles. These are shown in Fig. 2 for t = 3.35 s at the end of the time period in which maximum W_{dia} is reached. The abscissa in these plots, r, is a fluxsurface label, often referred to as the effective minor radius. The red data points depict n^{e} and T^{e}, the density and temperature of electrons, respectively, obtained from Thomson scattering, and electron temperatures obtained with the electron cyclotron emission (ECE) system are shown in black. Values of T^{i} determined with charge exchange recombination spectroscopy (CXRS) are given by the blue circles. See Methods for information concerning these measurements and the error bars associated with them. Profile fits to these experimental data are depicted in Fig. 2 by continuous curves in red for electrons and in blue for ions. The stored energy associated with these profiles is W = 1.01 MJ, in good agreement with the experimentally measured W_{dia}. To account for the uncertainty in experimental profiles, neoclassical results will also be determined here assuming certain variations of the plasma parameters, but, as will be seen, such sensitivity studies do not lead to fundamentally different conclusions.
Before proceeding to the analysis of these experimental results and their interpretation, it should be noted that two tenets underlie this process. First, that quantitative values of the neoclassical fluxes can be accurately determined for specified magnetic equilibria and plasma parameters, and second, that these fluxes represent the minimum value of transport that can be achieved for the given conditions. To substantiate the first of these claims, numerical means of solving the linearized drift kinetic equation and determining neoclassical transport coefficients have been thoroughly benchmarked as part of an international collaboration^{24}. For the second, although there are various examples in the literature of stellarator confinement being consistent with neoclassical theory in the plasma core^{25,28}, there are no published claims of experimental results for which the observed confinement exceeds neoclassical predictions.
The importance of the neoclassical optimization for attaining record fusiontripleproduct results in hightemperature stellarator plasmas can be demonstrated by calculating and comparing the neoclassical energy fluxes associated with the density and temperature profiles of Fig. 2 for configurations with different levels of optimization. This comparison is provided in Fig. 3a for the configurations (b) W7X standard, (c) W7X highmirror, (d) LHD R_{0} = 3.6 m, and (e) LHD R_{0} = 3.75 m, and depicts the sum of electron and ion neoclassical energy fluxes, normalized to the heating power P_{heat} = 4.5 MW of discharge 20180918.045, as a function of normalized minor radius. Of primary interest here is to isolate the influence of the magnetic field configuration on the neoclassical results and, to this end, effects of the plasma volume on confinement have been removed by a slight linear upscaling of both LHD configurations so as to have the same plasma volume as W7X; such scaling leaves the ϵ_{eff} profile and the aspect ratio unchanged. Influence of the magnetic field strength on the results is avoided by setting B_{0} = 2.5 T in all cases. This comparison demonstrates that the plasma profiles of discharge 20180918.045 in the W7X standard configuration are only conceivable in magnetic configurations that have a comparably high level of neoclassical optimization. Indeed, the temperatures attained are high enough to imply a peak neoclassical energy flux nearly commensurate with the heating power for the W7X highmirror configuration, and for the LHD examples these energy fluxes exceed P_{heat} over a portion of the plasma and are thus physically impossible.
Further details of the neoclassical results are provided for each individual configuration in the remaining four frames of Fig. 3. In these plots, constituent contributions to the neoclassical energy fluxes made by electrons, \({Q}_{{\rm{neo}}}^{{\rm{e}}}\), and ions, \({Q}_{{\rm{neo}}}^{{\rm{i}}}\), are shown by the red and blue curves, respectively; their sum is again given by the black ballandchain curve. To illustrate the sensitivity of this sum to variations in plasma parameters, results for the total neoclassical energy fluxes are also depicted with (n^{e}, T^{α}) replaced by (gn^{e}, g^{−1}T^{α}) and varying the scaling factor g in the range 0.9 ≤ g ≤ 1.1, where α = {e, i} indicates the values for electrons and ions, respectively). As the pressure profile is unaffected by this variation it is possible to perform all calculations using the same equilibrium. The values of collisionality remain sufficiently small in all these cases to ensure that the electron neoclassical fluxes are predominantly due to 1/ν transport. Given the strong temperature dependence of the losses in this regime, it is apparent that \({Q}_{{\rm{neo}}}^{{\rm{e}}}\) will attain its largest value for g = 0.9 and then decrease as g increases. The same behaviour is found for \({Q}_{{\rm{neo}}}^{{\rm{i}}}\) although the relative reduction is much weaker and can only be accounted for accurately by enforcing the ambipolarity constraint so as to obtain the correct value of the radial electric field. The sum \({Q}_{{\rm{neo}}}^{{\rm{e}}}+{Q}_{{\rm{neo}}}^{{\rm{i}}}\) is thus a monotonically decreasing function of g over the range of values considered and confirms the theoretical expectations mentioned in the previous section.
As expected qualitatively, the results of Fig. 3 confirm that \(V{}^{{\prime} }({Q}_{{\rm{n}}{\rm{e}}{\rm{o}}}^{{\rm{e}}}+{Q}_{{\rm{n}}{\rm{e}}{\rm{o}}}^{{\rm{i}}})\) is an increasing function of ϵ_{eff} but quantification of this statement is far more difficult. For example, comparing the results of the two W7X configurations, one finds \({Q}_{{\rm{neo}}}^{{\rm{e}}}\) larger in the highmirror case by a factor of five, which corresponds well with the \({{\epsilon }}_{{\rm{eff}}}^{3/2}\) dependence of 1/ν transport. However, the total neoclassical energy flux for W7X highmirror is increased by a factor only somewhat larger than two, owing to the much smaller increase in \({Q}_{{\rm{neo}}}^{{\rm{i}}}\). This example also demonstrates the nonlinearity of the neoclassical fluxes—quantitative accuracy in the determination of these fluxes cannot rely only on comparisons of ϵ_{eff}.
With regard to the higher levels of neoclassical energy transport for the LHD configurations, it should be emphasized that these are provided only to give the reader some idea of how large neoclassical energy fluxes can become at high T^{α} in stellaratortype devices that were not explicitly optimized to reduce neoclassical transport. Such observations have been accounted for in contemporary designs of heliotron reactors by shifting the plasma column further inwards and also adjusting the pitch modulation of the helical windings so as to obtain a tolerable level of neoclassical transport for the envisaged plasma conditions^{29}.
As explained in the previous section, the structure of B is critical in determining the level of electron neoclassical transport and the reduction of losses in the ion channel relies on the appearance of the radial electric field needed to satisfy the ambipolarity constraint. Experimental profiles of E_{r} obtained from Xray imaging crystal spectrometer (XICS) measurements have been published for other examples of pelletfuelled discharges in W7X and have been shown to conform with neoclassical predictions^{30}. The same claim can be made for discharge 20180918.045, which has E_{r} profiles obtained from CXRS in addition to those from XICS (see Extended Data Fig. 3). With this additional information, plasma profiles leading to large values of the triple product in W7X hightemperature plasmas may be taken as experimental evidence for the reduction of neoclassical losses through an appropriate optimization of the confining magnetic field. In particular, the very small values of ϵ_{eff}(ρ), which characterize the W7X standard configuration, are indispensable for achievement of the high electron and ion temperatures of discharge 20180918.045 with a heating power of only 4.5 MW. As the behaviour of the bootstrap current in W7X has already been shown to conform with neoclassical expectations^{31}, reduction of all neoclassical fluxes and flows in line with the W7X optimization goals is also substantiated.
Confirming the importance of the W7X neoclassical optimization for improving plasma performance is only a first step, however, as hightripleproduct phases have been of limited duration in the device so far. It is postulated that—owing to properties of the turbulent transport^{32}—such performance requires the establishment and sustainment of sufficiently steep density gradients. A similar tendency was also observed in the preceding Wendelstein experiment, W7AS, where the appearance of ‘optimum confinement’ conditions was always accompanied by such gradients^{25}. These W7AS discharges were heated using a combination of ECRH and neutral beam injection (NBI), the latter providing the plasma with a strong central particle source, while simultaneously the edge particle source due to recycling neutrals dropped to unusually low levels. Under such conditions, steep density gradients could be maintained throughout the 250ms heating pulse but, although long in comparison to τ_{E}, such pulses were insufficient to claim ‘steadystate’ conditions as the sinks due to neutralparticle pumping by plasmafacing components did not saturate during this time. At W7X, NBI was successfully commissioned during the second portion of the 2017–2018 campaign and will allow future investigations into whether optimum confinement conditions can also be realized in this device during the 10s duration of NBI pulses. To truly test the steadystate perspectives of the HELIAS concept, the ECRH at W7X has been designed to provide the plasma with 1,800 s of continuouswave power and a watercooled, highheatflux divertor is currently being installed in the device to provide the necessary particle and heat exhaust over this period of time. If steep density gradients are indeed the key to improving the confinement of W7X plasmas, density profile tailoring over such time scales will probably need to rely on the capabilities of a new steadystate pellet injector, which should also go into operation during the next campaign.
Methods
Theoretical details
In hightemperature stellarator plasmas, the radial components of the fluxsurfaceaveraged neoclassical particle and energy flux densities (Γ_{neo} and Q_{neo}, respectively) of particle species α may be expressed as^{24}
where n = n(r), T = T(r) and q are the density, temperature and charge of the given species, respectively, E_{r} = E_{r}(r) is the radial electric field and the δ_{ij} are normalized transport coefficients
comprising appropriate combinations of elements of the neoclassical transport matrix
where K ≡ κ/T = mv^{2}/(2T) is the normalized kinetic energy and D is the socalled monoenergetic radial transport coefficient. The terminology ‘radial’ is used here to denote quantities that are oriented perpendicularly to flux surfaces so that the radial coordinate r should be understood as a fluxsurface label.
To understand how the neoclassical energy transport scales with various plasma and configuration parameters, it is sufficient to consider a simple heuristic description of randomwalk diffusion processes having D ∝ ℱ(Δr)^{2}ν_{eff} where ℱ is the fraction of particles participating in the process, Δr is the characteristic step size of such particles and ν_{eff} is the frequency with which a step is taken. For the cases of interest here, the transport is due to the pitchangle scattering portion of the linearized collision operator, allowing one to express the ‘effective’ step frequency as ν_{eff} = ν/ℱ^{2} where ν is the 90°deflection frequency and the ℱ^{−2} enhancement accounts for the fact that scattering through the portion of phase space comprising ℱ occurs more often than scattering through 90°. The heuristic expression for the transport coefficient then simplifies to D ∝ (Δr)^{2}ν/ℱ leaving only Δr and ℱ to be determined.
The variation of B along field lines in toroidal devices causes a B × ∇B drift of guiding centre particle trajectories in this field, having a characteristic velocity v_{d} = κ/(qR_{0}B_{0}) (refs. ^{6,7}). For particles localized in a stellarator’s helical ripple, this drift will lead to a displacement from the flux surface by an amount Δr = min(v_{d}/ν_{eff}, v_{d}/Ω_{E}), depending on whether the effective collision frequency, ν_{eff}, or the E × B precessional frequency, Ω_{E} = E_{r}/(rB_{0}), is larger for the given particle. The first of these cases is typical for electrons in a hightemperature stellarator plasma and will be considered in the remainder of this paragraph. For the simplest of stellarator magnetic fields having a constant helicalripple amplitude over the entire flux surface, ϵ_{h} = ϵ_{h}(r), the fraction of localized particles scales as \( {\mathcal F} \propto {{\epsilon }}_{{\rm{h}}}^{1/2}\) and the heuristic expression for the radial transport coefficient then yields \(D\propto {{\epsilon }}_{{\rm{h}}}^{3/2}{v}_{{\rm{d}}}^{2}/\nu \). Substituting this result into the formula for the L_{ij} and taking the collision frequency proportional to nv^{−3}, one recovers the scaling \({L}_{ij}\propto {{\epsilon }}_{{\rm{h}}}^{3/2}{T}^{7/2}{n}^{1}{({R}_{0}{B}_{0})}^{2}\). This result can be generalized to arbitrary stellarator fields if one replaces ϵ_{h} with the appropriate value of ϵ_{eff}, which explains why this latter quantity is given the name ‘effective helical ripple for 1/ν transport’. To derive the scalings of the neoclassical energy flux in the 1/ν regime, one approximates V′ ∝ rR_{0} and the inverse gradient scale lengths appearing in the braces of equation (1) as a^{−1} to obtain \(V{}^{{\prime} }{Q}_{{\rm{n}}{\rm{e}}{\rm{o}}}\propto {{\epsilon }}_{{\rm{e}}{\rm{f}}{\rm{f}}}^{3/2}{T}^{9/2}{R}_{0}^{1}{B}_{0}^{2}\).
For a fusion plasma having equal species temperatures the collision frequency of electrons will exceed that of the ions by roughly two orders of magnitude. Consequently, for localized ions the characteristic radial step size of the diffusive randomwalk process will be Δr ∝ v_{d}/Ω_{E} and ℱ will combine the fractions of phase space in which particle orbits change their trapping states due to either collisions or drifts. For the stellarator field with ϵ_{h} = ϵ_{h}(r), the former case is characterized by a ‘collisional boundary layer’ of width^{13} \({{\epsilon }}_{{\rm{h}}}^{1/2}{({\nu }_{{\rm{eff}}}/{{\Omega }}_{E})}^{1/2}={(\nu /{{\Omega }}_{E})}^{1/2}\) whereas the latter is signified by ℱ_{tr}, the fraction of ‘transition’ orbits^{33}, which is a largely geometric factor reflecting the topology of the local maxima of B. Substitution of these quantities into the heuristic expression for the transport coefficient then yields for the ions \(D\propto {({v}_{{\rm{d}}}/{\varOmega }_{E})}^{2}\nu {(\sqrt{\nu /{\varOmega }_{E}}+{ {\mathcal F} }_{{\rm{t}}{\rm{r}}})}^{1}\) which is more commonly found in the literature as two separate results for the limiting cases in which one of the terms within the parentheses is far larger than its counterpart, the socalled √ν regime with \(D\propto {v}_{{\rm{d}}}^{2}\,{(\nu /{\varOmega }_{E}^{3})}^{1/2}\) and the ‘ν regime’ with \(D\propto {({v}_{{\rm{d}}}/{{\Omega }}_{E})}^{2}\nu /{ {\mathcal F} }_{{\rm{tr}}}\). Using \({\mathscr{A}}\) to signify those parameters in ν/Ω_{E} that are not of direct relevance to plasmaparameter scalings, the neoclassical ion energy losses are then seen to obey
Although this expression is conveniently compact, it leaves the complicated dependence of E_{r} on plasma and device parameters unspecified. This will be addressed next for the case of most relevance to the reduction of neoclassical transport in W7X.
Unlike the axisymmetric tokamak, neoclassical particle fluxes are not intrinsically ambipolar in a stellarator, and thus a theoretical means of determining the E_{r} profile is provided by enforcing the ambipolarity constraint. In a pure hydrogen plasma (using α = e to denote electrons and α = i for ions) for which n^{e} = n^{i} = n and for which q^{e} = −e and q^{i} = e, where e is the elementary charge, equating \({{\Gamma }}_{{\rm{neo}}}^{{\rm{e}}}\) and \({{\Gamma }}_{{\rm{neo}}}^{{\rm{i}}}\) will yield
Substituting this result back into the neoclassical expressions, one obtains for the particle flux density
and for the energy flux densities
where the species indices are chosen to be [α, β] = [e, i] or [i, e] as appropriate. However, one should recall that the \({L}_{ij}^{{\rm{i}}}\) are dependent on E_{r}, so that profitable use of these equations requires that special circumstances hold. One such example is the fusionrelevant case that has T^{e} = T^{i} = T for which the radial electric field equation becomes
and will yield E_{r} ∝ T for the limiting case in which electron 1/ν transport has been sufficiently reduced to satisfy \({L}_{11}^{{\rm{e}}}/{L}_{11}^{{\rm{i}}}\ll 1\). In the same limit, the particle flux density is well approximated by
making it evident that electrons are the ratecontrolling species. Additionally, one notices that the factor outside the braces in \({Q}_{{\rm{neo}}}^{\alpha }\) simplifies to \(nT{L}_{11}^{{\rm{e}}}\) and that the largest normalized transport coefficients associated with the temperature gradients are \({\delta }_{22}^{{\rm{e}}}+{\delta }_{12}^{{\rm{i}}}{\delta }_{21}^{{\rm{e}}}\) for electrons and \(({\delta }_{22}^{{\rm{i}}}{\delta }_{12}^{{\rm{i}}}{\delta }_{21}^{{\rm{i}}}){L}_{11}^{{\rm{i}}}/{L}_{11}^{{\rm{e}}}\) for ions. At first glance, the largeness of \({L}_{11}^{{\rm{i}}}/{L}_{11}^{{\rm{e}}}\) might lead one to expect a considerable difference in the size of these two terms, but the strong temperature dependence of 1/ν transport makes the \({\delta }_{ij}^{{\rm{e}}}\) considerably larger than their \({\delta }_{ij}^{{\rm{i}}}\) counterparts. Indeed, under the assumption of ‘pure’ regimes, the L_{ij} may be expressed in terms of the socalled gamma function, and one finds \({\delta }_{12}^{{\rm{e}}}=7/2\), \({\delta }_{21}^{{\rm{e}}}=5\) and \({\delta }_{22}^{{\rm{e}}}=45/2\) whereas \(1/2\le {\delta }_{12}^{{\rm{i}}}\le 5/4\), \(2\le {\delta }_{21}^{{\rm{i}}}\le 11/4\) and \(3\le {\delta }_{22}^{{\rm{i}}}\le 99/16\) where the smaller values hold in the ν regime and the larger for the √ν regime^{34}. When these limits apply one obtains \(25\le {\delta }_{22}^{{\rm{e}}}+{\delta }_{12}^{{\rm{i}}}{\delta }_{21}^{{\rm{e}}}\le 115/4\) and \(2\le {\delta }_{22}^{{\rm{i}}}{\delta }_{12}^{{\rm{i}}}{\delta }_{21}^{{\rm{i}}}\le 11/4\) so that the δ_{ij} combination of relevance for electrons is an order of magnitude larger than its counterpart for ions. The neoclassical energy transport of electrons and ions will thus be of similar magnitudes and the strong reduction of \({L}_{11}^{{\rm{e}}}\) is clearly seen to be of benefit to both species.
Details concerning the calculation of neoclassical energy fluxes
For the magnetic configurations considered in this work all equilibria have been determined with the Variational Moments Equilibrium Code (VMEC)^{35} and subsequently expressing the results in terms of Boozer flux coordinates^{36}. This provides all information concerning the magnetic field needed as input by the Drift Kinetic Equation Solver (DKES)^{37}, which is used to prepare a dataset of monoenergetic transport coefficients covering the entire range of ν* and E_{r} values relevant for determining the L_{ij} given any combination of density and temperature producing the pressure profile of the VMEC equilibrium. The radial electric field profile is determined selfconsistently by using numerical rootfinding techniques to determine solutions to the nonlinear ambipolarity constraint.
For the W7X calculations, equilibria for the pressure profile of discharge 20180918.045 at t = 3.35 s have been used, whereas the LHD results presented here have been obtained for vacuum equilibria. The LHD equilibria thus ignore the outward shift of the plasma column due to nonzero plasma pressure—the socalled Shafranov shift^{38}—which leads to an increase of the neoclassical energy losses^{39,40}. This deleterious effect can be counteracted to an extent by using the vertical field coils to shift the plasma axis back to its vacuum position, but a deformation in the shape of flux surfaces remains, which serves to degrade the neoclassical confinement. The neoclassical energy losses calculated here for LHD are thus ‘best case’ results. By contrast, the W7X optimization had the explicit goal of using finiteplasmapressure effects to its benefit; for the highmirror configuration, ϵ_{eff} decreases monotonically as the pressure increases, whereas for the standard configuration, smalltomodest pressure has little influence on this quantity.
The results plotted in Fig. 3 also allow a rough assessment of the reactor prospects for these four configurations with regard to their neoclassical energy confinement. For a deuteriumtritium plasma meant to produce a power of 3 gigawatts thermal (GWth) and that has 600 MW of αparticle heating, it is necessary to increase the temperatures of discharge 20180918.045 by a factor of six and the density by a factor of three, assuming that reactorsized versions of W7X and LHD have their dimensions increased by a factor of four and their magnetic field strength doubled. Scaling up the configurations leaves ϵ_{eff} unchanged so that the electron 1/ν energy fluxes will increase according to \({T}^{9/2}{R}_{0}^{1}{B}_{0}^{2}\), and thus by a factor of roughly 200. Scaling the maximum electron energy fluxes of Fig. 3 by this factor, one obtains 75 MW for W7X standard, 340 MW for W7X highmirror, 640 MW for LHD R_{0} = 3.6 m and 1,320 MW for LHD R_{0} = 3.75 m. For stellarator reactor plasmas, ion neoclassical energy fluxes are at least as large as those of the electrons but even the somewhat extreme assumption of \({Q}_{{\rm{neo}}}^{{\rm{i}}}=2{Q}_{{\rm{neo}}}^{{\rm{e}}}\) would make the total neoclassical energy fluxes only a fraction of the αparticle heating for W7X standard. At first glance, the reactor prospects of W7X highmirror appear questionable, but one must realize that this backoftheenvelope estimate ignores changes to the configuration due to increased plasma pressure. In the current example, the normalized plasma pressure, β ∝ nTB^{−2}, is increased by a factor of 4.5 and, when this is accounted for, ϵ_{eff} in the highmirror configuration is reduced by enough that \({Q}_{{\rm{neo}}}^{{\rm{e}}}\) would drop by a factor of two. For the LHD cases it is necessary to choose different scaling factors; the R_{0} = 3.6 m case becomes viable with a factor five larger size and the somewhat reduced temperature that this allows. (This assumes that the vertical field coils of the heliotron are used to compensate the Shafranov shift as the neoclassical transport in such a highpressure equilibrium would otherwise become intolerable).
It should also be pointed out that the optimization of W7X was undertaken roughly thirty years ago and that great improvements in the projected fastparticle confinement of HELIAS have been made in the intervening years. Such improvements have simultaneously reduced the ϵ_{eff} values from the per cent level of W7X to the per mille level in new reactor candidates, thereby decreasing neoclassical energy transport for prospective fusion plasmas to very small levels. This development, coupled with the predominance of turbulent transport in the W7X experiment, have focused recent theoretical and numerical efforts on the further optimization of the HELIAS concept to also contend with this transport channel. Such efforts are still in their infancy, however, and it is not yet possible to foresee what combination of optimized magnetic field and plasma conditions will best reduce turbulent transport and what influence such a combination will have on the other HELIAS optimization goals.
Background and experimental details
Since its inception, W7X has been conceived as a means of demonstrating the attractive properties that advanced stellarators offer as prospective fusion reactors. Ultimately, this envisions not only the achievement of hightripleproduct plasma operation but the ability to maintain such operation over time scales far in excess of those characterizing all physical processes of the plasma and the plasmafacing components. To enable such ‘steadystate’ scenarios, electron cyclotron resonance heating (ECRH) of W7X plasmas is available using ten gyrotrons, each of which is capable of providing continuouswave power over a time interval of 1,800 s (ref. ^{41}). The gyrotron frequency of 140 GHz is resonant at the second harmonic for a magnetic field strength of B = 2.5 T with cutoff densities of n_{c} = 1.2 × 10^{20} m^{−3} for waves with extraordinarymode (X2) polarization and 2n_{c} for ordinarymode (O2) polarization. At these densities and for W7X plasma volumes of nearly 30 m^{3} the collisional transfer of energy from electrons to ions should be excellent, so that T^{i} ≈ T^{e} can be expected in spite of lacking a direct means of heating the ions. Such operational conditions also mimic qualitatively those of a reactor where the heating power of fusion αparticles goes chiefly to electrons, which subsequently heat deuterium and tritium ions by collisional energy exchange.
With regard to plasmafacing components, steadystate operation of W7X will only become possible after the installation of a watercooled, highheatflux (HHF) divertor capable of providing the necessary particle and power exhaust throughout the course of such experiments. The installation of this divertor with all its intricate plumbing requirements was started at the end of 2018, with completion scheduled for the end of 2021. The experimental results considered here were obtained in the latter portion of 2018, during the second half of a campaign in which initial experience with divertor operation in W7X was acquired through the use of an uncooled test divertor unit (TDU) having the same geometry that its HHF successor will have. Pellet fuelling was already observed to improve plasma performance during the first portion of this TDU campaign^{3,26} and such results motivated further investigations once wall conditioning had been greatly improved following ‘boronization’ of the device, early in the second half of the TDU campaign. On many occasions, considerable increases of the plasma diamagnetic energy, W_{dia}, were measured in the aftermath of pellet fuelling and the discharge 20180918.045 is used here to illustrate such results.
Time traces for this discharge are plotted in Extended Data Fig. 2. Following a short startup phase in Xmode to achieve plasma breakdown, the ECRH power (first plot in the figure) is launched exclusively with Omode polarization so as to avoid large amounts of reflected power should pellet fuelling lead transiently to densities above the X2 cutoff. Also shown by the dotted curve is the radiated power measured with bolometer arrays; this power originates in roughly equal portions from the confined plasma and from the region outside the last closed flux surface, which includes the large magnetic islands that are the basis of the W7X divertor concept. Initially, the lineintegrated electron density, ∫dℓn^{e}, is held at a value of roughly 3 × 10^{19} m^{2}, which is too low to promote effective collisional energy transfer from electrons to ions (where the dℓ appearing in the integrand is the differential length along the line of sight). Consequently, the 3.2 MW of ECRH used during this phase produces plasmas that have T^{e} ≫ T^{i} centrally as the core temperature measurements from Thomson scattering and the Xray imaging crystal spectrometer (XICS) diagnostics show (third plot in the figure). A density ramp is initiated 1.86 s into the discharge by injecting a series of 28 frozen hydrogen pellets into the plasma at a frequency of 30 Hz using a socalled blower gun^{42}. Although the fuelling efficiency of the first few pellets is poor, noticeable improvement occurs thereafter^{26} and increments in the lineintegrated density caused by individual pellets become clearly discernible with ∫dℓn^{e} slightly exceeding 12 × 10^{19} m^{−2} when the series ends at t = 2.8 s. At this time point, the ECRH power is increased to 4.5 MW to maintain T^{e} at values sufficient for good O2 absorption. This is further aided at W7X by a multipass launch scheme which uses specially prepared reflecting surfaces to redirect unabsorbed ECRH power back into the plasma. With this scheme, raytracing simulations predict that more than 95% of the launched ECRH power will be deposited in the plasma during the highperformance phase of discharge 20180918.045, a value that has been confirmed by the analysis of strayradiation measurements from the experiment taken with ‘sniffer’ probes.
The pellets obviously provide the plasma with a particle source but their thermalization also introduces a considerable energy sink into the electron energy balance particularly as the thermalization process also strongly heats the ions for the given initial conditions^{43}. In response, T^{e} falls and the central electron temperature has dropped below 2 keV—nearly to the level of central T^{i}—when the increase in ECRH power occurs. To this time point, pellet fuelling has been accompanied by a rise in the diamagnetic energy from W_{dia} = 0.40 MJ at t = 1.8 s to W_{dia} = 0.68 MJ at t = 1.8 s. With the larger heating power, the time rate of change of W_{dia} then increases considerably and the diamagnetic energy attains its maximum value of 1.02 MJ after another 400 ms. In the absence of pellet fuelling, the lineintegrated density decreases throughout this time, making it evident that a temperature increase must also have taken place and, indeed, the central values exceed 2.5 keV for both electrons and ions. Thus, an experimental situation has arisen that has high density and high temperatures simultaneously, a situation ideal for testing the efficacy of the W7X optimization regarding reduction of the neoclassical transport. The discharge shown here is particularly attractive for such an investigation as the 1.02 MJ are maintained for 230 ms which corresponds to one energy confinement time given the 4.5 MW of ECRH used to heat the plasma. This simplifies considerations of the energy balance as the ∂W/∂t term appearing in this equation will be of negligible importance by the end of the highenergy phase.
ECRH power deposition was not highly localized for discharge 20180918.045, as the electron cyclotron waves were launched with ordinarymode polarization. Nevertheless, the great majority of the heating power is deposited within ρ < 0.4 at the time point of the analysis. At this radius, and in the neighbouring vicinity at least out to ρ = 0.55, the sum of electron and ion neoclassical energy fluxes exceeds 25% of the heating power and is thus of relevance when considering the energy balance. A larger share of this balance remains unaccounted for, however, even after deducting the observed energy losses due to radiation. This shortfall indicates that turbulent transport must still be of importance for explaining the experimental results, especially in the plasma periphery where low temperatures lead to negligible neoclassical energy transport. In spite of the level of turbulent transport needed to fully explain this discharge, the high temperatures achieved are nevertheless responsible for an experimentally relevant neoclassical energy flux, even though the W7X standard configuration has extremely small values of the effective helical ripple. Assuming such temperatures had been achieved in configurations without a comparably high level of neoclassical optimization quickly leads to the physically impossible result of neoclassical energy fluxes exceeding the heating power, as is demonstrated in the comparison of Fig. 3.
Diagnostics
The plasma energy has been measured with a diamagnetic loop, located in the (toroidal) symmetry plane of one of the five field periods^{44}. It encircles the plasma and is equipped with four compensation coils, which are also located inside the vacuum vessel and directly attached to the diamagnetic loop itself. These do not encircle the plasma and can therefore be used to compensate measurements of the main loop for errors due to eddy currents in the adjacent vacuum vessel as well as fluctuations of externally driven currents in the main superconducting magnetic field coils.
The lineintegrated electron density was determined by a singlechannel dispersion interferometer employing a CO_{2} laser measuring at a fundamental wavelength of 10.6 μm (ref. ^{45}). The probing beam passes through the plasma twice by making use of a corner cube reflector; the singlepass path length through the plasma is roughly 1.3 m. The statistical error for ∫dℓn^{e} is generally given as 10^{18} m^{−2} and there was an additional systematic error ≤4 × 10^{18} m^{−2} during the portion of the experimental campaign during which discharge 20180918.045 was performed.
The Thomson scattering system at W7X employs three YAG lasers to provide full profiles of electron density and temperature at a rate of 30 Hz. Scattered light is collected by two invessel optical systems and routed to polychromators outside the torus hall via optical fibres. Fortytwo spatial channels are available, each employing five interference filters to provide spectral resolution^{46,47}. Spectral calibration is done ex vessel by withdrawing the optical systems and illuminating a diffusereflecting screen with a supercontinuum light source^{46}. Absolute calibration for density measurements was obtained by performing antiStokes rotational Raman scattering in nitrogen. Bayesian analysis is used to determine the most probable electron temperature and density for each Thomson scattering volume, as well as their uncertainties and crosscorrelation. Error bars in Fig. 2 depict the width of the 95% confidence interval divided by four, which corresponds to the onestandarddeviation interval of a normal distribution; see ref. ^{47} for full details.
Electron cyclotron emission is used at W7X to determine the electron temperature using a 32channel heterodyne radiometer probing X2 emission^{48}. The T^{e} data plotted in Fig. 2 give ECE values of the electron temperature averaged over a centred 20ms time window, with the error bars depicting the standard deviation of these values and also accounting for systematic calibration uncertainties, largely due to the unknown thermal drifts of the radiometer sensitivity. An optimized lineofsight is realized across the 3D plasma shape by means of an in vessel Gauss telescope. The radiometer channels are calibrated relative to each other and absolute calibration is carried out using a second ‘identical’ Gauss telescope outside the torus that views a hotcold light source^{49}. For the results presented here, the absolute calibration has been scaled up by a factor of 1.58 to account for the relative mismatch with Thomson data on the day of the experiment. Direct interpretation of the radiometer signals as a local temperature measurement would require that ideal blackbody conditions be fulfilled in the plasma, whereas in reality, nonthermal ‘hot’ components of the electron distribution function appear in the spectrum as well. A further complication arises in semitransparent plasmas owing to a degradation of spatial resolution. The T^{e} profile is therefore inferred by applying Bayesian analysis with forward modelling of the radiation transport in the plasma; for full details, see ref. ^{50}.
The CXRS diagnostic at W7X makes use of a highétendue spectrometer to provide measurements of T^{i} at 51 spatial locations as well as E_{r} at 25 of these. Given the prevalence of carbon in the device’s plasmafacing components, this element is always present as an impurity in W7X plasmas and the visible charge exchange line C^{VI} at 529.07 nm provides the strongest possible signal. Active charge exchange measurements are enabled by using short ‘blips’ of neutral beam injection (NBI) for diagnostic purposes and subtracting the passive spectrum observed before and after the NBI from this data. A complete description may be found in ref. ^{51}.
The XICS diagnostic is based on spectroscopic analysis of emission from highly charged argon impurities that are seeded into the plasma in trace amounts for diagnostic purposes^{52}. The XICS system records a 1D image of lineintegrated spectra, from which the ion temperature is found by measuring the Dopplerbroadened width of the emission lines. The local plasma parameters are found by using a tomographic inversion based on a known VMEC equilibrium^{53}. Detailed descriptions of the XICS diagnostic on W7X can be found in ref. ^{54} although improvements in the diagnostic analysis have since been implemented, including compensation for spherical aberrations and the subpixel distribution of photons on the detector. At the time point analysed for this discharge, the XICS values of T^{i} exceed those of CXRS by 150 to 200 eV and indicate T^{i} > T^{e} for r > 0.2 m despite the fact that no direct heating of the ions occurs using ECRH. This is at odds with energy balance considerations, however, which argue for T^{i} ≈ T^{e} outside the region of power deposition. These expectations are better fulfilled by the CXRS data and, given the strong dependence of the neoclassical losses on temperature, it has therefore been decided to err on the side of caution by using only the CXRS values of T^{i} for the profile fits used in the calculations of neoclassical fluxes. The time evolution of the T^{i} profile measured by CXRS has also been confirmed through the measurements of the XICS system.
Plasma radiation is measured by two bolometer cameras—each with arrays of detectors that have a 5 μm blackened goldfoil absorber—installed in one of the midfieldperiod toroidalsymmetry planes of W7X, and with spatial resolution of between 3 and 4 cm (ref. ^{55}). Lineintegrated signals from 65 channels are used to obtain radiation intensity distributions by tomographic reconstruction with ‘relative gradient smoothing’ as regularization functional (to be published). Fluxsurfaceaveraged radial emissivity profiles are then derived by averaging these 2Demissivity distributions in the poloidal direction. The total radiated power loss is a linear interpolation of radiation from the observation volume to that of the entire plasma volume. Toroidal variations of the radiation strength are not considered, an assumption supported by the results of edge modelling.
Data availability
The data depicted in the plots of this paper and other findings of this study are available from the corresponding author upon reasonable request.
Change history
12 October 2021
A Correction to this paper has been published: https://doi.org/10.1038/s4158602104023y
References
 1.
Wesson, J. Tokamaks 2nd edn (Oxford Univ. Press, 1997).
 2.
Grieger, G. et al. Modular stellarator reactors and plans for Wendelstein 7X. Fusion Technol. 21, 1767–1778 (1992).
 3.
Pedersen, T. S. et al. First results from divertor operation in Wendelstein 7X. Plasma Phys. Contr. Fusion 61, 014035 (2019).
 4.
Wolf, R. C. et al. Performance of Wendelstein 7X stellarator plasmas during the first divertor operation phase. Phys. Plasmas 26, 082504 (2019).
 5.
Lawson, J. D. Some criteria for a power producing thermonuclear reactor. Proc. Phys. Soc. B 70, 6–10 (1957).
 6.
Alfvén, H. On the motion of a charged particle in a magnetic field. Ark. Mat. Astr. Fys. 27A, 1–20 (1940).
 7.
Northrop, T. G. The guiding center approximation to charged particle motion. Ann. Phys. 15, 79–101 (1961).
 8.
Galeev, A. A. & Sagdeev, R. Z. Transport phenomena in a collisionless plasma in a toroidal magnetic system. Sov. Phys. JETP 26, 233–240 (1968).
 9.
Hinton, F. L. & Hazeltine, R. D. Theory of plasma transport in toroidal confinement systems. Rev. Mod. Phys. 48, 239–308 (1976).
 10.
Connor, J. W. Invariance principles and plasma confinement. Plasma Phys. Contr. Fusion 30, 619–650 (1988).
 11.
Iiyoshi, A., Fujiwara, M., Motojima, O., Ohyabu, N. & Yamazaki, K. Design study for the Large Helical Device. Fusion Technol. 17, 169–187 (1990).
 12.
Wobig, H. & Rehker, S. A stellarator coil system without helical windings. In Proc. 7th Symp. Fusion Technology 345–357 (1972).
 13.
Galeev, A. A., Sagdeev, R. Z., Furth, H. P. & Rosenbluth, M. N. Plasma diffusion in a toroidal stellarator. Phys. Rev. Lett. 22, 511–514 (1969).
 14.
Ho, D. D. M. & Kulsrud, R. M. Neoclassical transport in stellarators. Phys. Fluids 30, 442–461 (1987).
 15.
Dommaschk, W., Lotz, W. & Nührenberg, J. MonteCarlo simulations of neoclassical transport in stellarators. Nucl. Fusion 24, 794–799 (1984).
 16.
Beidler, C. D. & Hitchon, W. N. G. Ripple transport in helicalaxis advanced stellarators: a comparison with classical stellarator/torsatrons. Plasma Phys. Contr. Fusion 36, 317–353 (1994).
 17.
Beidler, C. D. et al. Physics and engineering design for Wendelstein 7X. Fusion Technol. 17, 148–168 (1990).
 18.
Painter, S. L. & Lyon, J. F. Transport analysis of stellarator reactors. Nucl. Fusion 31, 2271–2290 (1991).
 19.
Beidler, C. D. et al. The Helias reactor. In 16th Int. Conf. Fusion Energy Vol. 3, 407–418 (International Atomic Energy Agency, 1996).
 20.
Beidler, C. D. et al. The Helias reactor HSR4/18. Nucl. Fusion 41, 1759–1766 (2001).
 21.
Sugama, H., Okamoto, M., Horton, W. & Wakatani, M. Transport processes and entropy production in toroidal plasmas with gyrokinetic electromagnetic turbulence. Phys. Plasmas 3, 2379–2394 (1996).
 22.
Helander, P. Theory of plasma confinement in nonaxisymmetric magnetic fields. Rep. Prog. Phys. 77, 087001 (2014); corrigendum 81, 099501 (2018).
 23.
Geiger, J. et al. Physics in the magnetic configuration space of W7X. Plasma Phys. Contr. Fusion 57, 014004 (2015).
 24.
Beidler, C. D. et al. Benchmarking of the monoenergetic transport coefficients — results from the International Collaboration on Neoclassical Transport in Stellarators (ICNTS). Nucl. Fusion 51, 076001 (2011).
 25.
Hirsch, M. et al. Major results from the stellarator W7AS. Plasma Phys. Contr. Fusion 50, 053001 (2008).
 26.
Baldzuhn, J. et al. Pellet fueling experiments in Wendelstein 7X. Plasma Phys. Contr. Fusion 61, 095012 (2019).
 27.
Bozhenkov, S. A. et al. Highperformance plasmas after pellet injections in Wendelstein 7X. Nucl. Fusion 60, 066011 (2020).
 28.
Dinklage, A. et al. Intermachine validation study of neoclassical transport modelling in medium to highdensity stellaratorheliotron plasmas. Nucl. Fusion 53, 063022 (2013).
 29.
Goto, T. et al. Conceptual design of a compact helical fusion reactor FFHRc1 for the early demonstration of yearlong electric power generation. Nucl. Fusion 59, 076030 (2019).
 30.
Pablant, N. A. et al. Investigation of the neoclassical ambipolar electric field in ionroot plasmas on W7X. Nucl. Fusion 60, 036021 (2020).
 31.
Dinklage, A. et al. Magnetic configuration effects on the Wendelstein 7X stellarator. Nat. Phys. 14, 855–860 (2018); correction 14, 867 (2018).
 32.
Alcusón, J. A. et al. Suppression of electrostatic microinstabilities in maximumJ stellarators. Plasma Phys. Contr. Fusion 62, 035005 (2020).
 33.
Mynick, H. E. Effect of collisionless detrapping on nonaxisymmetric transport in a stellarator with radial electric field. Phys. Fluids 26, 2609–2615 (1983).
 34.
Maaßberg, H., Beidler, C. D. & Simmet, E. E. Density control problems in large stellarators with neoclassical transport. Plasma Phys. Contr. Fusion 41, 1135–1153 (1999).
 35.
Hirshman, S. P., van Rij, W. I. & Merkel, P. Threedimensional free boundary calculations using a spectral Green’s function method. Comput. Phys. Commun. 43, 143–155 (1986).
 36.
Boozer, A. H. Guiding center drift equations. Phys. Fluids 23, 904–908 (1980).
 37.
van Rij, W. I. & Hirshman, S. P. Variational bounds for transport coefficients in threedimensional toroidal plasmas. Phys. Fluids B 1, 563–569 (1989).
 38.
Shafranov, V. D. Equilibrium of a toroidal pinch in a magnetic field. Sov. At. Energy 13, 1149–1158 (1963).
 39.
Mynick, H. E. The effect of finite β on stellarator transport. Phys. Fluids 28, 1139–1147 (1985).
 40.
Wakasa, A., Murakami, S. & Oikawa, S. Study of neoclassical transport in LHD plasmas by applying the DCOM/NNW neoclassical transport database. Plasma Fusion Res. 3, S1030 (2008).
 41.
Erckmann, V. et al. Electron cyclotron heating for W7X: physics and technology. Fus. Sci. Technol. 52, 291–312 (2007).
 42.
Dibon, M. et al. Blower gun pellet injection system for W7X. Fusion Eng. Des. 9899, 1759–1762 (2015).
 43.
Aleynikov, P., Breizman, B. N., Helander, P. & Turkin, Y. Plasma ion heating by cryogenic pellet injection. J. Plasma Phys. 85, 905850105 (2019).
 44.
Rahbarnia, K. et al. Diamagnetic energy measurement during the first operational phase at the Wendelstein 7X stellarator. Nucl. Fusion 58, 096010 (2018).
 45.
Brunner, K. J. et al. Realtime dispersion interferometry for density feedback in fusion devices. J. Instrum. 13, P09002 (2018).
 46.
Pasch, E. et al. The Thomson scattering system at Wendelstein 7X. Rev. Sci. Instrum. 87, 11E729 (2016).
 47.
Bozhenkov, S. A. et al. The Thomson scattering diagnostic at Wendelstein 7X and its performance in the first operation phase. JINST 12, P10004 (2017).
 48.
Hirsch, M. et al. ECE diagnostic for the initial operation of Wendelstein 7X. In 20th Joint Workshop on Electron Cyclotron Emission and Electron Cyclotron Resonance Heating (EC20) (eds Oosterbeek, J. et al.) 03007 (2019).
 49.
Höfel, U. et al. Bayesian modeling of microwave radiometer calibration on the example of the Wendelstein 7X electron cyclotron emission diagnostic. Rev. Sci. Instrum. 90, 043502 (2019).
 50.
Höfel, U. Bayesian Analysis of Electron Cyclotron Emission Measurements at Wendelstein 7X. PhD thesis, Technische Universität Berlin (2020); https://doi.org/10.14279/depositonce9621.
 51.
Ford, O. P. et al. Charge exchange recombination spectroscopy at Wendelstein 7X. Rev. Sci. Instrum. 91, 023507 (2020).
 52.
Bitter, M. et al. Objectives and layout of a highresolution Xray imaging crystal spectrometer for the Large Helical Device. Rev. Sci. Instrum. 81, 10E328 (2010).
 53.
Pablant, N. A. et al. Tomographic inversion techniques incorporating physical constraints for line integrated spectroscopy in stellarators and tokamaks. Rev. Sci. Instrum. 85, 11E424 (2014).
 54.
Langenberg, A. et al. Prospects of Xray imaging spectrometers for impurity transport: recent results from the stellarator Wendelstein 7X. Rev. Sci. Instrum. 89, 10G101 (2018).
 55.
Zhang, D. et al. First observation of a stable highly dissipative divertor plasma regime on the Wendelstein 7X stellarator. Phys. Rev. Lett. 123, 025002 (2019).
Acknowledgements
This work has been carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014–2018 and 2019–2020 under grant agreement no. 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission.
Funding
Open access funding provided by Max Planck Society.
Author information
Affiliations
Consortia
Contributions
J.B., S.A.B., Ye.O.K. and H.P.L. planned and coordinated the experiment described here. T.S. provided the ECRH power and its deposition profile. K.J.B. calculated the unabsorbed ECRH power using sniffer probe data. S.A.B., H.D., G.F., E.P. and E.R.S. supplied Thomson data. K.J.B. provided interferometer data. M.H. and U. Höfel supplied ECE data. A.L. and N.A.P. supplied XICS data. O.P.F. supplied CXRS data. U. Hergenhahn, M.K., S.K., A.P. and J. Svensson provided experimental data for the effective charge state or helped in analysis of this data. T.A., U.N., K.R., J. Schilling and H.T. supplied measurements of the diamagnetic energy. D.Z. provided results for the radiated power. A.A., O.P.F., N.A.P. and A.L. provided experimental profiles of E_{r}. M.N.A.B., U. Höfel, A.L., N.A.P. and F.W. performed data analysis. J.B. provided data concerning the pellet series. J.G. carried out VMEC calculations. H.M.S. created the DKES datasets and performed the calculations of neoclassical fluxes. Y.T. determined the profiles of effective helical ripple. C.D.B. created the figures (except for Extended Data Fig. 1) and is the principal author of this paper with contributions to the body of the text from P.H., T.S.P. and R.C.W. and with text concerning the diagnostics from K.J.B., M.H., N.A.P., K.R. and D.Z. M.D. provided the depiction of the W7X nonplanar coils. Particle orbit videos were prepared by M.B. and R.K. All authors provided feedback and contributed to improving the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Peer review information Nature thanks Dennis Whyte and the other, anonymous, reviewer(s) for their contribution to the peer review of this work. Peer reviewer reports are available.
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
Extended Data Fig. 1 Nonplanar coil system of W7X, viewed from above.
The complete set of superconducting coils also includes 20 planar coils, four in each field period, which are used to change the rotational transform and/or shift the plasma column; these are not shown as they remain currentfree for the configurations considered here.
Extended Data Fig. 2 Time traces for W7X discharge 20180918.045.
In the top frame, the launched ECRH power is shown, as well as the radiated power measured by bolometers (dotted curve). The second frame plots the lineintegrated density measured by an interferometer. The third frame depicts ‘core’ electron (red) and ion (blue) temperatures from Thomson scattering and XICS measurements, respectively. The diamagnetic energy trace during this discharge is given in the bottom frame. A series of 28 pellets is injected into the plasma at a frequency of 30 Hz during the time phase indicated in grey. W_{dia} values exceeding 1.02 MJ are recorded during the phase indicated in yellow.
Extended Data Fig. 3 Profiles of the radial electric field as a function of the normalized plasma radius.
Experimental results from CXRS (blue) and XICS (red) measurements are compared to the theoretical expectations from the ambipolarity constraint (black) for the plasma profiles of Fig. 2. Error bars depict one standard deviation in the evaluation of the measurements. The XICS analysis is described in ref. ^{30}; documentation of the methodology used to evaluate the CXRS measurements is in preparation.
Supplementary information
Supplementary Information
This file contains detailed notes regarding Supplementary Videos 1 and 2.
Video 1
Localized ion trajectory in the W7X standard configuration – see Supplementary Information document for detailed description.
Video 2
Localized ion trajectory in the LHD configuration having a major radius of 3.75 m  see Supplementary Information document for detailed description.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Beidler, C.D., Smith, H.M., Alonso, A. et al. Demonstration of reduced neoclassical energy transport in Wendelstein 7X. Nature 596, 221–226 (2021). https://doi.org/10.1038/s4158602103687w
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s4158602103687w
Further reading

On the brink of a new era in nuclear fusion R&D
Nature Reviews Physics (2022)

Observation of a reducedturbulence regime with boron powder injection in a stellarator
Nature Physics (2022)
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.