Sarkas: A fast pure-python molecular dynamics suite for plasma physics
We present an open-source, performant, pure-python molecular dynamics (MD) suite for non-ideal plasmas. The code, Sarkas, aims to accelerate the research process by providing an MD code complete with pre- and post-processing tools. Sarkas offers the ease of use of Python while employing the Numba library to obtain execution speeds comparable to that of compiled languages. The available tools in Sarkas include graphical displays of the equilibration process through a Jupyter interface and the ability to compute quantities such as, radial distribution functions, autocorrelation functions, and Green-Kubo relations. Many force laws used to simulate plasmas are included in Sarkas, namely, pure Coulomb, Yukawa and Molière pair-potentials. Sarkas also contains quantum statistical potentials and fast Ewald methods are included where necessary. An object-oriented approach allows for easy modification of Sarkas, such as adding new time integrators, boundary conditions and force laws.