In computational chemistry, orbital-free density functional theory (OFDFT) is a quantum mechanical approach to electronic structure determination which is based on functionals of the electronic density. It is most closely related to the Thomas–Fermi model. Orbital-free density functional theory is, at present, less accurate than Kohn–Sham density functional theory models, but it has the advantage of being fast, so that it can be applied to large systems.
Kinetic energy of electrons: an orbital-dependent functional The Hohenberg–Kohn theorems guarantee that, for a system of atoms, there exists a functional of the electron density that yields the total energy. Minimization of this functional with respect to the density gives the ground-state density from which all of the system's properties can be obtained. Although the Hohenberg–Kohn theorems tell us that such a functional exists, they do not give us guidance on how to find it. In practice, the density functional is known exactly except for two terms. These are the electronic kinetic energy and the exchange–correlation energy. The lack of the true exchange–correlation functional is a well known problem in DFT, and there exists a huge variety of approaches to approximate this crucial component. In general, there is no known form for the interacting kinetic energy in terms of electron density. In practice, instead of deriving approximations for interacting kinetic energy, much effort was devoted to deriving approximations for non-interacting (Kohn–Sham) kinetic energy, which is defined as (in atomic units)
T S [ { ϕ i } ] = ∑ i = 1 N ⟨ ϕ i | − 1 2 ∇ 2 | ϕ i ⟩ , {\displaystyle T_{S}[\{\phi _{i}\}]=\sum _{i=1}^{N}\langle \phi _{i}|-{\frac {1}{2}}\nabla ^{2}|\phi _{i}\rangle ,}
where | ϕ i ⟩ {\displaystyle |\phi _{i}\rangle } is the i-th Kohn–Sham orbital. The summation is performed over all the occupied Kohn–Sham orbitals.
Thomas-Fermi (TF) kinetic energy One of the first attempts to do this (even before the formulation of the Hohenberg–Kohn theorem) was the Thomas–Fermi model (1927), which wrote the kinetic energy as
T TF [ n ] = 3 10 ( 3 π 2 ) 2 3 ⏟ C T F ∫ [ n ( r ) ] 5 3 d 3 r . {\displaystyle T_{\text{TF}}[n]=\underbrace {{\frac {3}{10}}(3\pi ^{2})^{\frac {2}{3}}} _{C_{TF}}\int [n(\mathbf {r} )]^{\frac {5}{3}}\,d^{3}r.}
This expression is based on the homogeneous electron gas (HEG) and a Local Density Approximation (LDA), thus, is not very accurate for most physical systems. By formulating Kohn–Sham kinetic energy in terms of electron density, one avoids diagonalizing the Kohn–Sham Hamiltonian for solving for the Kohn–Sham orbitals, therefore saving the computational cost. Since no Kohn–Sham orbital is involved in orbital-free density functional theory, one only needs to minimize the system's energy with respect to the electron density. An important bound for the TF kinetic energy is the Lieb-Thirring inequality.
Von Weizsäcker (vW) kinetic energy A notable historical improvement of the Thomas-Fermi model is the von Weizsäcker (vW) kinetic energy (1935), which is exactly the kinetic energy for noninteracting bosons and can be regarded as a Generalized Gradient approximation (GGA).
T vW [ n ] = 1 8 ∫ ∇ n ( r ) ⋅ ∇ n ( r ) n ( r ) d 3 r = ∫ n ( r ) ( − 1 2 Δ ) n ( r ) d 3 r {\displaystyle T_{\text{vW}}[n]={\frac {1}{8}}\int {\frac {\nabla n(\mathbf {r} )\cdot \nabla n(\mathbf {r} )}{n(\mathbf {r} )}}d^{3}r=\int {\sqrt {n(\mathbf {r} )}}(-{\frac {1}{2}}\Delta ){\sqrt {n(\mathbf {r} )}}d^{3}r}
… excerpt ends here. Continue reading the full article.
