Full four-component relativistic calculations of NMR shielding and indirect spin-spin coupling tensors in hydrogen halides

On the geometry dependence of the nuclear magnetic resonance chemical shift of mercury in thiolate complexes: A relativistic density functional theory study

Thiolate containing mercury(II) complexes of the general formula [Hg(SR)n ]2 - n have been of great interest since the toxicity of mercury was recognized. 199Hg nuclear magnetic resonance spectroscopy (NMR) is a powerful tool for characterization of mercury complexes. In this work, the Hg shielding constants in a series of [Hg(SR)n ]2 - n complexes are therefore investigated computationally with particular emphasis on their geometry dependence. Geometry optimizations and NMR chemical shift calculations are performed at the density functional theory (DFT) level with both the zeroth-order regular approximation (ZORA) and four-component relativistic methods. The four exchange-correlation (XC) functionals PBE0, PBE, B3LYP, and BLYP are used in combination with either Dyall's Gaussian-type (GTO) or Slater-type orbitals (STOs) basis sets. Comparing ZORA and four-component calculations, one observes that the calculated shielding constants for a given molecular geometry have a constant difference of ∼ 1070 ppm. This confirms that ZORA is an acceptable relativistic method to compute NMR chemical shifts. The combinations of four-component/PBE0/v3z and ZORA/PBE0/QZ4P are applied to explore the geometry dependence of the isotropic shielding. For a given coordination number, the distance between mercury and sulfur is the key factor affecting the shielding constant, while changes in bond and dihedral angles and even different side groups have relatively little impact.

Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-Accelerated Computer Architectures

We present the development and implementation of relativistic coupled cluster linear response theory (CC-LR), which allows the determination of molecular properties arising from time-dependent or time-independent electric, magnetic, or mixed electric-magnetic perturbations (within a common gauge origin for the magnetic properties) as well as taking into account the finite lifetime of excited states in the framework of damped response theory. We showcase our implementation, which is capable to offload the computationally intensive tensor contractions characteristic of coupled cluster theory onto graphical processing units, in the calculation of (a) frequency-(in)dependent dipole-dipole polarizabilities of IIB atoms and selected diatomic molecules, with a particular emphasis on the calculation of valence absorption cross sections for the I2 molecule; (b) indirect spin-spin coupling constants for benchmark systems such as the hydrogen halides (HX, X = F-I) as well the H2Se-H2O dimer as a prototypical system containing hydrogen bonds; and (c) optical rotations at the sodium D line for hydrogen peroxide analogues (H2Y2, Y = O, S, Se, Te). Thanks to this implementation, we are able to show the similarities in performance, but often the significant discrepancies, between CC-LR and approximate methods such as density functional theory. Comparing standard CC response theory with the flavor based upon the equation of motion formalism, we find that for valence properties such as polarizabilities, the two frameworks yield very similar results across the periodic table as found elsewhere in the literature; for properties that probe the core region, such as spin-spin couplings, on the other hand, we show a progressive differentiation between the two as relativistic effects become more important. Our results also suggest that as one goes down the periodic table, it may become increasingly difficult to measure pure optical rotation at the sodium D line due to the appearance of absorbing states.

Analytical Second-Order Properties for the Random Phase Approximation: Nuclear Magnetic Resonance Shieldings

A method for the analytical computation of nuclear magnetic resonance (NMR) shieldings within the direct random phase approximation (RPA) is presented. As a starting point, we use the RPA ground-state energy expression within the resolution-of-the-identity approximation in the atomic-orbital formalism. As has been shown in a recent benchmark study using numerical second derivatives [Glasbrenner, M. J. Chem. Theory Comput. 2022, 18, 192], RPA based on a Hartree-Fock reference shows accuracies comparable to coupled cluster singles and doubles (CCSD) for NMR chemical shieldings. Together with the much lower computational cost of RPA, it has emerged as an accurate method for the computation of NMR shieldings. Therefore, we aim to extend the applicability of RPA NMR to larger systems by introducing analytical second-order derivatives, making it a viable method for the accurate and efficient computation of NMR chemical shieldings.

On the Utmost Importance of the Basis Set Choice for the Calculations of the Relativistic Corrections to NMR Shielding Constants

The investigation of the sensitivity of the relativistic corrections to the NMR shielding constants (σ) to the configuration of angular spaces of the basis sets used on the atoms of interest was carried out within the four-component density functional theory (DFT). Both types of relativistic effects were considered, namely the so-called heavy atom on light atom and heavy atom on heavy atom effects, though the main attention was paid to the former. As a main result, it was found that the dependence of the relativistic corrections to σ of light nuclei (exemplified here by 1H and 13C) located in close vicinity to a heavy atom (exemplified here by In, Sn, Sb, Te, and I) on the basis set used on the light spectator atom was very much in common with that of the Fermi-contact contribution to the corresponding nonrelativistic spin-spin coupling constant (J). In general, it has been shown that the nonrelativistic J-oriented and σ-oriented basis sets, artificially saturated in the tight s-region, provided much better accuracy than the standard nonrelativistic σ-oriented basis sets when calculating the relativistic corrections to the NMR shielding constants of light nuclei at the relativistic four-component level of the DFT theory.

Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy

The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.

Computational 199 Hg NMR

Theoretical background and fundamental results dealing with the computation of mercury chemical shifts and spin-spin coupling constants are reviewed with a special emphasis on their stereochemical behavior and applications.

Relativistic and QED corrections to one-bond indirect nuclear spin-spin couplings in X$_2^{2+}$ and X$_3^{2+}$ ions (X = Zn, Cd, Hg)

The indirect spin-spin coupling tensor, $\bm J$, between mercury nuclei in systems containing this element can be of the order of few kHz and one of the largest measured. We analysed the physics behind the electronic mechanisms that contribute to the one- and two-bond couplings $^n \bm{J}_{\mathrm{Hg}-\mathrm{Hg}}$ ($n=1, 2$). For doing so, we performed calculations for $J$-couplings in the ionized $X_2^{2+}$ and $X_3^{2+}$ linear molecules ($X$ = Zn, Cd, Hg) within polarization propagator theory, using the random phase approximation and the pure zeroth--order approximation with Dirac--Hartree--Fock and Dirac--Kohn--Sham orbitals, both at four-component and ZORA levels. We show that the ``paramagnetic-like' mechanism contribute with more than 99.98\% to the total isotropic value of the coupling tensor. By analyzing the molecular and atomic orbitals involved in the total value of the response function, we find that the $s$-type valence atomic orbitals have a predominant role in the description of the coupling. This fact allows us to develop an effective model from which QED effects on $J$-couplings in the aforementioned ions can be estimated. Those effects were found to be within the interval $(0.7;~1.7)$\% of the total relativistic effect on isotropic one-bond $^1\bm{J}$ coupling, though ranging those corrections between the interval $(-0.4;~-0.2)$\% in Zn-containing ions, to $(-1.2;~-0.8)$\% in Hg-containing ions, of the total isotropic coupling constant in the studied systems. The estimated QED corrections show a visible dependence on the nuclear charge $Z$ of each atom $X$ in the form of a power-law proportional to $Z_X^5$.

Quantum Chemical Approaches to the Calculation of NMR Parameters: From Fundamentals to Recent Advances

Quantum chemical methods for the calculation of indirect NMR spin–spin coupling constants and chemical shifts are always in progress. They never stay the same due to permanently developing computational facilities, which open new perspectives and create new challenges every now and then. This review starts from the fundamentals of the nonrelativistic and relativistic theory of nuclear magnetic resonance parameters, and gradually moves towards the discussion of the most popular common and newly developed methodologies for quantum chemical modeling of NMR spectra.

A new double-reference correction scheme for accurate and efficient computation of NMR chemical shieldings

Our novel correction procedure yields high-accuracy DFT predictions of absolute NMR shieldings and enables outliers due to relativistic effects or manifestly inadequate modelling of electron correlation to be easily and unambiguously identified.

Hyperfine Coupling Constants in Local Exact Two-Component Theory

We present a highly efficient implementation of the electron-nucleus hyperfine coupling matrix within the one-electron exact two-component (X2C) theory. The complete derivative of the X2C Hamiltonian is formed, that is, the derivatives of the unitary decoupling transformation are considered. This requires the solution of the response and Sylvester equations, consequently increasing the computational costs. Therefore, we apply the diagonal local approximation to the unitary decoupling transformation (DLU). The finite nucleus model is employed for both the scalar potential and the vector potential. Two-electron picture-change effects are modeled with the (modified) screened nuclear spin-orbit approach. Our implementation is fully integral direct and OpenMP-parallelized. An extensive benchmark study regarding the Hamiltonian, the basis set, and the density functional approximation is carried out for a set of 12-17 transition-metal compounds. The error introduced by DLU is negligible, and the DLU-X2C Hamiltonian accurately reproduces its four-component "fully" relativistic parent results. Functionals with a large amount of Hartree-Fock exchange such as CAM-QTP-02 and ωB97X-D are generally favorable. The pure density functional r²SCAN performs remarkably and even outperforms the common hybrid functionals TPSSh and CAM-B3LYP. Fully uncontracted basis sets or contracted quadruple-ζ bases are required for accurate results. The capability of our implementation is demonstrated for [Pt(C₆Cl₅)₄]^- with more than 4700 primitive basis functions and four rare-earth single-molecule magnets: [La(OAr*)₃]^-, [Lu(NR₂)₃]^-, [Lu(OAr*)₃]^-, and [TbPc₂]^-. Here, the results with the spin-orbit DLU-X2C Hamiltonian are in an excellent agreement with the experimental findings of all Pt, La, Lu, and Tb molecules.

Benchmarking the Accuracy of the Direct Random Phase Approximation and σ-Functionals for NMR Shieldings

A method for computing NMR shieldings with the direct random phase approximation (RPA) and the closely related σ-functionals [Trushin, E.; Thierbach, A.; Görling, A. Toward chemical accuracy at low computational cost: density functional theory with σ-functionals for the correlation energy. J. Chem. Phys. 2021, 154, 014104] is presented, which is based on a finite-difference approach. The accuracy is evaluated in benchmark calculations using high-quality coupled cluster values as a reference. Our results show that the accuracy of the computed NMR shieldings using direct RPA is strongly dependent on the density functional theory reference orbitals and improves with increasing amounts of exact Hartree-Fock exchange in the functional. NMR shieldings computed with direct RPA using a Hartree-Fock reference are significantly more accurate than MP2 shieldings and comparable to CCSD shieldings. Also, the basis set convergence is analyzed and it is shown that at least triple-zeta basis sets are required for reliable results.

Paramagnetic NMR Shielding Tensors and Ring Currents: Efficient Implementation and Application to Heavy Element Compounds

We present an efficient implementation of paramagnetic NMR shielding tensors and shifts in a nonrelativistic and scalar-relativistic density functional theory framework. For the latter, we make use of the scalar exact two-component Hamiltonian in its local approximation, and generally we apply the well established (multipole-accelerated) resolution of the identity approximation and the seminumerical exchange approximation. The perturbed density matrix of a paramagnetic NMR shielding calculation is further used to study the magnetically induced current density and ring currents of open-shell systems as illustrated for [U@Bi₁₂]^3-. [U@Bi₁₂]^3- features delocalized highest occupied molecular orbitals and sustains a net diatropic ring current of ca. 18 nA/T through the Bi₁₂ torus similar to the all-metal aromatic heavy-element cluster [Th@Bi₁₂]^4-.

NMR Indirect Spin-Spin Coupling Constants in a Modern Quasi-Relativistic Density Functional Framework

A quasi-relativistic implementation of NMR indirect spin-spin coupling constants is presented. The exact two-component (X2C) Hamiltonian and its diagonal local approximation to the unitary decoupling transformation (DLU) are utilized together with the (modified) screened nuclear spin-orbit approach. In a restricted kinetic balance, the finite nucleus model is available for both the scalar and vector potentials. The implementation supports density functionals up to the fourth rung of Jacob's ladder, i.e., (range-separated) hybrid and local hybrid functionals based on a seminumerical ansatz. We assess the quality of our quasi-relativistic X2C approach by comparison with "fully" relativistic four-component results for small main-group molecules and alkynyl compounds. The mean absolute error introduced by the DLU scheme is less than 0.05 × 10¹⁹ T J^-2 of the reduced coupling constant for the small main-group molecules and 0.5 Hz for the alkynyl compounds. Thus, the error is significantly smaller than finite nucleus size effects for heavy elements. The basis set convergence and the impact of different density functional approximations are further studied. We propose a simple scheme to develop segmented-contracted relativistic all-electron basis sets for NMR spin-spin couplings. Our implementation allows us to perform calculations of extended molecules with reasonable computational effort, which is illustrated for the ¹J(¹¹⁹Sn, ³¹P) coupling constant of a low-valent tin phosphinidenide complex. The corresponding results are in good agreement with the experimental findings.

Quantum chemical calculations of 77 Se and 125 Te nuclear magnetic resonance spectral parameters and their structural applications

An accurate quantum chemical (QC) modeling of ⁷⁷ Se and ¹²⁵ Te nuclear magnetic resonance (NMR) spectra is deeply involved in the NMR structural assignment for selenium and tellurium compounds that are of utmost importance both in organic and inorganic chemistry nowadays due to their huge application potential in many fields, like biology, medicine, and metallurgy. The main interest of this review is focused on the progress in QC computations of ⁷⁷ Se and ¹²⁵ Te NMR chemical shifts and indirect spin-spin coupling constants involving these nuclei. Different computational methodologies that have been used to simulate the NMR spectra of selenium and tellurium compounds since the middle of the 1990s are discussed with a strong emphasis on their accuracy. A special accent is placed on the calculations resorting to the relativistic methodologies, because taking into account the relativistic effects appreciably influences the precision of NMR calculations of selenium and, especially, tellurium compounds. Stereochemical applications of quantum chemical calculations of ⁷⁷ Se and ¹²⁵ Te NMR parameters are discussed so as to exemplify the importance of integrated approach of experimental and computational NMR techniques.

Delocalized relativistic effects, from the viewpoint of halogen bonding

The ability of organic and inorganic compounds bearing both iodine and astatine atoms to form halogen-bond interactions is theoretically investigated. Upon inclusion of the relativistic spin-orbit interaction, the I-mediated halogen bonds are more affected than the At-mediated ones in many cases. This unusual outcome is disconnected from the behavior of iodine's electrons. The significant decrease of astatine electronegativity with the spin-orbit coupling triggers a redistribution of the electron density, which propagates relativistic effects toward the distant iodine atom. This mechanism can be controlled by introducing suitable substituents and, in particular, strengthened by taking advantage of electron-withdrawing inductive and mesomeric effects. Noticeable relativistic effects can actually be transferred to light atoms properties, e.g., the halogen-bond basicity of bridgehead carbon atoms doubled in propellane derivatives.

The DIRAC code for relativistic molecular calculations

DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.

On the heavy atom on light atom relativistic effect in the NMR shielding constants of phosphine tellurides

The relativistic HALA effect has been shown to depend on the spatial deformation of the lone electron pairs of a heavy atom, as demonstrated for alkyl and alkene phosphine tellurides. It was found that HALA effect on phosphorous nuclear magnetic resonance shielding constant is strongly dependent on the spatial arrangements of light substituents on phosphorus, resulting in the deformation of the lone electron pairs of tellurium.

Relativistic heavy atom effect on the 31 P NMR parameters of phosphine chalcogenides. Part 1. Chemical shifts

Four-component density functional theory calculations of ³¹ P NMR chemical shifts have been performed for the representative series of 56 phosphine chalcogenides in order to investigate an influence of different functional groups on the heavy atom relativistic effect on the NMR chemical shifts of light phosphorous atoms (Heavy Atom on Light Atom effect). The validity of the 4-component density functional theory approach used for the wide-scale calculations of the phosphorous chemical shifts in a wide series of phosphine chalcogenides has been confirmed on a small series of 5 representative compounds with the aid of high-quality coupled cluster singles and doubles calculations taking into account solvent, vibrational, and the relativistic corrections in comparison with the experiment.

Theoretical study of the NMR chemical shift of Xe in supercritical condition

In this work we investigate the level of theory necessary for reproducing the non-linear variation of the ¹²⁹Xe nuclear magnetic resonance (NMR) chemical shift with the density of Xe in supercritical conditions. In detail we study how the ¹²⁹Xe chemical shift depends under supercritical conditions on electron correlation, relativistic and many-body effects. The latter are included using a sequential-QM/MM methodology, in which a classical MD simulation is performed first and the chemical shift is then obtained as an average of quantum calculations of 250 MD snapshots conformations carried out for Xe _n clusters (n = 2 - 8 depending on the density). The analysis of the relativistic effects is made at the level of 4-component Hartree-Fock calculations (4c-HF) and electron correlation effects are considered using second order Møller-Plesset perturbation theory (MP2). To simplify the calculations of the relativistic and electron correlation effects we adopted an additive scheme, where the calculations on the Xe _n clusters are carried out at the non-relativistic Hartree-Fock (HF) level, while electron correlation and relativistic corrections are added for all the pairs of Xe atoms in the clusters. Using this approach we obtain very good agreement with the experimental data, showing that the chemical shift of ¹²⁹Xe in supercritical conditions is very well described by cluster calculations at the HF level, with small contributions from relativistic and electron correlation effects.

On the calculation of second-order magnetic properties using subsystem approaches in a relativistic framework

The implementation of second-order magnetic properties in a frozen density embedding scheme in a four component relativistic framework is outlined and applied to model H₂X–H₂O systems (X = Se, Te, Po).

Normal halogen dependence of 13 C NMR chemical shifts of halogenomethanes revisited at the four-component relativistic level

The 'Normal Halogen Dependence' of ¹³ C NMR chemical shifts in the series of halogenomethanes is revisited at the four-component relativistic level. Calculations of ¹³ C NMR chemical shifts of 70 halogenomethanes have been carried out at the density functional theory (DFT) and MP2 levels with taking into account relativistic effects using the four-component relativistic theory of Dirac-Coulomb within the different computational methods (4RPA, 4OPW91) and hybrid computational schemes (MP2 + 4RPA, MP2 + 4OPW91). The most efficient computational protocols are derived for practical purposes. Relativistic shielding effect reaches as much as several hundreds of ppm for heavy halogenomethanes, and to account for this effect in comparison with experiment at the qualitative level, relativistic Dyall's basis sets of triple-zeta quality or higher are to be used within the framework of the four-component relativistic theory taking into account solvent effects. Relativistic geometrical optimization (as compared with the non-relativistic level) is essential for the molecules containing at least two iodines at one carbon atom. Copyright © 2016 John Wiley & Sons, Ltd.

Relativistic Calculations of Nuclear Magnetic Resonance Parameters

Relativistic effects are important for the accurate evaluation of the observables of nuclear magnetic resonance (NMR) spectroscopy, the nuclear magnetic shielding and the indirect spin–spin coupling tensors. Some of the most notable relativistic effects, in particular for light elements in the vicinity of heavy nuclei, are due to spin–orbit effects, an effect difficult to evaluate when starting from a non-relativistic wavefunction. Two- and four-component relativistic methods include spin–orbit effects variationally, and the recent improvements in the computational efficiency of these methods open new opportunities for accurate calculations of NMR parameters also for molecules with heavy elements. We here present an overview of the different approximations that have been introduced for calculating relativistic effects with two- and four-component methods and how these methods can be used to calculate the NMR parameters. We will also give some examples of systems that have been studied computationally with two- and four-component relativistic methods and discuss the importance of relativistic effects on the shielding and indirect spin–spin coupling constants.

Toward an absolute NMR shielding scale using the spin-rotation tensor within a relativistic framework

How can one extend Flygare's rule to the relativistic framework? Three models are proposed here. The best of them shows that σ is related with the spin-rotation tensor, the atomic shielding and a new term coined as σ^SO-S.

Scalar and Spin-Orbit Relativistic Corrections to the NICS and the Induced Magnetic Field: The case of the E12(2-) Spherenes (E = Ge, Sn, Pb)

Can relativistic effects modify the NICS and the B(ind) values? In this manuscript we evaluate the relativistic corrections incorporated via the zeroth-order regular approximation to the calculations of nucleus-independent chemical shifts and the induced magnetic field (B(ind)) in the E12(2-) spherenes (E = Ge, Sn, Pb). We found that both electron delocalization descriptors are strongly affected by the relativistic corrections. For instance, for plumbaspherene, the difference in values from the nonrelativistic to the relativity-included calculation is almost 40 ppm! Our results show that the changes observed in the NICS and B(ind) values in the title cages are a consequence of the treatment of the relativistic effects. If these effects are included as scalar or spin-orbit calculations, then we can establish the next trend: Ge12(2-) is a nonaromatic species, Sn12(2-) is a low aromatic species, and Pb12(2-) is strongly aromatic, according to calculated NICS and B(ind) values. Thus, any prediction of electron delocalization in molecules containing heavy elements without considering an adequate treatment for relativistic effects may lead to an erroneous chemical interpretation.

Four-component relativistic DFT calculations of (77) Se NMR chemical shifts: A gateway to a reliable computational scheme for the medium-sized organoselenium molecules

A versatile high-accuracy computational scheme for the (77) Se nuclear magnetic resonance (NMR) chemical shifts of the medium-sized organoselenium compounds is suggested within a framework of a full four-component relativistic density functional theory (DFT). The main accuracy factors (DFT functionals, relativistic geometry, vibrational corrections, and solvent effects) are addressed. The best result is achieved with NMR-oriented KT2 functional of Keal-Tozer characterized by a fairly small error of only 30 ppm for the span of about 1700 ppm (<2%).

MP2 calculation of (77) Se NMR chemical shifts taking into account relativistic corrections

The main factors affecting the accuracy and computational cost of the Second-order Möller-Plesset perturbation theory (MP2) calculation of (77) Se NMR chemical shifts (methods and basis sets, relativistic corrections, and solvent effects) are addressed with a special emphasis on relativistic effects. For the latter, paramagnetic contribution (390-466 ppm) dominates over diamagnetic term (192-198 ppm) resulting in a total shielding relativistic correction of about 230-260 ppm (some 15% of the total values of selenium absolute shielding constants). Diamagnetic term is practically constant, while paramagnetic contribution spans over 70-80 ppm. In the (77) Se NMR chemical shifts scale, relativistic corrections are about 20-30 ppm (some 5% of the total values of selenium chemical shifts). Solvent effects evaluated within the polarizable continuum solvation model are of the same order of magnitude as relativistic corrections (about 5%). For the practical calculations of (77) Se NMR chemical shifts of the medium-sized organoselenium compounds, the most efficient computational protocols employing relativistic Dyall's basis sets and taking into account relativistic and solvent corrections are suggested. The best result is characterized by a mean absolute error of 17 ppm for the span of (77) Se NMR chemical shifts reaching 2500 ppm resulting in a mean absolute percentage error of 0.7%.

Relativistic effects in the one-bond spin-spin coupling constants involving selenium

One-bond spin-spin coupling constants involving selenium of seven different types, (1)  J(Se,X), X = (1) H, (13) C, (15)  N, (19)  F, (29) Si, (31) P, and (77) Se, were calculated in the series of 14 representative compounds at the SOPPA(CCSD) level taking into account relativistic corrections evaluated both at the RPA and DFT levels of theory in comparison with experiment. Relativistic corrections were found to play a major role in the calculation of (1)  J(Se,X) reaching as much as almost 170% of the total value of (1)  J(Se,Se) and up to 60-70% for the rest of (1)  J(Se,X). Scalar relativistic effects (Darwin and mass-velocity corrections) by far dominate over spin-orbit coupling in the total relativistic effects for all (1)  J(Se,X). Taking into account relativistic corrections at both random phase approximation and density functional theory levels essentially improves the agreement of theoretical results with experiment. The most 'relativistic' (1)  J(Se,Se) demonstrates a marked Karplus-type dihedral angle dependence with respect to the mutual orientation of the selenium lone pairs providing a powerful tool for stereochemical analysis of selenoorganic compounds.

What factors influence the metal-proton spin-spin coupling constants in mercury- and cadmium-substutited rubredoxin?

The indirect metal-proton spin-spin coupling constants between protons in cysteine groups and the mercury or cadmium nucleus have been calculated for a small model of Me-rubredoxin complex (Me = Cd, Hg) by means of density functional theory with zeroth-order regular approximation Hamiltonian (DFT-ZORA). The calculated spin-spin coupling constants, in spite of the moderate size of the model system, are in good agreement with the values measured in NMR experiment, which are in the 0.29-0.56 Hz range for the Cd complex and in the 0.57-2.20 Hz range for the Hg complex. The robustness of the chosen method has been verified by calculations with a number of different exchange-correlation functionals and basis sets. Additionally, it has been shown that the short- and long-distance metal-proton coupling constants are affected mainly by the values of the metal-proton distance and the H-N-C-C dihedral angle.

The Absolute Shielding Constants of Heavy Nuclei: Resolving the Enigma of the (119)Sn Absolute Shielding

We demonstrate that the apparent disagreement between experimental determinations and four-component relativistic calculations of the absolute shielding constants of heavy nuclei is due to the breakdown of the commonly assumed relation between the electronic contribution to the nuclear spin-rotation constants and the paramagnetic contribution to the NMR shielding constants. We demonstrate that this breakdown has significant consequences for the absolute shielding constant of (119)Sn, leading to errors of about 1000 ppm. As a consequence, we expect that many absolute shielding constants of heavy nuclei will be in need of revision.