Fully relativistic calculation of nuclear magnetic shieldings and indirect nuclear spin-spin couplings in group-15 and -16 hydrides

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 $$ {}_n $$ ]2 - 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 $$ {}_n $$ ]2 - 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∼ $$ \sim $$ 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.

New pecJ-n (n = 1, 2) Basis Sets for Selenium Atom Purposed for the Calculations of NMR Spin-Spin Coupling Constants Involving Selenium

We present new compact pecJ-n (n = 1, 2) basis sets for the selenium atom developed for the quantum-chemical calculations of NMR spin-spin coupling constants (SSCCs) involving selenium nuclei. These basis sets were obtained at the second order polarization propagator approximation with coupled cluster singles and doubles amplitudes (SOPPA(CCSD)) level with the property-energy consistent (PEC) method, which was introduced in our previous papers. The existing SSCC-oriented selenium basis sets are rather large in size, while the PEC method gives more compact basis sets that are capable of providing accuracy comparable to that reached using the property-oriented basis sets of larger sizes generated with a standard even-tempered technique. This is due to the fact that the PEC method is very different in its essence from the even-tempered approaches. It generates new exponents through the total optimization of angular spaces of trial basis sets with respect to the property under consideration and the total molecular energy. New basis sets were tested on the coupled cluster singles and doubles (CCSD) calculations of SSCCs involving selenium in the representative series of molecules, taking into account relativistic, solvent, and vibrational corrections. The comparison with the experiment showed that the accuracy of the results obtained with the pecJ-2 basis set is almost the same as that provided by a significantly larger basis set, aug-cc-pVTZ-J, while that achieved with a very compact pecJ-1 basis set is only slightly inferior to the accuracy provided by the former.

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.

Vibrational Corrections to NMR Spin-Spin Coupling Constants from Relativistic Four-Component DFT Calculations

Zero-point vibrational (ZPV) corrections to the nuclear spin–spin coupling constants have been calculated using four-component Dirac–Kohn–Sham DFT for H₂X (where X = O, S, Se, Te, Po), XH₃ (where X = N, P, As, Sb, Bi), and XH₄ (where X = C, Si, Ge, Sn, and Pb) molecules and for HC≡CPbH₃. The main goal was to study the influence of relativistic effects on the ZPV corrections and thus results calculated at relativistic and nonrelativistic approaches have been compared. The effects of relativity become notable for the ZPV corrections to the spin–spin coupling constants for compounds with lighter elements (selenium and germanium) than for the spin–spin coupling constants themselves. In the case of molecules containing heavier atoms, for instance BiH₃ and PbH₄, relativistic effects play a crucial role on the results and approximating ZPV corrections by the nonrelativistic results may lead to larger errors than omitting ZPV corrections altogether.

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.

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.

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.

Calculation of 15N NMR chemical shifts: Recent advances and perspectives

Recent advances in computation of ¹⁵N NMR chemical shifts are reviewed, concentrating mainly on practical aspects of computational protocols and accuracy factors. The review includes the discussion of the level of theory, the choice of density functionals and basis sets together with taking into account solvent effects, rovibrational corrections and relativistic effects. Computational aspects of ¹⁵N NMR are illustrated for the series of neutral and protonated open-chain nitrogen-containing compounds and nitrogen heterocycles, coordination and intermolecular complexes.

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 and electron-correlation effects on the nuclear magnetic resonance shieldings of molecules containing tin and lead atoms

The reference values for NMR magnetic shieldings, σ(ref), are of the highest importance when theoretical analysis of chemical shifts are envisaged. The fact that the nonrelativistically valid relationship among spin-rotation constants and magnetic shieldings is not any longer valid for heavy atoms requires that the search for σ(ref) for such atoms needs new strategies to follow. We present here results of σ(ref) that were obtained by applying our own simple procedure which mixes accurate experimental chemical shifts (δ) and theoretical magnetic shieldings (σ). We calculated σ(Sn) and σ(Pb) in a family of heavy-halogen-containing molecules. We found out that σ(ref)[Sn;Sn(CH3)4] in gas phase should be close to 3864.11 ± 20.05 ppm (0.5%). For Pb atom, σ(ref)[Pb;Pb(CH3)4] should be close to 14475.1 ± 500.7 ppm. Such theoretical values correspond to calculations with the relativistic polarization propagator method, RelPPA, at the RPA level of approach. They are closer to experimental values as compared to those obtained applying few different functionals such as PBE0, B3LYP, BLYP, BP86, KT2, and KT3 of the density functional theory, DFT. We studied tin and lead shieldings of the XY(4-n)Z(n) (X = Sn, Pb; Y, Z = H, F, Cl, Br, I) and PbH(4-n)I(n) (n = 0, 1, 2, 3, 4) family of compounds with four-component functionals as implemented in the DIRAC code. For these systems results of calculations with RelPPA-RPA are more reliable than DFT ones. We argue about why those DFT functionals must be modified in order to obtain more accurate results of NMR magnetic shieldings within the relativistic regime: first, there is a dependence among both electron-correlation and relativistic effects that should be introduced in some way in the functionals; and second, the DIRAC code uses standard nonrelativistic functionals and the functionals B3LYP and PBE0 were parametrized only with data taken from light elements. It can explain why they are not able to properly introduce relativistic effects on nuclear magnetic shieldings. We finally show that in the analysis of magnetic shieldings for the family of compounds mentioned above, one must consider the newest and so-called heavy-atom effect on vicinal heavy atoms, HAVHA. Such effects are among the most important relativistic effects in these kind of compounds.

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.

Conformational analysis and diastereotopic assignments in the series of selenium-containing heterocycles by means of 77Se-1H spin-spin coupling constants: a combined theoretical and experimental study

A combined theoretical and experimental study on the stereochemical behavior of (77)Se-(1)H spin-spin coupling constants has been performed at the second-order polarization propagator approach level together with heteronuclear multiple-bond correlation technique in the series of selenium-containing four-, five- and six-membered heterocycles including the derivatives of thiaselenetane, selenasilole, thiaselenole, thiaselenolane and dihydrothiaselenine. Geminal and vicinal (77)Se-(1)H spin-spin couplings were shown to have the pronounced stereochemical dependences in respect with the topology of the coupling pathway, internal rotation of the side-chain substituents and ring inversion providing a straightforward tool for the conformational analysis and diastereotopic assignments in the chiral organoselenium compounds.

Structural trends of 77Se-1H spin-spin coupling constants and conformational behavior of 2-substituted selenophenes

Experimental measurements and second-order polarization propagator approach (SOPPA) calculations of (77)Se-(1)H spin-spin coupling constants together with theoretical energy-based conformational analysis in the series of 2-substituted selenophenes have been carried out. A new basis set optimized for the calculation of (77)Se-(1)H spin-spin coupling constants has been introduced by extending the aug-cc-pVTZ-J basis for selenium. Most of the spin-spin coupling constants under study, especially vicinal (77)Se-(1)H couplings, demonstrated a remarkable stereochemical behavior with respect to the internal rotation of the substituent in the 2-position of the selenophene ring, which is of major importance in the stereochemical studies of the related organoselenium compounds.

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.

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.

Quasirelativistic theory for the magnetic shielding constant. III. Quasirelativistic second-order Møller–Plesset perturbation theory and its application to tellurium compounds

The quasirelativistic (QR) generalized unrestricted Hartree-Fock method for the magnetic shielding constant [R. Fukuda, M. Hada, and H. Nakatsuji, J. Chem. Phys. 118, 1015 (2003); R. Fukuda, M. Hada, and H. Nakatsuji, J. Chem. Phys.118, 1027 (2003)] has been extended to include the electron correlation effect in the level of the second-order Møller-Plesset perturbation theory (MP2). We have implemented the energy gradient and finite-perturbation methods to calculate the magnetic shielding constant at the QR MP2 level and applied to the magnetic shielding constants and the NMR chemical shifts of 125Te nucleus in various tellurium compounds. The calculated magnetic shielding constants and NMR chemical shifts well reproduced the experimental values. The relations of the chemical shifts with the natures of ligands, and the tellurium oxidation states were investigated. The chemical shifts in different valence states were explained by the paramagnetic shielding and spin-orbit terms. The tellurium 5p electrons are the dominant origin of the chemical shifts in the Te I and Te II compounds and the chemical shifts were explained by the p-hole mechanism. The tellurium d electrons also play an important role in the chemical shifts of the hypervalent compounds.

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.