Extended Kohn-Sham method

STRUTINSKY S SHELL-CORRECTION METHOD IN THE EXTENDED KOHN-SHAM SCHEME APPLICATION TO THE IONIZATION POTENTIAL, ELECTRON AFFINITY, ELECTRONEGATIVITY AND CHEMICAL HARDNESS OF ATOMS... 

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. 

Delchev, Y.I., A.I. Kuleff, J. Maruani, T. Mineva, and F. Zahariev. 2006. Strutinsky s shell-correction method in the extended Kohn-Sham scheme Application to the ionization potential, electron affinity, electroneg ativity and chemical hardness of atoms. In Recent Advances in the Theory of Chemical and Physical Systems, edited by Jean-Pierre Julien, Jean Maruani, and Didier Mayou, 159-177. New York Springer-Verlag. 

Lately, the CP-MD approach has been combined with a mixed QM/MM scheme [10-12] which enables the treatment of chemical reactions in biological systems comprising tens of thousands of atoms [11, 26]. Furthermore, CP-MD and mixed QM/MM CP-MD simulations have also been extended to the treatment of excited states within a restricted open-shell Kohn-Sham approach [16, 17, 27] or within a linear response formulation of TDDFT [16, 18], enabling the study of biological photoreceptors [28] and the in situ design of optimal fluorescence probes with tailored optical properties [32]. Among the latest extensions of this method are also the calculation of NMR chemical shifts [14]. 

Electronegativity (y) or chemical potential (ji), ionization potential (7), and electron affinity (A) can be computed for electronic systems from the Kohn Sham (KS) equation, which has been extended by Janak [42] and others [43 45] using the Xu method [46], In this approach, one gets a meaning for orbital energy as ... 

The relativistic correction for the kinetic energy in the Dirac equation is naturally applicable to the Kohn-Sham equation. This relativistic Kohn-Sham equation is called the Dirac-KohnSham equation (Rajagopal 1978 MacDonald and Vosko 1979). The Dirac-Kohn-Sham equation is founded on the Rajagopal-Callaway theorem, which is the relativistic expansion of the Hohenberg-Kohn theorem on the basis of QED (Rajagopal and Callaway 1973). In this theorem, two theorems are contained The first theorem proves that the four-component external potential, which is the vector-potential-extended external potential, is determined by the four-component current density, which is the current-density-extended electron density. On the other hand, the second theorem establishes the variational principle for every four-component current density. See Sect. 6.5 for vector potential and current density. Consequently, the solution of the Dirac-Kohn-Sham equation is represented by the four-component orbital. This four-component orbital is often called a molecular spinor. However, this name includes no indication of orbital, which is the solution of one-electron SCF equations moreover, the targets of the calculations are not restricted to molecules. Therefore, in this book, this four-component orbital is called an orbital spinor. The Dirac-Kohn-Sham wavefunction is represented by the Slater determinant of orbital spinors (see Sect. 2.3). Following the Roothaan method (see Sect. 2.5), orbital spinors are represented by a linear combination of the four-component basis spinor functions, Xp, ... 

---