Temperature dependent deuterium quadrupole coupling constants of short hydrogen bonds

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.

Temperature Dependence of NMR Parameters Calculated from Path Integral Molecular Dynamics Simulations

The influence of temperature on NMR chemical shifts and quadrupolar couplings in model molecular organic solids is explored using path integral molecular dynamics (PIMD) and density functional theory (DFT) calculations of shielding and electric field gradient (EFG) tensors. An approach based on convoluting calculated shielding or EFG tensor components with probability distributions of selected bond distances and valence angles obtained from DFT-PIMD simulations at several temperatures is used to calculate the temperature effects. The probability distributions obtained from the quantum PIMD simulations, which includes nuclear quantum effects, are significantly broader and less temperature dependent than those obtained with conventional DFT molecular dynamics or with 1D scans through the potential energy surface. Predicted NMR observables for the model systems were in excellent agreement with experimental data.

Deuterium quadrupolar tensors of L-histidine hydrochloride monohydrate-d7

A deuterium NMR study at 14 T of a single crystal of L-histidine hydrochloride monohydrate has determined the deuteron quadrupole coupling constants CQ, asymmetry parameters eta, and electric field gradient orientation for the imidazolium and primary ammonium groups and for a water of crystallization. The imidazolium deuterons, which have very long relaxation times, have quite different coupling constants, reflecting different hydrogen bonding but nearly identical orientations, with the most distinct principal axis in both cases nearly parallel to the N-D vector. The -ND3 groups undergo 3-fold hops about the C-N bond axis and have typical quadrupole couplings; the D2O undergoes 2-fold hops, leading to a tensor with a large asymmetry parameter. With appropriate corrections for vibrational averaging, density functional cluster calculations give an excellent fit to the imidazolium tensor magnitudes and orientations.