Quantum nature of the hydrogen bond

Nuclear Quantum Effects in Water Reorientation and Hydrogen-Bond Dynamics

We combine classical and ring polymer molecular dynamics simulations with the molecular jump model to provide a molecular description of the nuclear quantum effects (NQEs) on water reorientation and hydrogen-bond dynamics in liquid H₂O and D₂O. We show that while the net NQE is negligible in D₂O, it leads to a ∼13% acceleration in H₂O dynamics compared to a classical description. Large angular jumps-exchanging hydrogen-bond partners-are the dominant reorientation pathway (just as in a classical description); the faster reorientation dynamics arise from the increased jump rate constant. NQEs do not change the jump amplitude distribution, and no significant tunneling is found. The faster jump dynamics are quantitatively related to decreased structuring of the OO radial distribution function when NQEs are included. This is explained, via a jump model analysis, by competition between the effects of water's librational and OH stretch mode zero-point energies on the hydrogen-bond strength.

The Rosetta All-Atom Energy Function for Macromolecular Modeling and Design

Over the past decade, the Rosetta biomolecular modeling suite has informed diverse biological questions and engineering challenges ranging from interpretation of low-resolution structural data to design of nanomaterials, protein therapeutics, and vaccines. Central to Rosetta's success is the energy function: a model parametrized from small-molecule and X-ray crystal structure data used to approximate the energy associated with each biomolecule conformation. This paper describes the mathematical models and physical concepts that underlie the latest Rosetta energy function, called the Rosetta Energy Function 2015 (REF15). Applying these concepts, we explain how to use Rosetta energies to identify and analyze the features of biomolecular models. Finally, we discuss the latest advances in the energy function that extend its capabilities from soluble proteins to also include membrane proteins, peptides containing noncanonical amino acids, small molecules, carbohydrates, nucleic acids, and other macromolecules.

Hidden role of intermolecular proton transfer in the anomalously diffuse vibrational spectrum of a trapped hydronium ion

We report the vibrational spectra of the hydronium and methyl-ammonium ions captured in the C_3v binding pocket of the 18-crown-6 ether ionophore. Although the NH stretching bands of the CH₃NH₃⁺ ion are consistent with harmonic expectations, the OH stretching bands of H₃O⁺ are surprisingly broad, appearing as a diffuse background absorption with little intensity modulation over 800 cm^-1 with an onset ∼400 cm^-1 below the harmonic prediction. This structure persists even when only a single OH group is present in the HD₂O⁺ isotopologue, while the OD stretching region displays a regular progression involving a soft mode at about 85 cm^-1 These results are rationalized in a vibrationally adiabatic (VA) model in which the motion of the H₃O⁺ ion in the crown pocket is strongly coupled with its OH stretches. In this picture, H₃O⁺ resides in the center of the crown in the vibrational zero-point level, while the minima in the VA potentials associated with the excited OH vibrational states are shifted away from the symmetrical configuration displayed by the ground state. Infrared excitation between these strongly H/D isotope-dependent VA potentials then accounts for most of the broadening in the OH stretching manifold. Specifically, low-frequency motions involving concerted motions of the crown scaffold and the H₃O⁺ ion are driven by a Franck-Condon-like mechanism. In essence, vibrational spectroscopy of these systems can be viewed from the perspective of photochemical interconversion between transient, isomeric forms of the complexes corresponding to the initial stage of intermolecular proton transfer.

Quantum Dynamics and Spectroscopy of Ab Initio Liquid Water: The Interplay of Nuclear and Electronic Quantum Effects

Understanding the reactivity and spectroscopy of aqueous solutions at the atomistic level is crucial for the elucidation and design of chemical processes. However, the simulation of these systems requires addressing the formidable challenges of treating the quantum nature of both the electrons and nuclei. Exploiting our recently developed methods that provide acceleration by up to 2 orders of magnitude, we combine path integral simulations with on-the-fly evaluation of the electronic structure at the hybrid density functional theory level to capture the interplay between nuclear quantum effects and the electronic surface. Here we show that this combination provides accurate structure and dynamics, including the full infrared and Raman spectra of liquid water. This allows us to demonstrate and explain the failings of lower-level density functionals for dynamics and vibrational spectroscopy when the nuclei are treated quantum mechanically. These insights thus provide a foundation for the reliable investigation of spectroscopy and reactivity in aqueous environments.

Weak Supramolecular Interactions Governing Parallel and Antiparallel DNA Quadruplexes: Insights from Large-Scale Quantum Mechanics Analysis of Experimentally Derived Models

The topology and energetics of guanine (G) quadruplexes is governed by supramolecular interactions within their strands. In this work, an extensive quantum mechanical (QM) study has been performed to analyze supramolecular interactions that shape the stems of (4+0) parallel (P) and (2+2) antiparallel (AP) quadruplex systems. The large-scale (≈400 atoms) models of P and AP were constructed from high-quality experimental structures. The results provide evidence that each of the P and AP structures is shaped by a distinct network of supramolecular interactions. Analysis of electron topological characteristics of hydrogen bonds in P and AP systems indicates that the P model benefits from stronger intratetrad hydrogen bonding. For intertetrad stacking interactions, both noncovalent interaction plot and energy decomposition analysis approaches suggest that the stem of the P quadruplex benefits more from stacking than that of the AP stem; the difference in energetic stabilization for the two topologies is about 10 %. Stronger hydrogen-bonding and stacking interactions in the stem of the P quadruplex, relative to those in the AP system, can be an important indicator to explain the experimental observations that guanine-rich oligonucleotides tend to form all-parallel stems with an all-anti orientation of nucleobases. However, in addition to intrinsic stabilization, partial desolvation effects, which affect the energetics and dynamics of the G-quadruplex folding process, call for further investigations.

Nuclear quantum effects in a HIV/cancer inhibitor: The case of ellipticine

Ellipticine is a natural product that is currently being actively investigated for its inhibitory cancer and HIV properties. Here we use path-integral molecular dynamics coupled with excited state calculations to characterize the role of nuclear quantum effects on the structural and electronic properties of ellipticine in water, a common biological solvent. Quantum effects collectively enhance the fluctuations of both light and heavy nuclei of the covalent and hydrogen bonds in ellipticine. In particular, for the ellipticine-water system, where the proton donor and acceptor have different proton affinities, we find that nuclear quantum effects (NQEs) strengthen both the strong and the weak H bonds. This is in contrast to what is observed for the cases where the proton affinity of the donors and acceptors is same. These structural fluctuations cause a significant red-shift in the absorption spectra and an increase in the broadening, bringing it into closer agreement with the experiments. Our work shows that nuclear quantum effects alter both qualitatively and quantitatively the optical properties of this biologically relevant system and highlights the importance of the inclusion of these effects in the microscopic understanding of their optical properties. We propose that isotopic substitution will produce a blue shift and a reduction in the broadening of the absorption peak.

Perspective: Structure and dynamics of water at surfaces probed by scanning tunneling microscopy and spectroscopy

The detailed and precise understanding of water-solid interaction largely relies on the development of atomic-scale experimental techniques, among which scanning tunneling microscopy (STM) has proven to be a noteworthy example. In this perspective, we review the recent advances of STM techniques in imaging, spectroscopy, and manipulation of water molecules. We discuss how those newly developed techniques are applied to probe the structure and dynamics of water at solid surfaces with single-molecule and even submolecular resolution, paying particular attention to the ability of accessing the degree of freedom of hydrogen. In the end, we present an outlook on the directions of future STM studies of water-solid interfaces as well as the challenges faced by this field. Some new scanning probe techniques beyond STM are also envisaged.

Collective proton transfer in ordinary ice: local environments, temperature dependence and deuteration effects

Ordinary ice at low temperature: what about collective nuclear quantum effects in its chiral six rings?

Symmetry and dynamics of FHF− anion in vacuum, in CD2Cl2 and in CCl4. Ab initio MD study of fluctuating solvent–solute hydrogen and halogen bonds

Asymmetric solvation of FHF⁻ by halogen- and hydrogen-bonding solvents breaks the symmetry of the anion.

Different Ways of Hydrogen Bonding in Water - Why Does Warm Water Freeze Faster than Cold Water?

The properties of liquid water are intimately related to the H-bond network among the individual water molecules. Utilizing vibrational spectroscopy and modeling water with DFT-optimized water clusters (6-mers and 50-mers), 16 out of a possible 36 different types of H-bonds are identified and ordered according to their intrinsic strength. The strongest H-bonds are obtained as a result of a concerted push-pull effect of four peripheral water molecules, which polarize the electron density in a way that supports charge transfer and partial covalent character of the targeted H-bond. For water molecules with tetra- and pentacoordinated O atoms, H-bonding is often associated with a geometrically unfavorable positioning of the acceptor lone pair and donor σ^*(OH) orbitals so that electrostatic rather than covalent interactions increasingly dominate H-bonding. There is a striking linear dependence between the intrinsic strength of H-bonding as measured by the local H-bond stretching force constant and the delocalization energy associated with charge transfer. Molecular dynamics simulations for 1000-mers reveal that with increasing temperature weak, preferentially electrostatic H-bonds are broken, whereas the number of strong H-bonds increases. An explanation for the question why warm water freezes faster than cold water is given on a molecular basis.

Proton Probability Distribution in the O···H···O Low-Barrier Hydrogen Bond: A Combined Solid-State NMR and Quantum Chemical Computational Study of Dibenzoylmethane and Curcumin

We report a combined solid-state (¹H, ²H, ¹³C, ¹⁷O) NMR and plane-wave density functional theory (DFT) computational study of the O···H···O low-barrier hydrogen bonds (LBHBs) in two 1,3-diketone compounds: dibenzoylmethane (1) and curcumin (2). In the solid state, both 1 and 2 exist in the cis-keto-enol tautomeric form, each exhibiting an intramolecular LBHB with a short O···O distance (2.435 Å in 1 and 2.455 Å in 2). Whereas numerous experimental (structural and spectroscopic) and computational studies have been reported for the enol isomers of 1,3-diketones, a unified picture about the proton location within an LBHB is still lacking. This work reports for the first time the solid-state ¹⁷O NMR data for the O···H···O LBHBs in 1,3-diketones. The central conclusion of this work is that detailed information about the probability density distribution of the proton (nuclear zero-point motion) across an LBHB can be obtained from a combination of solid-state NMR and plane-wave DFT computations (both NMR parameter calculations and ab initio molecular dynamics simulations). We propose that the precise proton probability distribution across an LBHB should provide a common basis on which different and sometimes seemingly contradicting experimental results obtained from complementary techniques, such as X-ray diffraction, neutron diffraction, and solid-state NMR, can be reconciled.

Exact Factorization-Based Density Functional Theory of Electrons and Nuclei

The ground state energy of a system of electrons (r=r_{1},r_{2},…) and nuclei (R=R_{1},R_{2},…) is proven to be a variational functional of the electronic density n(r,R) and paramagnetic current density j_{p}(r,R) conditional on R, the nuclear wave function χ(R), an induced vector potential A_{μ}(R) and a quantum geometric tensor T_{μν}(R). n, j_{p}, A_{μ} and T_{μν} are defined in terms of the conditional electronic wave function Φ_{R}(r). The ground state (n,j_{p},χ,A_{μ},T_{μν}) can be calculated by solving self-consistently (i) conditional Kohn-Sham equations containing effective scalar and vector potentials v_{s}(r) and A_{xc}(r) that depend parametrically on R, (ii) the Schrödinger equation for χ(R), and (iii) Euler-Lagrange equations that determine T_{μν}. The theory is applied to the E⊗e Jahn-Teller model.

Proton Conductivity in Phosphoric Acid: The Role of Quantum Effects

Phosphoric acid has one of the highest intrinsic proton conductivities of any known liquids, and the mechanism of this exceptional conductivity remains a puzzle. Our detailed experimental studies discovered a strong isotope effect in the conductivity of phosphoric acids caused by (i) a strong isotope shift of the glass transition temperature and (ii) a significant reduction of the energy barrier by zero-point quantum fluctuations. These results suggest that the high conductivity in phosphoric acids is caused by a very efficient proton transfer mechanism, which is strongly assisted by quantum effects.

Anharmonic and Quantum Fluctuations in Molecular Crystals: A First-Principles Study of the Stability of Paracetamol

Molecular crystals often exist in multiple competing polymorphs, showing significantly different physicochemical properties. Computational crystal structure prediction is key to interpret and guide the search for the most stable or useful form, a real challenge due to the combinatorial search space, and the complex interplay of subtle effects that work together to determine the relative stability of different structures. Here we take a comprehensive approach based on different flavors of thermodynamic integration in order to estimate all contributions to the free energies of these systems with density-functional theory, including the oft-neglected anharmonic contributions and nuclear quantum effects. We take the two main stable forms of paracetamol as a paradigmatic example. We find that anharmonic contributions, different descriptions of van der Waals interactions, and nuclear quantum effects all matter to quantitatively determine the stability of different phases. Our analysis highlights the many challenges inherent in the development of a quantitative and predictive framework to model molecular crystals. However, it also indicates which of the components of the free energy can benefit from a cancellation of errors that can redeem the predictive power of approximate models, and suggests simple steps that could be taken to improve the reliability of ab initio crystal structure prediction.

Protons and Hydroxide Ions in Aqueous Systems

Understanding the structure and dynamics of water's constituent ions, proton and hydroxide, has been a subject of numerous experimental and theoretical studies over the last century. Besides their obvious importance in acid-base chemistry, these ions play an important role in numerous applications ranging from enzyme catalysis to environmental chemistry. Despite a long history of research, many fundamental issues regarding their properties continue to be an active area of research. Here, we provide a review of the experimental and theoretical advances made in the last several decades in understanding the structure, dynamics, and transport of the proton and hydroxide ions in different aqueous environments, ranging from water clusters to the bulk liquid and its interfaces with hydrophobic surfaces. The propensity of these ions to accumulate at hydrophobic surfaces has been a subject of intense debate, and we highlight the open issues and challenges in this area. Biological applications reviewed include proton transport along the hydration layer of various membranes and through channel proteins, problems that are at the core of cellular bioenergetics.

Nuclear Quantum Effects in Water at the Triple Point: Using Theory as a Link Between Experiments

One of the most prominent consequences of the quantum nature of light atomic nuclei is that their kinetic energy does not follow a Maxwell-Boltzmann distribution. Deep inelastic neutron scattering (DINS) experiments can measure this effect. Thus, the nuclear quantum kinetic energy can be probed directly in both ordered and disordered samples. However, the relation between the quantum kinetic energy and the atomic environment is a very indirect one, and cross-validation with theoretical modeling is therefore urgently needed. Here, we use state of the art path integral molecular dynamics techniques to compute the kinetic energy of hydrogen and oxygen nuclei in liquid, solid, and gas-phase water close to the triple point, comparing three different interatomic potentials and validating our results against equilibrium isotope fractionation measurements. We will then show how accurate simulations can draw a link between extremely precise fractionation experiments and DINS, therefore establishing a reliable benchmark for future measurements and providing key insights to increase further the accuracy of interatomic potentials for water.

Inverse Temperature Dependence of Nuclear Quantum Effects in DNA Base Pairs

Despite the inherently quantum mechanical nature of hydrogen bonding, it is unclear how nuclear quantum effects (NQEs) alter the strengths of hydrogen bonds. With this in mind, we use ab initio path integral molecular dynamics to determine the absolute contribution of NQEs to the binding in DNA base pair complexes, arguably the most important hydrogen-bonded systems of all. We find that depending on the temperature, NQEs can either strengthen or weaken the binding within the hydrogen-bonded complexes. As a somewhat counterintuitive consequence, NQEs can have a smaller impact on hydrogen bond strengths at cryogenic temperatures than at room temperature. We rationalize this in terms of a competition of NQEs between low-frequency and high-frequency vibrational modes. Extending this idea, we also propose a simple model to predict the temperature dependence of NQEs on hydrogen bond strengths in general.

Nuclear quantum effects of hydrogen bonds probed by tip-enhanced inelastic electron tunneling

We report the quantitative assessment of nuclear quantum effects on the strength of a single hydrogen bond formed at a water-salt interface, using tip-enhanced inelastic electron tunneling spectroscopy based on a scanning tunneling microscope. The inelastic scattering cross section was resonantly enhanced by "gating" the frontier orbitals of water via a chlorine-terminated tip, so the hydrogen-bonding strength can be determined with high accuracy from the red shift in the oxygen-hydrogen stretching frequency of water. Isotopic substitution experiments combined with quantum simulations reveal that the anharmonic quantum fluctuations of hydrogen nuclei weaken the weak hydrogen bonds and strengthen the relatively strong ones. However, this trend can be completely reversed when a hydrogen bond is strongly coupled to the polar atomic sites of the surface.

Electrolytes induce long-range orientational order and free energy changes in the H-bond network of bulk water

Electrolytes interact with water in many ways: changing dipole orientation, inducing charge transfer, and distorting the hydrogen-bond network in the bulk and at interfaces. Numerous experiments and computations have detected short-range perturbations that extend up to three hydration shells around individual ions. We report a multiscale investigation of the bulk and surface of aqueous electrolyte solutions that extends from the atomic scale (using atomistic modeling) to nanoscopic length scales (using bulk and interfacial femtosecond second harmonic measurements) to the macroscopic scale (using surface tension experiments). Electrolytes induce orientational order at concentrations starting at 10 μM that causes nonspecific changes in the surface tension of dilute electrolyte solutions. Aside from ion-dipole interactions, collective hydrogen-bond interactions are crucial and explain the observed difference of a factor of 6 between light water and heavy water.

Textured fluorapatite bonded to calcium sulphate strengthen stomatopod raptorial appendages

Stomatopods are shallow-water crustaceans that employ powerful dactyl appendages to hunt their prey. Deployed at high velocities, these hammer-like clubs or spear-like devices are able to inflict substantial impact forces. Here we demonstrate that dactyl impact surfaces consist of a finely-tuned mineral gradient, with fluorapatite substituting amorphous apatite towards the outer surface. Raman spectroscopy measurements show that calcium sulphate, previously not reported in mechanically active biotools, is co-localized with fluorapatite. Ab initio computations suggest that fluorapatite/calcium sulphate interfaces provide binding stability and promote the disordered-to-ordered transition of fluorapatite. Nanomechanical measurements show that fluorapatite crystalline orientation correlates with an anisotropic stiffness response and indicate significant differences in the fracture tolerance between the two types of appendages. Our findings shed new light on the crystallochemical and microstructural strategies allowing these intriguing biotools to optimize impact forces, providing physicochemical information that could be translated towards the synthesis of impact-resistant functional materials and coatings.

Nuclear quantum effects in water exchange around lithium and fluoride ions

We employ classical and ring polymer molecular dynamics simulations to study the effect of nuclear quantum fluctuations on the structure and the water exchange dynamics of aqueous solutions of lithium and fluoride ions. While we obtain reasonably good agreement with experimental data for solutions of lithium by augmenting the Coulombic interactions between the ion and the water molecules with a standard Lennard-Jones ion-oxygen potential, the same is not true for solutions of fluoride, for which we find that a potential with a softer repulsive wall gives much better agreement. A small degree of destabilization of the first hydration shell is found in quantum simulations of both ions when compared with classical simulations, with the shell becoming less sharply defined and the mean residence time of the water molecules in the shell decreasing. In line with these modest differences, we find that the mechanisms of the exchange processes are unaffected by quantization, so a classical description of these reactions gives qualitatively correct and quantitatively reasonable results. We also find that the quantum effects in solutions of lithium are larger than in solutions of fluoride. This is partly due to the stronger interaction of lithium with water molecules, partly due to the lighter mass of lithium and partly due to competing quantum effects in the hydration of fluoride, which are absent in the hydration of lithium.

Intramolecular structure and energetics in supercooled water down to 255 K

We studied the structure and energetics of supercooled water by means of X-ray Raman and Compton scattering. Under supercooled conditions down to 255 K, the oxygen K-edge measured by X-ray Raman scattering suggests an increase of tetrahedral order similar to the conventional temperature effect observed in non-supercooled water. Compton profile differences indicate contributions beyond the theoretically predicted temperature effect and provide a deeper insight into local structural changes. These contributions suggest a decrease of the electron mean kinetic energy by 3.3 ± 0.7 kJ (mol K)(-1) that cannot be modeled within established water models. Our surprising results emphasize the need for water models that capture in detail the intramolecular structural changes and quantum effects to explain this complex liquid.

Intramolecular H⋯S interactions in metal di-(isopropyl)dithiocarbamate complexes

Intramolecular C–H⋯S interactions create restricted rotation of groups within di(isopropyl)dithiocarbate complexes.

Going clean: structure and dynamics of peptides in the gas phase and paths to solvation

The gas phase is an artificial environment for biomolecules that has gained much attention both experimentally and theoretically due to its unique characteristic of providing a clean room environment for the comparison between theory and experiment. In this review we give an overview mainly on first-principles simulations of isolated peptides and the initial steps of their interactions with ions and solvent molecules: a bottom up approach to the complexity of biological environments. We focus on the accuracy of different methods to explore the conformational space, the connections between theory and experiment regarding collision cross section evaluations and (anharmonic) vibrational spectra, and the challenges faced in this field.

Stability of Complex Biomolecular Structures: van der Waals, Hydrogen Bond Cooperativity, and Nuclear Quantum Effects

Biomolecules are complex systems stabilized by a delicate balance of weak interactions, making it important to assess all energetic contributions in an accurate manner. However, it is a priori unclear which contributions make more of an impact. Here, we examine stacked polyglutamine (polyQ) strands, a peptide repeat often found in amyloid aggregates. We investigate the role of hydrogen bond (HB) cooperativity, van der Waals (vdW) dispersion interactions, and quantum contributions to free energies, including anharmonicities through density functional theory and ab initio path integral simulations. Of these various factors, we find that the largest impact on structural stabilization comes from vdW interactions. HB cooperativity is the second largest contribution as the size of the stacked chain grows. Competing nuclear quantum effects make the net quantum contribution small but very sensitive to anharmonicities, vdW, and the number of HBs. Our results suggest that a reliable treatment of these systems can only be attained by considering all of these components.

Modeling quantum nuclei with perturbed path integral molecular dynamics

Here we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics for modeling quantum nuclear effects.

The quantum nature of nuclear motions plays a vital role in the structure, stability, and thermodynamics of molecules and materials. The standard approach to model nuclear quantum fluctuations in chemical and biological systems is to use path-integral molecular dynamics. Unfortunately, conventional path-integral simulations can have an exceedingly large computational cost due to the need to employ an excessive number of coupled classical subsystems (beads) for quantitative accuracy. Here, we combine perturbation theory with the Feynman–Kac imaginary-time path integral approach to quantum mechanics and derive an improved non-empirical partition function and estimators to calculate converged quantum observables. Our perturbed path-integral (PPI) method requires the same ingredients as the conventional approach, but increases the accuracy and efficiency of path integral simulations by an order of magnitude. Results are presented for the thermodynamics of fundamental model systems, an empirical water model containing 256 water molecules within periodic boundary conditions, and ab initio simulations of nitrogen and benzene molecules. For all of these examples, PPI simulations with 4 to 8 classical beads recover the nuclear quantum contribution to the total energy and heat capacity at room temperature within a 3% accuracy, paving the way toward seamless modeling of nuclear quantum effects in realistic molecules and materials.

Theoretical vibrational spectra of OH(-)(H2O)2: the effect of quantum distribution and vibrational coupling

We performed ab initio path integral molecular dynamics simulations for the hydroxide-water cluster, OH(-)(H2O)2, at 50 K, 100 K, and 150 K to investigate its flexible structure. From our simulations, we found that nuclear quantum effects enhance hydroxide hydrogen atom inversion and the conformational change between isomers occurs by simultaneous rotation of the free hydrogen atom. We propose the importance of including the transition state conformer with C2 symmetry, for the description of this system at temperatures realized in predissociation experiments. Temperature dependence of relative populations of each conformer along with multidimensional vibrational calculations were used to simulate the vibrational spectra and compare with the experimental spectra of Johnson and coworkers. We assign the doublet peaks seen in the experiment at 2500 to 3000 cm(-1), as the mixture of the ionic hydrogen bonded OH stretching overtone, ionic hydrogen bonded OH bending overtone, and the combination band of the ionic hydrogen bonded OH stretch and bend, which are modulated by the van der Waals OO vibrations. We concluded that for OH(-)(H2O)2, the vibrational couplings between the ionic hydrogen bonded motion and floppy modes contribute to the broadening of peaks observed in the 2500 to 3000 cm(-1) region.

Strong frequency dependence of vibrational relaxation in bulk and surface water reveals sub-picosecond structural heterogeneity

Because of strong hydrogen bonding in liquid water, intermolecular interactions between water molecules are highly delocalized. Previous two-dimensional infrared spectroscopy experiments have indicated that this delocalization smears out the structural heterogeneity of neat H2O. Here we report on a systematic investigation of the ultrafast vibrational relaxation of bulk and interfacial water using time-resolved infrared and sum-frequency generation spectroscopies. These experiments reveal a remarkably strong dependence of the vibrational relaxation time on the frequency of the OH stretching vibration of liquid water in the bulk and at the air/water interface. For bulk water, the vibrational relaxation time increases continuously from 250 to 550 fs when the frequency is increased from 3,100 to 3,700 cm(-1). For hydrogen-bonded water at the air/water interface, the frequency dependence is even stronger. These results directly demonstrate that liquid water possesses substantial structural heterogeneity, both in the bulk and at the surface.

Enhancing the hydrogen bond between the bridged hydrogen atom of diborane and ammonia

The character of the bridged hydrogen atom (Hb) of B2H6 has become a hot issue in recent years. In this work, the complexes B2H6 · · · NH3, B2H2X4 · · · nNH3 (n = 1, 2) and 2HF · · · B2H2X4 · · · 2NH3 (X = Cl, Br, I) were constructed and studied based on the M06-2X calculations to investigate how to enhance the Hb · · · N hydrogen-bonded interaction. When the terminal hydrogen atoms (Ht) of B2H6 were replaced by X (X = Cl, Br, I) atoms, the Hb · · · N hydrogen bond were strengthened. According to the electrostatic potentials in B2H2X4, two HF molecules were added to the interspace of the B-H-B-H four-membered ring of the B2H2X4 · · · 2NH3 complexes, and H · · · X hydrogen bond formed, resulting in further enhancing effect of Hb · · · N hydrogen bond. As a result, the positive cooperative effect of Hb · · · N hydrogen bond and H · · · X hydrogen bond do enhance the interactions of each other. The two measures not only enhance the strength of Hb · · · N hydrogen bond, but also achieve the goal to make the Hb · · · N hydrogen bond perpendicular to B · · · B direction.

Potential Paths for the Hydrogen-Bond Relaxing With (H2O)N Cluster Size

Relaxation of the inter- and intra-molecular interactions for the hydrogen bond (O:H-O) between undercoordinated molecules determines the unusual behavior of water nanodroplets and nanobubbles. However, probing such potentials remains unreality. Here we show that the Lagrangian solution [Huang et al., J. Phys. Chem. B, 2013. 117: 13639] transforms the observed H-O bond (x = H) and O:H nonbond (x = L) lengths and phonon frequencies (dx, x) [Sun et al., J. Phys. Chem. Lett., 2013. 4: 2565] into the respective force constants and bond energies (kx, Ex) and hence enables the mapping of the potential paths for the O:H-O bond relaxing with water cluster size. Results show that molecular undercoordination not only reduces the molecular size (dH) with enhanced H-O energy from the bulk value of 3.97 to 5.10 eV for a H2O monomer, but also enlarges the molecular separation (dL) with reduced O:H energy from 95 to 35 meV for a dimer. The H-O energy gain raises the melting point from bulk value 273 to 310 K for the skin and the O:H energy loss lowers the freezing temperature from bulk value 258 to 202 K for 1.4 nm sized droplet, by dispersing the quasisolid phase boundaries.

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.

Ab initio molecular dynamics simulation of liquid water by quantum Monte Carlo

Although liquid water is ubiquitous in chemical reactions at roots of life and climate on the earth, the prediction of its properties by high-level ab initio molecular dynamics simulations still represents a formidable task for quantum chemistry. In this article, we present a room temperature simulation of liquid water based on the potential energy surface obtained by a many-body wave function through quantum Monte Carlo (QMC) methods. The simulated properties are in good agreement with recent neutron scattering and X-ray experiments, particularly concerning the position of the oxygen-oxygen peak in the radial distribution function, at variance of previous density functional theory attempts. Given the excellent performances of QMC on large scale supercomputers, this work opens new perspectives for predictive and reliable ab initio simulations of complex chemical systems.

Strong isotope effects on melting dynamics and ice crystallisation processes in cryo vitrification solutions

The nucleation and growth of crystalline ice during cooling, and further crystallization processes during re-warming are considered to be key processes determining the success of low temperature storage of biological objects, as used in medical, agricultural and nature conservation applications. To avoid these problems a method, termed vitrification, is being developed to inhibit ice formation by use of high concentration of cryoprotectants and ultra-rapid cooling, but this is only successful across a limited number of biological objects and in small volume applications. This study explores physical processes of ice crystal formation in a model cryoprotective solution used previously in trials on vitrification of complex biological systems, to improve our understanding of the process and identify limiting biophysical factors. Here we present results of neutron scattering experiments which show that even if ice crystal formation has been suppressed during quench cooling, the water molecules, mobilised during warming, can crystallise as detectable ice. The crystallisation happens right after melting of the glass phase formed during quench cooling, whilst the sample is still transiting deep cryogenic temperatures. We also observe strong water isotope effects on ice crystallisation processes in the cryoprotectant mixture. In the neutron scattering experiment with a fully protiated water component, we observe ready crystallisation occurring just after the glass melting transition. On the contrary with a fully deuteriated water component, the process of crystallisation is either completely or substantially supressed. This behaviour might be explained by nuclear quantum effects in water. The strong isotope effect, observed here, may play an important role in development of new cryopreservation strategies.

“On-the-fly” coupled cluster path-integral molecular dynamics: impact of nuclear quantum effects on the protonated water dimer

“On-the-fly” coupled cluster-based path-integral molecular dynamics simulations predict that the effective potential of the protonated water–dimer has a single-well only.

Hydride ion formation in stoichiometric UO2

We investigated hydrogen solubility in UO₂ using DFT and predicted that hydrogen species energetically prefers to exist as a hydride ion rather than a proton in a hydroxyl group and on diffusion hydrogen's charge state will change.

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₃.

Intermolecular OHN hydrogen bond with a proton moving in 3-methylpyridinium 2,6-dichloro-4-nitrophenolate

Temperature-dependent changes in the strong OHN hydrogen bond in 3-methylpyridinium 2,6-dichloro-4-nitrophenolate are used to discuss the proton transfer mechanism.

Facet-dependent photocatalytic and antibacterial properties of α-Ag2WO4crystals: combining experimental data and theoretical insights

We have combined experimental results and calculations with new paths to explain the photocatalytic and antibacterial activities of α-Ag₂WO₄crystals.