Slow manifolds of classical Pauli particle enable structure-preserving geometric algorithms for guiding center dynamics
- Univ. of Science and Technology of China, Hefei (China). Dept. of Plasma Physics and Fusion Technology
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Since variational symplectic integrators for the guiding center was proposed, structure-preserving geometric algorithms have become an active research field in plasma physics. We found that the slow manifolds of the classical Pauli particle enable a family of structure-preserving geometric algorithms for guiding center dynamics with long-term stability and accuracy. This discovery overcomes the difficulty associated with the unstable parasitic modes for variational symplectic integrators when applied to the degenerate guiding center Lagrangian. It is a pleasant surprise that Pauli's Hamiltonian for electrons, which predated the Dirac equation and marks the beginning of particle physics, reappears in classical physics as an effective algorithm for solving an important plasma physics problem. This technique is applicable to other degenerate Lagrangians reduced from regular Lagrangians.
- Research Organization:
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- National MC Energy R&D Program; National Key Research and Development Program; National Natural Science Foundation of China (NSFC); USDOE
- Grant/Contract Number:
- AC02-09CH11466; 2018YFE0304100; 2016YFA0400600; 2016YFA0400601; 2016YFA0400602; NSFC-11905220; 11805273
- OSTI ID:
- 1814587
- Alternate ID(s):
- OSTI ID: 1815301
- Journal Information:
- Computer Physics Communications, Vol. 265; ISSN 0010-4655
- Publisher:
- ElsevierCopyright Statement
- Country of Publication:
- United States
- Language:
- English
Similar Records
Discovering exact, gauge-invariant, local energy-momentum conservation laws for the electromagnetic gyrokinetic system by high-order field theory on heterogeneous manifolds
Gauge-symmetrization method for energy-momentum tensors in high-order electromagnetic field theories