Reference state Kohn-Sham

The Kohn-Sham equation becomes a reference system for DFT because it yields the correct ground state via ... 

Further pragmatic moves are described in details in numerous books and reviews of which we cite the most concise and recent Ref. [82], Two further hypotheses are an important complement to the above cited theorems. One is the locality hypothesis, another is the Kohn-Sham representation of the single determinant reference state in terms of orbitals. The locality has been seriously questioned by Nesbet in recent papers [83,84], however, it remains the only practically implemented solution for the DFT. The single determinant form of the reference state in its turn guarantees that all the averages of the electron-electron interaction appearing in this context are in fact calculated with the two-electron density given by the determinant term in Eq. (5) with no cumulant. 

The KS equations are obtained by differentiating the energy with respect to the KS molecular orbitals, analogously to the derivation of the Hartree-Fock equations, where differentiation is with respect to wavefunction molecular orbitals (Section 5.2.3.4). We use the fact that the electron density distribution of the reference system, which is by decree exactly the same as that of the ground state of our real system (see the definition at the beginning of the discussion of the Kohn-Sham energy), is given by (reference [9])... 

In Kohn-Sham DFT based approaches, expressions that are of similar structure as Eqs. (9a) and (9b) are obtained, but in the form of contributions from all occupied Kohn-Sham MOs The excited-state wavefunctions are at the same time formally replaced by the unoccupied MOs, and the many-electron perturbation operators /T(M41, etc. by their one-electron counterparts //(M-41, etc. Orbital energies e and ea formally substitute the total energies of the states (see later). Thus, similar interpretations of NMR parameters can be worked out in which the highest occupied MO-lowest unoccupied MO gap (HLG) plays a highly important role. It must be emphasized, though, that there is no one-to-one correspondence between the excited states of the SOS equations and the unoccupied orbitals which enter the DFT expressions, nor between excitation energies and orbital energy differences, i.e., there are no one-determinantal wavefunctions in Kohn-Sham DFT perturbation theory which approximate the reference and excited states. 

If the Kohn-Sham orbitals [52] of density functional theory (DFT) [53] are used instead of Hartree-Fock orbitals in the reference state [54], the RI can become essential for the realization of electron propagator calculations. Modern implementations of Kohn-Sham DFT [55] use the variational approximation of the Coulomb potential [45,46] (which is mathematically equivalent to the RI as presented above), and four-index integrals are not used at all. A very interesting example of this combination is the use of the GW approximation [56] for molecular systems [54],... 

Y. Shigeta et al., Electron propagator calculations with Kohn-Sham reference states. Int. J. Quantum. Chem. 85, 411 120 (2001)... 

Assuming the existence of such a Frechet functional derivative [26,102] constitutes the locality hypothesis. If this hypothesis were valid, the OEL and Kohn-Sham equations would be equivalent, determining the same model or reference state. 

In the present chapter, we will focus on the simulation of the dynamics of photoexcited nucleobases, in particular on the investigation of radiationless decay dynamics and the determination of associated characteristic time constants. We use a nonadiabatic extension of ab initio molecular dynamics (AIMD) [15, 18, 21, 22] which is formulated entirely within the framework of density functional theory. This approach couples the restricted open-shell Kohn-Sham (ROKS) [26-28] first singlet excited state, Su to the Kohn-Sham ground state, S0, by means of the surface hopping method [15, 18, 94-97], The current implementation employs a plane-wave basis set in combination with periodic boundary conditions and is therefore ideally suited to condensed phase applications. Hence, in addition to gas phase reference simulations, we will also present nonadiabatic AIMD (na-AIMD) simulations of nucleobases and base pairs in aqueous solution. 

For variational methods, such as Hartree-Fock (HF), multi-configurational self-consistent field (MCSCF), and Kohn-Sham density functional theory (KS-DFT), the initial values of the parameters are equal to zero and 0) thus corresponds to the reference state in the absence of the perturbation. The A operators are the non-redundant state-transfer or orbital-transfer operators, and carries no time-dependence (the sole time-dependence lies in the complex A parameters). Furthermore, the operator A (t)A is anti-Hermitian, and tlie exponential operator is thus explicitly unitary so that the norm of the reference state is preserved. Perturbation theory is invoked in order to solve for the time-dependence of the parameters, and we expand the parameters in orders of the perturbation... 

Hyperfine couplings, in particular the isotropic part which measures the spin density at the nuclei, puts special demands on spin-restricted wave-functions. For example, complete active space (CAS) approaches are designed for a correlated treatment of the valence orbitals, while the core orbitals are doubly occupied. This leaves little flexibility in the wave function for calculating properties of this kind that depend on the spin polarization near the nucleus. This is equally true for self-consistent field methods, like restricted open-shell Hartree-Fock (ROHF) or Kohn-Sham (ROKS) methods. On the other hand, unrestricted methods introduce spin contamination in the reference (ground) state resulting in overestimation of the spin-polarization. 

---