Mapping of N-C Bond Formation from a Series of Crystalline Peri-Substituted Naphthalenes by Charge Density and Solid-State NMR Methodologies.

A combination of charge density studies and solid state nuclear magnetic resonance (NMR) 1 JNC coupling measurements supported by periodic density functional theory (DFT) calculations is used to characterise the transition from an n-π* interaction to bond formation between a nucleophilic nitrogen atom and an electrophilic sp2 carbon atom in a series of crystalline peri-substituted naphthalenes. As the N⋅⋅⋅C distance reduces there is a sharp decrease in the Laplacian derived from increasing charge density between the two groups at ca. N⋅⋅⋅C = 1.8 Å, with the periodic DFT calculations predicting, and heteronuclear spin-echo NMR measurements confirming, the 1 JNC couplings of ≈3-6 Hz for long C-N bonds (1.60-1.65 Å), and 1 JNC couplings of <1 Hz for N⋅⋅⋅C >2.1 Å.

Measuring multiple 17O-13C J-couplings in naphthalaldehydic acid: a combined solid state NMR and density functional theory approach.

A combined multinuclear solid state NMR and gauge included projected augmented wave, density functional theory (GIPAW DFT) computational approach is evaluated to determine the four heteronuclear 1J(13C,17O) couplings in solid 17O enriched naphthalaldehydic acid. Direct multi-field 17O magic angle spinning (MAS), triple quantum MAS (3QMAS) and double rotation (DOR) experiments are initially utilised to evaluate the accuracy of the DFT approximations used in the calculation of the isotropic chemical shifts (δiso), quadrupole coupling constants (CQ) and asymmetry (ηQ) parameters. These combined approaches give δiso values of 313, 200 and 66 ppm for the carbonyl (C[double bond, length as m-dash]O), ether (-O-) and hydroxyl (-OH) environments, respectively, with the corresponding measured quadrupole products (PQ) being 8.2, 9.0 and 10.6 MHz. The geometry optimised DFT structure derived using the CASTEP code gives firm agreement with the shifts observed for the ether (δiso = 223, PQ = 9.4 MHz) and hydroxyl (δiso = 62, PQ = 10.5 MHz) environments but the unoptimised experimental XRD structure has better agreement for the carbonyl group (δiso = 320, PQ = 8.3 MHz). The determined δiso and ηQ values are shown to be consistent with bond lengths closer to 1.222 Å (experimental length) rather than the geometry optimised length of 1.238 Å. The geometry optimised DFT 1J(13C,17O) coupling to the hydroxyl is calculated as 20 Hz and the couplings to the ether were calculated to be 37 (O-C[double bond, length as m-dash]O) and 32 (O-C-OH) Hz. The scalar coupling parameters for the unoptimised experimental carbonyl group predict a 1J(13C,17O) value of 28 Hz, whilst optimisation gives a value of 27 Hz. These calculated 1J(13C,17O) couplings, together with estimations of the probability of each O environment being isotopically labelled (determined by electrospray ionisation mass spectrometry) and the measured refocussable transverse dephasing (T2') behaviour, are combined to simulate the experimental decay behaviour. Good agreement between the measured and calculated decay behaviour is observed.

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, 15N/13C 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-[15N] cyano-4'-[13C2] ethynylbiphenyl (1) and [15N3,13C6]-2,4,6-triethynyl-1,3,5-triazine (2). The synthesis of 2-15N3,13C6 is reported here for the first time via a multistep procedure, where the key element is the reaction of [15N3]-2,4,6-trichloro-1,3,5-triazine (5) with [13C2]-[(trimethylsilyl)ethynyl]zinc chloride (8) to afford its immediate precursor [15N3,13C6]-2,4,6-tris[(trimethylsilyl)ethynyl]-1,3,5-triazine (9). Experimentally determined hydrogen-bond-mediated 2hJCN 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 2hKCN (29.0 × 1019 kg m-2 s-2 A-2 (1) and 27.9 × 1019 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 (2hKNO: 28.1 and 36.8 × 1019 kg m-2 s-2 A-2). Moreover, the GIPAW calculations show a clear correlation between increasing 2hJCN (and 3hJCN) 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 2hJCN coupling.

Regularized SCAN functional.

We propose modifications to the functional form of the Strongly Constrained and Appropriately Normed (SCAN) density functional to eliminate numerical instabilities. This is necessary to allow reliable, automatic generation of pseudopotentials (including projector augmented-wave potentials). The regularized SCAN is designed to match the original form very closely, and we show that its performance remains comparable.

An NMR crystallography investigation of furosemide.

This paper presents an NMR crystallography study of three polymorphs of furosemide. Experimental magic-angle spinning (MAS) solid-state NMR spectra are reported for form I of furosemide, and these are assigned using density-functional theory (DFT)-based gauge-including projector augmented wave (GIPAW) calculations. Focusing on the three known polymorphs, we examine the changes to the NMR parameters due to crystal packing effects. We use a recently developed formalism to visualise which regions are responsible for the chemical shielding of particular sites and hence understand the variation in NMR parameters between the three polymorphs.

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

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

Tin chemical shift anisotropy in tin dioxide: On ambiguity of CSA asymmetry derived from MAS spectra.

Two different axial symmetries of the 119Sn chemical shift anisotropy (CSA) in tin dioxide with the asymmetry parameter (η) of 0 and 0.27 were reported previously based on the analysis of MAS NMR spectra. By analyzing the static powder pattern, we show that the 119Sn CSA is axially symmetric. A nearly axial symmetry and the principal axis system of the 119Sn chemical shift tensor in SnO2 were deduced from periodic scalar-relativistic density functional theory (DFT) calculations of NMR parameters. The implications of fast small-angle motions on CSA parameters were also considered, which could potentially lead to a CSA symmetry in disagreement with a crystal symmetry. Our analysis of experimental spectra using spectral simulations and iterative fittings showed that MAS spectra recorded at relatively high frequencies do not show sufficiently distinct features in order to distinguish CSAs with η ≈ 0 and η ≈ 0.4. The example of SnO2 shows that both the MAS lineshape and spinning sideband analyses may overestimate the η value by as much as ∼0.3 and ∼0.4, respectively. The results confirm that a static powder pattern must be analysed in order to improve the accuracy of the CSA asymmetry measurements. The measurements on SnO2 nanoparticles showed that the asymmetry parameter of the 119Sn CSA increases for nm-sized particles with a larger surface area compared to µm-sized particles. The increase of the η value for tin atoms near the surface in SnO2 was also confirmed by DFT calculations.

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

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

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 (CQ) 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.

Visualising crystal packing interactions in solid-state NMR: Concepts and applications.

In this article, we introduce and apply a methodology, based on density functional theory and the gauge-including projector augmented wave approach, to explore the effects of packing interactions on solid-state nuclear magnetic resonance (NMR) parameters. A visual map derived from a so-termed "magnetic shielding contribution field" can be made of the contributions to the magnetic shielding of a specific site-partitioning the chemical shift to specific interactions. The relation to the established approaches of examining the molecule to crystal change in the chemical shift and the nuclear independent chemical shift is established. The results are applied to a large sample of 71 molecular crystals and three further specific examples from supermolecular chemistry and pharmaceuticals. This approach extends the NMR crystallography toolkit and provides insight into the development of both cluster based approaches to the predictions of chemical shifts and for empirical predictions of chemical shifts in solids.

Investigating Unusual Homonuclear Intermolecular "Through-Space" J Couplings in Organochalcogen Systems.

Although the electron-mediated spin-spin or J coupling is conventionally viewed as transmitted via covalent bonds, examples of J couplings between atoms that are not formally bonded but are in close proximity (termed "through-space" J couplings) have been reported. In this work, we investigate the observation of homonuclear 31P J couplings in organochalcogen heterocycles, which occur between 31P in two separate molecules, confirming without doubt their through-space nature. The presence of this interaction is even more surprising for one compound, where it occurs between crystallographically equivalent species. Although crystallographically equivalent species need not be magnetically equivalent in the solid state, owing to the presence of anisotropic interactions, we demonstrate that it is not the shielding anisotropy that lifts magnetic equivalence, in this case, but the presence of heteronuclear couplings to 77Se. We support our experimental observations with periodic scalar-relativistic density functional theory calculations and coupling density deformation plots to visualize the mechanism of these interesting interactions.

Visualization and processing of computed solid-state NMR parameters: MagresView and MagresPython.

We introduce two open source tools to aid the processing and visualisation of ab-initio computed solid-state NMR parameters. The Magres file format for computed NMR parameters (as implemented in CASTEP v8.0 and QuantumEspresso v5.0.0) is implemented. MagresView is built upon the widely used Jmol crystal viewer, and provides an intuitive environment to display computed NMR parameters. It can provide simple pictorial representation of one- and two-dimensional NMR spectra as well as output a selected spin-system for exact simulations with dedicated spin-dynamics software. MagresPython provides a simple scripting environment to manipulate large numbers of computed NMR parameters to search for structural correlations.

Unusual intermolecular "through-space" j couplings in p-se heterocycles.

Solid-state NMR spectra of new P-Se heterocycles based on peri-substituted naphthalene motifs show the presence of unusual J couplings between Se and P. These couplings are between atoms in adjacent molecules and occur "through space", rather than through conventional covalent bonds. Experimental measurements are supported by relativistic DFT calculations, which confirm the presence of couplings between nonbonded atoms, and provide information on the pathway of the interaction. This observation improves the understanding of J couplings and offers insight into the factors that affect crystal packing in solids, for future synthetic exploitation.

Relativistic nuclear magnetic resonance J-coupling with ultrasoft pseudopotentials and the zeroth-order regular approximation.

We present a method for the first-principles calculation of nuclear magnetic resonance (NMR) J-coupling in extended systems using state-of-the-art ultrasoft pseudopotentials and including scalar-relativistic effects. The use of ultrasoft pseudopotentials is allowed by extending the projector augmented wave (PAW) method of Joyce et al. [J. Chem. Phys. 127, 204107 (2007)]. We benchmark it against existing local-orbital quantum chemical calculations and experiments for small molecules containing light elements, with good agreement. Scalar-relativistic effects are included at the zeroth-order regular approximation level of theory and benchmarked against existing local-orbital quantum chemical calculations and experiments for a number of small molecules containing the heavy row six elements W, Pt, Hg, Tl, and Pb, with good agreement. Finally, (1)J(P-Ag) and (2)J(P-Ag-P) couplings are calculated in some larger molecular crystals and compared against solid-state NMR experiments. Some remarks are also made as to improving the numerical stability of dipole perturbations using PAW.

First-Principles Band Alignments at the Si:Anatase TiO2 Interface.

TiO2 has been identified as a promising electron transport layer in Si solar cells. Experiments have revealed that the Si:TiO2 interface undergoes structural changes depending on how it was fabricated. However, less is understood about the sensitivity of electronic properties, such as band alignments, to these changes. Here, we present first-principles calculations of band alignments between Si and anatase TiO2, investigating different surface orientations and terminations. By calculating vacuum-level alignments, we observe a large band offset reduction of 2.5 eV for the O-terminated Si slab compared to other terminations. Furthermore, a 0.5 eV increase is found for the anatase (101) surface compared to (001). We compare the band offsets obtained through vacuum alignment with four different heterostructure models. Even though the heterostructure models contain an excess of oxygen, their offsets agree well with vacuum-level alignments using stoichiometric or H-terminated slabs, and the reduction in band offsets seen for the O-terminated Si slab is not observed. Additionally, we have investigated different exchange-correlation treatments including PBE + U, postprocessing GW corrections, and the meta-GGA rSCAN functional. We find that rSCAN provides more accurate band offsets than PBE, but further corrections are still required to achieve <0.5 eV accuracy. Overall, our study quantifies the importance of surface termination and orientation for this interface.

Ultra-high resolution 17O solid-state NMR spectroscopy of biomolecules: a comprehensive spectral analysis of monosodium L-glutamate·monohydrate.

Monosodium L-glutamate monohydrate, a multiple oxygen site (eight) compound, is used to demonstrate that a combination of high-resolution solid-state NMR spectroscopic techniques opens up new possibilities for (17)O as a nuclear probe of biomolecules. Eight oxygen sites have been resolved by double rotation (DOR) and multiple quantum (MQ) NMR experiments, despite the (17)O chemical shifts lying within a narrow shift range of <50 ppm. (17)O DOR NMR not only provides high sensitivity and spectral resolution, but also allows a complete set of the NMR parameters (chemical shift anisotropy and electric-field gradient) to be determined from the DOR spinning-sideband manifold. These (17)O NMR parameters provide an important multi-parameter comparison with the results from the quantum chemical NMR calculations, and enable unambiguous oxygen-site assignment and allow the hydrogen positions to be refined in the crystal lattice. The difference in sensitivity between DOR and MQ NMR experiments of oxygen in bio/organic molecules is also discussed. The data presented here clearly illustrates that a high resolution (17)O solid-state NMR methodology is now available for the study of biomolecules, offering new opportunities for resolving structural information and hence new molecular insights.

Towards homonuclear J solid-state NMR correlation experiments for half-integer quadrupolar nuclei: experimental and simulated 11B MAS spin-echo dephasing and calculated 2J(BB) coupling constants for lithium diborate.

Magic-angle spinning (MAS) NMR spin-echo dephasing is systematically investigated for the spin I = 3/2 (11)B nucleus in lithium diborate, Li(2)O·2B(2)O(3). A clear dependence on the quadrupolar frequency (ω(Q)(PAS)/2π = 3C(Q)/[4I(2I- 1)]) is observed: the B3 (larger C(Q)) site dephases more slowly than the B4 site at all investigated MAS frequencies (5 to 20 kHz) at 14.1 T. Increasing the MAS frequency leads to markedly slower dephasing for the B3 site, while there is a much less evident effect for the B4 site. Considering samples at 5, 25, 80 (natural abundance) and 100% (11)B isotopic abundance, dephasing becomes faster for both sites as the (11)B isotopic abundance increases. The experimental behaviour is rationalised using density matrix simulations for two and three dipolar-coupled (11)B nuclei. The experimentally observed slower dephasing for the larger C(Q) (B3) site is reproduced in all simulations and is explained by the reintroduction of the dipolar coupling by the so-called "spontaneous quadrupolar-driven recoupling mechanism" having a different dependence on the MAS frequency for different quadrupolar frequencies. Specifically, isolated spin-pair simulations show that the spontaneous quadrupolar-driven recoupling mechanism is most efficient when the quadrupolar frequency is equal to twice the MAS frequency. While for isolated spin-pair simulations, increasing the MAS frequency leads to faster dephasing, agreement with experiment is observed for three-spin simulations which additionally include the homogeneous nature of the homonuclear dipolar coupling network. First-principles calculations, using the GIPAW approach, of the (2)J(11B-11B) couplings in lithium diborate, metaborate and triborate are presented: a clear trend is revealed whereby the (2)J(11B-11B) couplings increase with increasing B-O-B bond angle and B-B distance. However, the calculated (2)J(11B-11B) couplings are small (0.95, 1.20 and 2.65 Hz in lithium diborate), thus explaining why no zero crossing due to J modulation is observed experimentally, even for the sample at 25% (11)B where significant spin-echo intensity remains out to durations of ∼200 ms.

High-resolution (19)F MAS NMR spectroscopy: structural disorder and unusual J couplings in a fluorinated hydroxy-silicate.

High-resolution (19)F magic angle spinning (MAS) NMR spectroscopy is used to study disorder and bonding in a crystalline solid. (19)F MAS NMR reveals four distinct F sites in a 50% fluorine-substituted deuterated hydrous magnesium silicate (clinohumite, 4Mg(2)SiO(4)·Mg(OD(1-x)F(x))(2) with x = 0.5), indicating extensive structural disorder. The four (19)F peaks can be assigned using density functional theory (DFT) calculations of NMR parameters for a number of structural models with a range of possible local F environments generated by F(-)/OH(-) substitution. These assignments are supported by two-dimensional (19)F double-quantum MAS NMR experiments that correlate F sites based on either spatial proximity (via dipolar couplings) or through-bond connectivity (via scalar, or J, couplings). The observation of (19)F-(19)F J couplings is unexpected as the fluorines coordinate Mg atoms and the Mg-F interaction is normally considered to be ionic in character (i.e., there is no formal F-Mg-F covalent bonding arrangement). However, DFT calculations predict significant (19)F-(19)F J couplings, and these are in good agreement with the splittings observed in a (19)F J-resolved MAS NMR experiment. The existence of these J couplings is discussed in relation to both the nature of bonding in the solid state and the occurrence of so-called "through-space" (19)F-(19)F J couplings in solution. Finally, we note that we have found similar structural disorder and spin-spin interactions in both synthetic and naturally occurring clinohumite samples.

Time averaging of NMR chemical shifts in the MLF peptide in the solid state.

Since experimental measurements of NMR chemical shifts provide time and ensemble averaged values, we investigated how these effects should be included when chemical shifts are computed using density functional theory (DFT). We measured the chemical shifts of the N-formyl-L-methionyl-L-leucyl-L-phenylalanine-OMe (MLF) peptide in the solid state, and then used the X-ray structure to calculate the (13)C chemical shifts using the gauge including projector augmented wave (GIPAW) method, which accounts for the periodic nature of the crystal structure, obtaining an overall accuracy of 4.2 ppm. In order to understand the origin of the difference between experimental and calculated chemical shifts, we carried out first-principles molecular dynamics simulations to characterize the molecular motion of the MLF peptide on the picosecond time scale. We found that (13)C chemical shifts experience very rapid fluctuations of more than 20 ppm that are averaged out over less than 200 fs. Taking account of these fluctuations in the calculation of the chemical shifts resulted in an accuracy of 3.3 ppm. To investigate the effects of averaging over longer time scales we sampled the rotameric states populated by the MLF peptides in the solid state by performing a total of 5 micros classical molecular dynamics simulations. By averaging the chemical shifts over these rotameric states, we increased the accuracy of the chemical shift calculations to 3.0 ppm, with less than 1 ppm error in 10 out of 22 cases. These results suggests that better DFT-based predictions of chemical shifts of peptides and proteins will be achieved by developing improved computational strategies capable of taking into account the averaging process up to the millisecond time scale on which the chemical shift measurements report.