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.
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)
where is the i-th KohnâÂÂSham orbital. The summation is performed over all the occupied KohnâÂÂSham orbitals.
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
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.
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).
A conceptually really important quantity in OFDFT is the Pauli kinetic energy. As the Kohn-Sham correlation energy links the real system of interacting electrons to the artificial Kohn-Sham (KS) system of noninteracting electrons, the Pauli kinetic energy links the KS system to the fictitious system noninteracting model bosons. In the same way as the KS interacting energy it is highly KS-orbital dependent and must be in practice approximated.
The term Pauli comes from the fact, that the functional is related to the Pauli exclusion principle. for an electron number of .
The exchange energy in orbital-free density functional theory (OFDFT) is the Dirac exchange as a Local Density Approximation (LDA) (1930) It is related to the homogeneous electron gas (HEG).
An important bound for the LDA exchange energy is the Lieb-Oxford inequality.
State-of-the-art kinetic energy density functionals for orbital-free density functional theory, and still under active research, are the so-called nonlocal (NL) kinetic energy density functionals. They admit the general form
where is a nonlocal kernel depending on the electron density, and is some constant. The kernel is often determined based on the linear response function of the uniform electron gas, which makes the KEDF highly accurate for simple metals. Efforts have also been made to extend the KEDF to semiconductors by designing the kernel based on the linear response function in these systems.
The analogue to the Kohn-Sham (KS) equations (1965) in Orbital-free density functional theory (OFDFT) is the Levy-Perdew-Sahni (LPS) equation (1984), an effectively bosonic Schrödinger equation
where is the Kohn-Sham (KS) potential, the Pauli potential, the highest occupied KS orbital and the square root of the density. One big benefit of the LPS equation being so intimately related to the KS equations is that an existing KS code can be easily modified in an OF code with ejecting all orbitals except for one in the Self-Consistent-Field (SCF) cycle.
Starting from Euler-Lagrange equation of density functional theory , simultaneously adding and subtracting the von Weizsäcker potential, i.e. the functional derivative of the vW functional and acknowleding the definition of the Pauli kinetic energy, while the functional derivative of the Pauli kinetic energy with respect to the density is the Pauli potential . Expanding the functional derivative via chain rule and as a last step multiplying both sides by the square root of the density yields the LPS equation.
With the linear transformation and by defining the bosonic potential as the LPS equation evolves to the bosonic Schrödinger equation
.
Note that the normalization constraint in BraâÂÂket notation holds, since .
Quite recently also a time-dependent version of OFDFT has been developed. It is also implemented in DFTpy.
A free open-source software package for OFDFT DFTpy has been developed by the Pavanello Group. It was launched in 2020. The most recent version is 2.1.1.