Nuclear Quantum Effects and Hydrogen Bonding in Liquids

The Degradation of the 2-Hydroxyethylhydrazinium Nitrate Ionic Liquid System

The mechanism by which the 2-hydroxyethylhydrazinium nitrate (HEHN) ionic liquid (IL) dissociates via the formation and loss of nitric acid is investigated utilizing the effective fragment potential (EFP) method and various ab initio methods. Additionally, the interactions that dictate the stability of the HEHN IL and contribute to the formation of nitric acid are further assessed using the quasi-atomic orbital (QUAO) bonding analysis. From this work, it is suggested that the formation of nitric acid is dictated by the presence of hydrogen-bonding interactions which appear to become increasingly ionic in nature as the system approaches the bulk, limiting the formation of nitric acid.

Chemical features and machine learning assisted predictions of protein-ligand short hydrogen bonds

There are continuous efforts to elucidate the structure and biological functions of short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms reside more than 0.3 Å closer than the sum of their van der Waals radii. In this work, we evaluate 1070 atomic-resolution protein structures and characterize the common chemical features of SHBs formed between the side chains of amino acids and small molecule ligands. We then develop a machine learning assisted prediction of protein-ligand SHBs (MAPSHB-Ligand) model and reveal that the types of amino acids and ligand functional groups as well as the sequence of neighboring residues are essential factors that determine the class of protein-ligand hydrogen bonds. The MAPSHB-Ligand model and its implementation on our web server enable the effective identification of protein-ligand SHBs in proteins, which will facilitate the design of biomolecules and ligands that exploit these close contacts for enhanced functions.

Effective prediction of short hydrogen bonds in proteins via machine learning method

Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into account 21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.

Conditional Wave Function Theory: A Unified Treatment of Molecular Structure and Nonadiabatic Dynamics

We demonstrate that a conditional wave function theory enables a unified and efficient treatment of the equilibrium structure and nonadiabatic dynamics of correlated electron-ion systems. The conditional decomposition of the many-body wave function formally recasts the full interacting wave function of a closed system as a set of lower-dimensional (conditional) coupled "slices". We formulate a variational wave function ansatz based on a set of conditional wave function slices and demonstrate its accuracy by determining the structural and time-dependent response properties of the hydrogen molecule. We then extend this approach to include time-dependent conditional wave functions and address paradigmatic nonequilibrium processes including strong-field molecular ionization, laser-driven proton transfer, and nuclear quantum effects induced by a conical intersection. This work paves the road for the application of conditional wave function theory in equilibrium and out-of-equilibrium ab initio molecular simulations of finite and extended systems.

Artificial Neural Networks as Propagators in Quantum Dynamics

The utilization of artificial neural networks (ANNs) provides strategies for accelerating molecular simulations. Herein, ANNs are implemented as propagators of the time-dependent Schrödinger equation to simulate the quantum dynamics of systems with time-dependent potentials. These ANN propagators are trained to map nonstationary wavepackets from a given time to a future time within the discrete variable representation. Each propagator is trained for a specified time step, and iterative application of the propagator enables the propagation of wavepackets over long time scales. Such ANN propagators are developed and applied to one- and two-dimensional proton transfer systems, which exhibit nuclear quantum effects such as hydrogen tunneling. These ANN propagators are trained for either a specific time-independent potential or general potentials that can be time-dependent. Hierarchical, multiple time step algorithms enable parallelization, and the extension to higher dimensions is straightforward. This strategy is applicable to quantum dynamical simulations of diverse chemical and biological processes.

Artificial Neural Networks as Mappings between Proton Potentials, Wave Functions, Densities, and Energy Levels

Artificial neural networks (ANNs) have become important in quantum chemistry. Herein, applications to nuclear quantum effects, such as zero-point energy, vibrationally excited states, and hydrogen tunneling, are explored. ANNs are used to solve the time-independent Schrödinger equation for single- and double-well potentials representing hydrogen-bonded molecular systems capable of proton transfer. ANN mappings are trained to predict the lowest five proton vibrational energies, wave functions, and densities from the proton potentials and to predict the excited state proton vibrational energies and densities from the proton ground state density. For the inverse problem, ANN mappings are trained to predict the proton potential from the proton vibrational energy levels or the proton ground state density. This latter mapping is theoretically justified by the first Hohenberg-Kohn theorem establishing a one-to-one correspondence between the external potential and the ground state density. ANNs for two- and three-dimensional systems are also presented to illustrate the straightforward extension to higher dimensions.

Development of nuclear basis sets for multicomponent quantum chemistry methods

The nuclear-electronic orbital (NEO) framework provides a practical approach for directly incorporating nuclear quantum effects and non-Born-Oppenheimer effects of specified nuclei, typically protons, into quantum chemistry calculations. Multicomponent wave function based methods, such as NEO coupled cluster singles and doubles, and multicomponent density functional theory (DFT), such as NEO-DFT, require the appropriate selection of electronic and nuclear basis sets. Although a wide array of electronic basis sets are available, systematically developed nuclear basis sets that balance accuracy and efficiency have been lacking. Herein, a series of nuclear basis sets are developed and shown to be accurate and efficient for describing both ground and excited state properties of multicomponent systems in which electrons and specified protons are treated quantum mechanically. Three series of Gaussian-type nuclear basis sets, denoted PB4, PB5, and PB6, are developed with varying levels of angular momentum. A machine-learning optimization procedure relying on the Gaussian process regression method is utilized to accelerate the optimization process. The basis sets are validated in terms of predictions of ground state energies, proton densities, proton affinities, and proton vibrational excitation energies, allowing the user to select the desired balance between accuracy and efficiency for the properties of interest. These nuclear basis sets will enhance the tractability of NEO methods for applications to a wide range of chemical systems.

Quantum Stabilization of the Frustrated Hydrogen Bonding Structure in the Hydrogen Fluoride Trimer

We performed ab initio path integral molecular dynamics (PIMD) and molecular dynamics (MD) simulations to discuss the thermal and nuclear quantum effects on the stabilities of hydrogen bonding network in a hydrogen fluoride trimer (HF)₃ cluster. By the conventional molecular orbital calculation, the (HF)₃ cluster has an equilateral triangle shape, which has a frustration in the chemical structure of the hydrogen bonds, whereas the hydrogen bonding structure of a hydrogen fluoride dimer (HF)₂ cluster is nearly perpendicular to the acceptor molecule. The ratio of the triangular structures with the three hydrogen bondings in the PIMD simulation is larger than that in the MD one, whereas nonhydrogen bonding conformations such as a dimerlike structure are often found in MD simulation. The nuclear quantum effect stabilizes the frustrated hydrogen bonding network of the triangular (HF)₃ cluster.

Unraveling the structural and chemical features of biological short hydrogen bonds

Short hydrogen bonds are ubiquitous in biological macromolecules and exhibit distinctive proton potential energy surfaces and proton sharing properties.

The three-dimensional architecture of biomolecules often creates specialized structural elements, notably short hydrogen bonds that have donor–acceptor separations below 2.7 Å. In this work, we statistically analyze 1663 high-resolution biomolecular structures from the Protein Data Bank and demonstrate that short hydrogen bonds are prevalent in proteins, protein–ligand complexes and nucleic acids. From these biological macromolecules, we characterize the preferred location, connectivity and amino acid composition in short hydrogen bonds and hydrogen bond networks, and assess their possible functional importance. Using electronic structure calculations, we further uncover how the interplay of the structural and chemical features determines the proton potential energy surfaces and proton sharing conditions in biological short hydrogen bonds.

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.

Theoretical study of the H/D isotope effect on phase transition of hydrogen-bonded organic conductor κ-H3(Cat-EDT-TTF)2

κ-H₃(Cat-EDT-TTF)₂ (H-TTF) is a hydrogen-bonded π-electron system. Only its isotopologue, D-TTF, shows the phase transition. We obtained a symmetric single-well effective-PEC for H-TTF and low-barrier effective-PEC for D-TTF.

Theoretical analysis of geometry and NMR isotope shift in hydrogen-bonding center of photoactive yellow protein by combination of multicomponent quantum mechanics and ONIOM scheme

Multicomponent quantum mechanical (MC_QM) calculation has been extended with ONIOM (our own N-layered integrated molecular orbital + molecular mechanics) scheme [ONIOM(MC_QM:MM)] to take account of both the nuclear quantum effect and the surrounding environment effect. The authors have demonstrated the first implementation and application of ONIOM(MC_QM:MM) method for the analysis of the geometry and the isotope shift in hydrogen-bonding center of photoactive yellow protein. ONIOM(MC_QM:MM) calculation for a model with deprotonated Arg52 reproduced the elongation of O-H bond of Glu46 observed by neutron diffraction crystallography. Among the unique isotope shifts in different conditions, the model with protonated Arg52 with solvent effect reasonably provided the best agreement with the corresponding experimental values from liquid NMR measurement. Our results implied the availability of ONIOM(MC_QM:MM) to distinguish the local environment around hydrogen bonds in a biomolecule.

Ab initio molecular dynamics study of water at constant pressure using converged basis sets and empirical dispersion corrections

It is generally believed that studies of liquid water using the generalized gradient approximation to density functional theory require dispersion corrections in order to obtain reasonably accurate structural and dynamical properties. Here, we report on an ab initio molecular dynamics study of water in the isothermal-isobaric ensemble using a converged discrete variable representation basis set and an empirical dispersion correction due to Grimme [J. Comp. Chem. 27, 1787 (2006)]. At 300 K and an applied pressure of 1 bar, the density obtained without dispersion corrections is approximately 0.92 g/cm(3) while that obtained with dispersion corrections is 1.07 g/cm(3), indicating that the empirical dispersion correction overestimates the density by almost as much as it is underestimated without the correction for this converged basis. Radial distribution functions exhibit a loss of structure in the second solvation shell. Comparison of our results with other studies using the same empirical correction suggests the cause of the discrepancy: the Grimme dispersion correction is parameterized for use with a particular basis set; this parameterization is sensitive to this choice and, therefore, is not transferable to other basis sets.

Quantum nature of the hydrogen bond

Hydrogen bonds are weak, generally intermolecular bonds, which hold much of soft matter together as well as the condensed phases of water, network liquids, and many ferroelectric crystals. The small mass of hydrogen means that they are inherently quantum mechanical in nature, and effects such as zero-point motion and tunneling must be considered, though all too often these effects are not considered. As a prominent example, a clear picture for the impact of quantum nuclear effects on the strength of hydrogen bonds and consequently the structure of hydrogen bonded systems is still absent. Here, we report ab initio path integral molecular dynamics studies on the quantum nature of the hydrogen bond. Through a systematic examination of a wide range of hydrogen bonded systems we show that quantum nuclear effects weaken weak hydrogen bonds but strengthen relatively strong ones. This simple correlation arises from a competition between anharmonic intermolecular bond bending and intramolecular bond stretching. A simple rule of thumb is provided that enables predictions to be made for hydrogen bonded materials in general with merely classical knowledge (such as hydrogen bond strength or hydrogen bond length). Our work rationalizes the influence of quantum nuclear effects, which can result in either weakening or strengthening of the hydrogen bonds, and the corresponding structures, across a broad range of hydrogen bonded materials. Furthermore, it highlights the need to allow flexible molecules when anharmonic potentials are used in force field-based studies of quantum nuclear effects.

Effective force field for liquid hydrogen fluoride from ab initio molecular dynamics simulation using the force-matching method

A recently developed force-matching method for obtaining effective force fields for condensed matter systems from ab initio molecular dynamics (MD) simulations has been applied to fit a simple nonpolarizable two-site pairwise force field for liquid hydrogen fluoride. The ab initio MD in this case was a Car-Parrinello (CP) MD simulation of 64 HF molecules at nearly ambient conditions within the Becke-Lee-Yang-Parr approximation to the electronic density functional theory. The force-matching procedure included a fit of short-ranged nonbonded forces, bonded forces, and atomic partial charges. The performance of the force-match potential was examined for the gas-phase dimer and for the liquid phase at various temperatures. The model was able to reproduce correctly the bent structure and energetics of the gas-phase dimer, while the results for the structural properties, self-diffusion, vibrational spectra, density, and thermodynamic properties of liquid HF were compared to both experiment and the CP MD simulation. The force-matching model performs well in reproducing nearly all of the liquid properties as well as the aggregation behavior at different temperatures. The model is computationally cheap and compares favorably to many more computationally expensive potential energy functions for liquid HF.

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.

On the variation of the structure of liquid deuterium fluoride with temperature

The structure of liquid deuterium fluoride has been measured using pulsed neutron diffraction and high energy x-ray diffraction techniques as a function of temperature. The neutron experiments were performed at T=296+/-2 K, 246+/-2 K, and 193+/-2 K and the x-ray measurements carried out at 296+/-2 K and 195+/-2 K. The x-ray pair correlation functions, which are dominated by fluorine-fluorine interactions, show the first peak at approximately 2.53+/-0.05 A remains very nearly invariant with decreasing temperature. Peaks around 4.5 and 5.0 A also appear at both temperatures in the x-ray data. In contrast, the intermolecular peaks in the total neutron pair correlation function show that significant systematic local structural changes occur as the temperature is lowered. The first intermolecular peak position shortens from 1.64+/-0.05 A at 296 K to 1.56+/-0.05 A at 195 K. Although there are overlapping contributions from the intermolecular hydrogen-fluorine and hydrogen-hydrogen correlations, it is clear that the temperature dependent structural changes are largely due to a rearrangement of the deuterium atom positions in the fluid. By comparison with partial structure factor data the hydrogen bonds appear to become more linear at lower temperatures.

The World of Non-Covalent Interactions: 2006

The review focusses on the fundamental importance of non-covalent interactions in nature by illustrating specific examples from chemistry, physics and the biosciences. Laser spectroscopic methods and both ab initio and molecular modelling procedures used for the study of non-covalent interactions in molecular clusters are briefly outlined. The role of structure and geometry, stabilization energy, potential and free energy surfaces for molecular clusters is extensively discussed in the light of the most advanced ab initio computational results for the CCSD(T) method, extrapolated to the CBS limit. The most important types of non-covalent complexes are classified and several small and medium size non-covalent systems, including H-bonded and improper H-bonded complexes, nucleic acid base pairs, and peptides and proteins are discussed with some detail. Finally, we evaluate the interpretation of experimental results in comparison with state of the art theoretical models: this is illustrated for phenol...Ar, the benzene dimer and nucleic acid base pairs. A review with 270 references.

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.