The Optimum Maximum Allowed Displacement in Monte Carlo Simulation of One-Component Plasma

In this paper, a periodic One-Component Plasma (OCP) system of N-point particles is simulated by Monte Carlo (MC) technique in three dimensions. Because of the long range nature of the Coulomb potential, no cut-off distance is considered in calculations (i.e, for each particle i, the effect of the other N-1 particles on i, is taken into account). The maximum allowed displacement "dmax" used in MC simulation controls the convergence to the equilibrium state of the system. An optimum maximum allowed displacement, O-dmax, is found and is given by a function of the temperature and the density of the system. Obtaining this function is done statistically by fitting the calculated data resulted from simulating the system at different values of temperature and density. The obtained O-dmax gives fast convergence of the simulation.