An explicit, energy-conserving particle-in-cell scheme for relativistic plasmas

This paper extends an explicit, energy-conserving particle-in-cell scheme to relativistic plasmas by solving a local optimization problem that enforces exact energy conservation, demonstrating its compatibility with standard field solvers and superior performance on relativistic test problems despite rare instances of non-real solutions.

The Gist

Imagine you are trying to simulate a chaotic dance of billions of tiny, super-fast particles (like electrons) moving through space, interacting with invisible electric and magnetic fields. This is what scientists call a plasma. To do this on a computer, they use a method called Particle-in-Cell (PIC).

Think of the computer screen as a giant grid (like a chessboard). The particles are the pieces moving around, and the grid holds the map of the electric and magnetic forces.

The Problem: The "Leaky Bucket"

In traditional computer simulations, there's a major flaw. As the simulation runs, tiny mathematical errors pile up. It's like trying to carry water in a bucket with a slow, invisible leak. Over time, the water (energy) disappears from the bucket, or worse, the bucket starts filling up with water that wasn't there to begin with.

In physics simulations, this "leak" or "overflow" is called grid heating. It's a ghostly artifact where the computer thinks the plasma is getting hotter and more energetic just because of the math errors, not because of any real physical reason. This ruins the simulation, making it inaccurate.

The Solution: The "Perfect Balance"

The authors of this paper have developed a new, explicit (fast and straightforward) way to run these simulations that acts like a perfectly sealed bucket.

Here is how their new method works, using a simple analogy:

1. The Standard Step: Imagine you are pushing a shopping cart (a particle) through a store. You calculate where it should go next based on the current forces.
2. The Correction: In old methods, you just let the cart roll there. In this new method, after you calculate the new spot, you pause and ask: "Wait, did this move create or destroy any energy?"
3. The Optimization: If the answer is "yes," the computer performs a tiny, instant mathematical adjustment. It's like a very smart shopper who, realizing they spent a penny too much or too little, instantly adjusts their path by a microscopic amount to ensure the total cost (energy) remains exactly the same as it was before.
4. The Result: The simulation runs fast (it's "explicit," meaning it doesn't get bogged down in complex calculations), but it never loses or gains energy artificially.

The "Relativistic" Twist

The paper specifically tackles relativistic plasmas. This means the particles are moving so fast they are close to the speed of light. At these speeds, the rules of physics get weird (time slows down, mass seems to increase).

The authors took their "perfect balance" method, which was already good for slow-moving particles, and upgraded it to handle these super-fast, relativistic particles. They had to rewrite the math to account for these speed-of-light effects, but the core idea remains the same: force the energy to stay constant.

Does it work?

The authors tested their new method on four different "stress tests" involving high-speed particle beams and instabilities (chaotic behaviors).

• Accuracy: The new method predicted the behavior of the plasma just as well as the old, standard methods.
• Energy Conservation: This is the big win. While the old methods let energy drift by a noticeable amount over time, the new method kept the energy locked in with extreme precision (down to the level of tiny computer rounding errors).
• Rare Glitches: The math behind the "correction" step is so precise that, in extremely rare cases, it might suggest a mathematically impossible result (like an imaginary number). However, the authors found this happens so rarely (like finding a needle in a haystack) that it doesn't matter for practical use. They simply fix those few rare cases without breaking the simulation.

In Summary

This paper presents a new, faster, and more accurate way to simulate super-hot, fast-moving space plasma. It solves the age-old problem of simulations "leaking" energy by adding a smart, instant correction step that ensures the total energy of the system is perfectly preserved, all while running efficiently on modern computers.

Technical Summary

Technical Summary: An Explicit, Energy-Conserving Particle-in-Cell Scheme for Relativistic Plasmas

Problem Statement
Particle-in-cell (PIC) schemes are the standard for simulating kinetic plasma phenomena but historically suffer from a lack of discrete conservation properties. This deficiency leads to "grid heating," an artificial energy gain arising from finite grid instabilities. While fully implicit and semi-implicit schemes have been developed to conserve total energy exactly, they often incur high computational costs or require solving non-linear systems. Recent efforts have introduced explicit, energy-conserving schemes for non-relativistic plasmas, but extending these to the relativistic Vlasov-Maxwell system presents unique challenges. Relativistic effects are critical in applications such as tokamak runaway electrons, inertial confinement fusion, and astrophysical plasmas. Existing relativistic energy-conserving methods are either semi-implicit or fully explicit but rely on splitting schemes that require serial particle advancement within cells, limiting parallel scalability.

Methodology
The authors extend a recently developed explicit, energy-conserving PIC scheme (originally for non-relativistic systems) to the relativistic Vlasov-Maxwell system. The core methodology involves a standard time-integration step followed by a correction step at the end of each time-step to enforce exact energy conservation.

1. Particle Advance: The scheme utilizes the Higuera-Cary particle pusher, which preserves phase-space volume and accurately captures the $E \times B$ drift. The update for the proper velocity $u_p$ is nonlinearly implicit but analytically solvable, maintaining the scheme's effective explicit nature.
2. Field Solvers: The method is designed to be compatible with popular electromagnetic field solvers, specifically the Yee/Finite Difference Time Domain (FDTD) and the Pseudo-Spectral Analytic Time Domain (PSATD) methods. These solvers are chosen for their ability to accurately capture the light-wave dispersion relation, which is vital in relativistic applications.
3. Energy Correction via Optimization: The central innovation is a correction step where the updated proper velocity $u_p^{n+1}$ is determined by solving a constrained optimization problem.
   • Objective: Minimize the distance between the corrected velocity and a predicted velocity $u_p^\dagger$ (which is second-order accurate).
   • Constraint: Enforce exact conservation of total energy. The constraint equates the work done by the fields (derived from the discrete current and electric field) to the change in particle kinetic energy.
   • Solution: The optimization problem is analytically solvable. The solution scales the predicted velocity by a factor $\Gamma_p^n$, derived explicitly from the known velocities and Lorentz factors. This preserves the explicit nature of the scheme and ensures local, parallel scalability.
4. Handling Imaginary Solutions: In rare instances, the optimization constraint may yield an imaginary $\Gamma_p^n$. The authors note that for practical simulation parameters, these instances are extremely rare. When they occur, the scheme sets $\Gamma_p^n = 1$, which maintains second-order temporal accuracy and results in negligible energy errors (near machine precision).

Key Contributions

• Relativistic Extension: The paper successfully generalizes the explicit, energy-conserving PIC framework to the relativistic regime, addressing the non-trivial coupling between the Lorentz factor and velocity updates.
• Analytic Solvability: The authors demonstrate that the constrained optimization problem required for energy conservation admits an analytic solution, avoiding the need for iterative solvers or implicit time-stepping.
• Scalability: Unlike previous explicit relativistic energy-conserving schemes that require serial particle updates within cells, this method imposes no such requirement, offering improved parallel scalability.
• Solver Compatibility: The scheme is shown to be compatible with both FDTD and PSATD field solvers, with specific adaptations (such as time-averaging the electric field for PSATD) to maintain exact energy conservation.
• Accuracy Analysis: The paper provides a rigorous analysis showing the scheme is second-order accurate. It further clarifies that while the scaling of the correction factor $\Gamma_p^n$ differs slightly from the non-relativistic case (scaling as $O(\Delta t^3)$ rather than $O(\Delta t^4)$ for large velocities), this does not degrade the overall second-order accuracy of the scheme.

Numerical Results
The scheme is verified against four standard relativistic test problems:

1. Relativistic Two-Stream Instability: The scheme reproduces the theoretical relativistic growth rates and demonstrates significantly improved energy conservation compared to standard PIC, with negligible occurrences of imaginary $\Gamma$ values.
2. Relativistic Landau Damping: The method matches linear theory predictions for relativistic damping rates and achieves energy conservation improvements of approximately six orders of magnitude over standard PIC.
3. Relativistic Weibel Instability: The scheme correctly captures the linear growth and nonlinear saturation of the instability. It achieves double-precision energy conservation with zero problematic particles observed.
4. Filamentation Instability: In this challenging 2D test, the scheme reproduces theoretical growth rates and shows energy conservation improvements of up to seven orders of magnitude over standard methods, again with no problematic particles.

Significance and Claims
The authors claim that this work provides a robust, fully explicit, and energy-conserving alternative for relativistic plasma simulations. The primary significance lies in the ability to eliminate grid heating and mitigate finite grid instabilities without the computational overhead of implicit solvers or the scalability limitations of serial-splitting schemes. The paper modestly notes that while the optimization problem does not guarantee a real solution in all mathematical cases, such instances are sufficiently rare for practical parameters to permit "dramatic improvements in energy conservation." The method is presented as a direct extension of non-relativistic techniques, preserving the structure and efficiency of the original scheme while adapting it to the complexities of relativistic dynamics. Future work is suggested for application to more challenging problems and extension to relativistic collisional plasmas.