{"id":"a885728d-00fb-4e33-aa83-38b4c5e9ca74","arxiv_id":"2505.17757","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A constrained optimal mode method gives upper bounds on slab ITG linear growth that track the true growth rate and critical gradient, while the tightest possible bound is shown to be a definitional Case-Van Kampen energy.","lead":"This paper shows how to construct energetic upper bounds on linear gyrokinetic instabilities that are much tighter than earlier nonlinear bounds, and demonstrates them for the slab ion-temperature-gradient mode. The practical variant, called constrained optimal modes, reproduces the linear growth rate's dependence on key parameters, including the critical gradient, using only a few moment constraints.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Constrained-optimum bound lacks a proof that the global maximum of free-energy growth is attained and captured by the real roots of (4.16); if the maximum occurs at a non-stationary point or at infinity, the claimed rigorous upper bound could fail.","rationale":"The reader identified the same weakest assumption: the practical bound depends on the unproved global-maximization and real-root-selection step. I concur. The exact-bound result is a construction that is definitional once Case–Van Kampen completeness is accepted, and the authors flag this. The constrained-optimal-mode bound is the main novel claim, and its status as a rigorous upper bound requires that the algebraic reduction to (4.16) captures the global maximum of D/H over the feasible set. The paper supplies stationary equations and numerical examples, but no existence or compactness argument, so the concern is genuine and load-bearing. The proposed test—re-deriving the algebra and directly verifying the bound against the linear dispersion relation and a direct optimization—would settle whether the concern lands. No change to the conditional verdict is needed.","tokens_in":14091,"tokens_out":17243,"duration_ms":148888,"concrete_test":"Independently re-derive the reduction from the Lagrangian (4.13) to (4.16) using a computer algebra system, verifying that no division by κ1, γ', or λ1−1 eliminates valid cases. Then, for a parameter scan (for example, τ ∈ [0.5, 2], η ∈ [1, 10], κ∥ ∈ [0.1, 10], b ∈ [0, 2]), compute every linear eigenmode from (A.1) and check that each eigenfrequency (ω_r, γ) satisfies (4.16) with ω'_r = ω_r. In addition, directly maximize D/H over the three-moment feasible set with a global optimizer for a few challenging cases, including near-critical parameters, and verify that the direct maximum never exceeds Λ_max computed from the quartic. If any eigenmode is missed or any direct maximum exceeds the quartic bound, the upper-bound claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The practical upper bound rests on identifying Λ_max with the maximum over real ω'_r of real roots of the quartic (4.16). These roots are necessary conditions for interior extrema of the Lagrangian (4.13), not a demonstration that the constrained variational problem attains a global maximum on the feasible set defined by (4.9), (4.10) and (4.12). No compactness argument, constraint-qualification check, or boundary analysis is provided, and complex roots are discarded without showing that the true supremum cannot occur at a degenerate or non-stationary point. If the true supremum of D/H over the feasible set is larger than the largest real root—for example, at ω'_r → ∞ or at a boundary where some moment vanishes—then Λ_max could fall below the free-energy growth rate of a genuine linear eigenmode, invalidating the claimed rigorous bound. The numerical agreement in Figures 1–3 is suggestive but only samples a limited parameter range. The exact-bound result in Section 3 is definitional and self-flagged, and its reliance on Case (1959) completeness is acceptable; the load-bearing gap is in Section 4.2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","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.","tokens_in":14376,"tokens_out":6009,"duration_ms":50632,"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":[{"comment":"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.","section":"Section 4, Eq. (4.16)"},{"comment":"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.","section":"Appendix B and Appendix C"}],"minor_comments":[{"comment":"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":"Appendix B, final line"},{"comment":"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":"Section 4.1.1"},{"comment":"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":"Section 3.2"},{"comment":"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.","section":"Section 5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is within the scope of Journal of Plasma Physics, and the central idea is promising. The two major concerns are fixable: a rigorous attainment proof for the constrained variational problem and a complete derivation of the quartic. I would not recommend rejection on the current evidence, but the phrase 'rigorous upper bound' should not appear in the final version until the global-maximum issue is resolved. The authors' frank acknowledgement that the tightest-bound result is a somewhat trivial consequence of diagonalisation is commendable and should be retained."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The genuinely new thing here is the constrained optimal mode construction. Adding moment-equation constraints to the free-energy maximization captures the real-frequency effects that previous energetic bounds missed, and it reproduces the slab ITG critical gradient and density-gradient stabilization surprisingly well. That is a solid methodological advance, and the numerical agreement in Figures 1–3 is convincing. The Case–Van Kampen energy section is neat, but as the authors themselves say, it is a somewhat trivial consequence of diagonalizing the linear operator; it requires complete knowledge of the spectrum and so is more of a theoretical benchmark than a practical tool. I do not see that as a flaw, only a matter of emphasis.\n\nThe load-bearing gap is in Section 4.2. The bound is claimed to be rigorous, but the actual procedure maximizes over the real roots of a quartic (4.16) without showing that the true supremum of D/H on the constrained set is attained at one of those roots. No compactness argument, boundary analysis, or constraint-qualification check is given. The stress-test concern is fair: if the supremum occurs at a non-stationary point or at infinity, the computed Λ_max could understate the growth of a genuine eigenmode. I do not have a counterexample, and the numerics suggest the procedure is right, but the rigor claim currently outruns the proof. A referee should ask for either a compactness/attainment argument or a more careful statement that the bound is conjectured rather than proven.\n\nAlso, the reduction to the quartic is only sketched; the appendices give coefficients but not the elimination steps. That is probably fine for a specialist journal, but an independent check would be useful. The paper is honest about limitations, the literature is well covered, and the connection to gyrofluid closures is appropriately drawn.\n\nWho is this for? Plasma theorists working on gyrokinetic stability, especially stellarator optimization where cheap rigorous bounds would be valuable. It deserves a serious referee and likely publication after revision. I would send it to peer review without hesitation, but I would not accept the 'rigorous' language until the global-maximization step is either proven or explicitly qualified.","headline":"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.","tokens_in":14872,"tokens_out":2932,"would_cite":true,"duration_ms":27911,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"An energetic norm built from the complete linear-mode basis makes optimal-mode growth exactly equal eigenmode growth.","keywords":["energetic bounds","gyrokinetics","slab ITG","Case–Van Kampen modes","optimal modes","free energy","linear growth rate","critical gradient"],"falsifier":"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.","tokens_in":13868,"feed_emoji":"⚡","tokens_out":6689,"duration_ms":48657,"temperature":0.7,"pith_summary":"This paper asks how tightly linear gyrokinetic instability growth can be bounded from above by the growth of optimal modes of an energetic norm, and it sharpens the answer in two steps. The authors first construct the Case–Van Kampen energy, a positive norm built from the projection coefficients of the perturbed distribution onto the complete basis of discrete and continuum linear modes, and prove that its optimal-mode growth rates are exactly the linear eigenmode growth rates: $\\Lambda=\\gamma_n$ for each discrete mode and $\\Lambda=0$ for the continuum. They then abandon the need for full spectral knowledge and derive constrained optimal modes that maximise Helmholtz free-energy growth subject only to two moment constraints and an energy-balance constraint that all linear eigenmodes satisfy. Solving that variational problem reduces to a quartic polynomial whose maximum over the real frequency reproduces the slab ITG critical gradient $\\kappa_{\\parallel,\\mathrm{cr}}=\\sqrt{2\\tau(1+\\tau)}$, as well as density-gradient stabilisation and finite-Larmor-radius dependence. If these claims hold, energetic bounding becomes a tight and computationally cheap way to predict linear instability thresholds and growth rates without solving the eigenmode problem.","feed_headline":"Energetic bound now tracks slab ITG growth, critical gradient included","feed_subtitle":"Constrained optimal modes capture real-frequency resonance effects and yield a cheap, tight upper bound on linear growth.","key_machinery":"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.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"Supplies the completeness theorem for the discrete-plus-continuum eigenmode basis that makes the Case–Van Kampen energy a positive definite norm.","marker":"Case (1959)"},{"why":"Provides the theoretical treatment of singular Van Kampen modes and their normalization used in the continuum projection.","marker":"van Kampen & Felderhof (1967)"},{"why":"The previous generalized free-energy bound for the slab ITG that this paper shows is loose and then tightens with constrained optimal modes.","marker":"Plunk & Helander (2023)"},{"why":"Establishes the Helmholtz free-energy norm and its balance, the energetic quantity maximised by the constrained optimal modes.","marker":"Helander & Plunk (2022)"},{"why":"Source of the slab ITG critical gradient expression that the constrained bound reproduces.","marker":"Kadomtsev & Pogutse (1970)"},{"why":"Provides the linear dispersion relation and fluid-limit comparison used to validate the bound against eigenmode growth.","marker":"Plunk et al. (2014)"},{"why":"Example of a drift-kinetic system where the Case–Van Kampen energy is a nonlinear invariant, supporting the nonlinear-stability consequence.","marker":"Plunk (2015)"},{"why":"Related dynamically constrained free-energy approach that the authors compare their constrained optimal modes to.","marker":"Kotschenreuther et al. (2024)"}],"fun_headline_variants":["Tight linear bound for slab ITG instability","Constrained modes tighten gyrokinetic growth bound","Real-frequency effects sharpen energetic bound","Slab ITG growth now tightly bounded energetically","New energy norm matches linear ITG growth"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["Tight linear bound for slab ITG instability","Constrained modes tighten gyrokinetic growth bound","Real-frequency effects sharpen energetic bound","Slab ITG growth now tightly bounded energetically","New energy norm matches linear ITG growth"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000181,"raw_usage":{"total_tokens":1300,"prompt_tokens":931,"completion_tokens":369,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":547,"completion_tokens_details":{"reasoning_tokens":299}},"tokens_in":547,"tokens_out":369,"duration_ms":3222,"temperature":1.0,"reasoning_tokens":299,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T14:41:38.916694+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"M 1959 Plasma oscillations","cited_arxiv_id":null,"evidence_quote":"Supplies the completeness theorem for the discrete-plus-continuum eigenmode basis that makes the Case–Van Kampen energy a positive definite norm."},{"cited_title":"& Felderhof, B.U","cited_arxiv_id":null,"evidence_quote":"Provides the theoretical treatment of singular Van Kampen modes and their normalization used in the continuum projection."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"The previous generalized free-energy bound for the slab ITG that this paper shows is loose and then tightens with constrained optimal modes."},{"cited_title":"& Plunk, G","cited_arxiv_id":null,"evidence_quote":"Establishes the Helmholtz free-energy norm and its balance, the energetic quantity maximised by the constrained optimal modes."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Source of the slab ITG critical gradient expression that the constrained bound reproduces."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the linear dispersion relation and fluid-limit comparison used to validate the bound against eigenmode growth."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Example of a drift-kinetic system where the Case–Van Kampen energy is a nonlinear invariant, supporting the nonlinear-stability consequence."},{"cited_title":", Liu, X","cited_arxiv_id":null,"evidence_quote":"Related dynamically constrained free-energy approach that the authors compare their constrained optimal modes to."}],"review_version":1}