Search

The concept of guiding center is potent in astrophysics, space plasmas, fusion researches, and arc plasmas to solve the multi-scale dynamics of magnetized plasmas. In this letter, we rigorously prove that the guiding center dynamics can generally be described as a constrained canonical Hamiltonian system with two constraints in six dimensional phase space, and that the solution flow of the guiding center lies on a canonical symplectic sub-manifold. The guiding center can thus be modeled as a pseudo-particle with an intrinsic magnetic moment, which properly replaces the charged particle dynamics on time scales larger than the gyro-period. The complete dynamical behaviors, such as the velocity and force, of the guiding center pseudo-particle can be clearly deduced from the model. Furthermore, a series of related theories, such as symplectic numerical methods, the canonical gyro-kinetic theory, and canonical particle-in-cell algorithms can be systematically developed based on the canonical guiding center system. The canonical guiding center theory also provides an enlightenment for the origin of the intrinsic magnetic moment., Comment: 14pages, 4 figures

Gyrocenter dynamics of charged particles plays a fundamental role in plasma physics. In particular, accuracy and conservation of energy are important features for correctly performing long-time simulations. For this purpose, we here propose arbitrarily high-order energy conserving methods for its simulation. The analysis and the efficient implementation of the methods are fully described, and some numerical tests are reported., Comment: 23 pages, 4 figures

We prove that in finite dimensions, a Parity-Time (PT)-symmetric Hamiltonian is necessarily pseudo-Hermitian regardless of whether it is diagonalizable or not. This result is different from Mostafazadeh's, which requires the Hamiltonian to be diagonalizable. PT-symmetry breaking often occurs at exceptional points where the Hamiltonian is not diagonalizable. Our result implies that PT-symmetry breaking is equivalent to the onset of instabilities of pseudo-Hermitian systems, which was systematically studied by Krein et al. in 1950s. In particular, we show that the mechanism of PT-symmetry breaking is the resonance between eigenmodes with different Krein signatures., Comment: 11pages, 1 figure. arXiv admin note: text overlap with arXiv:1801.01676

Parity-Time (PT)-symmetry is being actively investigated as a fundamental property of observables in quantum physics. We show that the governing equations of the classical two-fluid interaction and the incompressible fluid system are PT-symmetric, and the well-known Kelvin-Helmholtz instability is the result of spontaneous PT-symmetry breaking. It is expected that all classical conservative systems governed by Newton's law admit PT-symmetry, and the spontaneous breaking thereof is a generic mechanism for classical instabilities. Discovering the PT-symmetry of systems in fluid dynamics and plasma physics and identifying the PT-symmetry breaking responsible for instabilities enable new techniques to classical physics and enrich the physics of PT-symmetry., Comment: 11 pages, 1 figure

It is generally believed that Parity-Time (PT)-symmetry breaking occurs when eigenvalues or both eigenvalues and eigenvectors coincide. However, we show that this well-accepted picture of PT-symmetry breaking is incorrect. Instead, we demonstrate that the physical mechanism of PT-symmetry breaking is the resonance between positive- and negative-action modes. It is proved that PT-symmetry breaking occurs when and only when this resonance condition is satisfied, and this mechanism applies to all known PT-symmetry breakings observed in different branches of physics. The result is achieved by proving a remarkable fact that in finite dimensions, a PT-symmetric Hamiltonian is necessarily pseudo-Hermitian, regardless whether it is diagonalizable or not., Comment: 15 pages, 3 figures

Three-wave instability is a fundamental process that has important applications in many branches of physics. It is widely accepted that the resonant condition $\omega_{z}\approx\omega_{x}+\omega_{y}$ for participating waves is the criteria for the onset of the instability. We show that this condition is neither sufficient nor necessary, instead, the exact criteria for the onset of the instability is that a positive-action mode resonates with a negative-action mode. This mechanism is imposed by the topology and geometry of the spectral space. Guided by this new theory, additional instability bands previously unknown are discovered., Comment: 14 pages, 4 figures

A lattice Maxwell system is developed with gauge-symmetry, symplectic structure and discrete space-time symmetry. Noether's theorem for Lie group symmetries is generalized to discrete symmetries for the lattice Maxwell system. As a result, the lattice Maxwell system is shown to admit a discrete local energy-momentum conservation law corresponding to the discrete space-time symmetry. These conservative properties make the discrete system an effective algorithm for numerically solving the governing differential equations on continuous space-time. Moreover, the lattice model, respecting all conservation laws and geometric structures, is as good as and probably more preferable than the continuous Maxwell model. Under the simulation hypothesis by Bostrom and in consistent with the discussion on lattice QCD by Beane et al., the two interpretations of physics laws on space-time lattice could be essentially the same., Comment: 12 pages

Structure-preserving geometric algorithm for the Vlasov-Maxwell (VM) equations is currently an active research topic. We show that spatially-discretized Hamiltonian systems for the VM equations admit a local energy conservation law in space-time. This is accomplished by proving that for a general spatially-discretized system, a global conservation law always implies a discrete local conservation law in space-time when the algorithm is local. This general result demonstrates that Hamiltonian discretizations can preserve local conservation laws, in addition to the symplectic structure, both of which are the intrinsic physical properties of infinite dimensional Hamiltonian systems in physics.

An infinite dimensional canonical symplectic structure and structure-preserving geometric algorithms are developed for the photon-matter interactions described by the Schr\"odinger-Maxwell equations. The algorithms preserve the symplectic structure of the system and the unitary nature of the wavefunctions, and bound the energy error of the simulation for all time-steps. This new numerical capability enables us to carry out first-principle based simulation study of important photon-matter interactions, such as the high harmonic generation and stabilization of ionization, with long-term accuracy and fidelity., Comment: 17 pages, 7 figures

Relativistic dynamics of a charged particle in time-dependent electromagnetic fields has theoretical significance and a wide range of applications. It is often multi-scale and requires accurate long-term numerical simulations using symplectic integrators. For modern large-scale particle simulations in complex, time-dependent electromagnetic field, explicit symplectic algorithms are much more preferable. In this paper, we treat the relativistic dynamics of a particle as a Hamiltonian system on the cotangent space of the space-time, and construct for the first time explicit symplectic algorithms for relativistic charged particles of order 2 and 3 using the sum-split technique and generating functions., Comment: 12 pages, 2 figures

We construct high order symmetric volume-preserving methods for the relativistic dynamics of a charged particle by the splitting technique with processing. Via expanding the phase space to include time $t$, we give a more general construction of volume-preserving methods that can be applied to systems with time-dependent electromagnetic fields. The newly derived methods provide numerical solutions with good accuracy and conservative properties over long time of simulation. Furthermore, because of the use of processing technique the high order methods are explicit, and cost less than the methods derived from standard compositions, thus are more efficient. The results are verified by the numerical experiments. Linear stability analysis of the methods show that the high order processed method allows larger time step size during integration.

An explicit high-order noncanonical symplectic algorithm for ideal two-fluid systems is developed. The fluid is discretized as particles in the Lagrangian description, while the electromagnetic fields and internal energy are treated as discrete differential form fields on a fixed mesh. With the assistance of Whitney interpolating forms, this scheme preserves the gauge symmetry of the electromagnetic field, and the pressure field is naturally derived from the discrete internal energy. The whole system is solved using the Hamiltonian splitting method discovered by He et al., which was been successfully adopted in constructing symplectic particle-in-cell schemes. Because of its structure preserving and explicit nature, this algorithm is especially suitable for large-scale simulations for physics problems that are multi-scale and require long-term fidelity and accuracy. The algorithm is verified via two tests: studies of the dispersion relation of waves in a two-fluid plasma system and the oscillating two-stream instability.

The two-stream instability is probably the most important elementary example of collective instabilities in plasma physics and beam-plasma systems. For a warm plasma with two charged particle species based on a 1D warm-fluid model, the instability diagram of the two-stream instability exhibits an interesting band structure that has not been explained. We show that the band structure for this instability is the consequence of the Hamiltonian nature of the warm two-fluid system. Interestingly, the Hamiltonian nature manifests as a complex G-Hamiltonian structure in wave-number space, which directly determines the instability diagram. Specifically, it is shown that the boundaries between the stable and unstable regions are locations for Krein collisions between eigenmodes with different Krein signatures. In terms of physics, this rigorously implies that the system is destabilized when a positive-action mode resonates with a negative-action mode, and that this is the only mechanism by which the system can be destabilized. It is anticipated that this physical mechanism of destabilization is valid for other collective instabilities in conservative systems in plasma physics, accelerator physics, and fluid dynamics systems, which admit infinite-dimensional Hamiltonian structures., Comment: 16 pages, 3 figures

Dynamics of a charged particle in the canonical coordinates is a Hamiltonian system, and the well-known symplectic algorithm has been regarded as the de facto method for numerical integration of Hamiltonian systems due to its long-term accuracy and fidelity. For long-term simulations with high efficiency, explicit symplectic algorithms are desirable. However, it is widely accepted that explicit symplectic algorithms are only available for sum-separable Hamiltonians, and that this restriction severely limits the application of explicit symplectic algorithms to charged particle dynamics. To overcome this difficulty, we combine the familiar sum-split method and a generating function method to construct second and third order explicit symplectic algorithms for dynamics of charged particle. The generating function method is designed to generate explicit symplectic algorithms for product-separable Hamiltonian with form of $H(\mathbf{p},\mathbf{q})=\mathbf{p}_{i}f(\mathbf{q})$ or $H(\mathbf{p},\mathbf{q})=\mathbf{q}_{i}f(\mathbf{p})$. Applied to the simulations of charged particle dynamics, the explicit symplectic algorithms based on generating functions demonstrate superiorities in conservation and efficiency., Comment: 17 pages, 3 figures

Explicit high-order non-canonical symplectic particle-in-cell algorithms for classical particle-field systems governed by the Vlasov-Maxwell equations are developed. The algorithm conserves a discrete non-canonical symplectic structure derived from the Lagrangian of the particle-field system, which is naturally discrete in particles. The electromagnetic field is spatially-discretized using the method of discrete exterior calculus with high-order interpolating differential forms for a cubic grid. The resulting time-domain Lagrangian assumes a non-canonical symplectic structure. It is also gauge invariant and conserves charge. The system is then solved using a splitting method discovered by He et al., which produces five exactly-soluable sub-systems, and high-order structure- preserving algorithms follow by combinations. The explicit, high-order, and conservative nature of the algorithms is especially suitable for long-term simulations of particle-field systems with extremely large number of degrees of freedom on massively parallel supercomputers. The algorithms have been tested and verified by the two physics problems, i.e., the nonlinear Landau damping and the electron Bernstein wave.

We show that, when applied to any non-canonical Hamiltonian system, any integrator that is symplectic for canonical Hamiltonian problems is actually conjugate symplectic for the non-canonical structure. This result is useful because it implies that canonically symplectic methods may be successfully applied to long-time integrations of non-canonical Hamiltonian problems, thus avoiding the need to construct ad hoc new methods. Numerical results for three non-canonical Hamiltonian systems demonstrate that (canonically) symplectic methods have significant advantages in numerical accuracy and near energy preservation over non-symplectic methods., Comment: This paper has been withdrawn by the author due to a error in the theorem

Hamiltonian integration methods for the Vlasov-Maxwell equations are developed by a Hamiltonian splitting technique. The Hamiltonian functional is split into five parts, i.e., the electrical energy, the magnetic energy, and the kinetic energy in three Cartesian components. Each of the subsystems is a Hamiltonian system with respect to the Morrison-Marsden-Weinstein Poisson bracket and can be solved exactly. Compositions of the exact solutions yield Poisson structure preserving, or Hamiltonian, integration methods for the Vlasov-Maxwell equations, which have superior long-term fidelity and accuracy.

The paper [1] by Crouseilles, Einkemmer, and Faou used an incorrect Poisson bracket for the Vlasov-Maxwell equations. If the correct Poisson bracket is used, the solution of one of the subsystems cannot be computed exactly in general. As a result, one cannot construct a symplectic scheme for the Vlasov-Maxwell equations using the splitting Hamiltonian method proposed in Ref [1]., Comment: 6 pages

Particle-in-Cell (PIC) simulation is the most important numerical tool in plasma physics. However, its long-term accuracy has not been established. To overcome this difficulty, we developed a canonical symplectic PIC method for the Vlasov-Maxwell system by discretizing its canonical Poisson bracket. A fast local algorithm to solve the symplectic implicit time advance is discovered without root searching or global matrix inversion, enabling applications of the proposed method to very large-scale plasma simulations with many, e.g., $10^{9}$, degrees of freedom. The long-term accuracy and fidelity of the algorithm enables us to numerically confirm Mouhot and Villani's theory and conjecture on nonlinear Landau damping over several orders of magnitude using the PIC method, and to calculate the nonlinear evolution of the reflectivity during the mode conversion process from extraordinary waves to Bernstein waves., Comment: Revised submission with more detailed derivation and a new application. 16 pages and 6 figures