Related papers: Comment on "Symplectic integration of magnetic sys…
Time-centered, hence second-order, methods for integrating the relativistic momentum of charged particles in an electromagnetic field are derived. A new method is found by averaging the momentum before use in the magnetic rotation term, and…
Dependable numerical results from long-time simulations require stable numerical integration schemes. For Hamiltonian systems, this is achieved with symplectic integrators, which conserve the symplectic condition and exactly solve for the…
Symplectic integration algorithms have become popular in recent years in long-term orbital integrations because these algorithms enforce certain conservation laws that are intrinsic to Hamiltonian systems. For problems with large variations…
We construct a particle integrator for nonrelativistic particles by means of the splitting method based on the exact flow of the equation of motion of particles in the presence of constant electric and magnetic field. This integrator is…
A modification of the standard Boris algorithm, called filtered Boris algorithm, is proposed for the numerical integration of the equations of motion of charged particles in a strong non-uniform magnetic field in the asymptotic scaling…
The Lorentz equations describe the motion of electrically charged particles in electric and magnetic fields and are used widely in plasma physics. The most popular numerical algorithm for solving them is the Boris method, a variant of the…
The Boris algorithm, a closely related variational integrator and a newly proposed filtered variational integrator are studied when they are used to numerically integrate the equations of motion of a charged particle in a non-uniform strong…
Two documents have recently been posted on the arXiv describing a numerical integration algorithm: "symplectic orbit/spin tracking code for all-electric storage rings" [1] and some computational results therefrom [2]. This note comments…
A basic leapfrog integrator and its energy-preserving and variational / symplectic variants are proposed and studied for the numerical integration of the equations of motion of relativistic charged particles in an electromagnetic field. The…
Xiao and Qin [Computer Physics Comm., 265:107981, 2021] recently proposed a remarkably simple modification of the Boris algorithm to compute the guiding centre of the highly oscillatory motion of a charged particle with step sizes that are…
In order to perform numerical studies of long-term stability in nonlinear Hamiltonian systems, one needs a numerical integration algorithm which is symplectic. Further, this algorithm should be fast and accurate. In this paper, we propose…
Numerical algorithms based on variational and symplectic integrators exhibit special features that make them promising candidates for application to general relativity and other constrained Hamiltonian systems. This paper lays part of the…
An improved Boris algorithm for simulating the motion of charged particles in electromagnetic fields has been developed. This enhancement addresses the issue of inaccurate fast-scale cyclotron phase calculations present in the original…
For a separable Hamiltonian, there are two fundamental, time-symmetric, second-order velocity-Verlet (VV) and position-Verlet (PV) symplectic integrators. Similarly, there are two VV and PV version of exact energy conserving algorithms for…
The equations of motion of a single particle subject to an arbitrary electric and a static magnetic field form a Poisson system. We present a second-order time integration method which preserves well the Poisson structure and compare it to…
A class of explicit numerical schemes is developed to solve for the relativistic dynamics and spin of particles in electromagnetic fields, using the Lorentz-BMT equation formulated in the Clifford algebra representation of Baylis. It is…
This work gives a Lie operator derivation of various Boris solvers via a detailed study of trajectory errors in a constant magnetic field. These errors in the gyrocenter location and the gyroradius are the foundational basis for why Boris…
A new approach is developed to integrate numerically the equations of motion for systems of interacting rigid polyatomic molecules. With the aid of a leapfrog framework, we directly involve principal angular velocities into the integration,…
Hamiltonian systems of ordinary and partial differential equations are fundamental mathematical models spanning virtually all physical scales. A critical property for the robustness and stability of computational methods in such systems is…
The interaction of electrically charged particles with magnetic fields is a fundamental problem in several areas of physics. An example is the motion of energetic particles through a magnetized plasma. The most accurate and reliable way to…