New perspectives in the PAW/GIPAW approach: J(P-O-Si) coupling constants, antisymmetric parts of shift tensors and NQR predictions

Solid-State Photo-NMR Study on Light-Induced Nitrosyl Linkage Isomers Uncovers Their Structural, Electronic, and Diamagnetic Nature

A light-induced linkage NO isomer (MS1) in trans-[Ru(¹⁵NO)(py)₄¹⁹F](ClO₄)₂ is detected and measured for the first time by solid-state MAS NMR. Chemical shift tensors of ¹⁵N and ¹⁹F, along with ⁿJ(¹⁵N-¹⁹F) spin-spin couplings and T₁ relaxation times of MS1, are compared with the ground state (GS) at temperatures T < 250 K. Isotropic chemical shifts (¹⁵N and ¹⁹F) are well resolved for two crystallographically independent cations (A and B) [Ru(¹⁵NO)(py)₄¹⁹F]²⁺, allowing to define separately both populations of MS1 isomers and thermal decay rates for two structural sites. The relaxation times T₁ of ¹⁹F in the case of GS (30/38.6 s for sites A/B) and MS1 (11.6/11.8 s for sites A/B) indicate that both isomers are diamagnetic, which is the first experimental evidence of diamagnetic properties of MS1 in ruthenium nitrosyl. After light irradiation (λ = 420 nm), the NO ligand rotates by nearly 180° from F-Ru-N-O to F-Ru-O-N, whereby the isotropic chemical shifts of δ_iso(¹⁵N) increase and those of δ_iso(¹⁹F) decrease. The ⁿJ(¹⁵N-¹⁹F) couplings increase from ²J(¹⁵N-Ru-¹⁹F)_GS = 71 Hz to ³J(¹⁵N-O-Ru-¹⁹F)_MS1 = 105 Hz. These results are interpreted on the basis of DFT-CASTEP calculations including Bader-, Mulliken-, and Hirshfeld-charge density distributions of both states.

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.

Probing 29Si-17O connectivities and proximities by solid-state NMR

The measurement of dipolar and J- couplings between ²⁹Si and ¹⁷O isotopes is challenging owing to (i) the low abundance of both isotopes and (ii) their close Larmor frequencies, which only differ by 19%. These issues are circumvented here by the use of isotopic enrichment and dedicated triple-resonance magic-angle spinning NMR probe. The surface of ²⁹Si-enriched silica was labelled with ¹⁷O isotope and heated at 80 and 200 °C. ²⁹Si-¹⁷O connectivities and proximities were probed using two-dimensional (2D) through-bond and through-space heteronuclear multiple-quantum coherences (J- and D-HMQC) experiments between ¹⁷O and ²⁹Si nuclei. The simulation of the build-up of the J- and D-HMQC signals allowed the first experimental measurement of J- and dipolar coupling constants between ¹⁷O and ²⁹Si nuclei. These HMQC experiments allow distinguishing two distinct siloxane (SiOSi) oxygen sites: (i) those covalently bonded to Q³ and Q⁴ groups, having a hydroxyl group as a second neighbour and (ii) those covalently bonded to two Q⁴ groups. The measured J- and dipolar coupling constants of siloxane ¹⁷O nucleus with Q^{4 29}Si nuclei differ from those with Q^{3 29}Si nuclei. These results indicate that the ²⁹Si-¹⁷O one-bond J-coupling and Si-O bond length depend on the second neighbours of the Si atoms.

31P MAS NMR and DFT study of crystalline phosphate matrices

We present the study of the phosphorus local environment by using ³¹P MAS NMR in a series of seven double monophosphates M^IIM^IV(PO₄)₂ (M^II and M^IV being divalent and tetravalent cations, respectively) of yavapaiite and low-yavapaiite type crystal structures. Solid-state and cluster DFT calculations were found to be efficient for predicting the ³¹P isotropic chemical shift and chemical shift anisotropy. To achieve this performance, however, a proper computational optimisation of the experimental structural data was required. From the three optimisation methods tested, the full optimisation provided the best reference structure for the calculation of the NMR parameters of the studied phosphates. Also, a better prediction of the chemical shifts was possible by using a correction to the GIPAW calculated shielding.

Weak Intermolecular CH···N Hydrogen Bonding: Determination of 13CH-15N Hydrogen-Bond Mediated J Couplings by Solid-State NMR Spectroscopy and First-Principles Calculations

Weak hydrogen bonds are increasingly hypothesized to play key roles in a wide range of chemistry from catalysis to gelation to polymer structure. Here, ¹⁵N/¹³C spin-echo magic-angle spinning (MAS) solid-state nuclear magnetic resonance (NMR) experiments are applied to "view" intermolecular CH···N hydrogen bonding in two selectively labeled organic compounds, 4-[¹⁵N] cyano-4'-[¹³C₂] ethynylbiphenyl (1) and [¹⁵N₃,¹³C₆]-2,4,6-triethynyl-1,3,5-triazine (2). The synthesis of 2-¹⁵N₃,¹³C₆ is reported here for the first time via a multistep procedure, where the key element is the reaction of [¹⁵N₃]-2,4,6-trichloro-1,3,5-triazine (5) with [¹³C₂]-[(trimethylsilyl)ethynyl]zinc chloride (8) to afford its immediate precursor [¹⁵N₃,¹³C₆]-2,4,6-tris[(trimethylsilyl)ethynyl]-1,3,5-triazine (9). Experimentally determined hydrogen-bond-mediated ^2hJ_CN couplings (4.7 ± 0.4 Hz (1) and 4.1 ± 0.3 Hz (2)) are compared with density functional theory (DFT) gauge-including projector augmented wave (GIPAW) calculations, whereby species-independent coupling values ^2hK_CN (29.0 × 10¹⁹ kg m^-2 s^-2 A^-2 (1) and 27.9 × 10¹⁹ kg m^-2 s^-2 A^-2 (2)) quantitatively demonstrate the J couplings for these "weak" CH···N hydrogen bonds to be of a similar magnitude to those for conventionally observed NH···O hydrogen-bonding interactions in uracil (^2hK_NO: 28.1 and 36.8 × 10¹⁹ kg m^-2 s^-2 A^-2). Moreover, the GIPAW calculations show a clear correlation between increasing ^2hJ_CN (and ^3hJ_CN) coupling and reducing C(H)···N and H···N hydrogen-bonding distances, with the Fermi contact term accounting for at least 98% of the isotropic ^2hJ_CN coupling.

121/123Sb Nuclear Quadrupole Resonance Spectroscopy: Characterization of Non-Covalent Pnictogen Bonds and NQR Crystallography

Pnictogen (or pnicogen) bonding is an attractive interaction between the electrophilic region of group 15 elements (N, P, As, Sb, Bi) and a nucleophile. This interaction for which unique applications in catalysis have recently been uncovered continues to gain popularity. Here, we investigate a series of pnictogen-bonded cocrystals based on SbF₃ and SbCl₃, prepared via mechanochemical ball milling, with ^121/123Sb ( I = ⁵/₂ and ⁷/₂, respectively) nuclear quadrupole resonance (NQR) spectroscopy. Observed NQR frequency shifts upon cocrystallization are on the order of 0.1 to 10 MHz and are clearly diagnostic of the formation of pnictogen bonds to antimony. Further evidence for pnictogen bonding is obtained by complementary ¹³C cross-polarization magic-angle spinning solid-state NMR experiments. DFT calculations of NMR parameters as well as natural localized molecular orbital analyses support the experimental findings and elucidate the electronic origins of the experimental NQR frequency shifts. This work provides insights into the changes in the antimony quadrupolar coupling constant upon pnictogen bonding: strikingly, the decreases noted here parallel those known for hydrogen bonds, but contrast with the increases reported for halogen bonds. The utility of the observed antimony nuclear quadrupolar coupling constants in constraining structural models of cocrystals for which diffraction-based structures are unavailable, i.e., a rudimentary implementation of NQR crystallography, is established. Overall, this work offers a new approach to understand emerging classes of electrophilic interactions and to contextualize them in the broader landscape of established chemical bonding paradigms.

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.

29Si NMR Chemical Shifts in Crystalline and Amorphous Silicon Nitrides

We investigate ²⁹Si nuclear magnetic resonance (NMR) chemical shifts, δ_iso, of silicon nitride. Our goal is to relate the local structure to the NMR signal and, thus, provide the means to extract more information from the experimental ²⁹Si NMR spectra in this family of compounds. We apply structural modeling and the gauge-included projector augmented wave (GIPAW) method within density functional theory (DFT) calculations. Our models comprise known and hypothetical crystalline Si₃N₄, as well as amorphous Si₃N₄ structures. We find good agreement with available experimental ²⁹Si NMR data for tetrahedral Si^[4] and octahedral Si^[6] in crystalline Si₃N₄, predict the chemical shift of a trigonal-bipyramidal Si^[5] to be about −120 ppm, and quantify the impact of Si-N bond lengths on ²⁹Si δ_iso. We show through computations that experimental ²⁹Si NMR data indicates that silicon dicarbodiimide, Si(NCN)₂ exhibits bent Si-N-C units with angles of about 143° in its structure. A detailed investigation of amorphous silicon nitride shows that an observed peak asymmetry relates to the proximity of a fifth N neighbor in non-bonding distance between 2.5 and 2.8 Å to Si. We reveal the impact of both Si-N(H)-Si bond angle and Si-N bond length on ²⁹Si δ_iso in hydrogenated silicon nitride structure, silicon diimide Si(NH)₂.

Solid-state NMR study on changes of phosphate and proton species in metal pyrophosphate composite (MP2 O7 -MO2 ) ceramics

In order to investigate the effects of composite production and the role of the counter cation for metal phosphate conductors, changes in the solid-state nuclear magnetic resonance (NMR) spectra and magnetic relaxation times caused by the removal of volatile products and water washing were examined for various metal pyrophosphate (MP₂ O₇ -MO₂ ; M = Sn, Si, Ti, and Zr). Acidic species could be detected by ¹ H MAS NMR spectra for all the composites except ZrP₂ O₇ -ZrO₂ that had the lowest conductivity. The ³¹ P DD-MAS NMR spectra for MP₂ O₇ -MO₂ composites showed different signal patterns depending on the counter cations participating in the ion exchange as a result of different microstructures. Combinational analysis of ³¹ P DD-MAS and ³¹ P CP-MAS NMR spectra of the composites indicated that protonic bulk phosphates were observed at slightly lower fields than non-protonic bulk phosphates in all of the MP₂ O₇ -MO₂ composites. After water washing, the acidic species and the protonic bulk phosphates of MP₂ O₇ -MO₂ composites disappeared or were reduced to trace amounts. The T₁ H values of the water-washed composites lengthened because of removal of orthophosphoric acid, although the T₁ P values remained almost unchanged. The results of the solid-state NMR studies suggest that the protonic bulk phosphates of MP₂ O₇ -MO₂ composites do not generally distribute in the bulk but exist in the interface between excess H₃ PO₄ and the bulk. Copyright © 2016 John Wiley & Sons, Ltd.

14N NQR lineshape in nanocrystals: An ab initio investigation of urea

¹⁴N nuclear quadrupole resonance (NQR) lineshapes mostly contain information of low interest, although in nanocrystals they may display some unexpected behaviour. In this work, we present an ab initio computational study of the ¹⁴N NQR lineshapes in urea nanocrystals as a function of the nanocrystal size and geometry, focusing on the surface induced broadening of the lineshapes. The lineshapes were obtained through a calculation of the electric field gradient for each nitrogen site in the nanocrystal separately, taking into account the individual crystal field by embedding the molecule of interest in a suitable lattice of point multipoles representing other urea molecules in the nanocrystal. The small influence of distant molecules is found with a series expansion, using the in-crystal Sternheimer shieldings which we also calculated ab initio. We have considered nanocrystals with two geometries: a sphere and a cube, with characteristic sizes between 5 and 100 nm. Our calculations suggest that there is a dramatic difference between the linewidths for the two geometries. For spheres, we find a steep drop in linewidths at ∼10 nm; at 5 nm the linewidth is ∼11 kHz, whereas for sizes above 20 nm the linewidth is practically negligible (<100 Hz). For cubes, on the other hand, we find a steady 1/size decrease, from 12 kHz at 10 nm to 1.2 kHz at 100 nm. This analysis is important for ¹⁴N NQR spectroscopy of crystalline pharmaceuticals, where nanoparticles are increasingly more often embedded in some sort of matrix. Although this is only a theoretical analysis, we believe that this work can serve as a guidance for the forthcoming experimental analysis.

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

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

Temperature dependence of X-ray absorption and nuclear magnetic resonance spectra: probing quantum vibrations of light elements in oxides

Probing the quantum thermal fluctuations of nuclei in light-element oxides using XANES and NMR spectroscopies.

From crystalline to amorphous calcium pyrophosphates: A solid state Nuclear Magnetic Resonance perspective

Hydrated calcium pyrophosphates (CPP, Ca2P2O7·nH2O) are a fundamental family of materials among osteoarticular pathologic calcifications. In this contribution, a comprehensive multinuclear NMR (Nuclear Magnetic Resonance) study of four crystalline and two amorphous phases of this family is presented. (1)H, (31)P and (43)Ca MAS (Magic Angle Spinning) NMR spectra were recorded, leading to informative fingerprints characterizing each compound. In particular, different (1)H and (43)Ca solid state NMR signatures were observed for the amorphous phases, depending on the synthetic procedure used. The NMR parameters of the crystalline phases were determined using the GIPAW (Gauge Including Projected Augmented Wave) DFT approach, based on first-principles calculations. In some cases, relaxed structures were found to improve the agreement between experimental and calculated values, demonstrating the importance of proton positions and pyrophosphate local geometry in this particular NMR crystallography approach. Such calculations serve as a basis for the future ab initio modeling of the amorphous CPP phases.

STATEMENT OF SIGNIFICANCE

The general concept of NMR crystallography is applied to the detailed study of calcium pyrophosphates (CPP), whether hydrated or not, and whether crystalline or amorphous. CPP are a fundamental family of materials among osteoarticular pathologic calcifications. Their prevalence increases with age, impacting on 17.5% of the population after the age of 80. They are frequently involved or associated with acute articular arthritis such as pseudogout. Current treatments are mainly directed at relieving the symptoms of joint inflammation but not at inhibiting CPP formation nor at dissolving these crystals. The combination of advanced NMR techniques, modeling and DFT based calculation of NMR parameters allows new original insights in the detailed structural description of this important class of biomaterials.

Selective measurements of long-range homonuclear J-couplings in solid-state NMR

We demonstrate here that the principle of frequency-selective spin-echoes can be extended to the measurements of long-range homonuclear scalar J-couplings in the solid-state. Singly or doubly frequency-selective pulses were used to generate either a J-modulated experiment (S) or a reference experiment (S0). The combination of these two distinct experiments provides experimental data that, in favorable cases, are insensitive to incoherent relaxation effects, and which can be used to estimate long-range homonuclear J-couplings in multiple spin-systems. The concept is illustrated in the case of a uniformly (13)C and (15)N labeled sample of L-histidine, where the absolute value of homonuclear J-couplings between two spins separated by one, two or three covalent bonds are measured. Moreover, we show that a (2)J((15)N-C-(15)N) coupling as small as 0.9 Hz can be precisely measured with the method presented here.

Investigation of the interface in silica-encapsulated liposomes by combining solid state NMR and first principles calculations

In the context of nanomedicine, liposils (liposomes and silica) have a strong potential for drug storage and release schemes: such materials combine the intrinsic properties of liposome (encapsulation) and silica (increased rigidity, protective coating, pH degradability). In this work, an original approach combining solid state NMR, molecular dynamics, first principles geometry optimization, and NMR parameters calculation allows the building of a precise representation of the organic/inorganic interface in liposils. {(1)H-(29)Si}(1)H and {(1)H-(31)P}(1)H Double Cross-Polarization (CP) MAS NMR experiments were implemented in order to explore the proton chemical environments around the silica and the phospholipids, respectively. Using VASP (Vienna Ab Initio Simulation Package), DFT calculations including molecular dynamics, and geometry optimization lead to the determination of energetically favorable configurations of a DPPC (dipalmitoylphosphatidylcholine) headgroup adsorbed onto a hydroxylated silica surface that corresponds to a realistic model of an amorphous silica slab. These data combined with first principles NMR parameters calculations by GIPAW (Gauge Included Projected Augmented Wave) show that the phosphate moieties are not directly interacting with silanols. The stabilization of the interface is achieved through the presence of water molecules located in-between the head groups of the phospholipids and the silica surface forming an interfacial H-bonded water layer. A detailed study of the (31)P chemical shift anisotropy (CSA) parameters allows us to interpret the local dynamics of DPPC in liposils. Finally, the VASP/solid state NMR/GIPAW combined approach can be extended to a large variety of organic-inorganic hybrid interfaces.

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

In 2001, Mauri and Pickard introduced the gauge including projected augmented wave (GIPAW) method that enabled for the first time the calculation of all-electron NMR parameters in solids, i.e. accounting for periodic boundary conditions. The GIPAW method roots in the plane wave pseudopotential formalism of the density functional theory (DFT), and avoids the use of the cluster approximation. This method has undoubtedly revitalized the interest in quantum chemical calculations in the solid-state NMR community. It has quickly evolved and improved so that the calculation of the key components of NMR interactions, namely the shielding and electric field gradient tensors, has now become a routine for most of the common nuclei studied in NMR. Availability of reliable implementations in several software packages (CASTEP, Quantum Espresso, PARATEC) make its usage more and more increasingly popular, maybe indispensable in near future for all material NMR studies. The majority of nuclei of the periodic table have already been investigated by GIPAW, and because of its high accuracy it is quickly becoming an essential tool for interpreting and understanding experimental NMR spectra, providing reliable assignments of the observed resonances to crystallographic sites or enabling a priori prediction of NMR data. The continuous increase of computing power makes ever larger (and thus more realistic) systems amenable to first-principles analysis. In the near future perspectives, as the incorporation of dynamical effects and/or disorder are still at their early developments, these areas will certainly be the prime target.