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₄).

Big picture of relativistic molecular quantum mechanics

Any quantum mechanical calculation on electronic structure ought to choose first an appropriate Hamiltonian H and then an Ansatz for parameterizing the wave function Ψ, from which the desired energy/property E(λ) can finally be calculated. Therefore, the very first question is: what is the most accurate many-electron Hamiltonian H? It is shown that such a Hamiltonian i.e. effective quantum electrodynamics (eQED) Hamiltonian, can be obtained naturally by incorporating properly the charge conjugation symmetry when normal ordering the second quantized fermion operators. Taking this eQED Hamiltonian as the basis, various approximate relativistic many-electron Hamiltonians can be obtained based entirely on physical arguments. All these Hamiltonians together form a complete and continuous ‘Hamiltonian ladder’, from which one can pick up the right one according to the target physics and accuracy. As for the many-electron wave function Ψ, the most intriguing questions are as follows. (i) How to do relativistic explicit correlation? (ii) How to handle strong correlation? Both general principles and practical strategies are outlined here to handle these issues. Among the electronic properties E(λ) that sample the electronic wave function nearby the nuclear region, nuclear magnetic resonance (NMR) shielding and nuclear spin-rotation (NSR) coupling constant are especially challenging: they require body-fixed molecular Hamiltonians that treat both the electrons and nuclei as relativistic quantum particles. Nevertheless, they have been formulated rigorously. In particular, a very robust ‘relativistic mapping’ between the two properties has been established, which can translate experimentally measured NSR coupling constants to very accurate absolute NMR shielding scales that otherwise cannot be obtained experimentally. Since the most general and fundamental issues pertinent to all the three components of the quantum mechanical equation HΨ = EΨ (i.e. Hamiltonian H, wave function Ψ, and energy/property E(λ)) have fully been understood, the big picture of relativistic molecular quantum mechanics can now be regarded as established.

Relativistic calculations of magnetic resonance parameters: background and some recent developments

This article outlines some basic concepts of relativistic quantum chemistry and recent developments of relativistic methods for the calculation of the molecular properties that define the basic parameters of magnetic resonance spectroscopic techniques, i.e. nuclear magnetic resonance shielding, indirect nuclear spin-spin coupling and electric field gradients (nuclear quadrupole coupling), as well as with electron paramagnetic resonance g-factors and electron-nucleus hyperfine coupling. Density functional theory (DFT) has been very successful in molecular property calculations, despite a number of problems related to approximations in the functionals. In particular, for heavy-element systems, the large electron count and the need for a relativistic treatment often render the application of correlated wave function ab initio methods impracticable. Selected applications of DFT in relativistic calculation of magnetic resonance parameters are reviewed.