Ion Pair Potentials-of-Mean-Force in Water

We demonstrate that with two small modifications, the popular dielectric continuum model is capable of predicting, with high accuracy, ion solvation thermodynamics in numerous polar solvents, and ion solvation free energies in water--co-solvent mixtures. The first modification involves perturbing the macroscopic dielectric-flux interface condition at the solute--solvent interface with a nonlinear function of the local electric field, giving what we have called a solvation-lay...

We formulate low-frequency charge-fluctuation forces between charged cylinders - parallel or skewed - in salt solution: forces from dipolar van der Waals fluctuations and those from the correlated monopolar fluctuations of mobile ions. At high salt concentrations forces are exponentially screened. In low-salt solutions dipolar energies go as $R^{-5}$ or $R^{-4}$; monopolar energies vary as $R^{-1}$ or $\ln{R}$, where $R$ is the minimal separation between cylinders. However, p...

We argue that the kosmotropes remain strongly hydrated in the vicinity of a hydrophobic surface, while the chaotropes lose their hydration shell and can become adsorbed to the interface. The mechanism of adsorption is still a subject of debate. We argue that there are two driving forces for anionic adsorption: the hydrophobic cavitational energy and the interfacial electrostatic surface potential of water. While the cavitational contribution to ionic adsorption is now well ac...

The AGBNP2 implicit solvent model, an evolution of the Analytical Generalized Born plus Non-Polar (AGBNP) model we have previously reported, is presented with the aim of modeling hydration effects beyond those described by conventional continuum dielectric representations. A new empirical hydration free energy component based on a procedure to locate and score hydration sites on the solute surface is introduced to model first solvation shell effects, such as hydrogen bonding,...

An `ab initio' molecular dynamics simulation of a Na+ ion in aqueous solution is presented and discussed. The calculation treats a Na+ ion and 32 water molecules with periodic boundary conditions on a cubic volume determined by an estimate of zero partial molar volume for this solute in water at normal density and at a temperature of 344+/-24 K. Analysis of the last half of the 12 ps trajectory shows 4.6 water molecules occupying the inner hydration shell of the Na+ ion on av...

In a floating water bridge the total radiation scattering of water stressed by a moderately strong electric field (1mV/nm) was compared to water without an applied electric field using X-ray and small angle neutron scattering. Structure refinement was carried out using the EPSR method and the TIP4P/2005 water model. These results did not reveal a significant difference in the local static structure of water however analysis of the simulation indicated that the modeled local p...

Theories of solvation free energies often involve electrostatic potentials at the position of a solute charge. Simulation calculations that apply cutoffs and periodic boundary conditions based on molecular centers result in center-dependent contributions to electrostatic energies due to a systematic sorting of charges in radial shells. This sorting of charges induces a surface-charge density at the cutoff sphere or simulation-box boundary that depends on the choice of molecul...

A new computational model for Sodium Chloride, the NaCl/{\epsilon}, is proposed. The Force Fields employed here for the description of the NaCl is based on a set of radial particle-particle pair potentials involving Lennard-Jones (LJ) and Coulombic forces. The parameterization is obtained fitting the density of the crystal and the density and the dielectric constant of the mixture of salt with water at diluted solution. Our model shows good agreement with the experimental val...

Determining the solvation free energies of single ions in water is one of the most fundamental problems in physical chemistry and yet many unresolved questions remain. In particular, the ability to decompose the solvation free energy into simple and intuitive contributions will have important implications for models of electrolyte solution. Here, we provide definitions of the various types of single ion solvation free energies based on different simulation protocols. We calcu...

It is a promising extension of the quantum mechanical/molecular mechanical (QM/MM) approach to incorporate the solvent molecules surrounding the QM solute into the QM region to ensure the adequate description of the electronic polarization of the solute. However, the solvent molecules in the QM region inevitably diffuse into the MM bulk during the QM/MM simulation. In this article we developed a simple and efficient method, referred to as boundary constraint with correction (...