10,079 research outputs found
Ion Sizes and Finite-Size Corrections for Ionic-Solvation Free Energies
Free energies of ionic solvation calculated from computer simulations exhibit
a strong system size dependence. We perform a finite-size analysis based on a
dielectric-continuum model with periodic boundary conditions. That analysis
results in an estimate of the Born ion size. Remarkably, the finite-size
correction applies to systems with only eight water molecules hydrating a
sodium ion and results in an estimate of the Born radius of sodium that agrees
with the experimental value.Comment: 2 EPS figure
Ion Pair Potentials-of-Mean-Force in Water
Recent molecular simulation and integral equation results alkali-halide ion
pair potentials-of-mean-force in water are discussed. Dielectric model
calculations are implemented to check that these models produce that
characteristic structure of contact and solvent-separated minima for oppositely
charged ions in water under physiological thermodynamic conditions. Comparison
of the dielectric model results with the most current molecular level
information indicates that the dielectric model does not, however, provide an
accurate description of these potentials-of-mean-force. We note that linear
dielectric models correspond to modelistic implementations of second-order
thermodynamic perturbation theory for the excess chemical potential of a
distinguished solute molecule. Therefore, the molecular theory corresponding to
the dielectric models is second-order thermodynamic perturbation theory for
that excess chemical potential. The second-order, or fluctuation, term raises a
technical computational issue of treatment of long-ranged interactions similar
to the one which arises in calculation of the dielectric constant of the
solvent. It is contended that the most important step for further development
of dielectric models would be a separate assessment of the first-order
perturbative term (equivalently the {\it potential at zero charge} ) which
vanishes in the dielectric models but is generally nonzero. Parameterization of
radii and molecular volumes should then be based of the second-order
perturbative term alone. Illustrative initial calculations are presented and
discussed.Comment: 37 pages and 8 figures. LA-UR-93-420
Balancing Local Order and Long-Ranged Interactions in the Molecular Theory of Liquid Water
A molecular theory of liquid water is identified and studied on the basis of
computer simulation of the TIP3P model of liquid water. This theory would be
exact for models of liquid water in which the intermolecular interactions
vanish outside a finite spatial range, and therefore provides a precise
analysis tool for investigating the effects of longer-ranged intermolecular
interactions. We show how local order can be introduced through quasi-chemical
theory. Long-ranged interactions are characterized generally by a conditional
distribution of binding energies, and this formulation is interpreted as a
regularization of the primitive statistical thermodynamic problem. These
binding-energy distributions for liquid water are observed to be unimodal. The
gaussian approximation proposed is remarkably successful in predicting the
Gibbs free energy and the molar entropy of liquid water, as judged by
comparison with numerically exact results. The remaining discrepancies are
subtle quantitative problems that do have significant consequences for the
thermodynamic properties that distinguish water from many other liquids. The
basic subtlety of liquid water is found then in the competition of several
effects which must be quantitatively balanced for realistic results.Comment: 8 pages, 6 figure
Quasi-chemical theory with a soft cutoff
In view of the wide success of molecular quasi-chemical theory of liquids,
this paper develops the soft-cutoff version of that theory. This development
has important practical consequences in the common cases that the packing
contribution dominates the solvation free energy of realistically-modeled
molecules because treatment of hard-core interactions usually requires special
purpose simulation methods. In contrast, treatment of smooth repulsive
interactions is typically straightforward on the basis of widely available
software. This development also shows how fluids composed of molecules with
smooth repulsive interactions can be treated analogously to the molecular-field
theory of the hard-sphere fluid. In the treatment of liquid water,
quasi-chemical theory with soft-cutoff conditioning doesn't change the
fundamental convergence characteristics of the theory using hard-cutoff
conditioning. In fact, hard cutoffs are found here to work better than softer
ones.Comment: 5 pages, 2 figure
- …