+ Site Statistics
+ Search Articles
+ PDF Full Text Service
How our service works
Request PDF Full Text
+ Follow Us
Follow on Facebook
Follow on Twitter
Follow on LinkedIn
+ Subscribe to Site Feeds
Most Shared
PDF Full Text
+ Translate
+ Recently Requested

Computation of methodology-independent single-ion solvation properties from molecular simulations. IV. Optimized Lennard-Jones interaction parameter sets for the alkali and halide ions in water



Computation of methodology-independent single-ion solvation properties from molecular simulations. IV. Optimized Lennard-Jones interaction parameter sets for the alkali and halide ions in water



Journal of Chemical Physics 134(14): 144104



The raw single-ion solvation free energies computed from atomistic (explicit-solvent) simulations are extremely sensitive to the boundary conditions and treatment of electrostatic interactions used during these simulations. However, as shown recently [M. A. Kastenholz and P. H. Hünenberger, J. Chem. Phys. 124, 224501 (2006); M. M. Reif and P. H. Hünenberger, J. Chem. Phys. 134, 144103 (2010)], the application of appropriate correction terms permits to obtain methodology-independent results. The corrected values are then exclusively characteristic of the underlying molecular model including in particular the ion-solvent van der Waals interaction parameters, determining the effective ion size and the magnitude of its dispersion interactions. In the present study, the comparison of calculated (corrected) hydration free energies with experimental data (along with the consideration of ionic polarizabilities) is used to calibrate new sets of ion-solvent van der Waals (Lennard-Jones) interaction parameters for the alkali (Li(+), Na(+), K(+), Rb(+), Cs(+)) and halide (F(-), Cl(-), Br(-), I(-)) ions along with either the SPC or the SPC/E water models. The experimental dataset is defined by conventional single-ion hydration free energies [Tissandier et al., J. Phys. Chem. A 102, 7787 (1998); Fawcett, J. Phys. Chem. B 103, 11181] along with three plausible choices for the (experimentally elusive) value of the absolute (intrinsic) hydration free energy of the proton, namely, ΔG(hyd)(⊖)[H(+)] = -1100, -1075 or -1050 kJ mol(-1), resulting in three sets L, M, and H for the SPC water model and three sets L(E), M(E), and H(E) for the SPC/E water model (alternative sets can easily be interpolated to intermediate ΔG(hyd)(⊖)[H(+)] values). The residual sensitivity of the calculated (corrected) hydration free energies on the volume-pressure boundary conditions and on the effective ionic radius entering into the calculation of the correction terms is also evaluated and found to be very limited. Ultimately, it is expected that comparison with other experimental ionic properties (e.g., derivative single-ion solvation properties, as well as data concerning ionic crystals, melts, solutions at finite concentrations, or nonaqueous solutions) will permit to validate one specific set and thus, the associated ΔG(hyd)(⊖)[H(+)] value (atomistic consistency assumption). Preliminary results (first-peak positions in the ion-water radial distribution functions, partial molar volumes of ionic salts in water, and structural properties of ionic crystals) support a value of ΔG(hyd)(⊖)[H(+)] close to -1100 kJ·mol(-1).

Please choose payment method:






(PDF emailed within 0-6 h: $19.90)

Accession: 052272795

Download citation: RISBibTeXText

PMID: 21495739

DOI: 10.1063/1.3567022


Related references

Computation of methodology-independent single-ion solvation properties from molecular simulations. III. Correction terms for the solvation free energies, enthalpies, entropies, heat capacities, volumes, compressibilities, and expansivities of solvated ions. Journal of Chemical Physics 134(14): 144103, 2011

Molecular dynamics simulations of the dynamic and energetic properties of alkali and halide ions using water-model-specific ion parameters. Journal of Physical Chemistry. B 113(40): 13279-13290, 2009

Computation of methodology-independent ionic solvation free energies from molecular simulations. I. The electrostatic potential in molecular liquids. Journal of Chemical Physics 124(12): 124106, 2006

Computation of methodology-independent ionic solvation free energies from molecular simulations. II. The hydration free energy of the sodium cation. Journal of Chemical Physics 124(22): 224501, 2006

Optimized ensemble Monte Carlo simulations of dense Lennard-Jones fluids. Journal of Chemical Physics 123(20): 204501, 2005

Molecular dynamics simulations of the interaction between polyhydroxylated compounds and Lennard-Jones walls: preferential affinity/exclusion effects and their relevance for bioprotection. Molecular Simulation 36(9): 708-728, 2010

Thermodynamics of small alkali metal halide cluster ions: comparison of classical molecular simulations with experiment and quantum chemistry. Journal of Physical Chemistry. a 119(3): 488-500, 2015

Lennard-Jones parameters for small diameter carbon nanotubes and water for molecular mechanics simulations from van der Waals density functional calculations. Journal of Computational Chemistry 33(6): 652-658, 2012

Computation of low pressures by molecular dynamics based on Lennard-Jones (12-6) potentials. Molecular Physics 55(6): 1339-1350, 1985

Solvation thermodynamics of alkali and halide ions in ionic liquids through integral equations. Journal of Chemical Physics 129(7): 074509, 2008

Simulated Solvation of Organic Ions II: Study of Linear Alkylated Carboxylate Ions in Water Nanodrops and in Liquid Water. Propensity for Air/Water Interface and Convergence to Bulk Solvation Properties. Journal of Physical Chemistry. B 119(36): 12094-12107, 2015

Crystal lattice properties fully determine short-range interaction parameters for alkali and halide ions. Journal of Chemical Physics 137(6): 064104, 2012

Thermodynamic, Structural and Transport Properties of Lennard-Jones Liquid Systems. A Molecular Dynamics Simulations of Liquid Helium, Neon, Methane and Nitrogen. International Journal of Molecular Sciences 4(12): 595-606, 2003

Langevin and molecular dynamics simulations of Lennard-Jones liquids. Molecular Physics 55(5): 1063-1074, 1985

Interfacing CRYSTAL/AMBER to Optimize QM/MM Lennard⁻Jones Parameters for Water and to Study Solvation of TiO₂ Nanoparticles. Molecules 23(11):, 2018