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.

This work opens a series of papers where we develop a general quasioptical theory for modeconverting electromagnetic beams in plasma and implement it in a numerical algorithm. Here, the basic theory is introduced. We consider a general quasimonochromatic multi-component wave in a weakly inhomogeneous linear medium with no sources. For any given dispersion operator that governs the wave field, we explicitly calculate the approximate operator that governs the wave envelope ψ to the second order in the geometrical-optics parameter. Then, we further simplify this envelope operator by assuming that the gradient of ψ transverse to the local group velocity is much larger than the corresponding parallel gradient. This leads to a parabolic differential equation for ψ ("quasioptical equation") in the basis of the geometrical-optics polarization vectors. Scalar and mode-converting vector beams are described on the same footing. We also explain how to apply this model to electromagnetic waves in general. In the next papers of this series, we report successful quasioptical modeling of radiofrequency wave beams in magnetized plasma based on this theory.

Inhomogeneous drift-wave turbulence can be modeled as an effective plasma where drift waves act as quantumlike particles and the zonal-flow velocity serves as a collective field through which they interact. This effective plasma can be described by a Wigner-Moyal equation (WME), which generalizes the quasilinear wave-kinetic equation (WKE) to the full-wave regime, i.e., resolves the wavelength scale. Unlike waves governed by manifestly quantumlike equations, whose WMEs can be borrowed from quantum mechanics and are commonly known, drift waves have Hamiltonians very different from those of conventional quantum particles. This causes unusual phase-space dynamics that is typically not captured by the WKE. We demonstrate how to correctly model this dynamics with the WME instead. Specifically, we report full-wave phase-space simulations of the zonal-flow formation (zonostrophic instability), deterioration (tertiary instability), and the so-called predator-prey oscillations. We also show how the WME facilitates analysis of these phenomena, namely, (i) we show that full-wave effects critically affect the zonostrophic instability, particularly its nonlinear stage and saturation; (ii) we derive the tertiary-instability growth rate; and (iii) we demonstrate that, with full-wave effects retained, the predator-prey oscillations do not require zonal-flow collisional damping, contrary to previous studies. We also show how the famous Rayleigh-Kuo criterion, which has been missing in wave-kinetic theories of drift-wave turbulence, emerges from the WME.

Newcomb's Lagrangian for ideal magnetohydrodynamics (MHD) in Lagrangian labeling is discretized using discrete exterior calculus. Variational integrators for ideal MHD are derived thereafter. Besides being symplectic and momentum-preserving, the schemes inherit built-in advection equations from Newcomb's formulation, and therefore avoid solving them and the accompanying error and dissipation. We implement the method in 2D and show that numerical reconnection does not take place when singular current sheets are present. We then apply it to studying the dynamics of the ideal coalescence instability with multiple islands. The relaxed equilibrium state with embedded current sheets is obtained numerically.