Search

We present an expanded Wolf formalism for direct summation of long-range dipole-dipole interactions and rule-of-thumbs how to choose optimal spherical cutoff (Rc) and damping parameter (α). This is done by comparing liquid radial distribution functions, dipole-dipole orientation correlations, particle energies, and dielectric constants, with Ewald sums and the Reaction field method. The resulting rule states that ασ < 1 and αRc > 3 for reduced densities around ρ(∗) = 1 where σ is the particle size. Being a pair potential, the presented approach scales linearly with system size and is applicable to simulations involving point dipoles such as the Stockmayer fluid and polarizable water models.

We present a new method for Monte Carlo or Molecular Dynamics numerical simulations of three-dimensional polar fluids. The simulation cell is defined to be the surface of the northern hemisphere of a four-dimensional (hyper)sphere. The point dipoles are constrained to remain tangent to the sphere and their interactions are derived from the basic laws of electrostatics in this geometry. The dipole-dipole potential has two singularities which correspond to the following boundary conditions: when a dipole leaves the northern hemisphere at some point of the equator, it reappears at the antipodal point bearing the same dipole moment. We derive all the formal expressions needed to obtain the thermodynamic and structural properties of a polar liquid at thermal equilibrium in actual numerical simulation. We notably establish the expression of the static dielectric constant of the fluid as well as the behavior of the pair correlation at large distances. We report and discuss the results of extensive numerical Monte Carlo simulations for two reference states of a fluid of dipolar hard spheres and compare these results with previous methods with a special emphasis on finite size effects.

A review of the literature on the calculation of electrostatic potentials, fields, and field gradients in systems consisting of charges and dipoles using the Ewald summation technique is presented. Discrepancies between the previous formulas are highlighted, and an error in the derivation of the reciprocal contributions to the electrostatic field and field gradient is corrected. The new formulas for the field and field gradient are shown to exhibit a termwise identity with the ones for the electrostatic energy.

The static dielectric properties of Stockmayer fluids are investigated in the hyperspherical geometry, S(3). Different methods of obtaining the static dielectric constant ε(r) are compared. Tested methods include the evaluation of the Kirkwood factor, fluctuations of the total dipole moment, and a two-center potential correlation formula to obtain the dielectric constant through effective interactions. With no coupling to the "surrounding," the different methods give consistent estimates of the dielectric constant. Adding a coupling to the surrounding gives large size dependencies and the two-center potential correlation formula breaks down. For low dipole moments, there is a good agreement in the dielectric constant with previous studies.