Search

A disaccharide is a challenging case for high-resolution (1)H solid-state NMR because of the 24 distinct protons (14 aliphatic and 10 OH) having (1)H chemical shifts that all fall within a narrow range of approximately 3 to 7 ppm. High-resolution (1)H (500 MHz) double-quantum (DQ) combined rotation and multiple pulse sequence (CRAMPS) solid-state NMR spectra of beta-maltose monohydrate are presented. (1)H-(1)H DQ-SQ CRAMPS spectra are presented together with (1)H (DQ)-(13)C correlation spectra obtained with a new pulse sequence that correlates a high-resolution (1)H DQ dimension with a (13)C single quantum (SQ) dimension using the refocused INEPT pulse-sequence element to transfer magnetization via one-bond (13)C-(1)H J couplings. Compared to the observation of only a single broad peak in a (1)H DQ spectrum recorded at 30 kHz magic-angle spinning (MAS), the use of DUMBO (1)H homonuclear decoupling in the (1)H DQ CRAMPS experiment allows the resolution of distinct DQ correlation peaks which, in combination with first-principles chemical shift calculations based on the GIPAW (Gauge Including Projector Augmented Waves) plane-wave pseudopotential approach, enables the assignment of the (1)H resonances to the 24 distinct protons. We believe this to be the first experimental solid-state NMR determination of the hydroxyl OH (1)H chemical shifts for a simple sugar. Variable-temperature (1)H-(1)H DQ CRAMPS spectra reveal small increases in the (1)H chemical shifts of the OH resonances upon decreasing the temperature from 348 K to 248 K.

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.

The NMR parameters obtained from solid-state DFT calculations within the GIPAW approach for (17)O- and (69/71)Ga-sites in a range of predominantly oxide-based (group II monoxides, SrTiO(3), BaZrO(3), BaSnO(3), BaTiO(3), LaAlO(3), LaGaO(3), SrZrO(3), MgSiO(3) and Ba(2)In(2)O(5)), and gallate (alpha- and beta-Ga(2)O(3), LiGaO(2), NaGaO(2), GaPO(4) and LaGaO(3)) materials are compared with experimental values, with a view to the future application of a similar approach to doped phases of interest as candidate intermediate temperature solid oxide fuel cell (ITSOFC) electrolytes. Isotropic and anisotropic chemical shift parameters, quadrupolar coupling constants, and associated asymmetries are presented and analyzed. The unusual GaO(5) site occurring in LaGaGe(2)O(7) is also fully characterised. In general, it is found that the theoretical results closely track the experimental trends, though some deviations are identified and discussed, particularly in regard to quadrupolar eta(Q)-values. The high quality of the computed results suggests that this approach can be extended to study more complex and disordered phases., (Copyright 2010 Elsevier Inc. All rights reserved.)

To efficiently obtain multiple-quantum magic-angle spinning (MQMAS) spectra of the nuclide 45Sc (I=7/2), we have combined several previously suggested techniques to enhance the signal-to-noise ratio and to improve spectral resolution for the test sample, scandium sulphate pentahydrate (ScSPH). Whereas the 45Sc-3QMAS spectrum of ScSPH does not offer sufficient resolution to clearly distinguish between the 3 scandium sites present in the crystal structure, these sites are well-resolved in the 5QMAS spectrum. The loss of sensitivity incurred by using MQMAS with 5Q coherence order is partly compensated for by using fast-amplitude modulated (FAM) sequences to improve the efficiency of both 5Q coherence excitation and conversion. Also, heteronuclear decoupling is employed to minimise dephasing of the 45Sc signal during the 5Q evolution period due to dipolar couplings with the water protons in the ScSPH sample. Application of multi-pulse decoupling schemes such as TPPM and SPINAL results in improved sensitivity and resolution in the F(1) (isotropic) dimension of the 5QMAS spectrum, the best results being achieved with the recently suggested SW(f)-TPPM sequence. By numerical fitting of the 45Sc-NMR spectra of ScSPH from 3QMAS, 5QMAS and single-quantum MAS at magnetic fields B(0)=9.4 T and 17.6 T, the isotropic chemical shift delta(iso), the quadrupolar coupling constant chi, and the asymmetry parameter eta were obtained. Averaging over all experiments, the NMR parameters determined for the 3 scandium sites, designated (a), (b) and (c) are: delta(iso)(a)=-15.5+/-0.5 ppm, chi(a)=5.60+/-0.10 MHz, eta(a)=0.06+/-0.05; delta(iso)(b)=-12.9+/-0.5 ppm, chi(b)=4.50+/-0.10 MHz, eta(b)=1.00+/-0.00; and delta(iso)(c)=-4.7+/-0.2 ppm, chi(c)=4.55+/-0.05 MHz, eta(c)=0.50+/-0.02. The NMR scandium species were assigned to the independent crystallographic sites by evaluating their experimental response to proton decoupling, and by density functional theory (DFT) calculations using the PAW and GIPAW approaches, in the following way: Sc(1) to (c), Sc(2) to (a), and Sc(3) to (b). The need to compute NMR parameters using an energy-optimised crystal structure is once again demonstrated., (2009 Elsevier Inc. All rights reserved.)

A fast method for crystal structure determination using crystal structure prediction and solid-state (1)H NMR is presented. This technique does not need any prior knowledge except the chemical formula; resonance assignment is not necessary. Starting from an ensemble of predicted crystal structures for powdered thymol, comparison between experimental and calculated (1)H solid-state isotropic NMR chemical shifts is sufficient to determine which predicted structure corresponds to the powder under study. The same approach using proton-proton spin-diffusion data is successful and can be used for cross-validation.

The application of periodic density functional theory-based methods to the calculation of (95)Mo electric field gradient (EFG) and chemical shift (CS) tensors in solid-state molybdenum compounds is presented. Calculations of EFG tensors are performed using the projector augmented-wave (PAW) method. Comparison of the results with those obtained using the augmented plane wave + local orbitals (APW+lo) method and with available experimental values shows the reliability of the approach for (95)Mo EFG tensor calculation. CS tensors are calculated using the recently developed gauge-including projector augmented-wave (GIPAW) method. This work is the first application of the GIPAW method to a 4d transition-metal nucleus. The effects of ultra-soft pseudo-potential parameters, exchange-correlation functionals and structural parameters are precisely examined. Comparison with experimental results allows the validation of this computational formalism.

We have used density-functional-theory methods to study phases of iron at terapascal (TPa) pressures. Hexagonal-close-packed (hcp) iron is stable from below 0.1 TPa to multi-TPa pressures, where we find a window of 7-21 TPa in which the face-centred-cubic (fcc) phase is slightly more stable. At 34 TPa we find a transition to a body-centred-tetragonal (bct) phase, which is a small distortion of the body-centred-cubic (bcc) structure. The bcc phase also becomes more stable than the hcp above 35 TPa. The bct and bcc phases become considerably more stable than hcp and fcc at even higher pressures, and we expect this result to hold at high temperatures.

We describe the results of a neutron powder diffraction study of perdeuterated ammonia monohydrate (AMH, ND(3).D(2)O) carried out in the range 102

We have determined the crystal structure of ammonia monohydrate phase II (AMH II) employing a combination of ab initio computational structure prediction and structure solution from neutron powder diffraction data using direct space methods. Neutron powder diffraction data were collected from perdeuterated AMH II using the D2B high-resolution diffractometer at the Institut Laue-Langevin. AMH II crystallizes in space-group Pbca with 16 formula units in a unit-cell of dimensions a = 18.8285(4) A, b = 6.9415(2) A, c = 6.8449(2) A, and V = 894.61(3) A3 [rho(calc)(deuterated) = 1187.56(4) kg m(-3)] at 502 MPa, 180 K. The structure is characterized by sheets of tessellated pentagons formed by orientationally ordered O-D...O, O-D...N, and N-D...O hydrogen-bonds; these sheets are stacked along the a-axis and connected by N-D...O hydrogen bonds alone. With the exception of the simple body-centered-cubic high-pressure phases of ammonia monohydrate and ammonia dihydrate, this is the first complex molecular structure of any of the high-pressure stoichiometric ammonia hydrates to be determined. The powder structure solution is complemented by an ab initio structure prediction using density functional theory which gives an almost identical hydrogen bonding network.

The (13)C NMR chemical shifts of fluorinated semiconducting single-walled carbon nanotubes (SWNTs) were computed using a gaugeincluding projectoraugmented plane wave (GIPAW) density functional method. The chemical shifts of the fluorinated carbons (C(alpha)) were rather insensitive to the degree and pattern of functionalization, as well as to the nanotube radius. The calculated shifts were typically between 82 and 84 ppm, which is in excellent agreement with a recent experimental value of 83.5 ppm. Because of the insensitivity of the shifts to details of the nanotube's electronic structure and diameter, the NMR signals of the C(alpha) carbons are a useful indicator of the degree of functionalization in a heterogeneous bulk sample. At high degrees of functionalization, the shifts of atoms neighboring C(alpha) might also be useful indicators of the functionalization pattern.

A protocol for the structure determination of powdered solids at natural abundance by NMR is presented and illustrated for the case of the small drug molecule thymol. The procedure uses proton spin-diffusion data from two-dimensional NMR experiments in combination with periodic DFT refinements incorporating (1)H and (13)C NMR chemical shifts. For thymol, the method yields a crystal structure for the powdered sample, which differs by an atomic root-mean-square-deviation (all atoms except methyl group protons) of only 0.07 A from the single crystal X-ray diffraction structure with DFT-optimized proton positions.

Ab initio density-functional-theory calculations and a structure-searching technique are used to identify candidate high-pressure phases of lithium (Li). We predict threefold coordinated structures to be stable in the pressure range 40-450 GPa and fourfold structures at higher pressures. We describe these low-coordination phases as elemental electrides. All of the stable phases are metallic but the Cmca-24 structure, and two distortions of it which are marginally the most stable in the pressure range 86-106 GPa, are nearly semiconducting with densities of electronic states at the Fermi energy of only a few percent of the free-electron value.

Density-functional-theory calculations and a structure-searching method are used to identify candidate high-pressure phases of nitrogen. We find six structures which are calculated to be more stable than previously studied structures at some pressures. Our four new molecular structures give insight into the most efficient packings of nitrogen molecules at high pressures, and we predict two new nonmolecular structures to be stable at very high pressures.

NMR chemical shifts were calculated from first principles for well defined crystalline organic solids. These density functional theory calculations were carried out within the plane-wave pseudopotential framework, in which truly extended systems are implicitly considered. The influence of motional effects was assessed by averaging over vibrational modes or over snapshots taken from ab initio molecular dynamics simulations. It is observed that the zero-point correction to chemical shifts can be significant, and that thermal effects are particularly noticeable for shielding anisotropies and for a temperature-dependent chemical shift. This study provides insight into the development of highly accurate first principles calculations of chemical shifts in solids, highlighting the role of motional effects on well defined systems.

A plane wave based method for the calculation of core-level spectra is presented. We provide details of the implementation of the method in the pseudopotential density functional code CASTEP, including technical issues concerning the calculations, and discuss the applicability and accuracy of the method. A number of examples are provided for comparing the results to both experiment and other density functional theory techniques.

Spectral features, chemical shifts, and absolute thresholds of electron energy loss near-edge structure (ELNES) and x-ray absorption near-edge structure (XANES) for selected compounds, i.e. TiO(2) (rutile), TiO(2) (anatase), SrTiO(3), Ti(2)O(3), Al(2)O(3), AlN and β-Ga(2)O(3), were calculated by a plane wave pseudopotential method. Experimental ELNES/XANES of those compounds were well reproduced when an excited pseudopotential, which includes a core hole, was used. In addition to the spectral features, it was found that chemical shifts among different compounds were also reproduced by correcting the contribution of the excited pseudopotentials to the energy of the core orbital.

A riddle solved! Despite its simple formula, the structure of the (SCN)(x) polymer has remained elusive since its first synthesis in 1929. From energetics as well as NMR chemical shifts, based on DFT calculations, we have strong evidence that it is indeed a tangle of linear chains, made up from N-linked S(2)C(2)N five-membered rings.Molecular fragments and crystal structures based on proposed structures for polythiocyanogen were studied by using molecular and solid-state electronic structure calculations at the density functional theory level. The energetics and chemical shifts from both types of calculations indicate that a planar N-linked chain consisting of 1,2,4-dithiazole five-membered rings with adjacent rings pointing in opposite directions is the most likely local structure of the (SCN)(x) polymer.

Heteronuclear solid-state magic-angle spinning (MAS) NMR experiments for probing (15)N-(17)O dipolar and J couplings are presented for [(2)H(NH(3)),1-(13)C,(15)N,(17)O(2)]glycine.(2)HCl and [(15)N(2),(17)O(2)]uracil. Two-dimensional (15)N-(17)O correlation spectra are obtained using the R(3)-HMQC experiment; for glycine.(2)HCl, the intensity of the resolved peaks for the CO and C-O(2)H (17)O resonances corresponds to the relative magnitude of the respective (15)N-(17)O dipolar couplings. (17)O-(15)N REDOR curves are presented for glycine.(2)HCl; fits of the initial buildup (DeltaS/S < 0.2) yield effective dipolar couplings in agreement with (+/-20%) the root-sum-squared dipolar couplings determined from the crystal structure. Experimental (15)N-(17)O REAPDOR curves for the (15)N resonances in glycine.(2)HCl and uracil fit well to the universal curve presented by Goldbourt et al. (J. Am. Chem. Soc. 2003, 125, 11194). Heteronuclear (13)C-(17)O and (15)N-(17)O J couplings were experimentally determined from fits of the quotient of the integrated intensity obtained in a heteronuclear and a homonuclear spin-echo experiment, S(Q)(tau) = S(HET)(tau)/S(HOM)(tau). For glycine.(2)HCl, (1)J(CO) was determined as 24.7 +/- 0.2 and 25.3 +/- 0.3 Hz for the CO and C-O(2)H resonances, respectively, while for uracil, the average of the two NH...O hydrogen-bond-mediated J couplings was determined as 5.1 +/- 0.6 Hz. In addition, two-bond intramolecular J couplings, (2)J(OO) = 8.8 +/- 0.9 Hz and (2)J(N1,N3) = 2.7 +/- 0.1 Hz, were determined for glycine.(2)HCl and uracil, respectively. Excellent agreement was found with J couplings calculated using the CASTEP code using geometrically optimized crystal structures for glycine.HCl [(1)J(CO)(CO) = 24.9 Hz, (1)J(CO)(COH) = 27.5 Hz, (2)J(OO) = 7.9 Hz] and uracil [(2h)J(N1,O4) = 6.1 Hz, (2h)J(N3,O4) = 4.6 Hz, (2)J(N1,N3) = 2.7 Hz].

A general approach for structural interpretation of local disorder in partially ordered solids is proposed, combining high-resolution two-dimensional (2D) nuclear magnetic resonance (NMR) and first principles calculations. We show that small chemical shift variations of the order of a ppm can be interpreted in detailed structural terms with advanced density functional theory methods. Focusing on a model system of bisphosphinoamine, we demonstrate that the existence and the spatial range of small amplitude disorder can be probed using quantitative statistical analyses of 2D NMR line shapes obtained from through-space correlation experiments collected using variable mixing times. We show how low-energy vibration modes calculated from first principles can be conveniently used not as a cause of disorder but, instead, to generate a basis set of physically plausible local distortions to describe candidate static distributions of local geometries. Calculations of (31)P NMR isotropic chemical shifts are then used for the first time to simulate 2D correlation lineshapes associated with these distortions, which permit their evaluation as a potential source of disorder by comparison to experimental 2D cross-peaks between phosphorus sites. This new type of structural constraints allows the identification of changes in the bonding geometry that most likely contribute to the local structural disorder. We thus identify at least one type of structural deformation that is compatible with the experimental 2D NMR data and is also within the order of magnitude of the "thermal ellipsoids" associated with the uncertainties on the atomic positions of the X-ray diffraction structure.

This article presents results of first-principles calculations of quadrupolar parameters measured by solid-state nuclear magnetic measurement (NMR) spectroscopy. Different computational methods based on density functional theory were used to calculate the quadrupolar parameters. Through a series of illustrations from different areas of solid state inorganic chemistry, it is shown how quadrupolar solid-state NMR properties can be tackled by a theoretical approach and can yield structural information., ((c) 2008 Wiley Periodicals, Inc.)

The high-resolution 27Al and 31P NMR spectra of two as-synthesized forms of the microporous aluminophosphate AlPO-14 and the corresponding calcined-dehydrated form were assigned using both "first-principles" calculations of NMR parameters (GIPAW, as implemented in NMR-CASTEP) and a 27Al-31P heteronuclear correlation NMR experiment (MQ-J-HETCOR) that exploits 27Al multiple-quantum coherences and J couplings to identify Al-O-P linkages. NMR parameters calculated from published AlPO-14 crystal structures, which are derived from powder X-ray diffraction (XRD) data, are in poor agreement with experiment and it was necessary to optimize the structure geometry using energy minimization before satisfactory agreement was obtained. Comparison of simulated powder XRD patterns from the experimental and the energy-minimized structures shows that the changes in relative atomic positions in the optimized structure are relatively small and yield only minor adjustments in the Bragg peak intensities. These results indicate that a combination of NMR spectroscopy and first-principles calculation of NMR parameters may soon be considered a generally useful step in the refinement of the structures of microporous materials derived from powder diffraction data.

Ammonia is an important compound with many uses, such as in the manufacture of fertilizers, explosives and pharmaceuticals. As an archetypal hydrogen-bonded system, the properties of ammonia under pressure are of fundamental interest, and compressed ammonia has a significant role in planetary physics. We predict new high-pressure crystalline phases of ammonia (NH(3)) through a computational search based on first-principles density-functional-theory calculations. Ammonia is known to form hydrogen-bonded solids, but we predict that at higher pressures it will form ammonium amide ionic solids consisting of alternate layers of NH(4)(+) and NH(2)(-) ions. These ionic phases are predicted to be stable over a wide range of pressures readily obtainable in laboratory experiments. The occurrence of ionic phases is rationalized in terms of the relative ease of forming ammonium and amide ions from ammonia molecules, and the volume reduction on doing so. We also predict that the ionic bonding cannot be sustained under extreme compression and that, at pressures beyond the reach of current static-loading experiments, ammonia will return to hydrogen-bonded structures consisting of neutral NH(3) molecules.

A recently developed method for calculating NMR J coupling in solid-state systems is applied to calculate hydrogen-bond-mediated (2h) J NN couplings across intra- or intermolecular N-H...N hydrogen bonds in two 6-aminofulvene-1-aldimine derivatives and the ribbon structure formed by a deoxyguanosine derivative. Excellent quantitative agreement is observed between the calculated solid-state J couplings and those previously determined experimentally in two recent spin-echo magic-angle-spinning NMR studies ( Brown, S. P. ; et al. Chem. Commun. 2002, 1852-1853 and Pham, T. N. ; et al. Phys. Chem. Chem. Phys. 2007, 9, 3416-3423 ). For the 6-aminofulvene-1-aldimines, the differences in (2h) J NN couplings in pyrrole and triazole derivatives are reproduced, while for the guanosine ribbons, an increase in (2h) J NN is correlated with a decrease in the N-H...N hydrogen-bond distance. J couplings are additionally calculated for isolated molecules of the 6-aminofulevene-1-aldimines extracted from the crystal with and without further geometry optimization. Importantly, it is shown that experimentally observed differences between J couplings determined by solution- and solid-state NMR are not solely due to differences in geometry; long-range electrostatic effects of the crystal lattice are shown to be significant also. J couplings that are yet to be experimentally measured are calculated. Notably, (2h) J NO couplings across N-H...O hydrogen bonds are found to be of a similar magnitude to (2h) J NN couplings, suggesting that their utilization and quantitative determination should be experimentally feasible.

Weak hydrogen bonding in uracil and 4-cyano-4'-ethynylbiphenyl, for which single-crystal diffraction structures reveal close CH...O=C and C[triple bond]CH...N[triple bond]C distances, is investigated in a study that combines the experimental determination of 1H, 13C, and 15N chemical shifts by magic-angle spinning (MAS) solid-state NMR with first-principles calculations using plane-wave basis sets. An optimized synthetic route, including the isolation and characterization of intermediates, to 4-cyano-4'-ethynylbiphenyl at natural abundance and with 13C[triple bond]13CH and 15N[triple bond]C labeling is described. The difference in chemical shifts calculated, on the one hand, for the full crystal structure and, on the other hand, for an isolated molecule depends on both intermolecular hydrogen bonding interactions and aromatic ring current effects. In this study, the two effects are separated computationally by, first, determining the difference in chemical shift between that calculated for a plane (uracil) or an isolated chain (4-cyano-4'-ethynylbiphenyl) and that calculated for an isolated molecule and by, second, calculating intraplane or intrachain nucleus-independent chemical shifts that quantify the ring current effects caused by neighboring molecules. For uracil, isolated molecule to plane changes in the 1H chemical shift of 2.0 and 2.2 ppm are determined for the CH protons involved in CH...O weak hydrogen bonding; this compares to changes of 5.1 and 5.4 ppm for the NH protons involved in conventional NH...O hydrogen bonding. A comparison of CH bond lengths for geometrically relaxed uracil molecules in the crystal structure and for geometrically relaxed isolated molecules reveals differences of no more than 0.002 A, which corresponds to changes in the calculated 1H chemical shifts of at most 0.1 ppm. For the C[triple bond]CH...N[triple bond]C weak hydrogen bonds in 4-cyano-4'-ethynylbiphenyl, the calculated molecule to chain changes are of similar magnitude but opposite sign for the donor 13C and acceptor 15N nuclei. In uracil and 4-cyano-4'-ethynylbiphenyl, the CH hydrogen-bonding donors are sp2 and sp hybridized, respectively; a comparison of the calculated changes in 1H chemical shift with those for the sp3 hybridized CH donors in maltose (Yates et al. J. Am. Chem. Soc. 2005, 127, 10216) reveals no marked dependence on hybridization for weak hydrogen-bonding strength.

We have combined a computational search strategy with first-principles density-functional-theory calculations to identify metastable phases of H(2)O under pressure. The most stable structures consist of water molecules, while the most energetic metastable phases consist of oxygen and hydrogen molecules. In between lie many other metastable phases, consisting of various combinations of a few small molecules. It may be possible to synthesize some of these metastable phases, and we use our results to understand the nature of the crystalline metastable phase of H(2)O recently synthesized by Mao et al. [Science 314, 636 (2006)].

17O, 15N, 13C, and 1H NMR chemical shieldings are calculated using density functional theory to differentiate among the three primarily helical forms, 310, alpha, and pi in polyalanine peptides under periodic boundary conditions. This study suggests 17O as the best observable, as it has been demonstrated to be sensitive to hydrogen bonding and highly affected by small changes in the polypeptide in helix conformations. This theoretical study seeks to characterize the subtle conformational differences of helical structures by NMR chemical shift observables which may lead to important questions in experimental structure determination on the basis of using chemical shifts to identify protein secondary structures.

Computations for chemical shifts of molecular organic compounds using the gauge-including projector augmented wave method and the NMR-CASTEP code are reviewed. The methods are briefly introduced, and some general aspects involving the sources of uncertainty in the results are explored. The limitations are outlined. Successful applications of the computations to problems of interpretation of NMR results are discussed and the range of areas in which useful information is obtained is illustrated by examples. The particular value of the computations for comparing shifts between resonances where the same chemical site is involved is emphasised. Such cases arise for shifts between different crystallographically independent molecules of the same chemical species, between polymorphs and for shift anisotropies and asymmetries., (Copyright © 2007 John Wiley & Sons, Ltd.)

A method to calculate NMR J-coupling constants from first principles in extended systems is presented. It is based on density functional theory and is formulated within a planewave-pseudopotential framework. The all-electron properties are recovered using the projector augmented wave approach. The method is validated by comparison with existing quantum chemical calculations of solution-state systems and with experimental data. The approach has also been applied to the silicophosphate, Si(5)O(PO(4))(6), giving (31)P-(29)Si-couplings which are in excellent agreement with experiment.

The 29Si and 17O NMR parameters of six polymorphs of MgSiO3 were determined through a combination of high-resolution solid-state NMR and first-principles gauge including projector augmented wave (GIPAW) formalism calculations using periodic boundary conditions. MgSiO3 is an important component of the Earth's mantle that undergoes structural changes as a function of pressure and temperature. For the lower pressure polymorphs (ortho-, clino-, and protoenstatite), all oxygen species in the 17O high-resolution triple-quantum magic angle spinning (MAS) NMR spectra were resolved and assigned. These assignments differ from those tentatively suggested in previous work on the basis of empirical experimental correlations. The higher pressure polymorphs of MgSiO3 (majorite, akimotoite, and perovskite) are stabilized at pressures corresponding to the Earth's transition zone and lower mantle, with perovskite being the major constituent at depths >660 km. We present the first 17O NMR data for these materials and confirm previous 29Si work in the literature. The use of high-resolution multiple-quantum MAS (MQMAS) and satellite-transition MAS (STMAS) experiments allows us to resolve distinct oxygen species, and full assignments are suggested. The six polymorphs exhibit a wide variety of structure types, providing an ideal opportunity to consider the variation of NMR parameters (both shielding and quadrupolar) with local structure, including changes in coordination number, local geometry (bond distances and angles), and bonding. For example, we find that, although there is a general correlation of increasing 17O chemical shift with increasing Si-O bond length, the shift observed also depends upon the exact coordination environment.

The 13C NMR chemical shifts for functionalized (7,0), (8,0), (9,0), and (10,0) single-walled carbon nanotubes (SWNTs) have been studied computationally using gauge-including projector-augmented plane-wave (GIPAW) density functional theory (DFT). The functional groups NH, NCH3, NCH2OH, and CH2NHCH2 have been considered, and different sites where covalent addition or substitution may occur have been examined. The shifts of the carbons directly attached to the group are sensitive to the bond which has been functionalized and may, therefore, be used to identify whether the group has reacted with a parallel or a diagonal C-C bond. The addition of NH to a parallel bond renders the functionalized carbons formally sp3-hybridized, yielding shifts of around 44 ppm, independent of the SWNT radius. Reaction with a diagonal bond retains the formal sp2 hybridization of the substituted carbons, and their shifts are slightly lower or higher than those of the unsubstituted carbon atoms. The calculated 1H NMR shifts of protons in the functional groups are also dependent upon the SWNT-group interaction. Upon decreasing the degree of functionalization for the systems where the group is added to a parallel bond, the average chemical shift of the unfunctionalized carbons approaches that of the pristine tube. At the same time, the shifts of the functionalized carbons remain independent upon the degree of functionalization. For the SWNTs where N-R attaches to a parallel bond, the average shift of the sp2 carbons was found to be insensitive to the substituent R. Moreover, the shifts of the functionalized sp3 carbons, as well as of the carbons within the group itself, are independent of the SWNT radius. The results indicate that a wealth of knowledge may be obtained from the 13C NMR of functionalized SWNTs.

The nuclear magnetic resonance (NMR) shielding and electric field gradient (EFG) tensors of three polymorphs of Mg(2)SiO(4), forsterite (alpha-Mg(2)SiO(4)), wadsleyite (beta-Mg(2)SiO(4)) and ringwoodite (gamma-Mg(2)SiO(4)), have been calculated using a density functional theory (DFT) approach with a planewave basis set and pseudopotential approximation. These Mg(2)SiO(4) polymorphs are the principal components of the Earth down to depths of 660 km and have been proposed as the hosts of water in the Earth's upper mantle and transition zone. A comparison of our calculations with single-crystal spectroscopic data in the literature for the alpha-polymorph, forsterite, shows that both the magnitude and orientation of the shielding and EFG tensors for O and Si can be obtained with sufficient accuracy to distinguish subtle differences in atomic positions between published structures. We compare calculated (17)O MAS NMR quadrupolar powder lineshapes directly with experimental lineshapes and show that we are able to reproduce them within the precision with which the NMR parameters may be determined from multi-parameter fitting. The relatively small amounts of sample available for the beta- and gamma-polymorphs, arising from the high pressures required for synthesis, has hindered the extraction of NMR parameters in previous work. The application of DFT calculations to these high-pressure polymorphs confirms previous spectral assignments, and provides deeper insight into the empirical correlations and observations reported in the literature. These first-principles methods are highly promising for the determination of local bonding in more complex materials, such as the hydrated forms of Mg(2)SiO(4), by aiding analysis of their multinuclear NMR spectra.

The (13)C CPMAS spectrum is presented for the polymorph of oxybuprocaine hydrochloride which is stable at room temperature, i.e. Mod. II degrees . It shows crystallographic splittings arising from the fact that there are two molecules, with substantially different conformations, in the asymmetric unit. An INADEQUATE two-dimensional experiment was used to link signals for the same independent molecule. The chemical shifts are discussed in relation to the crystal structure. Of the four ethyl groups attached to NH(+) nitrogens, one gives rise to unusually low chemical shifts, very different from those of the other three ethyl groups. This is attributed empirically to gamma-gauche conformational effects, as is confirmed by shielding computations. These considerations allow (13)C signals to be assigned to specific carbons in the two crystallographically inequivalent molecules in the crystal structure. Indeed, information about the conformations is inherent in the NMR spectrum, which thus provides data of crystallographic significance. A (13)C/(1)H HETCOR experiment enabled resolution to be obtained in the (1)H dimension and allowed (1)H and (13)C signals for the same independent molecule to be linked.

NMR chemical shifts were calculated for semiconducting (n,0) single-walled carbon nanotubes (SWNTs) with n ranging from 7 to 17. Infinite isolated SWNTs were calculated using a gauge-including projector-augmented plane-wave (GIPAW) approach with periodic boundary conditions and density functional theory (DFT). In order to minimize intertube interactions in the GIPAW computations, an intertube distance of 8 A was chosen. For the infinite tubes, we found a chemical shift range of over 20 ppm for the systems considered here. The SWNT family with lambda = mod(n, 3) = 0 has much smaller chemical shifts compared to the other two families with lambda = 1 and lambda = 2. For all three families, the chemical shifts decrease roughly inversely proportional to the tube's diameter. The results were compared to calculations of finite capped SWNT fragments using a gauge-including atomic orbital (GIAO) basis. Direct comparison of the two types of calculations could be made if benzene was used as the internal (computational) reference. The NMR chemical shifts of finite SWNTs were found to converge very slowly, if at all, to the infinite limit, indicating that capping has a strong effect (at least for the (9,0) tubes) on the calculated properties. Our results suggest that (13)C NMR has the potential for becoming a useful tool in characterizing SWNT samples.

We show how powder samples at natural isotopic abundance can be assigned to crystal structures by using high-resolution proton and carbon-13 solid-state NMR spectra in combination with first principles calculations. Homonuclear proton double-quantum spectra in combination with through-bond proton-carbon HSQC spectra are used to assign the NMR spectra. We then show that the proton chemical shifts can be included in the process of assigning the spectra to a crystal structure using first principles calculations. The method is demonstrated on the K salt of penicillin G.

The distorted perovskites NaTaO(3) and NaNbO(3) have been studied using (23)Na multiple-quantum (MQ) MAS NMR. NaTaO(3) was prepared by high temperature solid state synthesis and the NMR spectra are consistent with the expected room temperature structure of the material (space group Pbnm), with a single crystallographic sodium site. Two samples of NaNbO(3) were studied. The first, a commercially available sample which was annealed at 900 degrees C, showed two crystallographic sodium sites, as expected for the room temperature structure of the material (space group Pbcm). The second sample, prepared by a low temperature hydrothermal method, showed the presence of four sodium sites, two of which match the expected room temperature structure and the second pair, another polymorph of the material (space group P21ma). This is consistent with powder X-ray diffraction data which showed weak extra peaks which can be accounted for by the presence of this second polymorph. Density functional theory (DFT) calculations support our conclusions, and aid assignment of the NMR spectra. Finally, we discuss the measured NMR parameters in relation to other studies of sodium in high coordination sites in the solid state.

High-pressure phases of silane SiH4 are predicted using first-principles electronic structure methods. We search for low-enthalpy structures by relaxing from randomly chosen initial configurations, a strategy which is demonstrated to work well for unit cells containing up to at least ten atoms. We predict that silane will metallize at higher pressures than previously anticipated but might show high-temperature superconductivity at experimentally accessible pressures.

Studies on the solid-state structure of two polymorphs of 4-methyl-2-nitroacetanilide (MNA) were conducted using magic-angle spinning (13)C, (15)N and (1)H NMR spectroscopy, together with first-principles computations of NMR shielding (including use of a program that takes explicit account of the translational symmetry inherent in crystalline structures). The effects on (13)C chemical shifts of side-chain rotations have been explored. Information derived from these studies was then incorporated within a systematic space-search methodology for elucidation of trial crystallographic structures from powder XRD.

A (13)C CPMAS NMR experiment at high field (11.7 T) has produced significantly improved dispersion for the alpha form of testosterone, allowing revisions and extensions to be made to the assignments. Correlations shown by an INADEQUATE two-dimensional spectrum, recorded at 16.5 T, have allowed the components of most of the doublet signals to be grouped into two sets (for the two crystallographically independent molecules). First-principles computations, employing a fully solid-state approach, have been used to obtain values for the crystallographic splittings, which are discussed in relation to the experimental values. This procedure enables assignments to the two groups to be suggested for all but one of the remaining doublet signals. It also allows the two sets of signals to be identified specifically to the two independent molecules in the crystal structure. Computations were also carried out for the beta form of testosterone (a dihydrate). The shift differences between the alpha and beta forms were compared with the experimental data, with encouraging results. Comparisons were also made between computed and experimental shielding anisotropies and asymmetries for three of the carbons of the alpha form. The methodology has a high potential for future applications, though more examples need to be evaluated before general conclusions can be drawn.

The use of the spin-dependent pseudopotentials has been shown to markedly enhance the transferability of the commonly used spin-neutral pseudopotential method for the study of the structural and magnetic properties of transition-metal-containing materials. Unfortunately, because the method was based on the rather expensive norm-conserving pseudopotential formalism, the method was limited to the study of fairly small systems. Here we present an extension of the spin-dependent pseudopotential method for the far more computationally advantageous ultrasoft formalism and show that it is very easy to add such a feature to any preexisting computer code. We benchmark our new method by comparing to previously published results and then apply it to the study of several relevant test cases: bulk Ni, Fe, and Co, as well as a Pd atomic wire.

13C, 14N, 15N, 17O, and 35Cl NMR parameters, including chemical shift tensors and quadrupolar tensors for 14N, 17O, and 35Cl, are calculated for the crystalline forms of various amino acids under periodic boundary conditions and complemented by experiment where necessary. The 13C shift tensors and 14N electric field gradient (EFG) tensors are in excellent agreement with experiment. Similarly, static 17O NMR spectra could be precisely simulated using the calculation of the full chemical shift (CS) tensors and their relative orientation with the EFG tensors. This study allows correlations to be found between hydrogen bonding in the crystal structures and the 17O NMR shielding parameters and the 35Cl quadrupolar parameters, respectively. Calculations using the two experimental structures for L-alanine have shown that, while the calculated isotropic chemical shift values of 13C and 15N are relatively insensitive to small differences in the experimental structure, the 17O shift is markedly affected.

Two-dimensional (1)H-(13)C MAS-J-HMQC solid-state NMR spectra of the two anomeric forms of maltose at natural abundance are presented. The experimental (1)H chemical shifts of the CH and CH(2) protons are assigned using first-principles chemical shift calculations that employ a plane-wave pseudopotential approach. Further calculations show that the calculated change in the (1)H chemical shift when comparing the full crystal and an isolated molecule is a quantitative measure of intermolecular C-H...O weak hydrogen bonding. Notably, a clear correlation between a large chemical shift change (up to 2 ppm) and both a short H...O distance (<2.7 A) and a CHO bond angle greater than 130 degrees is observed, thus showing that directionality is important in C-H...O hydrogen bonding.

The spin-other-orbit (SOO) contribution to the g-tensor (DeltagSOO) of electron paramagnetic resonance arises due to the interaction of electron-spin magnetic moment with the magnetic field produced by the orbital motion of other electrons. A similar mechanism is responsible for the leading term in nuclear magnetic-shielding tensors sigma. We demonstrate that analogous to sigma, paramagnetic DeltagSOO contribution exhibits a pronounced dependence on the choice of the magnetic-field gauge. The gauge corrections to DeltagSOO are similar in magnitude, and opposite in sign, to the paramagnetic SOO term. We calculate gauge-invariant DeltagSOO values using gauge-including atomic orbitals and density-functional theory. For organic radicals, complete gauge-invariant DeltagSOO values typically amount to less than 500 parts per million (ppm), and are small compared to other g-tensor contributions. For the first-row transition-metal compounds, DeltagSOO may contribute several thousand ppm to the g-tensor, but are negligible compared to the remaining deviations from experiment. With popular choices for the magnetic-field gauge, the individual gauge-variant contributions may be an order of magnitude higher, and do not provide a reliable estimation of DeltagSOO.

The temperature dependence of (17)O and (25)Mg NMR chemical shifts in solid MgO have been calculated using a first-principles approach. Density functional theory, pseudopotentials, a plane-wave basis set, and periodic boundary conditions were used both to describe the motion of the nuclei and to compute the NMR chemical shifts. The chemical shifts were obtained using the gauge-including projector augmented wave method. In a crystalline solid, the temperature dependence is due to both (i) the variation of the averaged equilibrium structure and (ii) the fluctuation of the atoms around this structure. In MgO, the equilibrium structure at each temperature is uniquely defined by the cubic lattice parameters, which we take from experiment. We evaluate the effect of the fluctuations within a quasiharmonic approximation. In particular, the dynamical matrix, defining the harmonic Hamiltonian, has been computed for each equilibrium volume. This harmonic Hamiltonian was used to generate nuclear configurations that obey quantum statistical mechanics. The chemical shifts were averaged over these nuclear configurations. The results reproduce the previously published experimental NMR data measured on MgO between room temperature and 1000 degrees C. It is shown that the chemical shift behavior with temperature cannot be explained by thermal expansion alone. Vibrational corrections due to the fluctuations of atoms around their equilibrium position are crucial to reproduce the experimental results.

We have computed the (17)O NMR parameters of an amorphous calcium aluminosilicate (CAS) from first-principles. The atomic coordinates of a CAS glass of composition (CaO)(0.21)(Al(2)O (3))(0.12)(SiO (2))(0.67) were obtained by quenching a liquid to room temperature by the means of ab initio molecular-dynamics simulations of the Car-Parrinello type. The structure of the glass is found to be overall in good agreement with diffraction experiments. Some excess nonbridging O (NBO) atoms are found and are compensated by tricluster O atoms, i.e., by 3-fold coordinated O atoms to 4-fold coordinated Al or Si atoms. The glass coordinates were used to compute the (17)O NMR parameters using GGA-DFT and a correction of the Ca 3d orbital energy. The chemical shifts and the electric field gradients were obtained with the gauge including projector augmented-wave (GIPAW) and the projector augmented-wave (PAW) methods, respectively. The simulated 2D-3QMAS NMR spectrum of the CAS glass is in very good agreement with the available experimental data, notably because it takes into account the disorder present in the glass. This agreement further validates our CAS glass model. We show that the oxygen triclusters are not visible in a 2D-3QMAS NMR (17)O spectrum since their NMR parameters overlap with those of the Al-O-Si, Si-O-Si, or Al-O-Al sites. Finally, correlations between the structural characteristics and the values of the NMR parameters are extracted from the calculation with the aim of helping the interpretation of NMR spectra of glasses of similar compositions.

The 1H, 13C and 19F magic-angle spinning NMR spectra have been recorded for Form 1 of flurbiprofen. In the case of 19F, spinning sideband analysis has produced data for the components of the shielding tensor. The chemical shift of the hydrogen-bonded proton was found to be 14.0 ppm. Shielding parameters for all three nuclei have been calculated using Density Functional Theory (DFT) together with the Gauge Including Projector Augmented Wave (GIPAW) method which takes full allowance for the repetition inherent in crystalline structures. Such computations were made for the reported geometry, for a structure with all the atomic positions relaxed using DFT, and with only the hydrogen positions relaxed. The relationships of the computed shifts to those observed are discussed. In general, the correlations are good.

We apply density functional theory (DFT) to the calculation of the (17)O NMR parameters in Ca and Mg oxides and aluminosilicates. We study the accuracy of the Perdew, Burke, and Ernzerhof (PBE) generalized-gradient approximation to DFT in the description of these systems and the origin of the experimentally observed large dependence of the (17)O chemical shift on the alkaline earth ion. We find that (i) the partially covalent nature of the Ca-O bond has a huge impact on the O chemical shifts. The Ca-O covalence alone explains why in Ca oxides and aluminosilicates the (17)O chemical shifts are much more deshielded than those of the corresponding Mg compounds. (ii) The Ca-O covalence is overestimated by the PBE functional. Thus PBE-DFT is not able to reproduce the measured (17)O NMR parameters in Ca oxide and Ca aluminosilicates. (iii) It is possible to correct for the PBE-DFT deficiency in a simple and transferable way and to predict very accurate (17)O NMR parameters. Such accuracy allows us to assign the (17)O NMR spectra of two important model systems: the grossite aluminate (CaAl(4)O(7)) and the wollastonite (CaSiO(3)) silicate.

1H, 13C, 17O and 31P NMR parameters, including chemical shift tensors and quadrupolar parameters for 17O, were calculated for phenylphosphonic acid, C6H5PO(OH)2, under periodic boundary conditions. The results are in very good agreement with experimental data and permit the unambiguous assignment of all the sites present in the structure. In particular, the 17O NMR parameters of the P=O and P-OH environments were precisely determined, which should help in the characterization of the bonding mode of phosphonate molecules in hybrid solids. Moreover, the effect of intermolecular interactions on the NMR parameters were investigated by comparing the results of the calculations in the crystal and in an isolated molecule of phenylphosphonic acid., (Copyright 2004 John Wiley & Sons, Ltd.)

We show how to describe the coupling of electrons to nonuniform magnetic fields in the framework of the widely used norm-conserving pseudopotential approximation for electronic structure calculations. Our derivation applies to magnetic fields that are smooth on the scale of the core region. The method is validated by application to the calculation of the magnetic susceptibility of molecules within density functional theory (DFT) in the local density approximation. Our results are compared with high-quality all-electron DFT results obtained using Gaussian basis sets and another recently proposed pseudopotential formalism.