On the Implementation of Friedman Boundary Conditions in Liquid Water Simulations

We have applied the image approximation to the reaction field as suggested by H.L. Friedman [Mol. Phys., 29, 1533 (1975)] by investigating appropriate cavity sizes and system parameters for use in molecular simulations. The energy of and the structure around a central simple point charge (SPC) water molecule in a dielectric cavity was found to be in good agreement with the properties of a liquid sample. To confine the water molecules within the cavity, we introduced a short-range repulsion between a real charge and its image as the Lennard-Jones repulsive potential between oxygen atoms of the SPC potential. For a system of 65 water molecules a cavity radius of 10.45 Å is appropriate; this radius is altered to 12.00 Å for a cavity surrounding 113 molecules. The effect of the boundary is restricted to the outer-most water layer which is in contact with the dielectric continuum.