Kohn Sham methods

In practical calculations making use of the Kohn-Sham method, the Kohn-Sham equation is used. This equation is a one-electron SCF equation applying the Slater determinant to the wavefunction of the Hartree method, similarly to the Hartree-Fock method. Therefore, in the same manner as the Hartree-Fock equation, this equation is derived to determine the lowest energy by means of the Lagrange multiplier method, subject to the normalization of the wavefunction (Parr and Yang 1994). As a consequence, it gives a similar Fock operator for the nonlinear equation. 

The difference between this Fock operator and the Hartree-Fock counterpart in Eq. (2.51) is only the exchange-correlation potential functional, Exc, which substitutes for the exchange operator in the Hartree-Eock operator. That is, in the electron-electron interaction potential, only the exchange operator is replaced with the approximate potential density functionals of the exchange interactions and electron correlations, while the remaining Coulomb operator, Jj, which is represented as the interaction of electron densities, is used as is. The point is that the electron correlations, which are incorporated as the interactions between electron configurations in wavefunction theories (see Sect. 3.3), are simply included 

Differently to the Hartree-Fock method, the total electronic energy in the Kohn-Sham method is generally calculated using not the exchange-correlation potential functional but the exchange-correlation energy functional, as 

The exchange-correlation potential functional, Fxc, is the first derivative of this exchange-correlation energy functional with respect to electron density. 

Similarly to the Hartree-Fock method, the Kohn-Sham method is solved by the 

---

Filatov, M., Shaik, S., 1999, Application of Spin-Restricted Open-Shell Kohn-Sham Method to Atomic and Molecular Multiplet States , J. Chem. Phys., 110, 116. 

Note that the Kohn-Sham Hamiltonian hKS [Eq. (4.1)] is a local operator, uniquely determined by electron density15. This is the main difference with respect to the Hartree-Fock equations which contain a nonlocal operator, namely the exchange part of the potential operator. In addition, the KS equations incorporate the correlation effects through Vxc whereas they are lacking in the Hartree-Fock SCF scheme. Nevertheless, though the latter model cannot be considered a special case of the KS equations, there are some similarities between the Hartree-Fock and the Kohn-Sham methods, as both lead to a set of one-electron equations allowing to describe an n-electron system. 

It is important to keep in mind that all global quantities appearing in the equations above can be evaluated within the spin-polarized Kohn-Sham method. All the working equations to calculate them depend on the eigenvalues of the frontier Kohn-Sham spin-down and up orbitals. The expressions can be consulted in Refs. [12,18]. 

Filatov M, Shaik S (1999) A spin-restricted ensamble-referenced Kohn-Sham method and its application to diradical situations, Chem Phys Lett, 304 429-437... 

Density functional theory (DFT) provides an efficient method to include correlation energy in electronic structure calculations, namely the Kohn-Sham method 1 in addition, it constitutes a solid support to reactivity models.2 DFT framework has been used to formalize empirical reactivity descriptors, such as electronegativity,3 hardness4 and electrophilicity index.5 The frontier orbital theory was generalized by the introduction of Fukui function,6 and new reactivity parameters have also been proposed.7,8 Moreover, relationships between those parameters have been found, and general methods to relate new quantities exist.9... 

The purpose of this chapter will be to review the fundamentals of ab initio MD. We will consider here Density Functional Theory based ab initio MD, in particular in its Car-Parrinello version. We will start by introducing the basics of Density Functional Theory and the Kohn-Sham method, as the method chosen to perform electronic structure calculation. This will be followed by a rapid discussion on plane wave basis sets to solve the Kohn-Sham equations, including pseudopotentials for the core electrons. Then we will discuss the critical point of ab initio MD, i.e. coupling the electronic structure calculation to the ionic dynamics, using either the Born-Oppenheimer or the Car-Parrinello schemes. Finally, we will extend this presentation to the calculation of some electronic properties, in particular polarization through the modern theory of polarization in periodic systems. 

The key point of the Kohn-Sham method is to consider an auxiliary system of N non-interacting electrons for estimating the kinetic energy of the real (interacting) system [120-122,125]. 

With the development of GGA functionals, description of molecular systems with the Kohn-Sham method reached a precision similar to other quantum theory methods. It was quickly shown that the GGA s could also well reproduce the hydrogen bond properties. Short after, liquid water at ambient condition was first simulated by Car-Parrinello MD, with a sample of 32 water molecules with periodic boundary conditions [31]. Since then, many simulations of liquid water at different temperatures and pressures and of water solutions have been performed [32-39]. Nowadays, Car-Parrinello MD has become a major tool for the study of aqueous solutions [40-64]. 

Kinetic solvent effect (KSE), on antioxidant activity 877-881 Knoevenagel condensation 1405 Kohn-Sham method 68 Koilands 1427 Koilates 1427... 

To motivate the Kohn-Sham method, we return to molecular Hamiltonian [Eq. (2)] and note that, were it not for the electron-electron repulsion terms coupling the electrons, we could write the Hamiltonian operator as a sum of one-electron operators and solve Schrodinger equation by separation of variables. This motivates the idea of replacing the electron-electron repulsion operator by an average local representation thereof, w(r), which we may term the internal potential. The Hamiltonian operator becomes... 

The Kohn-Sham wave function, KS, is not expected to be a good approximation to the exact wave function indeed, it is a worse approximation to the exact wave function than the Hartree-Fock wave function. Flowever, unlike the electron density obtained from the Flartree-Fock equations, the Kohn-Sham method yields, in principle, the exact electron density. Thus we do not need to use the Kohn-Sham wave function to compute the properties of chemical systems. Rather, motivated by the first Hohenberg-Kohn theorem, we compute properties directly from the Kohn Sham electron density. How one does this, for any given system and for any property of interest, is an active topic of research. 

Because of its critical role in constructing the potential energy surface [cf. Eq. (3)] for a molecule, thence in the prediction of molecular structure and chemical reactivity, we mention how one may compute the electronic energy of a system using the Kohn-Sham method. In particular, one has... 

No useful explicit form for Exc[p] is known, but approximate exchange-correlation energy functionals often provide an excellent approximation to the energetic properties of the molecule. To explain how approximate exchange-correlation functionals work (and when they fail), recall the essence of the Kohn-Sham method the Kohn-Sham method constructs the model system with the energy functional... 

---