
==== Front
Nat Commun
Nat Commun
Nature Communications
2041-1723
Nature Publishing Group UK London

39245748
52182
10.1038/s41467-024-52182-z
Article
Stable Deuterium-Tritium plasmas with improved confinement in the presence of energetic-ion instabilities
http://orcid.org/0000-0003-0900-5564
Garcia Jeronimo jeronimo.garcia@cea.fr

1
Kazakov Yevgen 2
Coelho Rui 3
http://orcid.org/0000-0003-0482-0981
Dreval Mykola 4
de la Luna Elena 5
http://orcid.org/0000-0002-4815-3407
Solano Emilia R. 5
http://orcid.org/0000-0002-9608-280X
Štancar Žiga 6
http://orcid.org/0000-0002-6114-0539
Varela Jacobo 78
Baruzzo Matteo 9
http://orcid.org/0000-0001-7947-2841
Belli Emily 10
Bonofiglo Phillip J. 11
Candy Jeff 10
Maggi Costanza F. 6
Mailloux Joelle 6
http://orcid.org/0000-0001-6491-8759
Mazzi Samuele 1
Ongena Jef 2
Ruiz Juan R. 12
http://orcid.org/0000-0002-1858-4046
Poradzinski Michal 6
Sharapov Sergei 6
http://orcid.org/0000-0002-7220-8092
Zarzoso David 13
JET contributors
1 grid.457341.0 CEA, IRFM, Saint-Paul-lez-Durance, France
2 grid.270680.b Laboratory for Plasma Physics, LPP-ERM/KMS, EUROfusion Consortium member, TEC Partner, Brussels, Belgium
3 grid.9983.b 0000 0001 2181 4263 Instituto de Plasmas e Fusao Nuclear, Instituto Superior Técnico, Universidade de Lisboa, Lisboa, Portugal
4 https://ror.org/00183pc12 grid.425540.2 0000 0000 9526 3153 National Science Center Kharkiv Institute of Physics and Technology, Kharkiv, Ukraine
5 grid.420019.e 0000 0001 1959 5823 Laboratorio Nacional de Fusión, CIEMAT, Madrid, Spain
6 https://ror.org/0361bwx64 grid.9689.e 0000 0001 0683 2623 United Kingdom Atomic Energy Authority, Culham Campus, Abingdon, UK
7 https://ror.org/03ths8210 grid.7840.b 0000 0001 2168 9183 Universidad Carlos III de Madrid, Leganes, Madrid Spain
8 https://ror.org/00hj54h04 grid.89336.37 0000 0004 1936 9924 Institute for Fusion Studies, Department of Physics, University of Texas at Austin, Austin, TX USA
9 Dip.to Fusione e Tecnologie per la Sicurezza Nucleare, ENEA C. R. Frascati, via E. Fermi 45, Frascati (Roma), Italy
10 https://ror.org/03ngjpk76 grid.192673.8 0000 0004 0634 455X General Atomics, PO Box 85608, San Diego, CA USA
11 https://ror.org/03vn1ts68 grid.451320.1 0000 0001 2151 1350 Princeton Plasma Physics Laboratory, Princeton, NJ USA
12 https://ror.org/052gg0110 grid.4991.5 0000 0004 1936 8948 Rudolf Peierls Centre for Theoretical Physics, University of Oxford, Oxford, UK
13 https://ror.org/035xkbk20 grid.5399.6 0000 0001 2176 4817 Aix Marseille Univ, CNRS, Centrale Med, M2P2, Marseille, France
8 9 2024
8 9 2024
2024
15 784613 12 2023
29 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 licence, and indicate if you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence 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 licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Providing stable and clean energy sources is a necessity for the increasing demands of humanity. Energy produced by Deuterium (D) and Tritium (T) fusion reactions, in particular in tokamaks, is a promising path towards that goal. However, there is little experience with plasmas formed by D-T mixtures, since most of the experiments are currently performed in pure D. After more than 20 years, the Joint European Torus (JET) has carried out new D-T experiments with the aim of exploring some of the unique characteristics expected in future fusion reactors, such as the presence of highly energetic ions in low plasma rotation conditions. A new stable, high confinement and impurity-free D-T regime, with reduction of energy losses with respect to D, has been found. Multiscale physics mechanisms critically determine the thermal confinement. These crucial achievements importantly contribute to the establishment of fusion energy generation as an alternative to fossil fuels.

The use of clean energy sources is essential for the humankind. Here, authors show new experiments carried out with Deuterium-Tritium fusion reactions. Results show that energy production by such reactions can be more efficient than expected, confirming fusion as an alternative to fossil fuels.

Subject terms

Magnetically confined plasmas
Experimental nuclear physics
issue-copyright-statement© Springer Nature Limited 2024
==== Body
pmcIntroduction

Modern societies are eager to increase their energy resources, which are currently largely provided by fossil fuels. In this context, plasmas, i.e. the fourth state of matter, which are characterized by the presence of free-charged particles, can provide a carbon-free energy source through fusion reactions of light atom nuclei. For that purpose, plasmas must be well confined and reach high pressure in order to overcome the electrostatic repulsion of particles. Magnetic containment is one of the most promising routes toward this goal. The tokamak concept, in which plasmas are confined both by a high toroidal electric current, Ip, and a toroidal magnetic field, BT, is particularly promising. Tokamaks have made great progress from their initial proposal1,2. Plasmas with good thermal confinement have been achieved in high-performance tokamak discharges when turbulence-driven energy losses (one of the main physics mechanisms by which plasmas lose their thermal confinement) are strongly reduced. This is the case for the H-modes3 or internal transport barrier (ITB) regimes4, in which steep pressure gradients are formed, leading to high plasma pressure. Recently, the EAST5 and KSTAR6 tokamaks have reported significant advances. Long-sustained high-performance plasmas were obtained in conditions of low turbulent heat transport with simultaneous avoidance of dangerous magnetohydrodynamic (MHD) bursts, called edge localized modes (ELMs)7, which are characteristic of H-modes. Such bursts can lead to rapid expulsion of edge plasma and hence to high levels of heat and particle flux to the tokamak wall.

On the one hand, these advances have clarified the route towards a potential fusion energy commercial reactor. On the other hand, the results obtained are not enough to provide a clear insight into how future energy-producing reactor plasmas are expected to behave. The reason is that in significant contrast to plasmas produced nowadays, formed almost exclusively by pure Deuterium (D), energy-producing plasmas will use Deuterium and Tritium (T) to produce 14.1 MeV neutrons and 3.5 MeV 4He nuclei (alpha particles). The presence of T has been identified as a potential source of significant changes with respect to pure D plasmas, particularly in terms of turbulence or MHD characteristics8,9. Furthermore, the presence of a small and yet very energetic population of alpha particles, which will provide the main self-heating mechanism through collisions to the thermal plasma, leads to conditions that are mostly unexplored and have been identified as potentially detrimental. The highly energetic fusion-born alpha particles mainly transfer energy to electrons (rather than to thermal ions) by collisions, thus providing strong electron heating. Conditions in which the electron temperature, Te, is significantly higher than the ion temperature, Ti, have been identified as leading to unfavorable destabilization of turbulence and clamping of the ion temperatures at very low values, as shown in plasmas with pure electron heating by means of electromagnetic waves5,10. Moreover, alpha particles do not provide significant torque and therefore lead to low toroidal rotation, reducing rotation-driven suppression of turbulence11,12. Finally, they can resonantly destabilize Alfvén waves that can produce stochastic transport of alpha particles and hence reduction of fusion energy generation13–16.

These novel characteristics are in contrast to present-day experiments which are typically heated by neutral beam injection (NBI) at energies  ∼ 100 keV. NBI dominantly heats ions, not electrons, and provides significant torque. Unlike fusion-generated alpha particles, they produce a sizeable fast-ion density fraction with Ti/Te > 1. In such conditions, it is well known that high confinement can be obtained6.

The Joint European Torus (JET) tokamak17, which is the only tokamak in the world capable of operating with T, has undergone a new experimental campaign in D-T, DTE218, with the aim of providing solid evidence on the characteristics of D-T plasmas. A world record of fusion energy was produced in JET at the end of 202119 by developing H-mode plasmas with ELMs. However, such an impressive result was obtained under conditions of high NBI power. In this paper, we describe DTE2 plasmas that, rather than maximizing the production of fusion energy, aim to capture other important characteristics expected in future fusion reactors, such as the simultaneous development of conditions with dominant total heating fraction transferred to the electrons, low input torque, and the triggering of Alfvén wave instabilities. Since the simultaneous development of plasmas with dominant alpha particle heating and low NBI input power is not possible at JET, the expected conditions in future fusion tokamak reactors were reproduced by using external heating with electromagnetic waves resonating at the ion cyclotron frequency (ICRF). In these conditions, low external torque is applied while generating ions in the MeV energy range and significant energetic ions-related perturbations. These characteristics were not studied in the DTE2 fusion energy record, nor in the first D-T campaign in JET (DTE1)9 or TFTR8 in the 90s, since in both cases the maximization of the fusion power by using high NBI power was mostly explored and in such conditions plasmas with Ti/Te > 1 and high ion heating and rotation are obtained. Therefore, the study presented in this paper represents the first time that this path has been pursued in D-T plasmas.

We show that very good properties are achieved in terms of energy confinement and stability. The type of confinement expected in baseline reactor plasmas is obtained because of the better energy confinement in D-T with respect to the same conditions in D. In particular, for the electrons, the energy losses by transport are low, which allows temperatures of about 110 million K. For the ions, the core heat transport is significantly reduced in D-T compared to D when instabilities generated by the energetic ions are observed. This type of plasma provides an integrated solution for future tokamak reactors since, despite the development of an H-mode, deleterious ELMs are avoided. These results suggest improved energy confinement in future D-T plasma conditions.

Results

D-T plasmas development and comparison to D

Several plasma configurations and D-T concentrations were explored in JET to cover a wide range of possible configurations in future tokamak reactors. To minimize the external torque and hence the toroidal rotation, the plasmas presented in this study were mostly heated with ICRF, which ensures that a low external torque is applied. ICRF heating can accelerate ions to MeV energies, and hence it was used as a proxy for alpha particle-heated dominated plasmas. In these plasmas, 1% of H was used as a minority wave resonator. Additional lower heating power levels of NBI were used with the aim of reaching high temperatures. The discharge #99896, with major radius, R0 = 3 m, and minor radius, a = 0.91 m, shown in Fig. 1A, represents the type of plasma performed. The chosen D-T concentration was  ∼ 50%D–50%T as it is expected to deliver the maximum fusion power in future tokamak fusion reactors. A summary of the main global characteristics of this discharge is shown in Supplementary Table 1 in the Supplementary Information.Fig. 1 Main characteristics of the D-T discharge #99896.

A Time evolution of discharge #99896, with toroidal current Ip = 1.9 MA, magnetic field BT = 2.75 T, and q95 = 4.5, heated mainly with ICRF power, PICRF = 4.5 MW. The NBI power, PNBI ∼ 3.5 MW, was also injected with deuterium beams, before 9 s, and tritium beams, after 9 s. The radiated power, Prad represents 60% of the total input power. The power produced by D-T fusion reactions obtained reaches a maximum of Pfus ∼ 0.5 MW. B Time evolution of edge fluctuations as obtained from the BeII line emission from the inner divertor. C Time evolution of βN, defined as βN = βaBT/Ip[%] with β the ratio between magnetic and thermal pressure and a the plasma minor radius, H98(y, 2) and fGr=n¯e/nGr the Greenwald fraction with n¯e the average density and nGr the Greenwald density defined as nGr = Ip/πa2. D Accumulated NBI and ICRF radial power deposition for the discharge #99896. E Time evolution of magnetic perturbations detected by the Mirnov coils.

During the phase with ICRF and NBI heating, the thermal confinement time of this plasma, τ, calculated excluding the contribution of energetic ions, is the reference in the ITER baseline conditions as H98(y, 2) = τ/τIPB98 ≥ 1 (Fig. 1C) with τIPB98 the energy confinement time predicted by the IPB98(y,2) scaling20. This result is obtained under low rotation conditions, since the Mach number at ρ = 0.5, M = vtor/cs ∼ 0.15, with vtor the plasma rotation, cs=Te/mi the sound speed and mi the Deuterium ion mass, is lower than that expected in ITER D-T plasmas21.

As shown in Fig. 1B these results are obtained in the presence of small edge fluctuations rather than ELMs with βN = 1.2 and fGr ∼ 0.45 (Fig. 1C). Analysis using the TRANSP suite of modeling codes22 indicates that 56% of the total auxiliary heating is transferred to electrons as shown in Fig. 1D. The high electron heating is confirmed by the observation that Te > Ti, notably in the region ρ < 0.4 with Te/Ti ∼ 1.4 on the magnetic axis. On the contrary, during the NBI-only phase, H98(y, 2) = 0.7, indicating a lower confinement compared to the phase with injected ICRF heating. In terms of power produced by D-T fusion reactions, Pfus ∼ 0.5 MW was obtained. Finally, the radiated power reaches 60% of the total power injected and is located at the plasma separatrix close to the X-point, which indicates that the accumulation of core impurities, usually a concern in high confinement plasmas in a metallic wall environment such as in JET23, is avoided.

The good thermal confinement of discharge #99896 is obtained in the presence of generally-believed deleterious electromagnetic perturbations over a wide range of frequencies when ICRF heating is added, as shown in Fig. 1E. We can identify perturbations generated by the interplay between energetic ions and Alfvén waves, such as toroidal Alfvén eigenmodes (TAE)24,25 and reversed-shear Alfvén eigenmodes (RSAE)26. Fishbone instabilities27,28 related to the interaction between energetic ions and MHD are also present. All of this activity corroborates the presence of highly energetic ions in the plasma. MHD activity and in particular neoclassical tearing modes (NTM)29,30 are also found.

An equivalent discharge, #100871, in terms of input power, radiated power, and density, with pure D was also carried out to compare to the characteristics of D-T discharges. The match was successful, including a similar pattern of magnetic instabilities. The comparison between the temperature and electron density profiles between D-T and D is shown in Fig. 2A, B. Although the electron density is nearly identical in D-T and D, both Te and Ti are higher in D-T for the entire plasma radius. In particular, in the plasma core, an increase in Ti slope is obtained. Since the heating deposition profiles are very similar in D-T and D in the plasma inner core, as verified with TRANSP, such a change in Ti could have its origin in differences in thermal transport. Indeed, a significant change is found, as demonstrated by calculating the power balance heat diffusivity for ions, χi, and electrons, χe, (see Fig. 2C). In D-T, χi shows a drop in the plasma core starting from ρ = 0.4 and reaching χi,D−T/χi,D ∼ 0.5 at ρ = 0.2. The electron heat diffusivity, χe, remains very low for D-T and D, which leads to very peaked Te, similar to ITB plasmas observed in other high electron heated plasmas such as the super I-mode developed in EAST5. Importantly, χi,D−T/χe,D−T ∼ 2 while it is doubled for pure D indicating deteriorated ion thermal confinement.Fig. 2 Comparison between the D-T discharge #99896 and the D counterpart #100871.

A Comparison between Ti and Te for the D-T discharge #99896 and the D counterpart #100871. Ti is measured by the charge-exchange technique on impurity ions. Te is obtained by means of LIDAR and high-resolution Thomsom scattering (HRTS). An average over 8.5 s-8.7 s is performed. ρ is defined as the square root of the normalized toroidal magnetic flux. Shaded error bars represent the standard deviation of the time-averaged signals and the systematic diagnostic uncertainties. B Comparison between electron density, ne, for the D-T discharge #99896 and the D counterpart #100871. ne is measured with HRTS. Shaded error bars represent standard deviation of the time-averaged signals and the systematic diagnostic uncertainties. C Comparison between χi and χe obtained by power balance analysis for the D-T discharge #99896 and the D counterpart #100871. Shaded error bars represent standard deviation.

D-T density-ratio scans were performed to evaluate the impact of T on plasma characteristics, as shown in Supplementary Figs. 1 and 2. As an example, the discharge #99817 was performed in conditions similar to #99896 but with 85% T fraction and BT = 3.7 T. A summary of the main global characteristics of this discharge is shown in Supplementary Information Table 1. The total heating to electrons under these conditions was 70% as obtained from TRANSP, with Pfus ∼ 1 MW and Te reaching 110 million K and Ti 60 million K with only 8 MW of input power. Similar to discharge #99896, a broad range of energetic-ion instabilities was obtained. Compared to the case with  ∼ 50%D–50%T, χi,D−T/χe,D−T is reduced to  ∼ 1, which means that the presence of T in the plasma is a key player in reducing ion thermal energy losses due to heat transport.

Core instabilities analyses

Beyond large-scale MHD and Alfvén-driven fluctuations, small-scale microturbulence driven by the Ion Temperature Gradient (ITG) instability31 is one of the major threats to plasma thermal confinement. In JET plasmas, ITG modes, destabilized at ion-gyroradius scales,  ∼1 cm, are usually responsible for enhancing core radial energy transport even in strongly electron-heated plasmas32. Instead, energy transport driven by instabilities at electron-gyroradius scales,  ∼ 0.1 mm, i.e. by electron temperature gradient modes (ETG), is found not to contribute significantly to overall outward transport in the plasma core33. In this section, an analysis of the fluctuations found in plasma #99896 is given, and their role in producing good thermal confinement is clarified.

A Fourier analysis of the measured magnetic fluctuations is shown in the Supplementary Information Fig. 3. A full range of activity is found at both low and high frequencies. In the frequency range 1 kHz  < f < 40 kHz, modes with toroidal mode number n = 0 − 7 are detected, while for 120 kHz  < f < 200 kHz the modes detected cover n = − 5 to n = 6. However, nonlinear mode-mode interactions are at the origin of some of these fluctuations. This is demonstrated by performing a mode-mode bi-coherence analysis34. A nonlinear interplay between high-frequency TAE and low-frequency NTM is detected, as shown in Fig. 3A. In particular, toroidal mode numbers n = 1 and n = 2 among other interactions are a result of such nonlinear interplays (this is shown in detail in the Supplementary Information). A necessary condition for such an interaction to occur is that the radial locations of TAE and NTM are close to each other. This requirement is further verified by the analysis of the location of the different perturbations in real space using several techniques, as clarified in the methods section. RSAE and fishbones are located inside q = 1, at ρ ∼ 0.25, with q the safety factor defined as q ≡ dΨt/dΨp with Ψt the toroidal magnetic flux and Ψp the poloidal magnetic flux. RSAE are destabilized very close to the magnetic axis. TAE and NTM are located just outside q = 1, at 0.25 < ρ < 0.45.Fig. 3 Core plasma fluctuations characteristics.

A Bicoherence analysis of the perturbations found in discharge #99896. The analysis is performed at t = 7.7 s. B Logarithmic power of the density fluctuations as obtained from reflectometry at major radius R ∼ 3.36 m, ρ ∼ 0.35, and t = 8.4 s, for the D-T discharge #99896 with only NBI heating or with full NBI and ICRF heating and comparison to the pure D discharge #100871. C Electrostatic potential, Φ, fluctuations obtained for the D-T discharge #99896 by the global code FAR3D when considering two species of energetic ions, H and D, accelerated by the ICRF power. The yellow circle represents the q = 1 surface. Inside q = 1, a n = 1 perturbation is obtained which is identified as a fishbone instability. Outside q = 1, TAEs are obtained.

Further evidence of radial and temporal overlap of MHD and TAE beyond q = 1, located at R ∼ 3.25 m, is obtained from density fluctuations using a reflectometry diagnostic35 at R ∼ 3.36 m. It is shown in Fig. 3B that high-frequency fluctuations for both D-T and D are detected in the TAE and NTM frequency ranges. Furthermore, except for the frequencies corresponding to TAE and NTM, density fluctuations are lower for D-T than for D, notably in the typical range of ITG fluctuations f < 150 kHz. This supports the improvement in confinement for the thermal ions shown in Section 2 although further analyses are required in order to fully characterize turbulence reduction by means of radial correlation from reflectometry. Density fluctuations do not increase significantly when 4.5 MW of ICRF power is added on top of NBI power, except at the TAE and NTM frequencies. This is important because, in general, turbulence increases with increasing input power leading to the so-called thermal confinement degradation with input power12.

The origin of the destabilization of magnetic fluctuations by energetic ions has been analyzed with the global gyrofluid code FAR3D36. Two energetic ion species, H and D, were considered because they are both ICRF-accelerated by means of the first and second harmonic absorption. Their characteristics are obtained from TRANSP. The frequency and location of the fishbone and TAE perturbations obtained from linear simulations with FAR3D agree with experimental data as shown in the Supplementary Figs. 4 and 5 in the Supplementary Information, which clarifies that the perturbations are destabilized by ICRF-accelerated ions. Nonlinear simulations including both energetic ion species are shown in Fig. 3C. The radial extension of the electrostatic potential perturbation, 2.4 m  < R < 3.5 m, agrees well with the experimental location and clearly shows that the energetic-ion-induced perturbations extend up to mid-radius, coinciding with the radial extension of the decrease in thermal energy transport losses in D-T compared to D.

In addition to the nonlinear interplay between TAE and NTM, other complex and multi-scale interactions are of paramount interest and critically determine the performance. It is found with FAR3D that alpha particles do not destabilize any perturbation due to their low density, however, nonlinear interplay with fishbones, located inside q = 1 at ρ ∼ 0.25, can induce radial transport and losses of alpha particles, partially depleting the plasma axis of such particles, as shown in Fig. 4A. This is experimentally corroborated, as shown in Fig. 4B, by means of the fast ion loss detector (FILD)37. In the initial phase of the discharge, when the activity of the fishbone is especially strong, alpha particle losses are detected, as can be seen at t ∼ 8 s. In the later phase of the discharge, fishbone activity is reduced in intensity and no further losses are detected with origin in the fishbone perturbation, but rather in the NTM at f ∼ 15 kHz. Unlike fishbones, no alpha losses with origin in TAE are detected in both the FAR3D simulations and the experiment. Importantly, it is found in FAR3D that losses with origin on ICRF accelerated protons are much lower than those from alpha particles in the presence of fishbones. This is an important result, as it shows that alpha particles are very sensitive to magnetic perturbations even if they are not at the origin of such instabilities.Fig. 4 Alpha particle transport and losses in D-T.

A Alpha particle density profile as calculated with TRANSP assuming no alpha particle transport (nα,eq) and comparison to the profile obtained from the FAR3D code after the full development of the fishbone instability (nα,final). B Alpha particle loss frequency spectrum obtained by using the fast ion loss detector (FILD) with channels that are receptive to 3.5 MeV alpha particles.

Therefore, it becomes clear that although magnetic perturbations do not prevent the reach of high confinement, they can lead to loss of fusion power. These results show the critical interaction between magnetic perturbations and alpha particles, and we conclude that it is essential to control such interactions to produce high fusion power in future tokamak reactors.

Nevertheless, magnetic perturbations induced by energetic ions can also lead to beneficial effects that may have a positive impact on reducing energy losses by heat transport. This is the case of the interplay with the so-called zonal flows, i.e. thermal plasma flows with f ∼ 0 and poloidal and toroidal perturbation mode numbers, n, m = 0, 0. Zonal flows were theoretically predicted38 and also obtained in dedicated numerical simulations in the presence of energetic ions36,39,40. They are known to reduce transport driven by ITG turbulence41, in particular in the presence of energetic MeV ions as shown in D-3He plasmas32,42.

The generation of zonal flows is studied for discharge #99896 by analyzing the energy transfer between fishbone and TAE perturbations and the thermal plasma. The 2D pattern of n = 0 structures for zonal poloidal flows, as obtained from FAR3D nonlinear simulations, are shown in Fig. 5A, B, D, E considering two different energies for the energetic ions. Clearly, fishbones and TAE drive zonal flows, with higher intensity with increased energetic ion energy, and hence stronger instability drive. To study the influence of T on the zonal flow generation, an additional simulation with FAR3D is done artificially replacing T by D and keeping the total amount of thermal ions constant. As shown in Fig. 5C, F, although the generation of zonal flow by fishbone and TAE is also present in pure D, the intensity of zonal flow is lower than in D-T, which could explain the lower turbulent transport found in D-T compared to D, as hinted previously in numerical studies performed in D-T plasmas when turbulent transport is close to threshold43,44.Fig. 5 Zonal flow generation by energetic particle instabilities.

2D pattern of n = 0, m = 0 structures of zonal poloidal flows, Vth(0, 0), for the TAE and fishbone instabilities. Vth(0, 0) is defined as Vth(0, 0) = Er(0, 0) × BT with Er(0, 0) the n, m = 0, 0 component of the perturbed radial electric field. The dependence of the zonal flow intensity on the perturbation strength is studied by scanning the energetic ion equivalent temperature (Tf) using two values, Tf = 1 MeV and Tf = 500 keV in the FAR3D code. Zonal flow generation increases with increasing perturbation intensity for both TAEs (A, B) and fishbones (D, E). The radial extension of zonal flow activity coincides with the extension of the two perturbations. The zonal flow activity in D-T is compared to the one in pure D (C and F) by artificially replacing T by D in FAR3D. The intensity of the zonal flow is lower in D than in D-T.

The role of T and energetic ions on the good plasma confinement of discharge #99896 is further analyzed by performing simulations with CGYRO45. CGYRO solves the gyrokinetic-Maxwell equations46 to obtain the electrostatic and electromagnetic fluctuations and corresponding turbulent energy transport. Local simulations are performed at ρ = 0.31, as it is beyond the inverse radius of the sawtooth and is a location in which TAE is detected. The simulations are carried out by including and excluding energetic ions as separate species in addition to the electrons, D and T species. Due to the low energetic ion density compared to the electron density,  ∼3%, the growth rates obtained in linear simulations are nearly unaffected in the ITG scales by the presence of energetic ions as shown in Fig. 6A, B. However, low ky modes with TAE frequency are destabilized. These results show that the type of plasma found in JET is different to other plasmas dominated by energetic ions effects, such as the FIRE mode6, which is characterized by strong turbulence reduction with energetic ion dilution and linear effects47.Fig. 6 Gyrokinetic analysis of core turbulence reduction in D-T.

A Growth rate, γ, and B frequency, ω, spectrum obtained from linear simulations with the CGYRO code. ky is the binormal wavenumber normalized to the proton sound gyroradius ρs. C Energy flux obtained from gyrokinetic simulations performed with the CGYRO code for discharge #99896 at ρ = 0.31. The simulations are performed including and excluding the energetic ion component. Values of the ion thermal energy flux deduced from power balance in TRANSP (black horizontal dashed line) are only obtained when the energetic ion component is included in the simulations as a separate species. The total thermal ion energy flux obtained in D-T including energetic ions is compared to the one obtained assuming that all the thermal ions are D while keeping the rest of the parameters fixed. The energy flux in pure D is significantly higher than in D-T (a zoom of those bars is displayed in the inset on the top right). Simulations with a/LTFI=0 and TFI/Te = 5.25 are performed to linearly stabilize the low ky energetic ion mode while keeping the energetic ions in the simulations. Fluxes obtained with stabilized mode and energetic ions cannot reproduce the experimental fluxes.

Regarding non-linear effects, as shown in Fig. 6C, the energy flux, Q, obtained when including energetic ions is close to turbulence threshold and it is comparable to the values calculated from power balance analysis from TRANSP. Importantly, such a strong reduction in thermal energy flux compared to the case without energetic ions is accompanied by a high increase in zonal flow shearing activity, ωE×B, which is more than ten times higher when energetic ions are included, thus confirming the results obtained with FAR3D.

Two extra simulations with CGYRO were performed to reveal the physics mechanism by which turbulent transport is strongly reduced in the presence of energetic ions. In both cases, the aim is to keep the energetic ions in the simulations but to stabilize the energetic ion mode at low ky by assuming either a/LTFI=0 or by reducing the energy of the energetic ions four times with respect to the standard case, down to TFI/Te = 5.25 as shown in Fig. 6A, B. As shown in Fig. 6C, the presence of energetic ions without an energetic ion mode has some stabilizing effect on the thermal fluxes, but the matching of experimental fluxes is only obtained when the energetic ion mode is destabilized.

Furthermore, there is a clear asymmetry between the transport obtained in D and in T, with the T transport systematically lower than that for D, χi,T ∼ 0.8χi,D. Such a difference has an important consequence on the total flux in D-T compared to that of pure D. This is numerically analyzed by performing an alternative simulation in which the T ions are artificially considered as D, thus performing a pure D simulation. The turbulent energy flux ratio QD−T/QD = 0.67 is similar to the power balance obtained for the discharges #100871 in D and #99896 in D-T at the same radial location, Q99896/Q100871 = 0.71. This result confirms expectations from purely numerical simulations performed for D-T plasmas43,48,49. From the numerical point of view, global effects from profile shearing were investigated in CGYRO and found to negligibly affect the thermal fluxes at the radial location studied.

Regarding turbulence in the ETG scales, multiscale simulation capturing both energetic ion modes and ETG is very computationally challenging due to extreme spatial resolution requirements and thus beyond the scope of this work. However, for the case without energetic ions, preliminary multiscale simulations capturing both the ion and electron gyroradius scales indicate that the nonlinear ETG transport is suppressed by the ion-scale fluctuations.

In summary, the analyses of core plasma fluctuations indicate an optimum route towards the generation of fusion power in D-T tokamak plasmas whereby energetic ion instabilities remain in conditions of negligible or weak alpha particle transport, while they can induce thermal energy transport reduction by means of zonal flows, notably in the presence of T. Importantly, this is obtained under conditions of some energetic ions characteristics relevant to ITER burning plasmas, e.g. the energetic ion density is similar to the one expected in ITER,  ∼1%48, avoiding energetic ion dilution as usually happens in strong NBI heated plasmas.

Pedestal formation in D-T

As depicted in Fig. 7A and B, the D-T discharge #99896 shows the formation of an H-mode with a steep edge temperature gradient, i.e. a pedestal, at ρpol ∼ 0.95. The temperature of the D counterpart is lower at the same location, while the density is nearly identical. The fact that the pedestal pressure is higher with increasing isotope mass has been observed in H-mode plasmas with ELMs (ELMy H-mode). However, in ELMy H-mode regimes with typical type-I ELMs50–52, an increase in edge density rather than temperature is obtained with increasing isotope mass53.Fig. 7 Pedestal formation in D-T and comparison to D.

A Comparison between edge ne for the D-T discharge #99896 and D discharge #100871. ρpol is defined as the normalized poloidal flux. B Comparison between edge Te for the D-T discharge #99896 and D discharge #100871. The vertical dashed line represents the location of the plasma separatrix. The profiles are obtained from HRTS averaged in the time window 8.5 s-8.9 s for #99896 and 8.4 s-8.8 s for #100871. The evaluation of the error bars in panels A and B is done by deriving the expected signal levels at a given temperature and by calculating the standard deviation based on the photoelectron statistics, the plasma background light variation, and the detector noise.

With the aim of further investigating the origin of the pedestal found in discharge #99896, it is compared to two D-T discharges, one with a clear transition to H-mode with only NBI power and the development of typical type-I ELMy H-mode and another, with NBI and ICRF heating that remains in L-mode, i.e. without temperature pedestal. As shown in Fig. 8A, the edge density of the discharges #99896 and the one in L-mode is nearly identical, demonstrating that the density of discharge #99896 remains in L-mode, yet the pedestal formation is evident as the edge temperature nearly reaches that of the ELMy plasma (Fig. 8B). Importantly, discharge #99896 has no ELMs as shown in Fig. 8C in which the divertor oscillations from BeII divertor emission are compared to those obtained in L-mode and with ELMs. Clearly, the fluctuations are closer to those obtained in L-mode.Fig. 8 Pedestal characteristics in D-T.

Comparison between discharge #99896 and discharge #99502, with PNBI = 12.5 MW, in H-mode with type-I ELMs, and #99776, in L-mode, with PNBI = 5.4 MW and PICRF = 3.3 MW both obtained at Ip=2.5 MA, BT=3.7 T and q95=4.5. A Edge ne. B Edge Te. The profiles are obtained from HRTS averaged in the time window 8.5 s-8.9 s for #99896, 7.3 s-7.6 s for #99502, and 8.6 s-8.9 s for #99776. The evaluation of the error bars in panels A and B is done by deriving the expected signal levels at a given temperature and by calculating the standard deviation based on the photoelectron statistics, the plasma background light variation, and the detector noise. C Comparison of the BeII line emission from the inner divertor for the same discharges. The vertical dashed line represents the location of the plasma separatrix.

The spontaneous generation of plasmas with no ELMs and a pedestal for the temperature is systematic in these kinds of plasmas performed in D-T at different Ip and BT. It is apparent from the available data that, in those plasma conditions, the input power is close to the L to H-mode transition power threshold, but below the power required to fully develop ELMs. This is supported by the fact that during the phase where the temperature pedestal is sustained an n = 0 coherent mode can be observed in the magnetic sensors with a frequency of 5 kHz. This mode, known in JET as M-mode54, is typically detected in JET immediately after the L to H-mode transition. Long phases with an M-mode present are typically observed in the L-H transition at low density in JET, where pedestal dynamics similar to that described here have been identified55. This aspect was further studied in the particular case of discharge #99896, for which it was found that an additional 3 MW of input power lead to the formation of an ELMy H-mode regime. Regarding the magnetic configuration in which these results were obtained, it was used the standard B × ∇B in JET with the ion ∇B drift towards the dominant X-point. This configuration is known to be ‘favourable’ in terms of power requirements for the access to H-mode.

These plasmas show similarities with the I-mode regime56–58, found in D, which is characterized by the formation of a pedestal for the temperature, while the density remains in L-mode. Similar to the I-mode, the no-ELM regime described here has good impurity transport properties, with no impurity accumulation. However, no signs of the so-called weakly coherent mode (WCM), typically found in I-mode plasmas, have been detected.

A scenario towards D-T burning plasmas

Compared to the more typical H modes that develop ELMs and are heated with NBI power, the D-T discharge #99896 provides an attractive alternative with similar thermal energy, but obtained at lower input power, Ip and BT. Such a feature could lead to a more economical and simpler tokamak design. This is shown in Fig. 9A–C by comparing the discharge #99896 and D-T H-mode discharge #99501, at higher PNBI, Ip, and BT and heated with pure NBI power. Clearly, the density is higher at the edge for discharge #99501, since a pedestal is formed for both the density and the temperature, and yet the total thermal energy content of the two discharges, Wth, is similar, Wth ∼ 2.4 MJ. The reason is that the lower edge density is compensated by the higher core temperatures and lower core energy transport losses as shown in Fig. 9D by comparing χi, which is nearly ten times lower in the plasma core for discharge #99896.Fig. 9 Scenario performance in D-T.

Comparison between discharge #99896, obtained at PNBI = 3.5 MW, PICRF = 4.5 MW, Ip = 1.9 MA, BT = 2.75 T, and the D-T discharge #99501, obtained at PNBI = 9.5 MW, Ip = 2.5 MA, BT = 3.7 T. A ne. B Te. C Ti. Shaded error bars represent the standard deviation of the time-averaged signals and the systematic diagnostic uncertainties in panels A–C. D χi/χGB. χGB is the GyroBohm diffusivity defined as χGB=Te3/2mp1/2/(e2BT2a) with mp the proton mass, e the electron charge, and a the plasma minor radius. Shaded error bars represent standard deviation.

Discussion

The path towards commercial fusion reactors, although better understood in recent decades, still poses physics and technological uncertainties. The size, magnetic configuration, type of confinement, or power exhaust techniques expected in future tokamak devices are not fully established. Therefore, it is of fundamental importance to further clarify a safe and clear path toward the generation of efficient energy by means of fusion reactions. This is especially important for D-T plasmas. The presence of T, the generation of a high neutron rate at 14.1 MeV energy, or the presence of a significant population of alpha particles, are all characteristics of future fusion reactors in burning D-T plasmas that are not present in ubiquitous pure-D discharges. Studying the impact of such differences is critical in order to properly characterize how D-T fusion reactors might behave. In particular, T can have a strong impact on confinement, impurity generation, and stability, whereas alpha particles can lead to significant destabilization of magnetic perturbations and provide electron heating.

The JET tokamak has recently conducted a new D-T campaign after the first ones were developed at TFTR8 and JET9 in the 90’s. In view of clarifying key physics elements that will characterize future D-T plasmas, several scientific directions have been explored. High fusion power generation has been obtained by using NBI heating, which is the main heating mechanism at JET. Other specific experiments have focused on several important topics that were not previously studied in D-T. A particular emphasis has been put on the exploration of some of the unique features expected in future D-T plasmas, i.e. simultaneous dominant electron heating, low rotation, and fully destabilized energetic ions instabilities. Since dominant electron heating by alpha particles is not possible in the JET tokamak at low rotation, the electron heating and high energetic particle generation, able to destabilize magnetic perturbations in stationary conditions, has been externally provided by using the maximum ICRF power available in DTE2.

A major result has been obtained suggesting that some reactor-relevant plasma conditions may be very beneficial. The development of large-scale energetic particle perturbations in the presence of highly energetic ions significantly reduces the conductive-convective energy losses driven by microturbulence in the plasma core leading to a good core energy confinement. The strong impact of zonal flows has been found to play a key role in the reduction of turbulent energy transport in the conditions explored. Importantly, at such low levels of energy losses, there is a clear asymmetry between T and D, since the presence of T significantly enhances the zonal flow activity and further reduces the energy losses by transport, resulting in better global confinement in D-T compared to D.

Such results are obtained in a novel regime developed at JET, similar to the I-mode, consisting of the onset of a pedestal in the temperature while the density remains in L-mode and no damaging ELMs are detected. This new high confinement plasma regime is obtained close to the L-H transition threshold, and therefore requiring low input power. The presence of T is crucial also in the plasma edge, as at equivalent engineering parameters, a pedestal is developed in D-T but not in pure D.

Exploring a broad range of plasma conditions in present-day tokamaks is essential to evaluate how D-T plasmas might behave in the future. This is because the physics mechanisms expected to play a role in ITER and the future fusion reactor cannot be fully reproduced in an integrated way in existing tokamaks, and therefore specific studies must be performed. From this perspective, our studies complement the results obtained in DTE1 as they focus in different plasma regimes. DTE1 provided clear evidence of alpha heating, while differences in energy transport in D-T compared to D were weak and mostly with origin on the pedestal in H-mode plasmas heated with NBI. In DTE2, such results have been reproduced in H-mode plasmas when NBI is used as the main heating system and Ti/Te > 152,59. However, our results significantly expand the knowledge about D-T by exploring some other conditions expected in future D-T plasmas that cannot be obtained in high NBI heated conditions. We show that close to the turbulence threshold in the presence of energetic ions instabilities, as expected in ITER48, zonal flows can play a significant role and provide a route to lower core energy transport in D-T than in D. Importantly, the change from C-wall in DTE1, to Be and W wall and divertor in DTE2 has expanded the operational regime in which plasmas can be developed, significantly broadening the possibilities to perform experiments under conditions, notably at the plasma edge, not developed in DTE144.

Our findings pave the way for a more economical and simpler design of tokamaks, confirming that nuclear fusion by means of magnetically confined D-T plasmas is a promising source of clean energy. However, further studies, such as compatibility with power exhaust capabilities and exploration at higher density and power, are required to fully qualify these plasmas as a solid route toward tokamak reactors. Furthermore, it is necessary to perform more detailed modeling activities to analyze multiscale effects involving energetic ions and plasma perturbations at different spatial and temporal scales, including the non-linear interaction of alpha particles with magnetic perturbations. The most important aspect is to develop plasmas in which, unlike in the JET results shown in this paper, alpha particle heating is dominant. To this end, experimental and modeling efforts in D-T including a significant population of alpha particles, as expected in ITER60 or SPARC61, are essential.

Methods

Experimental design

The JET tokamak has investigated some of the most important fusion reactor conditions by conducting a new D-T campaign with Be/W wall. To reproduce the simultaneous high electron heating, low torque, and the destabilization of energetic ions-related instabilities expected in future tokamak reactors, JET has mostly used ICRF power rather than NBI. Several experimental conditions were explored to find an optimum plasma state in terms of confinement, energetic ion production, and D-T fusion power yield. One of the limitations of such exploratory work was the ICRF power, and hence the amount of electron heating available, which was limited to  ∼4.5 MW, whereas the NBI power was used up to 10 MW in pure NBI plasmas. The ICRF frequency used was 55 MHz at BT= 3.7 T and 42 MHz at BT= 2.75 T.

Different scans were performed for Ip and BT, e.g. Ip was explored in the range 1.9 MA  < Ip < 2.5 MA at two BT = 2.75 T, 3.7 T. The D-T concentration ratio was scanned from 40% T to 85% T by using different valves injecting neutral D and T gases.

Key diagnostics for D-T operation

The ion temperature profiles in this paper were obtained from charge exchange recombination spectroscopy (CXRS)62 measurements of impurity ions and electron temperature profiles from combined analysis of LIDAR Thomson scattering63 and high-resolution Thomson scattering (HRTS)64 diagnostics. The density profiles were taken from HRTS measurements, with the density normalized to match the line-average density measured by a far-infrared interferometer.

Mirnov coils are used as a standard MHD diagnostic on almost all tokamak devices. The coils are installed within the vacuum vessel close to the plasma boundary and provide a measurement of the time derivative of the magnetic field. Magnetic spectrograms (Fourier decomposition of the Mirnov coil signal) can then be used to identify relevant oscillation frequencies associated with MHD activity. In JET a number of coil arrays with high-frequency response are available, allowing activity in the Alfvén range to be observed.

The time-resolved neutron yield in JET is measured using three fission chambers, containing 235U and 238U, located outside the vacuum vessel.

The Alfvén eigenmode active diagnostic (AEAD)65 is characterized by six toroidally spaced antennas, each with independent power and phasing, whose aim is to actively excite marginally stable TAEs.

Alpha particle losses are detected by the fast ion loss detector (FILD) consisting of a Faraday cup array37. The Faraday cup array is composed of multiple cups that span a wide poloidal angle below the outboard midplane at a single toroidal location with a minor radial extent.

ELMs are characterized by the BeII emission signal from the inner divetor region.

The plasma isotopic composition is measured at the divertor comparing the relative amplitude of Balmer Dα and Tα spectral lines. The D and T ratio in the plasma core is assumed to be equal to the edge, as is usually the case in JET in the presence of multi-ion plasmas when turbulence is driven by ITG66.

The JET X-mode reflectometry diagnostic35 is composed of four distinct radial correlation reflectometers. All these reflectometers probe the mid-plane JET plasma. Plasma fluctuations can be obtained from the phase fluctuations of the reflectometer signal.

Magnetic perturbation spatial location and q profile verification

The q profile for discharge #99896 has been obtained by means of a loop between the EFIT code and the TRANSP code. The EFIT code calculates the magnetic equilibrium with the input data for the energetic ions content from TRANSP simulations. After a few iterations, a converged q profile is obtained. The validation of the q profile obtained from TRANSP and used for modeling with FAR3D was carried out against a series of diagnostics and MHD markers. As markers with the strongest signature in the diagnostics, the destabilised NTMs were identified and their toroidal mode number calculated using a toroidal array of Mirnov coils. The radial location of the associated rational surface q = 4/3 was inferred using two methods: the first uses as proxy the location of the phase inversion of the perturbed electron temperature derived from Electron Cyclotron Emission (ECE) at the NTM frequency, and the second matched the NTM frequency as derived from the Mirnov coils to the Doppler-shifted plasma rotation i.e. n × Vphi, where Vphi is the toroidal rotation of the main plasma ions as derived from CXRS diagnostic. The radial location of the q = 1/1 surface was inferred from the inversion radius of the ECE temperature profile during sawtooth crashes as well as from the fishbones signatures in the perturbed plasma temperature from ECE, evidencing a typical kink-like pattern inside the q = 1 surface. Lastly, the RSAEs and TAEs were located using Soft X-ray cameras, interferometry and reflectometry.

As shown in Supplementary Fig. 6, in the Supplementary Information, the agreement between the q profile from TRANSP and the MHD markers is good for discharge #99896. The q profile obtained for the D discharge #100871 is very similar compared to the #99896.

Experimental profiles fitting

The profile fitting algorithm makes use of a Gaussian process regression (GPR), which is not limited by a selection of specific fit functions and provides a statistically rigorous estimation of the confidence bounds of the fit. For more details, see the book on the topic written by Carl Edward Rasmussen and Christopher K. I. Williams http://gaussianprocess.org/gpml/chapters/.

TRANSP simulations

The pulses shown in this article were analyzed by interpretive simulations performed with the TRANSP modelling suite67 coupled with external heating modules NUBEAM (NBI)68 and TORIC (ICRF)69, and prepared with the OMFIT integrated modelling platform70. Interpretive analysis was based on the use of fitted profiles, including electron density and temperatures. The fitting of Te, ne and Ti were performed on data obtained from HRTS and CXRS. The fitting of experimental profiles consists on applying a global third-order polynomial fit in the range ρ ≲ 0.8 (with the additional constraint ∂Ti(0)/∂r = 0).

FAR3D description and simulations parameters

The gyrofluid FAR3D code solves the linear and nonlinear reduced resistive MHD equations describing the thermal plasma evolution coupled with the first two moments of the gyro-kinetic equation, the equations of the energetic particle density and parallel velocity moments71,72, introducing the wave-particle resonance effects required for Landau damping/growth. The correct model calibration requires performing gyrokinetic simulations to calculate the Landau closure coefficients in the gyrofluid simulations, matching the analytic TAE growth rates of the two-pole approximation of the plasma dispersion function with a Lorentzian energy distribution function for the energetic particles. The lowest-order Lorentzian is matched with a Maxwellian distribution by choosing an equivalent average energy. Further details of the model equations can be found in references73,74.

A set of linear simulations is performed to reproduce the Alfven eigenmode (AE) activity observed in the discharge, identifying the resonance induced by populations of energetic particles (EP) as passing D and trapped H. The analysis is based on a parametric study with respect to the EP beta (EP density in the plasma) and energy, calculating AEs consistent with the frequency range, plasma radial location, modes number and AE family observed in the experiment. Nonlinear simulations including passing D and trapped H populations are performed to analyze the saturation phase of the AE instabilities, particularly the energetic particle transport induced, the generation of zonal structure, and the nonlinear interaction between different EP populations. The simulations are performed using the EP model profiles obtained from TRANSP, the measured thermal plasma profiles, and the equilibrium calculated with VMEC code75.

CGYRO description and simulations parameters

The CGYRO code45 solves the electromagnetic gyrokinetic-Maxwell equations46. Local simulations were carried out at ρ = 0.31. Shaped, up-down symmetric flux-surface geometry was used and multi-species collisions were included using the Sugama collision operator76. Transverse and compressional electromagnetic fluctuations were retained. Rotation effects were assumed to be small in the core and not included. Kinetic electrons, D and T as separate species, were included in the simulations. Regarding the energetic-ion species, a lumped H-D species with effective mass averaged between H and D was assumed and modeled by fitting the energetic particle distribution to an equivalent high-temperature Maxwellian equilibrium distribution.

The simulations used a radial box length of Lx = 673ρs and a binormal box of length Ly = 628ρs. Nr = 768 radial modes and Nα = 64 complex toroidal modes were retained. Other resolution parameters were: Nθ = 24 (field line resolution), Nξ = 24 (pitch-angle resolution), Nu = 8 (energy resolution) with maximum energy umax2 = 8. The definitions of the CGYRO numerical resolution parameters can be found in45. The energy flux is provided in GB unit defined as QGB=neTecsρ*2, with cs=Te/mp, ρ* = ρs/a is the ratio of the proton sound gyroradius, ρs = cs/Ωc, to the system size, with Ωc = eB0/mp the ion gyrofrequency. Convergence tests have been performed indicating that the CGYRO results are well resolved.

The zonal flow shearing is defined in CGYRO as:1 ωE×B=∑kxkx2ρs2ϕ^(ky,kx)2tky=0

where kx is the radial wavenumber, ϕ^(ky,kx) is the fluctuating electrostatic potential and 〈〉 denotes the temporal average.

The equilibrium profile and geometry parameters are given in Supplementary Table 2 in the Suplementary Information.

Supplementary information

Supplementary Information

Peer Review File

Supplementary information

The online version contains supplementary material available at 10.1038/s41467-024-52182-z.

Acknowledgements

J. Garcia would like to thank Gerardo Giruzzi for fruitful discussions. This work has been carried out within the framework of the EUROfusion Consortium, funded by the European Union via the Euratom Research and Training Programme (Grant Agreement nos. 101052200 — EUROfusion). Views and opinions expressed are however those of the author(s) only and do not necessarily reflect those of the European Union or the European Commission. Neither the European Union nor the European Commission can be held responsible for them. This work was supported in part by Grants FIS2017-85252-R and PID2021-127727OB-I00 funded by the Spanish Ministry of Science, Innovation and Universities MICIU/AEI/10.13039/501100011033, by ERDF “A way of making Europe” and by ERDF/EU. An award of computer time was provided by the INCITE program and ALCC program. This research used resources from the Oak Ridge Leadership Computing Facility, which is an Office of Science User Facility supported under Contract DE-AC05-00OR22725. Computing resources were also provided by the National Energy Research Scientific Computing Center, which is an Office of Science User Facility supported under Contract DEAC02-05CH11231. This work was partially supported by the project US DOE under grant DE-FG02-04ER54742. D.Z. received financial support from the AIM4EP Project (ANR-21-CE30-0018), funded by the French National Research Agency (ANR).

Author contributions

The reported experiments were devised and jointly led by Y.K., J.O., S.S., J.G., and M.B., with the key coordination of E. de la L., C.F.M and J.M. The TRANSP simulations were performed by Ž.Š. and M.P. Gyrokinetic simulations and subsequent analyses were performed by E.B, J.C. and S.M. FAR3D simulations were performed by J.V. with the assistance of D.Z. Reflectometer analyses were performed by M.D. and J.R.R. MHD analyses were performed by R.C. and M.D. Pedestal analyses were performed by E.de la L. and E.S. Alpha particle losses were investigated by P.J.B. The manuscript was written by J.G. and E. de la L. with feedback by all the authors.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. A peer review file is available.

Data availability

The JET experimental data is stored in the PPF (Processed Pulse File) system which is a centralised data storage and retrieval system for data derived from raw measurements within the JET Torus, and from other sources such as simulation programs. These data are fully available for the EUROfusion consortium members and can be accessed by non-members under request to EUROfusion. Numerical data supporting the outcome of this study are available from the corresponding author upon request.

Code availability

The research codes cited in the paper require a prior detailed knowledge of the implemented physics models and are under continuous development. The corresponding author can be contacted for any further information.

Competing interests

The authors declare no competing interests.

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

A list of authors and their affiliations appears at the end of the paper.

A full list of members and their affiliations appears in the Supplementary Information
==== Refs
References

1. Smirnov VP Tokamak foundation in USSR/Russia 1950-1990 Nucl. Fusion 2009 50 014003 10.1088/0029-5515/50/1/014003
Smirnov, V. P. Tokamak foundation in USSR/Russia 1950-1990. Nucl. Fusion 50, 014003 (2009).10.1088/0029-5515/50/1/014003
2. Ongena J Koch R Wolf R Zohm H Magnetic-confinement fusion Nat. Phys. 2016 12 398 410 10.1038/nphys3745
Ongena, J., Koch, R., Wolf, R. & Zohm, H. Magnetic-confinement fusion. Nat. Phys. 12, 398–410 (2016).10.1038/nphys3745
3. Wagner F Development of an edge transport barrier at the H-mode transition of Asdex Phys. Rev. Lett. 1984 53 1453 1456 10.1103/PhysRevLett.53.1453
Wagner, F. et al. Development of an edge transport barrier at the H-mode transition of Asdex. Phys. Rev. Lett. 53, 1453–1456 (1984).10.1103/PhysRevLett.53.1453
4. Connor JW the ITB Database Group, and the ITPA Topical Group on transport and internal barrier physics Nucl. Fusion 2004 44 R1 R49 10.1088/0029-5515/44/4/R01
Connor, J. W. et al. the ITB Database Group, and the ITPA Topical Group on transport and internal barrier physics. Nucl. Fusion 44, R1–R49 (2004).10.1088/0029-5515/44/4/R01
5. Song Y Realization of thousand-second improved confinement plasma with super I-mode in tokamak EAST Sci. Adv. 2023 9 eabq5273 10.1126/sciadv.abq5273 36608124
Song, Y. et al. Realization of thousand-second improved confinement plasma with super I-mode in tokamak EAST. Sci. Adv. 9, eabq5273 (2023).36608124 10.1126/sciadv.abq5273
6. Han H A sustained high-temperature fusion plasma regime facilitated by fast ions Nature 2022 609 269 275 10.1038/s41586-022-05008-1 36071190
Han, H. et al. A sustained high-temperature fusion plasma regime facilitated by fast ions. Nature 609, 269–275 (2022).36071190 10.1038/s41586-022-05008-1
7. Zohm H Edge localized modes (ELMs) Plasma Phys. Control. Fusion 1996 38 105 128 10.1088/0741-3335/38/2/001
Zohm, H. Edge localized modes (ELMs). Plasma Phys. Control. Fusion 38, 105–128 (1996).10.1088/0741-3335/38/2/001
8. Strachan JD TFTR DT experiments Plasma Phys. Control. Fusion 1997 39 B103 10.1088/0741-3335/39/12B/008
Strachan, J. D. et al. TFTR DT experiments. Plasma Phys. Control. Fusion 39, B103 (1997).10.1088/0741-3335/39/12B/008
9. Jacquinot J Overview of ITER physics deuterium-tritium experiments in JET Nucl. Fusion 1999 39 235 10.1088/0029-5515/39/2/307
Jacquinot, J. et al. Overview of ITER physics deuterium-tritium experiments in JET. Nucl. Fusion 39, 235 (1999).10.1088/0029-5515/39/2/307
10. Beurskens MNA Confinement in electron heated plasmas in Wendelstein 7-X and Asdex Upgrade; the necessity to control turbulent transport Nucl. Fusion 2021 62 016015 10.1088/1741-4326/ac36f1
Beurskens, M. N. A. et al. Confinement in electron heated plasmas in Wendelstein 7-X and Asdex Upgrade; the necessity to control turbulent transport. Nucl. Fusion 62, 016015 (2021).10.1088/1741-4326/ac36f1
11. Hahm TS Burrell KH Flow shear induced fluctuation suppression in finite aspect ratio shaped tokamak plasma Phys. Plasmas 1995 2 1648 1651 10.1063/1.871313
Hahm, T. S. & Burrell, K. H. Flow shear induced fluctuation suppression in finite aspect ratio shaped tokamak plasma. Phys. Plasmas 2, 1648–1651 (1995).10.1063/1.871313
12. Doyle EJ Chapter 2: Plasma confinement and transport Nucl. Fusion 2007 47 S18 10.1088/0029-5515/47/6/S02
Doyle, E. J. et al. Chapter 2: Plasma confinement and transport. Nucl. Fusion 47, S18 (2007).10.1088/0029-5515/47/6/S02
13. Fasoli A Chapter 5: Physics of energetic ions Nucl. Fusion 2007 47 S264 10.1088/0029-5515/47/6/S05
Fasoli, A. et al. Chapter 5: Physics of energetic ions. Nucl. Fusion 47, S264 (2007).10.1088/0029-5515/47/6/S05
14. Gorelenkov N Pinches SD Toi K Energetic particle physics in fusion research in preparation for burning plasma experiments Nucl. Fusion 2014 54 125001 10.1088/0029-5515/54/12/125001
Gorelenkov, N., Pinches, S. D. & Toi, K. Energetic particle physics in fusion research in preparation for burning plasma experiments. Nucl. Fusion 54, 125001 (2014).10.1088/0029-5515/54/12/125001
15. Todo Y Introduction to the interaction between energetic particles and Alfvén eigenmodes in toroidal plasmas Rev. Mod. Plasma Phys. 2019 3 1 10.1007/s41614-018-0022-9
Todo, Y. Introduction to the interaction between energetic particles and Alfvén eigenmodes in toroidal plasmas. Rev. Mod. Plasma Phys. 3, 1 (2019).10.1007/s41614-018-0022-9
16. Heidbrink WW White RB Mechanisms of energetic-particle transport in magnetically confined plasmas Phys. Plasmas 2020 27 030901 10.1063/1.5136237
Heidbrink, W. W. & White, R. B. Mechanisms of energetic-particle transport in magnetically confined plasmas. Phys. Plasmas 27, 030901 (2020).10.1063/1.5136237
17. Mailloux J Overview of JET results for optimising ITER operation Nucl. Fusion 2022 62 042026 10.1088/1741-4326/ac47b4
Mailloux, J. et al. Overview of JET results for optimising ITER operation. Nucl. Fusion 62, 042026 (2022).10.1088/1741-4326/ac47b4
18. Maggi, C. F. et al. Overview of T and D-T results in jet with ITER-like wall. Nuclear Fusion 64, 112012 (2024).
19. Maslov, M. et al. Jet d-t scenario with optimized non-thermal fusion. Nuclear Fusion 63, 112002 (2023).
20. ITER Physics Expert Group on Confin Transport, ITER Physics Expert Group on Confin Database, and ITER Physics Basis Editors. Chapter 2: Plasma confinement and transport. Nuclear Fusion 39(12): 2175–2249 dec 1999. 10.1088/0029-5515/39/12/302
21. Polevoi AR PFPO plasma scenarios for exploration of long pulse operation in ITER Nucl. Fusion 2023 63 076003 10.1088/1741-4326/acd06f
Polevoi, A. R. et al. PFPO plasma scenarios for exploration of long pulse operation in ITER. Nucl. Fusion 63, 076003 (2023).10.1088/1741-4326/acd06f
22. Hawryluk, R. J. An empirical approach to tokamak transport. In Physics of plasmas close to thermonuclear conditions, pages 19–46. Elsevier, 1981.
23. Angioni C Tungsten transport in JET H-mode plasmas in hybrid scenario, experimental observations and modelling Nucl. Fusion 2014 54 083028 10.1088/0029-5515/54/8/083028
Angioni, C. et al. Tungsten transport in JET H-mode plasmas in hybrid scenario, experimental observations and modelling. Nucl. Fusion 54, 083028 (2014).10.1088/0029-5515/54/8/083028
24. Cheng CZ Chen L Chance MS High-n ideal and resistive shear Alfvén waves in tokamaks Ann. Phys. 1985 161 21 47 10.1016/0003-4916(85)90335-5
Cheng, C. Z., Chen, L. & Chance, M. S. High-n ideal and resistive shear Alfvén waves in tokamaks. Ann. Phys. 161, 21–47 (1985).10.1016/0003-4916(85)90335-5
25. Cheng CZ Chance MS Low-n shear Alfvén spectra in axisymmetric toroidal plasmas Phys. Fluids 1986 29 3695 10.1063/1.865801
Cheng, C. Z. & Chance, M. S. Low-n shear Alfvén spectra in axisymmetric toroidal plasmas. Phys. Fluids 29, 3695 (1986).10.1063/1.865801
26. Sharapov SE Alfvén wave cascades in a tokamak Phys. Plasmas 2002 9 2027 2036 10.1063/1.1448346
Sharapov, S. E. et al. Alfvén wave cascades in a tokamak. Phys. Plasmas 9, 2027–2036 (2002).10.1063/1.1448346
27. Chen L White RB Rosenbluth MN Excitation of internal kink modes by trapped energetic beam ions Phys. Rev. Lett. 1984 52 1122 1125 10.1103/PhysRevLett.52.1122
Chen, L., White, R. B. & Rosenbluth, M. N. Excitation of internal kink modes by trapped energetic beam ions. Phys. Rev. Lett. 52, 1122–1125 (1984).10.1103/PhysRevLett.52.1122
28. Coppi B Porcelli F Theoretical model of fishbone oscillations in magnetically confined plasmas Phys. Rev. Lett. 1986 57 2272 2275 10.1103/PhysRevLett.57.2272 10033680
Coppi, B. & Porcelli, F. Theoretical model of fishbone oscillations in magnetically confined plasmas. Phys. Rev. Lett. 57, 2272–2275 (1986).10033680 10.1103/PhysRevLett.57.2272
29. Chang Z Observation of nonlinear neoclassical pressure-gradient–driven tearing modes in tftr Phys. Rev. Lett. 1995 74 4663 4666 10.1103/PhysRevLett.74.4663 10058567
Chang, Z. et al. Observation of nonlinear neoclassical pressure-gradient–driven tearing modes in tftr. Phys. Rev. Lett. 74, 4663–4666 (1995).10058567 10.1103/PhysRevLett.74.4663
30. Sauter O Beta limits in long-pulse tokamak discharges Phys. Plasmas 1997 4 1654 1664 10.1063/1.872270
Sauter, O. et al. Beta limits in long-pulse tokamak discharges. Phys. Plasmas 4, 1654–1664 (1997).10.1063/1.872270
31. Romanelli F Ion temperature-gradient-driven modes and anomalous ion transport in tokamaks Phys. Fluids B: Plasma Phys. 1989 1 1018 1025 10.1063/1.859023
Romanelli, F. Ion temperature-gradient-driven modes and anomalous ion transport in tokamaks. Phys. Fluids B: Plasma Phys. 1, 1018–1025 (1989).10.1063/1.859023
32. Mazzi S Enhanced performance in fusion plasmas through turbulence suppression by megaelectronvolt ions Nat. Phys. 2022 18 776 782 10.1038/s41567-022-01626-8
Mazzi, S. et al. Enhanced performance in fusion plasmas through turbulence suppression by megaelectronvolt ions. Nat. Phys. 18, 776–782 (2022).10.1038/s41567-022-01626-8
33. Mantica P The role of electron-scale turbulence in the jet tokamak: Experiments and modelling Nucl. Fusion 2021 61 096014 10.1088/1741-4326/ac146e
Mantica, P. et al. The role of electron-scale turbulence in the jet tokamak: Experiments and modelling. Nucl. Fusion 61, 096014 (2021).10.1088/1741-4326/ac146e
34. Elgar S Guza RT Statistics of bicoherence IEEE Trans. Acoust. Speech Signal Process. 1988 36 1667 1668 10.1109/29.7555
Elgar, S. & Guza, R. T. Statistics of bicoherence. IEEE Trans. Acoust. Speech Signal Process. 36, 1667–1668 (1988).10.1109/29.7555
35. Hacquin S Localized X-mode reflectometry measurements of Alfvén eigenmodes on the JET tokamak Plasma Phys. Control. Fusion 2007 49 1371 10.1088/0741-3335/49/9/002
Hacquin, S. et al. Localized X-mode reflectometry measurements of Alfvén eigenmodes on the JET tokamak. Plasma Phys. Control. Fusion 49, 1371 (2007).10.1088/0741-3335/49/9/002
36. Spong DA Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-landau closure model Nucl. Fusion 2021 61 116061 10.1088/1741-4326/ac2990
Spong, D. A. et al. Nonlinear dynamics and transport driven by energetic particle instabilities using a gyro-landau closure model. Nucl. Fusion 61, 116061 (2021).10.1088/1741-4326/ac2990
37. Bonofiglo PJ Improvements to the Faraday cup fast ion loss detector and magnetohydrodynamic induced fast ion loss measurements in Joint European Torus plasmas Rev. Sci. Instrum. 2020 91 093502 10.1063/5.0014278 33003824
Bonofiglo, P. J. et al. Improvements to the Faraday cup fast ion loss detector and magnetohydrodynamic induced fast ion loss measurements in Joint European Torus plasmas. Rev. Sci. Instrum. 91, 093502 (2020).33003824 10.1063/5.0014278
38. Chen L Zonca F Nonlinear excitations of zonal structures by toroidal Alfvén eigenmodes Phys. Rev. Lett. 2012 109 145002 10.1103/PhysRevLett.109.145002 23083251
Chen, L. & Zonca, F. Nonlinear excitations of zonal structures by toroidal Alfvén eigenmodes. Phys. Rev. Lett. 109, 145002 (2012).23083251 10.1103/PhysRevLett.109.145002
39. Mishchenko A Numerical tools for burning plasmas Plasma Phys. Control. Fusion 2023 65 064001 10.1088/1361-6587/acce68
Mishchenko, A. et al. Numerical tools for burning plasmas. Plasma Phys. Control. Fusion 65, 064001 (2023).10.1088/1361-6587/acce68
40. Di Siena A Görler T Doerk H Poli E Bilato R Fast-ion stabilization of tokamak plasma turbulence Nucl. Fusion 2018 58 054002 10.1088/1741-4326/aaaf26
Di Siena, A., Görler, T., Doerk, H., Poli, E. & Bilato, R. Fast-ion stabilization of tokamak plasma turbulence. Nucl. Fusion 58, 054002 (2018).10.1088/1741-4326/aaaf26
41. Diamond PH Itoh SI Itoh K Hahm TS Zonal flows in plasmas - a review Plasma Phys. Control. Fusion 2005 47 R35 10.1088/0741-3335/47/5/R01
Diamond, P. H., Itoh, S. I., Itoh, K. & Hahm, T. S. Zonal flows in plasmas - a review. Plasma Phys. Control. Fusion 47, R35 (2005).10.1088/0741-3335/47/5/R01
42. Garcia J JET Contributors. Electromagnetic and fast ions effects as a key mechanism for turbulent transport suppression at JET Plasma Phys. Control. Fusion 2022 64 104002 10.1088/1361-6587/ac8613
Garcia, J. & JET Contributors. Electromagnetic and fast ions effects as a key mechanism for turbulent transport suppression at JET. Plasma Phys. Control. Fusion 64, 104002 (2022).10.1088/1361-6587/ac8613
43. Garcia J Görler T Jenko F Giruzzi G Gyrokinetic nonlinear isotope effects in tokamak plasmas Nucl. Fusion 2016 57 014007 10.1088/1741-4326/57/1/014007
Garcia, J., Görler, T., Jenko, F. & Giruzzi, G. Gyrokinetic nonlinear isotope effects in tokamak plasmas. Nucl. Fusion 57, 014007 (2016).10.1088/1741-4326/57/1/014007
44. Garcia J New H-mode regimes with small ELMs and high thermal confinement in the Joint European Torus Phys. Plasmas 2022 29 032505 10.1063/5.0072236
Garcia, J. et al. New H-mode regimes with small ELMs and high thermal confinement in the Joint European Torus. Phys. Plasmas 29, 032505 (2022).10.1063/5.0072236
45. Candy J Belli EA Bravenec RV A high-accuracy eulerian gyrokinetic solver for collisional plasmas J. Comput. Phys. 2016 324 73 93 10.1016/j.jcp.2016.07.039
Candy, J., Belli, E. A. & Bravenec, R. V. A high-accuracy eulerian gyrokinetic solver for collisional plasmas. J. Comput. Phys. 324, 73–93 (2016).10.1016/j.jcp.2016.07.039
46. Sugama H Horton W Nonlinear electromagnetic gyrokinetic equation for plasmas with large mean flows Phys. Plasmas 1998 5 2560 2573 10.1063/1.872941
Sugama, H. & Horton, W. Nonlinear electromagnetic gyrokinetic equation for plasmas with large mean flows. Phys. Plasmas 5, 2560–2573 (1998).10.1063/1.872941
47. Kim D Turbulence stabilization in tokamak plasmas with high population of fast ions Nucl. Fusion 2023 63 124001 10.1088/1741-4326/acffda
Kim, D. et al. Turbulence stabilization in tokamak plasmas with high population of fast ions. Nucl. Fusion 63, 124001 (2023).10.1088/1741-4326/acffda
48. Garcia J Görler T Jenko F Isotope and fast ions turbulence suppression effects: Consequences for high-β ITER plasmas Phys. Plasmas 2018 25 055902 10.1063/1.5016331
Garcia, J., Görler, T. & Jenko, F. Isotope and fast ions turbulence suppression effects: Consequences for high-β ITER plasmas. Phys. Plasmas 25, 055902 (2018).10.1063/1.5016331
49. Belli EA Candy J Waltz RE Reversal of simple hydrogenic isotope scaling laws in tokamak edge turbulence Phys. Rev. Lett. 2020 125 015001 10.1103/PhysRevLett.125.015001 32678657
Belli, E. A., Candy, J. & Waltz, R. E. Reversal of simple hydrogenic isotope scaling laws in tokamak edge turbulence. Phys. Rev. Lett. 125, 015001 (2020).32678657 10.1103/PhysRevLett.125.015001
50. Maggi CF Isotope effects on L-H threshold and confinement in tokamak plasmas Plasma Phys. Control. Fusion 2017 60 014045 10.1088/1361-6587/aa9901
Maggi, C. F. et al. Isotope effects on L-H threshold and confinement in tokamak plasmas. Plasma Phys. Control. Fusion 60, 014045 (2017).10.1088/1361-6587/aa9901
51. Frassinetti L Effect of the isotope mass on pedestal structure, transport and stability in D, D/T and t plasmas at similar βn and gas rate in JET-ILW type I ELMy H-modes Nucl. Fusion 2023 63 112009 10.1088/1741-4326/acf057
Frassinetti, L. et al. Effect of the isotope mass on pedestal structure, transport and stability in D, D/T and t plasmas at similar βn and gas rate in JET-ILW type I ELMy H-modes. Nucl. Fusion 63, 112009 (2023).10.1088/1741-4326/acf057
52. Schneider PA Isotope physics of heat and particle transport with tritium in JET-ILW type-I ELMy H-mode plasmas Nucl. Fusion 2023 63 112010 10.1088/1741-4326/acf560
Schneider, P. A. et al. Isotope physics of heat and particle transport with tritium in JET-ILW type-I ELMy H-mode plasmas. Nucl. Fusion 63, 112010 (2023).10.1088/1741-4326/acf560
53. Horvath L Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas Nucl. Fusion 2021 61 046015 10.1088/1741-4326/abdd77
Horvath, L. et al. Isotope dependence of the type I ELMy H-mode pedestal in JET-ILW hydrogen and deuterium plasmas. Nucl. Fusion 61, 046015 (2021).10.1088/1741-4326/abdd77
54. Solano ER Axisymmetric oscillations at l-h transitions in JET: M-mode Nucl. Fusion 2016 57 022021 10.1088/0029-5515/57/2/022021
Solano, E. R. et al. Axisymmetric oscillations at l-h transitions in JET: M-mode. Nucl. Fusion 57, 022021 (2016).10.1088/0029-5515/57/2/022021
55. Delabie, E. et al. The low density type III ELMy H-mode regime on JET-ILW: a low density H-mode compatible with a tungsten divertor? BP10.00052, 58th Annual Meeting of the APS Division of Plasma Physics, October 31-November 4 2016, (San Jose, California, 2016).
56. Ryter F H-mode power threshold and transition in Asdex Upgrade Plasma Phys. Control. Fusion 1998 40 725 10.1088/0741-3335/40/5/032
Ryter, F. et al. H-mode power threshold and transition in Asdex Upgrade. Plasma Phys. Control. Fusion 40, 725 (1998).10.1088/0741-3335/40/5/032
57. Greenwald M Transport phenomena in Alcator C-mod H-modes Plasma Phys. Control. Fusion 1998 40 789 10.1088/0741-3335/40/5/044
Greenwald, M. et al. Transport phenomena in Alcator C-mod H-modes. Plasma Phys. Control. Fusion 40, 789 (1998).10.1088/0741-3335/40/5/044
58. Whyte DG I-mode: an h-mode energy confinement regime with l-mode particle transport in Alcator C-mod Nucl. Fusion 2010 50 105005 10.1088/0029-5515/50/10/105005
Whyte, D. G. et al. I-mode: an h-mode energy confinement regime with l-mode particle transport in Alcator C-mod. Nucl. Fusion 50, 105005 (2010).10.1088/0029-5515/50/10/105005
59. Kiptily VG Evidence of electron heating by alpha particles in JET deuterium-tritium plasmas Phys. Rev. Lett. 2023 131 075101 10.1103/PhysRevLett.131.075101 37656860
Kiptily, V. G. et al. Evidence of electron heating by alpha particles in JET deuterium-tritium plasmas. Phys. Rev. Lett. 131, 075101 (2023).37656860 10.1103/PhysRevLett.131.075101
60. Shimada M Chapter 1: Overview and summary Nucl. Fusion 2007 47 S1 S17 10.1088/0029-5515/47/6/S01
Shimada, M. et al. Chapter 1: Overview and summary. Nucl. Fusion 47, S1–S17 (2007).10.1088/0029-5515/47/6/S01
61. Creely AJ Overview of the SPARC tokamak J. Plasma Phys. 2020 86 865860502 10.1017/S0022377820001257
Creely, A. J. et al. Overview of the SPARC tokamak. J. Plasma Phys. 86, 865860502 (2020).10.1017/S0022377820001257
62. Hawkes NC Instrumentation for the upgrade to the JET core charge-exchange spectrometers Rev. Sci. Instrum. 2018 89 10D113 10.1063/1.5037639 30399852
Hawkes, N. C. et al. Instrumentation for the upgrade to the JET core charge-exchange spectrometers. Rev. Sci. Instrum. 89, 10D113 (2018).30399852 10.1063/1.5037639
63. Maslov M Beurskens MNA Kempenaars M Flanagan J Status of the jet lidar Thomson scattering diagnostic J. Instrum. 2013 8 C11009 10.1088/1748-0221/8/11/C11009
Maslov, M., Beurskens, M. N. A., Kempenaars, M. & Flanagan, J. Status of the jet lidar Thomson scattering diagnostic. J. Instrum. 8, C11009 (2013).10.1088/1748-0221/8/11/C11009
64. Pasqualotto R High resolution Thomson scattering for Joint European Torus (JET) Rev. Sci. Instrum. 2004 75 3891 3893 10.1063/1.1787922
Pasqualotto, R. et al. High resolution Thomson scattering for Joint European Torus (JET). Rev. Sci. Instrum. 75, 3891–3893 (2004).10.1063/1.1787922
65. Tinguely RA Simultaneous measurements of unstable and stable Alfvén eigenmodes in JET Nucl. Fusion 2022 62 112008 10.1088/1741-4326/ac899e
Tinguely, R. A. et al. Simultaneous measurements of unstable and stable Alfvén eigenmodes in JET. Nucl. Fusion 62, 112008 (2022).10.1088/1741-4326/ac899e
66. Maslov M Observation of enhanced ion particle transport in mixed H/D isotope plasmas on JET Nucl. Fusion 2018 58 076022 10.1088/1741-4326/aac342
Maslov, M. et al. Observation of enhanced ion particle transport in mixed H/D isotope plasmas on JET. Nucl. Fusion 58, 076022 (2018).10.1088/1741-4326/aac342
67. Ongena J Voitsekhovitch I Evrard M McCune D Numerical transport codes Fusion Sci. Technol. 2012 61 180 189 10.13182/FST12-A13505
Ongena, J., Voitsekhovitch, I., Evrard, M. & McCune, D. Numerical transport codes. Fusion Sci. Technol. 61, 180–189 (2012).10.13182/FST12-A13505
68. Pankin A McCune D Andre R Bateman G Kritz A The tokamak Monte Carlo fast ion module NUBEAM in the National Transport Code Collaboration library Comput. Phys. Commun. 2004 159 157 184 10.1016/j.cpc.2003.11.002
Pankin, A., McCune, D., Andre, R., Bateman, G. & Kritz, A. The tokamak Monte Carlo fast ion module NUBEAM in the National Transport Code Collaboration library. Comput. Phys. Commun. 159, 157–184 (2004).10.1016/j.cpc.2003.11.002
69. Brambilla M Numerical simulation of ion cyclotron waves in tokamak plasmas Plasma Phys. Control. Fusion 1999 41 1 10.1088/0741-3335/41/1/002
Brambilla, M. Numerical simulation of ion cyclotron waves in tokamak plasmas. Plasma Phys. Control. Fusion 41, 1 (1999).10.1088/0741-3335/41/1/002
70. Grierson BA Orchestrating TRANSP simulations for interpretative and predictive tokamak modeling with OMFIT Fusion Sci. Technol. 2018 74 101 115 10.1080/15361055.2017.1398585
Grierson, B. A. et al. Orchestrating TRANSP simulations for interpretative and predictive tokamak modeling with OMFIT. Fusion Sci. Technol. 74, 101–115 (2018).10.1080/15361055.2017.1398585
71. Hedrick CL Leboeuf J-N Spong DA Alpha-Alfvén local dispersion relation and solutions Phys. Fluids B: Plasma Phys. 1992 4 3869 3882 10.1063/1.860344
Hedrick, C. L., Leboeuf, J.-N. & Spong, D. A. Alpha-Alfvén local dispersion relation and solutions. Phys. Fluids B: Plasma Phys. 4, 3869–3882 (1992).10.1063/1.860344
72. Spong DA Carreras BA Hedrick CL Linearized gyrofluid model of the alpha-destabilized toroidal Alfvén eigenmode with continuum damping effects Phys. Fluids B: Plasma Phys. 1992 4 3316 3328 10.1063/1.860386
Spong, D. A., Carreras, B. A. & Hedrick, C. L. Linearized gyrofluid model of the alpha-destabilized toroidal Alfvén eigenmode with continuum damping effects. Phys. Fluids B: Plasma Phys. 4, 3316–3328 (1992).10.1063/1.860386
73. Varela J Shimizu A Spong DA Garcia L Ghai Y Study of the Alfvén eigenmodes stability in CFQS plasma using a landau closure model Nucl. Fusion 2021 61 026023 10.1088/1741-4326/abd072
Varela, J., Shimizu, A., Spong, D. A., Garcia, L. & Ghai, Y. Study of the Alfvén eigenmodes stability in CFQS plasma using a landau closure model. Nucl. Fusion 61, 026023 (2021).10.1088/1741-4326/abd072
74. Varela J Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode Nucl. Fusion 2021 61 126016 10.1088/1741-4326/ac26a0
Varela, J. et al. Theoretical analysis of the saturation phase of the 1/1 energetic-ion-driven resistive interchange mode. Nucl. Fusion 61, 126016 (2021).10.1088/1741-4326/ac26a0
75. Hirshman SP Whitson JC Steepest-descent moment method for three-dimensional magnetohydrodynamic equilibria Phys. Fluids 1983 26 3553 3568 10.1063/1.864116
Hirshman, S. P. & Whitson, J. C. Steepest-descent moment method for three-dimensional magnetohydrodynamic equilibria. Phys. Fluids 26, 3553–3568 (1983).10.1063/1.864116
76. Sugama H Watanabe T-H Nunami M Linearized model collision operators for multiple ion species plasmas and gyrokinetic entropy balance equations Phys. Plasmas 2009 16 112503 10.1063/1.3257907
Sugama, H., Watanabe, T.-H. & Nunami, M. Linearized model collision operators for multiple ion species plasmas and gyrokinetic entropy balance equations. Phys. Plasmas 16, 112503 (2009).10.1063/1.3257907
