Computational Materials Science
Dynamical mean field theory and its extensions
Dynamical mean field theory (DMFT) and its numerous extensions represent a big step forward in our efforts to develop reliable methods for describing electronic correlations. Depending on the strength of the electronic correlation, these non-perturbative methods correctly yield a weakly correlated metal, a strongly correlated metal, or a Mott insulator. Applied to model Hamiltonians, DMFT was a big success and advanced our understanding of correlation effects like the Mott-Hubbard metal-insulator transition. More recently, extensions of DMFT such as the so-called quantum cluster theories as well as other extensions based on diagrammatic approaches have been quickly unveiling the missing aspects of the complete picture, where also non-local physics beyond DMFT is treated.
DMFT has been successfully merged with conventional bandstructure calculations in the local density approximation (LDA), thereby allowing for the realistic calculation of materials with strong electronic correlations such as transition metal oxides or heavy Fermion systems, which was hitherto not possible.
Electronic Structure Calculations with DMFT
One of the most important challenges of theoretical physics is the development of tools for the accurate calculation of material properties. In solid state theory, we already know our ''theory of everything'', which consists of three terms: the kinetic energy, the lattice potential, and the Coulomb interaction between electrons:
![]() | (1) |
Here, and
denote the position of electron
and ion
with charge
and
, respectively,
is the Laplace operator of the kinetic energy,
and
are the vacuum dielectric and Planck constant, see Fig. 1 for an illustration. Despite knowing the Hamiltonian, it is impossible to solve, even numerically, if more than a very few electrons are involved. This is due to the last term, the Coulomb interaction, which correlates the movement of every electron
with every other electron
.
Since electrons have to be described quantum-mechanically, the numerical effort grows exponentially with the number of electrons.
Solid State Hamiltonian To calculate a material, one has to take into account three terms: the kinetic energy favoring the electrons to move through the crystal (black), the lattice potential of the ions (green), and the Coulomb interaction between the electrons (red). Due to the latter, the first electron is repelled by the second, and it is energetically favorable to hop somewhere else, as depicted. Hence, the movement of every electron is correlated with that of every other, which prevents even a numerical solution. | LDA approximation LDA is an approximation which allows to calculate material properties but which dramatically simplifies the electronic correlations: every electron moves independently, i.e., uncorrelated, within a time-averaged local density of the other electrons. |
In this situation, one can either dramatically simplify the Hamiltonian (1), hoping that the simplified model allows for a qualitative understanding including correlation (many-body) effects, or employ equally dramatic approximations to deal with (1). These two paths have been followed by the two big communities of solid state theory, the many-body model-Hamiltonian and the density functional community.
Within density functional theory, the local density approximation (LDA) turned out to be unexpectedly successful, and established itself as the method for realistic solid state calculations in the last century. This is surprising because LDA is a serious approximation to the Coulomb interaction. Basically, an electron at sees a time-averaged density of the other electrons
, with a corresponding local LDA potential
, see Fig. 1. This reduces the problem to a single electron calculation.
The success of LDA shows, that this treatment is actually sufficient for many materials, both for calculating ground state energies and bandstructures, implying that electronic correlations are rather weak in these materials. But, there are important classes of materials where LDA fails, such as transition metal oxides or heavy fermion systems. In these materials the valence orbitals are the and
orbitals. For two electrons in these orbitals the distance
is particularly short, and electronic correlations particularly strong.
Many such transition metal oxides are Mott insulators, where, the on-(lattice-)site Coulomb repulsion splits the LDA bands into two sets of Hubbard bands. One can envisage the lower Hubbard band as consisting of all states with one electron on every lattice site and the upper Hubbard band as those states where two electrons are on the same lattice site. Since it costs an energy
to have two electrons on the same lattice sites, the latter states are completely empty and the former completely filled with a gap of size
in-between, see Fig. 2. Other transition metal oxides and heavy fermion systems are strongly correlated metals, with heavy quasiparticles at the Fermi energy, described by an effective mass or inverse weight
, see Fig. 2.
Weakly correlated metal The on-site Coulomb interaction | Strongly correlated metal In this intermediate regime, one has already Hubbard bands, like for | Mott insulator If the Coulomb interaction becomes large ( |
The entire parameter regime and physics is described by LDA |
Dynamical mean-field theory (DMFT) is a modern, non-perturbative many-body technique, which describes such kind of correlation physics. It was first applied to model Hamiltonians - with big success. Hence, DMFT was recently merged with the conventional LDA approach for calculating materials with strong electronic correlations realistically. The first LDADMFT calculations turned out to be a breakthrough in this respect, which might not be astonishing as it is able to describe the full range of materials from weakly to strongly correlated metals to Mott insulators, see Fig. 2. LDA
DMFT was successful in calculating, among others, the
-
transition in cerium, the metal-insulator transition in V
O
, and the
-phase of plutonium.
Outlook
Electronic structure calculations with DMFT have just begun. Many important physical quantities such as the optical and thermal conductivity, magnetic and orbital susceptibilities, and thermodynamic properties have not yet been calculated or only for a single material. Many important transitions and the physics beyond have neither yet been studied with LDADMFT nor are they understood. Besides LDA
DMFT, other methods employing DMFT to do electronic structure calculations are presently developed. Generalizations to extend DMFT towards cluster algorithms, which can describe, e.g., d-wave superconductivity, have been established and have been already successfully applied to model Hamiltonians.