Reaction Field Simulations of Monatomic and Diatomic Dipolar Fluids

Abstract

Firstly, the conditions are explored for which molecular dynamics simulations yield reliable thermodynamic results when the reaction field method is applied for dipolar fluids. The influence of the system size, of the scaling method for avoiding the total energy updrift, and of the reaction field dielectric constant are studied. With appropriate conditions, the results for the Stockmayer fluid when compared to previous results show excellent agreement up to very high dipole moments. In the next step the method is applied to two-centre Lennard-Jones plus ideal axial dipole (2CLJD) molecules of elongation L = 0.505. For one state point, the influence of the molecular shape as well as of the dipolar strength on energy and pressure is studied in detail. Finally, in the project of developing a physically based equation of state for polar fluids, thermodynamic results are presented for 2CLJD fluids of elongation L = 0.505 for four different dipole moments at 29 state points in a temperature and density grid. For these 116 runs which all were performed with vectorized codes on a CYBER 205 also the uncertainties in the results are given as have been estimated from the running averages.