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.     

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.     

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.     

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.     

Theoretical study of hyperfine interactions in small arsenic-containing radicals

Various density functional theory (DFT) and ab initio MP2 and CCSD methods were employed in calculations of arsenic isotropic and anisotropic hyperfine parameters in three small radicals, AsH(2), AsO(2), and H(2)AsO. Convergent basis sets for these calculations were specially derived starting from Dunning's correlation-consistent sets. DFT methods proved to be particularly appropriate choice, because the results obtained with the suitable functionals are in accordance with the available experimental values. Additionally, mechanisms of hyperfine couplings were investigated by examining individual orbital contributions to spin density. The spin polarization mechanism in AsH(2) was studied by evaluating one- and two-electron integrals in spin-restricted and unrestricted case. Apart from nonrelativistic, scalar-relativistic Douglas-Kroll-Hess DFT calculations were performed to examine relativistic impacts on spin density distribution.     

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.     

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.     

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.     

Relativistic spin-orbit effects on hyperfine coupling tensors by density-functional theory

A second-order perturbation theory treatment of spin-orbit corrections to hyperfine coupling tensors has been implemented within a density-functional framework. The method uses the all-electron atomic mean-field approximation and/or spin-orbit pseudopotentials in incorporating one- and two-electron spin-orbit interaction within a first-principles framework. Validation of the approach on a set of main-group radicals and transition metal complexes indicates good agreement between all-electron and pseudopotential results for hyperfine coupling constants of the lighter nuclei in the system, except for cases in which scalar relativistic effects become important. The nonrelativistic Fermi contact part of the isotropic hyperfine coupling constants is not always accurately reproduced by the exchange-correlation functionals employed, particularly for the triplet and pi-type doublet radicals in the present work. For this reason, ab initio coupled-cluster singles and doubles with perturbative triples results for the first-order contributions have been combined in the validation calculations with the density-functional results for the second-order spin-orbit contributions. In the cases where spin-orbit corrections are of significant magnitude relative to the nonrelativistic first-order terms, they improve the agreement with experiment. Antisymmetric contributions to the hyperfine tensor arise from the spin-orbit contributions and are discussed for the IO2 radical, whereas rovibrational effects have been evaluated for RhC, NBr, and NI.     

Electric field effects on the shielding constants of noble gases: a four-component relativistic Hartree-Fock study

Second derivatives of nuclear shielding constants with respect to an electric field, i.e., shielding polarizabilities, have been calculated for the noble gas atoms from helium to xenon. The calculations have been carried out using the four-component relativistic Hartree-Fock method. In order to assess the importance of the individual relativistic corrections, the shielding polarizabilities have also been calculated at the nonrelativistic Hartree-Fock level, with spin-orbit and scalar (Darwin and mass-velocity) effects having been established by perturbative methods. Electron correlation effects have been estimated using the second-order polarization propagator approach. The relativistic effects on the tensor components of the shielding polarizabilities are found to be larger and changing less regularly with the atomic number than for the shielding constant itself. However, there is a partial cancellation of the contributions to the parallel and perpendicular components of the shielding polarizability and as a consequence the mean shielding polarizability is far less affected than the individual components.     

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.     

A double quantum (129)Xe NMR experiment for probing xenon in multiply-occupied cavities of solid-state inclusion compounds

A method is presented for detecting multiple xenon atoms in cavities of solid-state inclusion compounds using (129)Xe double quantum NMR spectroscopy. Double quantum filtered (129)Xe NMR spectra, performed on the xenon clathrate of Dianin's compound were obtained under high-resolution Magic-Angle Spinning (MAS) conditions, by recoupling the weak (129)Xe-(129)Xe dipole-dipole couplings that exist between xenon atoms in close spatial proximity. Because the (129)Xe-(129)Xe dipole-dipole couplings are generally weak due to dynamics of the atoms and to large internuclear separations, and since the (129)Xe Chemical Shift Anisotropy (CSA) tends to be relatively large, a very robust dipolar recoupling sequence was necessary, with the symmetry-based SR26 dipolar recoupling sequence proving appropriate. We have also attempted to measure the (129)Xe-(129)Xe dipole-dipole coupling constant between xenon atoms in the cavities of the xenon-Dianin's compound clathrate and have found that the dynamics of the xenon atoms (as investigated with molecular dynamics simulations) as well as (129)Xe multiple spin effects complicate the analysis. The double quantum NMR method is useful for peak assignment in (129)Xe NMR spectra because peaks arising from different types of absorption/inclusion sites or from different levels of occupancy of single sites can be distinguished. The method can also help resolve ambiguities in diffraction experiments concerning the order/disorder in a material.     

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.     

The dependence on temperature and concentration of 129Xe NMR chemical shift values of organoxenon(II) salts [RXe][Y] and of the parent molecule XeF2

The dependence of the 129Xe NMR chemical shift value of XeF2 on temperature and concentration was determined in a variety of prototypic media: in acidic (anhydrous HF, aHF), nonprotic but polar (dichloromethane), and basic (CD3CN-EtCN, 1:3 v/v) solvents. The 129Xe NMR spectra of a representative series of organoxenon(II) salts [RXe][Y] (R = C6F5, heptafluoro-1,4-cyclohexadien-1-yl (cyclo-1,4-C6F7), pentafluoro-1,4-cyclohexadien-3-on-1-yl (cyclo-3-O-1,4-C6F5), CF2=C(CF3), (CF3)2CFC[triple bond]C, C4H9C[triple bond]C; Y = BF4, AsF6) in aHF showed, compared with XeF2-aHF, a quantitatively less distinct but qualitatively related dependence of delta(129Xe) vs temperature. The dependence of their delta(129Xe) values on concentration in aHF is negligible. An explanation for the different behavior of [RXe][Y] and XeF2 is offered.     

Density-functional calculations of relativistic spin-orbit effects on nuclear magnetic shielding in paramagnetic molecules

Terms arising from the relativistic spin-orbit effect on both hyperfine and Zeeman interactions are introduced to density-functional theory calculation of nuclear magnetic shielding in paramagnetic molecules. The theory is a generalization of the former nonrelativistic formulation for doublet systems and is consistent to O(alpha4), the fourth power of the fine structure constant, for the spin-orbit terms. The new temperature-dependent terms arise from the deviation of the electronic g tensor from the free-electron g value as well as spin-orbit corrections to hyperfine coupling tensor A, the latter introduced in the present work. In particular, the new contributions include a redefined isotropic pseudocontact contribution that consists of effects due to both the g tensor and spin-orbit corrections to hyperfine coupling. The implementation of the spin-orbit terms makes use of all-electron atomic mean-field operators and/or spin-orbit pseudopotentials. Sample results are given for group-9 metallocenes and a nitroxide radical. The new O(alpha4) corrections are found significant for the metallocene systems while they obtain small values for the nitroxide radical. For the isotropic shifts, none of the three beyond-leading-order hyperfine contributions are negligible.     

Interpretation of indirect nuclear spin-spin coupling tensors for polyatomic xenon fluorides and group 17 fluorides: results from relativistic density-functional calculations

Significant improvements have been made recently in the calculation of NMR indirect nuclear spin-spin coupling tensors (J). In particular, the relativistic zeroth-order regular approximation density-functional theory (ZORA-DFT) approach holds great promise for the calculation of spin-spin coupling constants for a variety of chemical systems containing heavy nuclei. In the present work, the ZORA-DFT method is applied to the calculation of the complete reduced coupling tensors, K, for a range of chlorine-, bromine-, iodine-, and xenon-containing species: K(Cl,F) for ClF(2)(+), ClF(3), ClF(4)(+), ClF(5), ClF(6)(-), and ClF(6)(+); K(Br,F) for BrF(3), BrF(6)(-), and BrF(6)(+); K(I,F) for IF(4)(+) and IF(6)(+); K(Xe,F) for XeF(+), XeF(2), XeF(3)(+), XeF(4), XeF(5)(-), XeF(5)(+), and XeF(7)(+). These species represent a wide variety of geometrical bonding arrangements. Agreement between the calculated coupling constants and available experimental data is excellent, and the absolute sign of the coupling constants is provided. It is shown that (1)K(iso) may be positive or negative even within the same molecule, e.g., K(Cl,F)(iso) may be of either sign, depending on the local environment. Periodic trends in (1)K(iso) for isovalent and isostructural molecules are evident. The spin-spin coupling anisotropies, Delta K, and the orientations of the K tensors are also determined. The success of the calculations is a direct result of employing reliable geometries and considering both scalar and spin-orbit relativistic effects. The dependence of K(Cl,F)(iso) and K(Xe,F)(iso) on the local molecular and electronic structure is discussed in terms of the paramagnetic spin-orbit (PSO) and combined Fermi-contact spin-dipolar (FC+SD) coupling mechanisms. The PSO term depends strongly on the number of valence shell electron lone pairs on the central heavy atom, and the FC+SD contribution increases with the Cl[bond]F or Xe[bond]F bond length for a given series of compounds. This interpretation allows for the successful rationalization of the existing experimental data.     

Renner-Teller and spin-orbit vibronic coupling effects in linear triatomic molecules with a half-filled pi shell

The vibronic and spin-orbit-induced interactions among the (3)Sigma(-), (1)Delta, and (1)Sigma(+) electronic states arising from a half-filled pi orbital of a linear triatomic molecule are considered, employing the microscopic (Breit-Pauli) spin-orbit coupling operator. The 6 x 6 Hamiltonian matrix is derived in a diabatic spin-orbital electronic basis set, including terms up to fourth order in the expansion of the molecular Hamiltonian in the bending normal coordinate about the linear geometry. The symmetry properties of the Hamiltonian are analyzed. Aside from the nonrelativistic fourth-order Renner-Teller vibronic coupling within the (1)Delta state and the second-order nonrelativistic vibronic coupling between the (1)Sigma(+) and (1)Delta states, there exist zeroth-order, first-order, as well as third-order vibronic coupling terms of spin-orbit origin. The latter are absent when the phenomenological expression for the spin-orbit coupling operator is used instead of the microscopic form. The effects of the nonrelativistic and spin-orbit-induced vibronic coupling mechanisms on the (3)Sigma(-), (1)Delta, and (1)Sigma(+) adiabatic potential energy surfaces as well as on the spin-vibronic energy levels are discussed for selected parameter values.     

Performance of nonrelativistic and quasi-relativistic hybrid DFT for the prediction of electric and magnetic hyperfine parameters in 57Fe Mössbauer spectra

57Fe electric and magnetic hyperfine parameters were calculated for a series of 10 iron model complexes, covering a wide range of oxidation and spin states. Employing the B3LYP hybrid method, results from nonrelativistic density functional theory (DFT) and quasi-relativistic DFT within the zero-order regular approximation (ZORA) were compared. Electron densities at the iron nuclei were calculated and correlated with experimental isomer shifts. It was shown that the fit parameters do not depend on a specific training set of iron complexes and are, therefore, more universal than might be expected. The nonrelativistic and quasi-relativistic electron densities gave fit parameters of similar quality; the ZORA densities are only shifted by a factor of 1.32, upward in the direction of the four-component Dirac-Fock value. From a correlation of calculated electric field gradients and experimental quadrupole splittings, the value of the 57Fe nuclear quadrupole moment was redetermined to a value of 0.16 barn, in good agreement with other studies. The ZORA approach gave no additional improvement of the calculated quadrupole splittings in comparison to the nonrelativistic approach. The comparison of the calculated and measured 57Fe isotropic hyperfine coupling constants (hfcc's) revealed that both the ZORA approach and the inclusion of spin-orbit contributions lead to better agreement between theory and experiment in comparison to the nonrelativistic results. For all iron complexes with small spin-orbit contributions (high-spin ferric and ferryl systems), a distinct underestimation of the isotropic hfcc's was found. Scaling factors of 1.81 (nonrelativistic DFT) and 1.69 (ZORA) are suggested. The calculated 57Fe isotropic hfcc's of the remaining model systems (low-spin ferric and high-spin ferrous systems) contain 10-50% second-order contributions and were found to be in reasonable agreement with the experimental results. This is assumed to be the consequence of error cancellation because g-tensor calculations for these systems are of poor quality with the existing DFT approaches. Excellent agreement between theory and experiment was found for the 57Fe anisotropic hfcc's. Finally, all of the obtained fit parameters were used for an application study of the [Fe(H2O)6]3+ ion. The calculated spectroscopic data are in good agreement with the Mossbauer and electron paramagnetic resonance results discussed in detail in a forthcoming paper.