Powder Crystallography by Proton Solid-State NMR Spectroscopy

Advanced solid-state NMR spectroscopy and its applications in zeolite chemistry

Solid-state NMR spectroscopy (ssNMR) can provide details about the structure, host-guest/guest-guest interactions and dynamic behavior of materials at atomic length scales. A crucial use of ssNMR is for the characterization of zeolite catalysts that are extensively employed in industrial catalytic processes. This review aims to spotlight the recent advancements in ssNMR spectroscopy and its application to zeolite chemistry. We first review the current ssNMR methods and techniques that are relevant to characterize zeolite catalysts, including advanced multinuclear and multidimensional experiments, in situ NMR techniques and hyperpolarization methods. Of these, the methodology development on half-integer quadrupolar nuclei is emphasized, which represent about two-thirds of stable NMR-active nuclei and are widely present in catalytic materials. Subsequently, we introduce the recent progress in understanding zeolite chemistry with the aid of these ssNMR methods and techniques, with a specific focus on the investigation of zeolite framework structures, zeolite crystallization mechanisms, surface active/acidic sites, host-guest/guest-guest interactions, and catalytic reaction mechanisms.

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.

Atomic-level structure determination of amorphous molecular solids by NMR

Structure determination of amorphous materials remains challenging, owing to the disorder inherent to these materials. Nuclear magnetic resonance (NMR) powder crystallography is a powerful method to determine the structure of molecular solids, but disorder leads to a high degree of overlap between measured signals, and prevents the unambiguous identification of a single modeled periodic structure as representative of the whole material. Here, we determine the atomic-level ensemble structure of the amorphous form of the drug AZD4625 by combining solid-state NMR experiments with molecular dynamics (MD) simulations and machine-learned chemical shifts. By considering the combined shifts of all 1H and 13C atomic sites in the molecule, we determine the structure of the amorphous form by identifying an ensemble of local molecular environments that are in agreement with experiment. We then extract and analyze preferred conformations and intermolecular interactions in the amorphous sample in terms of the stabilization of the amorphous form of the drug.

Nuclear spin diffusion under fast magic-angle spinning in solid-state NMR

Solid-state nuclear spin diffusion is the coherent and reversible process through which spin order is transferred via dipolar couplings. With the recent increases in magic-angle spinning (MAS) frequencies and magnetic fields becoming routinely applied in solid-state nuclear magnetic resonance, understanding how the increased 1H resolution obtained affects spin diffusion is necessary for interpretation of several common experiments. To investigate the coherent contributions to spin diffusion with fast MAS, we have developed a low-order correlation in Liouville space model based on the work of Dumez et al. (J. Chem. Phys. 33, 224501, 2010). Specifically, we introduce a new method for basis set selection, which accounts for the resonance-offset dependence at fast MAS. Furthermore, we consider the necessity of including chemical shift, both isotropic and anisotropic, in the modeling of spin diffusion. Using this model, we explore how different experimental factors change the nature of spin diffusion. Then, we show case studies to exemplify the issues that arise in using spin diffusion techniques at fast spinning. We show that the efficiency of polarization transfer via spin diffusion occurring within a deuterated and 100% back-exchanged protein sample at 60 kHz MAS is almost entirely dependent on resonance offset. We additionally identify temperature-dependent magnetization transfer in beta-aspartyl L-alanine, which could be explained by the influence of an incoherent relaxation-based nuclear Overhauser effect.

Mechanistic origins of methyl-driven Overhauser DNP

The Overhauser effect in the dynamic nuclear polarization (DNP) of non-conducting solids has drawn much attention due to the potential for efficient high-field DNP as well as a general interest in the underlying principles that enable the Overhauser effect in small molecules. We recently reported the observation of 1H and 2H Overhauser effects in H3C- or D3C-functionalized Blatter radical analogs, which we presumed to be caused by methyl rotation. In this work, we look at the mechanism for methyl-driven Overhauser DNP in greater detail, considering methyl librations and tunneling in addition to classical rotation. We predict the temperature dependence of these mechanisms using density functional theory and spin dynamics simulations. Comparisons with results from ultralow-temperature magic angle spinning-DNP experiments revealed that cross-relaxation at temperatures above 60 K originates from both libration and rotation, while librations dominate at lower temperatures. Due to the zero-point vibrational nature of these motions, they are not quenched by very low temperatures, and methyl-driven Overhauser DNP is expected to increase in efficiency down to 0 K, predominantly due to increases in nuclear relaxation times.

Probing a Hydrogen-π Interaction Involving a Trapped Water Molecule in the Solid State

The detection and characterization of trapped water molecules in chemical entities and biomacromolecules remains a challenging task for solid materials. We herein present proton-detected solid-state Nuclear Magnetic Resonance (NMR) experiments at 100 kHz magic-angle spinning and at high static magnetic-field strengths (28.2 T) enabling the detection of a single water molecule fixed in the calix[4]arene cavity of a lanthanide complex by a combination of three types of non-covalent interactions. The water proton resonances are detected at a chemical-shift value close to zero ppm, which we further confirm by quantum-chemical calculations. Density Functional Theory calculations pinpoint to the sensitivity of the proton chemical-shift value for hydrogen-π interactions. Our study highlights how proton-detected solid-state NMR is turning into the method-of-choice in probing weak non-covalent interactions driving a whole branch of molecular-recognition events in chemistry and biology.

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 ¹H-X (X = ¹⁴N, ³⁵Cl) interactions. This experiment can be used to acquire 2D ¹H{X}-HMQC filtered ¹H-¹H correlation as well as 2D ¹H-X HMQC spectra. Powder forms of dopamine·HCl and l-histidine·HCl·H₂O 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·H₂O, 2D peaks associated with ¹⁴N-¹H-¹H and ³⁵Cl-¹H-¹H distances of up to 4.4 and 3.9 Å have been detected. This is further corroborated by the observation of 2D peaks corresponding to ¹⁴N-¹H-¹H and ³⁵Cl-¹H-¹H 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.

Ambiguous structure determination from powder data: four different structural models of 4,11-di-fluoro-quinacridone with similar X-ray powder patterns, fit to the PDF, SSNMR and DFT-D

Four different structural models, which all fit the same X-ray powder pattern, were obtained in the structure determination of 4,11-di-fluoro-quinacridone (C₂₀H₁₀N₂O₂F₂) from unindexed X-ray powder data by a global fit. The models differ in their lattice parameters, space groups, Z, Z', molecular packing and hydrogen bond patterns. The molecules form a criss-cross pattern in models A and B, a layer structure built from chains in model C and a criss-cross arrangement of dimers in model D. Nevertheless, all models give a good Rietveld fit to the experimental powder pattern with acceptable R-values. All molecular geometries are reliable, except for model D, which is slightly distorted. All structures are crystallochemically plausible, concerning density, hydrogen bonds, intermolecular distances etc. All models passed the checkCIF test without major problems; only in model A a missed symmetry was detected. All structures could have probably been published, although 3 of the 4 structures were wrong. The investigation, which of the four structures is actually the correct one, was challenging. Six methods were used: (1) Rietveld refinements, (2) fit of the crystal structures to the pair distribution function (PDF) including the refinement of lattice parameters and atomic coordinates, (3) evaluation of the colour, (4) lattice-energy minimizations with force fields, (5) lattice-energy minimizations by two dispersion-corrected density functional theory methods, and (6) multinuclear CPMAS solid-state NMR spectroscopy (¹H, ¹³C, ¹⁹F) including the comparison of calculated and experimental chemical shifts. All in all, model B (perhaps with some disorder) can probably be considered to be the correct one. This work shows that a structure determination from limited-quality powder data may result in totally different structural models, which all may be correct or wrong, even if they are chemically sensible and give a good Rietveld refinement. Additionally, the work is an excellent example that the refinement of an organic crystal structure can be successfully performed by a fit to the PDF, and the combination of computed and experimental solid-state NMR chemical shifts can provide further information for the selection of the most reliable structure among several possibilities.

Solid-State NMR: Methods for Biological Solids

In the last two decades, solid-state nuclear magnetic resonance (ssNMR) spectroscopy has transformed from a spectroscopic technique investigating small molecules and industrial polymers to a potent tool decrypting structure and underlying dynamics of complex biological systems, such as membrane proteins, fibrils, and assemblies, in near-physiological environments and temperatures. This transformation can be ascribed to improvements in hardware design, sample preparation, pulsed methods, isotope labeling strategies, resolution, and sensitivity. The fundamental engagement between nuclear spins and radio-frequency pulses in the presence of a strong static magnetic field is identical between solution and ssNMR, but the experimental procedures vastly differ because of the absence of molecular tumbling in solids. This review discusses routinely employed state-of-the-art static and MAS pulsed NMR methods relevant for biological samples with rotational correlation times exceeding 100's of nanoseconds. Recent developments in signal filtering approaches, proton methodologies, and multiple acquisition techniques to boost sensitivity and speed up data acquisition at fast MAS are also discussed. Several examples of protein structures (globular, membrane, fibrils, and assemblies) solved with ssNMR spectroscopy have been considered. We also discuss integrated approaches to structurally characterize challenging biological systems and some newly emanating subdisciplines in ssNMR spectroscopy.

Proton-phosphorous connectivities revealed by high-resolution proton-detected solid-state NMR

Proton-detected solid-state NMR enables atomic-level insight in solid-state reactions, for instance in heterogeneous catalysis, which is fundamental for deciphering chemical reaction mechanisms. We herein introduce a phosphorus-31 radiofrequency channel in...

Determining the Three-Dimensional Structures of Silica-Supported Metal Complexes from the Ground Up

The immobilization of molecularly precise metal complexes to substrates, such as silica, provides an attractive platform for the design of active sites in heterogeneous catalysts. Specific steric and electronic variations of the ligand environment enable the development of structure-activity relationships and the knowledge-driven design of catalysts. At present, however, the three-dimensional environment of the precatalyst, much less the active site, is generally not known for heterogeneous single-site catalysts. We explored the degree to which NMR-based surface-to-complex interatomic distances could be used to solve the three-dimensional structures of three silica-supported metal complexes. The structure solution revealed unexpected features related to the environment around the metal that would be difficult to discern otherwise. This approach appears to be highly robust and, due to its simplicity, is readily applied to most single-site catalysts with little extra effort.

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.

NMR spectroscopy probes microstructure, dynamics and doping of metal halide perovskites

Solid-state magic-angle spinning NMR spectroscopy is a powerful technique to probe atomic-level microstructure and structural dynamics in metal halide perovskites. It can be used to measure dopant incorporation, phase segregation, halide mixing, decomposition pathways, passivation mechanisms, short-range and long-range dynamics, and other local properties. This Review describes practical aspects of recording solid-state NMR data on halide perovskites and how these afford unique insights into new compositions, dopants and passivation agents. We discuss the applicability, feasibility and limitations of ¹H, ¹³C, ¹⁵N, ¹⁴N, ¹³³Cs, ⁸⁷Rb, ³⁹K, ²⁰⁷Pb, ¹¹⁹Sn, ¹¹³Cd, ²⁰⁹Bi, ¹¹⁵In, ¹⁹F and ²H NMR in typical experimental scenarios. We highlight the pivotal complementary role of solid-state mechanosynthesis, which enables highly sensitive NMR studies by providing large quantities of high-purity materials of arbitrary complexity and of chemical shifts calculated using density functional theory. We examine the broader impact of solid-state NMR on materials research and how its evolution over seven decades has benefitted structural studies of contemporary materials such as halide perovskites. Finally, we summarize some of the open questions in perovskite optoelectronics that could be addressed using solid-state NMR. We, thereby, hope to stimulate wider use of this technique in materials and optoelectronics research.

Accurate fragment-based 51-V chemical shift predictions in molecular crystals

Nuclear magnetic resonance (NMR) spectroscopy plays a crucial role in determining molecular structure for complex biological and pharmaceutical compounds. NMR investigations are increasingly reliant on computation for mapping spectral features to chemical structures. Here we benchmark the accuracy of fragment-based ⁵¹V chemical shielding tensor calculations using a training set comprised of 10 biologically and pharmaceutically relevant oxovanadium complexes. Using our self-consistent reproduction of the Madelung potential (SCRMP) electrostatic embedding model, we demonstrate comparable performance between fragment methods and computationally demanding cluster-based techniques. Specifically, fragment methods employing hybrid density functionals are capable of reproducing the experimental ⁵¹V isotropic chemical shifts with a training set rms error of ~9 ​ppm, representing a 20% improvement over traditional plane wave techniques. We provide training set-derived linear regression models for mapping the absolute shieldings obtained from computation to the experimentally determined chemical shifts using four common density functionals; PBE0, B3LYP, PBE, and BLYP. Finally, we establish the utility of fragment methods and the reported regression parameters examining four oxovanadium structures excluded from the training set including the tetracoordinate oxovanadium silicate [Formula: see text] , VO¹⁵NGlySalbz which contains redox-active ligands, and the solid-state form of the common ⁵¹V NMR reference compound VOCl₃.

Recent progress in dipolar recoupling techniques under fast MAS in solid-state NMR spectroscopy

With the recent advances in NMR hardware and probe design technology, magic-angle spinning (MAS) rates over 100 ​kHz are accessible now, even on commercial solid NMR probes. Under such fast MAS conditions, excellent spectral resolution has been achieved by efficient suppression of anisotropic interactions, which also opens an avenue to the proton-detected NMR experiments in solids. Numerous methods have been developed to take full advantage of fast MAS during the last decades. Among them, dipolar recoupling techniques under fast MAS play vital roles in the determination of the molecular structure and dynamics, and are also key elements in multi-dimensional correlation NMR experiments. Herein, we review the dipolar recoupling techniques, especially those developed in the past two decades for fast-to-ultrafast MAS conditions. A major focus for our discussion is the ratio of RF field strength (in frequency) to MAS frequency, ν₁/ν_r, in different pulse sequences, which determines whether these dipolar recoupling techniques are suitable for NMR experiments under fast MAS conditions. Systematic comparisons are made among both heteronuclear and homonuclear dipolar recoupling schemes. In addition, the schemes developed specially for proton-detection NMR experiments under ultrafast MAS conditions are highlighted as well.

Recent advances in probing host–guest interactions with solid state nuclear magnetic resonance

A recent update on how solid state NMR has aided the interpretation and understanding of host–guest interactions in the field of supramolecular assemblies is provided.

Accuracy of 1H-1H distances measured using frequency selective recoupling and fast magic-angle spinning

Selective recoupling of protons (SERP) is a method to selectively and quantitatively measure magnetic dipole-dipole interaction between protons and, in turn, the proton-proton distance in solid-state samples at fast magic-angle spinning. We present a bimodal operator-based Floquet approach to describe the numerically optimized SERP recoupling sequence. The description calculates the allowed terms in the first-order effective Hamiltonian, explains the origin of selectivity during recoupling, and shows how different terms are modulated as a function of the radio frequency amplitude and the phase of the sequence. Analytical and numerical simulations have been used to evaluate the effect of higher-order terms and offsets on the polarization transfer efficiency and quantitative distance measurement. The experimentally measured ¹H-¹H distances on a fully protonated thymol sample are ∼10%-15% shorter than those reported from diffraction studies. A semi-quantitative model combined with extensive numerical simulations is used to rationalize the effect of the third-spin and the role of different parameters in the experimentally observed shorter distances. Measurements at high magnetic fields improve the match between experimental and diffraction distances. The measurement of ¹H-¹H couplings at offsets different from the SERP-offset has also been explored. Experiments were also performed on a perdeuterated ubiquitin sample to demonstrate the feasibility of simultaneously measuring multiple quantitative distances and to evaluate the accuracy of the measured distance in the absence of multispin effects. The estimation of proton-proton distances provides a boost to structural characterization of small pharmaceuticals and biomolecules, given that the positions of protons are generally not well defined in x-ray structures.

Hydrogen-Mediated Noncovalent Interactions in Solids: What Can NMR Crystallography Tell About?

Hydrogen atoms play a crucial role in the aggregation of organic (bio)molecules through diverse number of noncovalent interactions that they mediate, such as electrostatic in proton transfer systems, hydrogen bonding, and CH-π interactions, to mention only the most prominent. To identify and adequately describe such low-energy interactions, increasingly sensitive methods have been developed over time, among which quantum chemical computations have witnessed impressive advances in recent years. For reaching the present state-of-the-art, computations had to rely on a pool of relevant experimental data, needed at least for validation, if not also for other purposes. In the case of molecular crystals, the best illustration for the synergy between computations and experiment is given by the so-called NMR crystallography approach. Originally designed to increase the confidence level in crystal structure determination of organic compounds from powders, NMR crystallography is able now to offer also a wealth of information regarding the noncovalent interactions that drive molecules to pack in a given crystalline pattern or another. This is particularly true for the noncovalent interactions which depend on the exact location of labile hydrogen atoms in the system: in such cases, NMR crystallography represents a valuable characterization tool, in some cases complementing even the standard single-crystal X-ray diffraction technique. A concise introduction in the field is made in this mini-review, which is aimed at providing a comprehensive picture with respect to the current accuracy level reached by NMR crystallography in the characterization of hydrogen-mediated noncovalent interactions in organic solids. Different types of practical applications are illustrated with the example of molecular crystals studied by our research group, but references to other representative developments reported in the literature are also made. By summarizing the major concepts and methodological progresses, the present work is also intended to be a guide to the practical potential of this relatively recent analytical tool for the scientists working in areas where crystal engineering represents the main approach for rational design of novel materials.

Can proton-proton recoupling in fully protonated solids provide quantitative, selective and efficient polarization transfer?

Dipolar recoupling sequences have been used to probe spatial proximity of nuclear spins and were traditionally designed to probe rare spins such as ¹³C and/or ¹⁵N nuclei. The multi-spin dipolar-coupling network of the rare spins is weak due to smaller couplings and large chemical shift dispersion. Therefore, the recoupling approaches were tailored to design offset compensated or broadband sequences. In contrast, protons have a substantially stronger dipolar-coupling network and much narrower chemical shift range. Broadband recoupling sequences such as radio-frequency driven recoupling (RFDR), back-to-back (BABA), and lab frame proton-proton spin diffusion have been routinely used to characterize the structures of protein/macromolecules and small molecules. Recently selective ¹H-¹H recoupling sequences have been proposed that combine chemical shift offset of the resolved proton spectrum (at fast MAS) with first- and second-order dipolar recoupling Hamiltonians to obtain quantitative and qualitative proton distances, respectively. Herein, we evaluate the performances of broadband and selective proton recoupling sequences such as finite pulse RFDR (fp-RFDR), band-selective spectral spin diffusion (BASS-SD), second-order cross-polarization (SOCP), and selective recoupling of proton (SERP) in terms of the selectivity and efficiency of ¹H-¹H polarization transfers in a dense network of proton spins and explore the possibility of measuring ¹H-¹H distances. We use theoretical considerations, numerical simulations, and experiments to support the distinct advantages and disadvantages of each recoupling sequence. Experiments were performed on L-histidine.HCl.H₂O at a MAS frequency of 71.43 kHz. This study rationalizes the proper selection of ¹H-¹H recoupling sequences when working with fully protonated solids.

Molecular Mechanism of Crystalline-to-Amorphous Conversion of Pharmaceutical Solids from 19F Magic Angle Spinning NMR

Crystalline and amorphous materials usually possess distinct physicochemical properties due to major variations in long-range and local molecular packings. Enhanced fundamental knowledge of the molecular details of crystalline-to-amorphous interconversions is necessary to correlate the intermolecular structure to material properties and functions. While crystal structures can be readily obtained by X-ray crystallography, the microstructure of amorphous materials has rarely been explored due to a lack of high-resolution techniques capable of probing local molecular structures. Moreover, there is increasing interest in understanding the molecular nature of amorphous solids in pharmaceutical sciences due to the widespread utilization of amorphous active pharmaceutical ingredients (APIs) in pharmaceutical development for solubility and bioavailability enhancement. In this study, we explore multidimensional ¹³C and ¹⁹F magic angle spinning (MAS) NMR spectroscopy to study the molecular packing of amorphous posaconazole (POSA) in conjunction with the crystalline counterpart. Utilizing methods integrating homonuclear and heteronuclear ¹H, ¹³C, and ¹⁹F correlation spectroscopy and atomic ¹⁹F-to-¹³C distance measurements, we identified the major differences in molecular packing between crystalline and amorphous POSA. The intermolecular "head-to-head" interaction along the molecule's major axis, as well as the "head-to-tail" molecular packing perpendicular to the major axis in POSA crystals, was recapitulated by MAS NMR. Furthermore, critical intermolecular distances in the crystal lattice were determined. Most importantly, the head-to-tail contact of two neighboring molecules was found to be preserved in amorphous POSA, suggesting localized molecular order, whereas crucial interactions for head-to-head packing are absent in the amorphous form resulting in long-range disorder. Our study, likely one of the first documented examples, provides molecular-level structural details to understand the molecular mechanism of crystalline-to-amorphous conversion of fluorine-containing drug substances occurring in drug processing and development and establish a high-resolution experimental protocol for investigating amorphous materials.

A brief introduction to the basics of NMR spectroscopy and selected examples of its applications to materials characterization

Nuclear magnetic resonance (NMR) spectroscopy is an analytical technique that gives information on the local magnetic field around atomic nuclei. Since the local magnetic field of the nucleus is directly influenced by such features of the molecular structure as constitution, configuration, conformation, intermolecular interactions, etc., NMR can provide exhaustive information on the chemical structure, which is unrivaled by any other analytical method. Starting from the 1950s, NMR spectroscopy first revolutionized organic chemistry and became an indispensable tool for the structure elucidation of small, soluble molecules. As the technique evolved, NMR rapidly conquered other disciplines of chemical sciences. When the analysis of macromolecules and solids also became feasible, the technique turned into a staple in materials characterization, too. All aspects of NMR spectroscopy, including technical and technological development, as well as its applications in natural sciences, have been growing exponentially since its birth. Hence, it would be impossible to cover, or even touch on, all topics of importance related to this versatile analytical tool. In this tutorial, we aim to introduce the reader to the basic principles of NMR spectroscopy, instrumentation, historical development and currently available brands, practical cost aspects, sample preparation, and spectrum interpretation. We show a number of advanced techniques relevant to materials characterization. Through a limited number of examples from different fields of materials science, we illustrate the immense scope of the technique in the analysis of materials. Beyond our inherently limited introduction, an ample list of references should help the reader to navigate further in the field of NMR spectroscopy.

NMR crystallography of molecular organics

Developments of NMR methodology to characterise the structures of molecular organic structures are reviewed, concentrating on the previous decade of research in which density functional theory-based calculations of NMR parameters in periodic solids have become widespread. With a focus on demonstrating the new structural insights provided, it is shown how "NMR crystallography" has been used in a spectrum of applications from resolving ambiguities in diffraction-derived structures (such as hydrogen atom positioning) to deriving complete structures in the absence of diffraction data. As well as comprehensively reviewing applications, the different aspects of the experimental and computational techniques used in NMR crystallography are surveyed. NMR crystallography is seen to be a rapidly maturing subject area that is increasingly appreciated by the wider crystallographic community.

The origin of negative cross-peaks in proton-spin diffusion spectrum of fully protonated solids at fast MAS: Coherent or incoherent effect?

Spin-diffusion (SD) is amongst the first methods proposed to spatially transfer polarization between dipolar-coupled nuclear spins. Lab-frame SD has proved particularly useful in structural characterization of a large variety of molecules. During SD, the rate of magnetization transfer between the two nuclei depends on the square of the dipolar coupling and the zero-quantum lineshape of the two spins. The relative sign of the diagonal and cross-peaks is determined by the spin part of the dipolar Hamiltonian. Practically, SD experiments are used in two ways: (a) SD transfer amongst only the protons (known as proton spin-diffusion or PSD) and b) SD amongst rare nuclei, coupled to a strong proton bath, known as proton driven spin-diffusion (PDSD). It is well established that the diagonal and cross-peaks have the same sign during SD based polarization transfer. 2D PSD experiments recorded on Histidine.HCl.H₂O sample at fast magic angle spinning (MAS) show that some of the cross-peaks in the 2D spectrum are negative with respect to the diagonal peaks. Cross-relaxation due to stochastic motion is generally believed to give rise to such negative peaks. Herein, we use theoretical calculations, numerical simulations and experiments to show that the origin of the negative cross-peaks in PSD spectrum is due to coherent interactions. The origin of negative peaks can be specifically ascribed to a four spin, double-flip-double flop term, in the third-order Hamiltonian. These terms become the dominant terms at fast spinning when additional - conditions are satisfied.

Ultrafast 1H MAS NMR Crystallography for Natural Abundance Pharmaceutical Compounds

Magic angle spinning (MAS) NMR is a powerful method for the study of pharmaceutical compounds, and probes with spinning frequencies above 100 kHz enable an atomic-resolution analysis of sub-micromole quantities of fully protonated solids. Here, we present an ultrafast NMR crystallography approach for structural characterization of organic solids at MAS frequencies of 100-111 kHz. We assess the efficiency of ¹H-detected experiments in the solid state and demonstrate the utility of 2D and 3D homo- and heteronuclear correlation spectra for resonance assignments. These experiments are demonstrated for an amino acid, U-¹³C,¹⁵N-histidine, and also for the significantly larger, natural product Posaconazole, an antifungal compound investigated at natural abundance. Our results illustrate the power for characterizing organic molecules, enabled by exploiting the increased ¹H resolution and sensitivity at MAS frequencies above 100 kHz.

Crystal structure determination of an elusive methanol solvate – hydrate of catechin using crystal structure prediction and NMR crystallography

A useful short-cut was developed to limit the number of molecular conformations that need to be regarded in crystal structure prediction calculations, which led to the crystal structure determination of new methanol solvate – hydrate of catechin.

Understanding hydrogen-bonding structures of molecular crystals via electron and NMR nanocrystallography

Understanding hydrogen-bonding networks in nanocrystals and microcrystals that are too small for X-ray diffractometry is a challenge. Although electron diffraction (ED) or electron 3D crystallography are applicable to determining the structures of such nanocrystals owing to their strong scattering power, these techniques still lead to ambiguities in the hydrogen atom positions and misassignments of atoms with similar atomic numbers such as carbon, nitrogen, and oxygen. Here, we propose a technique combining ED, solid-state NMR (SSNMR), and first-principles quantum calculations to overcome these limitations. The rotational ED method is first used to determine the positions of the non-hydrogen atoms, and SSNMR is then applied to ascertain the hydrogen atom positions and assign the carbon, nitrogen, and oxygen atoms via the NMR signals for ¹H, ¹³C, ¹⁴N, and ¹⁵N with the aid of quantum computations. This approach elucidates the hydrogen-bonding networks in l-histidine and cimetidine form B whose structure was previously unknown.

Electron diffraction can be used to determine nanocrystal structures, but is unsuitable for locating hydrogen atoms. Here the authors combine electron diffraction, solid-state NMR and first-principles calculations to resolve the crystal structures and hydrogen-bonding networks of L-histidine and cimetidine form B.

Linear-scaling ab initio simulations of spin diffusion in rotating solids

We investigated the utility of locally restricting the basis sets involved in low-order correlations in Liouville space (LCL) calculations of spin diffusion. Using well-known classical models of spin diffusion, we describe a rationale for selecting the optimal basis set for such calculations. We then show that the use of these locally restricted basis sets provides the same computational accuracy as the full LCL set while reducing the computational time by several orders of magnitude. Speeding up the calculations also enables us to use higher maximum spin orders and increase the computational accuracy. Furthermore, unlike exact and full LCL calculations, locally restricted LCL calculations scale linearly with the system size and should thus enable the ab initio study of spin diffusion in spin systems containing several thousand spins.

Quantitative 1H-1H Distances in Protonated Solids by Frequency-Selective Recoupling at Fast Magic Angle Spinning NMR

Nuclear magnetic resonance (NMR) spectroscopy of protons in protonated solids is challenging. Fast magic angle spinning (MAS) and homonuclear decoupling schemes, in conjunction, with high magnetic fields have improved the proton resolution. However, experiments to quantitatively measure ¹H-¹H distances still remain elusive due to the dense proton-proton dipolar coupling network. A novel MAS solid-state NMR pulse sequence is proposed to selectively recouple and measure interproton distances in protonated samples. The phase-modulated sequence combined with a judicious choice of transmitter frequency is used to measure quantitative ¹H-¹H distances on the order of 3 Å in l-histidine·HCl·H₂O, despite the presence of other strongly coupled protons. This method provides a major boost to NMR crystallography approaches for structural determination of pharmaceutical molecules by directly measuring ¹H-¹H distances. The band-selective nature of the sequence also enables observation of selective ¹H-¹H correlations (e.g., H^N-H^N/H^N-H^α/Η^Ν-Η^Methyl) in peptides and proteins, which should serve as useful restraints in structure determination.

Structural Analyses of Alanine Trimer and Tetramer Crystals with Antiparallel and Parallel β-Sheet Structures Using Solid-State 1H Spin-Diffusion 2D Correlation NMR Spectroscopy

Poly-l-alanine (PLA) sequences are key elements of the crystalline domains of spider dragline and wild silkworm silks. In the present work, ¹H spin-diffusion two-dimensional (2D) correlation NMR spectra were observed for selectively deuterated (Ala)₃ and (Ala)₄ crystals to develop the analytical method for the structure of PLA sequences. The build-up curves of the cross peaks for three kinds of ¹H pairs in selectively deuterated (Ala)₃ and (Ala)₄ crystals were observed to obtain spin-diffusion rate constant k _{j, k} from relaxation master equations P _{i, j}(τ_m). The k _{j, k} values subsequently lead to effective interproton distance r _{j, k}^eff (obs) values for individual proton-proton pairs, which include intra- and intermolecular contributions. The r _{j, k}^eff (obs) values were compared to r _{j, k}^eff (calc) values obtained from the experimentally determined atomic coordinates of antiparallel (AP) β-sheet (Ala)₃ and (Ala)₄ and parallel (P) β-sheet of (Ala)₃ and (Ala)₄ crystals. The agreement between the r _{j, k}^eff (obs) and r _{j, k}^eff (calc) values was good for AP β-sheet (Ala)₃ and (Ala)₄ crystals but poor for P β-sheet (Ala)₃ and (Ala)₄ crystals. These deviations were obtained from the interproton distances of the interchain contributions due to different packing arrangements. The packing arrangements of the PLA region are important when considering the relevant structure and the mechanical properties of silks.

DNP-enhanced solid-state NMR spectroscopy of active pharmaceutical ingredients

Solid-state NMR spectroscopy has become a valuable tool for the characterization of both pure and formulated active pharmaceutical ingredients (APIs). However, NMR generally suffers from poor sensitivity that often restricts NMR experiments to nuclei with favorable properties, concentrated samples, and acquisition of one-dimensional (1D) NMR spectra. Here, we review how dynamic nuclear polarization (DNP) can be applied to routinely enhance the sensitivity of solid-state NMR experiments by one to two orders of magnitude for both pure and formulated APIs. Sample preparation protocols for relayed DNP experiments and experiments on directly doped APIs are detailed. Numerical spin diffusion models illustrate the dependence of relayed DNP enhancements on the relaxation properties and particle size of the solids and can be used for particle size determination when the other factors are known. We then describe the advanced solid-state NMR experiments that have been enabled by DNP and how they provide unique insight into the molecular and macroscopic structure of APIs. For example, with large sensitivity gains provided by DNP, natural isotopic abundance, ¹³ C-¹³ C double-quantum single-quantum homonuclear correlation NMR spectra of pure APIs can be routinely acquired. DNP also enables solid-state NMR experiments with unreceptive quadrupolar nuclei such as ² H, ¹⁴ N, and ³⁵ Cl that are commonly found in APIs. Applications of DNP-enhanced solid-state NMR spectroscopy for the molecular level characterization of low API load formulations such as commercial tablets and amorphous solid dispersions are described. Future perspectives for DNP-enhanced solid-state NMR experiments on APIs are briefly discussed.

Does Z' equal 1 or 2? Enhanced powder NMR crystallography verification of a disordered room temperature crystal structure of a p38 inhibitor for chronic obstructive pulmonary disease

The crystal structure of the Form A polymorph of N-cyclopropyl-3-fluoro-4-methyl-5-[3-[[1-[2-[2-(methylamino)ethoxy]phenyl]cyclopropyl]amino]-2-oxo-pyrazin-1-yl]benzamide (i.e., AZD7624), determined using single-crystal X-ray diffraction (scXRD) at 100 K, contains two molecules in the asymmetric unit (Z' = 2) and has regions of local static disorder. This substance has been in phase IIa drug development trials for the treatment of chronic obstructive pulmonary disease, a disease which affects over 300 million people and contributes to nearly 3 million deaths annually. While attempting to verify the crystal structure using nuclear magnetic resonance crystallography (NMRX), we measured ¹³C solid-state NMR (SSNMR) spectra at 295 K that appeared consistent with Z' = 1 rather than Z' = 2. To understand this surprising observation, we used multinuclear SSNMR (¹H, ¹³C, ¹⁵N), gauge-including projector augmented-wave density functional theory (GIPAW DFT) calculations, crystal structure prediction (CSP), and powder XRD (pXRD) to determine the room temperature crystal structure. Due to the large size of AZD7624 (ca. 500 amu, 54 distinct ¹³C environments for Z' = 2), static disorder at 100 K, and (as we show) dynamic disorder at ambient temperatures, NMR spectral assignment was a challenge. We introduce a method to enhance confidence in NMR assignments by comparing experimental ¹³C isotropic chemical shifts against site-specific DFT-calculated shift distributions established using CSP-generated crystal structures. The assignment and room temperature NMRX structure determination process also included measurements of ¹³C shift tensors and the observation of residual dipolar coupling between ¹³C and ¹⁴N. CSP generated ca. 90 reasonable candidate structures (Z' = 1 and Z' = 2), which when coupled with GIPAW DFT results, room temperature pXRD, and the assigned SSNMR data, establish Z' = 2 at room temperature. We find that the polymorphic Form A of AZD7624 is maintained at room temperature, although dynamic disorder is present on the NMR timescale. Of the CSP-generated structures, 2 are found to be fully consistent with the SSNMR and pXRD data; within this pair, they are found to be structurally very similar (RMSD₁₆ = 0.30 Å). We establish that the CSP structure in best agreement with the NMR data possesses the highest degree of structural similarity with the scXRD-determined structure (RMSD₁₆ = 0.17 Å), and has the lowest DFT-calculated energy amongst all CSP-generated structures with Z' = 2.

Improved resolution and simplification of the spin-diffusion-based NMR method for the structural analysis of mixed-linker MOFs

Nuclear magnetic resonance spectroscopy combined with modeling represents a powerful tool for the structural analysis of heterogeneous materials. In this contribution we describe an upgraded method, particularly suited for the structural analysis of mixed-linker metal-organic framework materials, which is based on the measurement and modeling of proton spin diffusion among constituents. We tested the method on a UiO-66-type metal-organic material, in which the organic building units were 1,4-benzenedicarboxylate and trans-1,4-cyclohexanedicarboxylate anions distributed within the framework in an unknown manner. We showed that resolution of the signals of different building units could be significantly enhanced by the carbon-detected version of the proton spin-diffusion measurement. Because this kind of measurement is much more time consuming than the proton-detected measurement and because one has to carry out several two-dimensional measurements to extract spin-diffusion curves, we inspected the possibility of reducing the number of such measurements. This could be done by limiting the analysis to short mixing times, for which, as shown in this contribution, linear approximation is valid. When working in the linear regime, only a few experimental points are needed to determine the slope of spin-diffusion curves. Usage of short spin-diffusion mixing times significantly shortened the total measurement time and also markedly simplified the modeling of spin-diffusion curves.

A 13C solid-state NMR investigation of four cocrystals of caffeine and theophylline

We report an analysis of the 13C solid-state NMR chemical shift data in a series of four cocrystals involving two active pharmaceutical ingredient (API) mimics (caffeine and theophylline) and two diacid coformers (malonic acid and glutaric acid). Within this controlled set, we make comparisons of the isotropic chemical shifts and the principal values of the chemical shift tensor. The dispersion at 14.1 T (600 MHz 1H) shows crystallographic splittings in some of the resonances in the magic angle spinning spectra. By comparing the isotropic chemical shifts of individual C atoms across the four cocrystals, we are able to identify pronounced effects on the local electronic structure at some sites. We perform a similar analysis of the principal values of the chemical shift tensors for the anisotropic C atoms (most of the ring C atoms for the API mimics and the carbonyl C atoms of the diacid coformers) and link them to differences in the known crystal structures. We discuss the future prospects for extending this type of study to incorporate the full chemical shift tensor, including its orientation in the crystal frame of reference.

Enhanced NMR Discrimination of Pharmaceutically Relevant Molecular Crystal Forms through Fragment-Based Ab Initio Chemical Shift Predictions

Chemical shift prediction plays an important role in the determination or validation of crystal structures with solid-state nuclear magnetic resonance (NMR) spectroscopy. One of the fundamental theoretical challenges lies in discriminating variations in chemical shifts resulting from different crystallographic environments. Fragment-based electronic structure methods provide an alternative to the widely used plane wave gauge-including projector augmented wave (GIPAW) density functional technique for chemical shift prediction. Fragment methods allow hybrid density functionals to be employed routinely in chemical shift prediction, and we have recently demonstrated appreciable improvements in the accuracy of the predicted shifts when using the hybrid PBE0 functional instead of generalized gradient approximation (GGA) functionals like PBE. Here, we investigate the solid-state ¹³C and ¹⁵N NMR spectra for multiple crystal forms of acetaminophen, phenobarbital, and testosterone. We demonstrate that the use of the hybrid density functional instead of a GGA provides both higher accuracy in the chemical shifts and increased discrimination among the different crystallographic environments. Finally, these results also provide compelling evidence for the transferability of the linear regression parameters mapping predicted chemical shieldings to chemical shifts that were derived in an earlier study.

Five decades of homonuclear dipolar decoupling in solid-state NMR: Status and outlook

It has been slightly more than fifty years since the first homonuclear spin decoupling scheme, Lee-Goldburg decoupling, was proposed for removing homonuclear dipolar interactions in solid-state nuclear magnetic resonance. A family of such schemes has made observation of high-resolution NMR spectra of abundant spins possible in various applications in solid state. This review outlines the strategies used in this field and the future prospects of homonuclear spin decoupling in solid-state NMR.

High resolution triple resonance micro magic angle spinning NMR spectroscopy of nanoliter sample volumes

To be able to study mass-limited samples and small single crystals, a triple resonance micro-magic angle spinning (μMAS) probehead for the application of high-resolution solid-state NMR of nanoliter samples was developed. Due to its excellent rf performance this allows us to explore the limits of proton NMR resolution in strongly coupled solids. Using homonuclear decoupling we obtain unprecedented (1)H linewidths for a single crystal of glycine (Δν(CH2) = 0.14 ppm) at high field (20 T) in a directly detected spectrum. The triple channel design allowed the recording of high-resolution μMAS (13)C-(15)N correlations of [U-(13)C-(15)N] arginine HCl and shows that the superior (1)H resolution opens the way for high-sensitivity inverse detection of heteronuclei even at moderate spinning speeds and rf-fields. Efficient decoupling leads to long coherence times which can be exploited in many correlation experiments.

Quantitative structure parameters from the NMR spectroscopy of quadrupolar nuclei

Nuclear magnetic resonance (NMR) spectroscopy is one of the most important characterization tools in chemistry, however, 3/4 of the NMR active nuclei are underutilized due to their quadrupolar nature. This short review centers on the development of methods that use solid-state NMR of quadrupolar nuclei for obtaining quantitative structural information. Namely, techniques using dipolar recoupling as well as the resolution afforded by double-rotation are presented for the measurement of spin–spin coupling between quadrupoles, enabling the measurement of internuclear distances and connectivities. Two-dimensional J-resolved-type experiments are then presented for the measurement of dipolar and J coupling, between spin-1/2 and quadrupolar nuclei as well as in pairs of quadrupolar nuclei. Select examples utilizing these techniques for the extraction of structural information are given. Techniques are then described that enable the fine refinement of crystalline structures using solely the electric field gradient tensor, measured using NMR, as a constraint. These approaches enable the solution of crystal structures, from polycrystalline compounds, that are of comparable quality to those solved using single-crystal diffraction.

Characterising the role of water in sildenafil citrate by NMR crystallography

A combination of solid-state NMR techniques, including ¹³C/¹H correlation, ²H magic-angle spinning NMR and first principles calculation are employed to characterise the role of water in different hydration states of sildenafil citrate.

Proton-Detected Solid-State NMR Spectroscopy of Bone with Ultrafast Magic Angle Spinning

While obtaining high-resolution structural details from bone is highly important to better understand its mechanical strength and the effects of aging and disease on bone ultrastructure, it has been a major challenge to do so with existing biophysical techniques. Though solid-state NMR spectroscopy has the potential to reveal the structural details of bone, it suffers from poor spectral resolution and sensitivity. Nonetheless, recent developments in magic angle spinning (MAS) NMR technology have made it possible to spin solid samples up to 110 kHz frequency. With such remarkable capabilities, (1)H-detected NMR experiments that have traditionally been challenging on rigid solids can now be implemented. Here, we report the first application of multidimensional (1)H-detected NMR measurements on bone under ultrafast MAS conditions to provide atomistic-level elucidation of the complex heterogeneous structure of bone. Our investigations demonstrate that two-dimensional (1)H/(1)H chemical shift correlation spectra for bone are obtainable using fp-RFDR (finite-pulse radio-frequency-driven dipolar recoupling) pulse sequence under ultrafast MAS. Our results infer that water exhibits distinct (1)H-(1)H dipolar coupling networks with the backbone and side-chain regions in collagen. These results show the promising potential of proton-detected ultrafast MAS NMR for monitoring structural and dynamic changes caused by mechanical loading and disease in bone.

Macroscopic nuclear spin diffusion constants of rotating polycrystalline solids from first-principles simulation

A method for quantitatively calculating nuclear spin diffusion constants directly from crystal structures is introduced. This approach uses the first-principles low-order correlations in Liouville space (LCL) method to simulate spin diffusion in a box, starting from atomic geometry and including both magic-angle spinning (MAS) and powder averaging. The LCL simulations are fit to the 3D diffusion equation to extract quantitative nuclear spin diffusion constants. We demonstrate this method for the case of (1)H spin diffusion in ice and L-histidine, obtaining diffusion constants that are consistent with literature values for (1)H spin diffusion in polymers and that follow the expected trends with respect to magic-angle spinning rate and the density of nuclear spins. In addition, we show that this method can be used to model (13)C spin diffusion in diamond and therefore has the potential to provide insight into applications such as the transport of polarization in non-protonated systems.