Solid-state NMR spectra and long intradimer bonds in the pi-[TCNE]22- dianion

The principal (13)C chemical-shift values for the pi-[TCNE](2)(2-) dimer anion within an array of counterions have been measured to understand better the electronic structure of these atypical chemical species in several related TCNE-based structures. The structure of pi-[TCNE](2)(2-) is unusual as it contains two very long C-C bond lengths (ca. 2.9 Angstroms) between the two monomeric units and has been found to exist as a singlet state, suggestive of a (1)A(1g) (b(2u)(2)b(1g)(0)) electronic configuration. A systematic study of several oxidation states of [TCNE](n) (n = 0, 1-, 2-) was conducted to determine how the NMR chemical-shift tensor values change as a function of electronic structure and to understand the interactions that lead to spin-pairing of the monomer units. The density functional theory (DFT) calculated nuclear shielding tensors are correlated with the experimentally determined principal chemical-shift values. Such theoretical methods provide information on the tensor magnitudes and orientations of their principal tensor components with respect to the molecular frame. Both theoretical and experimental ethylenic chemical-shielding tensors reveal high sensitivity in the component, delta(perpendicular), lying in the monomer molecular plane and perpendicular to the pi-electron plane. This largest shift dependence on charge density is observed to be about -111 ppm/e(-) for delta(perpendicular). The component in the molecular plane but parallel to the central C=C bond, delta(parallel), exhibits a sensitivity of approximately -43 ppm/e(-). However, the out-of-plane component delta'(perpendicular) shows a minimal dependence of -2.6 ppm/e(-) on the oxidation state (n) of [TCNE](n). These relative values support the claim that it is changes within the ethylenic pi-electrons and not the sigma-electrons that best account for the dramatic variations in bonding and shift tensors in this series of compounds. Concerning the intraion bonding, relatively weak Wiberg bond orders between the two monomeric components of the dimer correlate with the long bonds linking the two [TCNE(*)](-) monomers. The chemical-shift tensors for the cyano group, compared to the ethylene shifts, exhibit a reduced sensitivity on the TCNE oxidation state. The experimental principal chemical-shift components agree (within typical errors) with the calculated quantum mechanical shieldings used to correlate the bonding. The embedded ion model (EIM) was used to investigate the typically large electrostatic lattice potential in these ionic materials. Chemical-shielding principal values calculated with the EIM model differ from experiment by +/-3.82 ppm on average, whereas in the absence of an electrostatic field model, the experimental and theoretical results agree by +/-4.42 ppm, which is only a modest increase in error considering the overall ionic magnitudes associated with the tensor variations. Apparently, the effects of the sizable long-range electrostatic fields cancel when the shifts are computed because of lattice symmetry.     

Relationships between 31P chemical shift tensors and conformation of nucleic acid backbone: a DFT study

Density functional theory (DFT) has been applied to study the conformational dependence of 31P chemical shift tensors in B-DNA. The gg and gt conformations of backbone phosphate groups representing BI- and BII-DNA have been examined. Calculations have been carried out on static models of dimethyl phosphate (dmp) and dinucleoside-3',5'-monophosphate with bases replaced by hydrogen atoms in vacuo as well as in an explicit solvent. Trends in 31P chemical shift anisotropy (CSA) tensors with respect to the backbone torsion angles alpha, zeta, beta, and epsilon are presented. Although these trends do not change qualitatively upon solvation, quantitative changes result in the reduction of the chemical shift anisotropy. For alpha and zeta in the range from 270 degrees to 330 degrees and from 240 degrees to 300 degrees , respectively, the delta22 and delta33 principal components vary within as much as 30 ppm, showing a marked dependence on backbone conformation. The calculated 31P chemical shift tensor principal axes deviate from the axes of O-P-O bond angles by at most 5 degrees . For solvent models, our results are in a good agreement with experimental estimates of relative gg and gt isotropic chemical shifts. Solvation also brings the theoretical deltaiso of the gg conformation closer to the experimental gg data of barium diethyl phosphate.     

Experimental and computational characterization of the 17O quadrupole coupling and magnetic shielding tensors for p-nitrobenzaldehyde and formaldehyde

We have used solid-state 17O NMR experiments to determine the 17O quadrupole coupling (QC) tensor and chemical shift (CS) tensor for the carbonyl oxygen in p-nitro-[1-(17)O]benzaldehyde. Analyses of solid-state 17O NMR spectra obtained at 11.75 and 21.15 T under both magic-angle spinning (MAS) and stationary conditions yield the magnitude and relative orientation of these two tensors: CQ = 10.7 +/- 0.2 MHz, etaQ = 0.45 +/- 0.10, delta11 = 1050 +/- 10, delta22 = 620 +/- 10, delta33 = -35 +/- 10, alpha = 90 +/- 10, beta = 90 +/- 2, gamma = 90 +/- 10 degrees. The principal component of the 17O CS tensor with the most shielding, delta33, is perpendicular to the H-C=O plane, and the tensor component with the least shielding, delta11, lies along the C=O bond. For the 17O QC tensor, the largest (chi(zz)) and smallest (chi(xx)) components are both in the H-C=O plane being perpendicular and parallel to the C=O bond, respectively. This study represents the first time that these two fundamental 17O NMR tensors have been simultaneously determined for the carbonyl oxygen of an aldehyde functional group by solid-state 17O NMR. The reported experimental solid-state 17O NMR results provide the first set of reliable data to allow evaluation of the effect of electron correlation on individual CS tensor components. We found that the electron correlation effect exhibits significant influence on 17O chemical shielding in directions within the H-C=O plane. We have also carefully re-examined the existing experimental data on the 17O spin-rotation tensor for formaldehyde and proposed a new set of best "experimental" 17O chemical shielding tensor components: sigma11 = -1139 +/- 80, sigma22 = -533 +/- 80, sigma33 = 431 +/- 5, and sigma(iso) = -414 +/- 60 ppm. Using this new set of data, we have evaluated the accuracy of quantum chemical calculations of the 17O CS tensors for formaldehyde at the Hartree-Fock (HF), density-functional theory (DFT), M?ller-Plesset second-order perturbation (MP2), and coupled-cluster singles and doubles (CCSD) levels of theory. The conclusion is that, while results from HF and DFT tend to underestimate the electron correlation effect, the MP2 method overestimates its contribution. The CCSD results are in good agreement with the experimental data.     

A solid-state nitrogen-15 NMR and ab initio study of nitrobenzenes

Insight into the unexpectedly small range of isotropic nitrogen chemical shifts in nitrobenzene derivatives is gained through measurements of the chemical shift (CS) tensor by solid-state NMR experiments and ab initio molecular orbital (MO) and density functional theory (DFT) calculations. The principal components, delta(ii), of the (15)N CS tensors have been measured for nitrobenzene, 4-nitroaniline, 4-nitrotoluene, 4-nitroanisole, 4-nitroacetophenone, nitromesitylene, and 2,4,6-tri-tert-butylnitrobenzene. No obvious correlations of the delta(ii) values with traditional reactivity parameters were observed. The CS tensor components change significantly for the para-substituted nitrobenzenes, but these variations nearly cancel to yield isotropic shifts that fall in a range of only 3 ppm. Ab initio calculations of the delta(ii) values at the HF level are in poor agreement with the experimental values, whereas MP2 calculations and DFT calculations employing the B3LYP functional are in better agreement with experiment. The calculated (B3LYP/6-311G) delta(ii) values follow a trend in which delta(11) and delta(33) increase while delta(22) decreases with the accepted electron withdrawing ability of the para substituent. These changes tend to cancel yielding a variation in delta(iso) of only 4 ppm. These calculations indicate that the CS tensor has the same orientation as the carbon CS tensor in the isoelectronic benzoate anion: delta(11) bisects the O-N-O angle, delta(33) is perpendicular to the NO(2) plane, and delta(22) is in the NO(2) plane and perpendicular to delta(11).     

Probing local structure in zeolite frameworks: ultrahigh-field NMR measurements and accurate first-principles calculations of zeolite 29Si magnetic shielding tensors

The principal components of zeolite 29Si magnetic shielding tensors have been accurately measured and calculated for the first time. The experiments were performed at an ultrahigh magnetic field of 21.1 T in order to observe the small anisotropies of the 29Si shielding interactions that arise for Si atoms in near-tetrahedral geometries. A robust two-dimensional (2D) chemical shift anisotropy (CSA) recoupling pulse sequence was employed that enables quasi-static powder patterns to be resolved according to the isotropic chemical shifts. For the zeolites Sigma-2 and ZSM-12, it is demonstrated that the 29Si chemical shift (CS) tensor components measured by the recoupling experiment are in excellent agreement with those determined from spinning sidebands in slow magic-angle spinning (MAS) experiments. For the zeolite ZSM-5, the principal components of the 29Si CS tensors of 15 of the 24 Si sites were measured using the 2D CSA recoupling experiment, a feat that would not be possible with a slow MAS experiment due to the complexity of the spectrum. A simple empirical relationship between the 29Si CS tensors and local structural parameters could not be established. However, the 29Si magnetic shielding tensors calculated using Hartree-Fock ab initio calculations on clusters derived from the crystal structures are in excellent agreement with the experimental results. The accuracy of the calculations is strongly dependent on the quality of the crystal structure used in the calculation, indicating that the 29Si magnetic shielding interaction is extremely sensitive to the local structure around each Si atom. It is anticipated that the measurement and calculation of 29Si shielding tensors could be incorporated into the "NMR crystallography" of zeolites and other related silicate materials, possibly being used for structure refinements that may lead to crystal structures with very accurate Si and O atomic coordinates.     

A comparison of quantum chemical models for calculating NMR shielding parameters in peptides: mixed basis set and ONIOM methods combined with a complete basis set extrapolation

This article compares several quantum mechanical approaches to the computation of chemical shielding tensors in peptide fragments. First, we describe the effects of basis set quality up to the complete basis set (CBS) limit and level of theory (HF, MP2, and DFT) for four different atoms in trans N-methylacetamide. For both isotropic shielding and shielding anisotropy, the MP2 results in the CBS limit show the best agreement with experiment. The HF values show quite a different tendency to MP2, and even in the CBS limit they are far from experiment for not only the isotropic shielding of carbonyl carbon but also most shielding anisotropies. In most cases, the DFT values differ systematically from MP2, and small basis-set (double- or triple-zeta) results are often fortuitously in better agreement with the experiment than the CBS ones. Second, we compare the mixed basis set and ONIOM methods, combined with CBS extrapolation, for chemical shielding calculations at a DFT level using various model peptides. From the results, it is shown that the mixed basis set method provides better results than ONIOM, compared to CBS calculations using the nonpartitioned full systems. The information studied here will be useful in guiding the selection of proper quantum chemical models, which are in a tradeoff between accuracy and cost, for shielding studies of peptides and proteins.     

Acetylenic carbon-13 chemical shift tensors for diphenylacetylene and (eta2-Diphenylacetylene)Pt(PPh3)2: a solid-state NMR and theoretical study

The structure of (eta2-diphenylacetylene)Pt(PPh3) (2), as well as those of its dichloromethane and benzene solvates, is determined via X-ray crystallography. An investigation of the chemical shift (CS) tensors of the 13C-labeled carbons in Ph13C13CPh and (eta2-Ph13C13CPh)Pt(PPh3)2.(C6H6) is carried out via analysis of 13C NMR spectra from stationary solid samples. The principal components of the CS tensors as well as their orientations with respect to the 13C,13C internuclear vector are determined. DFT calculations of these CS tensors are in close agreement with the experimental values. For diphenylacetylene (tolane), the orientations and principal-component magnitudes of the alkynyl carbon CS tensors are comparable to those for other alkynyl carbons, although the CS tensor is not axially symmetric in this case. Coordination to platinum causes a change in the CS tensor orientation and a net increase in the isotropic chemical shift, resulting from a significant increase in two principal components (delta11 and delta33) while the third (delta22) decreases only slightly. The measured carbon CS tensors in the platinum complex bear a striking similarity to those of the alkenyl carbons in trans-Ph(H)C=C(H)Ph, and a short theoretical discussion of these observations is presented.     

A combined experimental and quantum chemistry study of selenium chemical shift tensors

A comprehensive investigation of selenium chemical shift tensors is presented. Experimentally determined chemical shift tensors were obtained from solid-state 77Se NMR spectra for several organic, organometallic, or inorganic selenium-containing compounds. The first reported indirect spin-spin coupling between selenium and chlorine is observed for Ph(2)SeCl(2) where 1J(77Se,35Cl)iso is 110 Hz. Selenium magnetic shielding tensors were calculated for all of the molecules investigated using zeroth-order regular approximation density functional theory, ZORA DFT. The computations provide the orientations of the chemical shift tensors, as well as a test of the theory for calculating the magnetic shielding interaction for heavier elements. The ZORA DFT calculations were performed with nonrelativistic, scalar relativistic, and scalar with spin-orbit relativistic levels of theory. Relativistic contributions to the magnetic shielding tensor were found to be significant for (NH4)2WSe4 and of less importance for organoselenium, organophosphine selenide, and inorganic selenium compounds containing lighter elements.     

Ring current effects in crystals. Evidence from 13C chemical shift tensors for intermolecular shielding in 4,7-di-t-butylacenaphthene versus 4,7-di-t-butylacenaphthylene

13C chemical shift tensor data from 2D FIREMAT spectra are reported for 4,7-di-t-butylacenaphthene and 4,7-di-t-butylacenaphthylene. In addition, calculations of the chemical shielding tensors were completed at the B3LYP/6-311G** level of theory. While the experimental tensor data on 4,7-di-t-butylacenaphthylene are in agreement with theory and with previous data on polycyclic aromatic hydrocarbons, the experimental and theoretical data on 4,7-di-t-butylacenaphthene lack agreement. Instead, larger than usual differences are observed between the experimental chemical shift components and the chemical shielding tensor components calculated on a single molecule of 4,7-di-t-butylacenaphthene, with a root mean square (rms) error of +/-7.0 ppm. The greatest deviation is concentrated in the component perpendicular to the aromatic plane, with the largest value being a 23 ppm difference between experiment and theory for the 13CH2 carbon delta11 component. These differences are attributed to an intermolecular chemical shift that arises from the graphitelike, stacked arrangement of molecules found in the crystal structure of 4,7-di-t-butylacenaphthene. This conclusion is supported by a calculation on a trimer of molecules, which improves the agreement between experiment and theory for this component by 14 ppm and reduces the overall rms error between experiment and theory to 4.0 ppm. This intermolecular effect may be modeled with the use of nuclei independent chemical shieldings (NICS) calculations and is also observed in the isotropic 1H chemical shift of the CH2 protons as a 4.2 ppm difference between the solution value and the solid-state chemical shift measured via a 13C-1H heteronuclear correlation experiment.     

TensorView: A software tool for displaying NMR tensors

The representation of nuclear magnetic resonance (NMR) tensors as surfaces on three-dimensional molecular models is an information-rich presentation that highlights the geometric relationship between tensor principal components and the underlying molecular and electronic structure. Here, we describe a new computational tool, TensorView, for depicting NMR tensors on the molecular framework. This package makes use of the graphical interface and built-in molecular display functionality present within the Mathematica programming environment and is robust for displaying tensor properties from a broad range of commercial and user-specific computational chemistry packages. Two mathematical forms for representing tensor interaction surfaces are presented, the popular ellipsoidal construct and the more technically correct “ovaloid” form. Examples are provided for chemical shielding and shift tensors, dipole–dipole and quadrupolar couplings, and atomic anisotropic displacement parameters (thermal ellipsoids) derived from NMR crystallography.     

First-principles calculations within periodic boundary conditions of the NMR shielding tensor for a transition metal nucleus in a solid state system: The example of 51V in AlVO4

We present the first density functional theory based calculations of NMR shielding parameters for a transition metal nucleus using periodic boundary conditions. These calculations employ the gauge-including projected augmented-wave pseudopotential approach. The quality of this method is discussed by comparing experimental and calculated chemical shift tensor eigenvalues for the quadrupolar 51V nucleus in the diamagnetic solid-state compound AlVO4. Furthermore, the combination of shielding tensor with fast and accurate projector augmented-wave electric field gradient tensor calculations allows us to determine the relative orientation of these two tensors.     

31P chemical shift tensors for canonical and non-canonical conformations of nucleic acids: a DFT study and NMR implications

31P chemical shift anisotropy (CSA) tensors have been calculated for a set of selected DNA and RNA backbone conformations using density functional theory. The set includes canonical A-RNA, A-DNA, BI-DNA, BII-DNA, ZI-DNA, and ZII-DNA as well as four A-RNA-type, seven non-A-RNA-type, and three non-canonical DNA conformations. Hexahydrated dimethyl phosphate has been employed as a model. The 31P chemical shift tensors obtained are discussed in terms of similarities in the behavior observed for gauche-gauche (gg) and gauche-trans (gt) conformations around the P-O bonds. We show that torsion angles alpha and zeta are major determinants of the isotropic chemical shift deltaiso and of the deltaCSA11 component of the traceless chemical shift tensor, which is revealed in separate ranges of both deltaiso and deltaCSA11 for gg- and gt-conformers, respectively. A clear distinction between the two conformation types has not been found for the deltaCSA22 and deltaCSA33 components, which is attributed to their different directional properties. The 31P CSA tensors exhibit considerable variations resulting in large spans of approximately 16 ppm for deltaCSA11 and approximately 22 ppm for deltaCSA22 and deltaCSA33. We examine the consequences of the CSA variations for predicting the chemical shift changes upon partial alignment deltacsa and for the values of CSA order parameters extracted from the analysis of 31P NMR relaxation data. The theoretical 31P CSA tensors as well as the experimental 31P CSA tensor of barium diethyl phosphate (BDEP) are used to calculate deltacsa for two eclipsed orientations of the CSA and molecular alignment tensors. Percentage differences between the CSA order parameters obtained using the theoretical 31P CSA tensors and the experimental 31P CSA tensor of BDEP, respectively, are also determined.     

Disiline-doped boron nitride nanotubes: A computational study

The properties of the electronic structure of the Disiline-doped boron nitride nanotubes (Disiline-BNNTs) are investigated by a density functional theory (DFT) calculation. The structural forms are firstly optimized and the CS tensors calculated. Subsequently, the chemical-shielding isotropic (CS^I) and chemical shielding anisotropic (CS^A) parameters are found. The shielding values of boron (B) and nitrogen (N) atoms were calculated by Gauge-Including Atomic Orbital (GIAO), Continuous Set of Gauge Transformations (CSGT) and Individual Gauges for Atoms in Molecules (IGAIM) methods, using B3LYP/6-311+G*. The B3LYP level of theory with IGAIM was the best method to evaluate the theoretical chemical shifts for studied models. The results reveal a significant effect of Disiline doping on the chemical shielding tensors at the sites of those ¹¹B and ¹⁵N nuclei located in the nearest neighborhood of the Disiline-doped ring. Furthermore, the values of dipole moments and HOMO-LUMO gaps change in the Disiline-doped models in comparison with the original pristine model.     

Study of potassium O,O'-Dibutyldithiophosphate combining DFT, 31P CP/MAS NMR and infrared spectroscopy

Dithiophosphates are used in many different industrial applications. To explain their functions and properties in these applications, a fundamental understanding on a molecular level is needed. Potassium O, O'-Dibutyldithiophosphate and its anion have been investigated by means of a combination of DFT and (31)P CP/MAS NMR and infrared spectroscopy. Several low-energy conformations were studied by DFT. Three different conformations with significantly different torsion angles of the O-C bond relative to the O-P-O plane were selected for further studies of infrared frequencies and (31)P NMR chemical-shift tensors. A good agreement between theoretical and experimental results was obtained, especially when the IR spectra or (31)P chemical shift tensor parameters of all three conformations were added, indicating that, because of the low energy difference between the conformations, the molecules are rapidly fluctuating between them.     

Investigation of the polymorphs of dimethyl-3,6-dichloro-2,5-dihydroxyterephthalate by (13)C solid-state NMR spectroscopy.

Two of the three conformational polymorphs of dimethyl-3,6-dichloro-2,5-dihydroxyterephthalate are studied by solid-state NMR techniques. The structural differences between the polymorphs have previously been studied by X-ray. In these two polymorphs named white and yellow due to their color, the major structural difference is the torsional angle between the ester group and the aromatic ring. The yellow form has a dihedral angle of 4 degrees between the plane of the aromatic ring and the plane of the ester group, while the white form has two different molecules per unit cell with dihedral angles of 70 degrees and 85 degrees. This change greatly affects the conjugation in the pi-electronic system. In addition, there are differences in the hydrogen-bonding patterns, with the white form having intermolecular hydrogen bonds and the yellow form having intramolecular hydrogen bonds. In this work, the carbon isotropic chemical shift values and the chlorine electric field gradient (EFG) tensor information are extracted from the (13)C MAS spectra, and the principal values of the chemical shift tensors of the carbons are obtained from 2D FIREMAT experiments. Quantum chemical calculations of the chemical shift tensor data as well as the EFG tensor are performed at the HF and DFT levels of theory on individual molecules and on stacks of three molecules to account for the important intermolecular interactions in the white form. The differences between the spectral data on the two polymorphs are discussed in terms of the known electronic and structural differences.     

13C carbonyl chemical shielding tensors: Comparing SCF,MBPT (2), and DFT predictions to experiment

In this work, we calculate the ¹³C nuclear magnetic resonance chemical shielding tensors for 18 carbonyl-containing compounds. The many-body perturbation theory (MBPT), self-consistent field (SCF), and density functional theory (DFT) formalisms were used with gauge including atomic orbitals (GIAO) to calculate the shielding tensors. Our data suggest that shielding tensors can be efficiently estimated by performing one MBPT(2) correlated calculation (e.g., at a reference geometry) and SCF-level calculations at other geometries and taking the SCF-to-correlated tensor element differences to be geometry independent. That is, the correlation contribution to the chemical shielding seems to be relatively constant over a considerable range of distortions. Treatment of correlation using DFT methods is shown to not be as systematically reliable as with MBPT(2). Data on 18 carbonyl compounds show that the single largest influence on the shielding tensor is the presence of nearby electron-withdrawing or electron-donating groups. Finally, although good agreement with powder or single-crystal experimental data is achieved for two or three tensor eigenvalues, systematic differences remain for one element; the origins of these differences are discussed. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 875–894, 1997     

Carbon-13 chemical shift tensors of disaccharides: measurement, computation and assignment

A recently developed chemical shift anisotropy amplification solid-state nuclear magnetic resonance (NMR) experiment is applied to the measurement of the chemical shift tensors in three disaccharides: sucrose, maltose, and trehalose. The measured tensor principal values are compared with those calculated from first principles using density functional theory within the planewave-pseudopotential approach. In addition, a method of assigning poorly dispersed NMR spectra, based on comparing experimental and calculated shift anisotropies as well as isotropic shifts, is demonstrated.