Structural Models and Molecular Thermodynamics of Hydration of Ions and Small Molecules

Self Cite

Propylene carbonate (PC) wets graphite with a contact angle of 31• at ambient conditions. Molecular dynamics simulations agree with this contact angle after 40% reduction of the strength of graphite-C atom Lennard-Jones interactions with the solvent, relative to the models used initially. A simulated nano-scale PC droplet on graphite displays a pronounced layering tendency and an Aztex pyramid structure for the droplet. Extrapolation of the computed tensions of PC liquid-vapor interface estimates the critical temperature of PC accurately to about 3%. PC molecules lie flat on the PC liquid-vapor surface, and tend to project the propyl carbon toward the vapor phase. For close PC neighbors in liquid PC, an important packing motif stacks carbonate planes with the outer oxygen of one molecule snuggled into the positively charged propyl end of another molecule so that neighboring molecule dipole moments are approximately antiparallel. The calculated thermal expansion coefficient and the dielectric constants for liquid PC agree well with experiment. The distribution of PC molecule binding energies is closely Gaussian. Evaluation of the density of the coexisting vapor then permits estimation of the packing contribution to the PC chemical potential, and that contribution is about 2/3rds of the magnitude of the contributions due to attractive interactions, with opposite sign.

Section: E Balance Of Excluded Volume and Attractive Interactionsmentioning

confidence: 99%

Interfaces of propylene carbonate

Chaudhari

Pratt

Pesika

et al. 2013

Self Cite

“…[81][82][83][84][85][86] For the calculation of single-ion hydration free energies, two main strategies have been explored to date for achieving such a combination. The most common approach is the cluster-continuum approach, 53,56,[87][88][89][90][91] often referred to as a particular instance of quasi-chemical theory, [92][93][94] in which QM energy calculations are performed on small ion-water clusters (different sizes and geometries) and the bulk intrinsic hydration free energy is calculated by subsequently embedding the cluster configurations into a CE solvent environment. The second approach relies on Car-Parrinello 95 (CPMD) or Born-Oppenheimer 96 (BOMD) molecular dynamics simulations, performed at the density-functional theory (DFT) level.…”

Section: Introductionmentioning

confidence: 99%

Absolute proton hydration free energy, surface potential of water, and redox potential of the hydrogen electrode from first principles: QM/MM MD free-energy simulations of sodium and potassium hydration

Hofer

Hünenberger

2018

104

The absolute intrinsic hydration free energy G • H + ,wat of the proton, the surface electric potential jump χ • wat upon entering bulk water, and the absolute redox potential V • H + ,wat of the reference hydrogen electrode are cornerstone quantities for formulating single-ion thermodynamics on absolute scales. They can be easily calculated from each other but remain fundamentally elusive, i.e., they cannot be determined experimentally without invoking some extra-thermodynamic assumption (ETA). The Born model provides a natural framework to formulate such an assumption (Born ETA), as it automatically factors out the contribution of crossing the water surface from the hydration free energy. However, this model describes the short-range solvation inaccurately and relies on the choice of arbitrary ion-size parameters. In the present study, both shortcomings are alleviated by performing first-principle calculations of the hydration free energies of the sodium (Na + ) and potassium (K + ) ions. The calculations rely on thermodynamic integration based on quantum-mechanical molecular-mechanical (QM/MM) molecular dynamics (MD) simulations involving the ion and 2000 water molecules. The ion and its first hydration shell are described using a correlated ab initio method, namely resolution-of-identity second-order Møller-Plesset perturbation (RIMP2). The next hydration shells are described using the extended simple point charge water model (SPC/E). The hydration free energy is first calculated at the MM level and subsequently increased by a quantization term accounting for the transformation to a QM/MM description. It is also corrected for finite-size, approximate-electrostatics, and potentialsummation errors, as well as standard-state definition. These computationally intensive simulations provide accurate first-principle estimates for G • H + ,wat , χ • wat , and V • H + ,wat , reported with statistical errors based on a confidence interval of 99%. The values obtained from the independent Na + and K + simulations are in excellent agreement. In particular, the difference between the two hydration free energies, which is not an elusive quantity, is 73.9 ± 5.4 kJ mol 1 (K + minus Na + ), to be compared with the experimental value of 71.7 ± 2.8 kJ mol 1 . The calculated values of G • H + ,wat , χ • wat , and V • H + ,wat ( 1096.7 ± 6.1 kJ mol 1 , 0.10 ± 0.10 V, and 4.32 ± 0.06 V, respectively, averaging over the two ions) are also in remarkable agreement with the values recommended by Reif and Hünenberger based on a thorough analysis of the experimental literature ( 1100 ± 5 kJ mol 1 , 0.13 ± 0.10 V, and 4.28 ± 0.13 V, respectively). The QM/MM MD simulations are also shown to provide an accurate description of the hydration structure, dynamics, and energetics. Published by AIP Publishing.

“…[63][64][65][66] The n values also closely mirror the GC values. This is expected, given that the switchover from a water being in the hydration shell to outside occurs near the first minimum in g(r).…”

Section: Resultsmentioning

confidence: 52%

Instantaneous, parameter-free methods to define a solute’s hydration shell

Chatterjee

Higham

Henchman

2015

Modeling the hydration of mono-atomic anions from the gas phase to the bulk phase: The case of the halide Instantaneous, parameter-free methods to define a solute's hydration shell A range of methods are presented to calculate a solute's hydration shell from computer simulations of dilute solutions of monatomic ions and noble gas atoms. The methods are designed to be parameter-free and instantaneous so as to make them more general, accurate, and consequently applicable to disordered systems. One method is a modified nearest-neighbor method, another considers solute-water Lennard-Jones overlap followed by hydrogen-bond rearrangement, while three methods compare various combinations of water-solute and water-water forces. The methods are tested on a series of monatomic ions and solutes and compared with the values from cutoffs in the radial distribution function, the nearest-neighbor distribution functions, and the strongest-acceptor hydrogen bond definition for anions. The Lennard-Jones overlap method and one of the force-comparison methods are found to give a hydration shell for cations which is in reasonable agreement with that using a cutoff in the radial distribution function. Further modifications would be required, though, to make them capture the neighboring water molecules of noble-gas solutes if these weakly interacting molecules are considered to constitute the hydration shell. C 2015 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution 3.0 Unported License.[http://dx