pith. sign in

arxiv: 2510.16542 · v2 · submitted 2025-10-18 · ⚛️ physics.comp-ph · physics.chem-ph· physics.plasm-ph

Extended phase-space symplectic integration for electron dynamics

Francois Mauger , Cristel Chandre This is my paper

Pith reviewed 2026-05-18 06:11 UTC · model grok-4.3

classification ⚛️ physics.comp-ph physics.chem-phphysics.plasm-ph
keywords symplectic integrationextended phase spacesplit-operator schemeselectron dynamicsplasma physicstime-dependent density functional theorynumerical stabilityaccuracy estimation
0
0 comments X

The pith

Extended phase space enables high-order symplectic split-operator integration for electron dynamics in plasma and TDDFT.

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

The paper shows how to introduce extended phase-space variables so that high-order symplectic split-operator schemes can integrate two classes of electron dynamics while preserving structure. One class is the classical motion of a charged particle in a uniform magnetic field plus turbulent electrostatic potential. The other is the time-dependent Kohn-Sham equations of density-functional theory. The authors supply the explicit extension steps, the stability condition required for the integrator, and a low-cost metric that monitors accuracy during a run. This framework is presented as a route to structure-preserving time stepping for both classical and quantum Hamiltonian systems that have finite or infinite degrees of freedom.

Core claim

By embedding the original electron dynamics in an extended phase space, the authors make it possible to apply high-order symplectic split-operator integrators to a one-and-a-half-degree-of-freedom plasma problem and to an infinite-degree-of-freedom TDDFT problem, derive the associated stability condition, and supply a simple on-the-fly accuracy metric that does not add significant computational cost.

What carries the argument

The extension procedure that adds auxiliary phase-space variables and later constrains them so the original dynamics are recovered while admitting a symplectic split-operator decomposition.

If this is right

  • High-order symplectic split-operator schemes become directly usable for long-time integration of the described plasma and TDDFT dynamics.
  • The derived stability condition sets the allowable time-step range for these integrators.
  • The inexpensive accuracy metric supplies a practical real-time diagnostic without extra force evaluations.
  • The same extension approach is claimed to apply to other classical and quantum Hamiltonian systems with finite or infinite degrees of freedom.

Where Pith is reading between the lines

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

  • Plasma turbulence simulations could run for longer times with less artificial dissipation of particle energy.
  • TDDFT propagations might maintain better conservation of total energy or other invariants over many optical cycles.
  • The stability condition could be checked on other perturbed Hamiltonian problems to test broader applicability.

Load-bearing premise

The added extended variables can be introduced and constrained without destroying the symplectic structure or injecting uncontrolled errors when the original system has one-and-a-half or infinitely many degrees of freedom.

What would settle it

A direct numerical check on either the plasma particle or the TDDFT system showing that the integrator violates symplecticity or that trajectories diverge uncontrollably once the extended variables are introduced and constrained.

Figures

Figures reproduced from arXiv: 2510.16542 by Cristel Chandre, Francois Mauger.

Figure 1
Figure 1. Figure 1: FIG. 1. Comparison of the error in energy given by Eq. (15) [PITH_FULL_IMAGE:figures/full_fig_p006_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Comparison of the error in energy given by Eq. (15) [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Comparison of the (a) accuracy and (b) efficacy of [PITH_FULL_IMAGE:figures/full_fig_p010_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Comparison of the (a) accuracy and (b) distance [PITH_FULL_IMAGE:figures/full_fig_p011_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Comparison of the (a) accuracy and (b) efficacy of [PITH_FULL_IMAGE:figures/full_fig_p014_6.png] view at source ↗
read the original abstract

We investigate the use of extended phase-space symplectic integration for simulating two different classes of electron dynamics. The first one, with one and a half degrees of freedom, comes from plasma physics and describes the classical dynamics of a charged particle in a strong, constant, and uniform magnetic field perturbed by a turbulent electrostatic potential. The second one, with an infinite number of degrees of freedom, comes from physical chemistry and corresponds to Kohn-Sham time-dependent density-functional theory. For both we lay out the extension procedure and stability condition for numerical integration of the dynamics using high-order symplectic split-operator schemes. We also identify a computationally inexpensive metric that can be used for on-the-fly estimation of the accuracy of simulations. Our work paves the way for broad application of symplectic split-operator integration of classical and quantum Hamiltonian systems with finite and infinite number of degrees of freedom by comparing different modes of implementation of extended phase space integration.

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 manuscript develops an extended phase-space formalism to enable high-order symplectic split-operator integration for electron dynamics in two settings: a classical 1.5-DOF plasma model (charged particle in uniform magnetic field plus turbulent electrostatic potential) and the infinite-DOF Kohn-Sham TDDFT system. It presents explicit extension procedures, derives associated stability conditions, and proposes a computationally inexpensive metric for on-the-fly accuracy monitoring.

Significance. If the stability conditions and metric are rigorously derived and shown to preserve the symplectic structure after constraint, the work would provide a practical route to long-time, structure-preserving simulations of Hamiltonian electron dynamics, which is valuable for both classical plasma problems and quantum TDDFT where artificial dissipation or energy drift must be avoided.

major comments (1)
  1. [TDDFT extension and stability section] The load-bearing claim for the infinite-DOF TDDFT case is that the phase-space extension can be introduced and subsequently constrained (or projected) while preserving the underlying symplectic form on the infinite-dimensional manifold and without generating uncontrolled truncation or constraint errors. The manuscript does not supply a derivation, machine-checked verification, or numerical test demonstrating that the proposed metric bounds deviation from symplecticity rather than merely a local residual; this assumption therefore remains unverified and directly affects the applicability statement for the TDDFT system.
minor comments (2)
  1. The abstract states that procedures and stability conditions are derived, yet the main text would benefit from explicit equation numbers or theorem statements that isolate the stability criterion for each system.
  2. Clarify whether the inexpensive metric is shown to be independent of the particular splitting order or whether it must be re-derived for each integrator.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of the manuscript and for highlighting the need for stronger justification of the infinite-dimensional TDDFT case. We address the single major comment below and will incorporate the requested clarifications in a revised version.

read point-by-point responses

  1. Referee: [TDDFT extension and stability section] The load-bearing claim for the infinite-DOF TDDFT case is that the phase-space extension can be introduced and subsequently constrained (or projected) while preserving the underlying symplectic form on the infinite-dimensional manifold and without generating uncontrolled truncation or constraint errors. The manuscript does not supply a derivation, machine-checked verification, or numerical test demonstrating that the proposed metric bounds deviation from symplecticity rather than merely a local residual; this assumption therefore remains unverified and directly affects the applicability statement for the TDDFT system.

    Authors: We agree that the manuscript would benefit from a more explicit derivation of symplectic preservation under the phase-space extension and subsequent constraint for the infinite-dimensional Kohn-Sham TDDFT system. The current text presents the extension procedure and stability condition by direct analogy with the finite-dimensional plasma model and identifies the inexpensive metric, but does not contain a self-contained proof that the projection step preserves the symplectic form on the infinite-dimensional manifold or a dedicated numerical demonstration that the metric bounds global deviation from symplecticity. In the revision we will add a dedicated subsection that (i) derives the preservation of the symplectic structure after constraint for the TDDFT case, (ii) states the resulting stability condition, and (iii) includes a numerical test on a representative infinite-dimensional model problem showing that the metric tracks deviation from symplecticity rather than only a local residual. This will directly support the applicability statement. revision: yes

Circularity Check

0 steps flagged

No circularity: extension procedure and stability condition derived from standard symplectic split-operator methods without reduction to inputs by construction.

full rationale

The paper explicitly lays out an extension procedure, stability condition, and inexpensive accuracy metric for symplectic integration applied to both the 1.5-DOF plasma system and infinite-DOF Kohn-Sham TDDFT. No quoted equations or claims show a result defined in terms of itself, a fitted parameter renamed as prediction, or a load-bearing step that collapses to a self-citation chain. The derivation chain remains self-contained against external benchmarks of symplectic integrators, with the central claims consisting of methodological extensions rather than tautological redefinitions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The approach rests on standard properties of Hamiltonian flows and symplectic integrators; no new free parameters or invented entities are introduced in the abstract.

axioms (1)
  • standard math Symplectic integrators preserve phase-space volume for Hamiltonian systems
    Invoked when claiming that the split-operator schemes remain symplectic after phase-space extension.

pith-pipeline@v0.9.0 · 5682 in / 1124 out tokens · 25011 ms · 2026-05-18T06:11:37.451975+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

35 extracted references · 35 canonical work pages

  1. [1]

    P. J. Morrison, Rev. Mod. Phys.70, 467 (1998)

    work page 1998

  2. [2]

    Mauger, C

    F. Mauger, C. Chandre, and T. Uzer, Phys. Rev. Lett. 102, 173002 (2009)

    work page 2009

  3. [3]

    Dubois, S

    J. Dubois, S. A. Berman, C. Chandre, and T. Uzer, Phys. Rev. Lett.121, 113202 (2018)

    work page 2018

  4. [4]

    N. W. Ashcroft and N. D. Mermin,Solid State Physics (Holt, Rinehart and Winston, New York, 1976)

    work page 1976

  5. [5]

    Altland and B

    A. Altland and B. D. Simons,Condensed Matter Field Theory, 2nd ed. (Cambridge University Press, 2010)

    work page 2010

  6. [6]

    W. H. Miller, Faraday Discussions110, 1 (1998)

    work page 1998

  7. [7]

    G´ omez Pueyo, M

    A. G´ omez Pueyo, M. A. Marques, A. Rubio, and A. Cas- tro, J. Chem. Theory Comput.14, 3040 (2018)

    work page 2018

  8. [8]

    Mauger, C

    F. Mauger, C. Chandre, M. Gaarde, K. Lopata, and K. Schafer, Communications in Nonlinear Science and Numerical Simulation129, 107685 (2024)

    work page 2024

  9. [9]

    Tabor,Chaos and Integrability in Nonlinear Dynam- ics: An Introduction(John Wiley & Sons, New York, 1989)

    M. Tabor,Chaos and Integrability in Nonlinear Dynam- ics: An Introduction(John Wiley & Sons, New York, 1989)

    work page 1989

  10. [10]

    Masoliver and A

    J. Masoliver and A. Ros, European Journal of Physics 32, 431 (2011). 14 FIG. 6. Comparison of the (a) accuracy and (b) efficacy of extended phase-space symplectic integration using different split-operator schemes – see legend and main text. 2S corre- sponds to the 2 nd-order Strang (a.k.a.Verlet) scheme, 4FR to the 4 th-order Forest-Ruth, and 4/6BM to th...

    work page 2011

  11. [11]

    V. I. Arnold,Mathematical Methods of Classical Mechan- ics, 2nd ed., Graduate Texts in Mathematics, Vol. 60 (Springer, New York, 1989)

    work page 1989

  12. [12]

    Hairer, C

    E. Hairer, C. Lubich, and G. Wanner,Geometric Nu- merical Integration: Structure-Preserving Algorithms for Ordinary Differential Equations, 2nd ed., Springer Series in Computational Mathematics (Springer, Berlin, 2006)

    work page 2006

  13. [13]

    J. E. Marsden and T. S. Ratiu,Introduction to Mechanics and Symmetry, 2nd ed., Texts in Applied Mathematics, Vol. 17 (Springer, New York, 1999)

    work page 1999

  14. [14]

    R. I. McLachlan, Commun. Comput. Phys.31, 987 (2022)

    work page 2022

  15. [15]

    Yoshida, Phys

    H. Yoshida, Phys. Lett. A150, 262 (1990)

    work page 1990

  16. [16]

    R. I. McLachlan and P. Atela, Nonlinearity5, 541 (1992)

    work page 1992

  17. [17]

    Kraus, K

    M. Kraus, K. Kormann, P. J. Morrison, and E. Son- nendr¨ ucker, J. Plasma Phys.83, 905830401 (2017)

    work page 2017

  18. [18]

    Laskar and P

    J. Laskar and P. Robutel, Celestial Mechanics and Dy- namical Astronomy80, 39 (2001)

    work page 2001

  19. [19]

    A. C. Petit, J. Laskar, G. Bou´ e, and M. Gastineau, A&A 628, A32 (2019)

    work page 2019

  20. [20]

    Blanes, F

    S. Blanes, F. Casas, A. Farr´ es, J. Laskar, J. Makazaga, and A. Murua, Applied Numerical Mathematics68, 58 (2013)

    work page 2013

  21. [21]

    Pihajoki, Celest Mech Dyn Astr121, 211 (2015)

    P. Pihajoki, Celest Mech Dyn Astr121, 211 (2015)

    work page 2015

  22. [22]

    Tao, Phys

    M. Tao, Phys. Rev. E94, 043303 (2016)

    work page 2016

  23. [23]

    T. F. Viscondi, I. L. Caldas, and P. J. Morrison, Physics of Plasmas24, 032102 (2017)

    work page 2017

  24. [24]

    Kohn and L

    W. Kohn and L. Sham, Phys. Rev.140, A1133 (1965)

    work page 1965

  25. [25]

    Runge and E

    E. Runge and E. K. U. Gross, Phys. Rev. Lett.52, 997 (1984)

    work page 1984

  26. [26]

    Python codes for theE×Bmodel:http://github.com/ cchandre/GC2D

    work page

  27. [27]

    Mauger and C

    F. Mauger and C. Chandre, SoftwareX28, 101968 (2024). [28]QMol-gridpackage repositories: GitHubhttps: //github.com/fmauger1/QMol-grid. MathWorkhttps: //www.mathworks.com/matlabcentral/fileexchange/ 171414-qmol-grid

    work page 2024

  28. [28]

    Moler and C

    C. Moler and C. Van Loan, SIAM Review45, 3 (2003)

    work page 2003

  29. [29]

    Blanes and P

    S. Blanes and P. Moan, J. Comp. Appl. Math.142, 313 (2002)

    work page 2002

  30. [30]

    J. R. Cary and A. J. Brizard, Reviews of Modern Physics 81, 693 (2009)

    work page 2009

  31. [31]

    Fass` o, Journal of Computational Physics189, 527 (2003)

    F. Fass` o, Journal of Computational Physics189, 527 (2003)

    work page 2003

  32. [32]

    Strang, SIAM J

    G. Strang, SIAM J. Numer. Analysis5, 506 (1968)

    work page 1968

  33. [33]

    Forest and R

    E. Forest and R. D. Ruth, Physica D43, 105 (1990)

    work page 1990

  34. [34]

    R. I. McLachlan, SIAM Journal on Scientific Computing 16, 151 (1995)

    work page 1995

  35. [35]

    Omelyan, I

    I. Omelyan, I. Mryglod, and R. Folk, Computer Physics Communications146, 188 (2002). [37]pyhamsyspython package:https://pypi.org/project/ pyhamsys/. GitHub:https://github.com/cchandre/ pyhamsys

    work page 2002