Solvation properties of Li+and CI−in water: molecular dynamics simulation with a non-rigid model

Surface Energies and Structure of Salt-Brine Interfaces

Permeability of salt formations is controlled by the equilibrium between the salt-brine and salt-salt interfaces described by the dihedral angle, which can change with the composition of the intergranular brine. Here, classical molecular dynamics (MD) simulations were used to investigate the structure and properties of the salt-brine interface to provide insight into the stability of salt systems. Mixed NaCl-KCl brines were investigated to explore differences in ion size on the surface energy and interface structure. Nonlinearity was noted in the salt-brine surface energy with increasing KCl concentration, and the addition of 10% KCl increased surface energies by 2-3 times (5.0 M systems). Size differences in Na⁺ and K⁺ ions altered the packing of dissolved ions and water molecules at the interface, impacting the surface energy. Additionally, ions at the interface had lower numbers of coordinating water molecules than those in the bulk and increased hydration for ions in systems with 100% NaCl or 100% KCl brines. Ultimately, small changes in brine composition away from pure NaCl altered the structure of the salt-brine interface, impacting the dihedral angle and the predicted equilibrium permeability of salt formations.

Effect of interionic interactions on the structure and dynamics of ionic solvation shells in aqueous electrolyte solutions

Molecular dynamics (MD) simulations to explore the structure and dynamics of the ionic solvation shell of alkali ions and halide ions in aqueous solution.

Absolute hydration free energy scale for alkali and halide ions established from simulations with a polarizable force field

A polarizable potential function for the hydration of alkali and halide ions is developed on the basis of the recent SWM4-DP water model [Lamoureux, G.; MacKerell, A. D., Jr.; Roux, B. J. Chem. Phys. 2003, 119, 5185]. Induced polarization is incorporated using classical Drude oscillators that are treated as auxiliary dynamical degrees of freedom. The ions are represented as polarizable Lennard-Jones centers, whose parameters are optimized to reproduce the binding energies of gas-phase monohydrates and the hydration free energies in the bulk liquid. Systematic exploration of the parameters shows that the monohydrate binding energies can be consistent with a unique hydration free energy scale if the computed hydration free energies incorporate the contribution from the air/water interfacial electrostatic potential (-540 mV for SWM4-DP). The final model, which can satisfyingly reproduce both gas and bulk-phase properties, corresponds to an absolute scale in which the intrinsic hydration free energy of the proton is -247 kcal/mol.

Absolute free energy calculations by thermodynamic integration in four spatial dimensions

An optimized technique for calculating the excess chemical potential of small molecules in dense liquids and the binding affinity of molecular ligands to biomolecules is reported. In this method, a molecular species is coupled to the system of interest via a nonphysical fourth spatial dimension w through which insertion or extraction can be carried out [R. Pomes, E. Eisenmesser, C. B. Post et al., J. Chem. Phys. 111, 3387 (1999)]. Molecular simulations are used to compute the potential of mean force (PMF) acting on the solute molecule in the fourth dimension. The excess chemical potential of that molecule is obtained as the difference in the PMF between fully coupled and fully decoupled systems. The simplicity, efficiency, and generality of the method are demonstrated for the calculation of the hydration free energies of water and methanol as well as sodium, cesium, and chloride ions. A significant advantage over other methods is that the 4D-PMF approach provides a single effective and general route for decoupling all nonbonded interactions (i.e., both Lennard-Jones and Coulombic) at once for both neutral and charged solutes. Direct calculation of the mean force from thermodynamic integration is shown to be more computationally efficient than calculating the PMF from umbrella sampling. Statistical error analysis suggests a simple strategy for optimizing sampling. The detailed analysis of systematic errors arising from the truncation of Coulombic interactions in a solvent droplet of finite size leads to straightforward corrections to ionic hydration free energies.