Douglas-Kroll transformation relativistic

Hess, B.A. and Kaldor, U. (2000) Relativistic all-electron coupled-cluster calculations on Au2 in the framework of the Douglas—Kroll transformation. Journal of Chemical Physics, 112, 1809-1813. 

It is clear from Ho that the Douglas-Kroll transformation makes use of a model space of relativistic free-particle spinors, and that it is defined by a perturbative expansion with the external potential as perturbation. Indeed, using the formulas given above, we get the familiar expressions for the second-order Douglas-Kroll-transformed Dirac operator, which is often dubbed Douglas-Kroll-Hess (DKH) operator... 

The Douglas-Kroll transformation [40] of the Dirac-Coulomb Hamiltonian in its implementation by HeB [41-45] leads to one of the currently most successful and popular forms of a relativistic no-pair Hamiltonian. The one-electron terms of the Douglas-Kroll-HeB (DKH) Hamiltonian have the form... 

Extended Douglas—Kroll transformations applied to the relativistic many-electron Hamiltonian... 

General two-component methods have been discussed in various chapters of the first part of this book, for instance in chapter 11 on Two-Component Methods and the Generalised Douglas-Kroll Transformation by Wolf, Reiher and Hess [165], in chapter 12 by Kutzelnigg on Perturbation Theory of Relativistic Effects [166] and in chapter 13 by Sundholm on Perturbation Theory Based on Quasi-Relativistic Hamiltonians [167]. 

A. Wolf, M. Reiher, B. Hess, Two-component methods and the generalized Douglas-Kroll transformation, in P. Schwerdtfeger (Ed.), Relativistic Electronic Structure Theory, Part 1, Fundamentals, Elsevier, Netherlands, 2002, pp. 627-668. 

Before leaving the theoretical formalism section, it is important to note that perturbation theory for relativistic effects can also be done at the fo n-con onent level, i.e. before elimination of the small component by a Foldy-Wouthitysen (FW) or Douglas-Kroll transformation. This is best done with direct perturbation theory (DPT) [71]. DPT involves a change of metric in the Dirac equation and an expansion of this modified Dirac eqtiation in powers of c . The four-component Levy-Leblond equation is the appropriate nonrelativistic limit. Kutzelnigg [72] has recently worked out in detail the simultaneous DPT for relativistic effects and magnetic fields (both external and... 

The reduction of the relativistic many-electron hamiltonian by expansion in powers of the external field is the second-order Douglas-Kroll transformation [29], and has been used with success by Hess and co-workers [30]. The operators which result from this transformation are non-singular, but the integrals over the resulting operators are complicated and have to be approximated, even for finite basis set expansions. The reduction of the Dirac-Coulomb-Breit equation to two-component form using direct perturbation theory has been described by Kutzelnigg and coworkers [26, 27, 31], Rutkowski [32], and van Lenthe et al. [33]. 

T. Nakajima, K. Hirao. Extended Douglas-Kroll transformations applied to the relativistic many-electron Hanrilto-nian. /. Chem. Phys., 119(8) (2003) 4105-4111. 

U. Kaldor and B. A. Hess, Ghent. Phys. Lett., 230, 1 (1994). Relativistic All-Electron Coupled-Cluster Calculations on the Gold Atom and Gold Hydride in the Framework of the Douglas-Kroll Transformation. 

This is the equivalent of the second-order Douglas-KroU operator, but it only involves operators that have been defined in the free-particle Foldy-Wouthuysen transformation. As for the Douglas-Kroll transformed Hamiltonian, spin separation may be achieved with the use of the Dirac relation to define a spin-fi ee relativistic Hamiltonian, and an approximation in which the transformation of the two-electron integrals is neglected, as in the Douglas-Kroll-Hess method, may also be defined. Implementation of this approximation can be carried out in the same way as for the Douglas-Kroll approximation both approximations involve the evaluation of kinematic factors, which may be done by matrix methods. 

Here we see a modified kinetic energy term that is cut off near the nuclei, a spin-free relativistic correction, and a spin-orbit term, both of which are regularized and behave as 1/r for small r. We may compare this with the regularization of the free-particle Foldy-Wouthuysen or Douglas-Kroll transformed Hamiltonian of section 16.3. The regularization clearly corrects the overestimation of relativistic effects that plagues the Pauli Hamiltonian. There is another consequence of the small r behavior. Since the relativistic correction operator behaves like 1 jr, it ought to be possible to use T zoRA variationally—and in fact we may demonstrate that there is a variational lower bound. 

We could of course proceed as we did for the Douglas-Kroll transformation and use a transformation that depends only on the nuclear potential. This would remove the awkwardness of having 1 /r,y in the denominator, but we still have the product of c f(2mc — Vi) with 1 /r,-y to deal with. If we are only interested in spin-free relativistic effects, we could neglect the transformation of the electron-electron interaction, as we did in the Douglas-Kroll-Hess approximation. This approximation yields the Hamiltonian... 

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. 

The Douglas-Kroll (DK) approach [153] can decouple the large and small components of the Dirac spinors in the presence of an external potential by repeating several unitary transformations. The DK transformation is a variant of the FW transformation [141] and adopts the external potential Vg t an expansion parameter instead of the speed of light, c, in the FW transformation. The DK transformation correct to second order in the external potential (DK2) has been extensively studied by Hess and co-workers [154], and has become one of the most familiar quasi-relativistic approaches. Recently, we have proposed the higher order DK method and applied the third-order DK (DK3) method to several systems containing heavy elements. 

These problems can be solved if one starts from the (untruncated ) Foldy-Wouthuysen transformation for a free particle, the only case for which the transformation is known anal3d ically, and incorporates the effects of the external potential on top. Along these lines, the so-called Douglas-Kroll-HeB (DKH) method [61-64] is constructed which is probably the most successful quasi-relativistic method in wave function based quantum chemistry. No details will be given here since this topic has been extensively discussed in volume 1 [34] of this series. Meanwhile several density functional implementations exist based on the Douglas-Kroll-HeJ3 approach [39-45]. In recent years. 

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]. 

A self-consistent scalar-relativistic (SR) version of the LCGTO-DF method has also been developed recently." "The SR variant employs a unitary second-order Douglas-Kroll-Hess (DKH) "" transformation for decoupling large and small components of the full four-component spinor solutions to the Dirac-Kohn-Sham equation. The approximate DKH transformation, very appropriate and efficient for molecular calculations, has been implemented this variant utilizes nuclear potential-based projectors and leaves the electron-electron interaction untransformed. 

---