Orbital contributions to relativistic corrections of the NMR nuclear magnetic shielding tensor originated in scalar field-dependent operators

The LRESC-Loc Model to Analyze Magnetic Shieldings with Localized Molecular Orbitals

The leading electronic mechanisms of relativistic effects in the NMR magnetic shieldings of heavy-atom (HA) containing molecules are well described by the linear response with elimination of small components model (LRESC). We show here first results from a new version of the LRESC model written in terms of localized molecular orbitals (LMOs) which is coined as LRESC-Loc. Those LMOs resemble "chemist's orbitals", representing lone-pairs, atomic cores, and bonds. The whole set of relativistic effects are expressed in terms of non-ligand-dependent and ligand-dependent contributions. We show the electronic origin of trends and behavior of different mechanisms in molecular systems which contain heavy elements that belong to any of the IB to VIIA groups of the periodic table. The SO mechanism has a well-defined dependence with the LPs (LPσ and LPπ) when the HAs have them, but the non-SO mechanisms mostly depend on other LMOs. In addition we propose here that the SO mechanism can be used to characterize interactions involving LPs and the non-SO mechanisms to characterize covalent and close-shell interactions. All our main results are in accord with previous findings, though we are now able to analyze them in a different manner.

Performance of the LRESC Model on top of DFT Functionals for Relativistic NMR Shielding Calculations

The linear response within the elimination of the small component model (LRESC) is an insightful and computationally efficient method for including relativistic effects on molecular properties like the nuclear magnetic shielding constants, spin-rotation constant, g-tensor, and electric field gradient of heavy atom containing molecules with atoms belonging up to the sixth row of the periodic table. One of its main advantages is its capacity to analyze the electronic origin of the different relativistic correcting terms. Until now, it was always applied on top of Hartree-Fock ground-state wave functions (LRESC/HF) to calculate and analyze NMR shieldings. In this work, we show the performance of the LRESC formalism on top of some density functional theory (DFT) functionals to compute tin shielding constants in SnX₄ (X = H, F, Cl, Br, I) molecular systems. We analyze the performance of each LRESC/DFT scheme on reproducing the electronic mechanisms of the shieldings, taking as a benchmark the results of relativistic calculations at the RPA level of approach (4c/RPA). As in previous works, we divide the LRESC relativistic correcting terms into two groups: core-dependent and ligand-dependent contributions. It is shown here that core-dependent corrections are well-reproduced for the selected DFT functionals, but some differences arise in the ligand-dependent ones. We focus on the performance of different functionals, including the same electron correlation part but containing different amounts of HF exchange. The best results are obtained for the BHandHLYP functional (50% of HF exchange) and the worst for BLYP (0%). When the percentage of HF exchange increases, ligand-dependent contributions are better described, and the final LRESC/DFT results are closer to those obtained with LRESC/HF and 4c/RPA methods. The spin-orbit correction to the shielding constant is one of the main ligand-dependent contributions (there are two more) with total value depending on the amount of HF exchange included in the functional. When the amount of HF exchange decreases, the spin-orbit contribution becomes larger, overestimating the shielding constant even when nonrelativisitc DFT values are much smaller than the nonrelativistic HF ones, as it happens for the heaviest molecular system studied here (SnI₄).

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.

Theory and computation of nuclear magnetic resonance parameters

The art of quantum chemical electronic structure calculation has over the last 15 years reached a point where systematic computational studies of magnetic response properties have become a routine procedure for molecular systems. One of their most prominent areas of application are the spectral parameters of nuclear magnetic resonance (NMR) spectroscopy, due to the immense importance of this experimental method in many scientific disciplines. This article attempts to give an overview on the theory and state-of-the-art of the practical computations in the field, in terms of the size of systems that can be treated, the accuracy that can be expected, and the various factors that would influence the agreement of even the most accurate imaginable electronic structure calculation with experiment. These factors include relativistic effects, thermal effects, as well as solvation/environmental influences, where my group has been active. The dependence of the NMR spectra on external magnetic and optical fields is also briefly touched on.

Relativistic heavy-atom effects on heavy-atom nuclear shieldings

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

Mean field linear response within the elimination of the small component formalism to evaluate relativistic effects on magnetic properties

The linear response within the elimination of the small component formalism is aimed at obtaining the leading order relativistic corrections to magnetic molecular properties in the context of the elimination of the small component approximation. In the present work we extend the method in order to include two-body effects in the form of a mean field one-body operator. To this end we consider the four-component Dirac-Hartree-Fock operator as the starting point in the evaluation of the second order relativistic expression of magnetic properties. The approach thus obtained is the fully consistent leading order approximation of the random phase approximation four-component formalism. The mean field effect on the relativistic corrections to both the diamagnetic and paramagnetic terms of magnetic properties taking into account both the Coulomb and Breit two-body interactions is considered.

Leading-order relativistic effects on nuclear magnetic resonance shielding tensors

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

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.

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.