  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\`ere 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.
