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.     

Relativistic effects on group-12 metal nuclear shieldings

The leading-order perturbation theory approach to relativistic effects on the nuclear magnetic shielding provides an economic method for obtaining the chemical shifts in heavy-element containing systems. The method features detailed analysis potential in terms of the different physical mechanisms affecting the shielding tensors of heavy nuclei. The perturbative nature, however, results in an increasing error with increasingly heavy elements in the system. In this work, we investigate the performance of the Breit-Pauli perturbation theory (BPPT) against fully relativistic four-component theory in computing the nuclear shielding constants as well as the chemical shifts with respect to corresponding atomic ions of group-12 metals, M = Zn, Cd, and Hg, in dimethyl M(CH(3))(2) and aqueous M(H(2)O)(6)(2+) complexes. It is shown that five out of the total of sixteen BPPT correction terms are responsible for most of the relativistic corrections for the chemical shift of studied metals. The relativity is important already for Cd and BPPT is proven to work well up to Hg for the chemical shift, as calibrated with the fully relativistic method.     

The Xe shielding surfaces for Xe interacting with linear molecules and spherical tops

The 129Xe nuclear magnetic resonance spectrum of xenon in gas mixtures of Xe with other molecules provides a test of the ab initio surfaces for the intermolecular shielding of Xe in the presence of the other molecule. We examine the electron correlation contributions to the Xe-CO2, Xe-N2, Xe-CO, Xe-CH4, and Xe-CF4 shielding surfaces and test the calculations against the experimental temperature dependence of the density coefficients of the Xe chemical shift in the gas mixtures at infinite dilution in Xe. Comparisons with the gas phase data permit the refinement of site-site potential functions for Xe-N2, Xe-CO, and Xe-CF4 especially for atom-Xe distances in the range 3.5-6 A. With the atom-atom shielding surfaces and potential parameters obtained in the present work, construction of shielding surfaces and potentials for applications such as molecular dynamics averaging of Xe chemical shifts in liquid solvents containing CH3, CH2, CF3, and CF2 groups is possible.     

DFT study of the NMR properties of xenon in covalent compounds and van der waals complexes-implications for the use of 129Xe as a molecular probe

The NMR properties (chemical shift and spin-spin coupling constants) of (129)Xe in covalent compounds and weakly bound complexes have been investigated by DFT methods including relativistic effects. For covalent species, a good agreement between experimental and calculated results is achieved without scalar relativistic effects, but their inclusion (with a triple-zeta, double-polarization basis set) leads to some improvement in the quality of the correlation. The spin-orbit coupling term has a significant effect on the shielding constant, but makes a small contribution to the chemical shift. Coupling constants contain substantial contributions from the Fermi contact and paramagnetic spin-orbit terms; unlike light nuclei the spin-dipole term is also large, whereas the diamagnetic spin-orbit term is negligible. For van der Waals dimers, the dependence of the xenon chemical shift and anisotropy is calculated as a function of the distance. Small (<1 Hz) but non-negligible through-space coupling constants between (129)Xe and (13)C or (1)H are predicted. Much larger couplings, of the order of few Hz, are calculated between xenon and (17)O in a model silicate residue.     

Effect of relativity on the ionization spectra of the xenon fluorides XeFn (n=2, 4, 6)

Noble gas compounds exhibit special chemical bonding situations and have been investigated by various spectroscopic and theoretical techniques. In this work we calculate the ionization spectra of the xenon fluorides (XeF2,XeF4, and XeF6) in the valence and subvalence (down to Xe 4d) areas by application of the recently developed Dirac-Hartree-Fock one-particle propagator technique. In this technique, the relativistic (four-component) and electron correlation effects are computed simultaneously. The xenon compounds show considerable spin-orbit splitting strongly influencing the photoelectron spectrum not reproducible in prior calculations. Comparison to one-component methods is made and the occurring satellite structures are interpreted. The satellite structures can be attributed either to the breakdown of the one-particle picture or to a reflection of intra-atomic and interatomic Auger decay processes within the molecule.     

Exploring new 129Xe chemical shift ranges in HXeY compounds: hydrogen more relativistic than xenon

Among rare gases, xenon features an unusually broad nuclear magnetic resonance (NMR) chemical shift range in its compounds and as a non-bonded Xe atom introduced into different environments. In this work we show that (129)Xe NMR chemical shifts in the recently prepared, matrix-isolated xenon compounds appear in new, so far unexplored (129)Xe chemical shift ranges. State-of-the-art theoretical predictions of NMR chemical shifts in compounds of general formula HXeY (Y = H, F, Cl, Br, I, -CN, -NC, -CCH, -CCCCH, -CCCN, -CCXeH, -OXeH, -OH, -SH) as well as in the recently prepared ClXeCN and ClXeNC species are reported. The bonding situation of Xe in the studied compounds is rather different from the previously characterized cases as Xe appears in the electronic state corresponding to a situation with a low formal oxidation state, between I and II in these compounds. Accordingly, the predicted (129)Xe chemical shifts occur in new NMR ranges for this nucleus: ca. 500-1000 ppm (wrt Xe gas) for HXeY species and ca. 1100-1600 ppm for ClXeCN and ClXeNC. These new ranges fall between those corresponding to the weakly-bonded Xe(0) atom in guest-host systems (δ < 300 ppm) and in the hitherto characterized Xe molecules (δ > 2000 ppm). The importance of relativistic effects is discussed. Relativistic effects only slightly modulate the (129)Xe chemical shift that is obtained already at the nonrelativistic CCSD(T) level. In contrast, spin-orbit-induced shielding effects on the (1)H chemical shifts of the H1 atom directly bonded to the Xe center largely overwhelm the nonrelativistic deshielding effects. This leads to an overall negative (1)H chemical shift in the range between -5 and -25 ppm (wrt CH(4)). Thus, the relativistic effects induced by the heavy Xe atom appear considerably more important for the chemical shift of the neighbouring, light hydrogen atom than that of the Xe nucleus itself. The predicted NMR parameters facilitate an unambiguous experimental identification of these novel compounds.     

Relativistic heavy-atom effects on heavy-atom nuclear shieldings

The principal relativistic heavy-atom effects on the nuclear magnetic resonance (NMR) shielding tensor of the heavy atom itself (HAHA effects) are calculated using ab initio methods at the level of the Breit-Pauli Hamiltonian. This is the first systematic study of the main HAHA effects on nuclear shielding and chemical shift by perturbational relativistic approach. The dependence of the HAHA effects on the chemical environment of the heavy atom is investigated for the closed-shell X(2+), X(4+), XH(2), and XH(3) (-) (X=Si-Pb) as well as X(3+), XH(3), and XF(3) (X=P-Bi) systems. Fully relativistic Dirac-Hartree-Fock calculations are carried out for comparison. It is necessary in the Breit-Pauli approach to include the second-order magnetic-field-dependent spin-orbit (SO) shielding contribution as it is the larger SO term in XH(3) (-), XH(3), and XF(3), and is equally large in XH(2) as the conventional, third-order field-independent spin-orbit contribution. Considering the chemical shift, the third-order SO mechanism contributes two-thirds of the difference of approximately 1500 ppm between BiH(3) and BiF(3). The second-order SO mechanism and the numerically largest relativistic effect, which arises from the cross-term contribution of the Fermi contact hyperfine interaction and the relativistically modified spin-Zeeman interaction (FC/SZ-KE), are isotropic and practically independent of electron correlation effects as well as the chemical environment of the heavy atom. The third-order SO terms depend on these factors and contribute both to heavy-atom shielding anisotropy and NMR chemical shifts. While a qualitative picture of heavy-atom chemical shifts is already obtained at the nonrelativistic level of theory, reliable shifts may be expected after including the third-order SO contributions only, especially when calculations are carried out at correlated level. The FC/SZ-KE contribution to shielding is almost completely produced in the s orbitals of the heavy atom, with values diminishing with the principal quantum number. The relative contributions converge to universal fractions for the core and subvalence ns shells. The valence shell contribution is negligible, which explains the HAHA characteristics of the FC/SZ-KE term. Although the nonrelativistic theory gives correct chemical shift trends in present systems, the third-order SO-I terms are necessary for more reliable predictions. All of the presently considered relativistic corrections provide significant HAHA contributions to absolute shielding in heavy atoms.     

Solvent effects and dynamic averaging of 195Pt NMR shielding in cisplatin derivatives

The influences of solvent effects and dynamic averaging on the (195)Pt NMR shielding and chemical shifts of cisplatin and three cisplatin derivatives in aqueous solution were computed using explicit and implicit solvation models. Within the density functional theory framework, these simulations were carried out by combining ab initio molecular dynamics (aiMD) simulations for the phase space sampling with all-electron relativistic NMR shielding tensor calculations using the zeroth-order regular approximation. Structural analyses support the presence of a solvent-assisted "inverse" or "anionic" hydration previously observed in similar square-planar transition-metal complexes. Comparisons with computationally less demanding implicit solvent models show that error cancellation is ubiquitous when dealing with liquid-state NMR simulations. After aiMD averaging, the calculated chemical shifts for the four complexes are in good agreement with experiment, with relative deviations between theory and experiment of about 5% on average (1% of the Pt(II) chemical shift range).     

NMR shielding constants in PH3, absolute shielding scale, and the nuclear magnetic moment of 31P

Ab initio values of the absolute shielding constants of phosphorus and hydrogen in PH(3) were determined, and their accuracy is discussed. In particular, we analyzed the relativistic corrections to nuclear magnetic resonance (NMR) shielding constants, comparing the constants computed using the four-component Dirac-Hartree-Fock approach, the four-component density functional theory (DFT), and the Breit-Pauli perturbation theory (BPPT) with nonrelativistic Hartree-Fock or DFT reference functions. For the equilibrium geometry, we obtained σ(P) = 624.309 ppm and σ(H) = 29.761 ppm. Resonance frequencies of both nuclei were measured in gas-phase NMR experiments, and the results were extrapolated to zero density to provide the frequency ratio for an isolated PH(3) molecule. This ratio, together with the computed shielding constants, was used to determine a new value of the nuclear magnetic dipole moment of (31)P: μ(P) = 1.1309246(50) μ(N).     

Ab initio calculations of NMR chemical shifts

The nuclear magnetic resonance chemical shift is one of the most powerful properties available for structure determination at the molecular level. A review of advances made in the ab initio calculation of chemical shielding during the past five years is presented. Specifically, progress in the areas including the effects of an unpaired electron, electron correlation, and relativistic effects into ab initio chemical shielding calculations, the tensor nature of the chemical shift, and intramolecular and intermolecular effects on the chemical shift will be covered.     

High-resolution electron spin resonance spectroscopy of XeF* in solid argon. The hyperfine structure constants as a probe of relativistic effects in the chemical bonding properties of a heavy noble gas atom

Xenon fluoride radicals were generated by solid-state chemical reactions of mobile fluorine atoms with xenon atoms trapped in Ar matrix. Highly resolved electron spin resonance spectra of XeF* were obtained in the temperature range of 5-25 K and the anisotropic hyperfine parameters were determined for magnetic nuclei 19F, 129Xe, and 131Xe using naturally occurring and isotopically enriched xenon. Signs of parallel and perpendicular hyperfine components were established from analysis of temperature changes in the spectra and from numerical solutions of the spin Hamiltonian for two nonequivalent magnetic nuclei. Thus, the complete set of components of hyperfine- and g-factor tensors of XeF* were obtained: 19F (Aiso=435, Adip=1249 MHz) and 129Xe (Aiso=-1340, Adip=-485 MHz); g(parallel)=1.9822 and g(perpendicular)=2.0570. Comparison of the measured hyperfine parameters with those predicted by density-functional theory (DFT) calculations indicates, that relativistic DFT gives true electron spin distribution in the 2Sigma+ ground-state, whereas nonrelativistic theory underestimates dramatically the electron-nuclear contact Fermi interaction (Aiso) on the Xe atom. Analysis of the obtained magnetic-dipole interaction constants (Adip) shows that fluorine 2p and xenon 5p atomic orbitals make a major contribution to the spin density distribution in XeF*. Both relativistic and nonrelativistic calculations give close magnetic-dipole interaction constants, which are in agreement with the measured values. The other relativistic feature is considerable anisotropy of g-tensor, which results from spin-orbit interaction. The orbital contribution appears due to mixing of the ionic 2Pi states with the 2Sigma+ ground state, and the spin-orbit interaction plays a significant role in the chemical bonding of XeF*.     

Leading-order relativistic effects on nuclear magnetic resonance shielding tensors

We present perturbational ab initio calculations of the nuclear-spin-dependent relativistic corrections to the nuclear magnetic resonance shielding tensors that constitute, together with the other relativistic terms reported by us earlier, the full leading-order perturbational set of results for the one-electron relativistic contributions to this observable, based on the (Breit-)Pauli Hamiltonian. These contributions are considered for the H(2)X (X = O,S,Se,Te,Po) and HX (X = F,Cl,Br,I,At) molecules, as well as the noble gas (Ne, Ar, Kr, Xe, Rn) atoms. The corrections are evaluated using the relativistic and magnetic operators as perturbations on an equal footing, calculated using analytical linear and quadratic response theory applied on top of a nonrelativistic reference state provided by self-consistent field calculations. The (1)H and heavy-atom nuclear magnetic shielding tensors are compared with four component, nearly basis-set-limit Dirac-Hartree-Fock calculations that include positronic excitations, as well as available literature data. Besides the easy interpretability of the different contributions in terms of familiar nonrelativistic concepts, the accuracy of the present perturbational scheme is striking for the isotropic part of the shielding tensor, for systems including elements up to Xe.     

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

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.     

Are ring currents still useful to rationalize the benzene proton magnetic shielding?

[structure: see text] The conventional interpretation of proton NMR chemical shifts is supported by large basis set ab initio quantum mechanical calculations. The benzene protons are predicted to lie within the deshielding zone defined in terms of the out-of-plane magnetic shielding domain. However, ring currents by themselves are not sufficient to account quantitatively for the observed benzene proton downfield chemical shift. sigma-Electron contributions must also be taken into account. The conventional explanation for the ethyne proton chemical shift is valid.     

Heats of formation of xenon fluorides and the fluxionality of XeF(6) from high level electronic structure calculations

Atomization energies at 0 K and heats of formation at 0 and 298 K are predicted for XeF(+), XeF(-), XeF(2), XeF(4), XeF(5)(-), and XeF(6) from coupled cluster theory (CCSD(T)) calculations with new correlation-consistent basis sets for Xe. To achieve near chemical accuracy (+/-1 kcal/mol), up to four corrections were added to the complete basis set binding energies based on frozen core coupled cluster theory energies: a correction for core-valence effects, a correction for scalar relativistic effects, a correction for first-order atomic spin-orbit effects, and in some cases, a second-order spin-orbit correction. Vibrational zero-point energies were computed at the coupled cluster level of theory. The structure of XeF(6) is difficult to obtain with the C(3)(v)() and O(h)() structures having essentially the same energy. The O(h)() structure is only 0.19 kcal/mol below the C(3)(v)() one at the CCSD(T)/CBS level using an approximate geometry for the C(3)(v)() structure. With an optimized C(3)(v)() geometry, the C(3)(v)() structure would probably become slightly lower in energy than the O(h)() one. The calculated heats of formation for the neutral XeF(n)() fluorides are less negative than the experimental values from the equilibrium measurements by 2.0, 7.7, and 12.2 kcal/mol for n = 2, 4, and 6, respectively. For the experimental values, derived from the photoionization measurements, this discrepancy becomes even larger, suggesting a need for a redetermination of the experimental values. Evidence is presented for the fluxionality of XeF(6) caused by the presence of a sterically active, free valence electron pair on Xe.     

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.     

Rapid calculation of protein chemical shifts using bond polarization theory and its application to protein structure refinement

Although difficult to analyze, NMR chemical shifts provide detailed information on protein structure. We have adapted the semi-empirical bond polarization theory (BPT) to protein chemical shift calculation and chemical shift driven protein structure refinement. A new parameterization for BPT amide nitrogen chemical shift calculation has been derived from MP2 ab initio calculations and successfully evaluated using crystalline tripeptides. We computed the chemical shifts of the small globular protein ubiquitin, demonstrating that BPT calculations can match the results obtained at the DFT level of theory at very low computational cost. In addition to the calculation of chemical shift tensors, BPT allows the calculation of chemical shift gradients and consequently chemical shift driven geometry optimizations. We applied chemical shift driven protein structure refinement to the conformational analysis of a set of Trypanosoma brucei (the causative agent of African sleeping sickness) tryparedoxin peroxidase Px III structures. We found that the interaction of Px III with its reaction partner Tpx seems to be governed by conformational selection rather than by induced fit.     

Novel method for the measurement of xenon gas solubility using 129Xe NMR spectroscopy

A novel method is presented for determining xenon partitioning between a gas phase and a liquid phase. An experimental setup which permits the simultaneous measurement of the 129Xe chemical shift in both the gas and the liquid phases, that is, under the same experimental conditions, has been designed. Xenon solubility is obtained via 129Xe chemical shift measurements in the gas phase. The method was validated against xenon solubility data from the literature; in general, the agreement is found to be within 3%. The solubility of xenon in three solvents for which data have not been previously reported (acetone, acetonitrile, and 1,1,2,2-tetrachloroethane) was determined using this novel method. 129Xe chemical shifts for dissolved xenon are also reported; it is found that xenon-xenon interactions may play a significant role in the liquid phase even at low equilibrium xenon pressures.