REVIEW 2 major objections 4 minor 38 references
Tightening energetic bounds on linear gyrokinetic instabilities
T0 review · 2 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read An energetic norm built from the complete linear-mode basis makes optimal-mode growth exactly equal eigenmode growth.
desk verdict The constrained optimal modes are the real contribution and look correct in slab geometry, but the paper's claim of a rigorous upper bound needs a proof that the quartic-root search actually finds the global maximum. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
Two objects carry the argument. The Case–Van Kampen energy is a diagonal norm in the complete set of linear modes; its balance law $dE/dt=2\sum_n\gamma_n|a_n|^2$ turns the optimal-mode variational problem (3.19) into the exact eigenvalue $\Lambda=\gamma_n$. The practical machinery is the constrained optimal-mode problem: maximize the free-energy drive $D$ subject to (4.9), (4.10), and (4.12), constraints that encode eigenmode phase relations between moments and free-energy balance; projecting the Euler–Lagrange equation onto the density and flow moments reduces it to the quartic (4.16), and the bound is $\Lambda_{\max}=\max_{\omega'_r}\Lambda$ over real roots.
What would settle it
A direct numerical search for a distribution function that satisfies constraints (4.9), (4.10), and (4.12) with free-energy growth above the quartic maximum would falsify the bound; equivalently, finding parameters where the slab ITG dispersion relation (A1) has a positive growth rate but the quartic (4.16) has no real positive root would show the constrained bound loses the instability.
Extended reading notes
Core claim
The central discovery is that loose energetic bounds are not a fundamental limitation: there exists a norm—the Case–Van Kampen energy $E=\sum_n |a_n|^2+\int |A(\omega)|^2 d\omega$—whose instantaneous optimal growth spectrum coincides with the linear eigenmode growth spectrum, giving equality in $\gamma \le \Lambda_{\max}$. Because the norm is diagonal in the eigenmode basis, transient growth is absent and the Landau-damped continuum does not decay the norm. For practical bounds, the paper shows that retaining just the real-frequency phase information through constrained optimal modes—maximising Helmholtz free-energy growth subject to eigenmode-like relations between density, flow, and higher moment plus free-energy balance $D=\gamma' H$—yields an upper bound that tracks the linear slab ITG growth rate, including a critical gradient and stabilization by density gradients that earlier nonlinear bounds missed.
Load-bearing premise
The load-bearing premise for the practical bound is that the feasible set defined by (4.9), (4.10), and (4.12) contains all linear eigenmodes, and that the maximum over real $\omega'_r$ of the real roots of (4.16) is the true supremum of free-energy growth on that set; the first part is plausible for eigenmodes, but the second is asserted rather than proved.
Editorial extensions
If this is right
- The slab ITG critical gradient $\kappa_{\parallel,\mathrm{cr}}=\sqrt{2\tau(1+\tau)}$ follows from a two-moment variational bound, with no need to solve the kinetic dispersion relation.
- Enlarging the number of moment constraints should tighten the bound toward the largest linear growth rate, because in the infinite-constraint limit the feasible set collapses to the linear eigenmodes.
- In any system where the Case–Van Kampen energy is a nonlinear invariant, linear stability implies nonlinear stability, since $\Lambda=0$ for all perpendicular wavenumbers means subcritical turbulence cannot occur.
- The constrained bound depends on the real frequency $\omega'_r$, so resonant stabilisation and density-gradient stabilisation enter through the optimization parameter rather than through an ad hoc closure.
Reading between the lines
- The same two-moment construction could be tried in toroidal geometry, where the continuum structure differs; a critical-gradient match there would be strong evidence the bound is capturing resonance physics rather than being a slab-specific accident.
- The gap between $\Lambda_{\max}$ and $\gamma_{\max}$ for simpler norms can now be read as a measure of non-normality, since the Case–Van Kampen result suggests the gap is a property of the chosen norm, not of the linear operator itself.
- A numerical implementation with an arbitrary number of Hermite–Laguerre moment constraints would give a tunable, rigorously valid upper bound that could be tested against linear gyrokinetic codes before being used in stellarator optimization.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper addresses the question of how tightly instantaneous energetic norms can bound the linear growth of gyrokinetic instabilities, using the slab ion-temperature-gradient (ITG) mode as a testbed. In Section 3, the authors reduce the linear gyrokinetic equation to a one-dimensional integral equation (3.1), invoke the Case-Van Kampen completeness theorem, and define a Case-Van Kampen energy E = sum_n |a_n|^2 + integral |A(omega)|^2 domega whose balance is dE/dt = 2 sum_n gamma_n |a_n|^2. The associated optimal-mode problem has solutions Lambda = gamma_n and Lambda = 0, so the tightest possible energetic bound coincides with the eigenmode growth spectrum, a result the authors themselves describe as a somewhat trivial consequence of diagonalisation. In Section 4, they construct constrained optimal modes that maximise Helmholtz free-energy growth subject to moment constraints and a free-energy balance constraint satisfied by linear eigenmodes. This leads to a quartic equation (4.16) for Lambda after maximising over the real frequency parameter omega'_r. Numerical comparisons in Figures 1-3 show that the constrained bound reproduces qualitative features of the linear growth rate, including the critical gradient kappa_{||,cr} = sqrt(2 tau (1+tau)), density-gradient stabilisation, and finite-Larmor-radius effects.
Significance. If the missing technical details are supplied, the paper would make two valuable contributions. First, it cleanly demonstrates that the looseness of previous energetic bounds is not a fundamental limitation: an energy norm built from eigenmode projection coefficients can make optimal growth and linear eigenmode growth coincide. Second, the constrained optimal mode construction is a low-dimensional, computationally efficient way to include real-frequency (phase) information in a variational bound, and the exact reproduction of the slab ITG critical gradient is a nontrivial and encouraging result. The paper also correctly situates the tightest-bound result as an existence statement about norms rather than a practical algorithm, since the Case-Van Kampen energy requires full spectral knowledge. The main load-bearing issue is that the 'rigorous upper bound' claim for the constrained optimal modes is not fully proven: the passage from the Lagrangian stationarity conditions to the global maximum over real omega'_r requires an attainment proof that is currently absent.
major comments (2)
- [Section 4, Eq. (4.16)] The identification of Lambda_max with the maximum over real omega'_r of the real roots of the quartic (4.16) is not justified as a rigorous upper bound. The Lagrangian (4.13) yields only first-order necessary conditions for interior extrema of the constrained problem; no compactness argument, constraint-qualification check, or analysis of boundaries (e.g., vanishing of one of the moments kappa_i or |omega'_r| tending to infinity) is provided, and complex roots are discarded without showing that the true supremum of D/H over the feasible set defined by (4.9), (4.10), and (4.12) is attained at a stationary point with a real Lambda. If the supremum is realised at a non-stationary point or at infinity, the computed Lambda_max could fall below the free-energy growth rate of a genuine linear eigenmode, invalidating the claimed bound. Please add an existence and attainment proof, or an explicit boundary and infinity analysis that rules out this failure.
- [Appendix B and Appendix C] The reduction from the Lagrangian (4.13) and the moment system (B 4)-(B 8) to the quartic (4.16) is only sketched; Appendix C lists the coefficients P, Q, and R but does not show the elimination steps. Since all numerical results in Figures 1-3 and the critical-gradient formula depend on equation (4.16), a full derivation (or a supplementary computer-algebra script) is needed for reproducibility and to verify that no non-generic cases, such as denominators involving 1 + tau - G_perp0 or eta vanishing, have been silently excluded.
minor comments (4)
- [Appendix B, final line] The definition 'tilde_lambda_3,4 = lambda_2,3 / (n T_i omega_* eta (1 - lambda_1))' appears to contain an index typo; presumably tilde_lambda_2 and tilde_lambda_3 are intended.
- [Section 4.1.1] The limit 'b -> 0 with eta -> infinity' is ambiguous because the plotted parameter kappa_|| = omega_* eta / (v_T k_||) is held fixed; please specify the ordering, e.g., eta -> infinity with omega_* -> 0 at fixed omega_* eta.
- [Section 3.2] The statement that the optimal modes are exactly the linear eigenmodes would benefit from a short direct verification that f_n satisfies the generalised eigenvalue problem (3.19); the projection argument shows the necessary value of Lambda but not explicitly that f_n is a solution.
- [Section 5] Because the abstract and introduction present the tightest possible bound as a main result, the conclusions should state explicitly that this bound is not practically computable without solving the linear problem, serving as an existence proof rather than a predictive tool.
Circularity Check
Section 3's 'tightest bound' is definitional: the Case–Van Kampen energy is built from eigenmode projection coefficients, so optimal growth equals eigenmode growth by construction; Section 4's constrained optimal modes are independent and non-circular.
-
self definitional
[Section 3.1–3.2, Eqs. (3.16)–(3.23)]
"by combining (3.14) and (3.15) we can construct what we will refer to as the Case–Van Kampen energy, which we define as, E = ∑_n |a_n(t)|^2 + ∫ |A(ω,t)|^2 dω. ... This may feel like a somewhat trivial consequence of the diagonalisation of the linear operator by the linear eigenmodes."
Equation (3.16) defines E directly from the coefficients a_n and A(ω) of the eigenfunction expansion (3.9), so the growth law (3.17) is merely the weighted sum of eigenvalues. The optimal-mode problem (3.19)–(3.21) is built from the same projection coefficients, and the orthogonality relations collapse the generalized eigenvalue problem to Λ = γ_n. Thus the 'tightest possible bound' and one-to-one correspondence are restatements of the diagonalization: the result is contained in the definition of E. The authors explicitly flag it as 'a somewhat trivial consequence,' confirming the reduction. This is transparent definitional circularity, not a fitted-parameter substitution.
full rationale
The Section 3 result is the only circular/definitional step. The Case–Van Kampen energy E is defined as the sum of squared projection coefficients of the linear eigenmode basis, so the balance dE/dt = 2∑γ_n|a_n|^2 and the optimal-mode solution Λ=γ_n are already contained in the construction. The authors themselves call it 'a somewhat trivial consequence of the diagonalisation,' which confirms the reduction. This does not invalidate the paper: it is an honest, explicitly disclosed identity, and the subsequent constrained-optimal-mode construction in Section 4 is genuinely non-circular—its constraints are necessary conditions for eigenmodes, not the eigenmode solution itself, and the resulting bounds are compared to the independent linear dispersion relation (A1). The skeptic's concern that Λ_max might not be attained by real roots of (4.16) is a rigor issue about the variational maximum, not a circularity. Overall score 6 reflects that one of the two headline claims is true by construction, while the practically useful bound is independent and self-contained.
Assumptions & free parameters
free parameters (2)
- omega'_r (constrained-mode real frequency)
- Arbitrary weights in Case-Van Kampen energy =
arbitrary
assumptions (4)
- standard math The discrete plus continuum Case-Van Kampen eigenmodes form a complete basis for the linear operator of Eq. (3.1), with orthogonality relations (3.10)-(3.11).
- domain assumption Slab geometry reduction: uniform magnetic field, single Fourier mode k_parallel, adiabatic electrons, single ion species.
- domain assumption Every linear eigenmode satisfies the moment constraints (4.9)-(4.10) and energy balance D = gamma H (4.12).
- ad hoc to paper The global maximum over omega'_r of the real roots of (4.16) yields the true supremum of Lambda over the constrained set.
invented entities (1)
-
Case-Van Kampen energy E = sum_n |a_n|^2 + integral |A(omega)|^2 domega
Cite this review
Pith. "Pith review of Tightening energetic bounds on linear gyrokinetic instabilities." pith.science (2026). https://pith.science/paper/2XNOVH4N
@misc{pith2026250517757,
author = {Pith},
title = {Pith review of: Tightening energetic bounds on linear gyrokinetic instabilities},
year = {2026},
howpublished = {\url{https://pith.science/paper/2XNOVH4N}},
note = {Machine review of arXiv:2505.17757}
}
read the original abstract
Bounding energetic growth of gyrokinetic instabilities is a complementary approach to linear instability analyses involving normal eigenmodes. Previous work has focused on upper bounds which are valid linearly and nonlinearly. However, if an upper bound on linear instability growth is desired, these nonlinearly valid bounds may be a poor predictor of the growth of the most unstable eigenmode. This is most evident for the simplest of instabilities: the ion-temperature-gradient (ITG) mode in slab geometry. In this work, we derive energetic upper bounds specifically for linear instability growth, focusing on the slab ITG. We show that there is no fundamental limitation on how tightly linear growth can be bounded by an energetic norm, with the tightest possible bound being given by a special energy comprised of projection coefficients of the linear eigenmode basis. Additionally, we consider `constrained optimal modes' that maximise energy growth subject to constraints that are also obeyed by the linear eigenmodes. This yields computationally efficient upper bounds that closely resemble the linear growth rate, capturing effects connected to the real frequency of instabilities, which have been absent in the energetic bounds considered thus far.
Figures
Reference graph
Works this paper leans on
-
[1]
, " * write output.state after.block = add.period write newline
ENTRY address author booktitle chapter edition editor howpublished institution journal key month note number organization pages publisher school series title type volume year eprint label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 'mid.sentence ...
-
[2]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in capitalize " " * FUNCT...
-
[3]
Barbour, Nathaniel , Dorland, William , Abel, Ian G. & Landreman, Matt 2025 Machine-learning Closure for Vlasov-Poisson Dynamics in Fourier-Hermite Space , arXiv:arXiv: 2504.13313
-
[4]
Barnes, M. , Parra, F. I. , Highcock, E. G. , Schekochihin, A. A. , Cowley, S. C. & Roach, C. M. 2011 Turbulent Transport in Tokamak Plasmas with Rotational Shear . Phys. Rev. Lett. 106 (17), 175004
work page 2011
-
[5]
Biglari, H. , Diamond, P. H. & Rosenbluth, M. N. 1989 Toroidal ion-pressure-gradient-driven drift instabilities and transport revisited . Physics of Fluids B: Plasma Physics 1 (1), 109--118
work page 1989
-
[6]
Bourdelle, C. , Garbet, X. , Imbeaux, F. , Casati, A. , Dubuit, N. , Guirlet, R. & Parisot, T. 2007 A new gyrokinetic quasilinear transport model applied to particle transport in tokamak plasmas . Physics of Plasmas 14 (11), 112501
work page 2007
-
[7]
Braginskii, S. I. 1965 Transport Processes in a Plasma . Reviews of Plasma Physics 1 , 205
work page 1965
-
[8]
Case, K. M 1959 Plasma oscillations . Annals of Physics 7 (3), 349--364
work page 1959
Show all 38 references
-
[9]
Costello, P. J. & Plunk, G. G. 2024 Energetic bounds on gyrokinetic instabilities. Part 4. Bounce-averaged electrons, arXiv:arXiv: 2404.06081
2024 arXiv
-
[10]
Physics of Plasmas 12 (7), 072309
Dannert, Tilman & Jenko, Frank 2005 Gyrokinetic simulation of collisionless trapped-electron mode turbulence . Physics of Plasmas 12 (7), 072309
2005
-
[11]
Journal of the Atmospheric Sciences 61 (9), 1086--1091
DelSole, Timothy 2004 The Necessity of Instantaneous Optimals in Stationary Turbulence . Journal of the Atmospheric Sciences 61 (9), 1086--1091
2004
-
[12]
Frei, B. J. , Hoffmann, A. C. D. , Ricci, P. , Brunner, S. & Tecchioll, Z. 2023 Moment-based approach to the flux-tube linear gyrokinetic model . Journal of Plasma Physics 89 (4), 905890414
2023
-
[13]
& Perkins, Francis W
Hammett, Gregory W. & Perkins, Francis W. 1990 Fluid moment models for Landau damping with application to the ion-temperature-gradient instability . Phys. Rev. Lett. 64 (25), 3019--3022
1990
-
[14]
Hatch, D R , Jenko, F , Navarro, A Ba \ n \'o n , Bratanov, V , Terry, P W & Pueschel, M J 2016 Linear signatures in nonlinear gyrokinetics: Interpreting turbulence with pseudospectra . New J. Phys. 18 (7), 075018
2016
-
[15]
Hatch, D. R. , Terry, P. W. , Jenko, F. , Merz, F. & Nevins, W. M. 2011 Saturation of Gyrokinetic Turbulence through Damped Eigenmodes . Phys. Rev. Lett. 106 (11), 115003
2011
-
[16]
& Plunk, G
Helander, P. & Plunk, G. G. 2021 Upper Bounds on Gyrokinetic Instabilities in Magnetized Plasmas . Phys. Rev. Lett. 127 (15), 155001
2021
-
[17]
& Plunk, G
Helander, P. & Plunk, G. G. 2022 Energetic bounds on gyrokinetic instabilities. Part 1. Fundamentals . Journal of Plasma Physics 88 (2), 905880207
2022
-
[18]
Heninger, J. M. & Morrison, P. J. 2018 An integral transform technique for kinetic systems with collisions . Physics of Plasmas 25 (8), 082118
2018
-
[19]
Huang, Ziyu , Dong, Chuanfei & Wang, Liang 2025 Machine-learning heat flux closure for multi-moment fluid modeling of nonlinear Landau damping . Proc. Natl. Acad. Sci. U.S.A. 122 (11), e2419073122 , arXiv:arXiv: 2503.11090
2025 arXiv
-
[20]
Kadomtsev, B. B. & Pogutse, O. P. 1970 Turbulence in Toroidal Systems . In Reviews of Plasma Physics \/ (ed. M. A. Leontovich ) , pp. 249--400 . Boston, MA: Springer US
1970
-
[21]
, Merz, F
Kammerer, M. , Merz, F. & Jenko, F. 2008 Exceptional points in linear gyrokinetics . Physics of Plasmas 15 (5), 052102 , arXiv:arXiv: https://pubs.aip.org/aip/pop/article-pdf/doi/10.1063/1.2909618/14083208/052102\_1\_online.pdf
2008 doi
-
[22]
, Liu, X
Kotschenreuther, M. , Liu, X. , Mahajan, S.M. , Hatch, D.R. & Merlo, G. 2024 Transport Barriers in magnetized plasmas- general theory with dynamical constraints . Nucl. Fusion 64 (7), 076033
2024
-
[23]
1997 Submarginal profiles and turbulent transport: An exactly solvable model
Krommes, John A. 1997 Submarginal profiles and turbulent transport: An exactly solvable model . Physics of Plasmas 4 (5), 1342--1356
1997
-
[24]
1946 On electron plasma oscillations
Landau, L.D. 1946 On electron plasma oscillations . Sov. Phys. JETP 16 , 574
1946
-
[25]
, Plunk, G
Landreman, M. , Plunk, G. G. & Dorland, W. 2015 Generalized universal instability: Transient linear amplification and subcritical turbulence . Journal of Plasma Physics 81 (5), 905810501
2015
-
[26]
Mandell, N. R. , Dorland, W. , Abel, I. , Gaur, R. , Kim, P. , Martin, M. & Qian, T. 2024 GX : A GPU-native gyrokinetic turbulence code for tokamak and stellarator design . Journal of Plasma Physics 90 (4), 905900402
2024
-
[27]
Mandell, N. R. , Dorland, W. & Landreman, M. 2018 Laguerre-- Hermite pseudo-spectral velocity formulation of gyrokinetics . Journal of Plasma Physics 84 (1), 905840108
2018
-
[28]
Morrison, P. J. & Pfirsch, D. 1992 Dielectric energy versus plasma energy, and Hamiltonian action-angle variables for the Vlasov equation . Physics of Fluids B: Plasma Physics 4 (10), 3038--3057
1992
-
[29]
Plunk, G. G. 2013 Landau damping in a turbulent setting . Physics of Plasmas 20 (3), 032304
2013
-
[30]
Plunk, G. G. 2015 On the nonlinear stability of a quasi-two-dimensional drift kinetic model for ion temperature gradient turbulence . Physics of Plasmas 22 (4), 042305
2015
-
[31]
Plunk, G. G. & Helander, P. 2022 Energetic bounds on gyrokinetic instabilities. Part 2. Modes of optimal growth . Journal of Plasma Physics 88 (3), 905880313
2022
-
[32]
Plunk, G. G. & Helander, P. 2023 Energetic bounds on gyrokinetic instabilities. Part 3. Generalized free energy . Journal of Plasma Physics 89 (4), 905890419
2023
-
[33]
Plunk, G. G. , Helander, P. , Xanthopoulos, P. & Connor, J. W. 2014 Collisionless microinstabilities in stellarators. III . The ion-temperature-gradient mode . Physics of Plasmas 21 (3), 032112
2014
-
[34]
, Plunk, G
Podavini, L. , Plunk, G. G. , Helander, P. & Zocco, A. 2024 Energetic bounds on gyrokinetic instabilities. Part V . Contrasting optimal and normal modes over the geometric landscape . To be published in Journal of Plasma Physics
2024
-
[35]
Pueschel, M. J. , Faber, B. J. , Citrin, J. , Hegna, C. C. , Terry, P. W. & Hatch, D. R. 2016 Stellarator Turbulence : Subdominant Eigenmodes and Quasilinear Modeling . Phys. Rev. Lett. 116 (8), 085001
2016
-
[36]
Roberg-Clark , G. T. , Plunk, G. G. , Xanthopoulos, P. , N \"u hrenberg, C. , Henneberg, S. A. & Smith, H. M. 2023 Critical gradient turbulence optimization toward a compact stellarator reactor concept . Phys. Rev. Res. 5 (3), L032030
2023
-
[37]
& Felderhof, B.U
van Kampen , N.G. & Felderhof, B.U. 1967 Theoretical Methods in Plasma Physics\/ . North-Holland Publishing Company
1967
-
[38]
, Highcock, E
van Wyk, F. , Highcock, E. G. , Schekochihin, A. A. , Roach, C. M. , Field, A. R. & Dorland, W. 2016 Transition to subcritical turbulence in a tokamak plasma . Journal of Plasma Physics 82 (6), 905820609
2016
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.