What is Density functional theory?

A theory that uses the fact that the ground-state energy is determined by the electron density, solving the Kohn–Sham equations to obtain the electronic states and energy of a material. The core of first-principles calculation.

Japanese name: 密度汎関数理論

Density functional theory appears in 7 practice questions across 1 topic, including Computational materials science.

Field: Computational materials science

The idea behind the theory

The properties of materials are determined by the states of their electrons, but the wave function of many mutually interacting electrons has too many variables to be solved as it is. Hohenberg and Kohn showed that the ground-state energy is determined by the electron density alone, a function of only three dimensions. Kohn and Sham replaced this with a problem in which fictitious non-interacting electrons move in an effective potential, arranged so that their electron density is the same as that of the real system. The complicated interactions between electrons are gathered into a term called the exchange-correlation energy.

Quantities obtained by calculation

DFT takes the kinds and arrangement of atoms as input and gives the total energy, the forces on the atoms, the electronic band structure and so on. Moving the atoms until the forces are zero gives the stable structure and lattice parameter; comparing energies of different structures gives phase stability; supercell calculations give the formation energies of defects such as vacancies; and calculations with applied strain give elastic constants. Because it uses no experimental data, it is also called first-principles calculation, and it is used to predict materials not yet made and to generate data for the CALPHAD method and for interatomic potentials.

Approximations and limits

The exact form of the exchange-correlation energy is unknown, so approximations such as the local density approximation (LDA) and generalized gradient approximation (GGA) are used. These have limits: band gaps come out smaller than experiment, van der Waals forces are described poorly, and the properties of oxides with localized d or f electrons are sometimes wrong; hybrid functionals and correction methods are used to address them. Ordinary calculations are also of the ground state at 0 K and handle a few hundred atoms. For finite-temperature properties, phonon calculations or molecular dynamics are combined.

Try it: questions about Density functional theory

Which statement about the ideas behind density functional theory (DFT) is correct?

  1. ABased on the electron density determining the ground-state energy, it solves the Kohn–Sham equations and treats electron interactions with an approximate exchange-correlation functional
  2. BIt fits the form of an interatomic potential to experimental data and solves Newton's equations of motion for the atoms, so the electronic states never appear anywhere in the calculation
  3. CIt obtains the wave functions of all electrons exactly, without any exchange-correlation approximation, so it predicts any property of any material with the same accuracy as experiment does
Show answer and explanation

Answer: A. Based on the electron density determining the ground-state energy, it solves the Kohn–Sham equations and treats electron interactions with an approximate exchange-correlation functional

The wave function of many electrons has too many variables to be solved as it is. DFT takes the electron density as its main variable, based on the Hohenberg–Kohn theorem that the ground-state energy is determined by the electron density, a function of only three dimensions. In practice, the Kohn–Sham equations, in which fictitious non-interacting electrons move in an effective potential, are solved repeatedly until the electron density no longer changes (a self-consistent calculation). The complicated interactions between electrons are gathered into the exchange-correlation functional, for which approximations such as LDA and GGA are used, so the results are not exact and their accuracy differs with the material and property. Solving equations of motion with an empirical interatomic potential is classical molecular dynamics.

Which is a correct tendency of the results of standard DFT calculations using LDA or GGA exchange-correlation functionals?

  1. AThey tend to overestimate the band gaps of semiconductors and insulators compared with experiment; for silicon GGA gives about 2 eV
  2. BThey tend to underestimate the band gaps of semiconductors and insulators compared with experiment; for silicon GGA gives about 0.6 eV
  3. CBesides the ground state at 0 K, they give finite-temperature free energies and melting points directly from a single structural optimization
Show answer and explanation

Answer: B. They tend to underestimate the band gaps of semiconductors and insulators compared with experiment; for silicon GGA gives about 0.6 eV

It is well known that with LDA or GGA, a band gap read from the difference of eigenvalues of the Kohn–Sham equations comes out smaller than experiment. For silicon, GGA gives about 0.6 eV, roughly half the experimental 1.12 eV, so the statement that they overestimate gets the trend backwards. When accurate values are needed, hybrid functionals or the GW approximation are used. On the other hand, quantities tied to the total energy of the ground state, such as lattice parameters, elastic moduli and formation energies, are obtained with good accuracy for many materials. Ordinary DFT is a calculation of the ground state at 0 K, so obtaining finite-temperature properties requires combining it with calculations of lattice vibrations (phonons) or molecular dynamics. Standard functionals also describe van der Waals forces poorly, so corrections are added for layered materials and adsorption.

FAQ

Are first-principles calculation and DFT the same thing?

First-principles calculation is a general term for methods that calculate from the basic laws of quantum mechanics without empirical parameters, and in materials science calculations by DFT are by far the most widely used. The two are therefore often used almost interchangeably, but more rigorous methods of quantum chemistry are also included in first-principles calculation.

Why do band gaps come out too small?

Because the LDA and GGA exchange-correlation energies retain an error in which an electron interacts with itself, and they do not correctly describe how the energy changes with the number of electrons. Another reason is that the eigenvalues of the Kohn–Sham equations do not, in principle, represent excitation energies. The GW approximation and hybrid functionals improve this.

How large a system can be calculated?

Ordinary calculations handle tens to hundreds of atoms, and the cost grows roughly as the cube of the number of atoms. Because periodic boundary conditions are used, a crystal can be represented by its unit cell or a supercell of a few hundred atoms, but problems needing large regions, such as dislocations and grain boundaries, are handled by combining DFT with molecular dynamics using interatomic potentials and similar methods.

Density functional theory in the explanations

Questions that mention Density functional theory (7)

  1. Which statement about the ideas behind density functional theory (DFT) is correct?
  2. Which is a correct tendency of the results of standard DFT calculations using LDA or GGA exchange-correlation functionals?
  3. A 108-atom supercell (with periodic boundary conditions) of a metal calculated with DFT has a total energy of −405.00 eV. Removing one atom to create one vacancy and relaxing the structure gives −400.10 eV. Which is closest to the vacancy formation energy calculated as E_v = E(N−1) − [(N−1)/N]E(N)?
  4. Which pairing of a materials computation method with the range of space and time it can usually handle is most appropriate?
  5. Which statement about the interatomic potentials used in molecular dynamics is correct?
  6. Which statement about the kinetic Monte Carlo (kMC) method is correct?
  7. Which statement about the CALPHAD method is correct?

Related terms

Molecular dynamics Interatomic potential Vacancy Phonon Phase-field method Finite element method Monte Carlo method CALPHAD method Gibbs energy Materials informatics

Back to the glossary