Kohn-Sham-like equations

Next, using the concept [2,64] of adiabatic connection, Kohn-Sham-like equations can be derived. We suppose the existence of a continuous path between the interacting and the noninteracting systems. The density , of the th electron state is the same along the path. 

To solve the Kohn-Sham-like equation (Equation 9.15), one has to find an approximation to the potential of the noninteracting system. The OPM [25] can be... 

Both the noninteracting Hamiltonian //(, and the Kohn-Sham-like potential w ( iij, / <,] r) v ofr) depend on / they are different for different excited states. The Kohn-Sham-like equations read... 

N Orthonormal Orbitals Kohn-Sham-Like Equations. 207... 

Let us finish this Section by discussing the relationship between the Kohn-Sham-like equations advanced above and the actual Kohn-Sham equations. From the perspective of local-scaling transformations, we can analyze this relationship as follows. First of all, we assume that, for an interacting system, we are able to select an orbit-generating wavefunction belonging to the... 

Euler-Lagrange equations for the intra-orbit optimization of N orthonormal orbitals Kohn-Sham-like equations... 

The single-particle orbitals are calculated from time-dependent Schrodinger- or Kohn-Sham-like equations. 

Using the concept of adiabatic connection Kohn-Sham-like equations can be derived. It is supposed that the density is the same for both the interacting and non-interacting systems, and there exists a continuous path between them. A coupling constant path is defined by the Schrodinger equation... 

The ground-state differential virial theorem for spherically symmetric Kohn-Sham potential for particles having zero angular momentum has been derived by March and Young and generalized by Nagy and March A similar derivation can be performed for a single excited state The Kohn-Sham-like equations take the form... 

These total energy expressions, after a variational treatment, lead to DFTB Kohn-Sham like equations of the form ... 

Here, C is the gauge constant, / is the boundary of the closed shells n > f indicating the vacant band and the upper continuum electron states matomic core and the states of a negative continuum (accounting for the electron vacuum polarization). The minimization of the functional ImSEninv leads to the Dirac-Kohn-Sham-like equations for the electron density that are numerically solved. Finally an optimal set of the IQP functions results. In concrete calculation it is sufficient to use the simplified procedure, which is reduced to the functional minimization using the variation of the correlation potential parameter b in Eq. 3.11 [20, 32]. The Dirac equations for the radial functions F and G (the large and small Dirac components respectively) are ... 

Section 11.2 presents the non-variational theory. In Sect. 11.3 Kohn-Sham-like equations are obtained through adiabatic connection. Density scaling is applied to obtain a generalized Kohn-Sham scheme in Sect. 11.4. The optimized potential and the KLI methods are generalized in Sect. 11.5. The last section is devoted to illustrative examples and discussion. 

N is always integer, but is generally non-integer. Therefore, the Kohn-Sham-like equations will differ from the ones corresponding to the Af-electron Kohn-Sham system (11.15). To construct another Kohn-Sham system we define the density... 

Applying the variational principle to the energy given by Eq. 1, Kohn and Sham reformulated the density functional theory by deriving a set of one-electron Hartree-like equations leading to the Kohn-Sham orbitals v().(r) involved in the calculation of p(r)15. The Kohn-Sham (KS) equations are written as follows ... 

The set of Kohn-Sham-like linear equations above represents the working equations of DFPT. They are usually solved by iterative linear algebra algorithms (conjugate-gradient minimization). 

In this short review, a brief overview of the underlying principles of TDDFT has been presented. The formal aspects for TDDFT in the presence of scalar potentials with periodic time dependence as well as TD electric and magnetic fields with arbitrary time dependence are discussed. This formalism is suitable for treatment of interaction with radiation in atomic and molecular systems. The Kohn-Sham-like TD equations are derived, and it is shown that the basic picture of the original Kohn-Sham theory in terms of a fictitious system of noninteracting particles is retained and a suitable expression for the effective potential is derived. 

The partial cancellation of correlation effects arising from the 1- and 2-matrices is a well-known fact in Kohn-Sham theory [1, 2, 72], As is to be expected, it also appears in the Kohn-Sham-type equations for an orbit cty. We would like to emphasize, nevertheless, that all terms in Eq. (82) can be explicitly calculated. In principle, therefore, within a particular orbit cty, the exchange-correlation energy term as well as the Kohn-Sham-type exchange-correlation potential can be explicitly obtained. The accuracy of the results depends, clearly, on our selection of the orbit-generating functions W 6 c Cn and W e c Sn. 

The Kohn-Sham-Dirac equation (27) introduced in the last section is the basis of most relativistic electronic structure calculations in solid state theory. There are certain aspects which make the numerical solution of this four-component equation more involved than its non-relativistic coimterpart The Hamiltonian of the Kohn-Sham-Dirac equation is, unlike its Schrodinger equivalent and unlike the field-theoretical Hamiltonian (7) with the properly chosen normal order, not bounded below. In the limit of free, non-interacting particles the solutions of the Kohn-Sham-Dirac equation are plane waves with energies e(k) = cVk -I- c, where positive energies correspond to electrons and states with negative energy can be interpreted as positrons. For numerical procedures, which preferably use variational techniques to find electronic solutions, this property of the Dirac operator causes a severe problem, which can be circumvented by certain techniques like the application of a squared Dirac operator or a projection onto the properly chosen electronic states according to their above definition after Eq. (19). 

In the embedding formalism introduced by Wesolowski and Warshel [3], the total electron density is partitioned into two components. One of them is not optimized (frozen) and the other is subject to optimization. The optimized component is treated in a Kohn-Sham-like way, i.e., by means of a reference system of non-interacting electrons. The multiplicative potential in one-electron equations for embedded orbitals, Eq. (1) or Eqs. (20) and (21) of Ref. [3], differs from the Kohn-Sham... 

The generalization of density-functional theory that allows different orbitals for electrons with different spins is called spin-density-functional theory (Parr and Yang, Chapter 8). In spin-DFT, one deals separately with the electron density p (r) due to the spin-a electrons and the density p (r) of the spin-j8 electrons, and functionals such as E c become functionals of these two quantities E c = Exc[p" ,p ]. One deals with separate Kohn-Sham eigenvalue equations for the spin-a orbitals and the spin-/3 orbitals, where these equations have v c = 8Exc[p , p ] /8p° and similarly for For species like CH3... 

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. 

Thus, once we know the various contributions in equation (5-15) we have a grip on the potential Vs which we need to insert into the one-particle equations, which in turn determine the orbitals and hence the ground state density and the ground state energy by employing the energy expression (5-13). It should be noted that Veff already depends on the density (and thus on the orbitals) through the Coulomb term as shown in equation (5-13). Therefore, just like the Hartree-Fock equations (1-24), the Kohn-Sham one-electron equations (5-14) also have to be solved iteratively. 

While the hydrodynamical scheme mentioned above involves the density quantities directly, an alternative second scheme based on their orbital partitioning along the lines of the Kohn-Sham [4] version of time-independent DFT has been derived by Ghosh and Dhara [14]. In this scheme, one obtains the exact densities p(r, t) and j(r, t) from the TD orbitals (///,(r, t) obtained by solving the effective one-particle TD Schrodinger-like equations given by... 

---