Juan I. Melo

Relativistic effects on the nuclear magnetic shielding tensor

A new approach for calculating relativistic corrections to the nuclear magnetic shieldings is presented. Starting from a full relativistic second order perturbation theory expression a two-component formalism is constructed by transforming matrix elements using the elimination of small component scheme and separating out the contributions from the no-virtual pair and the virtual pair part of the second order corrections to the energy. In this way we avoid a strong simplification used previously in the literature. We arrive at final expressions for the relativistic corrections which are equivalent to those of Fukui et al. [J. Chem Phys. 105, 3175 (1996)] and at some other additional terms correcting both the paramagnetic and the diamagnetic part of the nuclear magnetic shielding. Results for some relativistic corrections to the shieldings of the heavy and light nuclei in HX and CH3X (X=Br,I) at both random phase and second order polarization propagator approach levels are given.

Relativistic effects on nuclear magnetic shielding constants in HX and CH3X(X=Br, I) based on the linear response within the elimination of small component approach

Numerical calculations of relativistic effects on nuclear magnetic shielding constants sigma corresponding to all one-body operators obtained within a formalism developed in previous work were carried out. In this formalism, the elimination of small component scheme is applied to evaluate all quantities entering a four-component RSPT(2) expression of magnetic molecular properties. HX and CH3X (X=Br,I) were taken as model compounds. Calculations were carried out at the Hartree-Fock level for first-order quantities, and at the random-phase approximation (RPA) level for second- and third-order ones. It was found that values of sigma(X) are largely affected by several relativistic corrections not previously considered in the bibliography. sigma Values of the H nucleus are in close agreement with four-component RPA ones. Overall relativistic effects on the shift of sigma(X) from HX to CH3X are smaller than the nonrelativistic shifts.

Nuclear magnetic resonance shielding constants and chemical shifts in linear 199Hg compounds: A comparison of three relativistic computational methods

We investigate the importance of relativistic effects on NMR shielding constants and chemical shifts of linear HgL(2) (L = Cl, Br, I, CH(3)) compounds using three different relativistic methods: the fully relativistic four-component approach and the two-component approximations, linear response elimination of small component (LR-ESC) and zeroth-order regular approximation (ZORA). LR-ESC reproduces successfully the four-component results for the C shielding constant in Hg(CH(3))(2) within 6 ppm, but fails to reproduce the Hg shielding constants and chemical shifts. The latter is mainly due to an underestimation of the change in spin-orbit contribution. Even though ZORA underestimates the absolute Hg NMR shielding constants by ∼2100 ppm, the differences between Hg chemical shift values obtained using ZORA and the four-component approach without spin-density contribution to the exchange-correlation (XC) kernel are less than 60 ppm for all compounds using three different functionals, BP86, B3LYP, and PBE0. However, larger deviations (up to 366 ppm) occur for Hg chemical shifts in HgBr(2) and HgI(2) when ZORA results are compared with four-component calculations with non-collinear spin-density contribution to the XC kernel. For the ZORA calculations it is necessary to use large basis sets (QZ4P) and the TZ2P basis set may give errors of ∼500 ppm for the Hg chemical shifts, despite deceivingly good agreement with experimental data. A Gaussian nucleus model for the Coulomb potential reduces the Hg shielding constants by ∼100-500 ppm and the Hg chemical shifts by 1-143 ppm compared to the point nucleus model depending on the atomic number Z of the coordinating atom and the level of theory. The effect on the shielding constants of the lighter nuclei (C, Cl, Br, I) is, however, negligible.

Relativistic calculation of indirect NMR spin-spin couplings using the Douglas-Kroll-Hess approximation

We have employed the Douglas-Kroll-Hess approximation to derive the perturbative Hamiltonians involved in the calculation of NMR spin-spin couplings in molecules containing heavy elements. We have applied this two-component quasirelativistic approach using finite perturbation theory in combination with a generalized Kohn-Sham code that includes the spin-orbit interaction self-consistently and works with Hartree-Fock and both pure and hybrid density functionals. We present numerical results for one-bond spin-spin couplings in the series of tetrahydrides CH(4), SiH(4), GeH(4), and SnH(4). Our two-component Hartree-Fock results are in good agreement with four-component Dirac-Hartree-Fock calculations, although a density-functional treatment better reproduces the available experimental data.

Magnetic Exchange Couplings with Range-Separated Hybrid Density Functionals

We investigate the effect of Hartree-Fock range-separation on the calculation of magnetic exchange couplings in a set of nine bimetallic transition-metal complexes containing 3d elements (V, Cr, Mn, and Cu). To this end, we have compared magnetic exchange couplings calculated as self-consistent energy differences using two global hybrid functionals, B3LYP (Becke 3-parameter exchange and Lee-Yang-Parr correlation) and PBEh (hybrid Perdew-Burke-Ernzerhof) with the short-range separated HSE (Heyd-Scuseria-Ernzerhof) and the long-range corrected LC-ωPBE. Our results show that, although there is no clear superiority of any of these functionals when compared with experimental data, the LC-ωPBE provides a better description of the magnetization on the metallic centers, yielding self-consistent solutions that mimic more closely a Heisenberg-like behavior.

Relativistic corrections to the diamagnetic term of the nuclear magnetic shielding: Analysis of contributions from localized orbitals

We have calculated the relativistic corrections to the diamagnetic term of the nuclear magnetic shielding constants for a series of molecules containing heavy atoms. An analysis of the contributions from localized orbitals is performed. We establish quantitatively the relative importance of inner core and valence shell molecular orbitals in each correcting term. Contributions from the latter are much less important than those from the former. The calculated values of the correction sigma(L-PSO), first derived within the linear response elimination of small component formalism, show a power-law dependence on the nuclear charge approximately Z(3.5), in contrast with the approximately Z(3.1) behavior of the mass-velocity external-field correction to the paramagnetic term previously reported.

Formal relations connecting different approaches to calculate relativistic effects on molecular magnetic properties

In the present work a set of formal relations connecting different approaches to calculate relativistic effects on magnetic molecular properties are proven. The linear response (LR) within the elimination of the small component (ESC), Breit Pauli, and minimal-coupling approaches are compared. To this end, the leading order ESC reduction of operators within the minimal-coupling four-component approach is carried out. The equivalence of all three approaches within the ESC approximation is proven. It is numerically verified for the NMR nuclear-magnetic shielding tensor taking HX and CH3X (X=Br,I) as model compounds. Formal relations proving the gauge origin invariance of the full relativistic effect on the NMR nuclear-magnetic shielding tensor within the LR-ESC approach are presented.

Core-dependent and ligand-dependent relativistic corrections to the nuclear magnetic shieldings in MH4-n Y n (n = 0-4; M = Si, Ge, Sn, and Y = H, F, Cl, Br, I) model compounds.

The nuclear magnetic shieldings of Si, Ge, and Sn in MH4−nYn (M = Si, Ge, Sn; Y = F, Cl, Br, I and n = 1–4) molecular systems are highly influenced by the substitution of one or more hydrogens by heavy-halogen atoms. We applied the linear response elimination of small components (LRESC) formalism to calculate those shieldings and learn whether including only a few of the leading relativistic correction terms is sufficient to be able to quantitatively reproduce the full relativistic value. It was observed that the nuclear magnetic shieldings change as the number of heavy halogen substituents and their weights vary, and the pattern of σ(M) generally does not exhibit the normal halogen dependence (NHD) behavior that can be seen in similar molecular systems containing carbon atoms. We also analyzed each relativistic correction afforded by the LRESC method and split them in two: core-dependent and ligand-dependent contributions; we then looked for the electronic mechanisms involved in the different relativistic effects and in the total relativistic value. Based on this analysis, we were able to study the electronic mechanism involved in a recently proposed relativistic effect, the “heavy atom effect on vicinal heavy atom” (HAVHA), in more detail. We found that the main electronic mechanism is the spin–orbit or σpT(3) correction, although other corrections such as σpS(1) and σpS(3) are also important. Finally, we analyzed proton magnetic shieldings and found that, for molecules containing Sn as the central atom, σ(H) decreases as the number of heavy halogen substituents (of the same type: either F, Cl, or Br) increases, albeit at different rates for different halogens. σ(H) only increase as the number of halogen substituents increases if the halogen is iodine.

Relativistic two-component geometric approximation of the electron-positron contribution to magnetic properties in terms of Breit–Pauli spinors

An alternative approach for the calculation of the electron-positron (e-p) contribution to magnetic properties based on two-component Breit-Pauli spinors is presented. In it, the elimination of the small component scheme is applied to the inverse propagator matrix of e-p pairs. The effect of the positronic manifold is expressed as an operator acting on Breit-Pauli spinors. The operator form thus obtained sums up the relativistic correction as a geometric series and as a result a totally different behavior in the vicinity of a nucleus is obtained as compared to the one of the linear response approximation. This feature has deep influence in numerical values of the e-p contribution to the nuclear magnetic shielding of heavy atoms. Numerical calculations carried out for Kr, Xe, and I show that with this approach, the e-p contributions to this property are in good agreement with those of four-component methods.

NMR nuclear magnetic shielding anisotropy of linear molecules within the linear response within the elimination of the small component approach.

The influence of the spin-Zeeman (SZ) operator in the evaluation of the spin-orbit effect on the nuclear magnetic shielding tensor in the context of the linear response within the elimination of the small component approach is critically discussed. It is shown that such term yields no contribution to the isotropic nuclear magnetic shielding constant, but it may be of great importance in the determination of individual tensor components, and particularly of the tensor anisotropy. In particular, an interesting relation between the SZ and orbital Zeeman contributions to the spin-orbit effect for the case of linear molecules is shown to hold. Numerical examples for the BrH, IH, and XeF(2) molecules are presented which show that, provided the SZ term is taken into account, results of the individual shielding tensor components and the tensor anisotropy are in good agreement with those obtained by other theoretical methods, and particularly by the Dirac-Hartree-Fock approach.