Kohn-Sham equations solution

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

Bauernschmitt, R., Ahlrichs, R., 1996a, Stability Analysis for Solutions of the Closed Shell Kohn-Sham Equation , J. Chem. Phys., 104, 9047. 

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. 

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. 

Here the vector zs(r) is constructed from the kinetic energy tensor obtained by employing the solutions of the Kohn-Sham equation in Equation 7.42. Thus it is, in general, different from the vector z(r). A comparison of Equations 7.41 and 7.46 gives... 

The difficulty of this problem can be appreciated by noticing that in order to solve the Kohn-Sham equations exactly, one must have the exact exchange-correlation potential which, moreover, must be obtained from the exact exchange-correlation functional c[p( )] given by Eq. (160). As this functional is not known, the attempts to obtain a direct solution to the Kohn-Sham equations have had to rely on the use of approximate exchange-correlation functionals. This approximate direct method, however, does not satisfy the requirement of functional iV-representability,... 

These equations are superficially similar to Eq. (1.1). The main difference is that the Kohn-Sham equations are missing the summations that appear inside the full Schrodinger equation [Eq. (1.1)]. This is because the solution of the Kohn-Sham equations are single-electron wave functions that depend on only three spatial variables, ij ,(r). On the left-hand side of the Kohn-Sham equations there are three potentials, V, VH, and Vxc- The first... 

We have skipped over a whole series of important details in this process (How close do the two electron densities have to be before we consider them to be the same What is a good way to update the trial electron density How should we define the initial density ), but you should be able to see how this iterative method can lead to a solution of the Kohn-Sham equations that is self-consistent. 

As we discussed in Chapter 1, the main aim of a DFT calculation is to find the electron density that corresponds to the ground-state configuration of the system, p(r). The electron density is defined in terms of the solutions to the Kohn-Sham equations, i /-(r), by... 

QR Method. The first relativistic method is the so-called quasi-relativistic (QR) method. It has been developed by Snijders, Ziegler and co-workers (13). In this approach, a Pauli Hamiltonian is included into the self-consistent solution of the Kohn-Sham equations of DFT. The Pauli operator is in a DFT framework given by... 

Approximate solutions of the time-dependent Schrodinger equation can be obtained by using Frenkel variational principle within the PCM theoretical framework [17]. The restriction to a one-determinant wavefunction with orbital expansion over a finite atomic basis set leads to the following time-dependent Hartree-Fock or Kohn-Sham equation ... 

The first derivative of the density matrix with respect to the magnetic induction (dPfiv/dBi) is obtained by solving the coupled-perturbed Hartree-Fock (or Kohn-Sham) equations to which the first derivative of the effective Fock (or Kohn-Sham) operator with respect to the magnetic induction contributes. Due to the use of GIAOs, specific corrections arising from the effective operator Hcnv describing the environment effects will appear. We refer to Ref. [28] for the PCM model and to Ref. [29] for the DPM within either a HF or DFT description of the solute molecule. 

There are a number of band-structure methods that make varying approximations in the solution of the Kohn-Sham equations. They are described in detail by Godwal et al. (1983) and Srivastava and Weaire (1987), and we shall discuss them only briefly. For each method, one must eon-struct Bloch functions delocalized by symmetry over all the unit cells of the solid. The methods may be conveniently divided into (1) pesudopo-tential methods, (2) linear combination of atomic orbital (LCAO) methods (3) muffin-tin methods, and (4) linear band-structure methods. The pseudopotential method is described in detail by Yin and Cohen (1982) the linear muffin-tin orbital method (LMTO) is described by Skriver (1984) the most advanced of the linear methods, the full-potential linearized augmented-plane-wave (FLAPW) method, is described by Jansen... 

This paper gives a short overview of density functional calculations mainly based on the DV-Xa approach organized as follows. A short overview of Density Functional Theory, DFT, and Kohn-Sham equations is given in section II followed by a summary of different ways of solution of the Kohn-Sham equations in Sec. III. Comparisons of results from some old and some up-to-date density functional electronic structure calculations made by our group to show applications to clusters, surfaces, adsorbates on surfaces and Ceo are given in Sec. IV. Conclusions and outlook are summarized in Sec. V. 

---