skip to main content


Title: Calculating Nuclear Magnetic Resonance Chemical Shifts from Density Functional Theory: A Primer
Density functional theory (DFT) prediction of nuclear magnetic resonance (NMR) chemical shifts complements NMR experiments. Predicting chemical shifts accurately with DFT requires many different modeling decisions. Intended for novice modelers and nonexperts, this article discusses the considerations one should take in selecting a density functional, van der Waals dispersion correction, and basis set. It examines different strategies for handling systems in complex environments such as liquids, biomolecules, and crystals. Strategies include the use of cluster models, electrostatic embedding, continuum representations, periodic boundary conditions, and fragment‐based approaches. Finally, approaches for referencing the predicted absolute chemical shieldings for comparison against experimental chemical shifts are discussed.  more » « less
Award ID(s):
1665212
NSF-PAR ID:
10165326
Author(s) / Creator(s):
Date Published:
Journal Name:
eMagRes
Volume:
8
Issue:
3
ISSN:
2055-6101
Page Range / eLocation ID:
215-226
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Density functional theory (DFT) benchmark studies of 1H and 13C NMR chemical shifts often yield differing conclusions, likely due to non-optimal test molecules and non-standardized data acquisition. To address this issue, we carefully selected and measured 1H and 13C NMR chemical shifts for 50 structurally diverse small organic molecules containing atoms from only the first two rows of the periodic table. Our NMR dataset, DELTA50, was used to calculate linear scaling factors and to evaluate the accuracy of 73 density functionals, 40 basis sets, 3 solvent models, and 3 gauge-referencing schemes. The best performing DFT methodologies for 1H and 13C NMR chemical shift predictions were WP04/6-311++G(2d,p) and ωB97X-D/def2-SVP, respectively, when combined with the polarizable continuum solvent model (PCM) and gauge-independent atomic orbital (GIAO) method. Geometries should be optimized at the B3LYP-D3/6-311G(d,p) level including the PCM solvent model for the best accuracy. Predictions of 20 organic compounds and natural products from a separate probe set had root-mean-square deviations (RMSD) of 0.07 to 0.19 for 1H and 0.5 to 2.9 for 13C. Maximum deviations were less than 0.5 and 6.5 ppm for 1H and 13C, respectively.

     
    more » « less
  2. Abstract

    Cannabicitran is an important cannabinoid natural product produced byCannabis sativaand is often found at surprisingly high levels (up to ~10%) in “purified” commercial cannabidiol (CBD) extract preparations. Despite the prevalence of this molecule in CBD oil and other cannabinoid‐related products, and the rapidly expanding interest in cannabinoids for treatment of a wide range of physiological conditions, only unassigned1H NMR data and partial unambiguous13C assignments have been published. Herein, we report the complete1H and13C NMR assignments of cannabicitran and comparatively evaluate the performance of several density functional theory (DFT) methods with varying levels of theory for the calculation of NMR chemical shifts.

     
    more » « less
  3. Nuclear magnetic resonance (NMR) is one of the primary techniques used to elucidate the chemical structure, bonding, stereochemistry, and conformation of organic compounds. The distinct chemical shifts in an NMR spectrum depend upon each atom's local chemical environment and are influenced by both through-bond and through-space interactions with other atoms and functional groups. The in silico prediction of NMR chemical shifts using quantum mechanical (QM) calculations is now commonplace in aiding organic structural assignment since spectra can be computed for several candidate structures and then compared with experimental values to find the best possible match. However, the computational demands of calculating multiple structural- and stereo-isomers, each of which may typically exist as an ensemble of rapidly-interconverting conformations, are expensive. Additionally, the QM predictions themselves may lack sufficient accuracy to identify a correct structure. In this work, we address both of these shortcomings by developing a rapid machine learning (ML) protocol to predict 1 H and 13 C chemical shifts through an efficient graph neural network (GNN) using 3D structures as input. Transfer learning with experimental data is used to improve the final prediction accuracy of a model trained using QM calculations. When tested on the CHESHIRE dataset, the proposed model predicts observed 13 C chemical shifts with comparable accuracy to the best-performing DFT functionals (1.5 ppm) in around 1/6000 of the CPU time. An automated prediction webserver and graphical interface are accessible online at http://nova.chem.colostate.edu/cascade/. We further demonstrate the model in three applications: first, we use the model to decide the correct organic structure from candidates through experimental spectra, including complex stereoisomers; second, we automatically detect and revise incorrect chemical shift assignments in a popular NMR database, the NMRShiftDB; and third, we use NMR chemical shifts as descriptors for determination of the sites of electrophilic aromatic substitution. 
    more » « less
  4. null (Ed.)
    Rechargeable aluminum–graphite batteries are an emerging energy storage technology with great promise: they exhibit high rate performance, cyclability, and a discharge potential of approximately 2 V, while both electrodes are globally abundant, low cost, and inherently safe. The batteries use chloroaluminate-containing electrolytes and store charge in the graphite electrodes when molecular AlCl 4 − anions electrochemically intercalate within them. However, much remains to be understood regarding the ion intercalation mechanism, in part due to challenges associated with characterizing the chloroaluminate anions themselves. Here, we use solid-state 27 Al nuclear magnetic resonance (NMR) spectroscopy to probe the molecular-level electronic and magnetic environments of intercalated chloroaluminate anions at different states-of-charge. The results reveal broad 27 Al NMR signals associated with intercalated AlCl 4 − anions, reflecting high extents of local disorder. The intercalated anions experience a diversity of local environments, many of which are far from the ideal crystalline-like structures often depicted in graphite staging models. Density functional theory (DFT) calculations of the total 27 Al isotropic shifts enable the contributions of chemical shift, ring-current effects, and electric quadrupolar interactions to be disentangled quantitatively. In combination, the solid-state NMR and DFT results reveal the molecular geometries and environments of intercalated AlCl 4 − anions and capture the significant disorder present in intercalated graphite battery electrodes. 
    more » « less
  5. A fast, straightforward method for computing NMR chemical shieldings of crystalline solids is proposed. The method combines the advantages of both conventional approaches: periodic calculations using plane-wave basis sets and molecular computational approaches. The periodic calculations capture the periodic nature of crystalline solids, but the computational level of the electronic structure calculation is limited to general-gradient-approximation (GGA) density functionals. It is demonstrated that a correction to the GGA result calculated on an isolated molecule at a higher level of theory significantly improves the correlations between experimental and calculated chemical shifts while adding almost no additional computational cost. Corrections calculated with a hybrid density functional improved the accuracy of 13C, 15N and 17O chemical shift predictions significantly and allowed identifying errors in previously published experimental data. Applications of the approach to crystalline isocytosine, methacrylamide, and testosterone are presented. 
    more » « less