Calculation of binary magnetic properties and potential energy curve in xenon dimer: second virial coefficient of (129)Xe nuclear shielding

Quantum chemical calculations of the nuclear shielding tensor, the nuclear quadrupole coupling tensor, and the spin-rotation tensor are reported for the Xe dimer using ab initio quantum chemical methods. The binary chemical shift delta, the anisotropy of the shielding tensor Delta sigma, the nuclear quadrupole coupling tensor component along the internuclear axis chi( parallel ), and the spin-rotation constant C( perpendicular ) are presented as a function of internuclear distance. The basis set superposition error is approximately corrected for by using the counterpoise correction (CP) method. Electron correlation effects are systematically studied via the Hartree-Fock, complete active space self-consistent field, second-order M?ller-Plesset many-body perturbation, and coupled-cluster singles and doubles (CCSD) theories, the last one without and with noniterative triples, at the nonrelativistic all-electron level. We also report a high-quality theoretical interatomic potential for the Xe dimer, gained using the relativistic effective potential/core polarization potential scheme. These calculations used valence basis set of cc-pVQZ quality supplemented with a set of midbond functions. The second virial coefficient of Xe nuclear shielding, which is probably the experimentally best-characterized intermolecular interaction effect in nuclear magnetic resonance spectroscopy, is computed as a function of temperature, and compared to experiment and earlier theoretical results. The best results for the second virial coefficient, obtained using the CCSD(CP) binary chemical shift curve and either our best theoretical potential or the empirical potentials from the literature, are in good agreement with experiment. Zero-point vibrational corrections of delta, Delta sigma, chi (parallel), and C (perpendicular) in the nu=0, J=0 rovibrational ground state of the xenon dimer are also reported.     

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

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.     

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.     

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.     

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.     

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.     

Simulations of 129Xe NMR chemical shift of atomic xenon dissolved in liquid benzene

The isotropic ¹²⁹Xe NMR chemical shift of atomic Xe dissolved in liquid benzene was simulated by combining classical molecular dynamics and quantum chemical calculations of ¹²⁹Xe nuclear magnetic shielding. Snapshots from the molecular dynamics trajectory of xenon atom in a periodic box of benzene molecules were used for the quantum chemical calculations of isotropic ¹²⁹Xe chemical shift using nonrelativistic density functional theory as well as relativistic Breit?CPauli perturbation corrections. Thus, the correlation and relativistic effects as well as the temperature and dynamics effects could be included in the calculations. Theoretical results are in a very good agreement with the experimental data. The most of the experimentally observed isotropic ¹²⁹Xe shift was recovered in the nonrelativistic dynamical region, while the relativistic effects explain of about 8% of the total ¹²⁹Xe chemical shift.     

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.     

Noble gas-actinide complexes of the CUO molecule with multiple Ar,Kr, and Xe atoms in noble-gas matrices

Laser-ablated U atoms react with CO in excess argon to produce CUO, which is trapped in a triplet state in solid argon at 7 K, based on agreement between observed and relativistic density functional theory (DFT) calculated isotopic frequencies ((12)C(16)O, (13)C(16)O, (12)C(18)O). This observation contrasts a recent neon matrix investigation, which trapped CUO in a linear singlet state calculated to be about 1 kcal/mol lower in energy. Experiments with krypton and xenon give results analogous to those with argon. Similar work with dilute Kr and Xe in argon finds small frequency shifts in new four-band progressions for CUO in the same triplet states trapped in solid argon and provides evidence for four distinct CUO(Ar)(4-n)(Ng)(n) (Ng = Kr, Xe, n = 1, 2, 3, 4) complexes for each Ng. DFT calculations show that successively higher Ng complexes are responsible for the observed frequency progressions. This work provides the first evidence for noble gas-actinide complexes, and the first example of neutral complexes with four noble gas atoms bonded to one metal center.     

Rotational spectra of the Xe-(H2O)2 van der Waals trimer: xenon as a probe of electronic structure and dynamics

Rotational spectra of three isotopomers of the Xe-(H2O)2 van der Waals trimer were recorded using a pulsed-nozzle, Fourier transform microwave spectrometer. Nine [nine, four] a-type and twelve [eleven, seven] b-type transitions were measured for the 132Xe-(H2O)2 [129Xe-(H2O)2, 131Xe-(H2O)2] isotopomer. The determined rotational and centrifugal distortion constants were used to extract information about the structure and vibrational motions of the complex. The nuclear quadrupole hyperfine structures due to the 131Xe (nuclear spin quantum number I=3/2) nucleus were also detected. The large value of the off-diagonal nuclear quadrupole coupling constant chiab in particular provides detailed insight into the electronic environment of the xenon atom and the orientations of the water molecules within the complex. An effective structure that best reproduces the experimental 131Xe nuclear quadrupole coupling constants is rationalized by ab initio calculations. An overall goal of this line of work is to determine how the successive solvation of a xenon atom with water molecules affects the xenon electron distribution and its intermolecular interactions. The results may provide molecular level interpretations of 129Xe NMR data from, for example, imaging experiments.     

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.     

Solid-state 129Xe and 131Xe NMR study of the perxenate anion XeO(6)4-

Results of the first solid-state 131Xe NMR study of xenon-containing compounds are presented. The two NMR-active isotopes of xenon, 129Xe (I=1/2) and 131Xe (I=3/2), are exploited to characterize the xenon magnetic shielding and quadrupolar interactions for two sodium perxenate salts, Na4XeO6.xH2O (x=0, 2), at an applied magnetic field strength of 11.75 T. Solid-state 129/131Xe NMR line shapes indicate that the local xenon environment in anhydrous Na4XeO6 adopts octahedral symmetry, but upon hydration, the XeO6(4-) anion becomes noticeably distorted from octahedral symmetry. For stationary, anhydrous samples of Na4XeO6, the heteronuclear 129/131Xe-23Na dipolar interaction is the principal contributor to the breadth of the 129/131Xe NMR lines. For stationary and slow magic-angle-spinning samples of Na4XeO(6).2H2O, the anisotropic xenon shielding interaction dominates the 129Xe NMR line shape, whereas the 131Xe NMR line shape is completely dominated by the nuclear quadrupolar interaction. The xenon shielding tensor is approximately axially symmetric, with a skew of -0.7+/-0.3, an isotropic xenon chemical shift of -725.6+/-1.0 ppm, and a span of 95+/-5 ppm. The 131Xe quadrupolar coupling constant, 10.8+/-0.5 MHz, is large for a nucleus at a site of approximate Oh symmetry, and the quadrupolar asymmetry parameter indicates a lack of axial symmetry. This study demonstrates the extreme sensitivity of the 131Xe nuclear quadrupolar interaction to changes in the local xenon environment.     

129Xe NMR chemical shift in Xe@C60 calculated at experimental conditions: Essential role of the relativity,dynamics, and explicit solvent

The isotropic ¹²⁹Xe nuclear magnetic resonance (NMR) chemical shift (CS) in Xe@C₆₀ dissolved in liquid benzene was calculated by piecewise approximation to faithfully simulate the experimental conditions and to evaluate the role of different physical factors influencing the ¹²⁹Xe NMR CS. The ¹²⁹Xe shielding constant was obtained by averaging the ¹²⁹Xe nuclear magnetic shieldings calculated for snapshots obtained from the molecular dynamics trajectory of the Xe@C₆₀ system embedded in a periodic box of benzene molecules. Relativistic corrections were added at the Breit–Pauli perturbation theory (BPPT) level, included the solvent, and were dynamically averaged. It is demonstrated that the contribution of internal dynamics of the Xe@C₆₀ system represents about 8% of the total nonrelativistic NMR CS, whereas the effects of dynamical solvent add another 8%. The dynamically averaged relativistic effects contribute by 9% to the total calculated ¹²⁹Xe NMR CS. The final theoretical value of 172.7 ppm corresponds well to the experimental ¹²⁹Xe CS of 179.2 ppm and lies within the estimated errors of the model. The presented computational protocol serves as a prototype for calculations of ¹²⁹Xe NMR parameters in different Xe atom guest–host systems. © 2013 Wiley Periodicals, Inc.     

Calculations of nuclear quadrupole coupling in noble gas-noble metal fluorides: interplay of relativistic and electron correlation effects

The nuclear quadrupole coupling constants (NQCCs) of noble gas and noble metal nuclei in the recently found noble gas-noble metal fluorides (NgMF, where Ng=Ar,Kr,Xe and M=Cu,Ag,Au) are obtained theoretically by high-level ab initio calculations, where both relativistic and electron correlation effects are included, and compared to experimental results. Fully relativistic four-component Dirac-Hartree-Fock (DHF) calculations are carried out at the basis set limit for electric field gradient that couples with the electric quadrupole moment of the nucleus, and uncorrelated relativistic effects are extracted by comparing DHF results to nonrelativistic (NR) HF calculations. Electron correlation effects are investigated both at fully relativistic second-order Moller-Plesset (DMP2) and at NR MP2 levels of theory, as well as at the NR coupled-cluster singles and doubles with perturbational triples [CCSD(T)] level. The validity of the approximation where relativistic effects, on the one hand, and nonrelativistically obtained correlation effects, on the other hand, are evaluated separately and assumed to be additive, is investigated by comparison with the DMP2 results. Inclusion of relativistic effects is shown to be necessary for obtaining the correct NQCC trends as the nucleus of interest and/or its neighbors become heavier. Electron correlation treatment is needed for approaching quantitative agreement with the experimental NQCCs. The assumption of additive electron correlation and relativistic effects, corresponding to the NR correlation treatment added on top of relativistic DHF data, gives qualitatively correct noble gas NQCCs. For noble metal NQCCs, correlation treatment at the relativistic level of theory is mandatory for reaching agreement with experimental results. Current work also confirms the experimental trends of NQCCs, which have been taken as an indication of nearly covalent interaction between noble gas and noble metal in the heaviest present systems, especially in XeAuF.     

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.     

129Xe and 131Xe nuclear magnetic dipole moments from gas phase NMR spectra

³He, ¹²⁹Xe and ¹³¹Xe NMR measurements of resonance frequencies in the magnetic field B₀ = 11.7586 T in different gas phase mixtures have been reported. Precise radiofrequency values were extrapolated to the zero gas pressure limit. These results combined with new quantum chemical values of helium and xenon nuclear magnetic shielding constants were used to determine new accurate nuclear magnetic moments of ¹²⁹Xe and ¹³¹Xe in terms of that of the ³He nucleus. They are as follows: μ(¹²⁹Xe) = ?0.7779607(158)μ_N and μ(¹³¹Xe) = +0.6918451(70)μ_N. By this means, the new ‘helium method’ for estimations of nuclear dipole moments was successfully tested. Gas phase NMR spectra demonstrate the weak intermolecular interactions observed on the ³He and ¹²⁹Xe and ¹³¹Xe shielding in the gaseous mixtures with Xe, CO₂ and SF₆. Copyright © 2015 John Wiley & Sons, Ltd.     

Probing polymer colloids by Xe NMR

Model aqueous dispersions of polystyrene, poly(methyl methacrylate), poly(n-butyl acrylate) and a statistical copolymer poly(n-butyl acrylate-co-methyl methacrylate) were studied using xenon NMR spectroscopy. The ¹²⁹Xe NMR spectra of these various latexes reveal qualitative and quantitative differences in the number of peaks and in their line widths and chemical shifts. Above the glass transition temperature, exchange between xenon sorbed in the particle core and free xenon outside the particles is fast on the ¹²⁹Xe spectral time-scale and a single ¹²⁹Xe signal is observed. At temperatures below the glass transition temperature, the exchange between sorbed and free xenon is slow on the ¹²⁹Xe spectral time-scale and two ¹²⁹Xe NMR signals can be observed. If the signal of sorbed ¹²⁹Xe is observed, its chemical shift, line width and integral relative to the integral of free ¹²⁹Xe can be used for the characterization of the particle core. The line width of free ¹²⁹Xe provides the residence time of xenon outside the particles and can be used to determine the rate constant characterizing the kinetics of penetration of xenon in the particles. This rate constant emerges as promising parameter for the characterization of the polymer particle surface.