Density functionals relativistic corrections

The expressions (4.22)-(4.23) found in chap. 4 for the isomer shift 5 in nonrelativ-istic form may be applied to lighter elements up to iron without causing too much of an error. In heavier elements, however, the wave function j/ is subject to considerable modification by relativistic effects, particularly near the nucleus (remember that the spin-orbit coupling coefficient increases with Z ). Therefore, the electron density at the nucleus l /(o)P will be modified as well and the aforementioned equations for the isomer shift require relativistic correction. This has been considered [1] in a somewhat restricted approach by using Dirac wave functions and first-order perturbation theory in this approximation the relativistic correction simply consists of a dimensionless factor S (Z), which is introduced in the above equations for S,... 

On the other hand, high-level computational methods are limited, for obvious reasons, to very simple systems.122 Calculations are likely to have limited accuracy due to basis set effects, relativistic contributions, and spin orbit corrections, especially in the case of tin hydrides, but these concerns can be addressed. Given the computational economy of density functional theories and the excellent behavior of the hybrid-DFT B3LYP123 already demonstrated for calculations of radical energies,124 we anticipate good progress in the theoretical approach. We hope that this collection serves as a reference for computational work that we are certain will be forthcoming. 

Atomic units will be used throughout. The explicit density functionals representing the different contributions to the energy from the different terms of the hamiltonian are found performing expectation values taking Slater determinants of local plane waves as in the standard Fermi gas model. Those representing the first relativistic corrections are calculated in the Appendix. 

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. 

In subsection 3.1, we will present GGA and LDA calculations for Au clusters with 6first principles method outlined in section 2, which employs the same scalar-relativistic pseudo-potential for LDA and GGA (see Fig 1). These calculations show the crucial relevance of the level of density functional theory (DFT), namely the quality of the exchange-correlation functional, to predict the correct structures of Au clusters. Another, even more critical, example is presented in subsection 3.2, where we show that both approaches, LDA and GGA, predict the cage-like tetrahedral structure of Au2o as having lower energy than amorphous-like isomers, whereas for other Au clusters, namely Auig, Au ... 

The use of computational chemistry to address issues relative to process design was discussed in an article. The need for efficient software for massively parallel architectures was described. Methods to predict the electronic structure of molecules are described for the molecular orbital and density functional theory approaches. Two examples of electronic stracture calculations are given. The first shows that one can now make extremely accurate predictions of the thermochemistry of small molecules if one carefully considers all of the details such as zero-point energies, core-valence corrections, and relativistic corrections. The second example shows how more approximate computational methods, still based on high level electronic structure calculations, can be used to address a complex waste processing problem at a nuclear production facility (Dixon and Feller, 1999). 

In the previous section we presented the semi-classical electron-electron interaction we treated the electrons quantum mechanically but assumed that they interact via classical electromagnetic fields. The Breit retardation is only an approximate treatment of retardation and we shall now consider a more consistent treatment of the electron-electron interaction operator that also provides a bridge to relativistic DFT, which is current-density functional theory. For the correct description we have to take the quantization of electromagnetic fields into account (however, we will discuss only old, i.e., pre-1940 quantum electrodynamics). This means the two moving electrons interact via exchanged virtual photons with a specific angular frequency u>... 

The idea of density functional theory is that solely the electron density provides all necessary information about a system of electrons. As we have seen, when considering the interaction of moving charged fermions, it is the 4-current rather than the electron density p that has to be considered. In 1964, Hohenberg and Kohn (24) stated that in a non-relativistic framework the ground-state energy of a system can be calculated from its electron density and any trial density yields an energy that is an upper bound to the correct solution. 

All-electron DFT calculations were performed using the DMOL [24] code. These incorporated scalar relativistic corrections and employed the non-local exchange and correlation functional Perdew-Wang91 [25] denoted GGA in the rest of the paper, which is generally found to be superior to the local density approximation (EDA)... 

Density functional calculations (X) on the strengths of multiple metal-metal bonds in binuclear complexes such as M2Cl4(PR3)4 and M2X6, where the bonding primarily involves nd-orbitals rather than (n+l)s orbitals,reveal (XX) on the other hand only modest contributions from relativistic corrections (6-10 Kcal mol ) for the 5d homologues.The calculated order of stability for the M-M... 

---