Critical comparison of classical and quantum mechanical treatments of the phase equilibria of water

The Gibbs ensemble Monte Carlo simulation technique was used to compare the phase equilibria of the rigid TIP4P water model [Jorgensen et al., J. Chem. Phys. 79, 926 (1983)] utilizing classical and quantum statistical mechanics. The quantum statistical mechanical treatment generally resulted in lower liquid densities and higher vapor densities, narrowing the phase envelope. As a result, the calculated critical temperatures and normal boiling points were lower from the quantum simulations than the classical by 22 and 17 K, respectively, but the critical densities were equal within the estimated uncertainties. When the phase diagram from the quantum statistical mechanical treatment was increased by 22 K, it agreed with the classical results quite well throughout the entire simulated temperature range. A semiclassical treatment, involving a low order expansion in Planck's constant, resulted in good agreement with the path integral results for second virial coefficients, but gave densities and vapor pressures that fluctuated between the values for the classical and quantum statistical mechanics values, with no definite agreement with either.     

Global phase diagram for the honeycomb potential

We calculate the global phase diagram using classical statistical mechanics for an isotropic pair potential that has been previously [Rechtsman et al., Phys. Rev. Lett. 95, 228301 (2005)] shown to produce the low-coordinated two-dimensional honeycomb crystal as the ground-state structure. Low-coordinated crystals are of practical interest because they have desirable photonic band-gap properties. The phase diagram is obtained from Helmholtz free energies calculated using thermodynamic integration and Monte Carlo simulations. Our results show that the honeycomb crystal remains stable in the global phase diagram even after temperature effects are taken fully into account. Other stable phases in the phase diagram are high and low density triangular phases and a fluid phase. We find no evidence of gas-liquid or liquid-liquid phase coexistence.     

Self consistent tight binding model for dissociable water

We report results of development of a self consistent tight binding model for water. The model explicitly describes the electrons of the liquid self consistently, allows dissociation of the water and permits fast direct dynamics molecular dynamics calculations of the fluid properties. It is parameterized by fitting to first principles calculations on water monomers, dimers, and trimers. We report calculated radial distribution functions of the bulk liquid, a phase diagram and structure of solvated protons within the model as well as ac conductivity of a system of 96 water molecules of which one is dissociated. Structural properties and the phase diagram are in good agreement with experiment and first principles calculations. The estimated DC conductivity of a computational sample containing a dissociated water molecule was an order of magnitude larger than that reported from experiment though the calculated ratio of proton to hydroxyl contributions to the conductivity is very close to the experimental value. The conductivity results suggest a Grotthuss-like mechanism for the proton component of the conductivity.     

Free energy calculations for a flexible water model

In this work, we consider the problem of calculating the classical free energies of liquids and solids for molecular models with intramolecular flexibility. We show that thermodynamic integration from the fully-interacting solid of interest to a Debye crystal reference state, with anisotropic harmonic interactions derived from the Hessian of the original crystal, provides a straightforward route to calculating the Gibbs free energy of the solid. To calculate the molecular liquid free energy, it is essential to correctly account for contributions from both intermolecular and intramolecular motion; we employ thermodynamic integration to a Lennard-Jones reference fluid, coupled with direct evaluation of the molecular ro-vibrational partition function. These approaches are used to study the low-pressure classical phase diagram of the flexible q-TIP4P/F water model. We find that, while the experimental ice-I/liquid and ice-III/liquid coexistence lines are described reasonably well by this model, the ice-II phase is predicted to be metastable. In light of this finding, we go on to examine how the coupling between intramolecular flexibility and intermolecular interactions influences the computed phase diagram by comparing our results with those of the underlying rigid-body water model.     

An accurate and simple quantum model for liquid water

The path-integral molecular dynamics and centroid molecular dynamics methods have been applied to investigate the behavior of liquid water at ambient conditions starting from a recently developed simple point charge/flexible (SPC/Fw) model. Several quantum structural, thermodynamic, and dynamical properties have been computed and compared to the corresponding classical values, as well as to the available experimental data. The path-integral molecular dynamics simulations show that the inclusion of quantum effects results in a less structured liquid with a reduced amount of hydrogen bonding in comparison to its classical analog. The nuclear quantization also leads to a smaller dielectric constant and a larger diffusion coefficient relative to the corresponding classical values. Collective and single molecule time correlation functions show a faster decay than their classical counterparts. Good agreement with the experimental measurements in the low-frequency region is obtained for the quantum infrared spectrum, which also shows a higher intensity and a redshift relative to its classical analog. A modification of the original parametrization of the SPC/Fw model is suggested and tested in order to construct an accurate quantum model, called q-SPC/Fw, for liquid water. The quantum results for several thermodynamic and dynamical properties computed with the new model are shown to be in a significantly better agreement with the experimental data. Finally, a force-matching approach was applied to the q-SPC/Fw model to derive an effective quantum force field for liquid water in which the effects due to the nuclear quantization are explicitly distinguished from those due to the underlying molecular interactions. Thermodynamic and dynamical properties computed using standard classical simulations with this effective quantum potential are found in excellent agreement with those obtained from significantly more computationally demanding full centroid molecular dynamics simulations. The present results suggest that the inclusion of nuclear quantum effects into an empirical model for water enhances the ability of such model to faithfully represent experimental data, presumably through an increased ability of the model itself to capture realistic physical effects.     

A simple quantum statistical thermodynamics interpretation of an impressive phase diagram pressure shift upon (H/D) isotopic substitution in water + 3-methylpyridine

In a previous work (J. Phys. Chem. B 2003, 107, 9837), we reported liquid-liquid-phase splitting at negative pressures in mixtures of H2O + D2O + 3-methylpyridine (3-MP) at the limit of pure H2O as the solvent, thus extending for the first time the L-L phase diagrams to this metastable region. We showed that there is an intimate relation between pressure and solvent deuterium content. Isotopic substitution (H/D) in water provokes subtle entropic effects that, in turn, trigger a significant pressure shift, opening a pressure-wide miscibility window of as much as 1600 bar. Isotope effects are quantum in origin. Therefore, a model that is both pressure-dependent and considers quantization constitutes a necessary tool if one wishes to fully describe the p, T, x critical demixing in these systems. In the current work, the statistical-mechanical theory of isotope effects is combined with a compressible pressure-dependent model. This combination enabled us to predict successfully the overall L-L phase diagram via differences in the vibrational mode frequencies of water on its transfer from the pure state to that of dilution in 3-MP: each of the three librational modes undergo a calculated red-shift of -(250 +/- 30) cm(-1), while the overall internal frequencies contribution is estimated as a total +(400 +/- 25) cm(-1) blue-shift.     

Communication: Thermodynamics of water modeled using ab initio simulations

We regularize the potential distribution framework to calculate the excess free energy of liquid water simulated with the BLYP-D density functional. Assuming classical statistical mechanical simulations at 350 K model the liquid at 298 K, the calculated free energy is found in fair agreement with experiments, but the excess internal energy and hence also the excess entropy are not. The utility of thermodynamic characterization in understanding the role of high temperatures to mimic nuclear quantum effects and in evaluating ab initio simulations is noted.     

Analysis of nuclear quantum effects on hydrogen bonding

The impact of nuclear quantum effects on hydrogen bonding is investigated for a series of hydrogen fluoride (HF)n clusters and a partially solvated fluoride anion, F-(H2O). The nuclear quantum effects are included using the path integral formalism in conjunction with the Car-Parrinello molecular dynamics (PICPMD) method and using the second-order vibrational perturbation theory (VPT2) approach. For the HF clusters, a directional change in the impact of nuclear quantum effects on the hydrogen-bonding strength is observed as the clusters evolve toward the condensed phase. Specifically, the inclusion of nuclear quantum effects increases the F-F distances for the (HF)n=2-4 clusters and decreases the F-F distances for the (HF)n>4 clusters. This directional change occurs because the enhanced electrostatic interactions between the HF monomers become more dominant than the zero point energy effects of librational modes as the size of the HF clusters increases. For the F-(H2O) system, the inclusion of nuclear quantum effects decreases the F-O distance and strengthens the hydrogen bonding interaction between the fluoride anion and the water molecule because of enhanced electrostatic interactions. The vibrationally averaged 19F shielding constant for F-(H2O) is significantly lower than the value for the equilibrium geometry, indicating that the electronic density on the fluorine decreases as a result of the quantum delocalization of the shared hydrogen. Deuteration of this system leads to an increase in the vibrationally averaged F-O distance and nuclear magnetic shielding constant because of the smaller degree of quantum delocalization for deuterium.     

Efficient modeling of liquid phase photoemission spectra and reorganization energies: Difficult case of multiply charged anions

An efficient approach for quantitative modeling of liquid phase photoelectron spectra, reorganization energies, and redox potentials with DFT‐based molecular dynamics simulations is presented. The method is based on a large scale cluster‐continuum approach combined with the so‐called reflection principle (RP). Finite size clusters of solute molecules with solvating water molecules are at first generated using either classical molecular dynamics or molecular dynamics with a quantum thermostat which accounts for nuclear quantum effects. In the next step, the electron binding energies are calculated. Finite‐size corrections for (i) positions of electron binding energies and (ii) width of the spectrum are evaluated via a dielectric continuum approach. The performance of such a reflection principle with additional broadening approach (RP‐AB) for oxidation of multiply charged iron anions, [Fe(CN)₆]⁴⁻ and [Fe(CN)₆]³⁻ is demonstrated. The role of nuclear quantum effects is discussed as well as the relation between spectroscopic data and electrochemical quantities. Results are compared with recent liquid photoemission experiments, explaining the obstacles for applying liquid phase photoemission spectroscopy as a direct method for obtaining absolute redox potentials and suggesting a way to overcome them. © 2017 Wiley Periodicals, Inc.     

Contracted Schrödinger equation in quantum phase‐space

The phase space formulation of quantum mechanics is equivalent to standard quantum mechanics where averages are calculated by way of phase space integration as in the case of classical statistical mechanics. We derive the quantum hierarchy equations, often called the contracted Schrödinger equation, in the phase space representation of quantum mechanics which involves quasi‐distributions of position and momentum. We use the Wigner distribution for the phase space function and the Moyal phase space eigenvalue formulation to derive the hierarchy. We show that the hierarchy equations in the position, momentum, and position‐momentum representations are very similar in structure. © 2017 Wiley Periodicals, Inc.     

Measured and Predicted Solubility Phase Diagrams of Quaternary Systems LiBr-NaBr-MgBr2-H2O and LiBr-KBr-MgBr2-H2O at 298.15 K

The stable phase equilibria of quaternary systems LiBr-NaBr-MgBr₂-H₂O and LiBr-KBr-MgBr₂-H₂O at 298.15 K were studied by both experimental measurement(isothermal solution saturation method) and theoretical prediction(Pitzer model). The solubilities of the saturated solution have been determined experimentally and two stable phase diagrams and relevant water diagrams of the two quaternary systems were obtained. Results show that quaternary system LiBr-NaBr-MgBr₂-H₂O is hydrate II type as NaBr and NaBr·2H₂O coexistence. Its phase diagram consists of only one invariant point, four univariant curves, and five crystallization fields. The quaternary system LiBr-KBr-MgBr₂-H₂O is a complex type as the double salt KBr·MgBr₂·6H₂O formed. In addition to this double salt, the three single salts LiBr·2H₂O, KBr and MgBr₂·6H₂O also crystallize. In this paper, the solubilities of phase equilibria in above quaternary systems were also calculated by the Pitzer's electrolyte solution model. All the needed parameters can be obtained from the literature or be fitted by experimental data. On the Basis of the experimental and calculated results, the phase diagram of the quaternary system was plotted for comparison. It shows that the calculation results are consistent with the experimental ones.     

Impact of nuclear quantum effects on the molecular structure of bihalides and the hydrogen fluoride dimer

The structural impact of nuclear quantum effects is investigated for a set of bihalides, [XHX](-), X = F, Cl, and Br, and the hydrogen fluoride dimer. Structures are calculated with the vibrational self-consistent-field (VSCF) method, the second-order vibrational perturbation theory method (VPT2), and the nuclear-electronic orbital (NEO) approach. In the VSCF and VPT2 methods, the vibrationally averaged geometries are calculated for the Born-Oppenheimer electronic potential energy surface. In the NEO approach, the hydrogen nuclei are treated quantum mechanically on the same level as the electrons, and mixed nuclear-electronic wave functions are calculated variationally with molecular orbital methods. Electron-electron and electron-proton dynamical correlation effects are included in the NEO approach using second-order perturbation theory (NEO-MP2). The nuclear quantum effects are found to alter the distances between the heavy atoms by 0.02-0.05 A for the systems studied. These effects are of similar magnitude as the electron correlation effects. For the bihalides, inclusion of the nuclear quantum effects with the NEO-MP2 or the VSCF method increases the X-X distance. The bihalide X-X distances are similar for both methods and are consistent with two-dimensional grid calculations and experimental values, thereby validating the use of the computationally efficient NEO-MP2 method for these types of systems. For the hydrogen fluoride dimer, inclusion of nuclear quantum effects decreases the F-F distance with the NEO-MP2 method and increases the F-F distance with the VSCF and VPT2 methods. The VPT2 F-F distances for the hydrogen fluoride dimer and the deuterated form are consistent with the experimentally determined values. The NEO-MP2 F-F distance is in excellent agreement with the distance obtained experimentally for a model that removes the large amplitude bending motions. The analysis of these calculations provides insight into the significance of electron-electron and electron-proton correlation, anharmonicity of the vibrational modes, and nonadiabatic effects for hydrogen-bonded systems.     

Amphiphilic character of the hydrated proton in methanol-water solutions

The hydrated proton was studied in methanol-water solutions of varying methanol concentrations using the multistate empirical valence bond simulation method. Amphiphile-like behavior of the hydrated proton was noted from its anisotropic association with the methanol methyl groups. Molecular length immiscibility was also characterized through the enumeration of water and protonated water clusters. Excess proton diffusion was calculated across the varying methanol concentrations and found to be in good agreement with experiment after correcting for nuclear quantum effects.     

Efficient method to include nuclear quantum effects in the determination of phase boundaries

We developed a methodology to assess nuclear quantum effects in phase boundaries calculations that is based on the dynamical integration of Clausius-Clapeyron equation using path integral simulations. The technique employs non-equilibrium simulations that are very efficient. The approach was applied to the calculation of the melting line of Ne in an interval of pressures ranging from 1 to 3366 bar. Our results show a very good agreement with both experimental findings and results from previous calculations. The methodology can be applied to solid and liquid phases, without limitations regarding anharmonicities. The method allows the computation of coexistence lines for wide intervals of pressure and temperature using, in principle, a single simulation.     

Parameterization and Application of the General Amber Force Field to Model Fluoro Substituted Furanose Moieties and Nucleosides

Molecular mechanics force field calculations have historically shown significant limitations in modeling the energetic and conformational interconversions of highly substituted furanose rings. This is primarily due to the gauche effect that is not easily captured using pairwise energy potentials. In this study, we present a refinement to the set of torsional parameters in the General Amber Force Field (gaff) used to calculate the potential energy of mono, di-, and gem-fluorinated nucleosides. The parameters were optimized to reproduce the pseudorotation phase angle and relative energies of a diverse set of mono- and difluoro substituted furanose ring systems using quantum mechanics umbrella sampling techniques available in the IpolQ engine in the Amber suite of programs. The parameters were developed to be internally consistent with the gaff force field and the TIP3P water model. The new set of angle and dihedral parameters and partial charges were validated by comparing the calculated phase angle probability to those obtained from experimental nuclear magnetic resonance experiments.     

Thermodynamic properties of alloys of the binary In–La system

The thermochemical properties of melts of the binary In–La system were studied by the calorimetry method at 1250–1480 K over the whole concentration interval. It was shown that significant negative heat effects of mixing are characteristic features for these melts. Using the ideal associated solution (IAS) model, the activities of components, Gibbs energies and the entropies of mixing in the alloys, and the phase diagram of this system were calculated. They agree with the data from literature.     

Thermodynamic properties of alloys of the binary In–Yb system

The thermochemical properties of melts of the binary In–Yb system were studied by the calorimetry method at 1160–1380 K over the whole concentration interval. It was shown that significant negative heat effects of mixing are characteristic features for these melts. Using the ideal associated solution (IAS) model, the activities of components, Gibbs energies and the entropies of mixing in the alloys, and the phase diagram of this system were calculated. They agree with the data from literature.     

Communication: phase transitions, criticality, and three-phase coexistence in constrained cell models

In simulation studies of fluid-solid transitions, the solid phase is usually modeled as a constrained system in which each particle is confined to move in a single Wigner-Seitz cell. The constrained cell model has been used in the determination of fluid-solid coexistence via thermodynamic integration and other techniques. In the present work, the phase diagram of such a constrained system of Lennard-Jones particles is determined from constant-pressure simulations. The pressure-density isotherms exhibit inflection points which are interpreted as the mechanical stability limit of the solid phase. The phase diagram of the constrained system contains a critical and a triple point. The temperature and pressure at the critical and the triple point are both higher than those of the unconstrained system due to the reduction in the entropy caused by the single occupancy constraint.