THERMAL-ENERGY OF THE CRYSTALLINE ONE-COMPONENT PLASMA FROM DYNAMICAL SIMULATIONS

Molecular-dynamics simulations employing 1024 particles and a high-accuracy spline approximation of the Ewald potential have been used to measure the excess internal energy of the classical one-component plasma in both the fluid and solid (bcc) phases over a wide range (1 less-than-or-equal-to GAMMA less-than-or-equal-to 2000) of the coupling parameter. The energy data for GAMMA less-than-or-equal-to 300 are in excellent agreement with previous Monte Carlo calculations for this range, giving an independent corroboration of these calculations. Data in the crystalline regime 170 less-than-or-equal-to GAMMA less-than-or-equal-to 2000 are used to estimate the coefficients of the anharmonic energy component by least-squares fits. The first coefficient is in reasonable agreement with peturbation-theory predictions, although it is argued that determination of the coefficients in this manner is a rather ill-conditioned problem. The free-energy curves of the fluid and solid phases are found to cross at GAMMA(m) almost-equal-to 173, in agreement with recent estimates.