Track new results across pure and applied mathematics.
Search papers, labs, and topics across Lattice.
The authors construct a broad family of one-step exponential Boris-type integrators with matrix-valued filter functions for charged-particle dynamics where fast gyration is unresolved by the time step. Addressing the catastrophic resonance-induced error blow-up and field-dependent degradation found in standard and filtered Boris schemes, the framework enables the design of uniformly bounded filters. Using discrete variation-of-constants and summation-by-parts analyses, they prove field-independent error bounds that allow arbitrary step sizes while incurring only a controlled order reduction near resonant frequencies.
Resonance blow-ups in strong-field particle dynamics can be completely eliminated across arbitrary step sizes by switching to uniformly bounded matrix filters.
We study numerical time integration for the motion of a charged particle in a strong magnetic field, focusing on the regime in which the fast gyration is not resolved by the time step. Standard Boris-type methods are structure preserving and accurate when the gyration is resolved, but their error constants deteriorate with increasing field strength; filtered Boris methods improve this behavior but may suffer from singularities and resonance-induced error blow-up. In this work we consider the case of a constant strong magnetic field and derive a broad family of one-step exponential integrators from the variation-of-constants formula, formulated in terms of matrix-valued filter functions. This framework includes existing filtered Boris variants and permits the construction of schemes with uniformly bounded filters. Using a discrete variation-of-constants representation together with summation-by-parts arguments, we establish error bounds whose constants are independent of the magnetic field strength and characterize the filter conditions required for first- and second-order accuracy. In particular, we show that bounded filters avoid resonance blow-up and allow arbitrary time steps, at the cost of a controlled order reduction near resonant frequencies. Numerical experiments confirm the predicted convergence behavior, the different accuracy of position and velocity components, and the role of the leading error terms across a wide range of step-size-field-strength products.