Density functional Kohn-Sham equations

The total electron density is just the sum of the densities for the two types of electron. The exchange-correlation functional is typically different for the two cases, leading to a set of spin-polarised Kohn-Sham equations ... 

The application of density functional theory to isolated, organic molecules is still in relative infancy compared with the use of Hartree-Fock methods. There continues to be a steady stream of publications designed to assess the performance of the various approaches to DFT. As we have discussed there is a plethora of ways in which density functional theory can be implemented with different functional forms for the basis set (Gaussians, Slater type orbitals, or numerical), different expressions for the exchange and correlation contributions within the local density approximation, different expressions for the gradient corrections and different ways to solve the Kohn-Sham equations to achieve self-consistency. This contrasts with the situation for Hartree-Fock calculations, wlrich mostly use one of a series of tried and tested Gaussian basis sets and where there is a substantial body of literature to help choose the most appropriate method for incorporating post-Hartree-Fock methods, should that be desired. 

Kohn-Sham equations of the density functional theory then take on the following... 

Note that in all current implementations of TDDFT the so-called adiabatic approximation is employed. Here, the time-dependent exchange-correlation potential that occurs in the corresponding time-dependent Kohn-Sham equations and which is rigorously defined as the functional derivative of the exchange-correlation action Axc[p] with respect to the time-dependent electron-density is approximated as the functional derivative of the standard, time-independent Exc with respect to the charge density at time t, i. e.,... 

Morrison, R. C., Zhao, Q., 1995, Solution to the Kohn-Sham Equations Using Reference Densities from Accurate, Correlated Wave Functions for the Neutral Atoms Hehum Through Argon , Phys. Rev. A, 51, 1980. 

Scheiner, A. C., Baker, J., Andzelm, J. W., 1997, Molecular Energies and Properties from Density Functional Theory Exploring Basis Set Dependence of Kohn-Sham Equation Using Several Density Functionals , J. Comput. 

Wesolowski, T. A., H. Chermette, and J. Weber. Accuracy of Approximate Kinetic Energy Functionals in the Model of Kohn-Sham Equations with Constrained Electron Density the FH-NCH complex as a Test Case. J. Chem. Phys. In press. 

Morrison, R. C., Q. Zhao, R. C. Morrison, and R. G. Parr. 1995. Solution of the Kohn-Sham equations using reference densities from accurate, correlated wave functions for the neutral atoms helium through argon. Phys. Rev. A51, 1980. 

In the Kohn-Sham equation above, the Coulomb potential and the XC potential are obtained from their energy counterparts by taking the functional derivative of the latter with respect to the density. Thus... 

Most practical electronic structure calculations using density functional theory [1] involve solving the Kohn-Sham equations [2], The only unknown quantity in a Kohn-Sham spin-density functional calculation is the exchange-correlation energy (and its functional derivative) [2]... 

If you have a vague sense that there is something circular about our discussion of the Kohn-Sham equations you are exactly right. To solve the Kohn-Sham equations, we need to define the Hartree potential, and to define the Hartree potential we need to know the electron density. But to find the electron density, we must know the single-electron wave functions, and to know these wave functions we must solve the Kohn-Sham equations. To break this circle, the problem is usually treated in an iterative way as outlined in the following algorithm 1 2 3... 

Solve the Kohn-Sham equations defined using the trial electron density to find the single-particle wave functions, i i,(r). 

In fact, the true form of the exchange-correlation functional whose existence is guaranteed by the Hohenberg-Kohn theorem is simply not known. Fortunately, there is one case where this functional can be derived exactly the uniform electron gas. In this situation, the electron density is constant at all points in space that is, n(r) = constant. This situation may appear to be of limited value in any real material since it is variations in electron density that define chemical bonds and generally make materials interesting. But the uniform electron gas provides a practical way to actually use the Kohn-Sham equations. To do this, we set the exchange-correlation potential at each position to be the known exchange-correlation potential from the uniform electron gas at the electron density observed at that position ... 

This approximation uses only the local density to define the approximate exchange-correlation functional, so it is called the local density approximation (LDA). The LDA gives us a way to completely define the Kohn-Sham equations, but it is crucial to remember that the results from these equations do not exactly solve the true Schrodinger equation because we are not using the true exchange-correlation functional. 

We now have all the pieces in place to perform an HF calculation—a basis set in which the individual spin orbitals are expanded, the equations that the spin orbitals must satisfy, and a prescription for forming the final wave function once the spin orbitals are known. But there is one crucial complication left to deal with one that also appeared when we discussed the Kohn-Sham equations in Section 1.4. To find the spin orbitals we must solve the singleelectron equations. To define the Hartree potential in the single-electron equations, we must know the electron density. But to know the electron density, we must define the electron wave function, which is found using the individual spin orbitals To break this circle, an HF calculation is an iterative procedure that can be outlined as follows ... 

In the remainder of this section, we give a brief overview of some of the functionals that are most widely used in plane-wave DFT calculations by examining each of the different approaches identified in Fig. 10.2 in turn. The simplest approximation to the true Kohn-Sham functional is the local density approximation (LDA). In the LDA, the local exchange-correlation potential in the Kohn-Sham equations [Eq. (1.5)] is defined as the exchange potential for the spatially uniform electron gas with the same density as the local electron density ... 

A critical feature of this quantity is that it is nonlocal, that is, a functional based on this quantity cannot be evaluated at one particular spatial location unless the electron density is known for all spatial locations. If you look back at the Kohn-Sham equations in Chapter 1, you can see that introducing this nonlocality into the exchange-correlation functional creates new numerical complications that are not present if a local functional is used. Functionals that include contributions from the exact exchange energy with a GGA functional are classified as hyper-GGAs. 

Kohn-Sham Equations and Density Functional Models... 

Z is the nuclear charge, R-r is the distance between the nucleus and the electron, P is the density matrix (equation 16) and (qv Zo) are two-electron integrals (equation 17). f is an exchange/correlation functional, which depends on the electron density and perhaps as well the gradient of the density. Minimizing E with respect to the unknown orbital coefficients yields a set of matrix equations, the Kohn-Sham equations , analogous to the Roothaan-Hall equations (equation 11). 

Kohn-Sham Equations. The set of equations obtained by applying the Local Density Approximation to a general multi-electron system. An Exchange/Correlation Functional which depends on the electron density has replaced the Exchange Energy expression used in the Hartree-Fock Equations. The Kohn-Sham equations become the Roothaan-Hall Equations if this functional is set equal to the Hartree-Fock Exchange Energy expression. 

After the discovery of the relativistic wave equation for the electron by Dirac in 1928, it seems that all the problems in condensed-matter physics become a matter of mathematics. However, the theoretical calculations for surfaces were not practical until the discovery of the density-functional formalism by Hohenberg and Kohn (1964). Although it is already simpler than the Hartree-Fock formalism, the form of the exchange and correlation interactions in it is still too complicated for practical problems. Kohn and Sham (1965) then proposed the local density approximation, which assumes that the exchange and correlation interaction at a point is a universal function of the total electron density at the same point, and uses a semiempirical analytical formula to represent such universal interactions. The resulting equations, the Kohn-Sham equations, are much easier to handle, especially by using modern computers. This method has been the standard approach for first-principles calculations for solid surfaces. 

Density functional theory, 21, 31, 245-246 B3LYP functional, 246 Hartree-Fock-Slater exchange, 246 Kohn-Sham equations, 245 local density approximation, 246 nonlocal corrections, 246 Density matrix, 232 Determinantal wave function, 23 Dewar benzene, 290 from acetylene + cyclobutadiene, 290 interaction diagram, 297 rearrangement to benzene, 290, 296-297 DFT, see Density functional theory... 

In the Time Dependent Density Functional Theory (TDDFT) [16], the correlated many-electron problem is mapped into a set of coupled Schrodinger equations for each single electronic wavefunctions (o7 (r, t),j= 1, ), which yields the so-called Kohn-Sham equations (in atomic units)... 

---