Modeling Small Structural and Environmental Differences in Solids with 15 N NMR Chemical Shift Tensors

Predicting 35-Cl electric field gradient tensors in crystalline solids using cluster and fragment-corrected planewave density functional theory

Planewave-corrected methods have proven effective for accurately modeling nuclear magnetic resonance (NMR) parameters in crystalline systems. Recent work extended the application of planewave-corrected calculations beyond the second row, predicting EFG tensor parameters for 35Cl using a simple molecular correction to projector augmented-wave (PAW) density functional theory (DFT). Here we extend this work using fragment and cluster-based calculations coupled with polarizable continuum (PCM) methods to improve further the accuracy of planewave-corrected 35Cl EFG tensor calculations. Benchmark data from a test set comprised of 105 individual 35Cl EFG tensor principal components for chlorine-containing molecular crystals and crystalline chloride salts shows fragment-corrected planewave calculations using the PBE0 hybrid density functional improve the accuracy of predicted EFG tensor components by 30 % relative to traditional planewave calculations. We compare the influence of different geometry optimization methods and density functionals on the accuracy of predicted 35Cl EFG tensor parameters. Four cases of spectral assignment are presented to demonstrate the utility of improving the accuracy of predicted 35Cl EFG tensor parameters.

Predicting 51V nuclear magnetic resonance observables in molecular crystals

Solid-state nuclear magnetic resonance (NMR) spectroscopy and quantum chemical density functional theory (DFT) calculations are widely used to characterize vanadium centers in biological and pharmaceutically relevant compounds. Several techniques have been recently developed to improve the accuracy of predicted NMR parameters obtained from DFT. Fragment-based and planewave-corrected methods employing hybrid density functionals are particularly effective tools for solid-state applications. A recent benchmark study involving molecular crystal compounds found that fragment-based NMR calculations using hybrid density functionals improve the accuracy of predicted 51V chemical shieldings by 20% relative to traditional planewave methods. This work extends the previous study, including a careful analysis of 51V chemical shift anisotropy, electric field gradient calculations, and a more extensive test set. The accuracy of planewave-corrected techniques and recently developed fragment-based methods using electrostatic embedding based on the polarized continuum model (PCM) are found to be highly competitive with previous methods. Planewave-corrected methods achieve a 34% improvement in the errors of predicted 51V chemical shieldings relative to planewave. Additionally, planewave-corrected and fragment-based calculations were performed using PCM embedding, improving the accuracy of predicted 51V chemical shielding (CS) tensor principal values by 30% andC q values by 15% relative to traditional planewave methods. The performance of these methods is further examined using a redox-active oxovandium complex and a common 51V NMR reference compound.

Acceleration of indirect detection 195Pt solid-state NMR experiments by sideband selective excitation or alternative indirect sampling schemes

The measurement of the of chemical shift (CS) tensors via solid-state NMR (ssNMR) spectroscopy has proven to be a powerful probe of structure for organic molecules, biomolecules, and inorganic materials. However, when measuring the NMR spectra of heavy spin-1/2 isotopes the chemical shift anisotropy (CSA) is commonly on the order of thousands of parts per million, which makes acquisition of NMR spectra difficult due to the low NMR sensitivity imposed by the breadth of the signals and challenges in uniformly exciting the NMR spectrum. We have recently shown that complete ¹⁹⁵Pt NMR spectra could be rapidly measured by using ¹⁹⁵Pt saturation or excitation selective long pulses (SLP) with multiple rotor-cycle durations and RF fields less than 50 kHz into ¹H{¹⁹⁵Pt} or ¹H-³¹P{¹⁹⁵Pt} PE S-RESPDOR, TONE D-HMQC-4, J-resolved, and J-HMQC pulse sequences. The SLP only provide signal or dephasing when they are applied on resonance with a spinning sideband. The magic angle spinning ¹⁹⁵Pt NMR spectrum is reconstructed in the sideband selective NMR experiments by acquiring 1D NMR spectra at variable ¹⁹⁵Pt pulse offsets. In this work, we present a detailed investigation of the specific pulse conditions required for the ideal performance of sideband selective experiments. Sideband selective experiments are shown to be able to accurately reproduce MAS NMR spectra with minimal distortions of relative sideband intensities. It is also demonstrated that a ¹⁹⁵Pt NMR spectrum indirectly detected with HMQC can be rapidly obtained by acquiring a single rotor cycle of indirect dimension evolution points. We dub this method One Rotor Cycle of Acquisition (ORCA) HMQC. Sideband selective experiments and ORCA HMQC experiments are shown to provide a one order of magnitude improvement in experiment times as compared to conventional wideline HMQC experiments.

Benchmark accuracy of predicted NMR observables for quadrupolar 14 N and 17 O nuclei in molecular crystals

Nuclear quadrupole resonances for 14 N and ¹⁷ O nuclei are exquisitely sensitive to interactions with surrounding atoms. As a result, nitrogen and oxygen solid-state nuclear magnetic resonance (ssNMR) provides a powerful tool for investigating structure and dynamics in complex systems. First-principles calculations are increasingly used to facilitate spectral assignment and to evaluate and adjust crystal structures. Recent work combining the strengths of planewave density functional theory (DFT) calculations with a single molecule correction obtained using a higher level of theory has proven successful in improving the accuracy of predicted chemical shielding (CS) tensors and ¹⁷ O quadrupolar coupling constants ( C q ). Here we extend this work by examining the accuracy of predicted ¹⁴ N and ¹⁷ O electric field gradient (EFG) tensor components obtained using alternative planewave-corrections involving cluster and two-body fragment-based calculations. We benchmark the accuracy of CS and EFG tensor predictions for both nitrogen and oxygen using planewave, two-body fragment, and enhanced planewave-corrected techniques. Combining planewave and two-body fragment calculations reduces the error in predicted ¹⁷ O C q values by 35% relative to traditional planewave calculations. These enhanced planewave-correction methods improve the accuracy of ¹⁷ O and ¹⁴ N EFG tensor components by 15% relative to planewave DFT but yield minimal improvement relative to a simple molecular correction. However, in structural environments involving either high symmetry or strong intermolecular interactions, enhanced planewave-corrected methods provide a distinct advantage.

Improving the accuracy of GIPAW chemical shielding calculations with cluster and fragment corrections

Ab initio methods for predicting NMR parameters in the solid state are an essential tool for assigning experimental spectra and play an increasingly important role in structural characterizations. Recently, a molecular correction (MC) technique has been developed which combines the strengths of plane-wave methods (GIPAW) with single molecule calculations employing Gaussian basis sets. The GIPAW + MC method relies on a periodic calculation performed at a lower level of theory to model the crystalline environment. The GIPAW result is then corrected using a single molecule calculation performed at a higher level of theory. The success of the GIPAW + MC method in predicting a range of NMR parameters is a result of the highly local character of the tensors underlying the NMR observable. However, in applications involving strong intermolecular interactions we find that expanding the region treated at the higher level of theory more accurately captures local many-body contributions to the N15 NMR chemical shielding (CS) tensor. We propose alternative corrections to GIPAW which capture interactions between adjacent molecules at a higher level of theory using either fragment or cluster-based calculations. Benchmark calculations performed on N15 and C13 data sets show that these advanced GIPAW-corrected calculations improve the accuracy of chemical shielding tensor predictions relative to existing methods. Specifically, cluster-based N15 corrections show a 24% and 17% reduction in RMS error relative to GIPAW and GIPAW + MC calculations, respectively. Comparing the benchmark data sets using multiple computational models demonstrates that N15 CS tensor calculations are significantly more sensitive to intermolecular interactions relative to C13. However, fragment and cluster-based corrections that include direct hydrogen bond partners are sufficient for optimizing the accuracy of GIPAW-corrected methods. Finally, GIPAW-corrected methods are applied to the particularly challenging NMR spectral assignment of guanosine dihydrate which contains two guanosine molecules in the asymmetric unit.

Fast and Accurate Electric Field Gradient Calculations in Molecular Solids With Density Functional Theory

Modern approaches for calculating electric field gradient (EFF) tensors in molecular solids rely upon plane-wave calculations employing periodic boundary conditions (PBC). In practice, models employing PBCs are limited to generalized gradient approximation (GGA) density functionals. Hybrid density functionals applied in the context of gauge-including atomic orbital (GIAO) calculations have been shown to substantially improve the accuracy of predicted NMR parameters. Here we propose an efficient method that effectively combines the benefits of both periodic calculations and single-molecule techniques for predicting electric field gradient tensors in molecular solids. Periodic calculations using plane-wave basis sets were used to model the crystalline environment. We then introduce a molecular correction to the periodic result obtained from a single-molecule calculation performed with a hybrid density functional. Single-molecule calculations performed using hybrid density functionals were found to significantly improve the agreement of predicted 17O quadrupolar coupling constants (C q ) with experiment. We demonstrate a 31% reduction in the RMS error for the predicted 17O C q values relative to standard plane-wave methods using a carefully constructed test set comprised of 22 oxygen-containing molecular crystals. We show comparable improvements in accuracy using five different hybrid density functionals and find predicted C q values to be relatively insensitive to the choice of basis set used in the single molecule calculation. Finally, the utility of high-accuracy 17O C q predictions is demonstrated by examining the disordered 4-Nitrobenzaldehyde crystal structure.