Görling–Levy perturbation theory

Görling–Levy perturbation theory (GLPT) in Kohn–Sham (KS) density functional theory (DFT) is the analogue to what Møller–Plesset perturbation theory (MPPT) is in Hartree–Fock (HF) theory. Its basis is Rayleigh–Schrödinger perturbation theory (RSPT) and the adiabatic connection (AC). It describes electronic correlation effects. It is mostly used to second (GL2), rarely to third (GL3) or fourth (GL4) order, but becomes fast really increasingly computational expensive. It was published in 1993 and 1994 by Andreas Görling and Mel Levy.

Kohn–Sham correlation energy from Görling–Levy perturbation

The basis of GL perturbation theory is the adiabatic connection (AC) with the coupling constant $0\le \alpha \le 1$ connecting the artificial Kohn–Sham (KS) system of noninteracting electrons $\alpha =0$ to the real system of interacting electrons $\alpha =1$ with the AC Hamiltonian

${\hat{H}}_{\alpha }=\hat{T}+\alpha {\hat{V}}_{\text{ee}}+\sum _{i=1}^N{v}_{\alpha }(r_i)$

where $N$ is the number of electrons, $\hat{T}=-\frac{1}{2}{\sum }_i{\mathrm{\nabla }}_i^2$ the kinetic energy of the electrons, ${\hat{V}}_{\text{ee}}={\sum }_i{\sum }_{j>i}|r_i-r_j{|}^{-1}$ the electron-electron interaction. Görling and Levy expressed the coupling-strength dependent local multiplicative potential under the constraint, that the density $n(r)$ stays fixed along the AC as

$v_{\alpha }[n](r)=v_S[n](r)-\alpha v_{Hx}[n](r)-v_c^{\alpha }[n](r)$

where $v_S$ is the KS potential, $v_{Hx}$ the Hartree-exchange potential in first order, and the correlation potential for second order or higher $v_c^{\alpha }(r)={\alpha }^2{v}_c^{(2)}+{\alpha }^3{v}_c^{(3)}+{\alpha }^4{v}_c^{(4)}+\mathrm{...}$. As usual in perturbation theory we can express the correlation energy in a power series $E_c^{(0)}+\alpha E_c^{(1)}+{\alpha }^2{E}_c^{(2)}+\mathrm{...}$, where in GLPT the zeroth and first contribution vanish i.e. $E_c^{(0)}=E_c^{(1)}=0$. The second term is the Görling–Levy second order (GL2) correlation energy and can be evaluated with using the Slater–Condon rules and Brillouin's theorem in terms of occupied $i,j$ and unoccupied $a,b$ KS orbitals and eigenvalues Template:NumBlk where ${\mathrm{\Phi }}_S,{\mathrm{\Phi }}_k$ are ground state and excited KS determinants with their respective $E_0,E_k$ energies and $E_c^{\text{MP2}}$ is exactly the second order Møller–Plesset (MP2) correlation energy but evaluated with KS orbitals, $E_c^S$ the so called single excitation contribution to correlation which is missing in regular MPPT, but present in GLPT and ${\hat{v}}_x^{\text{NL}}$ is the nonlocal exchange operator from Hartree–Fock (HF) theory, ${\hat{v}}_x$ is the local Kohn–Sham (KS) exchange operator both evaluated with KS orbitals and lastly the notation $⟨ij||ab⟩=⟨ij|ab⟩-⟨ij|ba⟩$.

Hohenberg–Kohn functional from infinite Görling–Levy expansion

With GLPT up to infinite order one could in principle obtain the Hohenberg-Kohn (HK) functional exactly $F_{\text{HK}}[n]\equiv E[n]-\int drv(r)n(r)$ in terms of unoccupied and occupied KS orbitals $\{{\phi }_i\}$ and their eigenvalues $\{{\epsilon }_i\}$, where $E[n]$ is the electronic ground state energy and $v(r)$ the external potential. This is obviously only conceptually interesting since it is computational impossible. With the coupling constant expression Template:NumBlk

By setting $\alpha =1$ hence Template:NumBlk where in zeroth order $E_0=T_S=⟨{\mathrm{\Phi }}_S|\hat{T}|{\mathrm{\Phi }}_S⟩={\sum }_i\int dr{\varphi }_i^{\ast }(r)(-1/2{\mathrm{\nabla }}^2){\varphi }_i(r)$ is the KS kinetic energy with the KS potential $v_0(r)=v_S(r)$ and in first order $E_1=⟨{\mathrm{\Phi }}_S|{\hat{V}}_{\text{ee}}|{\mathrm{\Phi }}_S⟩=E_{Hx}[n]$ the Hartree-exchange (Hx) energy and its respective Hx potential $v_1(r)=v_{Hx}(r)$ and from second order the infinite GL$n$ correlation (c) energy with $n\to \mathrm{\infty }$, which is the exact Kohn–Sham (KS) correlation energy $\underset{n\to \mathrm{\infty }}{lim}{\sum }_{j=2}^n{E}_c^{\text{GLn}}[n]={\sum }_{j=2}^{\mathrm{\infty }}{E}_j=E_c[n]$ and the corresponding correlation potential $v_c(r)={\sum }_{j=2}^{\mathrm{\infty }}{v}_j(r)$. Similarly, if one would do Møller–Plesset perturbation theory up to infinite order one would obtain the exact Hartree–Fock (HF) correlation energy $\underset{n\to \mathrm{\infty }}{lim}{\sum }_{j=2}^n{E}_c^{\text{MPn}}=E_c^{\text{HF}}[\{{\mathrm{\Phi }}^{\text{HF}},{\mathrm{\Phi }}_i^a,{\mathrm{\Phi }}_{ij}^{ab},{\mathrm{\Phi }}_{ijk}^{abc},\mathrm{...}\}]$ where $i,j,k$ denote occupied and $a,b,c$ unoccupied HF orbitals and their respective singly, doubly, triply and so on excited Slater determinants. In this notation ${\mathrm{\Phi }}^{\text{HF}}$ is the HF determinant and ${\mathrm{\Phi }}_S$ the KS determinant.

Optimized Effective Potential (OEP) method

In the later half of their article Görling and Levy connect their perturbation theory to the optimized effective potential method.