Semi-empirical refinements of crystal structures using 17O quadrupolar-coupling tensors

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

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

Predicting 51V nuclear magnetic resonance observables in molecular crystals

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

An unusual ionic cocrystal of ponatinib hydrochloride: characterization by single-crystal X-ray diffraction and ultra-high field NMR spectroscopy

This study describes the discovery of a unique ionic cocrystal of the active pharmaceutical ingredient (API) ponatinib hydrochloride (pon·HCl), and characterization using single-crystal X-ray diffraction (SCXRD) and solid-state NMR (SSNMR) spectroscopy. Pon·HCl is a multicomponent crystal that features an unusual stoichiometry, with an asymmetric unit containing both monocations and dications of the ponatinib molecule, three water molecules, and three chloride ions. Structural features include (i) a charged imidazopyridazine moiety that forms a hydrogen bond between the ponatinib monocations and dications and (ii) a chloride ion that does not feature hydrogen bonds involving any organic moiety, instead being situated in a "square" arrangement with three water molecules. Multinuclear SSNMR, featuring high and ultra-high fields up to 35.2 T, provides the groundwork for structural interpretation of complex multicomponent crystals in the absence of diffraction data. A 13C CP/MAS spectrum confirms the presence of two crystallographically distinct ponatinib molecules, whereas 1D 1H and 2D 1H-1H DQ-SQ spectra identify and assign the unusually deshielded imidazopyridazine proton. 1D 35Cl spectra obtained at multiple fields confirm the presence of three distinct chloride ions, with density functional theory calculations providing key relationships between the SSNMR spectra and H⋯Cl- hydrogen bonding arrangements. A 2D 35Cl → 1H D-RINEPT spectrum confirms the spatial proximities between the chloride ions, water molecules, and amine moieties. This all suggests future application of multinuclear SSNMR at high and ultra-high fields to the study of complex API solid forms for which SCXRD data are unavailable, with potential application to heterogeneous mixtures or amorphous solid dispersions.

Resolving the Chemical Formula of Nesquehonite via NMR Crystallography, DFT Computation, and Complementary Neutron Diffraction

Nesquehonite is a magnesium carbonate mineral relevant to carbon sequestration envisioned for carbon capture and storage of CO₂ . Its chemical formula remains controversial today, assigned as either a hydrated magnesium carbonate [MgCO₃  ⋅ 3H₂ O], or a hydroxy- hydrated- magnesium bicarbonate [Mg(HCO₃ )OH ⋅ 2H₂ O]. The resolution of this controversy is central to understanding this material's thermodynamic, phase, and chemical behavior. In an NMR crystallography study, using rotational-echo double-resonance ¹³ C{¹ H} (REDOR), ¹³ C-¹ H distances are determined with precision, and the combination of ¹³ C static NMR lineshapes and density functional theory (DFT) calculations are used to model different H atomic coordinates. [MgCO₃  ⋅ 3H₂ O] is found to be accurate, and evidence from neutron powder diffraction bolsters these assignments. Refined H positions can help understand how H-bonding stabilizes this structure against dehydration to MgCO₃ . More broadly, these results illustrate the power of NMR crystallography as a technique for resolving questions where X-ray diffraction is inconclusive.

Hydrates of active pharmaceutical ingredients: A 35Cl and 2H solid-state NMR and DFT study

This study uses ³⁵Cl and ²H solid-state NMR (SSNMR) spectroscopy and dispersion-corrected plane-wave density functional theory (DFT) calculations to characterize the molecular-level structures and dynamics of hydrates of active pharmaceutical ingredients (APIs). We use ³⁵Cl SSNMR to measure the EFG tensors of the chloride ions to characterize hydrated forms of hydrochloride salts of APIs, along with two corresponding anhydrous forms. DFT calculations are used to refine the crystal structures of the APIs and determine relationships between the ³⁵Cl EFG tensors and the spatial arrangements of proximate hydrogen bonds, which are particularly influenced by interactions with water molecules. We find that the relationship between ³⁵Cl EFG tensors and local hydrogen bonding geometries is complex, but meaningful structure/property relationships can be garnered through use of DFT calculations. Specifically, for every case in which such a comparison could be made, we find that the hydrate has a smaller magnitude of C_Q than the corresponding anhydrous form, indicating a chloride ion environment with a ground-state electron density of higher spherical symmetry in the former. Finally, variable-temperature ³⁵Cl and ²H SSNMR experiments on a deuterium-exchanged sample of the API cimetidine hydrochloride monohydrate are used to monitor temperature-dependent influences on the spectra that may arise from motional influences on the ³⁵Cl and ²H EFG tensors. From the ²H SSNMR spectra, we determine that the motions of water molecules are characterized by jump-like motions about their C₂ rotational axes that occur on timescales that are unlikely to influence the ³⁵Cl central-transition (+1/2 ↔︎ -1/2) powder patterns (this is confirmed by ³⁵Cl SSNMR). Together, these methods show great promise for the future study of APIs in their bulk and dosage forms, especially variable hydrates in which crystallographic water content varies with external conditions such as humidity.

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

Nutraceuticals in Bulk and Dosage Forms: Analysis by 35Cl and 14N Solid-State NMR and DFT Calculations

This study uses ³⁵Cl and ¹⁴N solid-state NMR (SSNMR) spectroscopy and dispersion-corrected plane-wave density functional theory (DFT) calculations for the structural characterization of chloride salts of nutraceuticals in their bulk and dosage forms. For eight nutraceuticals, we measure the ³⁵Cl EFG tensor parameters of the chloride ions and use plane-wave DFT calculations to elucidate relationships between NMR parameters and molecular-level structure, which provide rapid NMR crystallographic assessments of structural features. We employ both ³⁵Cl direct excitation and ¹H→³⁵Cl cross-polarization methods to characterize a dosage form containing α-d-glucosamine HCl, observe possible impurity and/or adulterant phases, and quantify the weight percent of the active ingredient. To complement this, we also investigate ¹⁴N SSNMR spectroscopy and DFT calculations to characterize nitrogen atoms in the nutraceuticals. This includes a discussion of targeted acquisition experimental protocols (i.e., acquiring a select region of the overall pattern that features key discontinuities) that allow ultrawideline spectra to be acquired rapidly, even for unreceptive samples (i.e., those with long values of T₁(¹⁴N), short values of T₂^eff(¹⁴N), or very broad patterns). It is hoped that these experimental and computational protocols will be useful for the characterization of various solid forms of nutraceuticals (i.e., salts, polymorphs, hydrates, solvates, cocrystals, amorphous solid dispersions, etc.), help detect impurity and counterfeit solid phases in dosage forms, and serve as a foundation for future NMR crystallographic studies of nutraceutical solid forms, including studies using ab initio crystal structure prediction algorithms.

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

The ability to theoretically predict accurate NMR chemical shifts in solids is increasingly important due to the role such shifts play in selecting among proposed model structures. Herein, two theoretical methods are evaluated for their ability to assign ¹⁵ N shifts from guanosine dihydrate to one of the two independent molecules present in the lattice. The NMR data consist of ¹⁵ N shift tensors from 10 resonances. Analysis using periodic boundary or fragment methods consider a benchmark dataset to estimate errors and predict uncertainties of 5.6 and 6.2 ppm, respectively. Despite this high accuracy, only one of the five sites were confidently assigned to a specific molecule of the asymmetric unit. This limitation is not due to negligible differences in experimental data, as most sites exhibit differences of >6.0 ppm between pairs of resonances representing a given position. Instead, the theoretical methods are insufficiently accurate to make assignments at most positions.

Refining crystal structures using 13C NMR chemical shift tensors as a target function

A two-step process is described for refining crystal structures from any source.

Distinguishing the Structures of High-Pressure Hydrides with Nuclear Magnetic Resonance Spectroscopy

The structural characterization of high-pressure hydrides has encountered many difficulties mainly due to the weak X-ray scattering of hydrogen. Herein, we investigate the prospect of detecting the H₃S and LaH₁₀ structures with nuclear magnetic resonance (NMR) spectroscopy. Our calculations demonstrate that the different candidate structures of H₃S (or LaH₁₀) exhibit significant differences in the electric field gradient (EFG) tensor of the ³³S (or ¹³⁹La) sites, indicating that the NMR spectroscopy can well capture the structural differences, even the small changes in the atomic position, and hence can be used to effectively probe the structures and the phase transitions of H₃S and LaH₁₀. Our results clarify the relationship between the structures and the EFG tensor parameters and provide a potential means to detect the structures of high-pressure hydrides.

Towards Accurate Predictions of Proton NMR Spectroscopic Parameters in Molecular Solids

The factors contributing to the accuracy of quantum-chemical calculations for the prediction of proton NMR chemical shifts in molecular solids are systematically investigated. Proton chemical shifts of six solid amino acids with hydrogen atoms in various bonding environments (CH, CH₂ , CH₃ , OH, SH and NH₃ ) were determined experimentally using ultra-fast magic-angle spinning and proton-detected 2D NMR experiments. The standard DFT method commonly used for the calculations of NMR parameters of solids is shown to provide chemical shifts that deviate from experiment by up to 1.5 ppm. The effects of the computational level (hybrid DFT functional, coupled-cluster calculation, inclusion of relativistic spin-orbit coupling) are thoroughly discussed. The effect of molecular dynamics and nuclear quantum effects are investigated using path-integral molecular dynamics (PIMD) simulations. It is demonstrated that the accuracy of the calculated proton chemical shifts is significantly better when these effects are included in the calculations.

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.

Mechanical and electromechanical properties of functionalized hexagonal boron nitride nanosheet: A density functional theory study

Hydroxylation as a technique is mainly used to alter the chemical characteristics of hexagonal boron nitride (h-BN), affecting physical features as well as mechanical and electromechanical properties in the process, the extent of which remains unknown. In this study, effects of functionalization on the physical, mechanical, and electromechanical properties of h-BN, including the interlayer distance, Young's modulus, intrinsic strength, and bandgaps were investigated based on density functional theory. It was found that functionalized layers of h-BN have an average distance of about 5.48 Å. Analyzing mechanical properties of h-BN revealed great dependence on the degree of functionalization. For the amorphous hydroxylated hexagonal boron nitride nanosheets (OH-BNNS), the Young's modulus moves from 436 to 284 GPa as the coverage of -OH increases. The corresponding variations in the Young's modulus of the ordered OH-BNNS with analogous coverage are bigger at 460-290 GPa. The observed intrinsic strength suggested that mechanical properties are promising even after functionalization. Moreover, the resulted bandgap reduction drastically enhanced the electrical conductivity of this structure under imposed strains. The results from this work pave the way for future endeavors in h-BN nanocomposites research.

Predicting anisotropic thermal displacements for hydrogens from solid-state NMR: a study on hydrogen bonding in polymorphs of palmitic acid

The hydrogen-bonding environments at the COOH moiety in eight polycrystalline polymorphs of palmitic acid are explored using solid-state NMR. Although most phases have no previously reported crystal structure, measured 13C chemical shift tensors for COOH moieties, combined with DFT modeling establish that all phases crystallize with a cyclic dimer (R22(8)) hydrogen bonding arrangement. Phases A2, Bm and Em have localized OH hydrogens while phase C has a dynamically disordered OH hydrogen. The phase designated As is a mix of five forms, including 27.4% of Bm and four novel phases not fully characterized here due to insufficient sample mass. For phases A2, Bm, Em, and C the anisotropic uncertainties in the COOH hydrogen atom positions are established using a Monte Carlo sampling scheme. Sampled points are retained or rejected at the ±1σ level based upon agreement of DFT computed 13COOH tensors with experimental values. The collection of retained hydrogen positions bear a remarkable resemblance to the anisotropic displacement parameters (i.e. thermal ellipsoids) from diffraction studies. We posit that this similarity is no mere coincidence and that the two are fundamentally related. The volumes of NMR-derived anisotropic displacement ellipsoids for phases with localized OH hydrogens are 4.1 times smaller than those derived from single crystal X-ray diffraction and 1.8 times smaller than the volume of benchmark single crystal neutron diffraction values.

Mechanochemical syntheses and 35Cl solid-state NMR characterization of fluoxetine HCl cocrystals

Novel mechanochemical syntheses of cocrystals of fluoxetine HCl are presented, along with characterization of the molecular-level structures by ³⁵Cl solid-state NMR and DFT calculations.

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.