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.     

An experimental and theoretical study of the 13C and 31P chemical shielding tensors in solid O-phosphorylated amino acids

A series of l and dl forms of O-phosphorylated amino acids (serine, threonine, tyrosine) have been studied by using solid-state multinuclear NMR spectroscopy and ab initio calculations. Principal elements of the (13)C and (31)P chemical shielding tensors have been measured and discussed in relation to zwitterionic structures and intermolecular contacts. DFT calculations have been compared with experimental data showing their ability to reproduce experimentally obtained tensor values in this challenging class of compounds. The changes of orientation of (31)P chemical shielding tensor with respect to the molecular frame in the presence of hydrogen bonds have been revealed and discussed on the ground of theoretical calculations. The measurements of internuclear P...P distances, based on Zeeman magnetization exchange between (31)P spins with differing chemical shielding tensor orientations, were exploited for a clear distinction between enantiomers and racemates.     

Carbon-13 NMR shielding in the twenty common amino acids: comparisons with experimental results in proteins

We have used ab initio quantum chemical techniques to compute the (13)C(alpha) and (13)C(beta) shielding surfaces for the 14 amino acids not previously investigated (R. H. Havlin et al., J. Am. Chem. Soc. 1997, 119, 11951-11958) in their most popular conformations. The spans (Omega = sigma(33) - sigma(11)) of all the tensors reported here are large ( approximately 34 ppm) and there are only very minor differences between helical and sheet residues. This is in contrast to the previous report in which Val, Ile and Thr were reported to have large ( approximately 12 ppm) differences in Omega between helical and sheet geometries. Apparently, only the beta-branched (beta-disubstituted) amino acids have such large CSA span (Omega) differences; however, there are uniformly large differences in the solution-NMR-determined CSA (Deltasigma = sigma(orth) - sigma(par)) between helices and sheets in all amino acids considered. This effect is overwhelmingly due to a change in shielding tensor orientation. With the aid of such shielding tensor orientation information, we computed Deltasigma values for all of the amino acids in calmodulin/M13 and ubiquitin. For ubiquitin, we find only a 2.7 ppm rmsd between theory and experiment for Deltasigma over an approximately 45 ppm range, a 0.96 slope, and an R(2) = 0.94 value when using an average solution NMR structure. We also report C(beta) shielding tensor results for these same amino acids, which reflect the small isotropic chemical shift differences seen experimentally, together with similar C(beta) shielding tensor magnitudes and orientations. In addition, we describe the results of calculations of C(alpha), C(beta), C(gamma)1, C(gamma)2, and C(delta) shifts in the two isoleucine residues in bovine pancreatic trypsin inhibitor and the four isoleucines in a cytochrome c and demonstrate that the side chain chemical shifts are strongly influenced by chi(2) torsion angle effects. There is very good agreement between theory and experiment using either X-ray or average solution NMR structures. Overall, these results show that both C(alpha) backbone chemical shift anisotropy results as well as backbone and side chain (13)C isotropic shifts can now be predicted with good accuracy by using quantum chemical methods, which should facilitate solution structure determination/refinement using such shielding tensor surface information.     

Single-crystal studies of peptide prolyl and glycyl 15N shielding tensors

15N shielding tensors were determined for the central peptide groups in GGV, AGG, and APG by single-crystal NMR. We find that the angle between the downfield component (delta11) and the N-H or the N-C(delta) (pro) bonds is in the range of 20-23 degrees and in accord with previous solid-state NMR measurements. However, AGG, unlike APG or GGV, has a distorted peptide plane, and delta11 lies approximately in the plane of N, C(alpha), and H rather than in the peptide plane defined by heavy atoms. Accurate orientations of delta22 and delta33 were determined, and the usual assumption that delta22 is along the peptide normal was found only in APG which has a highly nonaxial tensor. More generally, delta22 and delta33 are rotated about the delta11 axis (36 degrees in GGV). These results are compared with DFT calculations to gain a structural understanding of the effects of intermolecular interactions on shielding tensor principal components and orientations. Trimeric clusters containing H-bonded neighbors predict the orientations of the principal components within 2-3 degrees, but calculated principal components are less quantitative. Possible reasons for this disagreement are explored.     

Experimental and theoretical investigation of the 13C and 15N chemical shift tensors in melanostatin-exploring the chemical shift tensor as a structural probe

The determination of backbone conformations in powdered peptides using 13C and 15N shift tensor information is explored. The 13C and 15N principal shift values in natural abundance 13C and 15N melanostatin (L-Pro-L-Leu-Gly amide) are measured using the FIREMAT technique. Furthermore, the orientation of the C-N bond in the 13C shift principal axis system for the backbone carbons is obtained from the presence of the 13C-14N dipolar coupling. The Ramachandran angles for the title compound are obtained from solid-state NMR data by comparing the experimentally determined shift tensor information to systematic theoretical shielding calculations on N-formyl-L-amino acid-amide models. The effects of geometry optimization and neglect of intermolecular interactions on the theoretical shielding values in the model compounds are investigated. The sets of NMR derived Ramachandran angles are assembled in a set of test structures that are compared to the available single-crystal X-ray structure. Shift tensor calculations on the test structures and the X-ray structure are used to further assess the importance of intermolecular interactions when the shift tensor is used as a structural probe in powdered peptides.     

A solid-state 95Mo NMR and computational investigation of dodecahedral and square antiprismatic octacyanomolybdate(IV) anions: is the point-charge approximation an accurate probe of local symmetry?

Solid-state 95Mo NMR spectroscopy is shown to be an efficient and effective tool for analyzing the diamagnetic octacyanomolybdate(IV) anions, Mo(CN)(8)4-, of approximate dodecahedral, D(2d), and square antiprismatic, D(4d), symmetry. The sensitivity of the Mo magnetic shielding (sigma) and electric field gradient (EFG) tensors to small changes in the local structure of these anions allows the approximate D(2d) and D(4d) Mo(CN)(8)4- anions to be readily distinguished. The use of high applied magnetic fields, 11.75, 17.63 and 21.1 T, amplifies the overall sensitivity of the NMR experiment and enables more accurate characterization of the Mo sigma and EFG tensors. Although the magnitudes of the Mo sigma and EFG interactions are comparable for the D(2d) and D(4d) Mo(CN)(8)4- anions, the relative values and orientations of the principal components of the Mo sigma and EFG tensors give rise to 95Mo NMR line shapes that are significantly different at the fields utilized here. Quantum chemical calculations of the Mo sigma and EFG tensors, using zeroth-order regular approximation density functional theory (ZORA DFT) and restricted Hartree-Fock (RHF) methods, have also been carried out and are in good agreement with experiment. The most significant and surprising result from the DFT and RHF calculations is a significant EFG at Mo for an isolated Mo(CN)(8)4- anion possessing an ideal square antiprismatic structure; this is contrary to the point-charge approximation, PCA, which predicts a zero EFG at Mo for this structure.     

Calculations of the contribution of ring currents to the chemical shielding anisotropy

Ring currents can exert a large effect upon the chemical shielding of NMR resonances. The analytical expression developed by Waugh and Fessenden (Waugh, J. S.; Fessenden, R. W. J. Am. Chem. Soc. 1957, 79, 846) and Johnson and Bovey (Johnson, C. E.; Bovey, F. A. J. Chem. Phys. 1957, 29, 1012) only quantifies the contribution of ring currents to the isotropic component of the shielding tensor. In the work described here an additional analytical expression is developed so that the contribution of ring currents to the full shielding tensor can be calculated, allowing an estimate of the influence of ring currents upon the chemical shielding anisotropy (CSA, Deltasigma). To test that this pair of analytical expressions can provide a reasonable estimate of the contribution of ring currents to the full shielding tensor a series of density functional calculations (DFT) were carried out. A shielding tensor in a model compound was calculated in two distinct ways. For the first series, DFT shielding calculations of the model compound were carried out in the presence of a benzene ring. For the second series a ring current contribution to the shielding tensor was calculated via the new expressions, and this was added to the result of a DFT shielding calculation which used in place of benzene the nonaromatic analogue 1,3 cyclohexadiene. The two series of results proved to be in excellent agreement. The pair of analytical expressions are used to calculate ring current contributions to the CSA (Deltasigma) of 1H(N) backbone amide resonances in a structure of the second type 2 module from the protein fibronectin. Significant CSA variations are predicted in particular for the 1H(N) of G42 which is most likely involved in a N-H...tpi aromatic hydrogen bond.     

Secondary structures of peptides and proteins via NMR chemical-shielding anisotropy (CSA) parameters

Complete nuclear magnetic resonance (NMR) chemical-shielding tensors, sigma, have been computed at different levels of density-functional theory (DFT), within the gauge-including atomic orbital (GIAO) formalism, for the atoms of the peptide model For-L-Ala-NH2 as a function of the backbone dihedral angles phi and psi by employing a dense grid of 10 degrees. A complete set of rigorously orthogonal symmetric tensor invariants, {sigma iso, rho, tau}, is introduced, where sigma iso is the usual isotropic chemical shielding, while the newly introduced rho and tau parameters describe the magnitude and the orientation/shape of the chemical-shielding anisotropy (CSA), respectively. The set {sigma iso, rho, tau} is unaffected by unitary transformations of the symmetric part of the shielding tensor. The mathematically and physically motivated {rho, tau} anisotropy pair is easily connected to more traditional shielding anisotropy measures, like span (Omega) and skew (kappa). The effectiveness of the different partitions of the CSA information in predicting conformations of peptides and proteins has been tested throughout the Ramachandran space by generating theoretical NMR anisotropy surfaces for our For-L-Ala-NH2 model. The CSA surfaces, including Omega(phi, psi), kappa(phi, psi), rho(phi, psi), and tau(phi, psi) are highly structured. Individually, none of these surfaces is able to distinguish unequivocally between the alpha-helix and beta-strand secondary structural types of proteins. However, two- and three-dimensional correlated plots, including Omega versus kappa, rho versus tau, and sigma iso versus rho versus tau, especially for 13Calpha, have considerable promise in distinguishing among all four of the major secondary structural elements.     

51V solid-state magic angle spinning NMR spectroscopy and DFT studies of oxovanadium(V) complexes mimicking the active site of vanadium haloperoxidases

A series of 11 oxovanadium(V) complexes mimicking the active site of vanadium haloperoxidases have been investigated by (51)V magic angle spinning NMR spectroscopy and density functional theory (DFT). The MAS spectra are dominated by the anisotropic quadrupolar and chemical shielding interactions; for these compounds, C(Q) ranges from 3 to 8 MHz, and delta(sigma) is in the range 340-730 ppm. The quadrupolar coupling and chemical shielding tensors as well as their relative orientations have been determined by numerical simulations of the spectra. The spectroscopic NMR observables appear to be very sensitive to the details of the electronic and geometric environment of the vanadium center in these complexes. For the four crystallographically characterized compounds from the series, the quadrupolar and chemical shielding anisotropies were computed at the DFT level using two different basis sets, and the calculated tensors were in general agreement with the experimental solid-state NMR data. A combination of (51)V solid-state NMR and computational methods is thus beneficial for investigation of the electrostatic and geometric environment in diamagnetic vanadium systems with moderate quadrupolar anisotropies.     

Proton magnetic shielding tensor in liquid water

The nuclear magnetic shielding tensor is a sensitive probe of the local electronic environment, providing information about molecular structure and intermolecular interactions. The magnetic shielding tensor of the water proton has been determined in hexagonal ice, but in liquid water, where the tensor is isotropically averaged by rapid molecular tumbling, only the trace of the tensor has been measured. We report here the first determination of the proton shielding anisotropy in liquid water, which, when combined with chemical shift data, yields the principal shielding components parallel (sigma(parallel)) and perpendicular (sigma(perpendicular)) to the O-H bond. We obtained the shielding anisotropy sigma(parallel)-sigma(perpendicular) by measuring the proton spin relaxation rate as a function of magnetic induction field in a water sample where dipole-dipole couplings are suppressed by H/D isotope dilution. The temperature dependence of the shielding components, determined from 0 to 80 degrees C, reflects vibrational averaging over a distribution of instantaneous hydrogen-bond geometries in the liquid and thus contains unique information about the temperature-dependent structure of liquid water. The temperature dependence of the shielding anisotropy is found to be 4 times stronger than that of the isotropic shielding. We analyze the liquid water shielding components in the light of previous NMR and theoretical results for vapor and ice. We show that a simple two-state model of water structure fails to give a consistent interpretation of the shielding data and we argue that a more detailed analysis is needed that quantitatively relates the shielding components to hydrogen bond geometry.     

Peptide 17O chemical shielding and electric field gradient tensors

Complete (17)O chemical shielding (CS) and quadrupole coupling (QC) tensors and their molecular orientations were determined for the central residues in two tripeptides Gly-Gly-Val (GGV) and Ala-Gly-Gly (AGG) by single-crystal NMR methods. Tensor orientations in the two peptides are very similar, however, principal components are different. The most shielded CS and smallest magnitude QC components are normal to the peptide plane, while the most deshielded CS and largest QC components are in the peptide plane either at an angle of 17 degrees (CS) or perpendicular (QC) to the C=O bond. Comparisons of principal components from experiment and DFT calculations indicate that the smaller shielding tensor span in GGV (549 ppm) compared to AGG (606 ppm) is likely due to two factors: a shorter "direct" H-bond distance to the peptide carbonyl oxygen and an "indirect" H bond of the peptide NH to a carboxylate rather than a carbonyl. We anticipate that (17)O NMR should be generally useful for probing H-bonding and local electrostatic interactions in proteins and polypeptides. Using the single-crystal data as an accurate reference, we show that a useful subset of the NMR parameters, QC and CS principal components and their relative orientation, can be obtained with reasonable accuracy from a very high-field (21.2 T), stationary sample powder spectrum.     

13C and (15)N chemical shift tensors in adenosine,guanosine dihydrate, 2'-deoxythymidine,and cytidine

The (13)C and (15)N chemical shift tensor principal values for adenosine, guanosine dihydrate, 2'-deoxythymidine, and cytidine are measured on natural abundance samples. Additionally, the (13)C and (15)N chemical shielding tensor principal values in these four nucleosides are calculated utilizing various theoretical approaches. Embedded ion method (EIM) calculations improve significantly the precision with which the experimental principal values are reproduced over calculations on the corresponding isolated molecules with proton-optimized geometries. The (13)C and (15)N chemical shift tensor orientations are reliably assigned in the molecular frames of the nucleosides based upon chemical shielding tensor calculations employing the EIM. The differences between principal values obtained in EIM calculations and in calculations on isolated molecules with proton positions optimized inside a point charge array are used to estimate the contributions to chemical shielding arising from intermolecular interactions. Moreover, the (13)C and (15)N chemical shift tensor orientations and principal values correlate with the molecular structure and the crystallographic environment for the nucleosides and agree with data obtained previously for related compounds. The effects of variations in certain EIM parameters on the accuracy of the shielding tensor calculations are investigated.     

Binding energies and 19F nuclear magnetic deshielding in paramagnetic halogen-bonded complexes of TEMPO with haloperfluorocarbons

19F NMR measurements and theoretical calculations were performed to study paramagnetic complexes of iodoperfluorocarbons with stable nitroxide radicals. Contrary to what is usually measured for diamagnetic halogen-bonded complexes involving iodoperfluorocarbons, it was found that the formation of complexes with the 2,2,6,6-tetramethyl(piperidin-1-yloxyl) (TEMPO) radical determines downfield shifts in the 19F NMR spectra. The experimental finding was confirmed by calculating nuclear shielding using density functional theory and correcting the isotropic diamagnetic (19)F chemical shift with contact interactions evaluated from the hyperfine coupling tensor. The computational analysis of the interaction between CF3I and TEMPO, by using DFT and MP2 theories, showed that the occurrence of the halogen bond between the interacting partners is associated with a significant charge transfer to CF3I and that the measured downfield shift is due to the occurring spin transfer.     

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     

Understanding sterol-membrane interactions part I: Hartree-Fock versus DFT calculations of 13C and 1H NMR isotropic chemical shifts of sterols in solution and analysis of hydrogen-bonding effects

1H and 13C NMR chemical shifts are exquisitely sensitive probes of the local environment of the corresponding nuclei. Ultimately, direct determination of the chemical shifts of sterols in their membrane environment has the potential to reveal their molecular interactions and dynamics, in particular concerning the hydrogen-bonding partners of their OH groups. However, this strategy requires an accurate and efficient means to quantify the influence of the various interactions on chemical shielding. Herein the validity of Hartree-Fock and DFT calculations of the 13C and 1H NMR chemical shifts of cholesterol and ergosterol are compared with one another and with experimental chemical shifts measured in solution at 500 MHz. A computational strategy (definition of basis set, simpler molecular models for the sterols themselves and their molecular complexes) is proposed and compared with experimental data in solution. It is shown in particular that the effects of hydrogen bonding with various functional groups (water as a hydrogen-bond donor and acceptor, acetone) on NMR chemical shifts in CDCl3 solution can be accurately reproduced with this computational approach.     

The isotropic nuclear magnetic shielding constants of acetone in supercritical water: a sequential Monte Carlo/quantum mechanics study including solute polarization

The nuclear isotropic shielding constants sigma((17)O) and sigma((13)C) of the carbonyl bond of acetone in water at supercritical (P=340.2 atm and T=673 K) and normal water conditions have been studied theoretically using Monte Carlo simulation and quantum mechanics calculations based on the B3LYP6-311++G(2d,2p) method. Statistically uncorrelated configurations have been obtained from Monte Carlo simulations with unpolarized and in-solution polarized solute. The results show that solvent effects on the shielding constants have a significant contribution of the electrostatic interactions and that quantitative estimates for solvent shifts of shielding constants can be obtained modeling the water molecules by point charges (electrostatic embedding). In supercritical water, there is a decrease in the magnitude of sigma((13)C) but a sizable increase in the magnitude of sigma((17)O) when compared with the results obtained in normal water. It is found that the influence of the solute polarization is mild in the supercritical regime but it is particularly important for sigma((17)O) in normal water and its shielding effect reflects the increase in the average number of hydrogen bonds between acetone and water. Changing the solvent environment from normal to supercritical water condition, the B3LYP6-311++G(2d,2p) calculations on the statistically uncorrelated configurations sampled from the Monte Carlo simulation give a (13)C chemical shift of 11.7+/-0.6 ppm for polarized acetone in good agreement with the experimentally inferred result of 9-11 ppm.     

95Mo nuclear magnetic resonance parameters of molybdenum hexacarbonyl from density functional theory: appraisal of computational and geometrical parameters

Solid-state (95)Mo nuclear magnetic resonance (NMR) properties of molybdenum hexacarbonyl have been computed using density functional theory (DFT) based methods. Both quadrupolar coupling and chemical shift parameters were evaluated and compared with parameters of high precision determined using single-crystal (95)Mo NMR experiments. Within a molecular approach, the effects of major computational parameters, i.e. basis set, exchange-correlation functional, treatment of relativity, have been evaluated. Except for the isotropic parameter of both chemical shift and chemical shielding, computed NMR parameters are more sensitive to geometrical variations than computational details. Relativistic effects do not play a crucial part in the calculations of such parameters for the 4d transition metal, in particular isotropic chemical shift. Periodic DFT calculations were tackled to measure the influence of neighbouring molecules on the crystal structure. These effects have to be taken into account to compute accurate solid-state (95)Mo NMR parameters even for such an inorganic molecular compound.     

Carbon chemical shift tensor components in quinolines and quinoline N-oxides

Chemical shift calculations are carried out for the quinoline carbons in 1,8-bis(2-isopropyl-4-quinolyl)naphthalene, 2-isopropylquinoline, amodiaquine, chloroquine, and quinine and the N-oxide of each compound. Ab initio calculations of the isotropic shielding values are in agreement with experimental chemical shifts. The calculations indicate that changes to the principal components of the shielding tensor upon N-oxidation are similar for each compound. Carbons 2, 4, 8, and 10 are largely shielded in each case as the nitrogen is oxidized. For C2, C4, and C10, this shielding is due to a large change in sigma11 and/or sigma22, indicating a change in pi-electron density. For C8, the large shielding change is due mainly to a change in sigma33, indicating a change in sigma-electron density. Upon examination and comparison of the calculated 13C shielding tensor components in the antimalarial drugs versus those in unsubstituted quinolines, it is found that amodiaquine and chloroquine have increased pi-electron density in the ring containing the amino side chain and quinine has increased pi-electron density in the opposite ring, containing the methoxy substituent.     

Investigating the vanadium environments in hydroxylamido V(V) dipicolinate complexes using 51V NMR spectroscopy and density functional theory

Using (51)V magic angle spinning solid-state NMR, SSNMR, spectroscopy and quantum chemical DFT calculations we have characterized the chemical shift and quadrupolar coupling parameters of a series of eight hydroxylamido vanadium(V) dipicolinate complexes of the general formula VO(dipic)(ONR1R2)(H2O) where R1 and R2 can be H, CH3, or CH2CH3. This class of vanadium compounds was chosen for investigation because of their seven-coordinate vanadium atom, a geometry for which there is limited (51)V SSNMR data. Furthermore, a systematic series of compounds with different electronic properties are available and allows for the effects of ligand substitution on the NMR parameters to be studied. The quadrupolar coupling constants, C(Q), are small, 3.0-3.9 MHz, but exhibit variations as a function of the ligand substitution. The chemical shift tensors in the solid state are sensitive to changes in both the hydroxylamide substituent and the dipic ligand, a sensitivity which is not observed for isotropic chemical shifts in solution. The chemical shift tensors span approximately 1000 ppm and are nearly axially symmetric. On the basis of DFT calculations of the chemical shift tensors, one of the largest contributors to the magnetic shielding anisotropy is an occupied molecular orbital with significant vanadium d(z)2 character along the V=O bond.     

Density functional theory investigation of hydrogen bonding effects on the oxygen, nitrogen and hydrogen electric field gradient and chemical shielding tensors of anhydrous chitosan crystalline structure

A systematic computational investigation was carried out to characterize the 17O, 14N and 2H electric field gradient, EFG, as well as 17O, 15N, 13C and 1H chemical shielding tensors in the anhydrous chitosan crystalline structure. To include the hydrogen-bonding effects in the calculations, the most probable interacting molecules with the target molecule in the crystalline phase were considered through a hexameric cluster. The computations were performed with the B3LYP method and 6-311++G(d,p) and 6-31++G(d,p) standard basis sets using the Gaussian 98 suite of programs. Calculated EFG and chemical shielding tensors were used to evaluate the 17O, 14N and 2H nuclear quadrupole resonance, NQR, and 17O, 15N, 13C and 1H nuclear magnetic resonance, NMR, parameters in the hexameric cluster, which are in good agreement with the available experimental data. The difference between the calculated NQR and NMR parameters of the monomer and hexamer cluster shows how much hydrogen bonding interactions affect the EFG and chemical shielding tensors of each nucleus. These results indicate that both O(3)-H(33)...O(5-3) and N-H(22)...O(6-4) hydrogen bonding have a major influence on NQR and NMR parameters. Also, the quantum chemical calculations indicate that the intra- and intermolecular hydrogen bonding interactions play an essential role in determining the relative orientation of EFG and chemical shielding principal components in the molecular frame axes.