pith. sign in

arxiv: 2605.16814 · v1 · pith:EMZRRXUDnew · submitted 2026-05-16 · ⚛️ physics.plasm-ph

Critical velocity-space mode scalings in linear and nonlinear Landau damping for the Vlasov--Poisson system

Noah K. Guberman , J. Coughlin , A. S. Joglekar This is my paper

Pith reviewed 2026-05-19 19:39 UTC · model grok-4.3

classification ⚛️ physics.plasm-ph
keywords Landau dampingVlasov-Poisson systemvelocity-space resolutioncollisional diffusionFourier-Hermite modesbounce frequencynonlinear plasma physics
0
0 comments X

The pith

Analytical scalings give the upper bound on velocity-space resolution needed for linear and nonlinear Landau damping in Vlasov-Poisson simulations.

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

The paper derives analytical scalings for the critical Fourier and Hermite velocity-space mode numbers required to capture Landau damping accurately when collisions are present. It employs a unified cascade-balance argument that equates the effects of wave-particle trapping and collisional diffusion to set these resolution limits. The scalings are expressed in terms of the bounce frequency, wavenumber, and collisional frequency. An ensemble of 800 Vlasov-Fokker-Planck simulations confirms the predicted dependencies on bounce frequency and collision frequency for both linear and nonlinear regimes.

Core claim

A unified cascade-balance argument yields analytical scalings for the critical Fourier and Hermite velocity-space mode numbers that bound the resolution needed for linear and nonlinear Landau damping mediated by collisional diffusion in the Vlasov-Poisson system; these scalings depend explicitly on the bounce frequency ω_b, wavenumber kλ_D, and electron-electron collisional frequency ν, with strong numerical agreement in the ω_b and ν dependencies.

What carries the argument

The cascade-balance argument that equates wave-particle bounce dynamics with collisional diffusion to identify the dominant balance setting critical mode numbers across regimes.

If this is right

  • Resolution requirements tighten as collisional frequency decreases because weaker diffusion allows finer velocity-space structures to persist longer.
  • Nonlinear regimes demand higher mode numbers than linear ones once bounce frequency exceeds a threshold set by the wavenumber.
  • The scalings provide an a priori estimate that removes the need for empirical resolution tests in Vlasov-Poisson codes.
  • Dependence on kλ_D implies that shorter-wavelength modes require proportionally finer velocity-space grids.

Where Pith is reading between the lines

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

  • These scalings could be used to design adaptive velocity-space grids that refine only where the local collisional diffusion length demands it.
  • Similar cascade-balance reasoning might extend to other velocity-space instabilities such as two-stream or bump-on-tail problems.
  • In multi-dimensional simulations the same logic would predict resolution needs along each velocity coordinate separately.

Load-bearing premise

The cascade-balance argument correctly identifies the dominant physical balance that sets the critical mode numbers across both linear and nonlinear regimes when collisions are present.

What would settle it

A high-resolution reference simulation run at mode numbers below the predicted critical values exhibits a measurable deviation in the measured damping rate or saturation level compared with runs that meet or exceed the scaling.

Figures

Figures reproduced from arXiv: 2605.16814 by A. S. Joglekar, J. Coughlin, Noah K. Guberman.

Figure 1
Figure 1. Figure 1: Left: linear (small-amplitude) phase mixing showing filamented structure in the mean-subtracted distribution function at t = 800 ω −1 p . Right: nonlinear (large-amplitude) phase mixing showing a pronounced vortex and its separatrix boundary in the total distribution function at t = 1200 ω −1 p . trapping. These two qualitatively distinct velocity-space dynamics motivate separate derivations of the critica… view at source ↗
Figure 2
Figure 2. Figure 2: Velocity-space cascade for a representative linear case (ωb = 0.009, kλD = 0.30, νˆ = 1.0×10−4 ). Left: log10 | ˆfk1 (t, s)| as a function of time and the dimensionless Fourier velocity wavenumber s. Right: log10 | ˆfk1 (t, m)| as a function of time and the square root of the Hermite mode number. In both panels, the cascade front advances linearly to higher mode numbers as phase mixing proceeds, before bei… view at source ↗
Figure 3
Figure 3. Figure 3: Same as [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Measured critical mode nc versus the theoretical prediction on log-log axes. The best-fit line (blue) corresponds to Eq. 4.1 evaluated using the best-fit C with fixed theoretical exponents. The four panels show, clockwise from top-left: sc (linear Fourier), mc (linear Hermite), mc (nonlinear Hermite), and sc (nonlinear Fourier). The corresponding R 2 values are given in table 1. Simulations where the casca… view at source ↗
read the original abstract

The velocity-space resolution required to accurately simulate kinetic phenomena in the 1D-1V Vlasov--Poisson system is generally not known a priori. In this work, we determine the upper bound on the resolution requirement for linear and nonlinear Landau damping mediated by collisional diffusion, deriving analytical scalings for the critical Fourier and Hermite velocity-space mode numbers using a unified cascade-balance argument. The resulting scalings depend on the bounce frequency $\omega_b$, wavenumber $k\lambda_D$, and electron-electron collisional frequency $\nu$. We validate these predictions against an ensemble of 800 Vlasov--Fokker--Planck simulations, finding strong agreement with the predicted $\omega_b$ and $\nu$ dependencies.

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

1 major / 2 minor

Summary. The paper derives analytical scalings for the critical Fourier and Hermite velocity-space mode numbers required to resolve linear and nonlinear Landau damping in the 1D-1V Vlasov-Poisson system. The scalings are obtained from a unified cascade-balance argument that equates the collisional diffusion time at the cutoff mode to the inverse bounce frequency ω_b (nonlinear) or linear Landau rate, with explicit dependence on wavenumber kλ_D and collisional frequency ν. These predictions are tested against an ensemble of 800 Vlasov-Fokker-Planck simulations, which reportedly show agreement with the predicted ω_b and ν dependencies.

Significance. If the scalings are robust, the work supplies a practical a priori estimate for velocity-space resolution in kinetic plasma simulations of Landau damping, a core process in collisionless and weakly collisional plasmas. The use of a large simulation ensemble to test the ν and ω_b scalings is a clear strength, as is the attempt at a unified treatment spanning linear and nonlinear regimes. The result could directly inform grid design in Vlasov codes and reduce unnecessary computational cost.

major comments (1)
  1. [Cascade-balance derivation] Cascade-balance argument (described in the abstract and the central derivation): the claim that equating collisional diffusion time at the critical Hermite/Fourier mode to 1/ω_b (or the linear rate) isolates the dominant cutoff mechanism is not yet demonstrated when phase mixing and trapping operate simultaneously. The 800-simulation ensemble confirms the overall ν and ω_b dependence but does not distinguish whether the assumed balance is realized or whether residual linear phase mixing or renormalized diffusion from trapped particles sets the cutoff instead. This assumption is load-bearing for the unified scaling across regimes.
minor comments (2)
  1. The abstract states 'strong agreement' with the ensemble but does not report quantitative measures (e.g., fit residuals or percentage deviation) for the predicted scalings; adding these would strengthen the validation claim.
  2. Notation for the critical mode numbers (Fourier and Hermite) should be introduced with explicit equations early in the text to aid readability for readers unfamiliar with the cascade argument.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful review and constructive comments on our manuscript. We address the single major comment below and indicate the revisions made.

read point-by-point responses

  1. Referee: Cascade-balance argument (described in the abstract and the central derivation): the claim that equating collisional diffusion time at the critical Hermite/Fourier mode to 1/ω_b (or the linear rate) isolates the dominant cutoff mechanism is not yet demonstrated when phase mixing and trapping operate simultaneously. The 800-simulation ensemble confirms the overall ν and ω_b dependence but does not distinguish whether the assumed balance is realized or whether residual linear phase mixing or renormalized diffusion from trapped particles sets the cutoff instead. This assumption is load-bearing for the unified scaling across regimes.

    Authors: We thank the referee for identifying this key point about the derivation. The cascade-balance argument is constructed by setting the collisional diffusion time at the cutoff mode equal to the inverse of the driving frequency scale (ω_b or the linear Landau rate), which identifies when dissipation overtakes the phase-space cascade. We agree that the 800 simulations primarily confirm the resulting scalings rather than furnishing direct evidence of the balance itself in regimes where phase mixing and trapping coexist. Nevertheless, the observed dependence on both ν and ω_b matches the specific functional form predicted by the balance and would be unlikely to arise from residual linear phase mixing or trapped-particle renormalized diffusion alone, which generally produce different parametric trends. In the revised manuscript we have expanded the central derivation section to include a clearer timescale-separation argument showing why the collisional cutoff dominates at the critical modes, together with representative plots of Hermite/Fourier mode energy evolution that illustrate the point at which collisional damping overtakes the cascade. revision: partial

Circularity Check

0 steps flagged

No circularity: scalings derived from independent physical balance and validated externally

full rationale

The paper presents a cascade-balance argument that equates collisional diffusion time at critical Hermite/Fourier modes to the inverse bounce frequency (or linear Landau rate), yielding explicit scalings in ω_b, kλ_D, and ν. This modeling step is a physical ansatz, not a redefinition or fit of the target mode numbers themselves. The resulting expressions are then compared to an independent ensemble of 800 Vlasov–Fokker–Planck simulations whose outcomes are not used to tune the balance. No self-citation chain, no fitted parameter renamed as prediction, and no reduction of the claimed scalings to their own inputs by construction appear in the derivation.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the validity of the cascade-balance argument as the mechanism that determines resolution requirements; no explicit free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption A unified cascade-balance argument captures the dominant physics setting the critical Fourier and Hermite mode numbers in both linear and nonlinear Landau damping when collisional diffusion is present.
    This assumption is invoked to derive the analytical scalings reported in the abstract.

pith-pipeline@v0.9.0 · 5666 in / 1336 out tokens · 60843 ms · 2026-05-19T19:39:44.489521+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

30 extracted references · 30 canonical work pages

  1. [1]

    Physics of Plasmas , author =

    Coulomb collision effects on linear. Physics of Plasmas , author =. 2014 , pages =. doi:10.1063/1.4875726 , abstract =

    work page doi:10.1063/1.4875726 2014

  2. [2]

    Machine learning of hidden variables in multiscale fluid simulation , volume =

    Joglekar, Archis S and Thomas, Alexander G R , month = sep, year =. Machine learning of hidden variables in multiscale fluid simulation , volume =. Machine Learning: Science and Technology , publisher =. doi:10.1088/2632-2153/acf81a , abstract =

    work page doi:10.1088/2632-2153/acf81a

  3. [3]

    Journal of Plasma Physics , author =

    Weakly collisional steady state linear and nonlinear plasma waves , volume =. Journal of Plasma Physics , author =. 2025 , keywords =. doi:10.1017/S0022377825000315 , abstract =

    work page doi:10.1017/s0022377825000315 2025

  4. [4]

    Collisionless

    O'Neil, Thomas , month = dec, year =. Collisionless. The Physics of Fluids , publisher =. doi:10.1063/1.1761193 , abstract =

    work page doi:10.1063/1.1761193

  5. [5]

    and Bernstein, Ira B

    Lenard, A. and Bernstein, Ira B. , month = dec, year =. Plasma. Physical Review , publisher =. doi:10.1103/PhysRev.112.1456 , abstract =

    work page doi:10.1103/physrev.112.1456

  6. [6]

    Journal of Plasma Physics , author =

    Phase mixing versus nonlinear advection in drift-kinetic plasma turbulence , volume =. Journal of Plasma Physics , author =. 2016 , pages =. doi:10.1017/S0022377816000374 , abstract =

    work page doi:10.1017/s0022377816000374 2016

  7. [7]

    Journal of Plasma Physics , author =

    Fluctuation-dissipation relations for a plasma-kinetic. Journal of Plasma Physics , author =. 2015 , pages =. doi:10.1017/S0022377814000622 , abstract =

    work page doi:10.1017/s0022377814000622 2015

  8. [8]

    Physics of Plasmas , author =

    Reduced fluid-kinetic equations for low-frequency dynamics, magnetic reconnection, and electron heating in low-beta plasmas , volume =. Physics of Plasmas , author =. 2011 , pages =. doi:10.1063/1.3628639 , abstract =

    work page doi:10.1063/1.3628639 2011

  9. [9]

    Journal of Plasma Physics , author =

    Laguerre–. Journal of Plasma Physics , author =. 2018 , keywords =. doi:10.1017/S0022377818000041 , abstract =

    work page doi:10.1017/s0022377818000041 2018

  10. [10]

    Plasma Physics and Controlled Fusion , author =

    Developments in the gyrofluid approach to. Plasma Physics and Controlled Fusion , author =. 1993 , pages =. doi:10.1088/0741-3335/35/8/006 , abstract =

    work page doi:10.1088/0741-3335/35/8/006 1993

  11. [11]

    Journal of Computational Physics , author =

    Asymptotic. Journal of Computational Physics , author =. 1984 , pages =. doi:10.1016/0021-9991(84)90124-4 , abstract =

    work page doi:10.1016/0021-9991(84)90124-4 1984

  12. [12]

    Joglekar, A. S. and Thomas, A. G. R. , month = mar, year =. On the evolution of a large-amplitude, weakly-collisional electron plasma wave , url =. doi:10.48550/arXiv.2603.11225 , abstract =

    work page doi:10.48550/arxiv.2603.11225

  13. [13]

    Su, C. H. and Oberman, C. , month = feb, year =. Collisional. Physical Review Letters , publisher =. doi:10.1103/PhysRevLett.20.427 , abstract =

    work page doi:10.1103/physrevlett.20.427

  14. [14]

    The Physics of Fluids , author =

    Dominant. The Physics of Fluids , author =. 1971 , pages =. doi:10.1063/1.1693397 , abstract =

    work page doi:10.1063/1.1693397 1971

  15. [15]

    The Physics of Fluids , author =

    Collisional damping of. The Physics of Fluids , author =. 1977 , pages =. doi:10.1063/1.861801 , abstract =

    work page doi:10.1063/1.861801 1977

  16. [16]

    The Physics of Fluids , author =

    Fourier‐. The Physics of Fluids , author =. 1967 , pages =. doi:10.1063/1.1762177 , abstract =

    work page doi:10.1063/1.1762177 1967

  17. [17]

    Physics of Plasmas , author =

    Kinetic simulation of steady states of ion temperature gradient driven turbulence with weak collisionality , volume =. Physics of Plasmas , author =. 2004 , pages =. doi:10.1063/1.1669393 , abstract =

    work page doi:10.1063/1.1669393 2004

  18. [18]

    Hatch, D. R. and Jenko, F. and Bañón Navarro, A. and Bratanov, V. , month = oct, year =. Transition between. Physical Review Letters , publisher =. doi:10.1103/PhysRevLett.111.175001 , abstract =

    work page doi:10.1103/physrevlett.111.175001

  19. [19]

    The Physics of Fluids , author =

    Numerical. The Physics of Fluids , author =. 1967 , pages =. doi:10.1063/1.1762272 , abstract =

    work page doi:10.1063/1.1762272 1967

  20. [20]

    Zakharov, V. E. and Karpman, V. I. , month = jan, year =. On the. Soviet Journal of Experimental and Theoretical Physics , publisher =

    work page

  21. [21]

    Sperling, J. L. , month = mar, year =. Electron. The Physics of Fluids , publisher =. doi:10.1063/1.862252 , abstract =

    work page doi:10.1063/1.862252

  22. [22]

    Townsend, Alex , year =. Fast

    work page

  23. [23]

    IMA Journal of Numerical Analysis , author =

    Fast computation of. IMA Journal of Numerical Analysis , author =. 2016 , pages =. doi:10.1093/imanum/drv002 , abstract =

    work page doi:10.1093/imanum/drv002 2016

  24. [24]

    Journal of Plasma Physics , author =

    Fourier–. Journal of Plasma Physics , author =. 2015 , pages =. doi:10.1017/S0022377814001287 , abstract =

    work page doi:10.1017/s0022377814001287 2015

  25. [25]

    and Ewart, Robert J

    Nastac, Michael L. and Ewart, Robert J. and Juno, James and Barnes, Michael and Schekochihin, Alexander A. , month = mar, year =. Universal fluctuation spectrum of. doi:10.48550/arXiv.2503.17278 , abstract =

    work page doi:10.48550/arxiv.2503.17278

  26. [26]

    Schekochihin, A. A. and Cowley, S. C. and Dorland, W. and Hammett, G. W. and Howes, G. G. and Quataert, E. and Tatsuno, T. , month = may, year =. Astrophysical gyrokinetics: kinetic and fluid turbulent cascades in magnetized weakly collisional plasmas , volume =. The Astrophysical Journal Supplement Series , publisher =. doi:10.1088/0067-0049/182/1/310 , ...

    work page doi:10.1088/0067-0049/182/1/310

  27. [27]

    The Physics of Fluids , author =

    Model. The Physics of Fluids , author =. 1964 , pages =. doi:10.1063/1.2746779 , abstract =

    work page doi:10.1063/1.2746779 1964

  28. [28]

    Physica , author =

    On the theory of stationary waves in plasmas , volume =. Physica , author =. 1955 , pages =. doi:10.1016/S0031-8914(55)93068-8 , abstract =

    work page doi:10.1016/s0031-8914(55)93068-8 1955

  29. [29]

    and Greene, John M

    Bernstein, Ira B. and Greene, John M. and Kruskal, Martin D. , month = nov, year =. Exact. Physical Review , publisher =. doi:10.1103/PhysRev.108.546 , abstract =

    work page doi:10.1103/physrev.108.546

  30. [30]

    Annals of Physics , author =

    Plasma oscillations , volume =. Annals of Physics , author =. 1959 , pages =. doi:10.1016/0003-4916(59)90029-6 , abstract =

    work page doi:10.1016/0003-4916(59)90029-6 1959