The PAW/GIPAW approach for computing NMR parameters: a new dimension added to NMR study of solids

First-principles investigation of positively charged and neutral oxygen vacancies in amorphous silica

The structural parameters, electron localization functions, electron paramagnetic resonance (EPR) parameters, formation energies, and thermodynamic transition levels of various oxygen vacancy defects in amorphous silica are comprehensively and integrally investigated by using density functional theory. The trends of changes in the oxygen vacancy defect structure and electron localization induced by the increase in distance between defective silicon atoms are clearly identified. It is shown that the dimer configuration may be the potential structure of the Eδ' center. For the back-projected unpuckered configuration and the puckered configuration, whose EPR parameters are more consistent with the experimental values of the Eγ' center, the unpaired electron localized on the sp3 hybridized silicon atom is a common feature. Due to the three-coordinated oxygen atom in the forward-oriented configuration, the EPR parameters are closest to those of the Eα' center. Transformations of oxygen vacancy defects under different charge states are studied by sequentially adding and removing electrons. The thermodynamic transition level analysis reveals that the dimer and forward configurations may behave as deep traps for electron accumulation. The back-projected puckered fourfold-coordinated and fivefold-coordinated configurations are comparatively stable and may be able to function as shallow traps for electron transport. The neutral double unpuckered, neutral back-projected puckered fourfold-coordinated, and neutral back-projected unpuckered configurations are more likely to lose electrons during hole trapping. As the bias voltage is repeatedly changed, the defect density of the puckered configuration may reduce, while that of the dimer and unpuckered configuration may take an opposite trend.

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

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

Can simple 'molecular' corrections outperform projector augmented-wave density functional theory in the prediction of 35 Cl electric field gradient tensor parameters for chlorine-containing crystalline systems?

Many-body expansion (MBE) fragment approaches have been applied to accurately compute nuclear magnetic resonance (NMR) parameters in crystalline systems. Recent examples demonstrate that electric field gradient (EFG) tensor parameters can be accurately calculated for 14 N and 17 O. A key additional development is the simple molecular correction (SMC) approach, which uses two one-body fragment (i.e., isolated molecule) calculations to adjust NMR parameter values established using 'benchmark' projector augmented-wave (PAW) density functional theory (DFT) values. Here, we apply a SMC using the hybrid PBE0 exchange-correlation (XC) functional to see if this can improve the accuracy of calculated 35 Cl EFG tensor parameters. We selected eight organic and two inorganic crystal structures and considered 15 chlorine sites. We find that this SMC improves the accuracy of computed values for both the 35 Cl quadrupolar coupling constant (CQ ) and the asymmetry parameter ( η Q ) by approximately 30% compared with benchmark PAW DFT values. We also assessed a SMC that offers local improvements not only in terms of the quality of the XC functional but simultaneously in the quality of the description of relativistic effects via the inclusion of spin-orbit effects. As the inorganic systems considered contain heavy atoms bonded to the chlorine atoms, we find further improvements in the accuracy of calculated 35 Cl EFG tensor parameters when both a hybrid functional and spin-orbit effects are included in the SMC. On the contrary, for chlorine-containing organics, the inclusion of spin-orbit relativistic effects using a SMC does not improve the accuracy of computed 35 Cl EFG tensor parameters.

Accurate and Transferable Machine Learning Potential for Molecular Dynamics Simulation of Sodium Silicate Glasses

An accurate and transferable machine learning (ML) potential for the simulation of binary sodium silicate glasses over a wide range of compositions (from 0 to 50% Na2O) was developed. The potential energy surface is approximated by the sum of atomic energy contributions mapped by a neural network algorithm from the local geometry comprising information on atomic distances and angles with neighboring atoms using the DeePMD code [Wang, H. Comput. Phys. Commun. 2018, 228, 178-184]. Our model was trained on a large data set of total energies and atomic forces computed at the density functional theory level on structures extracted from classical molecular dynamics (MD) simulations performed at several temperatures from 300 to 3000 K. This allows for the generation of a robust and transferable ML potential applicable over the full compositional range of glass formability at different temperatures that outperforms the empirical potentials available in the literature in reproducing structures and properties such as bond angle distribution, total distribution functions, and vibrational density of state. The generality of the approach enables the future training of a potential with other or more elements allowing for simulations of structures, properties, and behavior of ternary and multicomponent oxide glasses with nearly ab initio accuracy at a fraction of the computational cost.

MIM-ML: A Novel Quantum Chemical Fragment-Based Random Forest Model for Accurate Prediction of NMR Chemical Shifts of Nucleic Acids

We developed a random forest machine learning (ML) model for the prediction of 1H and 13C NMR chemical shifts of nucleic acids. Our ML model is trained entirely on reproducing computed chemical shifts obtained previously on 10 nucleic acids using a Molecules-in-Molecules (MIM) fragment-based density functional theory (DFT) protocol including microsolvation effects. Our ML model includes structural descriptors as well as electronic descriptors from an inexpensive low-level semiempirical calculation (GFN2-xTB) and trained on a relatively small number of DFT chemical shifts (2080 1H chemical shifts and 1780 13C chemical shifts on the 10 nucleic acids). The ML model is then used to make chemical shift predictions on 8 new nucleic acids ranging in size from 600 to 900 atoms and compared directly to experimental data. Though no experimental data was used in the training, the performance of our model is excellent (mean absolute deviation of 0.34 ppm for 1H chemical shifts and 2.52 ppm for 13C chemical shifts for the test set), despite having some nonstandard structures. A simple analysis suggests that both structural and electronic descriptors are critical for achieving reliable predictions. This is the first attempt to combine ML from fragment-based DFT calculations to predict experimental chemical shifts accurately, making the MIM-ML model a valuable tool for NMR predictions of nucleic acids.

Structural Modeling of O/F Correlated Disorder in TaOF3 and NbOF3-x(OH)x by Coupling Solid-State NMR and DFT Calculations

The structure of MOF3 (M = Nb, Ta) compounds was precisely modeled by combining powder X-ray diffraction, solid-state NMR spectroscopy, and semiempirical dispersion-corrected DFT calculations. It consists of stacked ∞(MOF3) layers along the c⃗ direction formed by heteroleptic corner-connected MX6 (X = O, F) octahedra. 19F NMR resonance assignments and occupancy rates of the anionic crystallographic sites have been revised. The bridging site is shared equally by the anions, and the terminal site is occupied by F only. An O/F correlated disorder is expected since cis-MO2F4 octahedra are favored, resulting in one-dimensional -F-M-O-M- strings along the <100> and <010> directions. Ten different 2 × 2 × 1 supercells per compound, fulfilling these characteristics, were built. Using DFT calculations and the GIPAW approach, the supercells were relaxed and the 19F isotropic chemical shift values were determined. The agreement between the experimental and calculated 19F spectra is excellent for TaOF3. The 1H and 19F experimental NMR spectra revealed that some of the bridging F atoms are substituted by OH groups, especially in NbOF3. New supercells involving OH groups were generated. Remarkably, the best agreement is obtained for the supercells with the composition closest to that estimated from the 19F NMR spectra, i.e., NbOF2.85(OH)0.15.

Unravelling the effect of paramagnetic Ni2+ on the 13C NMR shift tensor for carbonate in Mg2-xNixAl layered double hydroxides by quantum-chemical computations

Structural disorder and low crystallinity render it challenging to characterise the atomic-level structure of layered double hydroxides (LDH). We report a novel multi-step, first-principles computational workflow for the analysis of paramagnetic solid-state NMR of complex inorganic systems such as LDH, which are commonly used as catalysts and energy storage materials. A series of 13CO32--labelled Mg2-xNixAl-LDH, x ranging from 0 (Mg2Al-LDH) to 2 (Ni2Al-LDH), features three distinct eigenvalues δ11, δ22 and δ33 of the experimental 13C chemical shift tensor. The δii correlate directly with the concentration of the paramagnetic Ni2+ and span a range of |δ11 - δ33| ≈ 90 ppm at x = 0, increasing to 950 ppm at x = 2. In contrast, the isotropic shift, δiso(13C), only varies by -14 ppm in the series. Detailed insight is obtained by computing (1) the orbital shielding by periodic density-functional theory involving interlayer water, (2) the long-range pseudocontact contribution of the randomly distributed Ni2+ ions in the cation layers (characterised by an ab initio susceptibility tensor) by a lattice sum, and (3) the close-range hyperfine terms using a full first-principles shielding machinery. A pseudohydrogen-terminated two-layer cluster model is used to compute (3), particularly the contact terms. Due to negative spin density contribution at the 13C site arising from the close-by Ni2+ sites, this step is necessary to reach a semiquantitative agreement with experiment. These findings influence future NMR investigations of the formally closed-shell interlayer species within LDH, such as the anions or water. Furthermore, the workflow is applicable to a variety of complex materials.

Review: understanding the properties of amorphous materials with high-performance computing methods

Amorphous materials have no long-range order in their atomic structure. This makes much of the formalism for the study of crystalline materials irrelevant, and so elucidating their structure and properties is challenging. The use of computational methods is a powerful complement to experimental studies, and in this paper we review the use of high-performance computing methods in the simulation of amorphous materials. Five case studies are presented to showcase the wide range of materials and computational methods available to practitioners in this field. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

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.

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

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

3D to 0D cesium lead bromide: A 79/81Br NMR, NQR and theoretical investigation

Inorganic metal halides offer unprecedented tunability through elemental variation of simple three-element compositions, but can exhibit complicated phase behaviour, degradation, and microscopic phenomena (disorder/dynamics) that play an integral role for the bulk-level chemical and physical properties of these materials. Understanding the halogen chemical environment in such materials is crucial to addressing many of the concerns regarding implementing these materials in commercial applications. In this study, a combined solid-state nuclear magnetic resonance, nuclear quadrupole resonance and quantum chemical computation approach is used to interrogate the Br chemical environment in a series of related inorganic lead bromide materials: CsPbBr₃, CsPb₂Br₅, and Cs₄PbBr₆. The quadrupole coupling constants (C_Q) were determined to range from 61 to 114 MHz for ⁸¹Br, with CsPbBr₃ exhibiting the largest measured C_Q and Cs₄PbBr₆ the smallest. GIPAW DFT was shown to be an excellent pre-screening tool for estimating the EFG of Br materials and can increase experimental efficiency by providing good starting estimates for acquisition. Finally, the combination of theory and experiment to inform the best methods for expanding further to the other quadrupolar halogens is discussed.

Disorder to order: how halide mixing in MAPbI3-xBrx perovskites restricts MA dynamics

Mixed-halide lead perovskites are of particular interest for the design of tandem solar cells currently reaching record efficiencies. While halide phase segregation upon illumination of mixed perovskites is extensively studied, the effect of halide disorder on A cation dynamics is not well understood, despite its importance for charge carrier diffusion and lifetime. Here, we study the methylammonium (MA) reorientational dynamics in mixed halide MAPbI_3-xBr_x perovskites by a combined approach of experimental solid-state NMR spectroscopy and molecular dynamics (MD) simulations based on machine-learning force-fields (MLFF). ²⁰⁷Pb NMR spectra indicate the halides are randomly distributed over their lattice positions, whereas PXRD measurements show that all mixed MAPbI_3-xBr_x samples are cubic. The experimental ¹⁴N spectra and ¹H double-quantum (DQ) NMR data reveal anisotropic MA reorientations depending on the halide composition and thus associated disorder in the inorganic sublattice. MD calculations allow us to correlate these experimental results to restrictions of MA dynamics due to preferred MA orientations in their local Pb₈I_12-nBr_n "cages". Based on the experimental and simulated results, we develop a phenomenological model that correlates the ¹H dipolar coupling and thus the MA dynamics with the local composition and reproduces the experimental data over the whole composition range. We show that the dominant interaction between the MA cations and the Pb-X lattice that influences the cation dynamics is the local electrostatic potential being inhomogeneous in mixed halide systems. As such, we generate a fundamental understanding of the predominant interaction between the MA cations and the inorganic sublattice, as well as MA dynamics in asymmetric halide coordinations.

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.

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.

Quadrupolar 23Na+ NMR relaxation as a probe of subpicosecond collective dynamics in aqueous electrolyte solutions

Nuclear magnetic resonance relaxometry represents a powerful tool for extracting dynamic information. Yet, obtaining links to molecular motion is challenging for many ions that relax through the quadrupolar mechanism, which is mediated by electric field gradient fluctuations and lacks a detailed microscopic description. For sodium ions in aqueous electrolytes, we combine ab initio calculations to account for electron cloud effects with classical molecular dynamics to sample long-time fluctuations, and obtain relaxation rates in good agreement with experiments over broad concentration and temperature ranges. We demonstrate that quadrupolar nuclear relaxation is sensitive to subpicosecond dynamics not captured by previous models based on water reorientation or cluster rotation. While ions affect the overall water retardation, experimental trends are mainly explained by dynamics in the first two solvation shells of sodium, which contain mostly water. This work thus paves the way to the quantitative understanding of quadrupolar relaxation in electrolyte and bioelectrolyte systems.

Recent progress in solid-state NMR of spin-½ low-γ nuclei applied to inorganic materials

Significant technological and methodological advances in solid-state NMR techniques in recent years have increased the accessibility of nuclei with small magnetic moments (hereafter termed low-γ) underpinning an increased range of applications of such nuclei. These methodological advances are briefly summarised, including improvements in hardware and pulse sequences, as well as important developments in associated computational methods (e.g. first principles calculations, spectral simulation). Here spin-½ nuclei are the focus, with this Perspective complementing a very recent review that looked at half-integer spin low-γ quadrupolar nuclei. Reference is made to some of the original reports of such spin-½ nuclei, but recent progress in the relevant methodology and applications to inorganic materials (most within the last 10 years) of these nuclei are the focus. An overview of the current state-of-the-art of studying these nuclei is thereby provided for both NMR spectroscopists and materials researchers.

Dramatic Enhancement of Rare-Earth Metal-Organic Framework Stability Via Metal Cluster Fluorination

Rare-earth polynuclear metal-organic frameworks (RE-MOFs) have demonstrated high durability for caustic acid gas adsorption and separation based on gas adsorption to the metal clusters. The metal clusters in the RE-MOFs traditionally contain RE metals bound by μ₃-OH groups connected via organic linkers. Recent studies have suggested that these hydroxyl groups could be replaced by fluorine atoms during synthesis that includes a fluorine-containing modulator. Here, a combined modeling and experimental study was undertaken to elucidate the role of metal cluster fluorination on the thermodynamic stability, structure, and gas adsorption properties of RE-MOFs. Through systematic density-functional theory calculations, fluorinated clusters were found to be thermodynamically more stable than hydroxylated clusters by up to 8-16 kJ/mol per atom for 100% fluorination. The extent of fluorination in the metal clusters was validated through a ¹⁹F NMR characterization of 2,5-dihydroxyterepthalic acid (Y-DOBDC) MOF synthesized with a fluorine-containing modulator. ¹⁹F magic-angle spinning NMR identified two primary peaks in the isotropic chemical shift (δ_iso) spectra located at -64.2 and -69.6 ppm, matching calculated ¹⁹F NMR δ_iso peaks at -63.0 and -70.0 ppm for fluorinated systems. Calculations also indicate that fluorination of the Y-DOBDC MOF had negligible effects on the acid gas (SO₂, NO₂, H₂O) binding energies, which decreased by only ∼4 kJ/mol for the 100% fluorinated structure relative to the hydroxylated structure. Additionally, fluorination did not change the relative gas binding strengths (SO₂ > H₂O > NO₂). Therefore, for the first time the presence of fluorine in the metal clusters was found to significantly stabilize RE-MOFs without changing their acid-gas adsorption properties.

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

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

Leucine on Silica: A Combined Experimental and Modeling Study of a System Relevant for Origins of Life, and the Role of Water Coadsorption

Leucine on silica constitutes an interesting system from the point of view of origins of life studies since leucine coadsorbed on SiO₂ together with glutamic acid can give rise to rather long linear polypeptides upon activation. It is also an ideal system to test methods of molecular characterization of biomolecules deposited on mineral surfaces since it combines a small-scale model of peptides and proteins, which are among the most important components of biodevices, with one of the most widely used inorganic materials. We have deposited l-leucine on a high surface fumed silica in the submonolayer range and characterized it by a multipronged approach including macroscopic information (thermogravimetry, X-ray diffraction), in situ spectroscopic methods (IR, multinuclear solid-state NMR including single-pulse and CP-MAS, 2-D HETCOR), and molecular modeling by density functional theory (DFT), including calculation of NMR parameters. Specific information can be obtained on the adsorption and interaction mechanism. Leucine is rather strongly adsorbed without any covalent bonds, through the formation of a specific lattice of H-bonds that often involve coadsorbed water molecules. Its state is indeed strongly dependent on the drying procedure: insufficient drying results in liquid-like surroundings for the leucine functional groups, while vacuum drying only retains a limited number of waters (of the order of 5 per leucine molecule). The most stable models have zwitterionic leucine interacting directly with surface silanols through their ammonium group, while the carboxylate interacts through bridging waters. Experimental NMR chemical shifts are satisfactorily predicted for these models, and leucine can be viewed as a probe for specific groups of surface sites known as silanol nests.

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.

Regression Machine Learning Models Used to Predict DFT-Computed NMR Parameters of Zeolites

Machine learning approaches can drastically decrease the computational time for the predictions of spectroscopic properties in materials, while preserving the quality of the computational approaches. We studied the performance of kernel-ridge regression (KRR) and gradient boosting regressor (GBR) models trained on the isotropic shielding values, computed with density-functional theory (DFT), in a series of different known zeolites containing out-of-frame metal cations or fluorine anion and organic structure-directing cations. The smooth overlap of atomic position descriptors were computed from the DFT-optimised Cartesian coordinates of each atoms in the zeolite crystal cells. The use of these descriptors as inputs in both machine learning regression methods led to the prediction of the DFT isotropic shielding values with mean errors within 0.6 ppm. The results showed that the GBR model scales better than the KRR model.

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.

GIPAW Pseudopotentials of d Elements for Solid-State NMR

Computational methods are increasingly used to support interpreting, assigning and predicting the solid-state nuclear resonance magnetic spectra of materials. Currently, density functional theory is seen to achieve a good balance between efficiency and accuracy in solid-state chemistry. To be specific, density functional theory allows the assignment of signals in nuclear resonance magnetic spectra to specific sites and can help identify overlapped or missing signals from experimental nuclear resonance magnetic spectra. To avoid the difficulties correlated to all-electron calculations, a gauge including the projected augmented wave method was introduced to calculate nuclear resonance magnetic parameters with great success in organic crystals in the last decades. Thus, we developed a gauge including projected augmented pseudopotentials of 21 d elements and tested them on, respectively, oxides or nitrides (semiconductors), calculating chemical shift and quadrupolar coupling constant. This work can be considered the first step to improving the ab initio prediction of nuclear magnetic resonance parameters, and leaves open the possibility for inorganic compounds to constitute an alternative standard compound, with respect to tetramethylsilane, to calculate the chemical shift. Furthermore, this work represents the possibility to obtain results from first-principles calculations, to train a machine-learning model to solve or refine structures using predicted nuclear magnetic resonance spectra.

De Novo Crystal Structure Determination from Machine Learned Chemical Shifts

Determination of the three-dimensional atomic-level structure of powdered solids is one of the key goals in current chemistry. Solid-state NMR chemical shifts can be used to solve this problem, but they are limited by the high computational cost associated with crystal structure prediction methods and density functional theory chemical shift calculations. Here, we successfully determine the crystal structures of ampicillin, piroxicam, cocaine, and two polymorphs of the drug molecule AZD8329 using on-the-fly generated machine-learned isotropic chemical shifts to directly guide a Monte Carlo-based structure determination process starting from a random gas-phase conformation.

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.

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

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

89Y chemical shift anisotropy: a sensitive structural probe of layered yttrium hydroxides revealed by solid-state NMR spectroscopy and DFT calculations

Anion-exchangeable Y₂(OH)₅X·nH₂O (LYH-X, X = monovalent anions, n ≈ 1.5) materials are an ideal platform for incorporating the unique properties of layered metal hydroxides and rare-earth (RE) ions, and thus have exhibited promising prospects for various applications. To further improve the performance of LYH-X and related functional materials, their structure-property relationships must be explored. However, due to the intrinsic felxibility, extracting the local structural details of these materials is particularly challenging. In this work, we utilized a combined approach of ⁸⁹Y solid-state NMR (ssNMR) spectroscopy and density functional theory (DFT) calculations to reveal the response of ⁸⁹Y chemical shift anisotropy (CSA) in LYH-X to the structural changes including a small displacement of cationic yttrium hydroxide layers and intercalated anions. Such subtle structural changes are often associated with dehydration/rehydration, anion-exchange, exfoliation, and the self-assembly process of LYH-X and related functional materials, which are exceedingly difficult to detect using other techniques. The principal components of ⁸⁹Y CSA show a larger variation range than isotropic chemical shifts, making CSA a more sensitive probe. In addition, it is found that the response of ⁸⁹Y CSA to structural changes is distinct for Y sites with different local coordination environments, opening great opportunities to analyze each Y site within these materials. All these observations suggest that the strategy involving both experimental (⁸⁹Y ssNMR) and theoretical (DFT) approaches can be utilized to extract previously unavailable ultrafine structural information of LYH-X and related materials, and provide fruitful insights into their thorough structure-property relationships.

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.

Enhanced Fluoride Uptake by Layered Double Hydroxides under Alkaline Conditions: Solid-State NMR Evidence of the Role of Surface >MgOH Sites

Layered double hydroxides (LDHs) are potential low-cost filter materials for use in fluoride removal from drinking water, but molecular-scale defluoridation mechanisms are lacking. In this research, we employed ¹⁹F solid-state NMR spectroscopy to identify fluoride sorption products on 2:1 MgAl LDH and to reveal the relationship between fluoride sorption and the LDH structure. A set of six ¹⁹F NMR peaks centered at -140, -148, -156, -163, -176, and -183 ppm was resolved. Combining quantum chemical calculations based on density function theory (DFT) and ¹⁹F{²⁷Al} transfer of populations in double resonance (TRAPDOR) analysis, we could assign the peaks at -140, -148, -156, and -163 ppm to Al-F (F coordinated to surface Al) and those at -176 and -183 ppm to Mg-F (F coordinated to surface Mg only). Interestingly, the spectroscopic data reveal that the formation of Al-F is the predominant mode of F^- sorption at low pH, whereas the formation of Mg-F is predominant at high pH (or a higher Mg/Al ratio). This finding supports the fact that the F^- uptake of 2:1 MgAl LDH was nearly six times that of activated alumina at pH 9. Overall, we explicitly revealed the different roles of the surface >MgOH and >AlOH sites of LDHs in defluoridation, which explained why the use of classic activated alumina for defluoridation is limited at high pH. The findings from this research may also provide new insights into material screening for potential filters for F^- removal under alkaline conditions.

NMR Relaxation Rates of Quadrupolar Aqueous Ions from Classical Molecular Dynamics Using Force-Field Specific Sternheimer Factors

The nuclear magnetic resonance (NMR) relaxation of quadrupolar nuclei is governed by the electric field gradient (EFG) fluctuations at their position. In classical molecular dynamics (MD), the electron cloud contribution to the EFG can be included via the Sternheimer approximation, in which the full EFG at the nucleus that can be computed using quantum density functional theory (DFT) is considered to be proportional to that arising from the external, classical charge distribution. In this work, we systematically assess the quality of the Sternheimer approximation as well as the impact of the classical force field (FF) on the NMR relaxation rates of aqueous quadrupolar ions at infinite dilution. In particular, we compare the rates obtained using an ab initio parametrized polarizable FF, a recently developed empirical FF with scaled ionic charges and a simple empirical nonpolarizable FF with formal ionic charges. Surprisingly, all three FFs considered yield good values for the rates of smaller and less polarizable solutes (Li⁺, Na⁺, K⁺, Cl^-), provided that a model-specific Sternheimer parametrization is employed. Yet, the polarizable and scaled charge FFs yield better estimates for divalent and more polarizable species (Mg²⁺, Ca²⁺, Cs⁺). We find that a linear relationship between the quantum and classical EFGs holds well in all of the cases considered; however, such an approximation often leads to quite large errors in the resulting EFG variance, which is directly proportional to the computed rate. We attempted to reduce the errors by including first order nonlinear corrections to the EFG, yet no clear improvement for the resulting variance has been found. The latter result indicates that more refined methods for determining the EFG at the ion position, in particular those that take into account the instantaneous atomic environment around an ion, might be necessary to systematically improve the NMR relaxation rate estimates in classical MD.

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.

17 O solid-state NMR at ultrahigh magnetic field of 35.2 T: Resolution of inequivalent oxygen sites in different phases of MOF MIL-53(Al)

MIL-53(Al) is a member of the most extensively studied metal-organic framework (MOF) families owing to its "flexible" framework and superior stability. 17 O solid-state NMR (SSNMR) spectroscopy is an ideal site-specific characterization tool as it probes local oxygen environments. Because oxygen local structure is often altered during phase change, 17 O SSNMR can be used to follow phase transitions. However, 17 O is a challenging nucleus to study via SSNMR due to its low sensitivity and resolution arising from the very low natural abundance of 17 O isotope and its quadrupolar nature. In this work, we describe that by using 17 O isotopic enrichment and performing 17 O SSNMR experiments at an ultrahigh magnetic field of 35.2 T, all chemically and crystallographically inequivalent oxygen sites in two representative MIL-53(Al) (as-made and water adsorbed) phases can be completely resolved. The number of signals in each phase is consistent with that predicted from the space group refined from powder X-ray diffraction data. The 17 O 1D magic-angle spinning (MAS) and 2D triple-quantum MAS (3QMAS) spectra at 35.2 T furnish fine information about the host-guest interactions and the structural changes associated with phase transition. The ability to completely resolve multiple chemically and crystallographically inequivalent oxygen sites in MOFs at very high magnetic field, as illustrated in this work, significantly enhances the potential for using the NMR crystallography approach to determine crystal structures of new MOFs and verify the structures of existing MOFs obtained from refining powder X-ray diffraction data.

Recent progress in solid-state nuclear magnetic resonance of half-integer spin low-γ quadrupolar nuclei applied to inorganic materials

An overview is presented of recent progress in the solid-state nuclear magnetic resonance (NMR) observation of low-γ nuclei, with a focus on applications to inorganic materials. The technological and methodological advances in the last 20 years, which have underpinned the increased accessibility of low-γ nuclei for study by solid-state NMR techniques, are summarised, including improvements in hardware, pulse sequences and associated computational methods (e.g., first principles calculations and spectral simulation). Some of the key initial observations from inorganic materials of these nuclei are highlighted along with some recent (most within the last 10 years) illustrations of their application to such materials. A summary of other recent reviews of the study of low-γ nuclei by solid-state NMR is provided so that a comprehensive understanding of what has been achieved to date is available.

Solid-state nuclear magnetic resonance study of polymorphism in tris(8-hydroxyquinolinate)aluminium

Tris(8-hydroxyquinolinate)aluminium (Alq₃ ) is a metal-organic coordination complex, which is a widely used electroluminescent material in organic light-emitting diode technology. Crystalline Alq₃ is known to occur in five polymorphic forms (denoted α, β, γ, δ, and ε), although the structures of some of these polymorphs have been the subject of considerable debate. In particular, the structure of α-Alq₃ , which is a model for the local structure in amorphous films used in devices, is highly complex and has never been conclusively solved. In this work, we use solid-state nuclear magnetic resonance (NMR) and density functional theory (DFT) calculations to investigate the local structure of four Alq₃ samples. We find that the first structure proposed for α-Alq₃ is inconsistent with all of the samples studied, and DFT calculations further suggest that this structure is energetically unfavourable. Instead, samples containing the meridional (mer) isomeric form are found to contain local structures consistent with ε-Alq₃ , and a sample containing the facial (fac) isomeric form is consistent with a mixture of γ-Alq₃ and δ-Alq₃ . We also investigate the influence of different strategies for dispersion correction in DFT geometry optimisations. We find that a recently proposed modified semiempirical dispersion correction scheme gives good agreement with experiment. Furthermore, the DFT calculations also show that distinction between mer and fac isomers on the basis of η_Q that has been assumed in previous work is not always justified.

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

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

Analyzing Discrepancies in Chemical-Shift Predictions of Solid Pyridinium Fumarates

Highly accurate chemical-shift predictions in molecular solids are behind the success and rapid development of NMR crystallography. However, unusually large errors of predicted hydrogen and carbon chemical shifts are sometimes reported. An understanding of these deviations is crucial for the reliability of NMR crystallography. Here, recently reported large deviations of predicted hydrogen and carbon chemical shifts of a series of solid pyridinium fumarates are thoroughly analyzed. The influence of the geometry optimization protocol and of the computational level of NMR calculations on the accuracy of predicted chemical shifts is investigated. Periodic calculations with GGA, meta-GGA and hybrid functionals are employed. Furthermore, molecular corrections at the coupled-cluster singles-and-doubles (CCSD) level are calculated. The effect of nuclear delocalization on the structure and NMR shielding is also investigated. The geometry optimization with a computationally demanding hybrid functional leads to a substantial improvement in proton chemical-shift predictions.

Perfect and Defective 13C-Furan-Derived Nanothreads from Modest-Pressure Synthesis Analyzed by 13C NMR

The molecular structure of nanothreads produced by the slow compression of 13C4-furan was studied by advanced solid-state NMR. Spectral editing showed that >95% of carbon atoms were bonded to one hydrogen (C-H) and that there were 2-4% CH2, 0.6% C═O, and <0.3% CH3 groups. Alkenes accounted for 18% of the CH moieties, while trapped, unreacted furan made up 7%. Two-dimensional (2D) 13C-13C and 1H-13C NMR identified 12% of all carbon in asymmetric O-CH═CH-CH-CH- and 24% in symmetric O-CH-CH═CH-CH- rings. While the former represented defects or chain ends, some of the latter appeared to form repeating thread segments. Around 10% of carbon atoms were found in highly ordered, fully saturated nanothread segments. Unusually slow 13C spin-exchange with sites outside the perfect thread segments documented a length of at least 14 bonds; the small width of the perfect-thread signals also implied a fairly long, regular structure. Carbons in the perfect threads underwent relatively slow spin-lattice relaxation, indicating slow spin exchange with other threads and smaller amplitude motions. Through partial inversion recovery, the signals of the perfect threads were observed and analyzed selectively. Previously considered syn-threads with four different C-H bond orientations were ruled out by centerband-only detection of exchange NMR, which was, on the contrary, consistent with anti-threads. The observed 13C chemical shifts were matched well by quantum-chemical calculations for anti-threads but not for more complex structures like syn/anti-threads. These observations represent the first direct determination of the atomic-level structure of fully saturated nanothreads.

Structure determination of an amorphous drug through large-scale NMR predictions

Knowledge of the structure of amorphous solids can direct, for example, the optimization of pharmaceutical formulations, but atomic-level structure determination in amorphous molecular solids has so far not been possible. Solid-state nuclear magnetic resonance (NMR) is among the most popular methods to characterize amorphous materials, and molecular dynamics (MD) simulations can help describe the structure of disordered materials. However, directly relating MD to NMR experiments in molecular solids has been out of reach until now because of the large size of these simulations. Here, using a machine learning model of chemical shifts, we determine the atomic-level structure of the hydrated amorphous drug AZD5718 by combining dynamic nuclear polarization-enhanced solid-state NMR experiments with predicted chemical shifts for MD simulations of large systems. From these amorphous structures we then identify H-bonding motifs and relate them to local intermolecular complex formation energies.

A Review on Combination of Ab Initio Molecular Dynamics and NMR Parameters Calculations

This review focuses on a combination of ab initio molecular dynamics (aiMD) and NMR parameters calculations using quantum mechanical methods. The advantages of such an approach in comparison to the commonly applied computations for the structures optimized at 0 K are presented. This article was designed as a convenient overview of the applied parameters such as the aiMD type, DFT functional, time step, or total simulation time, as well as examples of previously studied systems. From the analysis of the published works describing the applications of such combinations, it was concluded that including fast, small-amplitude motions through aiMD has a noticeable effect on the accuracy of NMR parameters calculations.

Synthesis of the scandium chloride hydrates ScCl3·3H2O and Sc2Cl4(OH)2·12H2O and their characterisation by X-ray diffraction, 45Sc NMR spectroscopy and DFT calculations

The compounds ScCl3·3H2O (SCTH) and [{Sc(H2O)5(μ-OH)}2]Cl4·2H2O (SCOH), have been synthesised and characterised by single-crystal XRD, 45Sc NMR spectroscopy and DFT calculations, with the crystal structure of SCTH reported here for the first time. From 45Sc NMR measurements under static and MAS conditions, both chemical shift and quadrupolar coupling parameters have been determined. The quadrupolar coupling constants χ for the octahedrally coordinated scandium sites in SCTH are 2.0 ± 0.1 MHz for Sc(1) and 3.81 ± 0.05 MHz for Sc(2). For SCOH, where the hepta-coordination of the single scandium site constitutes a less symmetric electronic environment, 14.68 ± 0.05 MHz was found. DFT calculations for the static SCTH structure consistently overestimate the quadrupolar coupling constants, indicating the possible presence of crystal water dynamics on the NMR time scale.

Acidity and Local Confinement Effect in Mordenite Probed by Solid-State NMR Spectroscopy

Herein, utilizing acetonitrile as the probe molecule, the acidity and host-guest interactions of H-mordenite (H-MOR) zeolites are investigated comprehensively by solid-state NMR spectroscopy and theoretical calculation. The locations and local configurations of Brønsted acid sites (BASs) in H-MOR are revealed by multinuclear and multidimensional NMR experiments with adsorption/coadsorption of acetonitrile (CD₃CN) and trimethylphosphine (TMP). Moreover, the confinement effect of dual pores in MOR has been characterized via the quantitative determination of host-guest interactions between CH₃CN and BASs. The ¹H-¹⁵N dipolar measurement results and DFT calculations demonstrate that there are two kinds of acetonitrile molecules adsorbed in 12-membered ring (12MR) main channels with distinct mobility, where acetonitrile undergoes either partially restricted or highly flexible motion in the time scale of nanoseconds to microseconds. These two types of acetonitrile can exchange with temperature rising. In contrast, the mobility of acetonitrile in 8-membered ring (8MR) channels is very restricted due to the confinement of the framework.

Unveiling the stability of Sn/Si/graphite composites for Li-ion storage by physical, electrochemical and computational tools

Complex materials composed of two and three elements with high Li-ion storage capacity are investigated and tested as lithium-ion battery (LiB) negative electrodes. Namely, anodes containing tin, silicon, and graphite show very good performance because of the large gravimetric and volumetric capacity of silicon and structural support provided by tin and graphite. The performance of the composites during the first cycles was studied using ex situ magic angle spinning (MAS) ⁷Li Nuclear Magnetic Resonance (NMR), density functional theory (DFT) calculations, and electrochemical techniques. The best performance was obtained for Sn/Si/graphite in a 1 : 1 : 1 proportion, due to an emergent effect of the interaction between Sn and Si. The results suggest a stabilization effect of Sn over Si, providing a physical constraint that prevents Si pulverization. This mechanism ensures good cyclability over more than one hundred cycles, low capacity fading and high specific capacity.

Microscopic (Dis)order and Dynamics of Cations in Mixed FA/MA Lead Halide Perovskites

Recent developments in the field of high efficiency perovskite solar cells are based on stabilization of the perovskite crystal structure of FAPbI₃ while preserving its excellent optoelectronic properties. Compositional engineering of, for example, MA or Br mixed into FAPbI₃ results in the desired effects, but detailed knowledge of local structural features, such as local (dis)order or cation interactions of formamidinium (FA) and methylammonium (MA), is still limited. This knowledge is, however, crucial for their further development. Here, we shed light on the microscopic distribution of MA and FA in mixed perovskites MA_1-x FA _x PbI₃ and MA_0.15FA_0.85PbI_2.55Br_0.45 by combining high-resolution double-quantum ¹H solid-state nuclear magnetic resonance (NMR) spectroscopy with state-of-the-art near-first-principles accuracy molecular dynamics (MD) simulations using machine-learning force-fields (MLFFs). We show that on a small local scale, partial MA and FA clustering takes place over the whole MA/FA compositional range. A reasonable driving force for the clustering might be an increase of the dynamical freedom of FA cations in FA-rich regions. While MA_0.15FA_0.85PbI_2.55Br_0.45 displays similar MA and FA ordering as the MA_1-x FA _x PbI₃ systems, the average cation-cation interaction strength increased significantly in this double mixed material, indicating a restriction of the space accessible to the cations or their partial immobilization upon Br^- incorporation. Our results shed light on the heterogeneities in cation composition of mixed halide perovskites, helping to exploit their full optoelectronic potential.

Observing the three-dimensional dynamics of supported metal complexes

The dynamics of heterogeneous catalysts are linked to their activity and selectivity but are poorly understood. NMR enables for the determination of high-resolution dynamic structures for such sites and the mapping of accessible conformations.

Unveiling the Structure and Reactivity of Fatty-Acid Based (Nano)materials Thanks to Efficient and Scalable 17O and 18O-Isotopic Labeling Schemes

Fatty acids are ubiquitous in biological systems and widely used in materials science, including for the formulation of drugs and the surface-functionalization of nanoparticles. However, important questions regarding the structure and reactivity of these molecules are still to be elucidated, including their mode of binding to certain metal cations or materials surfaces. In this context, we have developed novel, efficient, user-friendly, and cost-effective synthetic protocols based on ball-milling, for the 17O and 18O isotopic labeling of two key fatty acids which are widely used in (nano)materials science, namely stearic and oleic acid. Labeled molecules were analyzed by 1H and 13C solution NMR, IR spectroscopy, and mass spectrometry (ESI-TOF and LC-MS), as well as 17O solid state NMR (for the 17O labeled species). In both cases, the labeling procedures were scaled-up to produce up to gram quantities of 17O- or 18O-enriched molecules in just half-a-day, with very good synthetic yields (all ≥84%) and enrichment levels (up to an average of 46% per carboxylic oxygen). The 17O-labeled oleic acid was then used for the synthesis of a metal soap (Zn-oleate) and the surface-functionalization of ZnO nanoparticles (NPs), which were characterized for the first time by high-resolution 17O NMR (at 14.1 and 35.2 T). This allowed very detailed insight into (i) the coordination mode of the oleate ligand in Zn-oleate to be achieved (including information on Zn···O distances) and (ii) the mode of attachment of oleic-acid at the surface of ZnO (including novel information on its photoreactivity upon UV-irradiation). Overall, this work demonstrates the high interest of these fatty acid-enrichment protocols for understanding the structure and reactivity of a variety of functional (nano)materials systems using high resolution analyses like 17O NMR.