First principles electronic structure: density functional theory - PowerPoint PPT Presentation

First principles electronic structure density functional theory
1 / 24

  • Uploaded on
  • Presentation posted in: General

First principles electronic structure: density functional theory. Electronic Schroedinger equation. Atomic units, classical nuclei, Born-Oppenheimer approximation. Describes 99% of condensed matter, materials physics, chemistry, biology, psychology, sociology, ....

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

Download Presentation

First principles electronic structure: density functional theory

An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -

Presentation Transcript

First principles electronic structure density functional theory

First principles electronic structure: density functional theory

Electronic schroedinger equation

Electronic Schroedinger equation

Atomic units, classical nuclei, Born-Oppenheimer approximation.

Describes 99% of condensed matter, materials physics, chemistry, biology, psychology, sociology, ....

3 basic parameters: electronic mass, proton mass, electric charge

First principles electronic structure density functional theory

  • Prediction: in principle, just need to know what atoms are, can predict all properties of quantum system to very high accuracy

No parameters, no ‘garbage in, garbage out’

Solving the schroedinger equation




Solving the Schroedinger Equation

  • For example, consider that we have n electrons populating a 3D space. Let’s divide 3D space into NxNxN=2x2x2 grid points. To reconstruct Ψ, how many points must we keep track of?

Two solution philosophies approximate

Two solution philosophies (approximate)

Density functional theory

Density functional theory

  • Hohenberg and Kohn (Phys. Rev. 1964) introduced two theorems:

    • Electron density <-> the external potential + constant

    • Total energy of any density is an upper bound to the ground state energy, if we know the functional

    • But mapping is unknown

Very clear and detailed proofs of these theorems can be found in “Electronic Structure,” Richard M. Martin (Cambridge University Press 2004 ,ISBN 0 521 78285 6).

First theorem

First theorem

  • For each external potential, there is a unique ground state electronic density

  • proof outline:

    • Suppose there are two different potentials

    • Suppose they have the same ground state density

    • Show that this leads to an inconsistency

Hohenberg kohn i

Hohenberg-Kohn I

- The external potential corresponds to a unique ground state electron density.

- A given ground state electron density corresponds to a unique external potential

- In particular, there is a one to one correspondence between the external potential and the ground state electron density

Second theorem

Second theorem

  • Variational principle with respect to the density

  • instead of finding an unknown function of 1023 variables, we need only find an unknown function of 3, the charge density.

  • incredible simplification:

  • Still need to find the functional!

What might f r look like

What might F[ρ(r)] look like?

From simple inspection:


Naively, we might expect the functional to contain terms that resemble the kinetic energy of the electrons and the coulomb interaction of the electrons

Ingredients of density functional theory


Ingredients of Density Functional Theory

  • Ingredients:

    • Note that what differs from one electronic system to another is the external potential of the ions

    • Hohenberg-Kohn I: one to one correpondence between the external potential and a ground state density

    • Hohenberg-Kohn II: Existence of a universal functional such that the ground state energy is minimized at the true ground state density

      • The universality is important. This functional is exactly the same for any electron problem. If I evaluate F for a given trial orbital, it will always be the same for that orbital - regardless of the system of particles.

    • Kohn-Sham: a way to approximate the functional F

Euler lagrange system

Euler-Lagrange System

The Hohenberg-Kohn theorems give us a variational statement about the ground state density:

“the exact density makes the functional derivative of F exactly equal to the negative of the external potential (to within a constant)”

If we knew how to evaluate F, we could solve all Coulombic problems exactly.

However, we do not know how to do this. We must, instead, approximate this functional. This is where Kohn-Sham comes in.

Kohn sham approach

Kohn-Sham Approach

Kohn and Sham said:

Separate kinetic energy coulomb energy, and other

Kinetic energy of the system of non-interacting electrons at the same density.

Coulomb is the electrostatic term (Hartree)

Exchange-correlation is everything else

Kohn sham approach1

Kohn-Sham Approach

The next step to solving for the energy is to introduce a set of one-electron orthonormal orbitals.

Now the variational condition can be applied, and one obtains the one-electron Kohn-Sham equations.

Where VXC is the exchange correlation functional, related to the xc energy as:

The exchange correlation functional

The Exchange-Correlation Functional

The exchange-correlation functional is clearly the key to success of DFT.

One of the great appealing aspects of DFT is that even relatively simple approximations to VXC can give quite accurate results.

The local density approximation (LDA) is by far the simplest and used to be the most widely used functional.

Approximate as the xc energy of homogeneous electron gas

Ceperley and Alder* performed accurate quantum Monte Carlo calculations for the electron gas

Fit energies(e.g., Perdew-Zunger) to get f(p)

Solving the kohn sham system

Solving the Kohn-Sham System

To solve the Kohn-Sham equations, a number of different methods exist.

Basis set expansion of the orbitals

-Localized orbitals : molecules, etc

-Plane waves: solids, metals, liquids

Meaning of the orbitals:

-Kohn-Sham: In principle meaningless, only representation of the density

-Hartree-Fock: Electron distributions in the non-interacting approximation

There is some ad-hoc justification for using K-S orbitals as approximate quasiparticle distributions, but it’s qualitative at best (see Fermi liquid theory).

Self consistent approach

self-consistent approach

  • like Hartree, LDA-DFT equations must be solved self-consistently

  • great effort to develop efficient and scalable algorithms

    • remarkably successful

    • widely available

      • can download computer codes that perform these calculations

Local density approximation

  • For example, here are various energy contributions for a Mn atom:

  • Hartree (ECV, EVV)

  • Kinetic (T0,V)

  • Exchange (EX)

  • Correlation energy is about EC ~ 0.1EX

Local Density Approximation

In the original Kohn-Sham paper, the authors themselves cast doubt on its accuracy for many properties. “We do not expect an accurate description of chemical bonding.”

And yet, not until at least 10 years later (the 70’s), time and time again it was shown that LDA provided remarkably accurate results.

LDA was shown to give excellent agreement with experiment for, e.g., lattice constants, bulk moduli, vibrational spectra, structure factors, and much more.

One of the reasons for its huge success is that, in the end, only a very small part of the energy is approximated.

Local density approximation1

Local Density Approximation

LDA also works well because errors in the approximation of exchange and correlation tend to cancel.

For example, in a typical LDA atom, there’s a ~10% underestimate in the exchange energy.

This error in exchange is compensated by a ~100-200% overestimate of the correlation energy.

Local density approximation2

Slow varying

Faster varying

Local Density Approximation

In theory, the LDA method should work best for systems with slowly varying densities (i.e., as close to a homogeneous electron gas as possible).

However, it is interesting that even for many systems where the density varies considerably, the LDA approach performs well!

Good bad local density approximation

Good & Bad - Local Density Approximation

Total energies

Structures of highly correlated systems (transition metals, FeO, NiO, predicts the non-magnetic phase of iron to be ground state)

Doesn’t describe weak interactions well.

Makes hydrogen bonds stronger than they should be.

Band gaps (shape and position is pretty good, but will underestimate gaps by roughly a factor of two; will predict metallic structure for some semiconductors)

  • Ground state densities well-represented

Cohesive energies are pretty good; LDA tends to overbind a system (whereas HF tends to underbind)

Bond lengths are good, tend to be underestimated by 1-2%

Good for geometries, vibrations, etc.

Local density approximation3

Local Density Approximation

One place where LDA performs poorly is in the calculation of excited states.

Beyond lda

Beyond LDA

First principles electronic structure density functional theory

  • There are many different density functionals

  • Each works best in different situations

  • Difficult to know if it worked for the right reasons

  • Practically, very good accuracy/computational cost

  • Login