Gauge Including Projector Augmented Wave Research Articles

Polymorph Identification for Flexible Molecules: Linear Regression Analysis of Experimental and Calculated Solution- and Solid-State NMR Data.

The Δδ regression approach of Blade et al. [ J. Phys. Chem. A 2020, 124(43), 8959-8977] for accurately discriminating between solid forms using a combination of experimental solution- and solid-state NMR data with density functional theory (DFT) calculation is here extended to molecules with multiple conformational degrees of freedom, using furosemide polymorphs as an exemplar. As before, the differences in measured 1H and 13C chemical shifts between solution-state NMR and solid-state magic-angle spinning (MAS) NMR (Δδexperimental) are compared to those determined by gauge-including projector augmented wave (GIPAW) calculations (Δδcalculated) by regression analysis and a t-test, allowing the correct furosemide polymorph to be precisely identified. Monte Carlo random sampling is used to calculate solution-state NMR chemical shifts, reducing computation times by avoiding the need to systematically sample the multidimensional conformational landscape that furosemide occupies in solution. The solvent conditions should be chosen to match the molecule's charge state between the solution and solid states. The Δδ regression approach indicates whether or not correlations between Δδexperimental and Δδcalculated are statistically significant; the approach is differently sensitive to the popular root mean squared error (RMSE) method, being shown to exhibit a much greater dynamic range. An alternative method for estimating solution-state NMR chemical shifts by approximating the measured solution-state dynamic 3D behavior with an ensemble of 54 furosemide crystal structures (polymorphs and cocrystals) from the Cambridge Structural Database (CSD) was also successful in this case, suggesting new avenues for this method that may overcome its current dependency on the prior determination of solution dynamic 3D structures.

1H and 13C chemical shift-structure effects in anhydrous β-caffeine and four caffeine-diacid cocrystals probed by solid-state NMR experiments and DFT calculations.

By using density functional theory (DFT) calculations, we refined the H atom positions in the structures of β-caffeine (C), α-oxalic acid (OA; (COOH)2), α-(COOH)2·2H2O, β-malonic acid (MA), β-glutaric acid (GA), and I-maleic acid (ME), along with their corresponding cocrystals of 2 : 1 (2C-OA, 2C-MA) or 1 : 1 (C-GA, C-ME) stoichiometry. The corresponding 13C/1H chemical shifts obtained by gauge including projector augmented wave (GIPAW) calculations agreed overall very well with results from magic-angle-spinning (MAS) nuclear magnetic resonance (NMR) spectroscopy experiments. Chemical-shift/structure trends of the precursors and cocrystals were examined, where good linear correlations resulted for all COO1H sites against the H⋯O and/or H⋯N H-bond distance, whereas a general correlation was neither found for the aliphatic/caffeine-stemming 1H sites nor any 13C chemical shift against either the intermolecular hydrogen- or tetrel-bond distance, except for the 13COOH sites of the 2C-OA, 2C-MA, and C-GA cocrystals, which are involved in a strong COOH⋯N bond with caffeine that is responsible for the main supramolecular stabilization of the cocrystal. We provide the first complete 13C NMR spectral assignment of the structurally disordered anhydrous β-caffeine polymorph. The results are discussed in relation to previous literature on the disordered α-caffeine polymorph and the ordered hydrated counterpart, along with recommendations for NMR experimentation that will secure sufficient 13C signal-resolution for reliable resonance/site assignments.

Do Models beyond Hybrid Density Functionals Increase the Agreement with Experiment for Predicted NMR Chemical Shifts or Electric Field Gradient Tensors in Organic Solids?

Ab initio predictions of chemical shifts and electric field gradient (EFG) tensor components are frequently used to help interpret solid-state nuclear magnetic resonance (NMR) experiments. Typically, these predictions employ density functional theory (DFT) with generalized gradient approximation (GGA) functionals, though hybrid functionals have been shown to improve accuracy relative to experiment. Here, the performance of a dozen models beyond the GGA approximation are examined for the prediction of solid-state NMR observables, including meta-GGA, hybrid, and double-hybrid density functionals and second-order Møller-Plesset perturbation theory (MP2). These models are tested on organic molecular crystal data sets containing 169 experimental 13C and 15N chemical shifts and 114 17O and 14N EFG tensor components. To make these calculations affordable, gauge-including projector augmented wave (GIPAW) Perdew-Burke-Ernzerhof (PBE) calculations with periodic boundary conditions are combined with a local intramolecular correction computed at the higher level of theory. Within the context of typical NMR property calculations performed on a static, DFT-optimized crystal structure, the benchmarking finds that the double-hybrid DFT functionals produce errors versus experiment that are no smaller than those of hybrid functionals in the best cases, and they can be larger. MP2 errors versus experiment are even bigger. Overall, no practical advantages are found for using any of the tested double-hybrid functionals or MP2 to predict experimental solid-state NMR chemical shifts and EFG tensor components for routine organic crystals, especially given the higher computational cost of those methods. This finding likely reflects error cancellation benefiting the hybrid functionals. Improving the accuracy of the predicted chemical shifts and EFG tensors relative to experiment would probably require more robust treatments of the crystal structures, their dynamics, and other factors.

Discovering the Solid-State Secrets of Lorlatinib by NMR Crystallography: To Hydrogen Bond or not to Hydrogen Bond

Lorlatinib is an active pharmaceutical ingredient (API) used in the treatment of lung cancer. Here, an NMR crystallography analysis is presented whereby the single-crystal X-ray diffraction structure (CSD: 2205098) determination is complemented by multinuclear (1H, 13C, 14/15N, 19F) magic-angle spinning (MAS) solid-state NMR and gauge-including projector augmented wave (GIPAW) calculation of NMR chemical shifts. Lorlatinib crystallises in the P21 space group, with two distinct molecules in the asymmetric unit cell, Z′ = 2. Three of the four NH2 hydrogen atoms form intermolecular hydrogen bonds, N30-H…N15 between the two distinct molecules and N30-H…O2 between two equivalent molecules. This is reflected in one of the NH21H chemical shifts being significantly lower, 4.0 ppm compared to 7.0 ppm. Two-dimensional 1H–13C, 14N-1H and 1H (double-quantum, DQ)-1H (single-quantum, SQ) MAS NMR spectra are presented. The 1H resonances are assigned and specific HH proximities corresponding to the observed DQ peaks are identified. The resolution enhancement at a 1H Larmor frequency of 1 GHz as compared to 500 or 600 MHz is demonstrated.

New Approach To Understanding the Experimental 133Cs NMR Chemical Shift of Clay Minerals via Machine Learning and DFT-GIPAW Calculations.

Structural determination of adsorbed atoms on layered structures such as clay minerals is a complex subject. Radioactive cesium (Cs) is an important element for environmental conservation, so it is vital to understand its adsorption structure on clay. The nuclear magnetic resonance (NMR) parameters of 133Cs, which can be determined from solid-state NMR experiments, are sensitive to the local neighboring structures of adsorbed Cs. However, determining the Cs positions from NMR data alone is difficult. This paper describes an approach for identifying the expected atomic positions on clay minerals by combining machine learning (ML) with experimentally observed chemical shifts. A linear ridge regression model for ML is constructed from the smooth overlap of atomic position descriptor and gauge-including projector augmented wave (GIPAW) ab initio data. The constructed ML model predicts the GIPAW data to within a 3 ppm root-mean-squared error. At this stage, the 133Cs chemical shifts can be instantaneously calculated from the Cs positions on any clay layers using ML. The inverse analysis, which derives the atomic positions from experimentally observed chemical shifts, is developed from the ML model. The input data for the inverse analysis are the layer structure and the experimentally observed chemical shifts. The Cs positions for the targeted chemical shifts are then output. Inverse analysis is applied to montmorillonite, and the resultant Cs positions are found to be consistent with previous results (Ohkubo, T.; et al. J. Phys. Chem. A 2018, 122, 9326-9337). The Cs positions on saponite clay are also clarified from experimentally observed chemical shifts and inverse analysis.

Calculated electron paramagnetic resonance g tensor and hyperfine parameters for zinc vacancy and N related defects in ZnO

Various defects in ZnO, focused on substitutional N$_O$ and N$_2$ in various sites, O-site, interstitial and Zn-site are studied using first-principles calculations with the goal of understanding the electron paramagnetic resonance (EPR) center reported for N$_2$ in ZnO and substitutional N on the O-site. The $g$ tensors are calculated using the gauge including projector augmented wave (GIPAW) method and compared with experiments. The $g$-tensor of the free N$_2^+$ and N$_2^-$ radicals and their various contributions within the GIPAW theory are analyzed first to provide a baseline reference for the accuracy of the method and for understanding the N$_2$ behavior in ZnO. Previous controversies on the site location of N$_2$ in ZnO for this EPR center and on the shallow or deep nature and donor or acceptor nature of this center are resolved. We find that the N$_2$ on the Zn site is mostly zinc-vacancy like in its spin density and $g$-tensor, while for the O-site, a model with the N$_2$ axis lying in-the basal plane and the singly occupied $\pi_g$-orbital along the {\bf c} axis provides good agreement with experiment. For the interstitial location, if the N$_2$ is not strongly interacting with the surroundings, no levels in the gap are found and hence also no possible EPR center. The calculated $g$-tensors for N$_O$ and $V_ Zn$ are also found to be in good agreement with experiment. The effects of different functionals affecting the localization of the spin density are shown to affect the $g$-tensor values.

Combining heteronuclear correlation NMR with spin-diffusion to detect relayed Cl–H–H and N–H–H proximities in molecular solids

Analysis of short-to-intermediate range intermolecular interactions offers a great way of characterizing the solid-state organization of small molecules and materials. This can be achieved by two-dimensional (2D) homo- and heteronuclear correlation NMR spectroscopy, for example, by carrying out experiments at high magnetic fields in conjunction with fast magic-angle spinning (MAS) techniques. But, detecting 2D peaks for heteronuclear dipolar coupled spin pairs separated by greater than 3 Å is not always straightforward, particularly when low-gamma quadrupolar nuclei are involved. Here, we present a 2D correlation NMR experiment that combines the advantages of heteronuclear-multiple quantum coherence (HMQC) and proton-based spin-diffusion (SD) pulse sequences using radio-frequency-driven-recouping (RFDR) to probe inter and intramolecular 1H-X (X = 14N, 35Cl) interactions. This experiment can be used to acquire 2D 1H{X}-HMQC filtered 1H–1H correlation as well as 2D 1H-X HMQC spectra. Powder forms of dopamine·HCl and l-histidine·HCl·H2O are characterized at high fields (21.1 T and 18.8 T) with fast MAS (60 kHz) using the 2D HMQC-SD-RFDR approach. Solid-state NMR results are complemented with NMR crystallography analyses using the gauge-including projector augmented wave (GIPAW) approach. For histidine·HCl·H2O, 2D peaks associated with 14N–1H–1H and 35Cl–1H–1H distances of up to 4.4 and 3.9 Å have been detected. This is further corroborated by the observation of 2D peaks corresponding to 14N–1H–1H and 35Cl–1H–1H distances of up to 4.2 and 3.7 Å in dopamine·HCl, indicating the suitability of the HMQC-SD-RFDR experiments for detecting medium-range proximities in molecular solids.

To what extent do bond length and angle govern the 13C and 1H NMR response to weak CH⋯O hydrogen bonds? A case study of caffeine and theophylline cocrystals

Weak hydrogen bonds are important structure-directing elements in supramolecular chemistry and biochemistry. We consider here weak CH⋯O hydrogen bonds in a series of cocrystals of theophylline and caffeine and assess to what extent the CH⋯O distance and angle govern the observed 13C and 1H isotropic chemical shifts. Gauge-including projector-augmented wave density functional theory (GIPAW DFT) calculations consistently predict a decrease in the 13C and 1H magnetic shielding constants upon hydrogen bond formation on the order of 2–5 ppm (13C) and 1–2 ppm (1H). These trends are reproduced using the machine-learning approach implemented in ShiftML. Experimental 13C and 1H chemical shifts obtained for powdered samples using one-dimensional NMR spectroscopy as well as heteronuclear correlation (HETCOR) spectroscopy correlate well with the GIPAW DFT results. However, the experimental 13C NMR response only correlates moderately well with the hydrogen bond length and angle, while the experimental 1H chemical shifts only show very weak correlations to these local structural elements. DFT computations on isolated imidazole-formaldehyde models show that the 13C and 1H chemical shifts generally decrease with the C⋯O distance but show no clear dependence on the CH⋯O angle. These results demonstrate that the 13C and 1H response to weak CH⋯O hydrogen bonding is influenced significantly by additional weak contacts within cocrystal heterodimeric units.

Comparing GIPAW with numerically exact chemical shieldings: The role of two-center contributions to the induced current.

In this study, we benchmark density functional theory gauge-including projector-augmented-wave (GIPAW) chemical shieldings against molecular shieldings for which basis set completeness has been achieved [Jensen et al., Phys. Chem. Chem. Phys. 18, 21145 (2016)]. We demonstrate the importance of two-center corrections for GIPAW hydrogen shieldings. For the other nuclei studied, standard GIPAW is sufficiently accurate. We find that GIPAW can be pushed to closely approach the basis set limit. The only source of small inaccuracies lies in the contribution to the shielding that is caused by surface currents, which we estimate comparing GIPAW susceptibilities to converged molecular magnetizabilities.

Refined Structures of O-Phospho-l-serine and Its Calcium Salt by New Multinuclear Solid-State NMR Crystallography Methods.

O-phospho-l-serine (Pser) and its Ca salt, Ca[O-phospho-l-serine]·H2O (CaPser), play important roles for bone mineralization and were recently also proposed to account for the markedly improved bone-adhesive properties of Pser-doped calcium phosphate-based cements for biomedical implants. However, the hitherto few proposed structural models of Pser and CaPser were obtained by X-ray diffraction, thereby leaving the proton positions poorly defined. Herein, we refine the Pser and CaPser structures by density functional theory (DFT) calculations and contrast them with direct interatomic-distance constraints from two-dimensional (2D) nuclear magnetic resonance (NMR) correlation experimentation at fast magic-angle spinning (MAS), encompassing double-quantum–single-quantum (2Q–1Q) 1H NMR along with heteronuclear 13C{1H} and 31P{1H} correlation NMR experiments. The Pser and CaPser structures before and after refinements by DFT were validated against sets of NMR-derived effective 1H–1H, 1H–31P, and 1H–13C distances, which confirmed the improved accuracy of the refined structures. Each distance set was derived from one sole 2D NMR experiment applied to a powder without isotopic enrichment. The distances were extracted without invoking numerical spin-dynamics simulations or approximate phenomenological models. We highlight the advantages and limitations of the new distance-extraction procedure. Isotropic 1H, 13C, and 31P chemical shifts obtained by DFT calculations using the gauge including projector augmented wave (GIPAW) method agreed very well with the experimental results. We discuss the isotropic and anisotropic 13C and 31P chemical-shift parameters in relation to the previous literature, where most data on CaPser are reported herein for the first time.

35 Cl-1 H Heteronuclear correlation magic-angle spinning nuclear magnetic resonance experiments for probing pharmaceutical salts.

Heteronuclear multiple-quantum coherence (HMQC) pulse sequences for establishing heteronuclear correlation in solid-state nuclear magnetic resonance (NMR) between 35 Cl and 1 H nuclei in chloride salts under fast (60 kHz) magic-angle spinning (MAS) and at high magnetic field (a 1 H Larmor frequency of 850 MHz) are investigated. Specifically, recoupling of the 35 Cl-1 H dipolar interaction using rotary resonance recoupling with phase inversion every rotor period or the symmetry-based SR42 1 pulse sequences are compared. In our implementation of the population transfer (PT) dipolar (D) HMQC experiment, the satellite transitions of the 35 Cl nuclei are saturated with an off-resonance WURST sweep, at a low nutation frequency, over the second spinning sideband, whereby the WURST pulse must be of the same duration as the recoupling time. Numerical simulations of the 35 Cl-1 H MAS D-HMQC experiment performed separately for each crystallite orientation in a powder provide insight into the orientation dependence of changes in the second-order quadrupolar-broadened 35 Cl MAS NMR lineshape under the application of dipolar recoupling. Two-dimensional 35 Cl-1 H PT-D-HMQC MAS NMR spectra are presented for the amino acids glycine·HCl and l-tyrosine·HCl and the pharmaceuticals cimetidine·HCl, amitriptyline·HCl and lidocaine·HCl·H2 O. Experimentally observed 35 Cl lineshapes are compared with those simulated for 35 Cl chemical shift and quadrupolar parameters as calculated using the gauge-including projector-augmented wave (GIPAW) method: the calculated quadrupolar product (PQ ) values exceed those measured experimentally by a factor of between 1.3 and 1.9.

Synergy of Solid-State NMR, Single-Crystal X-ray Diffraction, and Crystal Structure Prediction Methods: A Case Study of Teriflunomide (TFM).

In this work, for the first time, we present the X-ray diffraction crystal structure and spectral properties of a new, room-temperature polymorph of teriflunomide (TFM), CSD code 1969989. As revealed by DSC, the low-temperature TFM polymorph recently reported by Gunnam et al. undergoes a reversible thermal transition at −40 °C. This reversible process is related to a change in Z’ value, from 2 to 1, as observed by variable-temperature 1H–13C cross-polarization (CP) magic-angle spinning (MAS) solid-state NMR, while the crystallographic system is preserved (triclinic). Two-dimensional 13C–1H and 1H–1H double-quantum MAS NMR spectra are consistent with the new room-temperature structure, including comparison with GIPAW (gauge-including projector augmented waves) calculated NMR chemical shifts. A crystal structure prediction procedure found both experimental teriflunomide polymorphs in the energetic global minimum region. Differences between the polymorphs are seen for the torsional angle describing the orientation of the phenyl ring relative to the planarity of the TFM molecule. In the low-temperature structure, there are two torsion angles of 4.5 and 31.9° for the two Z’ = 2 molecules, while in the room-temperature structure, there is disorder that is modeled with ∼50% occupancy between torsion angles of −7.8 and 28.6°. These observations are consistent with a broad energy minimum as revealed by DFT calculations. PISEMA solid-state NMR experiments show a reduction in the C–H dipolar coupling in comparison to the static limit for the aromatic CH moieties of 75% and 51% at 20 and 40 °C, respectively, that is indicative of ring flips at the higher temperature. Our study shows the power of combining experiments, namely DSC, X-ray diffraction, and MAS NMR, with DFT calculations and CSP to probe and understand the solid-state landscape, and in particular the role of dynamics, for pharmaceutical molecules.

New Pharmaceutical Salts of Trazodone.

New pharmaceutically acceptable salts of trazodone (trazodone hydrogen bromide and trazodone 1-hydroxy-2-naphthonic acid) for the treatment of central nervous system disorders are synthesized and described. Although trazodone salts are poorly crystalline, single-crystal X-ray diffraction data for trazodone 1-hydroxy-2-naphthonic acid were collected and analyzed as well as compared to the previously described crystal structure of commercially available trazodone hydrochloride. The powder samples of all new salts were characterized by Fourier transform infrared spectroscopy, X-ray diffraction and 13C solid-state nuclear magnetic resonance spectroscopy. Spectroscopic studies were supported by gauge including projector augmented wave (GIPAW) calculations of carbon chemical shielding constants. The main goal of our research was to find salts with better physicochemical properties and to make an attempt to associate them with both the anion structure and the most prominent interactions exhibited by the protonated trazodone cation. The dissolution profiles of trazodone from tablets prepared from various salts with lactose monohydrate were investigated. The studies revealed that salts with simple anions show a fast release of the drug while the presence of more complex anion, more strongly interacting with the cation, effects a slow-release profile of the active substance and can be used for the preparation of the tables with a delay or prolonged mode of action.

Influence of Polymorphs and Local Defect Structures on NMR Parameters of Graphite Fluorides

The role of local molecular structure on calculated 13C and 19F NMR chemical shifts for graphite fluoride materials was explored by using gauge-including projector augmented wave (GIPAW) computatio...

A structural analysis of 2,5-diaryl-4H-2,4-dihydro-3H-1,2,4-triazol-3-ones: NMR in the solid state, X-ray crystallography, and GIPAW calculations.

The 1 H, 13 C, 15 N, and 19 F nuclear magnetic resonance (NMR) spectra of 11 2,5-diaryl-2,4-dihydro-3H-1,2,4-triazol-3-ones have been acquired in DMSO-d6 solution and the 13 C, 15 N, and 19 F NMR spectra have also been acquired in the solid state (solid-state nuclear magnetic resonance [SSNMR] and magic angle spinning [MAS]). The X-ray structures of Compounds 3, 5, and 6 have been determined by X-ray diffraction. Theoretical calculations at the gauge-independent atomic orbital (GIAO)/B3LYP/6-311++G(d,p) level have provided a set of 321 chemical shifts that were compared with 310 experimental values in DMSO-d6 . To obtain good agreements, some effects need to be included. The SSNMR chemical shifts have been compared with gauge-including projector-augmented wave (GIPAW) calculations and with the heavy atom-light atom (HALA) effects.

Investigating discrepancies between experimental solid-state NMR and GIPAW calculation: N[dbnd]C–N 13C and OH⋯O 1H chemical shifts in pyridinium fumarates and their cocrystals

An NMR crystallography analysis is presented for four solid-state structures of pyridine fumarates and their cocrystals, using crystal structures deposited in the Cambridge Crystallographic Data Centre, CCDC. Experimental one-dimensional one-pulse 1H and 13C cross-polarisation (CP) magic-angle spinning (MAS) nuclear magnetic resonance (NMR) and two-dimensional 14N–1H heteronuclear multiple-quantum coherence MAS NMR spectra are compared with gauge-including projector augmented wave (GIPAW) calculations of the 1H and 13C chemical shifts and the 14N shifts that additionally depend on the quadrupolar interaction. Considering the high ppm (>10 ​ppm) 1H resonances, while there is good agreement (within 0.4 ​ppm) between experiment and GIPAW calculation for the hydrogen-bonded NH moieties, the hydrogen-bonded fumaric acid OH resonances are 1.2–1.9 ​ppm higher in GIPAW calculation as compared to experiment. For the cocrystals of a salt and a salt formed by 2-amino-5-methylpyridinium and 2-amino-6-methylpyridinium ions, a large discrepancy of 4.2 and 5.9 ​ppm between experiment and GIPAW calculation is observed for the quaternary ring carbon 13C resonance that is directly bonded to two nitrogens (in the ring and in the amino group). By comparison, there is excellent agreement (within 0.2 ​ppm) for the quaternary ring carbon 13C resonance directly bonded to the ring nitrogen for the salt and cocrystal of a salt formed by 2,6-lutidinium and 2,5-lutidinium, respectively.

5-amino-2-methylpyridinium hydrogen fumarate: An XRD and NMR crystallography analysis.

Single-crystal X-ray diffraction structures of the 5-amino-2-methylpyridinium hydrogen fumarate salt have been solved at 150 and 300 K (CCDC 1952142 and 1952143). A base-acid-base-acid ring is formed through pyridinium-carboxylate and amine-carboxylate hydrogen bonds that hold together chains formed from hydrogen-bonded hydrogen fumarate ions. 1 H and 13 C chemical shifts as well as 14 N shifts that additionally depend on the quadrupolar interaction are determined by experimental magic angle spinning (MAS) solid-state nuclear magnetic resonance (NMR) and gauge-including projector-augmented wave (GIPAW) calculation. Two-dimensional homonuclear 1 H-1 H double-quantum (DQ) MAS and heteronuclear 1 H-13 C and 14 N-1 H spectra are presented. Only small differences of up to 0.1 and 0.6 ppm for 1 H and 13 C are observed between GIPAW calculations starting with the two structures solved at 150 and 300 K (after geometry optimisation of atomic positions, but not unit cell parameters). A comparison of GIPAW-calculated 1 H chemical shifts for isolated molecules and the full crystal structures is indicative of hydrogen bonding strength.

GIAO versus GIPAW: Comparison of Methods To Calculate 11B NMR Shifts of Icosahedral Closo-Heteroboranes toward Boron-Rich Borides.

In this work, we perform first-principle density functional theory calculations with the Perdew-Burke-Ernzerhof (PBE) exchange correlation functional to compare the results of the gauge-including atomic orbital (GIAO) method with the gauge-including projector-augmented wave (GIPAW) approach for isotropic 11B nuclear magnetic resonance shifts. GIPAW had been used successfully for the theoretical calculation of nuclear magnetic parameters of 11B species in strong ionic solid-phase compounds such as borates but had been applied very rarely to structures where boron is mainly involved in complex covalent bonding situations, for example, in icosahedra of boron-rich borides. Thus, we investigate the accuracy of both well-known methods and reliability of the effective treatment of core electrons on a test set containing 16 experimentally known closo-(hetero)dodecaboranes. In general, we find very good agreement between GIAO and GIPAW when compared to experimental observations. However, accidental degeneracies of the shift values are better predicted by GIPAW. The optimized molecular geometries on the PBE level agree well with gaseous electron diffraction data and lead to theoretical isotropic chemical 11B shifts with root-mean-square errors of 2.1 and 1.0 ppm depending on the used model of converting absolute shieldings to chemical shifts. The comparison with results from hybrid functionals (B3LYP, B3LYP-D2, and PBE0) shows a minor improvement in accuracy, which is in agreement with 13C shifts of sp3-hybridized species. In order to prove the reliability of the conversion parameters obtained by PBE, we report the calculated 11B shifts of 1,2-, 1,7-, and 1,12-PCB10H11 with GIAO and GIPAW to our knowledge for the first time. Additionally, Bader's analysis is carried out on the converged electron density for all boron species within the molecular test set, yielding no simple direct relation between charge and isotropic shifts.

Effects of Li Confined Motion on NMR Quadrupolar Interactions: A Combined 7 Li NMR and DFT-MD Study of LiR2 (PO4 )3 (R=Ti or Zr) Phases.

Ab initio molecular dynamics (MD) simulations and NMR GIPAW (gauge including projector augmented wave) calculations have been used to analyze the coordination and mobility of Li ions in LiTi2(PO4)3 (rhombohedral), LiZr2(PO4)3 (triclinic), and LiZr2(PO4)3 (rhombohedral) phases. Significant discrepancies are observed between static calculations of 7Li quadrupolar parameters and experimental values. The dynamical origin of this disagreement is demonstrated by incorporating in the calculations thermal vibrations and local motion of atoms with MD simulations. For LiTi2(PO4)3, the quadrupolar constant associated with Li ions grows with temperature because the local symmetry of the system decreases, whereas for the Zr phases, the quadrupolar constant decreases because thermal vibrations reduce the anisotropy of the interaction. Finally, for both Zr phases, MD yields Li distributions that compare well with disorder reported from diffraction studies.

Weak Intermolecular CH···N Hydrogen Bonding: Determination of 13CH-15N Hydrogen-Bond Mediated J Couplings by Solid-State NMR Spectroscopy and First-Principles Calculations.

Weak hydrogen bonds are increasingly hypothesized to play key roles in a wide range of chemistry from catalysis to gelation to polymer structure. Here, 15N/13C spin-echo magic-angle spinning (MAS) solid-state nuclear magnetic resonance (NMR) experiments are applied to "view" intermolecular CH···N hydrogen bonding in two selectively labeled organic compounds, 4-[15N] cyano-4'-[13C2] ethynylbiphenyl (1) and [15N3,13C6]-2,4,6-triethynyl-1,3,5-triazine (2). The synthesis of 2-15N3,13C6 is reported here for the first time via a multistep procedure, where the key element is the reaction of [15N3]-2,4,6-trichloro-1,3,5-triazine (5) with [13C2]-[(trimethylsilyl)ethynyl]zinc chloride (8) to afford its immediate precursor [15N3,13C6]-2,4,6-tris[(trimethylsilyl)ethynyl]-1,3,5-triazine (9). Experimentally determined hydrogen-bond-mediated 2hJCN couplings (4.7 ± 0.4 Hz (1) and 4.1 ± 0.3 Hz (2)) are compared with density functional theory (DFT) gauge-including projector augmented wave (GIPAW) calculations, whereby species-independent coupling values 2hKCN (29.0 × 1019 kg m-2 s-2 A-2 (1) and 27.9 × 1019 kg m-2 s-2 A-2 (2)) quantitatively demonstrate the J couplings for these "weak" CH···N hydrogen bonds to be of a similar magnitude to those for conventionally observed NH···O hydrogen-bonding interactions in uracil (2hKNO: 28.1 and 36.8 × 1019 kg m-2 s-2 A-2). Moreover, the GIPAW calculations show a clear correlation between increasing 2hJCN (and 3hJCN) coupling and reducing C(H)···N and H···N hydrogen-bonding distances, with the Fermi contact term accounting for at least 98% of the isotropic 2hJCN coupling.