Kohn Sham

Kohn-Sham or Slater exchange was more accurate for realistic systems [H]. Slater suggested that a parameter be introduced that would allow one to vary the exchange between the Slater and Kolm-Sham values [19]. The parameter, a, was often... 

Once a solution of the Kohn-Sham equation is obtained, the total energy can be computed from... 

The wavevector is a good quantum number e.g., the orbitals of the Kohn-Sham equations [21] can be rigorously labelled by k and spin. In tln-ee dimensions, four quantum numbers are required to characterize an eigenstate. In spherically syimnetric atoms, the numbers correspond to n, /, m., s, the principal, angular momentum, azimuthal and spin quantum numbers, respectively. Bloch s theorem states that the equivalent... 

Kleinman L 1997 Significance of the highest occupied Kohn-Sham eigenvalue Phys. Rev. B 56 12 042-5... 

By introducing this expression for the electron density and applying the appropriate variational condition the following one-electron Kohn-Sham equations result ... 

To. solve the Kohn-Sham equations a self-consistent approach is taken. An initial guess of the density is fed into Equation (3.47) from which a set of orbitals can be derived, leading to an improved value for the density, which is then used in the second iteration, and so on until convergence is achieved. 

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

In addition to the energy terms for the exchange-correlation contribution (which enables the total energy to be determined) it is necessary to have corresponding terms for the potential, Vxc[p(i )]/ which are used to solve the Kohn-Sham equations. These are obtained as the appropriate first derivatives using Equation (3.52). 

To solve the Kohn-Sham equations a number of different approaches and strategies have been proposed. One important way in which these can differ is in the choice of basis set for expanding the Kohn-Sham orbitals. In most (but not all) DPT programs for calculating the properties of molecular systems (rather than for solid-state materials) the Kohn-Sham orbitals are expressed as a linear combination of atomic-centred basis functions ... 

Equation (3.74) is the exact exchange energy (obtained from the Slater determinant the Kohn-Sham orbitals), is the exchange energy under the local spin densit) ... 

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

Density Eunctional Methods. The Kohn-Sham equations are... 

Properties can be computed by finding the expectation value of the property operator with the natural orbitals weighted by the occupation number of each orbital. This is a much faster way to compute properties than trying to use the expectation value of a multiple-determinant wave function. Natural orbitals are not equivalent to HF or Kohn-Sham orbitals, although the same symmetry properties are present. 

In this formulation, the electron density is expressed as a linear combination of basis functions similar in mathematical form to HF orbitals. A determinant is then formed from these functions, called Kohn-Sham orbitals. It is the electron density from this determinant of orbitals that is used to compute the energy. This procedure is necessary because Fermion systems can only have electron densities that arise from an antisymmetric wave function. There has been some debate over the interpretation of Kohn-Sham orbitals. It is certain that they are not mathematically equivalent to either HF orbitals or natural orbitals from correlated calculations. However, Kohn-Sham orbitals do describe the behavior of electrons in a molecule, just as the other orbitals mentioned do. DFT orbital eigenvalues do not match the energies obtained from photoelectron spectroscopy experiments as well as HF orbital energies do. The questions still being debated are how to assign similarities and how to physically interpret the differences. 

Kohn-Sham orbitals functions for describing the electron density in density functional theory calculations... 

In this equation Exc is the exchange correlation functional [46], is the partial charge of an atom in the classical region, Z, is the nuclear charge of an atom in the quantum region, is the distance between an electron and quantum atom q, r, is the distance between an electron and a classical atom c is the distance between two quantum nuclei, and r is the coordinate of a second electron. Once the Kohn-Sham equations have been solved, the various energy terms of the DF-MM method are evaluated as... 

In actual practice, self-consistent Kohn-Sham DFT calculations are performed in an iterative manner that is analogous to an SCF computation. This simiBarity to the methodology of Hartree-Fock theory was pointed out by Kohn and Sham. 

The Kohn-Sham equations look like standard HF equations, except that the exchange term is replaced with an exchange-correlation potential whose form is unknown. 

Implementation of the Kohn-Sham-LCAO procedure is quite simple we replace the standard exchange term in the HF-LCAO expression by an appropriate Vxc that will depend on the local electron density and perhaps also its gradient. The new integrals involved contain fractional powers of the electron density and cannot be evaluated analytically. There are various ways forward, all of which... 

As mentioned above, a KS-LCAO calculation adds one additional step to each iteration of a standard HF-LCAO calculation a quadrature to calculate the exchange and correlation functionals. The accuracy of such calculations therefore depends on the number of grid points used, and this has a memory resource implication. The Kohn-Sham equations are very similar to the HF-LCAO ones and most cases converge readily. 

---