Dirac-Fock

The Dirac equation can be readily adapted to the description of one electron in the held of the other electrons (Hartree-Fock theory). This is called a Dirac-Fock or Dirac-Hartree-Fock (DHF) calculation. 

Figure4.7 Relativistic bond contractions A re for Au2 calculated in the years from 1989 to 2001 using different quantum chemical methods. Electron correlation effects Acte = te(corn) — /"e(HF) at the relativistic level are shown on the right hand side of each bar if available. From the left to the right in chronological order Hartree-Fock-Slater results from Ziegler et al. [147] AIMP coupled pair functional results from Stbmberg and Wahlgren [148] EC-ARPP results from Schwerdtfeger [5] EDA results from Haberlen and Rdsch [149] Dirac-Fock-Slater...

Desclaux, J.P. (1973) Relativistic Dirac-Fock expectation values for atoms with Z = 1 to, Z = 120. Atomic Data and Nuclear Data Tables, 12, 311—406. 

Vilkas, M.J., Ishikawa, Y. and Koc, K. (1998) Quadratically convergent multiconfiguration Dirac-Fock and multireference relativistic configuration-interaction calculations for many-electron systems. Physical Review E, 58, 5096-5110. 

Bastug, T, Fleinemann, D., Sepp, W.-D., Kolb, D. and Fricke, B. (1993) All-electron Dirac-Fock-Slater SCF calculations of the Au2 molecule. Chemical Physics Letters, 211, 119-124. 

Pempointner, M., Seth, M. and Schwerdtfeger, P. (1998) A Point-Charge Model for the Nuclear Quadmpole Moment. Accurate Coupled-Cluster, Dirac-Fock, Douglas-KroD and Nomelativistic Hartree-Fock Calculations... 

Fig. 7.57 Plot of the isomer shift 5 of the 36.2 keV Mdssbauer transition of Os versus Dirac-Fock values for the electron density differences at the Os nuclei in free ion 5d configurations. The numbers of the data points refer to the numbering of the compounds in Table 7.9 (from [258])...

Fig. 7.78 Linear relation of the quadmpole splitting A q = ( jl)eqQ (1 + j /3)l/2 and the isomer shift b for aurous (a) and auric (b) compounds. Also included is a correlation with the relative change in electron density at the gold nucleus, Ali/r(o)P, as derived from Dirac-Fock atomic structure calculations for several electron configurations of gold. An approximate scale of the EFG (in the principal axes system) is given on the right-hand ordinate (from [341])...

The twin facts that heavy-atom compounds like BaF, T1F, and YbF contain many electrons and that the behavior of these electrons must be treated relati-vistically introduce severe impediments to theoretical treatments, that is, to the inclusion of sufficient electron correlation in this kind of molecule. Due to this computational complexity, calculations of P,T-odd interaction constants have been carried out with relativistic matching of nonrelativistic wavefunctions (approximate relativistic spinors) [42], relativistic effective core potentials (RECP) [43, 34], or at the all-electron Dirac-Fock (DF) level [35, 44]. For example, the first calculation of P,T-odd interactions in T1F was carried out in 1980 by Hinds and Sandars [42] using approximate relativistic wavefunctions generated from nonrelativistic single particle orbitals. 

Dirac-Fock (DF) level, ab initio calculations, P,T-odd interactions, 253 Dirac matrices, P,T-odd interaction operator, 251-253... 

There are two approaches to map crystal charge density from the measured structure factors by inverse Fourier transform or by the multipole method [32]. Direct Fourier transform of experimental structure factors was not useful due to the missing reflections in the collected data set, so a multipole refinement is a better approach to map charge density from the measured structure factors. In the multipole method, the crystal charge density is expanded as a sum of non-spherical pseudo-atomic densities. These consist of a spherical-atom (or ion) charge density obtained from multi-configuration Dirac-Fock (MCDF) calculations [33] with variable orbital occupation factors to allow for charge transfer, and a small non-spherical part in which local symmetry-adapted spherical harmonic functions were used. 

The relativistic coupled cluster method starts from the four-component solutions of the Drrac-Fock or Dirac-Fock-Breit equations, and correlates them by the coupled-cluster approach. The Fock-space coupled-cluster method yields atomic transition energies in good agreement (usually better than 0.1 eV) with known experimental values. This is demonstrated here by the electron affinities of group-13 atoms. Properties of superheavy atoms which are not known experimentally can be predicted. Here we show that the rare gas eka-radon (element 118) will have a positive electron affinity. One-, two-, and four-components methods are described and applied to several states of CdH and its ions. Methods for calculating properties other than energy are discussed, and the electric field gradients of Cl, Br, and I, required to extract nuclear quadrupoles from experimental data, are calculated. 

Heavy atoms exhibit large relativistic effects, often too large to be treated perturba-tively. The Schrodinger equation must be supplanted by an appropriate relativistic wave equation such as Dirac-Coulomb or Dirac-Coulomb-Breit. Approximate one-electron solutions to these equations may be obtained by the self-consistent-field procedure. The resulting Dirac-Fock or Dirac-Fock-Breit functions are conceptually similar to the familiar Hartree-Fock functions the Hartree-Fock orbitals are replaced, however, by four-component spinors. Correlation is no less important in the relativistic regime than it is for the lighter elements, and may be included in a similar manner. 

The no-pair DCB Hamiltonian (6) is used as a starting point for variational or many-body relativistic calculations [9], The procedure is similar to the nonrelativistic case, with the Hartree-Fock orbitals replaced by the four-component Dirac-Fock-Breit (DFB) functions. The spherical symmetry of atoms leads to the separation of the one-electron equation into radial and spin-angular parts [10], The radial four-spinor has the so-called large component the upper two places and the small component Q, in the lower two. The quantum number k (with k =j+ 1/2) comes from the spin-angular equation, and n is the principal quantum number, which counts the solutions of the radial equation with the same k. Defining... 

Here frs and (ri-l tM> are, respectively, elements of one-electron Dirac-Fock and antisymmetrized two-electron Coulomb-Breit interaction matrices over Dirac four-component spinors. The effect of the projection operators is now taken over by the normal ordering, denoted by the curly braces in (15), which requires annihilation operators to be moved to the right of creation operators as if all anticommutation relations vanish. The Fermi level is set at the top of the highest occupied positive-energy state, and the negative-energy states are ignored. 

Of the five group-13 elements, only B and A1 have experimentally well characterized electron affinities. Lists of recommended EAs [50,51] show errors ranging from 50% to 100% for Ga, In, and T1. Very few calculations have appeared for the latter atoms. These include the multireference configuration interaction (MRCI) ofAmau etal. using pseudopotentials [52], our relativistic coupled cluster work on T1 [45], and the multiconfiguration Dirac-Fock (MCDF) computation of Wijesundera [53]. 

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 molecular orbitals in the nonrelativistic and one-component calculations and the large component in the Dirac-Fock functions were spanned in the Cd s Ap9d)l[9slp6d basis of [63] and the H (5s 2p)/[35 l/>] set [61]. Contraction coefficients were taken from corresponding atomic SCF calculations. The basis for the small components in the Dirac-Fock calculations is derived by the MOLFDIR program from the large-component basis. The basis set superposition error is corrected by the counterpoise method [64]. The Breit interaction was found to have a very small effect and is therefore not included in the results. 

Table 1 Energies (in KeV) of single positive ions evaluated with (AH) a full relativistic kinetic energy functional without exchange [15] the c -order semi-relativistic functional (Eq. 46) without (1) and with (2) the relativistic exchange correction ((f-term), all using near-nuclear corrections, compared to Dirac-Fock (DF) values.

The practical implementation of this method leads to an integro-differential equation which has been solved for atoms [23] leading to the most consistent relativistic Thomas-Fermi estimations known up to now. The energy values for some atoms are displayed in Table 2 compared to Dirac-Fock ones. 

Table 2 Values of relativistic energies (E) and differences among relativistic and non-relativistic energies (AE) for neutral atoms in atomic units with the present approach using thefunctional given by Eq. (46) not including (1) or including (2) the term, compared to the results of Engel and Dreizler (ED) [23] using the relativistic Thomas-Fermi-Dirac- Weirsacker approach described in Section 2.6, and to Dirac-Fock values...

Results of similar accuracy as relativistic TFDW are found with a simple procedure based on near-nuclear correction which leave space for further improvements. For the reasons mentioned at the end of previous section the direct way to improve the present approach seems to be the refinement of the near nuclear corrections, a problem that we have just tackled with success in the non-relativistic framework [31,32]. The aim was to describe the near-nuclear region accurately by means of using the quantum mechanical exact asymptotic expression up to of the different ns eigenstates of Schodinger equation with a fit of the semiclassical potential at short distancies to the exact asymptotic behaviour (with four terms) of the potential near the nucleus. The result is that the density below Tq becomes very close to Hartree-Fock values and the improvement of the energy values is large (as an example, the energy of Cs+ is improved from the Ashby-Holzman result of-189.5 keV up to -205.6, very close to the HF value of -204.6 keV). This result makes us expect that a similar procedure in the relativistic framework may provide results comparable to Dirac-Fock ones. 

In this paper we present the first application of the ZORA (Zeroth Order Regular Approximation of the Dirac Fock equation) formalism in Ab Initio electronic structure calculations. The ZORA method, which has been tested previously in the context of Density Functional Theory, has been implemented in the GAMESS-UK package. As was shown earlier we can split off a scalar part from the two component ZORA Hamiltonian. In the present work only the one component part is considered. We introduce a separate internal basis to represent the extra matrix elements, needed for the ZORA corrections. This leads to different options for the computation of the Coulomb matrix in this internal basis. The performance of this Hamiltonian and the effect of the different Coulomb matrix alternatives is tested in calculations on the radon en xenon atoms and the AuH molecule. In the atomic cases we compare with numerical Dirac Fock and numerical ZORA methods and with non relativistic and full Dirac basis set calculations. It is shown that ZORA recovers the bulk of the relativistic effect and that ZORA and Dirac Fock perform equally well in medium size basis set calculations. For AuH we have calculated the equilibrium bond length with the non relativistic Hartree Fock and ZORA methods and compare with the Dirac Fock result and the experimental value. Again the ZORA and Dirac Fock errors are of the same order of magnitude. 

There are many problems in e.g. catalysis in which relativity may play a deciding role in the chemical reactivity. These problems generally involve large organic molecules which cannot be handled within the Dirac Fock framework. It is therefore necessary to reduce the work by making additional approximations. Generally used approaches are based on the Pauli expansion or on the Douglas Kroll transformation [3]. 

In this paper we use a regular approximation of the Dirac Fock formalism known as... 

It was previously shown [7] that one can, starting from the Dirac-Fock equation, derive the two component scaled ZORA equation... 

Equation (1) is obtained by using an expansion in E/ 2c - Vc) on the Dirac Fock equation. This expansion is valid even for a singular Coulombic potential near the nucleus, hence the name regular approximation. This is in contrast with the Pauli method, which uses an expansion in (E — V)I2(. Everything is written in terms of the two component ZORA orbitals, instead of using the large and small component Dirac spinors. This is an extra approximation with respect to the original formalism. 

Table 5 Comparison of equilibrium bondlengths for non relativistic ZORA and Dirac Fock calculations with the experimental value...

---