pith. sign in

arxiv: 2604.15024 · v1 · submitted 2026-04-16 · ⚛️ physics.plasm-ph

On nonlinear saturation of toroidal Alfv\'en eigenmode due to thermal plasma nonlinearities

Ningfei Chen , Thomas Hayward-Schneider , Fulvio Zonca , Zhiyong Qiu , Zhixin Lu , Xin Wang , Alessandro Biancalani , Alexey Mishchenko
show 2 more authors
Alberto Bottino Philipp Lauber
This is my paper

Pith reviewed 2026-05-10 09:23 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords toroidal Alfvén eigenmodenonlinear saturationthermal plasma nonlinearitiesgyrokinetic simulationphase-space zonal structureAlfvén continuumzonal fields
0
0 comments X

The pith

Thermal plasma nonlinearities govern toroidal Alfvén eigenmode saturation with stiffness above a 0.47 percent drive threshold.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper uses gyrokinetic particle-in-cell simulations to examine how toroidal Alfvén eigenmodes reach nonlinear saturation. It shows that thermal plasma effects alone can set the saturation level when the linear growth rate is sufficiently large relative to the mode frequency. The saturation occurs through a frequency downshift caused by phase-space structures in the thermal particles, leading the mode to merge with the continuum. Including zonal fields roughly doubles the saturation amplitude. Understanding these processes matters for predicting how energetic particles interact with waves in fusion devices.

Core claim

In single toroidal mode number simulations with zonal fields filtered out, the saturation level of TAE is governed by thermal plasma nonlinearities for gamma_L/omega_n > 0.47%, exhibiting weak dependence on the linear drive gamma_L. The TAE frequency decreases with increasing amplitude due to phase-space zonal structures of thermal plasmas, and saturation is reached when the mode merges into the continuum. Zonal fields counteract the PSZS effects, enhancing the saturation level by a factor of about two.

What carries the argument

Phase-space zonal structures (PSZS) formed by thermal particles, which induce a nonlinear frequency downshift in the TAE until it merges with the Alfvén continuum.

If this is right

  • The saturation amplitude shows little variation with changes in linear drive strength once above the threshold.
  • The TAE frequency downshifts progressively with growing amplitude.
  • At saturation, the mode transitions by merging into the continuum, accompanied by separation of poloidal harmonics.
  • Simulations including zonal fields yield approximately twice the saturation level compared to filtered cases.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • In full multi-mode simulations without filtering, the effective saturation might balance thermal nonlinearities against zonal field effects differently.
  • This stiffness could allow simpler estimates of TAE amplitudes in varying plasma conditions without detailed drive calculations.

Load-bearing premise

Filtering out zonal fields in single-toroidal-mode-number simulations isolates the thermal-plasma nonlinearity without altering the saturation physics or phase-space structures.

What would settle it

A simulation or experiment showing strong dependence of saturation level on linear drive even above 0.47% threshold, or no frequency downshift correlating with amplitude.

Figures

Figures reproduced from arXiv: 2604.15024 by Alberto Bottino, Alessandro Biancalani, Alexey Mishchenko, Fulvio Zonca, Ningfei Chen, Philipp Lauber, Thomas Hayward-Schneider, Xin Wang, Zhixin Lu, Zhiyong Qiu.

Figure 1
Figure 1. Figure 1: (a) The Alfvén continuum used in this paper, [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: The dependence of the (a) growth rate, and [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: The dependence of the (a) growth rate, and [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: The temporal evolution of the maximum of [PITH_FULL_IMAGE:figures/full_fig_p005_5.png] view at source ↗
Figure 7
Figure 7. Figure 7: The ϕ (s, ω) diagram of n = 6 TAE at ωA0t = 700 with only EP nonlinearity for nEP/ne = 0.002 and TEP = 400keV. The black solid curves represent the n = 6 Alfvén continuum shown in [PITH_FULL_IMAGE:figures/full_fig_p005_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: The temporal evolution of the maximum of [PITH_FULL_IMAGE:figures/full_fig_p007_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: The temporal evolution of the frequency of [PITH_FULL_IMAGE:figures/full_fig_p007_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: The ϕ (s, ω) diagram of n = 6 TAE at (a) ωA0t = 200, and (b) ωA0t = 400 for the simulation with all nonlinearities. The black solid curves represent the n = 6 Alfvén continuum shown in [PITH_FULL_IMAGE:figures/full_fig_p008_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: The temporal evolution of the ratio between [PITH_FULL_IMAGE:figures/full_fig_p008_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The temporal evolution of the maximum of [PITH_FULL_IMAGE:figures/full_fig_p008_12.png] view at source ↗
Figure 15
Figure 15. Figure 15: The temporal evolution of the maximum of [PITH_FULL_IMAGE:figures/full_fig_p009_15.png] view at source ↗
Figure 18
Figure 18. Figure 18: The radial mode structure of (a) radial electric field of zonal fields, i.e., ∂sδϕZ, and (b) radial derivative of zonal vector potential for nEP/ne = 0.002, and TEP = 400keV. δϕT Z and δϕE Z denote the theoretical zonal scalar potential generated by thermal plasma nonlinearity given by Equation (4) and resonant EPs contribution given by Equation (31) in Ref. [28], respectively. The subscripts “sim” and “a… view at source ↗
Figure 19
Figure 19. Figure 19: The radial mode structure of (a) radial electric field of zonal fields, i.e., ∂sδϕZ, and (b) modified radial electric field normalized by maximum of TAE intensity, i.e., ∂sδϕZ/ |δϕn| 2 for four values of Lni. dicular nonlinearity and quasi-neutrality condition [29]. The nonlinear nonadiabatic particle responses to zonal fields, i.e., PSZS, for thermal plasmas in the drift-kinetic limit can be derived as δ… view at source ↗
Figure 20
Figure 20. Figure 20: The dependence of saturation level of TAE [PITH_FULL_IMAGE:figures/full_fig_p012_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: The temporal evolution of eδϕn/Te for different combinations of time step and mass ratio. Here, the EP concentration nEP/ne = 0.002. Acknowledgement The authors acknowledge Professor Liu Chen (Zhe￾jiang University, PRC) for inspiration on physical under￾standing and theoretical analysis. This work has been carried out within the framework of the EUROfusion Consortium, funded by the European Union via the … view at source ↗
Figure 22
Figure 22. Figure 22: The radial mode structure of ∂sδϕZ for different electron marker numbers. Here, the mass ratio is mi/me = 1836, time step is ωci∆t = 12, and the EP concentration is nEP/ne = 0.002. with the electron marker number being 2 × 108 . Above all, the default settings for the cases with zonal fields are 2 × 107 , 2 × 108 , 4 × 107 for thermal ions, electrons, EPs, mi/me = 1836, ωci∆t = 12 for nEP/ne = 0.002. The … view at source ↗
read the original abstract

The nonlinear saturation of toroidal Alfven eigenmode (TAE) due to thermal plasma nonlinearities is investigated using gyrokinetic particle-in-cell simulations and theoretical analysis. In the single toroidal mode number simulations with zonal fields filtered out, we find that the saturation level of TAE is governed by thermal plasma nonlinearities for gamma_L/omega_n > 0.47%, which has weak dependence on the linear drive gamma_L, i.e., "stiffness" in saturation level. We find that the frequency of TAE decreases as the amplitude of it increases, which is induced by the phase-space zonal structure (PSZS) of thermal plasmas universally existed in particle-in-cell simulations. The saturation of TAE can be finally reached when the mode merges into the continuum. Following this process, the separation of neighboring poloidal harmonics and mode transition to energetic particle modes can be observed. In simulations with zonal fields, zonal fields can essentially counteract the effects of PSZS of thermal plasmas, leading to roughly a factor of 2 enhancement of the TAE saturation level compared to the single toroidal mode number simulation, implying the necessity of including zonal modes in evaluating the saturation level of TAE.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

3 major / 3 minor

Summary. The manuscript investigates the nonlinear saturation of toroidal Alfvén eigenmodes (TAEs) due to thermal plasma nonlinearities via gyrokinetic particle-in-cell simulations and theoretical analysis. In single-toroidal-mode-number runs with zonal fields filtered out, saturation is reported to be governed by thermal nonlinearities for gamma_L/omega_n > 0.47%, exhibiting stiffness (weak dependence on linear drive gamma_L). The frequency downshift is attributed to phase-space zonal structures (PSZS) of thermal ions, leading to continuum merging, poloidal harmonic separation, and transition to energetic-particle modes. Restoring zonal fields counteracts PSZS effects and increases saturation amplitude by a factor of approximately 2, implying the importance of zonal modes for accurate saturation estimates.

Significance. If the central claims hold, the work would advance understanding of TAE saturation in fusion-relevant plasmas by isolating thermal-plasma nonlinearities and PSZS as key mechanisms, with the stiffness result potentially simplifying amplitude predictions above threshold. The explicit contrast between filtered and unfiltered simulations usefully illustrates competing nonlinear channels. The combination of PIC evidence with interpretive analysis of frequency shifts and mode transitions is a strength, though the absence of analytic derivations or multi-mode benchmarks limits immediate generality.

major comments (3)
  1. [Abstract and simulation results] Abstract and simulation results: The claim that saturation is governed by thermal plasma nonlinearities with stiffness for gamma_L/omega_n > 0.47% rests exclusively on single-n simulations in which zonal fields are filtered. The manuscript states that restoring zonal fields counteracts PSZS and raises saturation amplitude by a factor of ~2; this indicates that the filtered setup omits a dominant saturation channel (possibly via E×B shearing or modified continuum damping), so the reported threshold and stiffness may be artifacts of the numerical isolation rather than intrinsic to thermal nonlinearities. Self-consistent multi-mode benchmarks are needed to confirm the result survives when zonal dynamics are retained.
  2. [Numerical setup and results sections] Numerical setup and results sections: The specific 0.47% threshold and stiffness are extracted from simulation time series without reported convergence tests (particle number, spatial resolution, time step), error bars on saturation levels, or systematic scans in gamma_L. This makes the quantitative claims provisional; small numerical variations could shift the apparent threshold or eliminate the weak gamma_L dependence.
  3. [Discussion of PSZS and frequency shift] Discussion of PSZS and frequency shift: The frequency decrease with amplitude and saturation via continuum merging are presented as induced by PSZS, but the link is observational from time series rather than derived from an equation relating PSZS amplitude to the observed shift or merging criterion. A quantitative model would be required to elevate this from correlation to causal mechanism supporting the saturation picture.
minor comments (3)
  1. Ensure the normalized drive strength gamma_L/omega_n is defined at first use in the main text, and clarify whether omega_n refers to a local or global frequency.
  2. Figure captions and text should explicitly state which runs include zonal-field filtering to avoid ambiguity when comparing saturation levels.
  3. Add references to prior gyrokinetic literature on phase-space zonal structures to place the PSZS observations in context.

Simulated Author's Rebuttal

3 responses · 1 unresolved

We thank the referee for the constructive comments, which help clarify the scope of our findings. We address each major comment point by point below. Revisions have been made to add caveats, numerical details, and expanded discussion, while acknowledging limitations in the current study.

read point-by-point responses

  1. Referee: The claim that saturation is governed by thermal plasma nonlinearities with stiffness for gamma_L/omega_n > 0.47% rests exclusively on single-n simulations in which zonal fields are filtered. ... this indicates that the filtered setup omits a dominant saturation channel ... Self-consistent multi-mode benchmarks are needed to confirm the result survives when zonal dynamics are retained.

    Authors: We agree that the filtered single-n setup isolates thermal nonlinearities and that zonal fields provide an important additional saturation channel, as demonstrated by the factor-of-two enhancement in our simulations. The stiffness and 0.47% threshold are explicitly reported for the filtered case to highlight the role of thermal plasma effects without zonal-mode interference. We do not claim these quantitative values apply unchanged to full multi-mode scenarios. In the revised manuscript we have added a dedicated paragraph in the discussion clarifying this scope and noting that self-consistent multi-mode runs would be valuable for assessing combined effects. revision: partial

  2. Referee: The specific 0.47% threshold and stiffness are extracted from simulation time series without reported convergence tests (particle number, spatial resolution, time step), error bars on saturation levels, or systematic scans in gamma_L. This makes the quantitative claims provisional.

    Authors: We acknowledge that a full set of convergence tests, error bars, and exhaustive gamma_L scans were not reported in the original submission. The presented runs used standard PIC resolutions for this class of problem, and the weak gamma_L dependence was observed consistently across the available drive values. The revised manuscript now includes additional details on the numerical parameters employed and a note on observed sensitivity to resolution. A complete systematic convergence study with error bars would require a new series of simulations. revision: partial

  3. Referee: The frequency decrease with amplitude and saturation via continuum merging are presented as induced by PSZS, but the link is observational from time series rather than derived from an equation relating PSZS amplitude to the observed shift or merging criterion. A quantitative model would be required.

    Authors: The frequency evolution and continuum-merging saturation are interpreted using the theoretical analysis already present in the manuscript, which connects PSZS-induced modifications of the effective potential to the observed downshift and the standard continuum-damping criterion for merging. While a closed-form analytic expression for the exact frequency shift magnitude is not derived, the simulations exhibit a consistent causal sequence. We have expanded the relevant section with a step-by-step interpretive model linking PSZS amplitude to the frequency trajectory and mode transition. revision: yes

standing simulated objections not resolved
  • Self-consistent multi-mode benchmarks with zonal dynamics retained

Circularity Check

0 steps flagged

No significant circularity; results are direct simulation observations

full rationale

The paper reports saturation levels, frequency shifts, and stiffness directly from gyrokinetic PIC simulation time series in single-n runs with zonal fields filtered. No analytic derivation chain is presented that reduces the saturation amplitude or the gamma_L/omega_n threshold to a fitted parameter, self-cited uniqueness theorem, or self-referential definition. The PSZS explanation and continuum-merging description are interpretive accounts of observed phase-space structures rather than equations that loop back to the inputs. Zonal-field effects are explicitly contrasted in separate runs, and the filtering choice is stated as a deliberate isolation step without being smuggled in via prior self-citation as an unverified ansatz. The central claim therefore remains self-contained against external benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 1 invented entities

The central claims rest on standard gyrokinetic ordering and the assumption that observed PSZS structures are physical rather than numerical artifacts; no new free parameters are introduced beyond simulation resolution choices.

axioms (1)
  • domain assumption Gyrokinetic ordering remains valid for the thermal and energetic particle distributions studied
    Invoked implicitly by the choice of gyrokinetic PIC method throughout the simulations.
invented entities (1)
  • phase-space zonal structure (PSZS) no independent evidence
    purpose: Explains the observed frequency downshift of the TAE with increasing amplitude
    Identified in the simulations as the cause of the frequency decrease; no independent experimental or analytic confirmation is provided in the abstract.

pith-pipeline@v0.9.0 · 5545 in / 1431 out tokens · 53204 ms · 2026-05-10T09:23:19.834581+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

37 extracted references · 37 canonical work pages

  1. [1]

    stiffness

    The saturation level of TAE shows two regimes. For the linear drive below γL/ωn ≃ 0.47%, the nonlin- ear saturation is dominated by EP nonlinearity. Oth- erwise, the saturation of TAE is dominated by thermal plasma nonlinearity, i.e., PSZS of thermal species, which is eδϕn/Te ∼ 0.1 and almost independent of linear drive, showing feature of “stiffness” . He...

    work page

  2. [2]

    Upon mode saturation, the decrease of mode frequency and the sep- aration of m = 10 and m = 11 poloidal harmonics can be observed

    The saturation of TAE can happen when the mode amplitude reaches eδϕn/Te ∼ 0.1 even if EPs evolve lin- early, which implies a threshold on the amplitude of TAE imposed by thermal plasma nonlinearity. Upon mode saturation, the decrease of mode frequency and the sep- aration of m = 10 and m = 11 poloidal harmonics can be observed

    work page

  3. [3]

    As the mode merges to the continuum, the mode saturation and separation of neighboring poloidal harmonics can be identified

    The decrease of mode frequency results from the modification of the potential well of TAE by PSZS of thermal species, as indicated by gyrokinetic theory devel- oped herein. As the mode merges to the continuum, the mode saturation and separation of neighboring poloidal harmonics can be identified. Meanwhile, the saturation level can be quantitatively derived...

    work page

  4. [4]

    However, in this case, the effect of PSZS on TAE is significantly counteracted, leading to enhance- ment of the saturation level of TAE from eδϕn/Te ∼ 0.1 to eδϕn/Te ∼ 0.2

    For reasonable linear drive γL/ωn > 1%, the satu- ration of TAE is also dominated by thermal plasma non- linearities. However, in this case, the effect of PSZS on TAE is significantly counteracted, leading to enhance- ment of the saturation level of TAE from eδϕn/Te ∼ 0.1 to eδϕn/Te ∼ 0.2

    work page

  5. [5]

    nonlin- ear

    For both cases with/without zonal fields, the satu- ration level of TAE is found to be proportionate to the square root of inverse aspect ratio, which is predicted by the theory. This implies a stronger TAE activity in devices with a larger inverse aspect ratio. The results given above imply that for future toka- maks with much stronger EP drive, the satur...

    work page

  6. [6]

    linear” and “nonlinear

    It is found that both the growth rate and frequency of TAE increase with increasing EP concentration. The mode structures of poloidal harmonics of TAE are shown in Figure 4. It is found that for both scalar and vector potentials, the m = 10 and m = 11 poloidal harmonics are dominant, which is a typical characteristic of TAE. Besides, we can also find that ...

    work page

  7. [7]

    The temporal evolution of eδϕn/Te of n = 6 TAE is shown in Figure 5

    Nonlinear saturation due to EP nonlinearity Before investigating the nonlinear saturation of TAE with thermal plasma nonlinearities, we need to identify the saturation due to EP nonlinearity [26], i.e., case D in Table I. The temporal evolution of eδϕn/Te of n = 6 TAE is shown in Figure 5. Here, we can observe a lin- ear growth before ωA0t ∼ 400. After th...

    work page

  8. [8]

    stiffness

    Nonlinear saturation with thermal plasma nonlinearity With the nonlinear saturation of TAE due to EP nonlinearity investigated above, the nonlinear saturation 2 4 6 8 10 20 L (103s-1) 10-3 10-2 10-1 100 sat Asat ALL Asat EP Asat T L,C2 L,C1 sat L 2 sat L Figure 6: The scaling law of the saturation of TAE versus the linear growth rate γL in single- n = 6 s...

    work page 2000

  9. [9]

    Chen and F

    L. Chen and F. Zonca, Review of Modern Physics 88, 015008 (2016)

    work page 2016

  10. [10]

    Chen, Physics of Plasmas 1, 1519 (1994)

    L. Chen, Physics of Plasmas 1, 1519 (1994)

    work page 1994

  11. [11]

    Fasoli, C

    A. Fasoli, C. Gormenzano, H. Berk, B. Breizman, S. Briguglio, D. Darrow, N. Gorelenkov, W. Heidbrink, A. Jaun, S. Konovalov, et al., Nuclear Fusion 47, S264 (2007)

    work page 2007

  12. [12]

    Cheng, L

    C. Cheng, L. Chen, and M. Chance, Ann. Phys. 161, 21 (1985)

    work page 1985

  13. [13]

    G. Y. Fu and J. W. Van Dam, Physics of Fluids B 1, 1949 (1989)

    work page 1949

  14. [14]

    H. L. Berk and B. N. Breizman, Physics of Fluids B 2, 2226 (1990)

    work page 1990

  15. [15]

    H. L. Berk and B. N. Breizman, Physics of Fluids B 2, 2235 (1990)

    work page 1990

  16. [16]

    Barberis, V

    T. Barberis, V. N. Duarte, E. J. Hartigan-OConnor, and N. N. Gorelenkov, Nuclear Fusion 65, 112005 (2025)

    work page 2025

  17. [17]

    Zonca, F

    F. Zonca, F. Romanelli, G. Vlad, and C. Kar, Phys. Rev. Lett. 74, 698 (1995)

    work page 1995

  18. [18]

    L. Chen, F. Zonca, R. Santoro, and G. Hu, Plasma physics and controlled fusion 40, 1823 (1998)

    work page 1998

  19. [19]

    Z. Qiu, L. Chen, and F. Zonca, Reviews of Modern Plasma Physics 7 (2023)

    work page 2023

  20. [20]

    Zonca, L

    F. Zonca, L. Chen, S. Briguglio, G. Fogaccia, G. Vlad, and X. Wang, New Journal of Physics 17, 013052 (2015)

    work page 2015

  21. [21]

    Chen and F

    L. Chen and F. Zonca, Phys. Rev. Lett. 109, 145002 (2012)

    work page 2012

  22. [22]

    L. Chen, P. Liu, R. Ma, Z. Lin, Z. Qiu, W. Wang, and F. Zonca, Nuclear Fusion 65, 016018 (2024)

    work page 2024

  23. [23]

    Q. Fang, G. Wei, N. Chen, L. Chen, F. Zonca, and Z. Qiu, Nuclear Fusion 65, 066004 (2025)

    work page 2025

  24. [24]

    Chen and F

    L. Chen and F. Zonca, Nuclear Fusion 47, 886 (2007)

    work page 2007

  25. [25]

    S. Wang, Z. Wang, and T. Wu, Phys. Rev. Lett. 132, 065106 (2024)

    work page 2024

  26. [26]

    L. Chen, Z. Qiu, and F. Zonca, Physics of Plasmas 29, 050701 (2022)

    work page 2022

  27. [27]

    Cheng, K

    Z. Cheng, K. Shen, and Z. Qiu, Nuclear Fusion 64, 066031 (2024)

    work page 2024

  28. [28]

    Cheng, K

    Z. Cheng, K. Shen, L. Chen, F. Zonca, and Z. Qiu, Nu- clear Fusion 65, 082002 (2025)

    work page 2025

  29. [29]

    R. Ma, P. Liu, L. Chen, F. Zonca, and Z. Qiu (2025), 2512.09243

    work page arXiv 2025

  30. [30]

    Lanti, N

    E. Lanti, N. Ohana, N. Tronko, T. Hayward-Schneider, A. Bottino, B. McMillan, A. Mishchenko, A. Scheinberg, A. Biancalani, P. Angelino, et al., Computer Physics Communications 251, 107072 (2020)

    work page 2020

  31. [31]

    Mishchenko, A

    A. Mishchenko, A. Kies, R. Kleiber, and M. Cole, Physics of Plasmas 21, 092110 (2014)

    work page 2014

  32. [32]

    Mishchenko, A

    A. Mishchenko, A. Bottino, A. Biancalani, R. Hatzky, T. Hayward-Schneider, N. Ohana, E. Lanti, S. Brunner, L. Villard, M. Borchardt, et al., Computer Physics Com- munications 238, 194 (2019)

    work page 2019

  33. [33]

    Könies, S

    A. Könies, S. Briguglio, N. Gorelenkov, T. Fehér, M. Isaev, P. Lauber, A. Mishchenko, D. A. Spong, Y. Todo, W. A. Cooper, et al., Nuclear Fusion 58, 126027 (2018)

    work page 2018

  34. [34]

    Biancalani, A

    A. Biancalani, A. Bottino, M. Cole, C. D. Troia, P. Lauber, A. Mishchenko, B. Scott, and F. Zonca, Plasma Physics and Controlled Fusion 59, 054004 (2017)

    work page 2017

  35. [35]

    Y. Todo, H. Berk, and B. Breizman, Nuclear Fusion 50, 084016 (2010)

    work page 2010

  36. [36]

    Z. Qiu, L. Chen, and F. Zonca, Nuclear Fusion 57, 056017 (2017)

    work page 2017

  37. [37]

    E. A. Frieman and L. Chen, Physics of Fluids 25, 502 (1982)

    work page 1982