Hydrogen bonds and van der waals forces in ice at ambient and high pressures

Hydrogen-bond potential for ice VIII-X phase transition

Repulsive force between the O-H bonding electrons and the O:H nonbonding pair within hydrogen bond (O-H:O) is an often overlooked interaction which dictates the extraordinary recoverability and sensitivity of water and ice. Here, we present a potential model for this hidden force opposing ice compression of ice VIII-X phase transition based on the density functional theory (DFT) and neutron scattering observations. We consider the H-O bond covalent force, the O:H nonbond dispersion force, and the hidden force to approach equilibrium under compression. Due to the charge polarization within the O:H-O bond, the curvatures of the H-O bond and the O:H nonbond potentials show opposite sign before transition, resulting in the asymmetric relaxation of H-O and O:H (O:H contraction and H-O elongation) and the H+ proton centralization towards phase X. When cross the VIII-X phase boundary, both H-O and O:H contract slightly. The potential model reproduces the VIII-X phase transition as observed in experiment. Development of the potential model may provide a choice for further calculations of water anomalies.

Raman and IR Spectra of Ice Ih and Ice XI with an Assessment of DFT Methods

IR and Raman spectroscopic technology can be directly used to identify the occurrence of ferroelectric ice XI in laboratory or extraterrestrial settings. The performance of 16 different DFT methods applied on the ice Ih, VIII, IX, and XI crystal phases is evaluated. Based on a selected DFT computational scheme, the IR and Raman spectra of ice Ih and XI are derived and compared. When the spectra, both IR and Raman, of ice Ih and ice XI are compared, the librational vibrations are found to be the most affected by the proton ordering. The spectroscopic fingerprint of ice XI can be used to distinguish ferroelectric ice XI from ice Ih in the universe. Furthermore, the existence of only one kind of H-bond in ice Ih is demonstrated from the overlapping subspectra for different types of H-bonded pair configurations in 16 isomers of ice Ih, which provides an illustration to the historic debate on whether one or two kinds of H-bonds existed in ice.

The fate of carbon dioxide in water-rich fluids under extreme conditions

Investigating the fate of dissolved carbon dioxide under extreme conditions is critical to understanding the deep carbon cycle in Earth, a process that ultimately influences global climate change. We used first-principles molecular dynamics simulations to study carbonates and carbon dioxide dissolved in water at pressures (P) and temperatures (T) approximating the conditions of Earth's upper mantle. Contrary to popular geochemical models assuming that molecular CO₂(aq) is the major carbon species present in water under deep Earth conditions, we found that at 11 GPa and 1000 K, carbon exists almost entirely in the forms of solvated carbonate ([Formula: see text]) and bicarbonate ([Formula: see text]) ions and that even carbonic acid [H₂CO₃(aq)] is more abundant than CO₂(aq). Furthermore, our simulations revealed that ion pairing between Na⁺ and [Formula: see text]/[Formula: see text] is greatly affected by P-T conditions, decreasing with increasing pressure at 800 to 1000 K. Our results suggest that in Earth's upper mantle, water-rich geofluids transport a majority of carbon in the form of rapidly interconverting [Formula: see text] and [Formula: see text] ions, not solvated CO₂(aq) molecules.

Getting the Right Answers for the Right Reasons: Toward Predictive Molecular Simulations of Water with Many-Body Potential Energy Functions

The central role played by water in fundamental processes relevant to different disciplines, including chemistry, physics, biology, materials science, geology, and climate research, cannot be overemphasized. It is thus not surprising that, since the pioneering work by Stillinger and Rahman, many theoretical and computational studies have attempted to develop a microscopic description of the unique properties of water under different thermodynamic conditions. Consequently, numerous molecular models based on either molecular mechanics or ab initio approaches have been proposed over the years. However, despite continued progress, the correct prediction of the properties of water from small gas-phase clusters to the liquid phase and ice through a single molecular model remains challenging. To large extent, this is due to the difficulties encountered in the accurate modeling of the underlying hydrogen-bond network in which both number and strength of the hydrogen bonds vary continuously as a result of a subtle interplay between energetic, entropic, and nuclear quantum effects. In the past decade, the development of efficient algorithms for correlated electronic structure calculations of small molecular complexes, accompanied by tremendous progress in the analytical representation of multidimensional potential energy surfaces, opened the doors to the design of highly accurate potential energy functions built upon rigorous representations of the many-body expansion (MBE) of the interaction energies. This Account provides a critical overview of the performance of the MB-pol many-body potential energy function through a systematic analysis of energetic, structural, thermodynamic, and dynamical properties as well as of vibrational spectra of water from the gas to the condensed phase. It is shown that MB-pol achieves unprecedented accuracy across all phases of water through a quantitative description of each individual term of the MBE, with a physically correct representation of both short- and long-range many-body contributions. Comparisons with experimental data probing different regions of the water potential energy surface from clusters to bulk demonstrate that MB-pol represents a major step toward the long-sought-after "universal model" capable of accurately describing the molecular properties of water under different conditions and in different environments. Along this path, future challenges include the extension of the many-body scheme adopted by MB-pol to the description of generic solutes as well as the integration of MB-pol in an efficient theoretical and computational framework to model acid-base reactions in aqueous environments. In this context, given the nontraditional form of the MB-pol energy and force expressions, synergistic efforts by theoretical/computational chemists/physicists and computer scientists will be critical for the development of high-performance software for many-body molecular dynamics simulations.

How van der Waals interactions determine the unique properties of water

Whereas the interactions between water molecules are dominated by strongly directional hydrogen bonds (HBs), it was recently proposed that relatively weak, isotropic van der Waals (vdW) forces are essential for understanding the properties of liquid water and ice. This insight was derived from ab initio computer simulations, which provide an unbiased description of water at the atomic level and yield information on the underlying molecular forces. However, the high computational cost of such simulations prevents the systematic investigation of the influence of vdW forces on the thermodynamic anomalies of water. Here, we develop efficient ab initio-quality neural network potentials and use them to demonstrate that vdW interactions are crucial for the formation of water's density maximum and its negative volume of melting. Both phenomena can be explained by the flexibility of the HB network, which is the result of a delicate balance of weak vdW forces, causing, e.g., a pronounced expansion of the second solvation shell upon cooling that induces the density maximum.

Modeling Polymorphic Molecular Crystals with Electronic Structure Theory

Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.

Benchmarking DFT and semiempirical methods on structures and lattice energies for ten ice polymorphs

Water in different phases under various external conditions is very important in bio-chemical systems and for material science at surfaces. Density functional theory methods and approximations thereof have to be tested system specifically to benchmark their accuracy regarding computed structures and interaction energies. In this study, we present and test a set of ten ice polymorphs in comparison to experimental data with mass densities ranging from 0.9 to 1.5 g/cm(3) and including explicit corrections for zero-point vibrational and thermal effects. London dispersion inclusive density functionals at the generalized gradient approximation (GGA), meta-GGA, and hybrid level as well as alternative low-cost molecular orbital methods are considered. The widely used functional of Perdew, Burke and Ernzerhof (PBE) systematically overbinds and overall provides inconsistent results. All other tested methods yield reasonable to very good accuracy. BLYP-D3(atm) gives excellent results with mean absolute errors for the lattice energy below 1 kcal/mol (7% relative deviation). The corresponding optimized structures are very accurate with mean absolute relative deviations (MARDs) from the reference unit cell volume below 1%. The impact of Axilrod-Teller-Muto (atm) type three-body dispersion and of non-local Fock exchange is small but on average their inclusion improves the results. While the density functional tight-binding model DFTB3-D3 performs well for low density phases, it does not yield good high density structures. As low-cost alternative for structure related problems, we recommend the recently introduced minimal basis Hartree-Fock method HF-3c with a MARD of about 3%.

Two Dimensional Ice from First Principles: Structures and Phase Transitions

Despite relevance to disparate areas such as cloud microphysics and tribology, major gaps in the understanding of the structures and phase transitions of low-dimensional water ice remain. Here, we report a first principles study of confined 2D ice as a function of pressure. We find that at ambient pressure hexagonal and pentagonal monolayer structures are the two lowest enthalpy phases identified. Upon mild compression, the pentagonal structure becomes the most stable and persists up to ∼2  GPa, at which point the square and rhombic phases are stable. The square phase agrees with recent experimental observations of square ice confined within graphene sheets. This work provides a fresh perspective on 2D confined ice, highlighting the sensitivity of the structures observed to both the confining pressure and the width.

Energy benchmarks for methane-water systems from quantum Monte Carlo and second-order Møller-Plesset calculations

The quantum Monte Carlo (QMC) technique is used to generate accurate energy benchmarks for methane-water clusters containing a single methane monomer and up to 20 water monomers. The benchmarks for each type of cluster are computed for a set of geometries drawn from molecular dynamics simulations. The accuracy of QMC is expected to be comparable with that of coupled-cluster calculations, and this is confirmed by comparisons for the CH4-H2O dimer. The benchmarks are used to assess the accuracy of the second-order Møller-Plesset (MP2) approximation close to the complete basis-set limit. A recently developed embedded many-body technique is shown to give an efficient procedure for computing basis-set converged MP2 energies for the large clusters. It is found that MP2 values for the methane binding energies and the cohesive energies of the water clusters without methane are in close agreement with the QMC benchmarks, but the agreement is aided by partial cancelation between 2-body and beyond-2-body errors of MP2. The embedding approach allows MP2 to be applied without loss of accuracy to the methane hydrate crystal, and it is shown that the resulting methane binding energy and the cohesive energy of the water lattice agree almost exactly with recently reported QMC values.

Singles correlation energy contributions in solids

The random phase approximation to the correlation energy often yields highly accurate results for condensed matter systems. However, ways how to improve its accuracy are being sought and here we explore the relevance of singles contributions for prototypical solid state systems. We set out with a derivation of the random phase approximation using the adiabatic connection and fluctuation dissipation theorem, but contrary to the most commonly used derivation, the density is allowed to vary along the coupling constant integral. This yields results closely paralleling standard perturbation theory. We re-derive the standard singles of Görling-Levy perturbation theory [A. Görling and M. Levy, Phys. Rev. A 50, 196 (1994)], highlight the analogy of our expression to the renormalized singles introduced by Ren and coworkers [Phys. Rev. Lett. 106, 153003 (2011)], and introduce a new approximation for the singles using the density matrix in the random phase approximation. We discuss the physical relevance and importance of singles alongside illustrative examples of simple weakly bonded systems, including rare gas solids (Ne, Ar, Xe), ice, adsorption of water on NaCl, and solid benzene. The effect of singles on covalently and metallically bonded systems is also discussed.

Water on BN doped benzene: a hard test for exchange-correlation functionals and the impact of exact exchange on weak binding

Density functional theory (DFT) studies of weakly interacting complexes have recently focused on the importance of van der Waals dispersion forces, whereas the role of exchange has received far less attention. Here, by exploiting the subtle binding between water and a boron and nitrogen doped benzene derivative (1,2-azaborine) we show how exact exchange can alter the binding conformation within a complex. Benchmark values have been calculated for three orientations of the water monomer on 1,2-azaborine from explicitly correlated quantum chemical methods, and we have also used diffusion quantum Monte Carlo. For a host of popular DFT exchange-correlation functionals we show that the lack of exact exchange leads to the wrong lowest energy orientation of water on 1,2-azaborine. As such, we suggest that a high proportion of exact exchange and the associated improvement in the electronic structure could be needed for the accurate prediction of physisorption sites on doped surfaces and in complex organic molecules. Meanwhile to predict correct absolute interaction energies an accurate description of exchange needs to be augmented by dispersion inclusive functionals, and certain non-local van der Waals functionals (optB88- and optB86b-vdW) perform very well for absolute interaction energies. Through a comparison with water on benzene and borazine (B3N3H6) we show that these results could have implications for the interaction of water with doped graphene surfaces, and suggest a possible way of tuning the interaction energy.

Communication: On the stability of ice 0, ice i, and I(h)

Using ab initio methods, we examine the stability of ice 0, a recently proposed tetragonal form of ice implicated in the homogeneous freezing of water [J. Russo, F. Romano, and H. Tanaka, Nat. Mater. 13, 670 (2014)]. Vibrational frequencies are computed across the complete Brillouin Zone using Density Functional Theory (DFT), to confirm mechanical stability and quantify the free energy of ice 0 relative to ice I(h). The robustness of this result is tested via dispersion corrected semi-local and hybrid DFT, and Quantum Monte-Carlo calculation of lattice energies. Results indicate that popular molecular models only slightly overestimate the stability of ice zero. In addition, we study all possible realisations of proton disorder within the ice zero unit cell, and identify the ground state as ferroelectric. Comparisons are made to other low density metastable forms of ice, suggesting that the ice i structure [C. J. Fennel and J. D. Gezelter, J. Chem. Theory Comput. 1, 662 (2005)] may be equally relevant to ice formation.

Comparison of x-ray absorption spectra between water and ice: new ice data with low pre-edge absorption cross-section

The effect of crystal growth conditions on the O K-edge x-ray absorption spectra of ice is investigated through detailed analysis of the spectral features. The amount of ice defects is found to be minimized on hydrophobic surfaces, such as BaF2(111), with low concentration of nucleation centers. This is manifested through a reduction of the absorption cross-section at 535 eV, which is associated with distorted hydrogen bonds. Furthermore, a connection is made between the observed increase in spectral intensity between 544 and 548 eV and high-symmetry points in the electronic band structure, suggesting a more extended hydrogen-bond network as compared to ices prepared differently. The spectral differences for various ice preparations are compared to the temperature dependence of spectra of liquid water upon supercooling. A double-peak feature in the absorption cross-section between 540 and 543 eV is identified as a characteristic of the crystalline phase. The connection to the interpretation of the liquid phase O K-edge x-ray absorption spectrum is extensively discussed.

On the room-temperature phase diagram of high pressure hydrogen: an ab initio molecular dynamics perspective and a diffusion Monte Carlo study

The finite-temperature phase diagram of hydrogen in the region of phase IV and its neighborhood was studied using the ab initio molecular dynamics (MD) and the ab initio path-integral molecular dynamics (PIMD). The electronic structures were analyzed using the density-functional theory (DFT), the random-phase approximation, and the diffusion Monte Carlo (DMC) methods. Taking the state-of-the-art DMC results as benchmark, comparisons of the energy differences between structures generated from the MD and PIMD simulations, with molecular and dissociated hydrogens, respectively, in the weak molecular layers of phase IV, indicate that standard functionals in DFT tend to underestimate the dissociation barrier of the weak molecular layers in this mixed phase. Because of this underestimation, inclusion of the quantum nuclear effects (QNEs) in PIMD using electronic structures generated with these functionals leads to artificially dissociated hydrogen layers in phase IV and an error compensation between the neglect of QNEs and the deficiencies of these functionals in standard ab initio MD simulations exists. This analysis partly rationalizes why earlier ab initio MD simulations complement so well the experimental observations. The temperature and pressure dependencies for the stability of phase IV were also studied in the end and compared with earlier results.

The individual and collective effects of exact exchange and dispersion interactions on the ab initio structure of liquid water

In this work, we report the results of a series of density functional theory (DFT) based ab initio molecular dynamics (AIMD) simulations of ambient liquid water using a hierarchy of exchange-correlation (XC) functionals to investigate the individual and collective effects of exact exchange (Exx), via the PBE0 hybrid functional, non-local van der Waals/dispersion (vdW) interactions, via a fully self-consistent density-dependent dispersion correction, and an approximate treatment of nuclear quantum effects, via a 30 K increase in the simulation temperature, on the microscopic structure of liquid water. Based on these AIMD simulations, we found that the collective inclusion of Exx and vdW as resulting from a large-scale AIMD simulation of (H2O)128 significantly softens the structure of ambient liquid water and yields an oxygen-oxygen structure factor, SOO(Q), and corresponding oxygen-oxygen radial distribution function, gOO(r), that are now in quantitative agreement with the best available experimental data. This level of agreement between simulation and experiment demonstrated herein originates from an increase in the relative population of water molecules in the interstitial region between the first and second coordination shells, a collective reorganization in the liquid phase which is facilitated by a weakening of the hydrogen bond strength by the use of a hybrid XC functional, coupled with a relative stabilization of the resultant disordered liquid water configurations by the inclusion of non-local vdW/dispersion interactions. This increasingly more accurate description of the underlying hydrogen bond network in liquid water also yields higher-order correlation functions, such as the oxygen-oxygen-oxygen triplet angular distribution, POOO(θ), and therefore the degree of local tetrahedrality, as well as electrostatic properties, such as the effective molecular dipole moment, that are in much better agreement with experiment.

Benchmarking the performance of density functional theory and point charge force fields in their description of sI methane hydrate against diffusion Monte Carlo

High quality reference data from diffusion Monte Carlo calculations are presented for bulk sI methane hydrate, a complex crystal exhibiting both hydrogen-bond and dispersion dominated interactions. The performance of some commonly used exchange-correlation functionals and all-atom point charge force fields is evaluated. Our results show that none of the exchange-correlation functionals tested are sufficient to describe both the energetics and the structure of methane hydrate accurately, while the point charge force fields perform badly in their description of the cohesive energy but fair well for the dissociation energetics. By comparing to ice Ih, we show that a good prediction of the volume and cohesive energies for the hydrate relies primarily on an accurate description of the hydrogen bonded water framework, but that to correctly predict stability of the hydrate with respect to dissociation to ice Ih and methane gas, accuracy in the water-methane interaction is also required. Our results highlight the difficulty that density functional theory faces in describing both the hydrogen bonded water framework and the dispersion bound methane.

Molecular to atomic phase transition in hydrogen under high pressure

The metallization of high-pressure hydrogen, together with the associated molecular to atomic transition, is one of the most important problems in the field of high-pressure physics. It is also currently a matter of intense debate due to the existence of conflicting experimental reports on the observation of metallic hydrogen on a diamond-anvil cell. Theoretical calculations have typically relied on a mean-field description of electronic correlation through density functional theory, a theory with well-known limitations in the description of metal-insulator transitions. In fact, the predictions of the pressure-driven dissociation of molecules in high-pressure hydrogen by density functional theory is strongly affected by the chosen exchange-correlation functional. In this Letter, we use quantum Monte Carlo calculations to study the molecular to atomic transition in hydrogen. We obtain a transition pressure of 447(3) GPa, in excellent agreement with the best experimental estimate of the transition 450 GPa based on an extrapolation to zero band gap from experimental measurements. Additionally, we find that C2/c is stable almost up to the molecular to atomic transition, in contrast to previous density functional theory (DFT) and DFT+quantum Monte Carlo studies which predict large stability regimes for intermediary molecular phases.

Predicting crystal structures and properties of matter under extreme conditions via quantum mechanics: the pressure is on

The role of quantum mechanical calculations in understanding and predicting the behavior of matter at extreme pressures is discussed in this feature contribution.

Role of hydrogen-bonding and its interplay with octahedral tilting in CH3NH3PbI3

Computed Kohn–Sham energies as a function of torsion angle for three rotational modes of the methylammonium group in orthorhombic CH₃NH₃PbI₃.

The accurate calculation of the band gap of liquid water by means of GW corrections applied to plane-wave density functional theory molecular dynamics simulations

Knowledge about the intrinsic electronic properties of water is imperative for understanding the behaviour of aqueous solutions that are used throughout biology, chemistry, physics, and industry.

Exploring the conformational preferences of 20-residue peptides in isolation: Ac-Ala19-Lys + H+vs. Ac-Lys-Ala19 + H+ and the current reach of DFT

Using a high-level density functional and an exhaustive search of conformation space, the predicted conformation of a 20-amino acid peptide explains two seemingly contradictory experiments.

Periodic MP2, RPA, and Boundary Condition Assessment of Hydrogen Ordering in Ice XV

Ice XV is the hydrogen-ordered form of the ice VI phase whose structure was predicted to be Cc and ferroelectric using periodic DFT approaches. However, neutron diffraction and Raman spectroscopy data show the structure to have P1̅ symmetry and to be antiferroelectric. Recent work1 using fragment-based MP2 and CCSD(T) approaches predicts the experimental structure as the ground state. We have reconsidered this problem using fully periodic MP2 and RPA approaches and find that the ferroelectric Cc structure is the lowest energy configuration. However, ubiquitously employed tinfoil boundary conditions stabilize polar structures. We suggest that ferroelectric Cc crystals can grow within a paraelectric ice VI matrix but may become unstable once a fraction of the matrix has become hydrogen-ordered. The reduction in dielectric constant causes P1̅ and other structures with small polarization to become favored, providing a possible resolution between observation and theoretical predictions.

Understanding molecular crystals with dispersion-inclusive density functional theory: pairwise corrections and beyond

CONSPECTUS: Molecular crystals are ubiquitous in many areas of science and engineering, including biology and medicine. Until recently, our ability to understand and predict their structure and properties using density functional theory was severely limited by the lack of approximate exchange-correlation functionals able to achieve sufficient accuracy. Here we show that there are many cases where the simple, minimally empirical pairwise correction scheme of Tkatchenko and Scheffler provides a useful prediction of the structure and properties of molecular crystals. After a brief introduction of the approach, we demonstrate its strength through some examples taken from our recent work. First, we show the accuracy of the approach using benchmark data sets of molecular complexes. Then we show its efficacy for structural determination using the hemozoin crystal, a challenging system possessing a wide range of strong and weak binding scenarios. Next, we show that it is equally useful for response properties by considering the elastic constants exhibited by the supramolecular diphenylalanine peptide solid and the infrared signature of water libration movements in brushite. Throughout, we emphasize lessons learned not only for the methodology but also for the chemistry and physics of the crystals in question. We further show that in many other scenarios where the simple pairwise correction scheme is not sufficiently accurate, one can go beyond it by employing a computationally inexpensive many-body dispersive approach that results in useful, quantitative accuracy, even in the presence of significant screening and/or multibody contributions to the dispersive energy. We explain the principles of the many-body approach and demonstrate its accuracy for benchmark data sets of small and large molecular complexes and molecular solids.

Thermodynamics and the hydrophobic effect in a core-softened model and comparison with experiments

A simple and computationally inexpensive core-softened model, originally proposed by Franzese [G. Franzese, J. Mol. Liq. 136, 267 (2007)], was adopted to show that it exhibits properties of waterlike fluid and hydrophobic effect. The potential used between particles is spherically symmetric with two characteristic lengths. Thermodynamics of nonpolar solvation were modeled as an insertion of a modified Lennard-Jones particle. It was investigated how the anomalous predictions of the model as well as the nonpolar solvation compare with the experimental data for water anomalies and the temperature dependence of noble gases hydration. It was shown that the model qualitatively follows the same trends as water. The model is able to reproduce waterlike anomalous properties (density maximum, heat capacity minimum, isothermal compressibility, etc.) and hydrophobic effect (minimum solubility for nonpolar solutes near ambient conditions, increased solubility of larger noble gases, etc.). It is argued that the model yields similar results as more complex and computationally expensive models.

Novel hydrogen hydrate structures under pressure

Gas hydrates are systems of prime importance. In particular, hydrogen hydrates are potential materials of icy satellites and comets, and may be used for hydrogen storage. We explore the H₂O–H₂ system at pressures in the range 0–100 GPa with ab initio variable-composition evolutionary simulations. According to our calculation and previous experiments, the H₂O–H₂ system undergoes a series of transformations with pressure, and adopts the known open-network clathrate structures (sII, C₀), dense “filled ice” structures (C₁, C₂) and two novel hydrate phases. One of these is based on the hexagonal ice framework and has the same H₂O:H₂ ratio (2:1) as the C₀ phase at low pressures and similar enthalpy (we name this phase Ih-C₀). The other newly predicted hydrate phase has a 1:2 H₂O:H₂ ratio and structure based on cubic ice. This phase (which we name C₃) is predicted to be thermodynamically stable above 38 GPa when including van der Waals interactions and zero-point vibrational energy, and explains previously mysterious experimental X-ray diffraction and Raman measurements. This is the hydrogen-richest hydrate and this phase has a remarkable gravimetric density (18 wt.%) of easily extractable hydrogen.