Relativistic computations, Douglas-Kroll

Molecules are more difficult to treat accurately than atoms, because of the reduced symmetry. An additional complication arises in relativistic calculations the Dirac-Fock-(-Breit) orbitals will in general be complex. One way to circumvent this difficulty is by the Douglas-Kroll-Hess transformation [57], which yields a one-component function with computational effort essentially equal to that of a nonrelativistic calculation. Spin-orbit interaction may then be added as a perturbation, implementation to AuH and Au2 has been reported [58]. Progress has also been made in the four-component formulation [59], and the MOLFDIR package [60] has been extended to include the CC method. Application to SnH4 has been described [61] here we present a recent calculation of several states of CdH and its ions [62], with one-, two-, and four-component methods. 

Accounting for relativistic effects in computational organotin studies becomes complicated, because Hartree-Fock (HF), density functional theory (DFT), and post-HF methods such as n-th order Mpller-Plesset perturbation (MPn), coupled cluster (CC), and quadratic configuration interaction (QCI) methods are non-relativistic. Relativistic effects can be incorporated in quantum chemical methods with Dirac-Hartree-Fock theory, which is based on the four-component Dirac equation. " Unformnately the four-component Flamiltonian in the all-electron relativistic Dirac-Fock method makes calculations time consuming, with calculations becoming 100 times more expensive. The four-component Dirac equation can be approximated by a two-component form, as seen in the Douglas-Kroll (DK) Hamiltonian or by the zero-order regular approximation To address the electron cor-... 

Operators that result from a DK transformation are directly given in the momentum representation. Hess et al. [29,31] developed a very efficient strategy to evaluate the corresponding matrix elements in a basis set representation it employs the eigenvectors of the operator as approximate momentum representation [29,31]. In practice, the two-component DK Hamiltonian is built of matrix representations of the three operators p, V, pVp + id(pV x p). This Douglas-Kroll-Hess (DKH) approach became one of the most successful two-component tools of relativistic computational chemistry [16,74]. In particular, many applications showed that the second-order operator 2 Is variationally stable [10,13,14,31,75,76,87]. 

The scalar relativistic (SR) corrections were calculated by the second-order Douglas-Kroll-Hess (DKH2) method [53-57] at the (U/R)CCSD(T) or MRCI level of theory in conjunction with the all-electron aug-cc-pVQZ-DK2 basis sets that had been recently developed for iodine [58]. The SR contributions, as computed here, account for the scalar relativistic effects on carbon as well as corrections for the PP approximation for iodine. Note, however, that the Stuttgart-Koln PPs that are used in this work include Breit corrections that are absent in the Douglas-Kroll-Hess approach [58]. 

Very recently McGrath [113] presented torsional ab initio DMSs of H2O2 computed using high levels of theory, CCSD(T)/aug-cc-pVNZ (N = D,T,Q) approximated to the CBS limit. The core electrons correlations were estimated with aug-cc-pwCVTZ and the relativistic contributions with Douglas-Kroll-Hess (DKH) theory. Only the torsional part of the DMSs was investigated and represented with the following Fourier series analytical function... 

---