Connection between the regular approximation and the normalized elimination of the small component in relativistic quantum theory

Generalized Foldy-Wouthuysen transformation for relativistic two-component methods: Systematic analysis of two-component Hamiltonians

The generalized Foldy-Wouthuysen (GFW) transformation was proposed as a generic form that unifies four types of transformations in relativistic two-component methods: unnormalized GFW(UN), and normalized form 1, form 2, and form 3 (GFW(N1), GFW(N2), and GFW(N3)). The GFW transformation covers a wide range of transformations beyond the simple unitary transformation of the Dirac Hamiltonian, allowing for the systematic classification of all existing two-component methods. New two-component methods were also systematically derived based on the GFW transformation. These various two-component methods were applied to hydrogen-like and helium-like ions. Numerical errors in energy were evaluated and classified into four types: the one-electron Hamiltonian approximation, the two-electron operator approximation, the newly defined "picture difference error (PDE)," and the error in determining the transformation, and errors in multi-electron systems were discussed based on this classification.

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.

Douglas-Kroll and infinite order two-component transformations of Dirac-Fock operator

We extended the conventional Douglas-Kroll (DK) and infinite order two-component (IOTC) methods to a technique applicable to Fock matrices, called extended DK (EDK) and extended IOTC (EIOTC), respectively. First, we defined a strategy to divide the Dirac-Fock operator into zero- and first-order terms. We then demonstrated that the first-order extended DK transformation, which is the Foldy-Wouthuysen transformation for the zero-order term, as well as the second- and third-order EDK and EIOTC, could be well defined. The EDK- and EIOTC-transformed Fock matrix, kinetic energy operator, nuclear attraction operator, and density matrix were derived. These equations were numerically evaluated, and it was found that these methods were accurate. In particular, EIOTC was consistent with the four-component approach. Four-component and extended two-component calculations are more expensive than non-relativistic calculations due to small-component-type two-electron integrals. We developed a new approximation formula, RIS-V, for small-component-type two-electron integrals, including the spin-orbit interaction between electrons. These results suggest that the RIS-V formula effectively accelerates the four-component and extended two-component methods.

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.

Exact two-component Hamiltonians for relativistic quantum chemistry: Two-electron picture-change corrections made simple

Based on self-consistent field (SCF) atomic mean-field (amf) quantities, we present two simple yet computationally efficient and numerically accurate matrix-algebraic approaches to correct both scalar-relativistic and spin-orbit two-electron picture-change effects (PCEs) arising within an exact two-component (X2C) Hamiltonian framework. Both approaches, dubbed amfX2C and e(xtended)amfX2C, allow us to uniquely tailor PCE corrections to mean-field models, viz. Hartree-Fock or Kohn-Sham DFT, in the latter case also avoiding the need for a point-wise calculation of exchange-correlation PCE corrections. We assess the numerical performance of these PCE correction models on spinor energies of group 18 (closed-shell) and group 16 (open-shell) diatomic molecules, achieving a consistent ≈10^-5 Hartree accuracy compared to reference four-component data. Additional tests include SCF calculations of molecular properties such as absolute contact density and contact density shifts in copernicium fluoride compounds (CnF_n, n = 2,4,6), as well as equation-of-motion coupled-cluster calculations of x-ray core-ionization energies of 5d- and 6d-containing molecules, where we observe an excellent agreement with reference data. To conclude, we are confident that our (e)amfX2C PCE correction models constitute a fundamental milestone toward a universal and reliable relativistic two-component quantum-chemical approach, maintaining the accuracy of the parent four-component one at a fraction of its computational cost.

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.

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.

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.

An efficient implementation of two-component relativistic exact-decoupling methods for large molecules

We present an efficient algorithm for one- and two-component relativistic exact-decoupling calculations. Spin-orbit coupling is thus taken into account for the evaluation of relativistically transformed (one-electron) Hamiltonian. As the relativistic decoupling transformation has to be evaluated with primitive functions, the construction of the relativistic one-electron Hamiltonian becomes the bottleneck of the whole calculation for large molecules. For the established exact-decoupling protocols, a minimal matrix operation count is established and discussed in detail. Furthermore, we apply our recently developed local DLU scheme [D. Peng and M. Reiher, J. Chem. Phys. 136, 244108 (2012)] to accelerate this step. With our new implementation two-component relativistic density functional calculations can be performed invoking the resolution-of-identity density-fitting approximation and (Abelian as well as non-Abelian) point group symmetry to accelerate both the exact-decoupling and the two-electron part. The capability of our implementation is illustrated at the example of silver clusters with up to 309 atoms, for which the cohesive energy is calculated and extrapolated to the bulk.

Development and application of the analytical energy gradient for the normalized elimination of the small component method

The analytical energy gradient of the normalized elimination of the small component (NESC) method is derived for the first time and implemented for the routine calculation of NESC geometries and other first order molecular properties. Essential for the derivation is the correct calculation of the transformation matrix U relating the small component to the pseudolarge component of the wavefunction. The exact form of ∂U/∂λ is derived and its contribution to the analytical energy gradient is investigated. The influence of a finite nucleus model and that of the picture change is determined. Different ways of speeding up the calculation of the NESC gradient are tested. It is shown that first order properties can routinely be calculated in combination with Hartree-Fock, density functional theory (DFT), coupled cluster theory, or any electron correlation corrected quantum chemical method, provided the NESC Hamiltonian is determined in an efficient, but nevertheless accurate way. The general applicability of the analytical NESC gradient is demonstrated by benchmark calculations for NESC/CCSD (coupled cluster with all single and double excitation) and NESC/DFT involving up to 800 basis functions.

Bonding in mercury-alkali molecules: Orbital-driven van der Waals complexes

The bonding situation in mercury-alkali diatomics HgA (²Σ⁺) (A = Li, Na, K, Rb) has been investigated employing the relativistic all-electron method Normalized Elimination of the Small Component (NESC), CCSD(T), and augmented VTZ basis sets. Although Hg,A interactions are typical of van der Waals complexes, trends in calculated D_e values can be explained on the basis of a 3-electron 2-orbital model utilizing calculated ionization potentials and the D_e values of HgA⁺(¹Σ⁺) diatomics. HgA molecules are identified as orbital-driven van der Waals complexes. The relevance of results for the understanding of the properties of liquid alkali metal amalgams is discussed.

Decoupling of the Dirac equation correct to the third order for the magnetic perturbation

A two-component relativistic theory accurately decoupling the positive and negative states of the Dirac Hamiltonian that includes magnetic perturbations is derived. The derived theory eliminates all of the odd terms originating from the nuclear attraction potential V and the first-order odd terms originating from the magnetic vector potential A, which connect the positive states to the negative states. The electronic energy obtained by the decoupling is correct to the third order with respect to A due to the (2n+1) rule. The decoupling is exact for the magnetic shielding calculation. However, the calculation of the diamagnetic property requires both the positive and negative states of the unperturbed (A=0) Hamiltonian. The derived theory is applied to the relativistic calculation of nuclear magnetic shielding tensors of HX (X=F,Cl,Br,I) systems at the Hartree-Fock level. The results indicate that such a substantially exact decoupling calculation well reproduces the four-component Dirac-Hartree-Fock results.

Relativistic calculation of nuclear magnetic shielding tensor using the regular approximation to the normalized elimination of the small component. II. Consideration of perturbations in the metric operator

A previous relativistic shielding calculation theory based on the regular approximation to the normalized elimination of the small component approach is improved by the inclusion of the magnetic interaction term contained in the metric operator. In order to consider effects of the metric perturbation, the self-consistent perturbation theory is used for the case of perturbation-dependent overlap integrals. The calculation results show that the second-order regular approximation results obtained for the isotropic shielding constants of halogen nuclei are well improved by the inclusion of the metric perturbation to reproduce the fully relativistic four-component Dirac-Hartree-Fock results. However, it is shown that the metric perturbation hardly or does not affect the anisotropy of the halogen shielding tensors and the proton magnetic shieldings.

Quasirelativistic theory. II. Theory at matrix level

The Dirac operator in a matrix representation in a kinetically balanced basis is transformed to the matrix representation of a quasirelativistic Hamiltonian that has the same electronic eigenstates as the original Dirac matrix (but no positronic eigenstates). This transformation involves a matrix X, for which an exact identity is derived and which can be constructed either in a noniterative way or by various iteration schemes, not requiring an expansion parameter. Both linearly convergent and quadratically convergent iteration schemes are discussed and compared numerically. The authors present three rather different schemes, for each of which even in unfavorable cases convergence is reached within three or four iterations, for all electronic eigenstates of the Dirac operator. The authors present the theory both in terms of a non-Hermitian and a Hermitian quasirelativistic Hamiltonian. Quasirelativistic approaches at the matrix level known from the literature are critically analyzed in the frame of the general theory.

An infinite-order two-component relativistic Hamiltonian by a simple one-step transformation

The authors report the implementation of a simple one-step method for obtaining an infinite-order two-component (IOTC) relativistic Hamiltonian using matrix algebra. They apply the IOTC Hamiltonian to calculations of excitation and ionization energies as well as electric and magnetic properties of the radon atom. The results are compared to corresponding calculations using identical basis sets and based on the four-component Dirac-Coulomb Hamiltonian as well as Douglas-Kroll-Hess and zeroth-order regular approximation Hamiltonians, all implemented in the DIRAC program package, thus allowing a comprehensive comparison of relativistic Hamiltonians within the finite basis approximation.

Relativistic calculation of nuclear magnetic shielding tensor including two-electron spin-orbit interactions

A relativistic calculation of nuclear magnetic shielding tensor including two-electron spin-orbit interactions is performed. In order to reduce the computational load in evaluating the two-electron relativistic integrals, the charge density is approximated by a linear combination of the squares of s-type spatial basis functions. Including the two-electron spin-orbit interaction effect is found to improve the calculation results.

Comment on “Quasirelativistic theory equivalent to fully relativistic theory” [J. Chem. Phys. 123, 241102 (2005)]

The connection between the exact quasirelativistic approach developed in the title reference [W. Kutzelnigg and W. Liu, J. Chem. Phys. 123, 241102 (2005)] and the method of elimination of the small component in matrix form developed previously by Dyall is explicitly worked out. An equation that links Hermitian and non-Hermitian formulations of the exact quasirelativistic theory is derived. Besides establishing a kinship between the existing formulations, the proposed equation can be employed for the derivation of new formulations of the exact quasirelativistic theory.

Infinite-order quasirelativistic density functional method based on the exact matrix quasirelativistic theory

The exact one-electron matrix quasirelativistic theory [Kutzelnigg and Liu, J. Chem. Phys. 123, 241102 (2005)] is extended to the effective one-particle Kohn-Sham scheme of density functional theory. Several variants of the resultant theory are discussed. Although they are in principle equivalent, consideration of computational efficiency strongly favors the one (F(+)) in which the effective potential remains untransformed. Further combined with the atomic approximation for the matrix X relating the small and large components of the Dirac spinors as well as a simple ansatz for correcting the two-electron picture change errors, a very elegant, accurate, and efficient infinite-order quasirelativistic approach is obtained, which is far simpler than all existing quasirelativistic theories and must hence be regarded as a breakthrough in relativistic quantum chemistry. In passing, it is also shown that the Dirac-Kohn-Sham scheme can be made as efficient as two-component approaches without compromising the accuracy. To demonstrate the performance of the new methods, atomic calculations on Hg and E117 are first carried out. The spectroscopic constants (bond length, vibrational frequency, and dissociation energy) of E117(2) are then reported. All the results are in excellent agreement with those of the Dirac-Kohn-Sham calculations.

Relativistic calculation of nuclear magnetic shielding using normalized elimination of the small component

The normalized elimination of the small component (NESC) theory, recently proposed by Filatov and Cremer, is extended to include magnetic interactions and applied to the calculation of the nuclear magnetic shielding in HX (X=F, Cl, Br, I) systems. The NESC calculations are performed at the levels of the zeroth-order regular approximation (ZORA) and the second-order regular approximation (SORA). The calculations show that the NESC-ZORA results are very close to the NESC-SORA results, except for the shielding of the I nucleus. Both the NESC-ZORA and NESC-SORA calculations yield very similar results to the previously reported values obtained using the relativistic infinite-order two-component coupled Hartree-Fock method. The difference between NESC-ZORA and NESC-SORA results is significant for the shieldings of iodine.

Quasirelativistic theory equivalent to fully relativistic theory

The Dirac operator in a matrix representation in a kinetically balanced basis is transformed to a quasirelativistic Hamiltonian matrix, that has the same electronic eigenstates as the original Dirac matrix. This transformation involves a matrix X, for which an exact identity is derived, and which can be constructed either in a noniterative way or by various iteration schemes, without requiring an expansion parameter. The convergence behavior of five different iteration schemes is studied numerically, with very promising results.

Infinite-Order Regular Approximation by the Metric Perturbation

The regular approximation methods for the reduction of the Dirac equation to a fully equivalent two-component form are considered in the framework of the perturbation theory. The usual Dirac hamiltonian is first transformed with the change of metric. Then, the change of metric is considered as a perturbation to the zeroth-order (ZORA) problem. General formulae for perturbation corrections to the ZORA wave function and energy are expressed solely in terms of the two-component solutions. The method presented in this paper gives the energy- independent scheme for the step-by-step generation of the infinite-order results which are equivalent to solutions of the Dirac equations. Several formal and computational aspects of the infinite-order regular approximation are discussed. It is concluded that, because of the use of well-behaved operators, the high-order regular approximation methods can be considered as competitive to high-order Douglas-Kroll approaches.