The effect of the spatial nonlocality of the Kirkwood g-factor on the determination of the long wavelength dielectric functions in dipolar fluids.

The Kirkwood g-factor that determines the long wavelength dielectric constant of a simple, isotropic, translationally invariant dipolar fluid is given by an integral of a dipole-dipole correlation function over a spherical region of a nonzero radius R(K) chosen such that any further increase in the radius leads to no change in the value of the integral, thereby defining a Kirkwood correlation length R(K). For radii less than the correlation length the integral defines a radius dependent (nonlocal) Kirkwood g-factor, implying a nonlocal dielectric function. The nonlocal nature of these quantities has important consequences for the determination of the long wavelength dielectric function from dipole fluctuations via the Kirkwood-Fröhlich connection. The dipole-dipole correlation function (the volume dipole auto-correlation function) commonly used in this determination involves particles residing solely within a sphere of radius R, unlike the correct correlation function which involves either a single particle with those particles in a spherical volume of radius R(K) or those particles in a spherical volume of radius R with those residing within a spherical volume of radius R+R(K). A procedure is suggested for extracting the infinite system dipole-dipole correlation function from results of simulations performed on finite spherical samples. Using some results reported in the recent literature, relative to the accurate correlation function the commonly used correlation function ranges from 27% too small for a sphere having a radius comparable to the Kirkwood correlation length to 4% too small at a radius of seven times that correlation length. As a result, the apparent dielectric constants, as determined by the conventional procedure of using the fluctuations of the sum of dipoles in a finite fixed volume, are also too small. This suggests that a dielectric constant extracted from computer simulations using a total dipole-total dipole correlation function in a given volume with other geometries and/or boundary conditions will result in similar errors.