# Plasma and Fusion Research

## Volume 11, 1303003 (2016)

# Letters

- Graduate School of Engineering, Tottori University, Minami 4-101, Koyama-cho, Tottori-shi, Tottori 680-8552, Japan
- 1)Japan Atomic Energy Agency, 2-166 Omotedate, Obuchi, Rokkasho, Aomori 039-3212, Japan

### Abstract

A high-accuracy numerical integration algorithm for a charged particle motion is developed. The algorithm is based on the Hamiltonian mechanics and the operator decomposition. The algorithm is made to be time-reversal symmetric, and its order of accuracy can be increased to any order by using a recurrence formula. One of the advantages is that it is an explicit method. An effective way to decompose the time evolution operator is examined; the Poisson tensor is decomposed and non-canonical variables are adopted. The algorithm is extended to a time dependent fields' case by introducing the extended phase space. Numerical tests showing the performance of the algorithm are presented. One is the pure cyclotron motion for a long time period, and the other is a charged particle motion in a rapidly oscillating field.

### Keywords

numerical integration, charged particle motion, time reversal symmetry, high accuracy, recurrence formula, time dependent field, ponderomotive force

### Full Text

### References

- [1] M. Suzuki, Phys. Lett. A 146, 319 (1990).
- [2] M. Suzuki, Phys. Lett. A 165, 387 (1992).
- [3] H. Yoshida, Phys. Lett. A 150, 262 (1990).
- [4] M. Furukawa, Proceedings of PLASMA2014 (Niigata, Japan, Nov. 2014) 18PB-107.
- [5] Y. He, Y. Sun, J. Liu and H. Qin, J. Comput. Phys. 281, 135 (2015).
- [6] A.J. Lichtenberg and M.A. Lieberman, Regular and Chaotic Dynamics, Second Ed. (Springer, 1992) pp.15–16.
- [7] Y. Ohkawa, Batchelor thesis (Department of Applied Mathematics and Physics, Faculty of Engineering, Tottori University, 2014), in Japanese.
- [8] R.D. Hazeltine and F.L. Waelbroeck, The Framework of Plasma Physics (Perseus Books, 1998) pp.42–44.