Reflecting boundary conditions for classical molecular dynamics simulations of nonideal plasmas

The influence of boundary conditions on results of the classical molecular dynamics simulations of nonideal electron-ion plasma is analyzed. A comprehensive study is performed for the convergence of per-particle potential energy and pressure with the number of particles using both the nearest image method (periodic boundaries) and harmonic reflective boundaries. As a result an error caused by finiteness of the simulation box is estimated. Moreover the electron oscillations given by the spectra of the current autocorrelation function are analyzed for both types of the boundary conditions. Nonideal plasmas with the nonideality parameter in range 0.26-2.6 is considered. To speed up the classical molecular dynamics simulations the graphics accelerators code is used.