pith. sign in

arxiv: 2605.17820 · v1 · pith:N27QU7XDnew · submitted 2026-05-18 · 🧮 math.NA · cs.NA

Solving Vlasov-Poisson system with an adaptive Hermite spectral method

Sihong Shao , Yanli Wang , Jie Wu This is my paper

Pith reviewed 2026-05-20 01:06 UTC · model grok-4.3

classification 🧮 math.NA cs.NA
keywords Vlasov-Poisson systemadaptive Hermite spectral methodfrequency indicatorconservative projectionfilamentationplasma simulationspectral method
0
0 comments X

The pith

An adaptive Hermite spectral method solves the Vlasov-Poisson system by dynamically adjusting resolution via a frequency indicator while preserving invariants through two-step projection.

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

The paper develops an adaptive Hermite spectral method for the Vlasov-Poisson system that describes collisionless plasma dynamics. A frequency indicator based on the size of high-order coefficients in the expansion detects when filamentation creates fine-scale structures that demand higher resolution. The method then changes the scaling factor in the symmetrically weighted Hermite basis and applies a two-step projection: a constrained optimization step enforces conservation of mass, momentum, energy, and the L2 norm, followed by an ODE approximation that updates the coefficients at linear cost. Experiments in one- and two-dimensional velocity settings demonstrate that the approach remains accurate and more efficient than non-adaptive alternatives. A reader would care because Vlasov-Poisson simulations routinely encounter filamentation that forces either prohibitive fixed resolution or loss of accuracy, and this technique offers a practical way to adapt locally without breaking conservation laws.

Core claim

The paper establishes that a frequency indicator computed from high-order Hermite coefficients can be used to adjust the scaling factor in the basis functions, and that a two-step projection operator—first formulated as a constrained optimization problem to hold the invariants fixed, then approximated by an ODE step—computes the updated coefficients efficiently while keeping mass, momentum, energy, and the L2 norm of the distribution function sufficiently accurate, thereby enabling the adaptive Hermite method to handle filamentation in the Vlasov-Poisson system, as confirmed by numerical tests in 1D1V and 2D2V regimes.

What carries the argument

The frequency indicator that quantifies the contribution of high-order expansion coefficients in the symmetrically weighted Hermite basis with scaling factor, together with the fast conservative projection operator that first solves a constrained optimization problem and then applies an ODE-based approximation to update coefficients.

If this is right

  • The adaptive adjustment allows the method to capture fine filamentary structures without committing to a uniformly high resolution from the outset.
  • Preservation of mass, momentum, energy, and L2 norm through the projection step supports stable integration over long times.
  • The linear complexity of the ODE approximation step keeps the cost of each adaptation modest even when changes occur frequently.
  • The approach extends at least to two-dimensional velocity spaces, as shown by the 2D2V experiments.

Where Pith is reading between the lines

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

  • Similar frequency-based indicators could be tested on other spectral bases for kinetic equations that develop filamentation.
  • Linking the indicator threshold to a physical scale such as the local Debye length might make adaptation more robust across different plasma regimes.
  • Direct comparison against a fixed high-order non-adaptive run on a problem with a known analytic solution would quantify the efficiency gain from adaptation.

Load-bearing premise

The frequency indicator based on high-order expansion coefficients reliably detects when resolution must increase due to filamentation, and the two-step projection preserves the invariants to sufficient accuracy for the target simulations.

What would settle it

A long-time simulation of the Vlasov-Poisson system in which strong filamentation develops yet the frequency indicator fails to increase the scaling factor, producing visible loss of accuracy or drift in the conserved quantities beyond the level tolerated in the reported tests.

Figures

Figures reproduced from arXiv: 2605.17820 by Jie Wu, Sihong Shao, Yanli Wang.

Figure 1
Figure 1. Figure 1: (Adaptive algorithm in Sec. 2.2) Illustration of the scaling adaptive adjustment procedures. When the frequency indicator exits the reference range, the scaling factor β is updated via a bidirectional line search on the discrete set S to minimize the indicator. Accordingly, the distribution function is approximated as f(t, x, v) ≈ X N k=0 ¯f β k (t, x)H¯ β k (v), ¯f β k (t, x) = Z R H¯ β k (v)f(t, x, v) ωβ… view at source ↗
Figure 2
Figure 2. Figure 2: (1D1V linear Landau damping in Sec. 5.1.1) Comparison between the scaling adaptive and non￾adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by the non-adaptive method (NA). (b) Evolution of the potential energy W E by the scaling adaptive method (SA). 5.1 1D1V examples We first consider several classical 1D1V benchmark problems, including linear Landau damping, … view at source ↗
Figure 3
Figure 3. Figure 3: (1D1V linear Landau damping in Sec. 5.1.1) (a) Evolution of the scaling factor β. (b) L 2 error of the scaling adaptive and non-adaptive methods with different expansion orders N at t = 20 and 30. up to approximately t ≈ 35, whereas the non-adaptive method exhibits recurrence at approximately t ≈ 18 for N = 64 and t ≈ 26 for N = 128. This demonstrates that the scaling adaptive algorithm achieves higher eff… view at source ↗
Figure 4
Figure 4. Figure 4: (1D1V linear Landau damping in Sec. 5.1.1) Evolution of the relative error for the total mass, energy, and L 2 norm with N = 64, and N = 128. Here, the solid lines are the numerical results with N = 128, while the dashed lines are those with N = 64. The blue lines are the error evolution of the total mass, the red lines are those of the total energy, while the purple lines are those of the L 2 norm. (a) Co… view at source ↗
Figure 5
Figure 5. Figure 5: (1D1V nonlinear Landau damping in Sec. 5.1.2) Comparison between the scaling adaptive and non-adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by the non-adaptive method. (b) Evolution of the potential energy W E by the scaling adaptive method. 0 5 10 15 20 25 30 35 40 t 0 1 2 3 4 5 6 7 8 9 10 - N = 128 N = 256 N = 512 N = 1024 (a) Scaling factor β 0 5 10 15 20 … view at source ↗
Figure 6
Figure 6. Figure 6: (1D1V nonlinear Landau damping in Sec. 5.1.2) Results of the scaling adaptive method with different expansion orders N. (a) Evolution of the scaling factor β. (b) Evolution of the errors in total mass, energy, and the L 2 norm, computed by the scaling adaptive method with N = 128. The computational time for the non-adaptive and scaling adaptive methods is listed in Tab. 3, where the time spent on obtaining… view at source ↗
Figure 7
Figure 7. Figure 7: (1D1V nonlinear Landau damping in Sec. 5.1.2) Comparison of the distribution function f between the scaling adaptive and non-adaptive method with different expansion orders. (a–d) t = 20. (e–h) t = 40 [PITH_FULL_IMAGE:figures/full_fig_p018_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: (1D1V two-stream instability in Sec. 5.1.3) Comparison between the scaling adaptive and non￾adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by both methods. (b) Evolution of the scaling factor β of the scaling adaptive method. (c) Evolution of the errors in total mass, energy, and the L 2 norm, computed by the scaling adaptive method with N = 128, ∆t = 0.001. T… view at source ↗
Figure 9
Figure 9. Figure 9: (1D1V two-stream instability in Sec. 5.1.3) Comparison of the distribution functions f between the scaling adaptive and non-adaptive method with different expansion orders N, t = 50. (a) Scaling adaptive method with N = 256. (b–d) Non-adaptive method. 0 5 10 15 20 25 30 t -6 -5 -4 -3 -2 -1 0 1 lo g p W E SA, N = 128 SA, N = 256 NA, N = 128 NA, N = 256 (a) WE 0 5 10 15 20 25 30 t 0 1 2 3 4 5 6 - N = 128 N =… view at source ↗
Figure 10
Figure 10. Figure 10: (1D1V bump-on-tail instability in Sec. 5.1.4) Comparison between the scaling adaptive and non￾adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by two methods. (b) Evolution of the scaling factor β. (c) Evolution of the errors in total mass, energy, and the L 2 norm, computed by the scaling adaptive method with N = 128 and ∆t = 0.00025. 20 [PITH_FULL_IMAGE:figu… view at source ↗
Figure 11
Figure 11. Figure 11: (1D1V bump-on-tail instability in Sec. 5.1.4) Comparison of the distribution functions f between the scaling adaptive and non-adaptive method with different expansion orders N, t = 30. (a–c) Scaling adaptive method. (d–f) Non-adaptive method. N = 128, 256, 512, and 1024 is presented in Fig. 10b, showing an increase over time as expected. The relative error for the mass, total energy, and the L 2 norm with… view at source ↗
Figure 12
Figure 12. Figure 12: (2D2V linear Landau damping in Sec. 5.2.1) Comparison between the scaling adaptive and non￾adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by the non-adaptive method. (b) Evolution of the potential energy W E by the scaling method. 5.2 2D2V examples In this subsection, three 2D2V numerical examples are studied, including the linear Landau damping, two￾stream i… view at source ↗
Figure 13
Figure 13. Figure 13: (2D2V linear Landau damping in Sec. 5.2.1) Results of the scaling adaptive method with different expansion orders N. (a) Evolution of the scaling factor β. (b) Evolution of the errors in total mass, energy, and the L 2 norm, computed by the scaling adaptive method with N1 = N2 = 64. problems, doubling the velocity expansion orders N1 and N2 leads to an approximately fourfold increase in computational cost… view at source ↗
Figure 14
Figure 14. Figure 14: (2D2V two-stream instability in Sec. 5.2.2) Comparison between the scaling adaptive and non￾adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by two methods. (b) Evolution of the scaling factor β. (c) Evolution of the relative errors in total mass, energy, and the L 2 norm, computed by the scaling adaptive method with N1 = N2 = 32. the anisotropic spreading of t… view at source ↗
Figure 15
Figure 15. Figure 15: (2D2V two-stream instability in Sec. 5.2.2) Comparison of the distribution f between the scal￾ing adaptive and non-adaptive method with different expansion orders N at t = 20. (a–d) Distribution f(v1, 0, x, L/2). (e–h) Distribution f(0, v2, L/2, y). accurately at moderate resolution. The evolution of the scaling factors is presented in Fig. 16b. Both β1 and β2 increase over time. A slight separation betwe… view at source ↗
Figure 16
Figure 16. Figure 16: (2D2V bump-on-tail instability in Sec. 5.2.3) Comparison between the scaling adaptive and non￾adaptive methods with different expansion orders N. (a) Evolution of the potential energy W E by two methods. (b) Evolution of the scaling factor β. (c) Evolution of the relative errors in total mass, energy, and the L 2 norm, computed by the scaling adaptive method with N1 = N2 = 32. 0 10 20 30 40 x -6 -4 -2 0 2… view at source ↗
Figure 17
Figure 17. Figure 17: (2D2V bump-on-tail instability in Sec. 5.2.3) Comparison of the distribution f between the scaling adaptive and non-adaptive method with different expansion orders N, t = 20. (a–c) Distribution f(v1, 0, x, L/2). (d–f) Distribution f(0, v2, L/2, y). 26 [PITH_FULL_IMAGE:figures/full_fig_p026_17.png] view at source ↗
read the original abstract

We propose an adaptive Hermite spectral method for the Vlasov-Poisson system based on a recently developed frequency indicator that measures the contribution of the high-order expansion coefficients. Precisely, the symmetrically weighted Hermite basis with a scaling factor is utilized to approximate the distribution function to satisfy the increasing resolution requirement, which, for example, is induced by filamentation. To implement the scaling adjustment, a fast conservative projection operator is constructed in two steps. The first step is to formulate the projection as a constrained optimization problem to preserve key invariants, including mass, momentum, energy, and the $L^2$ norm of the distribution function. The second step is an ODE-based approximation developed to compute the updated expansion coefficients with linear complexity. Numerical experiments with 1D1V and 2D2V settings validate the feasibility and efficiency of this proposed adaptive Hermite method.

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

2 major / 2 minor

Summary. The manuscript proposes an adaptive Hermite spectral method for the Vlasov-Poisson system. It uses a symmetrically weighted Hermite basis whose scaling factor is adjusted via a frequency indicator derived from the magnitude of high-order expansion coefficients, in order to resolve filamentation and other fine-scale structures. Scaling changes are implemented through a two-step projection: the update is first posed as a constrained optimization problem that enforces conservation of mass, momentum, energy, and the L² norm of the distribution function; this is then replaced by an ODE-based approximation that reduces the cost to linear in the number of modes. Numerical experiments in one- and two-dimensional velocity settings are reported to illustrate feasibility and efficiency.

Significance. If the method preserves the stated invariants to high accuracy under repeated adaptation and if the frequency indicator reliably detects the need for rescaling, the approach could offer a practical, structure-preserving spectral technique for long-time Vlasov-Poisson simulations at moderate cost. The explicit construction of a conservative projection and the linear-complexity ODE surrogate are concrete strengths that, if rigorously validated, would be of interest to the numerical plasma-physics community.

major comments (2)
  1. [§3.2] §3.2 (ODE approximation to the projection): the manuscript replaces the solution of the constrained optimization problem with an ODE integrator whose local truncation error is not bounded relative to the invariant constraints. Because adaptation may occur many times during a simulation, it is essential to show that the accumulated drift in mass, momentum, energy, and L² norm remains below a prescribed tolerance (e.g., 10^{-10}) over the full run; the current numerical tests do not report such long-term drift statistics.
  2. [§4] §4 (numerical experiments): the 1D1V and 2D2V results demonstrate qualitative agreement with reference solutions, yet no quantitative convergence study with respect to the number of Hermite modes or the adaptation threshold is provided. In particular, the observed order of accuracy after each rescaling step is not measured, leaving open whether the projection step degrades the spectral accuracy of the underlying Hermite expansion.
minor comments (2)
  1. [§2.3] The frequency indicator is defined in terms of the ratio of the sum of squared high-order coefficients to the total L² norm; a precise statement of the threshold value and its dependence on the scaling factor should be given explicitly rather than left as a tunable parameter.
  2. [Figures 3–6] Figure captions should state the number of modes, the adaptation frequency, and the final invariant errors for each plotted run.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the thorough review and valuable suggestions. We address the major comments point by point below and will incorporate the requested additions in the revised manuscript.

read point-by-point responses

  1. Referee: [§3.2] §3.2 (ODE approximation to the projection): the manuscript replaces the solution of the constrained optimization problem with an ODE integrator whose local truncation error is not bounded relative to the invariant constraints. Because adaptation may occur many times during a simulation, it is essential to show that the accumulated drift in mass, momentum, energy, and L² norm remains below a prescribed tolerance (e.g., 10^{-10}) over the full run; the current numerical tests do not report such long-term drift statistics.

    Authors: We agree that quantifying the accumulated drift in the invariants under repeated adaptations is important for demonstrating the robustness of the ODE-based projection. The present experiments illustrate short-term behavior and conservation properties during individual runs, but do not tabulate long-term drift statistics across many rescalings. In the revision we will add a dedicated numerical test that performs a large number of adaptations over an extended time interval and reports the maximum observed drift in mass, momentum, energy, and the L² norm, confirming that these remain below 10^{-10}. revision: yes

  2. Referee: [§4] §4 (numerical experiments): the 1D1V and 2D2V results demonstrate qualitative agreement with reference solutions, yet no quantitative convergence study with respect to the number of Hermite modes or the adaptation threshold is provided. In particular, the observed order of accuracy after each rescaling step is not measured, leaving open whether the projection step degrades the spectral accuracy of the underlying Hermite expansion.

    Authors: We acknowledge that the current numerical section emphasizes qualitative agreement and computational efficiency rather than systematic convergence rates. To address the concern about possible degradation of spectral accuracy by the projection, the revised manuscript will include quantitative studies: error norms versus number of Hermite modes for several fixed adaptation thresholds, together with measured convergence orders computed both immediately before and after rescaling events. These additions will clarify that the conservative projection preserves the expected spectral accuracy. revision: yes

Circularity Check

1 steps flagged

Minor self-citation on frequency indicator; central projection derivation remains independent

specific steps
  1. self citation load bearing [Abstract]
    "We propose an adaptive Hermite spectral method for the Vlasov-Poisson system based on a recently developed frequency indicator that measures the contribution of the high-order expansion coefficients."

    The adaptive framework is introduced as based on this indicator; if the indicator originates from prior overlapping-author work without independent re-derivation or external verification in the present manuscript, the central premise of reliable adaptation for filamentation inherits its justification from self-citation rather than standalone derivation.

full rationale

The paper constructs an adaptive Hermite method by combining a referenced frequency indicator with a newly formulated two-step projection (constrained optimization followed by ODE approximation) to handle scaling changes while preserving invariants. This construction does not reduce the claimed numerical validation or method feasibility to a tautological fit or self-referential definition. The indicator reference represents a single minor self-citation that is not load-bearing for the core contribution, as the projection operator and numerical experiments in 1D1V/2D2V settings provide independent content. No self-definitional, fitted-input, or ansatz-smuggling reductions are present.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The method relies on the standard properties of Hermite polynomials and the existence of a well-posed projection that preserves the listed invariants; no new physical entities are introduced. The frequency indicator and the ODE approximation step are algorithmic choices whose accuracy is asserted but not derived from first principles in the abstract.

free parameters (1)
  • scaling factor adjustment threshold
    The frequency indicator threshold that triggers rescaling is chosen to balance resolution and cost; its specific value is not derived and must be set for each problem.
axioms (2)
  • domain assumption The symmetrically weighted Hermite basis with scaling satisfies the increasing resolution requirement induced by filamentation.
    Invoked to justify the adaptive strategy for the Vlasov-Poisson distribution function.
  • standard math The constrained optimization problem admits a unique solution that exactly preserves mass, momentum, energy, and L2 norm.
    Used to formulate the projection operator.

pith-pipeline@v0.9.0 · 5682 in / 1424 out tokens · 40090 ms · 2026-05-20T01:06:13.699028+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

60 extracted references · 60 canonical work pages

  1. [1]

    J. W. Banks, A. G. Odu, R. Berger, T. Chapman, W. Arrighi, and S. Brunner. High-order accu- rate conservative finite difference methods for Vlasov equations in 2D+2V.SIAM J. Sci. Comput., 41(5):B953–B982, 2019. 29

    work page 2019

  2. [2]

    Birdsall and A

    C. Birdsall and A. Langdon.Plasma Physics via Computer Simulation. CRC Press, 2018

    work page 2018

  3. [3]

    Blaustein and F

    A. Blaustein and F. Filbet. A structure and asymptotic preserving scheme for the Vlasov–Poisson– Fokker–Planck model.J. Comput. Phys., 498:112693, 2024

    work page 2024

  4. [4]

    J. C. Butcher.Numerical Methods for Ordinary Differential Equations. Wiley, 2016

    work page 2016

  5. [5]

    Cai and R

    Z. Cai and R. Li. Numerical regularized moment method of arbitrary order for Boltzmann-BGK equa- tion.SIAM J. Sci. Comput., 32(5):2875–2907, 2010

    work page 2010

  6. [6]

    Z. Cai, B. Lin, and M. Lin. A positive and moment-preserving Fourier spectral method.SIAM J. Numer. Anal., 62(1):273–294, 2024

    work page 2024

  7. [7]

    Camporeale, G

    E. Camporeale, G. Delzanno, B. Bergen, and J. Moulton. On the velocity space discretization for the Vlasov–Poisson system: Comparison between implicit Hermite spectral and particle-in-cell methods. Comput. Phys. Commun., 198:47–58, 2016

    work page 2016

  8. [8]

    Camporeale, G

    E. Camporeale, G. L. Delzanno, G. Lapenta, and W. Daughton. New approach for the study of linear Vlasov stability of inhomogeneous systems.Phys. Plasmas, 13(9):092110, 2006

    work page 2006

  9. [9]

    J. Canosa. Numerical solution of Landau’s dispersion equation.J. Comput. Phys., 13(1):158–160, 1973

    work page 1973

  10. [10]

    Cassini and L

    F. Cassini and L. Einkemmer. Efficient 6D Vlasov simulation using the dynamical low-rank framework ensign.Comput. Phys. Commun., 280:108489, 2022

    work page 2022

  11. [11]

    Cheng and G

    C. Cheng and G. Knorr. The integration of the Vlasov equation in configuration space.J. Comput. Phys., 22(3):330–351, 1976

    work page 1976

  12. [12]

    T. Chou, S. Shao, and M. Xia. Adaptive Hermite spectral methods in unbounded domains.Appl. Numer. Math., 183:201–220, 2023

    work page 2023

  13. [13]

    Coughlin, J

    J. Coughlin, J. Hu, and U. Shumlak. Robust and conservative dynamical low-rank methods for the Vlasov equation via a novel macro-micro decomposition.J. Comput. Phys., 509:113055, 2024

    work page 2024

  14. [14]

    Degond, F

    P. Degond, F. Deluzet, L. Navoret, A.-B. Sun, and M.-H. Vignal. Asymptotic-Preserving Particle-In- Cell method for the Vlasov–Poisson system near quasineutrality.J. Comput. Phys., 229(16):5630–5652, 2010

    work page 2010

  15. [15]

    Delzanno

    G. Delzanno. Multi-dimensional, fully-implicit, spectral method for the Vlasov–Maxwell equations with exact conservation laws in discrete form.J. Comput. Phys., 301:338–356, 2015

    work page 2015

  16. [16]

    J. Denavit. Numerical simulation of plasmas with periodic smoothing in phase space.J. Comput. Phys., 9(1):75–98, 1972

    work page 1972

  17. [17]

    Einkemmer, K

    L. Einkemmer, K. Kormann, J. Kusch, R. G. McClarren, and J.-M. Qiu. A review of low-rank methods for time-dependent kinetic simulations.J. Comput. Phys., 538:114191, 2025

    work page 2025

  18. [18]

    Eliasson

    B. Eliasson. Outflow boundary conditions for the Fourier transformed three-dimensional Vlasov– Maxwell system.J. Comput. Phys., 225(2):1508–1532, 2007

    work page 2007

  19. [19]

    Filbet and T

    F. Filbet and T. Xiong. Conservative discontinuous Galerkin/Hermite spectral method for the Vlasov– Poisson system.Commun. Appl. Math. Comput., 4(1):34–59, 2022

    work page 2022

  20. [20]

    I. M. Gamba and S. H. Tharkabhushanam. Spectral-Lagrangian methods for collisional models of non-equilibrium statistical states.J. Comput. Phys., 228(6):2012–2036, 2009

    work page 2012

  21. [21]

    H. Grad. On the kinetic theory of rarefied gases.Commun. Pure Appl. Math., 2(4):331–407, 1949

    work page 1949

  22. [22]

    Guo and J.-M

    W. Guo and J.-M. Qiu. A local macroscopic conservative (LoMaC) low rank tensor method for the Vlasov dynamics.J. Sci. Comput., 101(3):61, 2024. 30

    work page 2024

  23. [23]

    J. P. Holloway. Spectral velocity discretizations for the Vlasov–Maxwell equations.Transport Theory Statist. Phys., 25(1):1–32, 1996

    work page 1996

  24. [24]

    Hu and H

    H. Hu and H. Yu. Scaling optimized Hermite approximation methods.SIAM J. Numer. Anal., 64(1):125– 147, 2026

    work page 2026

  25. [25]

    Z. Hu, Z. Cai, and Y. Wang. Numerical simulation of microflows using Hermite spectral methods.SIAM J. Sci. Comput., 42(1):B105–B134, 2020

    work page 2020

  26. [26]

    Huang and Z

    Y. Huang and Z. Yang. An analytically-solvable, wave-form asymptotic-preserving and energy- conserving time-splitting scheme for Vlasov–Poisson equations in the quasi-neutral regime.J. Comput. Phys., 544:114450, 2026

    work page 2026

  27. [27]

    Issan, O

    O. Issan, O. Chapurin, O. Koshkarov, and G. L. Delzanno. Effects of artificial collisions, filter- ing, and nonlocal closure approaches on Hermite-based Vlasov–Poisson simulations.Phys. Plasmas, 32(3):033906, 2025

    work page 2025

  28. [28]

    Issan, O

    O. Issan, O. Koshkarov, F. D. Halpern, B. Kramer, and G. L. Delzanno. Anti-symmetric and posi- tivity preserving formulation of a spectral method for Vlasov–Poisson equations.J. Comput. Phys., 514:113263, 2024

    work page 2024

  29. [29]

    Issan, O

    O. Issan, O. Koshkarov, F. D. Halpern, B. Kramer, and G. L. Delzanno. Conservative closures of the Vlasov–Poisson equations based on symmetrically weighted Hermite spectral expansion.J. Comput. Phys., 524:113741, 2025

    work page 2025

  30. [30]

    J. Juno, A. Hakim, J. TenBarge, E. Shi, and W. Dorland. Discontinuous Galerkin algorithms for fully kinetic plasmas.J. Comput. Phys., 353:110–147, 2018

    work page 2018

  31. [31]

    Kincaid and W

    D. Kincaid and W. Cheney.Numerical Analysis: Mathematics of Scientific Computing. American Mathematical Society, Providence, Rhode Island, third edition, 2009

    work page 2009

  32. [32]

    A. J. Klimas. A numerical method based on the Fourier–Fourier transform approach for modeling 1-D electron plasma evolution.J. Comput. Phys., 50(2):270–306, 1983

    work page 1983

  33. [33]

    A. J. Klimas. A method for overcoming the velocity space filamentation problem in collisionless plasma model solutions.J. Comput. Phys., 68(1):202–226, 1987

    work page 1987

  34. [34]

    Kormann, K

    K. Kormann, K. Reuter, and M. Rampp. A massively parallel semi-Lagrangian solver for the six- dimensional Vlasov–Poisson equation.Int. J. High Perform. Comput. Appl., 33(5):924–947, 2019

    work page 2019

  35. [35]

    Kormann and A

    K. Kormann and A. Yurova. A generalized Fourier–Hermite method for the Vlasov–Poisson system. BIT Numer. Math., 61(3):881–909, 2021

    work page 2021

  36. [36]

    H. Liu, X. Cai, Y. Cao, and G. Lapenta. An asymptotic-preserving conservative semi-Lagrangian scheme for the Vlasov–Maxwell system in the quasi-neutral limit.J. Comput. Phys., 528:113840, 2025

    work page 2025

  37. [37]

    Manzini, G

    G. Manzini, G. L. Delzanno, J. Vencels, and S. Markidis. A Legendre–Fourier spectral method with exact conservation laws for the Vlasov–Poisson system.J. Comput. Phys., 317:82–107, 2016

    work page 2016

  38. [38]

    Manzini, D

    G. Manzini, D. Funaro, and G. L. Delzanno. Convergence of spectral discretizations of the Vlasov– Poisson system.SIAM J. Numer. Anal., 55(5):2312–2335, 2017

    work page 2017

  39. [39]

    Moler and C

    C. Moler and C. V. Loan. Nineteen dubious ways to compute the exponential of a matrix.SIAM Rev., 20(4):801–836, 1978

    work page 1978

  40. [40]

    Mouhot and C

    C. Mouhot and C. Villani. On Landau damping.Acta Math., 207(1):29–201, 2011

    work page 2011

  41. [41]

    Pagliantini, G

    C. Pagliantini, G. L. Delzanno, and S. Markidis. Physics-based adaptivity of a spectral method for the Vlasov–Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space. J. Comput. Phys., 488:112252, 2023. 31

    work page 2023

  42. [42]

    Palmroth, U

    M. Palmroth, U. Ganse, Y. Pfau-Kempf, M. Battarbee, L. Turc, T. Brito, M. Grandin, S. Hoilijoki, A. Sandroos, and S. von Alfthan. Vlasov methods in space physics and astrophysics.Living Rev. Comput. Astrophys., 4(1), 2018

    work page 2018

  43. [43]

    Pareschi and T

    L. Pareschi and T. Rey. Moment preserving Fourier–Galerkin spectral methods and application to the Boltzmann equation.SIAM J. Numer. Anal., 60(6):3216–3240, 2022

    work page 2022

  44. [44]

    J. T. Parker and P. J. Dellar. Fourier–Hermite spectral representation for the Vlasov–Poisson system in the weakly collisional limit.J. Plasma Phys., 81(2), 2015

    work page 2015

  45. [45]

    C. Rampf. Cosmological Vlasov–Poisson equations for dark matter: Recent developments and connec- tions to selected plasma problems.Rev. Mod. Plasma Phys., 5(1), 2021

    work page 2021

  46. [46]

    Roytershteyn and G

    V. Roytershteyn and G. L. Delzanno. Spectral approach to plasma kinetic simulations based on Hermite decomposition in the velocity space.Front. Astron. Space Sci., 5:27, 2018

    work page 2018

  47. [47]

    J. W. Schumer and J. P. Holloway. Vlasov simulations using velocity-scaled Hermite representations. J. Comput. Phys., 144(2):626–661, 1998

    work page 1998

  48. [48]

    Serre.Matrices: Theory and Applications

    D. Serre.Matrices: Theory and Applications. Springer New York, 2010

    work page 2010

  49. [49]

    S. Shao, Y. Wang, and J. Wu. An adaptive Hermite spectral method for the Boltzmann equation, 2025. arXiv:2509.17981

    work page arXiv 2025

  50. [50]

    J. Shen, T. Tang, and L.-L. Wang.Spectral Methods: Algorithms, Analysis and Applications. Springer, Berlin, Heidelberg, 2011

    work page 2011

  51. [51]

    D. G. Swanson.Plasma kinetic theory. Series in Plasma Physics. CRC Press, Florida, 2008

    work page 2008

  52. [52]

    Tskhakaya, K

    D. Tskhakaya, K. Matyash, R. Schneider, and F. Taccogna. The Particle-In-Cell method.Contrib. Plasm. Phys., 47(8-9):563–594, 2007

    work page 2007

  53. [53]

    Vencels, G

    J. Vencels, G. L. Delzanno, A. Johnson, I. B. Peng, E. Laure, and S. Markidis. Spectral solver for multi- scale plasma physics simulations with dynamically adaptive number of moments.Procedia Comput. Sci., 51:1148–1157, 2015

    work page 2015

  54. [54]

    Vogman, U

    G. Vogman, U. Shumlak, and P. Colella. Conservative fourth-order finite-volume Vlasov–Poisson solver for axisymmetric plasmas in cylindrical (r, v r, vθ) phase space coordinates.J. Comput. Phys., 373:877– 899, 2018

    work page 2018

  55. [55]

    Wang and L

    H. Wang and L. Zhang. Convergence analysis of Hermite approximations for analytic functions.Math. Comput., 2025

    work page 2025

  56. [56]

    M. Xia, S. Shao, and T. Chou. Efficient scaling and moving techniques for spectral methods in un- bounded domains.SIAM J. Sci. Comput., 43(5):A3244–A3268, 2021

    work page 2021

  57. [57]

    M. Xia, S. Shao, and T. Chou. A frequency-dependentp-adaptive technique for spectral methods.J. Comput. Phys., 446:110627, 2021

    work page 2021

  58. [58]

    Xiong, G

    T. Xiong, G. Russo, and J.-M. Qiu. Conservative multi-dimensional semi-Lagrangian finite difference scheme: Stability and applications to the kinetic and fluid simulations.J. Sci. Comput., 79(2):1241– 1270, 2018

    work page 2018

  59. [59]

    S. Zaki, L. Gardner, and T. Boyd. A finite element code for the simulation of one-dimensional Vlasov plasmas. I. Theory.J. Comput. Phys., 79(1):184–199, 1988

    work page 1988

  60. [60]

    Zhang and I

    C. Zhang and I. M. Gamba. A conservative scheme for Vlasov–Poisson–Landau modeling collisional plasmas.J. Comput. Phys., 340:470–497, 2017. 32

    work page 2017