Plasma and Fusion Research

Volume 16, 1203100 (2021)

Rapid Communications

Double Leap-Frog Method for Large-Time-Step Particle Simulation to Keep Larmor Radius Small
Tomonori TAKIZUKA, Kenzo IBANO and Satoshi TOGO1)
Graduate School of Engineering, Osaka University, Suita 565-0871, Japan
Plasma Research Center, University of Tsukuba, Tsukuba 305-8577, Japan
(Received 6 September 2021 / Accepted 13 October 2021 / Published 1 November 2021)


A modified leap-frog (LF) scheme is presented that keeps the correct Larmor radius even in case of a large time step Δt compared to the cyclotron period Ω−1, ΩΔt ≫ 1, for the particle simulation of a plasma in the strong magnetic field. The Larmor radius simulated by the conventional LF method becomes very large for ΩΔt ≫ 1, and such a numerical condition has been avoided in general. If the LF method is applicable to such situations, new particle simulation codes can be more easily developed for a wide area of plasma physics. By repeating the LF steps doubly and adopting the averaged velocity to advance the particle position, the Larmor radius is kept real independently of the ΩΔt value. Proper nature on the energy conservation, magnetic moment conservation and drift-velocity realization is safely inherited from the LF method.


particle simulation, leap-frog method, time step, Larmor radius

DOI: 10.1585/pfr.16.1203100


  • [1] J.P. Verboncoeur, Plasma Phys. Control. Fusion 47, A231 (2005).
  • [2] J.P. Boris, Proc. 4th Conf. on Numerical Simulation of Plasmas, p.3 (1970).
  • [3] S. Zenitani and T. Umeda, Phys. Plasmas 25, 112110 (2018).
  • [4] T. Takizuka, Plasma Phys. Control. Fusion 59, 034008 (2017).
  • [5] T. Takizuka and H. Abe, J. Comput. Phys. 25, 205 (1977).
  • [6] T. Takizuka, K. Ibano and M. Yagi, Plasma Fusion Res. 13, 1203088 (2018).
  • [7] T. Takizuka, K. Ibano and M. Yagi, Plasma Fusion Res. 14, 1203091 (2019).
  • [8] K. Ibano, Y. Kikuchi, S. Togo, Y. Ueda and T. Takizuka, Nucl. Fusion 59, 076001 (2019).
  • [9] X. Garbet, Y. Idomura, L. Villard and T.H. Watanabe, Nucl. Fusion 50, 043002 (2010).