Adaptive timestepping Monte Carlo integration of Coulomb collisions
Here, we report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell–Jüttner statistics. The implementation is based on the Beliaev–Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the threedimensional particle momentum space and the fivedimensional guiding center phase space.
 Authors:

^{[1]}
;
^{[2]};
^{[1]}
 Aalto Univ., Espoo (Finland)
 Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
 Publication Date:
 Grant/Contract Number:
 259675; 298126; AC02–09CH11466
 Type:
 Accepted Manuscript
 Journal Name:
 Computer Physics Communications
 Additional Journal Information:
 Journal Volume: 222; Journal Issue: C; Journal ID: ISSN 00104655
 Publisher:
 Elsevier
 Research Org:
 Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
 Sponsoring Org:
 USDOE
 Country of Publication:
 United States
 Language:
 English
 Subject:
 72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS; 97 MATHEMATICS AND COMPUTING; Coulomb collision; Monte Carlo; Fokker–Planck equation; Milstein method
 OSTI Identifier:
 1441004
Sarkimaki, Konsta, Hirvijoki, E., and Terava, J.. Adaptive timestepping Monte Carlo integration of Coulomb collisions. United States: N. p.,
Web. doi:10.1016/j.cpc.2017.09.025.
Sarkimaki, Konsta, Hirvijoki, E., & Terava, J.. Adaptive timestepping Monte Carlo integration of Coulomb collisions. United States. doi:10.1016/j.cpc.2017.09.025.
Sarkimaki, Konsta, Hirvijoki, E., and Terava, J.. 2017.
"Adaptive timestepping Monte Carlo integration of Coulomb collisions". United States.
doi:10.1016/j.cpc.2017.09.025. https://www.osti.gov/servlets/purl/1441004.
@article{osti_1441004,
title = {Adaptive timestepping Monte Carlo integration of Coulomb collisions},
author = {Sarkimaki, Konsta and Hirvijoki, E. and Terava, J.},
abstractNote = {Here, we report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell–Jüttner statistics. The implementation is based on the Beliaev–Budker collision integral which allows both the test particle and the background plasma to be relativistic. The integration method supports adaptive time stepping, which is shown to greatly improve the computational efficiency. The Monte Carlo method is implemented for both the threedimensional particle momentum space and the fivedimensional guiding center phase space.},
doi = {10.1016/j.cpc.2017.09.025},
journal = {Computer Physics Communications},
number = C,
volume = 222,
place = {United States},
year = {2017},
month = {10}
}