Sarkas: A Fast Pure-Python Molecular Dynamics Suite for Plasma Physics
arXiv:2105.14955 · doi:10.1016/j.cpc.2021.108245
Abstract
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 but also 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.
References in corpus (11)
- Comments on PPPM, FMM, and the Ewald Method for Large Periodic Coulombic Systems
- How to mesh up Ewald sums (II): An accurate error estimate for the P3M algorithm
- Molecular Dynamics for Fermions
- Molecular Dynamics for Low Temperature Plasma-Surface Interaction Studies
- Kinetic modelling and molecular dynamics simulation of ultracold neutral plasmas including ionic correlations
- Crystallization of Carbon Oxygen Mixtures in White Dwarf Stars
- "The numerical accuracy of truncated Ewald sums for periodic systems with long-range Coulomb interactions"
- Temperature Relaxation in Hot Dense Hydrogen
- Efficacy of the Radial Pair Potential Approximation for Molecular Dynamics Simulations of Dense Plasmas
- Depinning dynamics of two-dimensional dusty plasmas on a one-dimensional periodic substrate
- Constrained Molecular Dynamics simulation of the Quark-Gluon Plasma