Density functional theory

Density functional theory is a computational quantum mechanical modelling method used to investigate the electronic structure of many-body systems, particularly atoms, molecules, and the condensed phases. It stands as one of the most popular and versatile techniques in computational physics, computational chemistry, and materials science. The theory reformulates the intractable many-body problem of interacting electrons in terms of the electron density, a simpler quantity, rather than the complex many-electron wavefunction. This approach provides a powerful framework for calculating ground-state properties of systems ranging from small molecules to extended solids and complex materials.

Theoretical foundations

The rigorous foundation was established by Pierre Hohenberg and Walter Kohn in 1964 with the Hohenberg–Kohn theorem. This theorem proves that the ground-state electron density uniquely determines all properties of a quantum system, including the wave function and total energy. A subsequent practical scheme was introduced by Walter Kohn and Lu Jeu Sham in 1965, leading to the Kohn–Sham equations. These equations map the complex interacting system onto a fictitious system of non-interacting electrons moving in an effective potential, known as the Kohn–Sham potential. This potential includes contributions from the electron-nucleus attraction, the classical electron-electron repulsion, and a critical term called the exchange-correlation functional, which encapsulates all quantum mechanical many-body effects.

Formulation and methodology

In practice, solving the Kohn–Sham equations involves a self-consistent iterative procedure. An initial guess for the electron density is used to construct the effective potential, from which the Kohn–Sham orbitals are computed by solving a set of one-electron Schrödinger-like equations. The new density calculated from these orbitals is then used to update the potential, and the cycle repeats until convergence is achieved. This process is computationally implemented within various software packages like VASP, Quantum ESPRESSO, and Gaussian (software). The methodology requires a basis set to represent the orbitals, with common choices including plane waves for periodic systems and Gaussian-type orbitals for molecular calculations, often used in conjunction with pseudopotentials to treat core electrons efficiently.

Approximations and functionals

The accuracy of calculations is entirely dependent on the approximation chosen for the unknown exchange-correlation functional. The simplest is the Local-density approximation, which assumes the functional depends only on the local value of the density. A major advancement was the introduction of the Generalized gradient approximation, which also includes the density gradient, improving accuracy for molecules and surfaces. More sophisticated meta-GGAs and hybrid functionals, such as B3LYP and PBE0, incorporate a portion of exact exchange from Hartree–Fock theory. The development of new functionals is an active area of research pursued by groups like those at the University of California, Irvine and the Max Planck Institute.

Applications in chemistry and physics

The method is extensively applied across scientific disciplines. In chemistry, it is used to calculate molecular structures, vibrational spectra, reaction mechanisms, and activation energies, aiding in the design of new catalysts and pharmaceuticals. In condensed matter physics and materials science, it is indispensable for predicting properties of semiconductors, superconductors, magnets, and nanomaterials. It has been crucial in studies of graphene, perovskites, and high-temperature superconductivity. Research institutions like MIT, Stanford University, and national labs such as Lawrence Berkeley National Laboratory routinely employ it for materials discovery.

Limitations and challenges

Despite its success, the theory has well-known limitations. The most significant is the approximate nature of exchange-correlation functionals, which can lead to errors in predicting band gaps, dissociation energies, and properties of strongly correlated systems like those containing transition metals or lanthanides. Describing van der Waals interactions accurately requires specialized dispersion-corrected functionals. Furthermore, the standard formalism is a ground-state theory, making the treatment of excited states, such as those probed in photoelectron spectroscopy, more challenging, though extensions like time-dependent density functional theory address this. Ongoing work at institutions like the Fritz Haber Institute focuses on developing more accurate functionals and methods for dynamical and non-equilibrium processes.

Category:Computational chemistry Category:Condensed matter physics Category:Quantum chemistry