Multiple Boris integrators for particle-in-cell simulation

The particle-in-cell (PIC) simulation has been extensively used in space plasma physics. In PIC simulation, the Boris method (1970) is a standard method to advance charged particles, due to its accuracy, robustness, and stability. Meanwhile, there is a growing demand for better particle solvers, because numerical errors might be accumulated in long-term PIC simulations. In this contribution, we propose a new family of Boris-type schemes for integrating gyration of charged particles in PIC simulation. The new solvers combine the popular Boris solver arbitrary n times, and therefore we call them the multiple Boris solvers. Using Chebyshev polynomials, a one-step numerical procedure is provided. The new solvers give n**2 times smaller errors and allow larger timesteps. We will also present benchmark results of the multiple Boris solvers and other particle solvers.