next up previous contents index
Next: Implementation Up: Orbital Dependent Kohn-Sham Density Previous: Orbital Dependent Kohn-Sham Density   Contents   Index

Theoretical Background

Approximations to the exchange-correlation (XC) functional of the Kohn-Sham (KS) Density Functional Theory (DFT) can be classified by the so-called “Jacob’s ladder.” The ground on which the ladder lies is the Hartree approximation (XC energy is zero), and the first rung is the local density approximation (LDA) in which the XC energy density is a simple local function of the density. The second rung of the Jacob’s ladder is the generalized gradient approximation (GGA): in this case the XC energy density depends also on the gradient of the density. In the third rung (meta-GGA) an additional variable is used, the Kohn-Sham kinetic energy density which allows, e.g., to construct self-correlation-free functionals. Functionals in the above rungs can have high accuracy for different class of problems in chemistry and solid-state physics, but their main limitation is the self-interaction error (SIE)[156,157,158,159]. To avoid the SIE the exchange must be treated exactly and this can be achieved by functionals in the fourth rung which depend explicitly on all the occupied KS orbitals. In the KS formalism the EXX (exact-exchange) energy is (for closed-shell systems ns = 2 )[156,157,158,159]:

ExEXX = - $\displaystyle {\frac{{n_s}}{{2}}}$$\displaystyle \sum_{a}^{{occ.}}$$\displaystyle \sum_{b}^{{occ.}}$$\displaystyle \int$$\displaystyle \int$d$\displaystyle \bf r$d$\displaystyle \bf r{^\prime}$$\displaystyle {\frac{{\phi^{\text{KS}}_a({\bf r}) \phi^{\text{KS}}_b({\bf r}) \...
...a({\bf r}') \phi^{\text{KS}}_b({\bf r}') }}{{\left\Vert{\bf r-r'}\right\Vert}}}$  . (16.1)

i.e. the same functional form of the Hartree-Fock (HF) exchange but computed with KS orbitals which are obtained using a self-consistent local EXX potential. At this point we should recall that hybrid DFT functionals (including HF exchange), doesn't belong to the KS formalism: in hybrid DFT, in fact, the non-local HF exchange operator $ \hat{{v}}^{\text}_{}$NLx($ \bf r,r'$) = - $ \sum_{a}^{{occ.}}$$ {\frac{{ \phi_a({\bf r})\phi_a({\bf r}')}}{{\left\Vert{\bf r-r'}\right\Vert}}}$ is employed in the self-consistent (Generalized Kohn-Sham) equations determining the orbitals.

While LDA, GGA, meta-GGA and hybrid functionals are implemented (for ground-state calculations) in the dscf and ridft, the odft module considers functionals of the fourth rung. Currently exchange-only orbital-dependent approaches are implemented in the odft module. The EXX KS local potential ( vxEXX($ \bf r$) ) can be obtained using the optimized effective potential (OEP) method (in each self-consistent step):[157,160,158,159]:

$\displaystyle \int$d$\displaystyle \bf r{^\prime}$χs($\displaystyle \bf r,r'$)vxEXX($\displaystyle \bf r$) = $\displaystyle \sum_{a}^{{occ.}}$$\displaystyle \sum_{s}^{{vir.}}$2ns$\displaystyle \left\langle\vphantom{{\phi_a}}\right.$φa$\displaystyle \left.\vphantom{{\phi_a}}\right\vert$$\displaystyle \hat{{v}}^{\text}_{}$NLx$\displaystyle \left\vert\vphantom{{\phi_s}}\right.$φs$\displaystyle \left.\vphantom{{\phi_s}}\right\rangle$$\displaystyle {\frac{{ \phi_s({\bf r})\phi_a({\bf r})}}{{\epsilon_a-\epsilon_s}}}$ (16.2)

where χs($ \bf r,r'$) = $ \sum_{a}^{{occ.}}$2ns$ \sum_{s}^{{vir.}}$$ {\frac{{ \phi_a({\bf r})\phi_s({\bf r})\phi_s({\bf r}')\phi_a({\bf r}')}}{{\epsilon_a-\epsilon_s}}}$ is the non-interacting density response.

An effective approximation to the OEP-EXX potential is given by the Localized Hartree-Fock (LHF) potential[156] which is given by

vxLHF($\displaystyle \bf r$) = - $\displaystyle \sum_{{ij}}^{{occ.}}$ns$\displaystyle {\frac{{ \phi_i({\bf r})\phi_j({\bf r})}}{{ \rho({\bf r})}}}$$\displaystyle \int$d$\displaystyle \bf r{^\prime}$$\displaystyle {\frac{{\phi_i({\bf r}')\phi_j({\bf r}')}}{{\left\Vert{\bf r-r'}\right\Vert }}}$  
  + $\displaystyle \sum_{{ij}}^{{occ.}}$ns$\displaystyle {\frac{{ \phi_i({\bf r})\phi_j({\bf r})}}{{ \rho({\bf r})}}}$$\displaystyle \left\langle\vphantom{{\phi_i}}\right.$φi$\displaystyle \left.\vphantom{{\phi_i}}\right\vert$vxLHF-$\displaystyle \hat{{v}}^{\text}_{}$NLx$\displaystyle \left\vert\vphantom{{\phi_j}}\right.$φj$\displaystyle \left.\vphantom{{\phi_j}}\right\rangle$ (16.3)

where the first term is called Slater potential and the second term correction term. If terms ij are neglected in the correction term, the Krieger-Li-Iafrate (KLI) potential[161] is obtained. Note that the Eq. (16.3) depends only on occupied orbitals, whereas Eq. (16.2) depends also on virtual orbitals. The LHF total energy is assumed to be the EXX total energy, even if LHF is not variational (although the deviation from the EXX energy is very small, usually below 0.01% ). The LHF potential is equivalent to the Common Energy Denominator Approximation (CEDA) [162] and to the Effective Local Potential (ELP) [163].

Both OEP-EXX and LHF (in contrast to functionals of the first three rungs) satisfy the HOMO condition[161]

$\displaystyle \left\langle\vphantom{{\phi_{HOMO}}}\right.$φHOMO$\displaystyle \left.\vphantom{{\phi_{HOMO}}}\right\vert$vx$\displaystyle \left\vert\vphantom{{\phi_{HOMO}}}\right.$φHOMO$\displaystyle \left.\vphantom{{\phi_{HOMO}}}\right\rangle$ = $\displaystyle \left\langle\vphantom{{\phi_{HOMO}}}\right.$φHOMO$\displaystyle \left.\vphantom{{\phi_{HOMO}}}\right\vert$$\displaystyle \hat{{v}}^{\text}_{}$NLx$\displaystyle \left\vert\vphantom{{\phi_{HOMO}}}\right.$φHOMO$\displaystyle \left.\vphantom{{\phi_{HOMO}}}\right\rangle$  , (16.4)

and the asymptotic relation[164,165]

vx($\displaystyle \bf r_{l}^{}$)$\displaystyle \;\stackrel{{r_l\rightarrow\infty}}{{\longrightarrow}}\;$$\displaystyle \left\langle\vphantom{{\phi_M}}\right.$φM$\displaystyle \left.\vphantom{{\phi_M}}\right\vert$vx-$\displaystyle \hat{{v}}^{\text}_{}$NLx$\displaystyle \left\vert\vphantom{{\phi_M}}\right.$φM$\displaystyle \left.\vphantom{{\phi_M}}\right\rangle$ - $\displaystyle {\frac{{1}}{{r_l}}}$ . (16.5)

where φM is the highest occupied orbital which do not have a nodal surface in the asymptotic region along direction $ \bf r_{l}^{}$ . Considering together with condition (16.4), we finally obtain that:
vx($ \bf r$) will approach -1/r along all directions where φHOMO($ \bf r$) does not have a nodal surface in the asymptotic region (e.g. this is the case of atoms);
on directions which belong to the nodal surface of the HOMO, the vx($ \bf r$) will approach $ \left\langle\vphantom{{\phi_M}}\right.$φM$ \left.\vphantom{{\phi_M}}\right\vert$vx-$ \hat{{v}}^{\text}_{}$NLx$ \left\vert\vphantom{{\phi_M}}\right.$φM$ \left.\vphantom{{\phi_M}}\right\rangle$ - 1/r .

Both OEP-EXX and LHF gives total energies very close to the Hartree-Fock one (actually ELHF > EEXX > EHF ), thus, without an appropriate correlation functional, these methods are not suitable for thermochemistry. On the other hand OEP-EXX and LHF give very good KS orbital spectra. In fact the eigenvalues of the HOMO is very close to the Hartree-Fock and to exact ionization potential (I.P): this is in contrast to functional of the first three rungs which underestimate the HOMO energy by several eVs. In addition a continuum set of bound unoccupied orbitals are obtained. Thus OEP-EXX or LHF KS orbitals are very good input quantities for computing NMR shielding constants[166], energy-levels in hybrid interfaces[167] and TD-DFT excitation energies [168] (the latter using LDA/GGA kernels, not the hybrid ones).

next up previous contents index
Next: Implementation Up: Orbital Dependent Kohn-Sham Density Previous: Orbital Dependent Kohn-Sham Density   Contents   Index