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.     

Ab initio study of the ground state properties of molecular oxygen

The electric and magnetic properties of the ground state of oxygen molecule are calculated by multiconfiguration self-consisted field (MCSCF) method and compared with experimental data: the quadrupole moment, polarizability, the 17O nuclear quadrupole coupling constant, magnetizability tensor, nuclear spin-rotation coupling constant and rotational g factor. The last two constants are calculated for all possible isotope modifications. The rotational, ESR and NMR spectra are discussed. Fermi-contact hyperfine coupling parameter is additionally estimated by different methods. The NMR chemical shielding tensor for 17O16O species at high temperature limit (without electron spin contribution) is predicted. Potential energy curves for 10 excited bound states and the internuclear distance dependence of the studied properties are also presented.     

The nuclear magnetic resonance line shapes of Xe in the cages of clathrate hydrates

We report, for the first time, a prediction of the line shapes that would be observed in the (129)Xe nuclear magnetic resonance (NMR) spectrum of xenon in the cages of clathrate hydrates. We use the dimer tensor model to represent pairwise contributions to the intermolecular magnetic shielding tensor for Xe at a specific location in a clathrate cage. The individual tensor components from quantum mechanical calculations in clathrate hydrate structure I are represented by contributions from parallel and perpendicular tensor components of Xe-O and Xe-H dimers. Subsequently these dimer tensor components are used to reconstruct the full magnetic shielding tensor for Xe at an arbitrary location in a clathrate cage. The reconstructed tensors are employed in canonical Monte Carlo simulations to find the Xe shielding tensor component along a particular magnetic field direction. The shielding tensor component weighted according to the probability of finding a crystal fragment oriented along this direction in a polycrystalline sample leads to a predicted line shape. Using the same set of Xe-O and Xe-H shielding functions and the same Xe-O and Xe-H potential functions we calculate the Xe NMR spectra of Xe atom in 12 distinct cage types in clathrate hydrates structures I, II, H, and bromine hydrate. Agreement with experimental spectra in terms of the number of unique tensor components and their relative magnitudes is excellent. Agreement with absolute magnitudes of chemical shifts relative to free Xe atom is very good. We predict the Xe line shapes in two cages in which Xe has not yet been observed.     

Xe chemical shift tensor in silicalite and SSZ-24

We report, for the first time, a theoretical prediction of the (129)Xe nuclear magnetic resonance chemical shift tensor of xenon atom in a single crystal of silicalite at near-zero occupancy and the temperature dependence of the Xe NMR chemical shift tensor for the polycrystalline silicalite at maximum occupancy. The former is a measure of the sensitivity of the Xe tensor components to the local structure of the channels without Xe-Xe contributions. The latter is a measure of the sensitivity of the Xe-Xe tensor components to the Xe-Xe distributions, as determined by the Xe-Xe potential function in competition with the Xe-silicalite potential function. Both theoretical predictions can be compared against Xe NMR experiments: the first against the Xe spectra collected as a function of rotation of the single crystal about the three crystalline axes in a magnetic field, and the second against variable temperature Xe NMR studies (below room temperature) of polycrystalline silicalite at maximum Xe occupancy. With the same parameter set (Xe-O potential and shielding functions), we predict the line shapes of Xe in SSZ-24 zeolite under various conditions of occupancy and temperature.     

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.     

Theoretical study of the nuclear spin-molecular rotation coupling for relativistic electrons and non-relativistic nuclei

A theoretical study of the relation between the relativistic formulation of the nuclear magnetic shielding and spin-rotation tensors is presented. To this end a theoretical expression of the relativistic spin-rotation tensor is formulated, considering a molecular Hamiltonian of relativistic electrons and non-relativistic nuclei. Molecular rotation effects are introduced considering the terms of the Born-Oppenheimer decomposition, which couple the electrons and nuclei dynamics. The loss of the simple relation linking both spectral parameters in the non-relativistic formulation is further analyzed carrying out a perturbative expansion of relativistic effects by means of the linear response within the elimination of the small component approach. It is concluded that relativistic effects on the spin-rotation tensor are less important than those of the nuclear magnetic shielding tensor.     

Theoretical predictions of nuclear magnetic resonance parameters in a novel organo-xenon species: chemical shifts and nuclear quadrupole couplings in HXeCCH

We calibrate the methodology for the calculation of nuclear magnetic resonance (NMR) properties in novel organo-xenon compounds. The available state-of-the-art quantum-chemical approaches are combined and applied to the HXeCCH molecule as the model system. The studied properties are (129)Xe, (1)H, and (13)C chemical shifts and shielding anisotropies, as well as (131)Xe and (2)H nuclear quadrupole coupling constants. The aim is to obtain, as accurately as currently possible, converged results with respect to the basis set, electron correlation, and relativistic effects, including the coupling of relativity and correlation. This is done, on one hand, by nonrelativistic correlated ab initio calculations up to the CCSD(T) level and, on the other hand, for chemical shifts and shielding anisotropies by the leading-order relativistic Breit-Pauli perturbation theory (BPPT) with correlated ab initio and density-functional theory (DFT) reference states. BPPT at the uncorrelated Hartree-Fock level as well as the corresponding fully relativistic Dirac-Hartree-Fock method are found to be inapplicable due to a dramatic overestimation of relativistic effects, implying the influence of triplet instability in this multiply bonded system. In contrast, the fully relativistic second-order Moller-Plesset perturbation theory method can be applied for the quadrupole coupling, which is a ground-state electric property. The performance of DFT with various exchange-correlation functionals is found to be inadequate for the nonrelativistic shifts and shielding anisotropies as compared to the CCSD(T) results. The relativistic BPPT corrections to these quantities can, however, be reasonably predicted by DFT, due to the improved triplet excitation spectrum as compared to the Hartree-Fock method, as well as error cancellation within the five main BPPT contributions. We establish three computationally feasible models with characteristic error margins for future calculations of larger organo-xenon compounds to guide forthcoming experimental NMR efforts. The predicted (129)Xe chemical shift in HXeCCH is in a novel range for this nucleus, between weakly bonded or solvated atomic xenon and xenon in the hitherto characterized molecules.     

Spin‐orbit ZORA and four‐component Dirac–Coulomb estimation of relativistic corrections to isotropic nuclear shieldings and chemical shifts of noble gas dimers

Hartree–Fock and density functional theory with the hybrid B3LYP and general gradient KT2 exchange‐correlation functionals were used for nonrelativistic and relativistic nuclear magnetic shielding calculations of helium, neon, argon, krypton, and xenon dimers and free atoms. Relativistic corrections were calculated with the scalar and spin‐orbit zeroth‐order regular approximation Hamiltonian in combination with the large Slater‐type basis set QZ4P as well as with the four‐component Dirac–Coulomb Hamiltonian using Dyall's acv4z basis sets. The relativistic corrections to the nuclear magnetic shieldings and chemical shifts are combined with nonrelativistic coupled cluster singles and doubles with noniterative triple excitations [CCSD(T)] calculations using the very large polarization‐consistent basis sets aug‐pcSseg‐4 for He, Ne and Ar, aug‐pcSseg‐3 for Kr, and the AQZP basis set for Xe. For the dimers also, zero‐point vibrational (ZPV) corrections are obtained at the CCSD(T) level with the same basis sets were added. Best estimates of the dimer chemical shifts are generated from these nuclear magnetic shieldings and the relative importance of electron correlation, ZPV, and relativistic corrections for the shieldings and chemical shifts is analyzed. © 2015 Wiley Periodicals, Inc.     

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.     

Nuclear spin relaxation due to chemical shift anisotropy of gas-phase 129Xe

Nuclear spin relaxation provides detailed dynamical information on molecular systems and materials. Here, first-principles modeling of the chemical shift anisotropy (CSA) relaxation time for the prototypic monoatomic (129)Xe gas is carried out, both complementing and predicting the results of NMR measurements. Our approach is based on molecular dynamics simulations combined with pre-parametrized ab initio binary nuclear shielding tensors, an "NMR force field". By using the Redfield relaxation formalism, the simulated CSA time correlation functions lead to spectral density functions that, for the first time, quantitatively determine the experimental spin-lattice relaxation times T(1). The quality requirements on both the Xe-Xe interaction potential and binary shielding tensor are investigated in the context of CSA T(1). Persistent dimers Xe(2) are found to be responsible for the CSA relaxation mechanism in the low-density limit of the gas, completely in line with the earlier experimental findings.     

The chemical shifts of Xe in the cages of clathrate hydrate Structures I and II

We report, for the first time, a calculation of the isotropic NMR chemical shift of 129Xe in the cages of clathrate hydrates Structures I and II. We generate a shielding surface for Xe in the clathrate cages by quantum mechanical calculations. Subsequently this shielding surface is employed in canonical Monte Carlo simulations to find the average isotropic Xe shielding values in the various cages. For the two types of cages in clathrate hydrate Structure I, we find the intermolecular shielding values [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-214.0 ppm, and [sigma(Xe@5(12)6(2) cage)-sigma(Xe atom)]=-146.9 ppm, in reasonable agreement with the values -242 and -152 ppm, respectively, observed experimentally by Ripmeester and co-workers between 263 and 293 K. For the 5(12) and 5(12)6(4) cages of Structure II we find [sigma(Xe@5(12) cage)-sigma(Xe atom)]=-206.7 ppm, and [sigma(Xe@5(12)6(4) cage)-sigma(Xe atom)]=-104.7 ppm, also in reasonable agreement with the values -225 and -80 ppm, respectively, measured in a Xe-propane type II mixed clathrate hydrate at 77 and 220-240 K by Ripmeester et al.     

Hyperfine spectrum of RbCl

The molecular beam electric resonance technique has been used to conduct a high precision examination of the hyperfine spectrum of the four isotopomers of RbCl. Coupling constants for the nuclear electric quadrupole interactions, the spin-rotation interactions, the tensor and scalar spin-spin interactions, and a rubidium nuclear octupole interaction, and their dependence on vibrational and rotational states have been determined. The dominant interaction, the rubidium nuclear electric quadrupole interaction, shows a small shift with substitution of the chlorine isotope.     

Grand canonical Monte Carlo simulations of the 129Xe NMR line shapes of xenon adsorbed in (+/-)-[Co(en)3]Cl3

The 129Xe NMR line shapes of xenon adsorbed in the nanochannels of the (+/-)-[Co(en)3]Cl3 ionic crystal have been calculated by grand canonical Monte Carlo (GCMC) simulations. The results of our GCMC simulations illustrate their utility in predicting 129Xe NMR chemical shifts in systems containing a transition metal. In particular, the nanochannels of (+/-)-[Co(en)3]Cl3 provide a simple, yet interesting, model system that serves as a building block toward understanding xenon chemical shifts in more complex porous materials containing transition metals. Using only the Xe-C and Xe-H potentials and shielding response functions derived from the Xe@CH4 van der Waals complex to model the interior of the channel, the GCMC simulations correctly predict the 129Xe NMR line shapes observed experimentally (Ueda, T.; Eguchi, T.; Nakamura, N.; Wasylishen, R. E. J. Phys. Chem. B 2003, 107, 180-185). At low xenon loading, the simulated 129Xe NMR line shape is axially symmetric with chemical-shift tensor components delta(parallel) = 379 ppm and delta(perpendicular) = 274 ppm. Although the simulated isotropic chemical shift, delta(iso) = 309 ppm, is overestimated, the anisotropy of the chemical-shift tensor is correctly predicted. The simulations provide an explanation for the observed trend in the 129Xe NMR line shapes as a function of the overhead xenon pressure: delta(perpendicular) increased from 274 to 292 ppm, while delta(parallel) changed by only 3 ppm over the entire xenon loading range. The overestimation of the isotropic chemical shifts is explained based upon the results of quantum mechanical 129Xe shielding calculations of xenon interacting with an isolated (+/-)-[Co(en)3]Cl3 molecule. The xenon chemical shift is shown to be reduced by about 12% going from the Xe@[Co(en)3]Cl3 van der Waals complex to the Xe@C2H6 fragment.     

C CP MAS NMR and DFT study of vascular-selective drugs felodipine and amlodipine

Amlodipine besylate (AM) and felodipine (FL) have been studied in solid by ¹³C CP MAS NMR and DFT. The spectra have been successfully reproduced from the theoretical calculations of the nuclear magnetic shielding tensors, using the GIAO method, after the SCF partial geometry optimisation. A very good correlation between the ¹³C chemical shifts and quadrupole coupling constants for chlorobenzene, AM and FL was obtained. ¹³C NMR as well as ³⁵Cl NQR do not differentiate between S and R enantiomers of AM and FL due to the symmetry of the 4-aryl ring comprising one (AM) or two (FL) chlorine atoms oriented in a perpendicular fashion over the 1,4-DHP ring.     

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.     

Theoretical predictions of the spectroscopic parameters in noble-gas molecules: HXeOH and its complex with water

We employ state-of-the-art methods and basis sets to study the effect of inserting the Xe atom into the water molecule and the water dimer on their NMR parameters. Our aim is to obtain predictions for the future experimental investigation of novel xenon complexes by NMR spectroscopy. Properties such as molecular structure and energetics have been studied by supermolecular approaches using HF, MP2, CCSD, CCSD(T) and MP4 methods. The bonding in HXeOH···H(2)O complexes has been analyzed by Symmetry-Adapted Perturbation Theory to provide the intricate insight into the nature of the interaction. We focus on vibrational spectra, NMR shielding and spin-spin coupling constants-experimental signals that reflect the electronic structures of the compounds. The parameters have been calculated at electron-correlated and Dirac-Hartree-Fock relativistic levels. This study has elucidated that the insertion of the Xe atom greatly modifies the NMR properties, including both the electron correlation and relativistic effects, the (129)Xe shielding constants decrease in HXeOH and HXeOH···H(2)O in comparison to Xe atom; the (17)O, as a neighbour of Xe, is deshielded too. The HXeOH···H(2)O complex in its most stable form is stabilized mainly by induction and dispersion energies.     

Solid-State Nuclear Magnetic Resonance Techniques for the Structural Characterization of Geminal Alane-Phosphane Frustrated Lewis Pairs and Secondary Adducts

The first comprehensive solid-state nuclear magnetic resonance (NMR) characterization of geminal alane-phosphane frustrated Lewis pairs (Al/P FLPs) is reported. Their relevant NMR parameters (isotropic chemical shifts, direct and indirect ²⁷Al-³¹P spin-spin coupling constants, and ²⁷Al nuclear electric quadrupole coupling tensor components) have been determined by numerical analysis of the experimental NMR line shapes and compared with values computed from the known crystal structures by using density functional theory (DFT) methods. Our work demonstrates that the ³¹P NMR chemical shifts for the studied Al/P FLPs are very sensitive to slight structural inequivalences. The ²⁷Al NMR central transition signals are spread out over a broad frequency range (>200 kHz), owing to the presence of strong nuclear electric quadrupolar interactions that can be well-reproduced by the static ²⁷Al wideband uniform rate smooth truncation (WURST) Carr-Purcell-Meiboom-Gill (WCPMG) NMR experiment. ²⁷Al chemical shifts and quadrupole tensor components offer a facile and clear distinction between three- and four-coordinate aluminum environments. For measuring internuclear Al⋅⋅⋅P distances a new resonance-echo saturation-pulse double-resonance (RESPDOR) experiment was developed by using efficient saturation via frequency-swept WURST pulses. The successful implementation of this widely applicable technique indicates that internuclear Al⋅⋅⋅P distances in these compounds can be measured within a precision of ±0.1 Å.     

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.