NMR crystallography: structure and properties of materials from solid-state nuclear magnetic resonance observables

Computational Protocol for the Spectral Assignment of NMR Resonances in Covalent Organic Frameworks

Solid-state nuclear magnetic resonance spectroscopy is routinely used in the field of covalent organic frameworks to elucidate or confirm the structure of the synthesized samples and to understand dynamic phenomena. Typically this involves the interpretation and simulation of the spectra through the assumption of symmetry elements of the building units, hinging on the correct assignment of each line shape. To avoid misinterpretation resulting from library-based assignment without a theoretical basis incorporating the impact of the framework, this work proposes a first-principles computational protocol for the assignment of experimental spectra, which exploits the symmetry of the underlying building blocks for computational feasibility. In this way, this protocol accommodates the validation of previous experimental assignments and can serve to complement new NMR measurements.

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.

Carbon-13 chemical shift tensor measurements for nitrogen-dense compounds

This paper reports the principal values of the 13 C chemical shift tensors for five nitrogen-dense compounds (i.e., cytosine, uracil, imidazole, guanidine hydrochloride, and aminoguanidine hydrochloride). Although these are all fundamentally important compounds, the majority do not have 13 C chemical shift tensors reported in the literature. The chemical shift tensors are obtained from 1 H→13 C cross-polarization magic-angle spinning (CP/MAS) experiments that were conducted at a high field of 18.8 T to suppress the effects of 14 N-13 C residual dipolar coupling. Quantum chemical calculations using density functional theory are used to obtain the 13 C magnetic shielding tensors for these compounds. The best agreement with experiment arises from calculations using the hybrid functional PBE0 or the double-hybrid functional PBE0-DH, along with the triple-zeta basis sets TZ2P or pc-3, respectively, and intermolecular effects modeled using large clusters of molecules with electrostatic embedding through the COSMO approach. These measurements are part of an ongoing effort to expand the catalog of accurate 13 C chemical shift tensor measurements, with the aim of creating a database that may be useful for benchmarking the accuracy of quantum chemical calculations, developing nuclear magnetic resonance (NMR) crystallography protocols, or aiding in applications involving machine learning or data mining. This work was conducted at the National High Magnetic Field Laboratory as part of a 2-week school for introducing undergraduate students to practical laboratory experience that will prepare them for scientific careers or postgraduate studies.

Unveiling the topology of partially disordered micro-crystalline nitro-perylenediimide with X-aggregate stacking: an integrated approach

Profound knowledge of the molecular structure and supramolecular organization of organic molecules is essential to understand their structure-property relationships. Herein we demonstrate the packing arrangement of partially disordered nitro-perylenediimide (NO2-PDI), revealing that the perylenediimide units exhibit an X-shaped packing pattern. The packing of NO2-PDI is derived using a complementary approach that utilises solid-state NMR (ssNMR) and 3D electron diffraction (3D ED) techniques. Perylenediimide (PDI) molecules are captivating due to their high luminescence efficiency and optoelectronic properties, which are related to supramolecular self-assembly. Increasing the alkyl chain length on the imide substituent poses a more significant challenge in crystallizing the resulting molecule. In addition to the alkyl tails, other functional groups, like the nitro group attached as a bay substituent, can also cause disorder. Such heterogeneity could lead to diffuse scattering, which then complicates the interpretation of diffraction experiment data, where perfect periodicity is expected. As a result, there is an unmet need to develop a methodology for solving the structures of difficult-to-crystallize materials. A synergistic approach is utilised in this manuscript to understand the packing arrangement of the disordered material NO2-PDI by making use of 3D ED, ssNMR and density functional theory calculations (DFT). The combination of these experimental and theoretical approaches provides great promise in enabling the structural investigation of novel materials with customized properties across various applications, which are, due to the internal disorder, very difficult to study by diffraction techniques. By effectively addressing these challenges, our methodology opens up new avenues for material characterization, thereby driving exciting advancements in the field.

Simple and User-Friendly Methodology for Crystal Water Determination by Quantitative Proton NMR Spectroscopy in Deuterium Oxide

In drug research and development, knowledge of the precise structure of an active ingredient is crucial. However, it is equally important to know the water content of the drug molecule, particularly the number of crystal waters present in its structure. Such knowledge ensures the avoidance of drug dosage and formulation errors since the number of water molecules affects the physicochemical and pharmaceutical properties of the molecule. Several methods have been used for crystal water measurements of organic compounds, of which thermogravimetry and crystallography may be the most common ones. To the best of our knowledge, solution-state NMR spectroscopy has not been used for crystal water determination in deuterium oxide. Quantitative NMR (qNMR) method will be presented in the paper with a comparison of single-crystal X-ray diffraction and thermogravimetric analysis results. The qNMR method for water content measurement is straightforward, reproducible, and accurate, including measurement of 1H NMR spectrum before and after the addition of the analyte compound, and the result can be calculated after integration of the reference compound, analyte, and HDO signals using the given equation. In practical terms, there is no need for weighing the samples under study, which makes it simple and is a clear advantage to the current determination methods. In addition, the crystal structures of two model bisphosphonates used herein are reported: that of monopotassium etidronate dihydrate and monosodium zoledronate trihydrate.

A focus on detection of polymorphs by dynamic nuclear polarization solid-state nuclear magnetic resonance spectroscopy

Solid-state nuclear magnetic resonance (ssNMR) spectroscopy has found increasing application as a method for quantification and structure determination of solid forms (polymorphs) of organic solids and active pharmaceutical ingredients (APIs). However, ssNMR spectroscopy suffers from low sensitivity and resolution, making it challenging to detect dilute solid forms that may be present after recrystallization or reaction with co-formers. Cousin et al. (S. F. Cousin et al., Chem. Sci., 2023, https://doi.org/10.1039/D3SC02063K) have demonstrated that dynamic nuclear polarization (DNP) enhanced 13C cross-polarization (CP) saturation recovery experiments can be used to detect dilute polymorphic forms that are present within a mixture of solid forms. Enhancement of the NMR signal by DNP and differences in signal build-up rates for different polymorphs provide the sensitivity and contrast needed to resolve NMR signals from minor polymorphic forms. This method demonstrated by Cousin et al. should aid the discovery of solid drug forms.

Modulation of Rotational Dynamics in Halogen-Bonded Cocrystalline Solids

Dynamic processes are responsible for the functionality of a range of materials, biomolecules, and catalysts. We report a detailed systematic study of the modulation of methyl rotational dynamics via the direct and the indirect influence of noncovalent halogen bonds. For this purpose, a novel series of cocrystalline architectures featuring halogen bonds (XB) to tetramethylpyrazine (TMP) is designed and prepared using gas-phase, solution, and solid-state mechanochemical methods. Single-crystal X-ray diffraction reveals the capacity of molecular bromine as well as weak chloro-XB donors to act as robust directional structure-directing elements. Methyl rotational barriers (Ea) measured using variable-temperature deuterium solid-state NMR range from 3.75 ± 0.04 kJ mol-1 in 1,3,5-trichloro-2,4,6-trifluorobenzene·TMP to 7.08 ± 0.15 kJ mol-1 in 1,4-dichlorotetrafluorobenzene·TMP. Ea data for a larger series of TMP cocrystals featuring chloro-, bromo-, and iodo-XB donors are shown to be governed by a combination of steric and electronic factors. The average number of carbon-carbon close contacts to the methyl group is found to be a key steric metric capable of rationalizing the observed trends within each of the Cl, Br, and I series. Differences between each series are accounted for by considering the strength of the σ-hole on the XB donor. One possible route to modulating dynamics is therefore via designer cocrystals of variable stoichiometry, maintaining the core chemical features of interest between a given donor and acceptor while simultaneously modifying the number of carbon close contacts affecting dynamics. These principles may provide design opportunities to modulate more complex geared or cascade dynamics involving larger functional groups.

The Structure of Boron Monoxide

Boron monoxide (BO), prepared by the thermal condensation of tetrahydroxydiboron, was first reported in 1955; however, its structure could not be determined. With the recent attention on boron-based two-dimensional materials, such as borophene and hexagonal boron nitride, there is renewed interest in BO. A large number of stable BO structures have been computationally identified, but none are supported by experiments. The consensus is that the material likely forms a boroxine-based two-dimensional material. Herein, we apply advanced 11B NMR experiments to determine the relative orientations of B(B)O2 centers in BO. We find that the material is composed of D2h-symmetric O2B-BO2 units that organize to form larger B4O2 rings. Further, powder diffraction experiments additionally reveal that these units organize to form two-dimensional layers with a random stacking pattern. This observation is in agreement with earlier density functional theory (DFT) studies that showed B4O2-based structures to be the most stable.

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.

Solid-state NMR Spectroscopy of Iodine(I) Complexes

Solid-state NMR has been applied to a series of Barluenga-type iodine(I) [L-I-L]PF₆ (L=pyridine, 4-ethylpyridine, 4-dimethylaminopyridine, isoquinoline) complexes as their hexafluorophosphate salts, as well as their respective non-liquid ligands (L), their precursor silver(I) complexes, and the respective N-methylated pyridinium and quinolinium hexafluorophoshate salts. These results are compared and contrasted to the corresponding solution studies and single-crystal X-ray structures. As the first study of its kind on the solid-state NMR behavior of halogen(I) complexes, practical considerations are also discussed to encourage wider utilization of this technique in the future.

Zwitterionic or Not? Fast and Reliable Structure Determination by Combining Crystal Structure Prediction and Solid-State NMR

When it comes to crystal structure determination, computational approaches such as Crystal Structure Prediction (CSP) have gained more and more attention since they offer some insight on how atoms and molecules are packed in the solid state, starting from only very basic information without diffraction data. Furthermore, it is well known that the coupling of CSP with solid-state NMR (SSNMR) greatly enhances the performance and the accuracy of the predictive method, leading to the so-called CSP-NMR crystallography (CSP-NMRX). In this paper, we present the successful application of CSP-NMRX to determine the crystal structure of three structural isomers of pyridine dicarboxylic acid, namely quinolinic, dipicolinic and dinicotinic acids, which can be in a zwitterionic form, or not, in the solid state. In a first step, mono- and bidimensional SSNMR spectra, i.e., ¹H Magic-Angle Spinning (MAS), ¹³C and ¹⁵N Cross Polarisation Magic-Angle Spinning (CPMAS), ¹H Double Quantum (DQ) MAS, ¹H-¹³C HETeronuclear CORrelation (HETCOR), were used to determine the correct molecular structure (i.e., zwitterionic or not) and the local molecular arrangement; at the end, the RMSEs between experimental and computed ¹H and ¹³C chemical shifts allowed the selection of the correct predicted structure for each system. Interestingly, while quinolinic and dipicolinic acids are zwitterionic and non-zwitterionic, respectively, in the solid state, dinicotinic acid exhibits in its crystal structure a "zwitterionic-non-zwitterionic continuum state" in which the proton is shared between the carboxylic moiety and the pyridinic nitrogen. Very refined SSNMR experiments were carried out, i.e., ¹⁴N-¹H Phase-Modulated (PM) pulse and Rotational-Echo Saturation-Pulse Double-Resonance (RESPDOR), to provide an accurate N-H distance value confirming the hybrid nature of the molecule. The CSP-NMRX method showed a remarkable match between the selected structures and the experimental ones. The correct molecular input provided by SSNMR reduced the number of CSP calculations to be performed, leading to different predicted structures, while RMSEs provided an independent parameter with respect to the computed energy for the selection of the best candidate.

Fast magic angle spinning for the characterization of milligram quantities of organic and biological solids at natural isotopic abundance by 13C-13C correlation DNP-enhanced NMR

We show that multidimensional solid-state NMR ¹³C-¹³C correlation spectra of biomolecular assemblies and microcrystalline organic molecules can be acquired at natural isotopic abundance with only milligram quantities of sample. These experiments combine fast Magic Angle Spinning of the sample, low-power dipolar recoupling, and dynamic nuclear polarization performed with AsymPol biradicals, a recently introduced family of polarizing agents. Such experiments are essential for structural characterization as they provide short- and long-range distance information. This approach is demonstrated on diverse sample types, including polyglutamine fibrils implicated in Huntington's disease and microcrystalline ampicillin, a small antibiotic molecule.

Exploring the Potential of Multinuclear Solid-State 1 H, 13 C, and 35 Cl Magnetic Resonance To Characterize Static and Dynamic Disorder in Pharmaceutical Hydrochlorides

Crystallographic disorder, whether static or dynamic, can be detrimental to the physical and chemical stability, ease of crystallization and dissolution rate of an active pharmaceutical ingredient. Disorder can result in a loss of manufacturing control leading to batch-to-batch variability and can lengthen the process of structural characterization. The range of NMR active nuclei makes solid-state NMR a unique technique for gaining nucleus-specific information about crystallographic disorder. Here, we explore the use of high-field ³⁵ Cl solid-state NMR at 23.5 T to characterize both static and dynamic crystallographic disorder: specifically, dynamic disorder occurring in duloxetine hydrochloride (1), static disorder in promethazine hydrochloride (2), and trifluoperazine dihydrochloride (3). In all structures, the presence of crystallographic disorder was confirmed by ¹³ C cross-polarization magic-angle spinning (CPMAS) NMR and supported by GIPAW-DFT calculations, and in the case of 3, ¹ H solid-state NMR provided additional confirmation. Applying ³⁵ Cl solid-state NMR to these compounds, we show that higher magnetic fields are beneficial for resolving the crystallographic disorder in 1 and 3, while broad spectral features were observed in 2 even at higher fields. Combining the data obtained from ¹ H, ¹³ C, and ³⁵ Cl NMR, we show that 3 exhibits a unique case of disorder involving the ⁺ N-H hydrogen positions of the piperazinium ring, driving the chloride anions to occupy three distinct sites.

Narrowing down the conformational space with solid-state NMR in crystal structure prediction of linezolid cocrystals

Many solids crystallize as microcrystalline powders, thus precluding the application of single crystal X-Ray diffraction in structural elucidation. In such cases, a joint use of high-resolution solid-state NMR and crystal structure prediction (CSP) calculations can be successful. However, for molecules showing significant conformational freedom, the CSP-NMR protocol can meet serious obstacles, including ambiguities in NMR signal assignment and too wide conformational search space to be covered by computational methods in reasonable time. Here, we demonstrate a possible way of avoiding these obstacles and making as much use of the two methods as possible in difficult circumstances. In a simple case, our experiments led to crystal structure elucidation of a cocrystal of linezolid (LIN), a wide-range antibiotic, with 2,3-dihydroxybenzoic acid, while a significantly more challenging case of a cocrystal of LIN with 2,4-dihydroxybenzoic acid led to the identification of the most probable conformations of LIN inside the crystal. Having four rotatable bonds, some of which can assume many discreet values, LIN molecule poses a challenge in establishing its conformation in a solid phase. In our work, a set of 27 conformations were used in CSP calculations to yield model crystal structures to be examined against experimental solid-state NMR data, leading to a reliable identification of the most probable molecular arrangements.

Combining X-ray and NMR Crystallography to Explore the Crystallographic Disorder in Salbutamol Oxalate

Salbutamol is an active pharmaceutical ingredient commonly used to treat respiratory distress and is listed by the World Health Organization as an essential medicine. Here, we establish the crystal structure of its oxalate form, salbutamol oxalate, and explore the nature of its crystallographic disorder by combined X-ray crystallography and ¹³C cross-polarization (CP) magic-angle spinning (MAS) solid-state NMR. The *C-OH chiral center of salbutamol (note that the crystal structures are a racemic mixture of the two enantiomers of salbutamol) is disordered over two positions, and the tert-butyl group is rotating rapidly, as revealed by ¹³C solid-state NMR. The impact of crystallization conditions on the disorder was investigated, finding variations in the occupancy ratio of the *C-OH chiral center between single crystals and a consistency across samples in the bulk powder. Overall, this work highlights the contrast between investigating crystallographic disorder by X-ray diffraction and solid-state NMR experiment, and gauge-including projector-augmented-wave (GIPAW) density functional theory (DFT) calculations, with their combined use, yielding an improved understanding of the nature of the crystallographic disorder between the local (i.e., as viewed by NMR) and longer-range periodic (i.e., as viewed by diffraction) scale.

Molecular and Electronic Structure of Isolated Platinum Sites Enabled by the Expedient Measurement of 195Pt Chemical Shift Anisotropy

Techniques that can characterize the molecular structures of dilute surface species are required to facilitate the rational synthesis and improvement of Pt-based heterogeneous catalysts. ¹⁹⁵Pt solid-state NMR spectroscopy could be an ideal tool for this task because ¹⁹⁵Pt isotropic chemical shifts and chemical shift anisotropy (CSA) are highly sensitive probes of the local chemical environment and electronic structure. However, the characterization of Pt surface-sites is complicated by the typical low Pt loadings that are between 0.2 and 5 wt% and broadening of ¹⁹⁵Pt solid-state NMR spectra by CSA. Here, we introduce a set of solid-state NMR methods that exploit fast MAS and indirect detection using a sensitive spy nucleus (¹H or ³¹P) to enable the rapid acquisition of ¹⁹⁵Pt MAS NMR spectra. We demonstrate that high-resolution wideline ¹⁹⁵Pt MAS NMR spectra can be acquired in minutes to a few hours for a series of molecular and single-site Pt species grafted on silica with Pt loading of only 3-5 wt%. Low-power, long-duration, sideband-selective excitation, and saturation pulses are incorporated into t₁-noise eliminated dipolar heteronuclear multiple quantum coherence, perfect echo resonance echo saturation pulse double resonance, or J-resolved pulse sequences. The complete ¹⁹⁵Pt MAS NMR spectrum is then reconstructed by recording a series of 1D NMR spectra where the offset of the ¹⁹⁵Pt pulses is varied in increments of the MAS frequency. Analysis of the ¹⁹⁵Pt MAS NMR spectra yields the ¹⁹⁵Pt chemical shift tensor parameters. Zeroth order approximation density functional theory calculations accurately predict ¹⁹⁵Pt CS tensor parameters. Simple and predictive orbital models relate the CS tensor parameters to the Pt electronic structure and coordination environment. The methodology developed here paves the way for the detailed structural and electronic analysis of dilute platinum surface-sites.

Half-integer-spin quadrupolar nuclei in magic-angle spinning paramagnetic NMR: The case of NaMnO2

A combination of solid-state NMR methods for the extraction of ²³Na shift and quadrupolar parameters in the as-synthesized, structurally complex NaMnO₂ Na-ion cathode material, under magic-angle spinning (MAS) is presented. We show that the integration of the Magic-Angle Turning experiment with Rotor-Assisted Population transfer (RAPT) can be used both to identify shifts and to extract a range of magnitudes for their quadrupolar couplings. We also demonstrate the applicability of the two-dimensional one pulse (TOP) based double-sheared Satellite Transition Magic-Angle Spinning (TOP-STMAS) showing how it can yield a spectrum with separated shift and second-order quadrupolar anisotropies, which in turn can be used to analyze a quadrupolar lineshape free of anisotropic bulk magnetic susceptibility (ABMS) induced shift dispersion and determine both isotropic shift and quadrupolar products. Combining all these experiments, the shift and quadrupolar parameters for all observed Na environments were extracted and yielded excellent agreement with the density functional theory (DFT) based models that were reported in previous literature. We expect these methods to open the door for new possibilities for solid-state NMR to probe half-integer quadrupolar nuclei in paramagnetic materials and other systems exhibiting large shift dispersion.

Complete resonance assignment of a pharmaceutical drug at natural isotopic abundance from DNP-Enhanced solid-state NMR

Solid-state dynamic nuclear polarization enhanced magic angle spinning (DNP-MAS) NMR measurements coupled with density functional theory (DFT) calculations enable the full resonance assignment of a complex pharmaceutical drug molecule without the need for isotopic enrichment. DNP dramatically enhances the NMR signals, thereby making possible previously intractable two-dimensional correlation NMR spectra at natural abundance. Using inputs from DFT calculations, herein we describe a significant improvement to the structure elucidation process for complex organic molecules. Further, we demonstrate that a series of two-dimensional correlation experiments, including ¹⁵N-¹³C TEDOR, ¹³C-¹³C INADEQUATE/SARCOSY, ¹⁹F-¹³C HETCOR, and ¹H-¹³C HETCOR, can be obtained at natural isotopic abundance within reasonable experiment times, thus enabling a complete resonance assignment of sitagliptin, a pharmaceutical used for the treatment of type 2 diabetes.

Resolving alternative structure determinations of indapamide using 13C solid-state NMR

The conflict between alternative crystal structures in the Cambridge Structural Database for the diuretic drug indapamide hemihydrate (IND) has been resolved with the aid of ¹³C solid-state NMR. IND is seen to contain multiple distinct molecules in the asymmetric unit (Z' = 4) rather than exhibiting disorder in the orientation of sulfonamide groups. The NMR crystallographic approach is a more effective tool for distinguishing between alternative structures than naïve judgements of quality based on crystallographic refinement agreement factors.

On the influence of water on THz vibrational spectral features of molecular crystals

The nanoscale structure of molecular assemblies plays a major role in many (μ)-biological mechanisms. Molecular crystals are one of the most simple of these assemblies and are widely used in a variety of applications from pharmaceuticals and agrochemicals, to nutraceuticals and cosmetics. The collective vibrations in such molecular crystals can be probed using terahertz spectroscopy, providing unique characteristic spectral fingerprints. However, the association of the spectral features to the crystal conformation, crystal phase and its environment is a difficult task. We present a combined computational-experimental study on the incorporation of water in lactose molecular crystals, and show how simulations can be used to associate spectral features in the THz region to crystal conformations and phases. Using periodic DFT simulations of lactose molecular crystals, the role of water in the observed lactose THz spectrum is clarified, presenting both direct and indirect contributions. A specific experimental setup is built to allow the controlled heating and corresponding dehydration of the sample, providing the monitoring of the crystal phase transformation dynamics. Besides the observation that lactose phases and phase transformation appear to be more complex than previously thought - including several crystal forms in a single phase and a non-negligible water content in the so-called anhydrous phase - we draw two main conclusions from this study. Firstly, THz modes are spread over more than one molecule and require periodic computation rather than a gas-phase one. Secondly, hydration water does not only play a perturbative role but also participates in the facilitation of the THz vibrations.

From Operando Raman Mechanochemistry to "NMR Crystallography": Understanding the Structures and Interconversion of Zn-Terephthalate Networks Using Selective 17O-Labeling

The description of the formation, structure, and reactivity of coordination networks and metal-organic frameworks (MOFs) remains a real challenge in a number of cases. This is notably true for compounds composed of Zn²⁺ ions and terephthalate ligands (benzene-1,4-dicarboxylate, BDC) because of the difficulties in isolating them as pure phases and/or because of the presence of structural defects. Here, using mechanochemistry in combination with operando Raman spectroscopy, the observation of the formation of various zinc terephthalate compounds was rendered possible, allowing the distinction and isolation of three intermediates during the ball-milling synthesis of Zn₃(OH)₄(BDC). An "NMR crystallography" approach was then used, combining solid-state NMR (¹H, ¹³C, and ¹⁷O) and density functional theory (DFT) calculations to refine the poorly described crystallographic structures of these phases. Particularly noteworthy are the high-resolution ¹⁷O NMR analyses, which were made possible in a highly efficient and cost-effective way, thanks to the selective ¹⁷O-enrichment of either hydroxyl or terephthalate groups by ball-milling. This allowed the presence of defect sites to be identified for the first time in one of the phases, and the nature of the H-bonding network of the hydroxyls to be established in another. Lastly, the possibility of using deuterated precursors (e.g., D₂O and d ₄-BDC) during ball-milling is also introduced as a means for observing specific transformations during operando Raman spectroscopy studies, which would not have been possible with hydrogenated equivalents. Overall, the synthetic and spectroscopic approaches developed herein are expected to push forward the understanding of the structure and reactivity of other complex coordination networks and MOFs.

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.

Atomic-resolution chemical characterization of (2x)72-kDa tryptophan synthase via four- and five-dimensional 1H-detected solid-state NMR

NMR chemical shifts provide detailed information on the chemical properties of molecules, thereby complementing structural data from techniques like X-ray crystallography and electron microscopy. Detailed analysis of protein NMR data, however, often hinges on comprehensive, site-specific assignment of backbone resonances, which becomes a bottleneck for molecular weights beyond 40 to 45 kDa. Here, we show that assignments for the (2x)72-kDa protein tryptophan synthase (665 amino acids per asymmetric unit) can be achieved via higher-dimensional, proton-detected, solid-state NMR using a single, 1-mg, uniformly labeled, microcrystalline sample. This framework grants access to atom-specific characterization of chemical properties and relaxation for the backbone and side chains, including those residues important for the catalytic turnover. Combined with first-principles calculations, the chemical shifts in the β-subunit active site suggest a connection between active-site chemistry, the electrostatic environment, and catalytically important dynamics of the portal to the β-subunit from solution.

Imaging active site chemistry and protonation states: NMR crystallography of the tryptophan synthase α-aminoacrylate intermediate

NMR-assisted crystallography-the integrated application of solid-state NMR, X-ray crystallography, and first-principles computational chemistry-holds significant promise for mechanistic enzymology: by providing atomic-resolution characterization of stable intermediates in enzyme active sites, including hydrogen atom locations and tautomeric equilibria, NMR crystallography offers insight into both structure and chemical dynamics. Here, this integrated approach is used to characterize the tryptophan synthase α-aminoacrylate intermediate, a defining species for pyridoxal-5'-phosphate-dependent enzymes that catalyze β-elimination and replacement reactions. For this intermediate, NMR-assisted crystallography is able to identify the protonation states of the ionizable sites on the cofactor, substrate, and catalytic side chains as well as the location and orientation of crystallographic waters within the active site. Most notable is the water molecule immediately adjacent to the substrate β-carbon, which serves as a hydrogen bond donor to the ε-amino group of the acid-base catalytic residue βLys87. From this analysis, a detailed three-dimensional picture of structure and reactivity emerges, highlighting the fate of the L-serine hydroxyl leaving group and the reaction pathway back to the preceding transition state. Reaction of the α-aminoacrylate intermediate with benzimidazole, an isostere of the natural substrate indole, shows benzimidazole bound in the active site and poised for, but unable to initiate, the subsequent bond formation step. When modeled into the benzimidazole position, indole is positioned with C3 in contact with the α-aminoacrylate Cβ and aligned for nucleophilic attack. Here, the chemically detailed, three-dimensional structure from NMR-assisted crystallography is key to understanding why benzimidazole does not react, while indole does.

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...

Multi-Component Crystals Containing Urea: Mechanochemical Synthesis and Characterization by 35Cl Solid-State NMR Spectroscopy and DFT Calculations

Mechanochemical synthesis provides new pathways for the rational design of multi-component crystals (MCCs) involving anionic or cationic components, which offer molecular-level architectures unavailable to MCCs comprised of strictly neutral components....

A toolbox for improving the workflow of NMR crystallography

NMR crystallography is a powerful tool with applications in structural characterization and crystal structure verification, to name two. However, applying this tool presents several challenges, especially for industrial users, in terms of consistency, workflow, time consumption, and the requirement for a high level of understanding of experimental solid-state NMR and GIPAW-DFT calculations. Here, we have developed a series of fully parameterized scripts for use in Materials Studio and TopSpin, based on the .magres file format, with a focus on organic molecules (e.g. pharmaceuticals), improving efficiency, robustness, and workflow. We separate these tools into three major categories: performing the DFT calculations, extracting & visualizing the results, and crystallographic modelling. These scripts will rapidly submit fully parameterized CASTEP jobs, extract data from the calculations, assist in visualizing the results, and expedite the process of structural modelling. Accompanied with these tools is a description on their functionality, documentation on how to get started and use the scripts, and links to video tutorials for guiding new users. Through the use of these tools, we hope to facilitate NMR crystallography and to harmonize the process across users.

Field-stepped ultra-wideline NMR at up to 36 T: On the inequivalence between field and frequency stepping

Field-stepped NMR spectroscopy at up to 36 T using the series-connected hybrid (SCH) magnet at the U.S. National High Magnetic Field Laboratory is demonstrated for acquiring ultra-wideline powder spectra of nuclei with very large quadrupolar interactions. Historically, NMR evolved from the continuous-wave (cw) field-swept method in the early days to the pulsed Fourier-transform method in the modern era. Spectra acquired using field sweeping are generally considered to be equivalent to those acquired using the pulsed method. Here, it is shown that field-stepped wideline spectra of half-integer spin quadrupolar nuclei acquired using WURST/CPMG methods can be significantly different from those acquired with the frequency-stepped method commonly used with superconducting magnets. The inequivalence arises from magnetic field-dependent NMR interactions such as the anisotropic chemical shift and second-order quadrupolar interactions; the latter is often the main interaction leading to ultra-wideline powder patterns of half-integer spin quadrupolar nuclei. This inequivalence needs be taken into account to accurately and correctly determine the quadrupolar coupling and chemical shift parameters. A simulation protocol is developed for spectral fitting to facilitate analysis of field-stepped ultra-wideline NMR spectra acquired using powered magnets. A MATLAB program which implements this protocol is available on request.

Structure Solution of Nano-Crystalline Small Molecules Using MicroED and Solid-State NMR Dipolar-Based Experiments

Three-dimensional electron diffraction crystallography (microED) can solve structures of sub-micrometer crystals, which are too small for single crystal X-ray crystallography. However, R factors for the microED-based structures are generally high because of dynamic scattering. That means R factor may not be reliable provided that kinetic analysis is used. Consequently, there remains ambiguity to locate hydrogens and to assign nuclei with close atomic numbers, like carbon, nitrogen, and oxygen. Herein, we employed microED and ssNMR dipolar-based experiments together with spin dynamics numerical simulations. The NMR dipolar-based experiments were ¹H-¹⁴N phase-modulated rotational-echo saturation-pulse double-resonance (PM-S-RESPDOR) and ¹H-¹H selective recoupling of proton (SERP) experiments. The former examined the dephasing effect of a specific ¹H resonance under multiple ¹H-¹⁴N dipolar couplings. The latter examined the selective polarization transfer between a ¹H-¹H pair. The structure was solved by microED and then validated by evaluating the agreement between experimental and calculated dipolar-based NMR results. As the measurements were performed on ¹H and ¹⁴N, the method can be employed for natural abundance samples. Furthermore, the whole validation procedure was conducted at 293 K unlike widely used chemical shift calculation at 0 K using the GIPAW method. This combined method was demonstrated on monoclinic l-histidine.

Mechanochemistry: A Green Approach in the Preparation of Pharmaceutical Cocrystals

Mechanochemistry is considered an alternative attractive greener approach to prepare diverse molecular compounds and has become an important synthetic tool in different fields (e.g., physics, chemistry, and material science) since is considered an ecofriendly procedure that can be carried out under solvent free conditions or in the presence of minimal quantities of solvent (catalytic amounts). Being able to substitute, in many cases, classical solution reactions often requiring significant amounts of solvents. These sustainable methods have had an enormous impact on a great variety of chemistry fields, including catalysis, organic synthesis, metal complexes formation, preparation of multicomponent pharmaceutical solid forms, etc. In this sense, we are interested in highlighting the advantages of mechanochemical methods on the obtaining of pharmaceutical cocrystals. Hence, in this review, we describe and discuss the relevance of mechanochemical procedures in the formation of multicomponent solid forms focusing on pharmaceutical cocrystals. Additionally, at the end of this paper, we collect a chronological survey of the most representative scientific papers reporting the mechanochemical synthesis of cocrystals.

Combining X-rays, neutrons and electrons, and NMR, for precision and accuracy in structure-function studies

The distinctive features of the physics-based probes used in understanding the structure of matter focusing on biological sciences, but not exclusively, are described in the modern context. This is set in a wider scope of holistic biology and the scepticism about `reductionism', what is called the `molecular level', and how to respond constructively. These topics will be set alongside the principles of accuracy and precision, and their boundaries. The combination of probes and their application together is the usual way of realizing accuracy. The distinction between precision and accuracy can be blurred by the predictive force of a precise structure, thereby lending confidence in its potential accuracy. These descriptions will be applied to the comparison of cryo and room-temperature protein crystal structures as well as the solid state of a crystal and the same molecules studied by small-angle X-ray scattering in solution and by electron microscopy on a sample grid. Examples will include: time-resolved X-ray Laue crystallography of an enzyme Michaelis complex formed directly in a crystal equivalent to in vivo; a new iodoplatin for radiation therapy predicted from studies of platin crystal structures; and the field of colouration of carotenoids, as an effective assay of function, i.e. their colouration, when unbound and bound to a protein. The complementarity of probes, as well as their combinatory use, is then at the foundation of real (biologically relevant), probe-artefacts-free, structure-function studies. The foundations of our methodologies are being transformed by colossal improvements in technologies of X-ray and neutron sources and their beamline instruments, as well as improved electron microscopes and NMR spectrometers. The success of protein structure prediction from gene sequence recently reported by CASP14 also opens new doors to change and extend the foundations of the structural sciences.

Dynamics & Spectroscopy with Neutrons-Recent Developments & Emerging Opportunities

This work provides an up-to-date overview of recent developments in neutron spectroscopic techniques and associated computational tools to interrogate the structural properties and dynamical behavior of complex and disordered materials, with a focus on those of a soft and polymeric nature. These have and continue to pave the way for new scientific opportunities simply thought unthinkable not so long ago, and have particularly benefited from advances in high-resolution, broadband techniques spanning energy transfers from the meV to the eV. Topical areas include the identification and robust assignment of low-energy modes underpinning functionality in soft solids and supramolecular frameworks, or the quantification in the laboratory of hitherto unexplored nuclear quantum effects dictating thermodynamic properties. In addition to novel classes of materials, we also discuss recent discoveries around water and its phase diagram, which continue to surprise us. All throughout, emphasis is placed on linking these ongoing and exciting experimental and computational developments to specific scientific questions in the context of the discovery of new materials for sustainable technologies.

Natural abundance oxygen-17 solid-state NMR of metal organic frameworks enhanced by dynamic nuclear polarization

The 17O resonances of zirconium-oxo clusters that can be found in porous Zr carboxylate metal-organic frameworks (MOFs) have been investigated by magic-angle spinning (MAS) NMR spectroscopy enhanced by dynamic nuclear polarization (DNP). High-resolution 17O spectra at 0.037% natural abundance could be obtained in 48 hours, thanks to DNP enhancement of the 1H polarization by factors ε(1H) = Swith/Swithout = 28, followed by 1H → 17O cross-polarization, allowing a saving in experimental time by a factor of ca. 800. The distinct 17O sites from the oxo-clusters can be resolved at 18.8 T. Their assignment is supported by density functional theory (DFT) calculations of chemical shifts and quadrupolar parameters. Protonation of 17O sites seems to be leading to large characteristic shifts. Hence, natural abundance 17O NMR spectra of diamagnetic MOFs can thus be used to probe and characterize the local environment of different 17O sites on an atomic scale.

Resolving the Acid Site Distribution in Zn-Exchanged ZSM-5 with Stimulated Raman Scattering Microscopy

Zeolites are widely used acid catalysts in research and in industrial processes. The catalytic performance of these materials is affected by the nature and concentration of Brønsted and Lewis acid sites. The balance between these types of active sites—and thus the activity and selectivity of the zeolite—can be altered by the introduction of metal species, e.g., by ion exchange. Although the acidic properties of zeolites are routinely characterized by bulk-scale techniques, this ensemble-averaged approach neglects the local variations in the material. Insights into the distribution of active sites at the single-particle level are thus critical to better understand the impact of post-synthetic modifications on the zeolite acidity. In this contribution, we spatially resolve Brønsted and Lewis acid sites in protonated and Zn-exchanged ZSM-5 crystals. To this end, the vibrational modes of pyridine chemisorbed on active sites are mapped with stimulated Raman scattering (SRS) microscopy. The SRS images reveal sharp inter- and intra-particle heterogeneities in the distribution of Lewis acid sites introduced upon ion exchange, ascribed to local variations in the Al content. Besides assessing the impact of Zn exchange on the active site distribution in ZSM-5 crystals, this approach enables uniquely to map the distribution of Lewis acid sites in catalysts at the single-particle level.

Bridging photochemistry and photomechanics with NMR crystallography: the molecular basis for the macroscopic expansion of an anthracene ester nanorod

Crystals composed of photoreactive molecules represent a new class of photomechanical materials with the potential to generate large forces on fast timescales. An example is the photodimerization of 9-tert-butyl-anthracene ester (9TBAE) in molecular crystal nanorods that leads to an average elongation of 8%. Previous work showed that this expansion results from the formation of a metastable crystalline product. In this article, it is shown how a novel combination of ensemble oriented-crystal solid-state NMR, X-ray diffraction, and first principles computational modeling can be used to establish the absolute unit cell orientations relative to the shape change, revealing the atomic-resolution mechanism for the photomechanical response and enabling the construction of a model that predicts an elongation of 7.4%, in good agreement with the experimental value. According to this model, the nanorod expansion does not result from an overall change in the volume of the unit cell, but rather from an anisotropic rearrangement of the molecular contents. The ability to understand quantitatively how molecular-level photochemistry generates mechanical displacements allows us to predict that the expansion could be tuned from +9% to -9.5% by controlling the initial orientation of the unit cell with respect to the nanorod axis. This application of NMR-assisted crystallography provides a new tool capable of tying the atomic-level structural rearrangement of the reacting molecular species to the mechanical response of a nanostructured sample.

Pharmaceutical Hydrates Analysis-Overview of Methods and Recent Advances

This review discusses a set of instrumental and computational methods that are used to characterize hydrated forms of APIs (active pharmaceutical ingredients). The focus has been put on highlighting advantages as well as on presenting some limitations of the selected analytical approaches. This has been performed in order to facilitate the choice of an appropriate method depending on the type of the structural feature that is to be analyzed, that is, degree of hydration, crystal structure and dynamics, and (de)hydration kinetics. The presented techniques include X-ray diffraction (single crystal X-ray diffraction (SCXRD), powder X-ray diffraction (PXRD)), spectroscopic (solid state nuclear magnetic resonance spectroscopy (ssNMR), Fourier-transformed infrared spectroscopy (FT-IR), Raman spectroscopy), thermal (differential scanning calorimetry (DSC), thermogravimetric analysis (TGA)), gravimetric (dynamic vapour sorption (DVS)), and computational (molecular mechanics (MM), Quantum Mechanics (QM), molecular dynamics (MD)) methods. Further, the successful applications of the presented methods in the studies of hydrated APIs as well as studies on the excipients' influence on these processes have been described in many examples.

Chlorine-35 Solid-State Nuclear Magnetic Resonance Spectroscopy as an Indirect Probe of the Oxidation Number of Tin in Tin Chlorides

Ultrawideline ³⁵Cl solid-state nuclear magnetic resonance (SSNMR) spectra of a series of 12 tin chlorides were recorded. The magnitude of the ³⁵Cl quadrupolar coupling constant (C_Q) was shown to consistently indicate the chemical state (oxidation number) of the bound Sn center. The chemical state of the Sn center was independently verified by tin Mössbauer spectroscopy. C_Q(³⁵Cl) values of >30 MHz correspond to Sn(IV), while C_Q(³⁵Cl) readings of <30 MHz indicate that Sn(II) is present. Tin-119 SSNMR experiments would seem to be the most direct and effective route to interrogating tin in these systems, yet we show that ambiguous results can emerge from this method, which may lead to an incorrect interpretation of the Sn oxidation number. The accumulated ³⁵Cl NMR data are used as a guide to assign the Sn oxidation number in the mixed-valent metal complex Ph₃PPd^ImSnCl₂. The synthesis and crystal structure of the related Ph₃PPt^ImSnCl₂ are reported, and ¹⁹⁵Pt and ³⁵Cl SSNMR experiments were also used to investigate its Pt-Sn bonding. Plane-wave DFT calculations of ³⁵Cl, ¹¹⁹Sn, and ¹⁹⁵Pt NMR parameters are used to model and interpret experimental data, supported by computed ¹¹⁹Sn and ¹⁹⁵Pt chemical shift tensor orientations. Given the ubiquity of directly bound Cl centers in organometallic and inorganic systems, there is tremendous potential for widespread usage of ³⁵Cl SSNMR parameters to provide a reliable indication of the chemical state in metal chlorides.

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.

Anisotropy and NMR spectroscopy

Abstract

In this paper, different aspects concerning anisotropy in Nuclear Magnetic Resonance (NMR) spectroscopy have been reviewed. In particular, the relevant theory has been presented, showing how anisotropy stems from the dependence of internal nuclear spin interactions on the molecular orientation with respect to the external magnetic field direction. The consequences of anisotropy in the use of NMR spectroscopy have been critically discussed: on one side, the availability of very detailed structural and dynamic information, and on the other side, the loss of spectral resolution. The experiments used to measure the anisotropic properties in solid and soft materials, where, in contrast to liquids, such properties are not averaged out by the molecular tumbling, have been described. Such experiments can be based either on static low-resolution techniques or on one- and two-dimensional pulse sequences exploiting Magic Angle Spinning (MAS). Examples of applications of NMR spectroscopy have been shown, which exploit anisotropy to obtain important physico-chemical information on several categories of systems, including pharmaceuticals, inorganic materials, polymers, liquid crystals, and self-assembling amphiphiles in water. Solid-state NMR spectroscopy can be considered, nowadays, one of the most powerful characterization techniques for all kinds of solid, either amorphous or crystalline, and semi-solid systems for the obtainment of both structural and dynamic properties on a molecular and supra-molecular scale.

Graphic abstract

Application of NMR Crystallography to Highly Disordered Templated Materials: Extensive Local Structural Disorder in the Gallophosphate GaPO-34A

We present an NMR crystallographic investigation of two as-made forms of the recently characterized gallophosphate GaPO-34A, which has an unusual framework composition with a Ga:P ratio of 7:6 and contains both hydroxide and fluoride anions and either 1-methylimidazolium or pyridinium as the structure-directing agent. We combine previously reported X-ray crystallographic data with solid-state NMR spectroscopy and periodic density functional theory (DFT) calculations to show that the structure contains at least three distinct types of disorder (occupational, compositional, and dynamic). The occupational disorder arises from the presence of six anion sites per unit cell, but a total occupancy of five of these, leading to full occupancy of four sites and partial occupancy of the fifth and sixth (which are related by symmetry). The mixture of OH and F present leads to compositional disorder on the occupied anion sites, although the occupancy of some sites by F is calculated to be energetically unfavorable and signals relating to F on these sites are not observed by NMR spectroscopy, confirming that the compositional disorder is not random. Finally, a combination of high-field ⁷¹Ga NMR spectroscopy and variable-temperature ¹³C and ³¹P NMR experiments shows that the structure directing agents are dynamic on the microsecond time scale, which can be supported by averaging the ³¹P chemical shifts calculated with the SDA in different orientations. This demonstrates the value of an NMR crystallographic approach, particularly in the case of highly disordered crystalline materials, where the growth of large single crystals for conventional structure determination may not be possible owing to the extent of disorder present.

Investigating discrepancies between experimental solid-state NMR and GIPAW calculation: NC-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 ¹H and ¹³C cross-polarisation (CP) magic-angle spinning (MAS) nuclear magnetic resonance (NMR) and two-dimensional ¹⁴N-¹H heteronuclear multiple-quantum coherence MAS NMR spectra are compared with gauge-including projector augmented wave (GIPAW) calculations of the ¹H and ¹³C chemical shifts and the ¹⁴N shifts that additionally depend on the quadrupolar interaction. Considering the high ppm (>10 ​ppm) ¹H 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 ¹³C 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 ¹³C 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.

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.

Multimodal Structure Solution with 19F NMR Crystallography of Spin Singlet Molybdenum Oxyfluorides

Complex crystal structures with subtle atomic-scale details are now routinely solved using complementary tools such as X-ray and/or neutron scattering combined with electron diffraction and imaging. Identifying unambiguous atomic models for oxyfluorides, needed for materials design and structure-property control, is often still a considerable challenge despite their advantageous optical responses and applications in energy storage systems. In this work, NMR crystallography and single-crystal X-ray diffraction are combined for the complete structure solution of three new compounds featuring a rare triangular early transition metal oxyfluoride cluster, [Mo₃O₄F₉]^5-. After framework identification by single-crystal X-ray diffraction, 1D and 2D solid-state ¹⁹F NMR spectroscopy supported by ab initio calculations are used to solve the structures of K₅[Mo₃O₄F₉]·3H₂O (1), K₅[Mo₃O₄F₉]·2H₂O (2), and K₁₆[Mo₃O₄F₉]₂[TiF₆]₃·2H₂O (3) and to assign the nine distinct fluorine sites in the oxyfluoride clusters. Furthermore, ¹⁹F NMR identifies selective fluorine dynamics in K₁₆[Mo₃O₄F₉]₂[TiF₆]₃·2H₂O. These dual scattering and spectroscopy methods are used to demonstrate the generality and sensitivity of ¹⁹F shielding to small changes in bond length, on the order of 0.01 Å or less, even in the presence of hydrogen bonding, metal-metal bonding, and electrostatic interactions. Starting from the structure models, the nature of chemical bonding in the molybdates is explained by molecular orbital theory and electronic structure calculations. The average Mo-Mo distance of 2.505 Å and diamagnetism in 1, 2, and 3 are attributed to a metal-metal bond order of unity along with a 1a²1e⁴ electronic ground state configuration for the [Mo₃O₄F₉]^5- cluster, leading to a rare trimeric spin singlet involving d² Mo⁴⁺ ions. The approach to structure solution and bonding analysis is a powerful strategy for understanding the structures and chemical properties of complex fluorides and oxyfluorides.