skip to main content


Title: Symmetry and 1 H NMR chemical shifts of short hydrogen bonds: impact of electronic and nuclear quantum effects
Short hydrogen bonds (SHBs), which have donor and acceptor separations below 2.7 Å, occur extensively in small molecules and proteins. Due to their compact structures, SHBs exhibit prominent covalent characters with elongated Donor–H bonds and highly downfield (>14 ppm) 1 H NMR chemical shifts. In this work, we carry out first principles simulations on a set of model molecules to assess how quantum effects determine the symmetry and chemical shift of their SHBs. From simulations that incorporate the quantum mechanical nature of both the electrons and nuclei, we reveal a universal relation between the chemical shift and the position of the proton in a SHB, and unravel the origin of the observed downfield spectral signatures. We further develop a metric that allows one to accurately and efficiently determine the proton position directly from its 1 H chemical shift, which will facilitate the experimental examination of SHBs in both small molecules and biological macromolecules.  more » « less
Award ID(s):
1904800
NSF-PAR ID:
10147864
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Physical Chemistry Chemical Physics
Volume:
22
Issue:
9
ISSN:
1463-9076
Page Range / eLocation ID:
4884 to 4895
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Despite its importance in electron transfer reactions and radiation chemistry, there has been disagreement over the fundamental nature of the hydrated electron, such as whether or not it resides in a cavity. Mixed quantum/classical simulations of the hydrated electron give different structures depending on the pseudopotential employed, and ab initio models of computational necessity use small numbers of water molecules and/or provide insufficient statistics to compare to experimental observables. A few years ago, Kumar et al. (J. Phys. Chem. A 2015, 119, 9148) proposed a minimalist ab initio model of the hydrated electron with only a small number of explicitly treated water molecules plus a polarizable continuum model (PCM). They found that the optimized geometry had four waters arranged tetrahedrally around a central cavity, and that the calculated vertical detachment energy and radius of gyration agreed well with experiment, results that were largely independent of the level of theory employed. The model, however, is based on a fixed structure at 0 K and does not explicitly incorporate entropic contributions or the thermal fluctuations that should be associated with the room-temperature hydrated electron. Thus, in this paper, we extend the model of Kumar et al. by running Born−Oppenheimer molecular dynamics (BOMD) of a small number of water molecules with an excess electron plus PCM at room temperature. We find that when thermal fluctuations are introduced, the level of theory chosen becomes critical enough when only four waters are used that one of the waters dissociates from the cluster with certain density functionals. Moreover, even with an optimally tuned range-separated hybrid functional, at room temperature the tetrahedral orientation of the 0 K first-shell waters is entirely lost and the central cavity collapses, a process driven by the fact that the explicit water molecules prefer to make H-bonds with each other more than with the excess electron. The resulting average structure is quite similar to that produced by a noncavity mixed quantum/classical model, so that the minimalist 4-water BOMD models suffer from problems similar to those of noncavity models, such as predicting the wrong sign of the hydrated electron’s molar solvation volume. We also performed BOMD with 16 explicit water molecules plus an extra electron and PCM. We find that the inclusion of an entire second solvation shell of explicit water leads to little change in the outcome from when only four waters were used. In fact, the 16-water simulations behave much like those of water cluster anions, in which the electron localizes at the cluster surface, showing that PCM is not acceptable for use in minimalist models to describe the behavior of the bulk hydrated electron. For both the 4- and 16-water models, we investigate how the introduction of thermal motions alters the predicted absorption spectrum, vertical detachment energy, and resonance Raman spectrum of the simulated hydrated electron. We also present a set of structural criteria that can be used to numerically determine how cavity-like (or not) a particular hydrated electron model is. All of the results emphasize that the hydrated electron is a statistical object whose properties are inadequately captured using only a small number of explicit waters, and that a proper treatment of thermal fluctuations is critical to understanding the hydrated electron’s chemical and physical behavior. 
    more » « less
  2. Molecular dynamics (MD) simulations were used to investigate the structure and lifetimes of hydrogen bonds and auto dissociation via proton transfer in bulk water using a reactive and dissociative all-atom potential that has previously been shown to match a variety of water properties and proton transfer. Using the topological model, each molecule's donated and accepted hydrogen bonds were labeled relative to the other hydrogen bonds on neighboring waters, providing a description of the effect of these details on the structure, dynamics and autoionization of water molecules. In agreement with prior data, asymmetric bonding at the sub-100 femtosecond timescale is observed, as well as the existence of linear, bifurcated, and dangling hydrogen bonds. The lifetime of the H-bond, 2.1 ps, is consistent with experimental data, with short time librations on the order of femtoseconds. The angular correlation functions, the presence of a second shell water entering the first shell, and OH vibrational stretch frequencies were all consistent with experiment or ab initio calculations. The simulations show short-lived (femtoseconds) dissociation of a small fraction of water molecules followed by rapid recombination. The role of the other H-bonds to the acceptor and on the donor plays an important part in proton transfer between the molecules in auto dissociation and is consistent with the role of a strong electric field caused by local (first and second shell) waters on initiating dissociation. The number of H-bonds to the donor water is 4.3 per molecule in the simulations, consistent with previous data regarding the number of hydrogen bonds required to generate this strong local electric field that enhances dissociation. The continuous lifetime autocorrelation function of the H-bond for those molecules that experience dissociation is considerably longer than that for all molecules that show no proton transfer. 
    more » « less
  3. Nuclear spin relaxation mechanisms are often difficult to isolate and identify, especially in molecules with internal flexibility. Here we combine experimental work with computation in order to determine the major mechanisms responsible for 31 P spin–lattice and singlet order (SO) relaxation in pyrophosphate, a physiologically relevant molecule. Using field-shuttling relaxation measurements (from 2 μT to 9.4 T) and rates calculated from molecular dynamics (MD) trajectories, we identified chemical shift anisotropy (CSA) and spin–rotation as the major mechanisms, with minor contributions from intra- and intermolecular coupling. The significant spin–rotation interaction is a consequence of the relatively rapid rotation of the –PO 3 2− entities around the bridging P–O bonds, and is treated by a combination of MD simulations and quantum chemistry calculations. Spin–lattice relaxation was predicted well without adjustable parameters, and for SO relaxation one parameter was extracted from the comparison between experiment and computation (a correlation coefficient between the rotational motion of the groups). 
    more » « less
  4. Abstract

    This study employs ab initio calculations based on density functional theory (DFT) to investigate the structural properties,1H-NMR spectra, and vibrational spectra of methane sulfonic acid (MSA) at low degree of hydration. The findings reveal that energetically stable structures are formed by small clusters consisting of one or two MSA molecules (m = 1 and 2) and one or two water molecules in (MSA)m·(H2O)n(m = 1–2 and n = 1–5).These stable structures arise from the formation of strong cyclic hydrogen bonds between the proton of the hydroxyl (OH) group in MSA and the water molecules. However, clusters containing three or more water molecules (n > 2) exhibit proton transfer from MSA to water, resulting in the formation of ion-pairs composed of CH3SO3and H3O+species. The measured1H-NMR spectra demonstrate the presence of hydrogen-bonded interactions between MSA and water, with a single MSA molecule interacting with water molecules. This interaction model accurately represents the hydrogen bonding network, as supported by the agreement between the experimental and calculated NMR chemical shift results.

     
    more » « less
  5. The absorption spectra of molecular organic chromophores in aqueous media are of considerable importance in environmental chemistry. In this work, the UV-vis spectra of benzoic acid (BA), the simplest aromatic carboxylic acid, in aqueous solutions at varying pH and in the presence of salts are measured experimentally. The solutions of different pH provide insights into the contributions from both the non-dissociated acid molecule and the deprotonated anionic species. The microscopic interpretation of these spectra is then provided by quantum chemical calculations for small cluster models of benzoic species (benzoic acid and benzoate anion) with water molecules. Calculations of the UV-vis absorbance spectra are then carried out for different clusters such as C 6 H 5 COOH·(H 2 O) n and C 6 H 5 COO − ·(H 2 O) n , where n = 0–8. The following main conclusions from these calculations and the comparison to experimental results can be made: (i) the small water cluster yields good quantitative agreement with observed solution experiments; (ii) the main peak position is found to be very similar at different levels of theory and is in excellent agreement with the experimental value, however, a weaker feature about 1 eV to lower energy (red shift) of the main peak is correctly reproduced only by using high level of theory, such as Algebraic Diagrammatic Construction (ADC); (iii) dissociation of the BA into ions is found to occur with a minimum of water molecules of n = 8; (iv) the deprotonation of BA has an influence on the computed spectrum and the energetics of the lowest energy electronic transitions; (v) the effect of the water on the spectra is much larger for the deprotonated species than for the non-dissociated acid. It was found that to reproduce experimental spectrum at pH 8.0, additional continuum representation for the extended solvent environment must be included in combination with explicit solvent molecules ( n ≥ 3); (vi) salts (NaCl and CaCl 2 ) have minimal effect on the absorption spectrum and; (vii) experimental results showed that B-band of neutral BA is not sensitive to the solvent effects whereas the effect of the water on the C-band is significant. The water effects blue-shift this band up to ∼0.2 eV. Overall, the results demonstrate the ability to further our understanding of the microscopic interpretation of the electronic structure and absorption spectra of BA in aqueous media through calculations restricted to small cluster models. 
    more » « less