Chemical shift computations on a crystallographic basis: some reflections and comments

NMR crystallography of amino acids

The development of NMR crystallography methods requires a reliable database of chemical shifts measured for systems with known crystal structure. We measured and assigned carbon and hydrogen chemical shifts of twenty solid natural amino acids of known polymorphic structure, meticulously determined using powder X-ray diffraction. We then correlated the experimental data with DFT-calculated isotropic shieldings. The small size of the unit cell of most amino acids allowed for advanced computations using various families of DFT functionals, including generalized gradient approximation (GGA), meta-GGA and hybrid DFT functionals. We tested several combinations of functionals for geometry optimizations and NMR calculations. For carbon shieldings, the widely used GGA functional PBE performed very well, although an improvement could be achieved by adding shielding corrections calculated for isolated molecules using a hybrid functional. For hydrogen nuclei, we observed the best performance for NMR calculations carried out with structures optimized at the hybrid DFT level. The high fidelity of the calculations made it possible to assign additional signals that could not be assigned based on experiments alone, for example signals of two non-equivalent molecules in the unit cell of some of the amino acids.

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.

Exploiting solid-state dynamic nuclear polarization NMR spectroscopy to establish the spatial distribution of polymorphic phases in a solid material

Solid-state DNP NMR can enhance the ability to detect minor amounts of solid phases within heterogenous materials. Here we demonstrate that NMR contrast based on the transport of DNP-enhanced polarization can be exploited in the challenging case of early detection of a small amount of a minor polymorphic phase within a major polymorph, and we show that this approach can yield quantitative information on the spatial distribution of the two polymorphs. We focus on the detection of a minor amount (<4%) of polymorph III of m-aminobenzoic acid within a powder sample of polymorph I at natural isotopic abundance. Based on proposed models of the spatial distribution of the two polymorphs, simulations of 1H spin diffusion allow NMR data to be calculated for each model as a function of particle size and the relative amounts of the polymorphs. A comparison between simulated and experimental NMR data allows the model(s) best representing the spatial distribution of the polymorphs in the system to be established.

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.

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.

Monitoring the influence of additives on the crystallization processes of glycine with dynamic nuclear polarization solid-state NMR

Crystallization is fundamental in many domains, and the investigation of the sequence of solid phases produced as a function of crystallization time is thus key to understand and control crystallization processes. Here, we used a solid-state nuclear magnetic resonance strategy to monitor the crystallization process of glycine, which is a model compound in polymorphism, under the influence of crystallizing additives, such as methanol or sodium chloride. More specifically, our strategy is based on a combination of low-temperatures and dynamic nuclear polarization (DNP) to trap and detect transient crystallizing forms, which may be present only in low quantities. Interestingly, our results show that these additives yield valuable DNP signal enhancements even in the absence of glycerol within the crystallizing solution.

A Machine Learning Model of Chemical Shifts for Chemically and Structurally Diverse Molecular Solids

Nuclear magnetic resonance (NMR) chemical shifts are a direct probe of local atomic environments and can be used to determine the structure of solid materials. However, the substantial computational cost required to predict accurate chemical shifts is a key bottleneck for NMR crystallography. We recently introduced ShiftML, a machine-learning model of chemical shifts in molecular solids, trained on minimum-energy geometries of materials composed of C, H, N, O, and S that provides rapid chemical shift predictions with density functional theory (DFT) accuracy. Here, we extend the capabilities of ShiftML to predict chemical shifts for both finite temperature structures and more chemically diverse compounds, while retaining the same speed and accuracy. For a benchmark set of 13 molecular solids, we find a root-mean-squared error of 0.47 ppm with respect to experiment for ¹H shift predictions (compared to 0.35 ppm for explicit DFT calculations), while reducing the computational cost by over four orders of magnitude.

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.

A structure determination protocol based on combined analysis of 3D-ED data, powder XRD data, solid-state NMR data and DFT-D calculations reveals the structure of a new polymorph of l-tyrosine

We report the crystal structure of a new polymorph of l-tyrosine (denoted the β polymorph), prepared by crystallization from the gas phase following vacuum sublimation. Structure determination was carried out by combined analysis of three-dimensional electron diffraction (3D-ED) data and powder X-ray diffraction (XRD) data. Specifically, 3D-ED data were required for reliable unit cell determination and space group assignment, with structure solution carried out independently from both 3D-ED data and powder XRD data, using the direct-space strategy for structure solution implemented using a genetic algorithm. Structure refinement was carried out both from powder XRD data, using the Rietveld profile refinement technique, and from 3D-ED data. The final refined structure was validated both by periodic DFT-D calculations, which confirm that the structure corresponds to an energy minimum on the energy landscape, and by the fact that the values of isotropic ¹³C NMR chemical shifts calculated for the crystal structure using DFT-D methodology are in good agreement with the experimental high-resolution solid-state ¹³C NMR spectrum. Based on DFT-D calculations using the PBE0-MBD method, the β polymorph is meta-stable with respect to the previously reported crystal structure of l-tyrosine (now denoted the α polymorph). Crystal structure prediction calculations using the AIRSS approach suggest that there are three other plausible crystalline polymorphs of l-tyrosine, with higher energy than the α and β polymorphs.

Resolving Atomic-Scale Interactions in Nonfullerene Acceptor Organic Solar Cells with Solid-State NMR Spectroscopy, Crystallographic Modelling, and Molecular Dynamics Simulations

Fused-ring core nonfullerene acceptors (NFAs), designated "Y-series," have enabled high-performance organic solar cells (OSCs) achieving over 18% power conversion efficiency (PCE). Since the introduction of these NFAs, much effort has been expended to understand the reasons for their exceptional performance. While several studies have identified key optoelectronic properties that govern high PCEs, little is known about the molecular level origins of large variations in performance, spanning from 5% to 18% PCE, for example, in the case of PM6:Y6 OSCs. Here, a combined solid-state NMR, crystallography, and molecular modeling approach to elucidate the atomic-scale interactions in Y6 crystals, thin films, and PM6:Y6 bulk heterojunction (BHJ) blends is introduced. It is shown that the Y6 morphologies in BHJ blends are not governed by the morphology in neat films or single crystals. Notably, PM6:Y6 blends processed from different solvents self-assemble into different structures and morphologies, whereby the relative orientations of the sidechains and end groups of the Y6 molecules to their fused-ring cores play a crucial role in determining the resulting morphology and overall performance of the solar cells. The molecular-level understanding of BHJs enabled by this approach will guide the engineering of next-generation NFAs for stable and efficient OSCs.

Solid-State Structural Properties of Alloxazine Determined from Powder XRD Data in Conjunction with DFT-D Calculations and Solid-State NMR Spectroscopy: Unraveling the Tautomeric Identity and Pathways for Tautomeric Interconversion

We report the solid-state structural properties of alloxazine, a tricyclic ring system found in many biologically important molecules, with structure determination carried out directly from powder X-ray diffraction (XRD) data. As the crystal structures containing the alloxazine and isoalloxazine tautomers both give a high-quality fit to the powder XRD data in Rietveld refinement, other techniques are required to establish the tautomeric form in the solid state. In particular, high-resolution solid-state ¹⁵N NMR data support the presence of the alloxazine tautomer, based on comparison between isotropic chemical shifts in the experimental ¹⁵N NMR spectrum and the corresponding values calculated for the crystal structures containing the alloxazine and isoalloxazine tautomers. Furthermore, periodic DFT-D calculations at the PBE0-MBD level indicate that the crystal structure containing the alloxazine tautomer has significantly lower energy. We also report computational investigations of the interconversion between the tautomeric forms in the crystal structure via proton transfer along two intermolecular N-H···N hydrogen bonds; DFT-D calculations at the PBE0-MBD level indicate that the tautomeric interconversion is associated with a lower energy transition state for a mechanism involving concerted (rather than sequential) proton transfer along the two hydrogen bonds. However, based on the relative energies of the crystal structures containing the alloxazine and isoalloxazine tautomers, it is estimated that under conditions of thermal equilibrium at ambient temperature, more than 99.9% of the molecules in the crystal structure will exist as the alloxazine tautomer.

DFT approach to predict 13C NMR chemical shifts of hydrocarbon species adsorbed on Zn-modified zeolite

13C MAS NMR spectroscopy is a powerful technique to study the mechanisms of hydrocarbon transformations on heterogeneous catalysts. It can reliably identify the surface intermediates and the adsorbed products based...

Solid-State NMR-Driven Crystal Structure Prediction of Molecular Crystals: The Case of Mebendazole

Among all possible NMR crystallography approaches for crystal-structure determination, crystal structure prediction - NMR crystallography (CSP-NMRX) has recently turned out to be a powerful method. In the latter, the original procedure exploited solid-state NMR (SSNMR) information during the final steps of the prediction. In particular, it used the comparison of computed and experimental chemical shifts for the selection of the correct crystal packing. Still, the prediction procedure, generally carried out with DFT methods, may require important computational resources and be quite time-consuming, especially if there are no available constraints to use at the initial stage. Herein, the successful application of this combined prediction method, which exploits NMR information also in the input step to reduce the search space of the predictive algorithm, is presented. Herein, this method was applied on mebendazole, which is characterized by desmotropism. The use of SSNMR data as constraints for the selection of the right tautomer and the determination of the number of independent molecules in the unit cell led to a considerably faster process, reducing the number of calculations to be performed. In this way, the crystal packing was successfully predicted for the three known phases of mebendazole. To evaluate the quality of the predicted structures, these were compared to the experimental ones. The crystal structure of phase B of mebendazole, in particular, was determined de novo by powder diffraction and is presented for the first time in this paper.

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.

Bayesian probabilistic assignment of chemical shifts in organic solids

A prerequisite for NMR studies of organic materials is assigning each experimental chemical shift to a set of geometrically equivalent nuclei. Obtaining the assignment experimentally can be challenging and typically requires time-consuming multidimensional correlation experiments. An alternative solution for determining the assignment involves statistical analysis of experimental chemical shift databases, but no such database exists for molecular solids. Here, by combining the Cambridge Structural Database with a machine learning model of chemical shifts, we construct a statistical basis for probabilistic chemical shift assignment of organic crystals by calculating shifts for more than 200,000 compounds, enabling the probabilistic assignment of organic crystals directly from their two-dimensional chemical structure. The approach is demonstrated with the ¹³C and ¹H assignment of 11 molecular solids with experimental shifts and benchmarked on 100 crystals using predicted shifts. The correct assignment was found among the two most probable assignments in more than 80% of cases.

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.

Accurate location of hydrogen atoms in hydrogen bonds of tizoxanide from the combination of experimental and theoretical models

To obtain detailed information about the position of hydrogen atoms in hydrogen bonds, HBs, of crystalline organic molecular compounds is not an easy task. In this work we propose a combination of ssNMR experimental data with theoretical procedures to get such information. Furthermore, the combination of experimental and theoretical models provides us with well-defined grounds to analyse the strength of π-stacking interactions between layers of hydrogen bonded molecules. Two different theoretical models were considered, both approaches being quite different. The first one is a solid-state model, so that the periodicity of a crystalline system underlies calculations of the electronic energy, the electronic density and NMR parameters. The other one is a molecular model in which molecules are taken as isolated monomers, dimers and tetramers. These two models were applied to the tizoxanide, TIZ, molecular crystal though it can widely be applied to any other molecular crystal. By the application of the quantum molecular model it was possible to learn about the way the intermolecular HBs affect the position of hydrogen atoms that belong to HBs in TIZ. This molecule has two intermolecular HBs that stabilize the structure of a basic dimer, but it also has an intramolecular HB in each monomer whose position should be optimized together with the other ones. We found that by doing this it is possible to obtain reliable results of calculations of NMR spectroscopic parameters. Working with the solid-state model we found that any local variation of the TIZ crystalline structure is correlated with the variation of the values of the NMR parameters of each nucleus. The excellent agreement between experimental and calculated chemical shifts leads to the conclusion that the N10-H10 bond distance should be (1.00 ± 0.02) Å.

A curious case of dynamic disorder in pyrrolidine rings elucidated by NMR crystallography

A pharmaceutical exhibits differing dynamics in crystallographically distinct pyrrolidine rings despite being nearly related by symmetry, with one performing ring inversions while the other is constrained to torsional librations. Using ¹³C solid-state magic-angle spinning (MAS) NMR and DFT calculations, we show that this contrast originates from C-HH-C close contacts and less efficient C-Hπ intermolecular interactions observed in the transition state of the constrained pyrrolidine ring, highlighting the influence of the crystallographic environment on the molecular motion.

Structural Study of MgyAl(8+x-2y)/3O4-xNx (0 < x < 0.5, 0 < y < 1) Spinel Probed by X-ray Diffraction, 27Al MAS NMR, and First-Principles Calculations

The deep understanding of the crystal structure and composition-structure relationship is important for modifying and designing solids to obtain functional materials with customized properties. However, because of multiple compositions and complex structures in some spinel solid solutions, the composition-structure relationship is unknown, or becomes very complicated and difficult to be controlled. In this work, the solid-state magic-angle spinning nuclear magnetic resonance (MAS NMR) technique and powder X-ray diffraction (XRD) Rietveld refinement were combined to characterize the crystal structure of quaternary disordered Mg_yAl_(8+x-2y)/3O_4-xN_x solid solutions in detail, which was supported by the first-principles density functional theory calculations. Diffraction data indicate that Mg ions preferentially enter the tetrahedron structure in Mg_yAl_(8+x-2y)/3O_4-xN_x solid solutions because the sum of the bond valence in tetrahedron is closer to the atomic valence of Mg. With the compositional change, the coordination polyhedra in the crystal structure will also adjust the volume according to the changes in lattice parameter, anion parameter, and inversion parameter. In addition, the populations, chemical shifts, and quadrupole coupling constants of Al located in different coordinated environments of the solid solutions were detected through the simulation and integration of ²⁷Al MAS NMR spectra, which were related to the structural parameters by the bond valence method. It turns out that ²⁷Al MAS NMR parameters are highly sensitive to the subtle changes in the local environment of Al caused by the preferential occupation of Mg in the tetrahedron. These results provide deep insights into the crystal structural details of novel spinel materials with multiple disorders.

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

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

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.

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.

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

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

Exploring mechanochemical parameters using a DoE approach: Crystal structure solution from synchrotron XRPD and characterization of a new praziquantel polymorph

A rotated Doehlert matrix was utilized to explore the experimental design space around the milling parameters of Praziquantel (PZQ) polymorph B formation in terms of frequency and milling time. Three experimental responses were evaluated on the resulting ground samples: two quantitative responses, i.e. median particle size by Laser Light scattering (LLS) and drug recovery by HPLC, and one qualitative dependent variable, i.e. the obtained PZQ crystalline form, characterized through X-Ray Powder Diffraction (XRPD) and confirmed by Differential Scanning Calorimetry (DSC) and Thermogravimetric analysis (TGA). Temperature inside the jars was kept under constant control during the milling process by using temperature sensor equipped jars (thermojars), thus allowing evaluation of the obtained solid states at each experimental point, considering the specific temperature of the process. This explorative analysis led to the finding of a novel PZQ polymorph, named "Form C", produced without degradation, then fully characterized, including by means of Synchrotron XRPD, Polarimetric, FT-IR, SS-NMR, ESEM and saturation solubility. Crystal structure was solved from XRPD data and its geometry was optimized by DFT calculations (CASTEP). Finally, Form C and Form A activity against adult schistosoma mansoni were compared through in vitro testing, and Form C's physical stability checked. The new polymorph, crystallizing in space group I2/c, physically stable for approximately 2 months, showed a m.p. of 106.84 °C and displayed excellent biopharmaceutical properties (water solubility of 382.69±9.26 mg/l), while preserving excellent activity levels against adult schistosoma mansoni.

Reducing the computational cost of NMR crystallography of organic powders at natural isotopic abundance with the help of 13 C-13 C dipolar couplings

Structure determination of functional organic compounds remains a formidable challenge when the sample exists as a powder. Nuclear magnetic resonance crystallography approaches based on the comparison of experimental and Density Functional Theory (DFT)-computed ¹ H chemical shifts have already demonstrated great potential for structure determination of organic powders, but limitations still persist. In this study, we discuss the possibility of using ¹³ C-¹³ C dipolar couplings quantified on powdered theophylline at natural isotopic abundance with the help of dynamic nuclear polarization, to realize a DFT-free, rapid screening of a pool of structures predicted by ab initio random structure search. We show that although ¹³ C-¹³ C dipolar couplings can identify structures possessing long range structural motifs and unit cell parameters close to those of the true structure, it must be complemented with other data to recover information about the presence and the chemical nature of the supramolecular interactions.

Improving the accuracy of solid-state nuclear magnetic resonance chemical shift prediction with a simple molecular correction

A fast, straightforward method for computing NMR chemical shieldings of crystalline solids is proposed. The method combines the advantages of both conventional approaches: periodic calculations using plane-wave basis sets and molecular computational approaches. The periodic calculations capture the periodic nature of crystalline solids, but the computational level of the electronic structure calculation is limited to general-gradient-approximation (GGA) density functionals. It is demonstrated that a correction to the GGA result calculated on an isolated molecule at a higher level of theory significantly improves the correlations between experimental and calculated chemical shifts while adding almost no additional computational cost. Corrections calculated with a hybrid density functional improved the accuracy of 13C, 15N and 17O chemical shift predictions significantly and allowed identifying errors in previously published experimental data. Applications of the approach to crystalline isocytosine, methacrylamide, and testosterone are presented.

An XRD and NMR crystallographic investigation of the structure of 2,6-lutidinium hydrogen fumarate

A crystallographic study highlighting the benefits of a combined XRD and NMR approach in investigating both stability and variation within an organic multicomponent crystal.

Ab initio computation for solid-state 31P NMR of inorganic phosphates: revisiting X-ray structures

The complete 31P NMR chemical shift tensors for 22 inorganic phosphates obtained from ab initio computation are found to correspond closely to experimentally obtained parameters. Further improvement was found when structures determined by diffraction were geometry optimized. Besides aiding in spectral assignment, the cases where correspondence is significantly improved upon geometry optimization point to the crystal structures requiring correction.

Chemical shift reference scale for Li solid state NMR derived by first-principles DFT calculations

For studying electrode and electrolyte materials for lithium ion batteries, solid-state (SS) nuclear magnetic resonance (NMR) of lithium moves into focus of current research. Theoretical simulations of magnetic resonance parameters facilitate the analysis and interpretation of experimental Li SS-NMR spectra and provide unique insight into physical and chemical processes that are determining the spectral profile. In the present paper, the accuracy and reliability of the theoretical simulation methods of Li chemical shielding values is benchmarked by establishing a reference scale for Li SS-NMR of diamagnetic compounds. The impact of geometry, ionic mobility and relativity are discussed. Eventually, the simulation methods are applied to the more complex lithium titanate spinel (Li₄Ti₅O₁₂, LTO), which is a widely discussed battery anode material. Simulation of the Li SS-NMR spectrum shows that the commonly adopted approach of assigning the resonances to individual crystallographic sites is not unambiguous.

Theoretical investigations into the variability of the 15N solid-state NMR parameters within an antimicrobial peptide ampullosporin A

The solid-state NMR measurements play an indispensable role in studies of interactions between biological membranes and peptaibols, which are amphipathic oligopeptides with a high abundance of alpha-aminobutyric acid (Aib). The solid-state NMR investigations are important in establishing the molecular models of the pore forming and antimicrobial properties of peptaibols, but rely on certain simplifications. Some of the underlying assumptions concern the parameters describing the 15N NMR chemical shielding tensor (CST) of the amide nitrogens in Aib and in conventional amino acids. Here the density functional theory (DFT) based calculations were applied to the known crystal structure of one of peptaibols, Ampullosporin A, in order to explicitly describe the variation of the 15N NMR parameters within its backbone. Based on the DFT computational data it was possible to verify the validity of the assumptions previously made about the differences between Aib and other amino acids in the isotropic part of the CST. Also the trends in the magnitudes and orientations of the anisotropic components of the CST, as revealed by the DFT calculations of the full periodic structure of Ampullosporin A, were thoroughly analyzed, and may be employed in future studies of peptaibols.

Perspective: Current advances in solid-state NMR spectroscopy

In contrast to the rapid and revolutionary impact of solution-state Nuclear Magnetic Resonance (NMR) on modern chemistry, the field of solid-state NMR has matured more slowly. This reflects the major technical challenges of much reduced spectral resolution and sensitivity in solid-state as compared to solution-state spectra, as well as the relative complexity of the solid state. In this perspective, we outline the technique developments that have pushed resolution to intrinsic limits and the approaches, including ongoing major developments in the field of Dynamic Nuclear Polarisation, that have enhanced spectral sensitivity. The information on local structure and dynamics that can be obtained using these gains in sensitivity and resolution is illustrated with a diverse range of examples from large biomolecules to energy materials and pharmaceuticals and from both ordered and highly disordered materials. We discuss how parallel developments in quantum chemical calculation, particularly density functional theory, have enabled experimental data to be translated directly into information on local structure and dynamics, giving rise to the developing field of "NMR crystallography."

Recent Advances in Solid-State Nuclear Magnetic Resonance Spectroscopy

The sensitivity of nuclear magnetic resonance (NMR) spectroscopy to the local atomic-scale environment offers great potential for the characterization of a diverse range of solid materials. Despite offering more information than its solution-state counterpart, solid-state NMR has not yet achieved a similar level of recognition, owing to the anisotropic interactions that broaden the spectral lines and hinder the extraction of structural information. Here, we describe the methods available to improve the resolution of solid-state NMR spectra and the continuing research in this area. We also highlight areas of exciting new and future development, including recent interest in combining experiment with theoretical calculations, the rise of a range of polarization transfer techniques that provide significant sensitivity enhancements, and the progress of in situ measurements. We demonstrate the detailed information available when studying dynamic and disordered solids and discuss the future applications of solid-state NMR spectroscopy across the chemical sciences.

On the structural aspects of solid solutions of enantiomers: an intriguing case study of enantiomer recognition in the solid state

The structural aspects of type 1 and type 2 solid solutions have been revised.

Full spectroscopic characterization of two crystal pseudopolymorphic forms of the antiandrogen cortexolone 17α-propionate for topic application

Cortexolone-17α-propionate (CP) is a topically active antiandrogen useful in the treatment of skin disorders. In the solid state, three anhydrous forms of this drug (CPI, CPII and CPIII) occur, together with one hydrated crystal (CPW). The single crystal structure of the monohydrated phase, CPW, compared with that of the anhydrous form CPIII, shows a markedly different solid state behavior. These latter pseudopolymorphic forms have also been fully characterized by spectroscopic methods.

New Polymorph Form of Dexamethasone Acetate

A new monohydrated polymorph of dexamethasone acetate was crystallized and its crystal structure characterized. The different analytical techniques used for describing its structural and vibrational properties were: single crystal and polycrystal X-ray diffraction, solid state nuclear magnetic resonance, infrared spectroscopy. A Hirshfeld surface analysis was carried out through self-arrangement cemented by H-bonds observed in this new polymorph. This new polymorph form appeared because of self-arrangement via classical hydrogen bonds around the water molecule.

Ab initio random structure searching of organic molecular solids: assessment and validation against experimental data

This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated.

Probing intermolecular interactions in a diethylcarbamazine citrate salt by fast MAS 1H solid-state NMR spectroscopy and GIPAW calculations

Fast magic-angle spinning (MAS) NMR is used to probe intermolecular interactions in a diethylcarbamazine salt, that is widely used as a treatment against adult worms of Wuchereria bancrofti which cause a common disease in tropical countries named filariasis. Specifically, a dihydrogen citrate salt that has improved thermal stability and solubility as compared to the free form is studied. One-dimensional ¹H, ¹³C and ¹⁵N and two-dimensional ¹H-¹³C and ¹⁴N-¹H heteronuclear correlation NMR experiments under moderate and fast MAS together with GIPAW (CASTEP) calculations enable the assignment of the ¹H, ¹³C and ¹⁴N/¹⁵N resonances. A two-dimensional ¹H-¹H double-quantum (DQ) -single-quantum (SQ) MAS spectrum recorded with BaBa recoupling at 60kHz MAS identifies specific proton-proton proximities associated with citrate-citrate and citrate-diethylcarbamazine intermolecular interactions.

A new polymorph of ciprofloxacin saccharinate: Structural characterization and pharmaceutical profile

In this study, a new polymorph of ciprofloxacin saccharinate has been thoroughly evaluated with respect to structural as well as biopharmaceutical properties. Preliminary characterization of the new polymorph was performed by differential scanning calorimetry and thermogravimetric analysis, later confirmed by Fourier transform infra-red spectroscopy. The crystal structure was determined from the powder X-ray diffraction pattern using the direct-space algorithm. It lies in the triclinic P-1 space group. It is having lattice parameters different from previously reported forms. The solid-state ¹³C NMR data calculated from the crystal structure by exploiting density functional theory were found to be in excellent agreement with corresponding experimental ¹³C NMR data, thus providing a robust validation of the authenticity of the structure. The prepared polymorph shows enhanced aqueous solubility and dissolution rate in contrast to previously reported polymorph. The new form demonstrated improved oral bioavailability and inhibition of bacterial species. The enhanced biopharmaceutical properties are attributed to the higher extent of its absorption and distribution with respect to the individual components. The new form is a potential candidate for the pharmaceutical industry.

Pharmaceutical cocrystals, salts and polymorphs: Advanced characterization techniques

The main goal of a novel drug development is to obtain it with optimal physiochemical, pharmaceutical and biological properties. Pharmaceutical companies and scientists modify active pharmaceutical ingredients (APIs), which often are cocrystals, salts or carefully selected polymorphs, to improve the properties of a parent drug. To find the best form of a drug, various advanced characterization methods should be used. In this review, we have described such analytical methods, dedicated to solid drug forms. Thus, diffraction, spectroscopic, thermal and also pharmaceutical characterization methods are discussed. They all are necessary to study a solid API in its intrinsic complexity from bulk down to the molecular level, gain information on its structure, properties, purity and possible transformations, and make the characterization efficient, comprehensive and complete. Furthermore, these methods can be used to monitor and investigate physical processes, involved in the drug development, in situ and in real time. The main aim of this paper is to gather information on the current advancements in the analytical methods and highlight their pharmaceutical relevance.

Exploring Systematic Discrepancies in DFT Calculations of Chlorine Nuclear Quadrupole Couplings

Previous studies have revealed significant discrepancies between density functional theory (DFT)-calculated and experimental nuclear quadrupolar coupling constants (C_Q) for chlorine atoms, particularly in ionic solids. Various aspects of the computations are systematically investigated here, including the choice of the DFT functional, basis set convergence, and geometry optimization protocol. The effects of fast (fs) time-scale dynamics are probed using molecular dynamics (MD) and nuclear quantum effects (NQEs) are considered using path-integral MD calculations. It is shown that the functional choice is the most important factor related to improving the accuracy of the quadrupolar coupling calculations, and that functionals beyond the generalized gradient approximation (GGA) level, such as hybrid and meta-GGA functionals, are required for good correlations with experiment. The influence of molecular dynamics and NQEs is less important than the functional choice in the studied systems. A method which involves scaling the calculated quadrupolar coupling constant is proposed here; its application leads to good agreement with experimental data.

Determination of a complex crystal structure in the absence of single crystals: analysis of powder X-ray diffraction data, guided by solid-state NMR and periodic DFT calculations, reveals a new 2'-deoxyguanosine structural motif

Derivatives of guanine exhibit diverse supramolecular chemistry, with a variety of distinct hydrogen-bonding motifs reported in the solid state, including ribbons and quartets, which resemble the G-quadruplex found in nucleic acids with sequences rich in guanine. Reflecting this diversity, the solid-state structural properties of 3',5'-bis-O-decanoyl-2'-deoxyguanosine, reported in this paper, reveal a hydrogen-bonded guanine ribbon motif that has not been observed previously for 2'-deoxyguanosine derivatives. In this case, structure determination was carried out directly from powder XRD data, representing one of the most challenging organic molecular structures (a 90-atom molecule) that has been solved to date by this technique. While specific challenges were encountered in the structure determination process, a successful outcome was achieved by augmenting the powder XRD analysis with information derived from solid-state NMR data and with dispersion-corrected periodic DFT calculations for structure optimization. The synergy of experimental and computational methodologies demonstrated in the present work is likely to be an essential feature of strategies to further expand the application of powder XRD as a technique for structure determination of organic molecular materials of even greater complexity in the future.

The application of tailor-made force fields and molecular dynamics for NMR crystallography: a case study of free base cocaine

Motional averaging has been proven to be significant in predicting the chemical shifts in ab initio solid-state NMR calculations, and the applicability of motional averaging with molecular dynamics has been shown to depend on the accuracy of the molecular mechanical force field. The performance of a fully automatically generated tailor-made force field (TMFF) for the dynamic aspects of NMR crystallography is evaluated and compared with existing benchmarks, including static dispersion-corrected density functional theory calculations and the COMPASS force field. The crystal structure of free base cocaine is used as an example. The results reveal that, even though the TMFF outperforms the COMPASS force field for representing the energies and conformations of predicted structures, it does not give significant improvement in the accuracy of NMR calculations. Further studies should direct more attention to anisotropic chemical shifts and development of the method of solid-state NMR calculations.