Dynamical mean field theory

Dynamical mean field theory

Dynamical Mean Field Theory (DMFT) is a method to determine the electronic structure of strongly correlated materials. In such materials, the approximation of independent electrons, which is used in Density Functional Theory and usual band structure calculations, breaks down. Dynamical Mean-Field Theory, a non-perturbative treatment of local interactions between electrons, bridges the gap between the nearly-free electron gas limit and the atomic limit of condensed-matter physics.[1]

DMFT consists in mapping a many-body lattice problem to a many-body local problem, called an impurity model.[2] While the lattice problem is in general intractable, the impurity model is usually solvable through various schemes. The mapping in itself does not constitute an approximation. The only approximation made in ordinary DMFT schemes is to assume the lattice self-energy to be a momentum-independent (local) quantity. This approximation becomes exact in the limit of lattices with an infinite coordination.

One of DMFT's main successes is to describe the phase transition between a metal and a Mott insulator when the strength of electronic correlations is increased. It has been successfully applied to real materials, in combination to the local density approximation of Density Functional Theory.

Contents

Relation to Mean-Field Theory

The DMFT treatment of lattice quantum models is similar to the mean-field theory (MFT) treatment of classical models such as the Ising model.[3] In the Ising model, the lattice problem is mapped onto an effective single site problem, whose magnetization is to reproduce the lattice magnetization through an effective "mean-field". This condition is called the self-consistency condition. It stipulates that the single-site observables should reproduce the lattice "local" observables by means of an effective field. While the N-site Ising Hamiltonian is hard to solve analytically (to date, analytical solutions exist only for the 1D and 2D case), the single-site problem is easily solved.

Likewise, DMFT maps a lattice problem (e.g the Hubbard model) onto a single-site problem. In DMFT, the local observable is the local Green's function. Thus, the self-consistency condition for DMFT is for the impurity Green's function to reproduce the lattice local Green's function through an effective mean-field which, in DMFT, is the hybridization function Δ(τ) of the impurity model. DMFT owes its name to the fact that the mean-field Δ(τ) is time-dependent, or dynamical. This also points to the major difference between the Ising MFT and DMFT: Ising MFT maps the N-spin problem into a single-site, single-spin problem. DMFT maps the lattice problem onto a single-site problem, but the latter fundamentally remains a N-body problem which captures the temporal fluctuations due to electron-electron correlations.

Description of DMFT for the Hubbard Model

The DMFT mapping

Single-orbital Hubbard model

The Hubbard model [4] describes the onsite interaction between electrons of opposite spin by a single parameter, U. The Hubbard Hamiltonian may take the following form:

 H_{\text{Hubbard}}=t \sum_{\langle ij \rangle \sigma}  c_{i\sigma}^{\dagger}c_{j\sigma} + U\sum_{i}n_{i \uparrow} n_{i\downarrow}

where c_i^{\dagger},c_i denote the creation and annihilation operators of an electron on a localized orbital on site i, and n_i=c_i^{\dagger}c_i.

The following assumptions have been made:

  • only one orbital contributes to the electronic properties (as might be the case of copper atoms in superconducting cuprates, whose d-bands are non-degenerate),
  • the orbitals are so localized that only nearest-neighbor hopping t is taken into account

The auxiliary problem: the Anderson impurity model

The Hubbard model is in general intractable under usual perturbation expansion techniques. DMFT maps this lattice model onto the so-called Anderson impurity model (AIM). This model describes the interaction of one site (the impurity) with a "bath" of electronic levels (described by the annihilation and creation operators ap and a_p^{\dagger}) through an hybridization function. The Anderson model corresponding to our single-site model is a single-orbital Anderson impurity model, whose hamiltonian formulation is the following:

H_{\text{AIM}}=\underbrace{\sum_{p}\epsilon_p a_p^{\dagger}a_p}_{H_{\text{bath}}} + \underbrace{\sum_{p\sigma}\left(V_{p}^{\sigma}c_{\sigma}^{\dagger}a_{p\sigma}+h.c\right)}_{H_{\text{mix}}}+\underbrace{U n_{\uparrow} n_{\downarrow}-\mu \left(n_{\uparrow}+n_{\downarrow}\right)}_{H_{\text{loc}}}

where

  • Hbath describes the non-correlated electronic levels \epsilon_p of the bath
  • Hloc describes the impurity, where two electrons interact with the energetical cost U
  • Hmix describes the hybridization (or coupling) between the impurity and the bath through hybridization terms V_p^{\sigma}

The Matsubara Green's function of this model, defined by  G_{\text{imp}}(\tau) = - \langle T c(\tau) c^{\dagger}(\tau)\rangle , is entirely determined by the parameters U and the so-called hybridization function  \Delta_\sigma(i\omega_n) = \sum_{p}\frac{|V_p^\sigma|^2}{i\omega_n-\epsilon_p}, which is the imaginary-time Fourier-transform of Δσ(τ).

This hybridization function describes the dynamics of electrons hopping in and out of the bath. It should reproduce the lattice dynamics such that the impurity Green's function is the same as the local lattice Green's function. It is related to the non-interacting Green's function by the relation:

(G0) − 1(iωn) = iωn + μ − Δ(iωn) (1)

Solving the Anderson impurity model consists in computing observables such as the interacting Green's function G(iωn) for a given hybridization function Δ(iωn) and U. It is a difficult but not intractable problem. There exists a number of ways to solve the AIM, such as

  • Exact diagonalization
  • Hirsch-Fye algorithm[5]
  • Iterative Perturbation Theory
  • Non-Crossing Approximation
  • Continuous-Time Quantum Monte Carlo algorithms[6]

Self-consistency equations

The self-consistency condition requires the impurity Green's function Gimp(τ) to coincide with the local lattice Green's function G_{ii}(\tau) = -\langle T c_i(\tau)c_i^{\dagger}(0)\rangle :

 G(i\omega_n) = G_{ii}(i\omega_n) = \sum_k \frac {1}{i\omega_n +\mu - \epsilon(k) - \Sigma(k,i\omega_n)}

where Σ(k,iωn) denotes the lattice self-energy.

DMFT approximation: locality of the lattice self-energy

The only DMFT approximations (apart from the approximation that can be made in order to solve the Anderson model) consists in neglecting the spatial fluctuations of the lattice self-energy, by equating it to the impurity self-energy:

 \Sigma(k,i\omega_n) \approx \Sigma_{imp}(i\omega_n)

This approximation becomes exact in the limit of lattices with infinite coordination, that is when the number of neighbors of each site is infinite. Indeed, one can show that in the diagrammatic expansion of the lattice self-energy, only local diagrams survive when one goes into the infinite coordination limit.

Thus, as in classical mean-field theories, DMFT is supposed to get more accurate as the dimensionality (and thus the number of neighbors) increases. Other put, for low dimensions, spatial fluctuations will render the DMFT approximation less reliable.

The DMFT Loop

In order to find the local lattice Green's function, one has to determine the hybridization function such that the corresponding impurity Green's function will coincide with the sought-after local lattice Green's function. An intuitive method would be the following: for a given U , μ and temperature T

  1. First, compute the non-interacting lattice Green's function G0(k,iωn), and extract its local part Gloc
  2. The self-consistency condition requires it to be equal to the impurity Green's function,  \mathcal{G}^0(\tau) = G_{loc}
  3. Compute the corresponding hybridization function through (1)
  4. Solve the AIM for a new impurity Green's function G_{imp}^0(\tau), extract its self-energy: \Sigma_{imp}(i\omega_n) = (G_{imp}^0)^{-1}(i\omega_n) -  (\mathcal{G}^0)^{-1}(i\omega_n)
  5. Make the DMFT approximation:  \Sigma(k,i\omega_n) \approx \Sigma_{imp}(i\omega_n)
  6. Compute the new lattice Green's function, extract its local part and go back to step 2 with a new  \mathcal{G}^1

Self-consistency is reached when G_{imp}^n = G_{imp}^{n+1}.

Applications

The local lattice Green's function and other impurity observables can be used to calculate a number of physical quantities as a function of correlations U, bandwidth, filling (chemical potential μ), and temperature T:

  • the spectral function (which gives the band structure)
  • the kinetic energy
  • the double occupancy of a site
  • response functions (compressibility, optical conductivity, specific heat)

In particular, the drop of the double-occupancy as U increases is a signature of the Mott transition.

Extensions of DMFT

DMFT has several extensions, extending the above formalism to multi-orbital, multi-site problems.[7]

Multi-orbital extension

DMFT can be extended to Hubbard models with multiple orbitals, namely with electron-electron interactions of the form Uαβnαnβ where α and β denote different orbitals. This is especially relevant for compounds whose d-orbitals are degenerate, such as iron in the newly discovered high-temperature iron-based superconductors (pnictides).

Cluster DMFT

In order to improve on the DMFT approximation, the Hubbard model can be mapped on a multi-site impurity (cluster) problem, which allows one to add some spatial dependence to the impurity self-energy. Typical clusters contain 4 to 8 sites.

Extended DMFT

DMFT can be applied to more general models such as the t-J model.

References and notes

  1. ^ Antoine Georges et al. (1996). "Dynamical mean-field theory of strongly correlated fermion systems and the limit of infinite dimensions". Reviews of Modern Physics 68 (1): 13. doi:10.1103/RevModPhys.68.13. 
  2. ^ A. Georges and G.Kotliar (1992). "Hubbard model in infinite dimensions". Physical Review B 45 (12): 6479. doi:10.1103/PhysRevB.45.6479. 
  3. ^ Antoine Georges (2004). Strongly Correlated Electron Materials: Dynamical Mean-Field Theory and Electronic Structure. "AIP Conference Proceedings". Lectures on the Physics of Highly Correlated Electron Systems VIII , American Institute of Physics Conference Proceedings Vol. 715 (715): 3–74. arXiv:cond-mat/0403123. doi:10.1063/1.1800733. 
  4. ^ John Hubbard (1963). "Electron Correlations in Narrow Energy Bands". Proceedings of the Royal Society A 276 (1365): 238. doi:10.1098/rspa.1963.0204. 
  5. ^ J.E. Hirsch and R.M. Fye (1986). "Monte Carlo Method for Magnetic Impurities in Metals". Physical Review Letters 56 (23): 2521. doi:10.1103/PhysRevLett.56.2521. 
  6. ^ A. N. Rubtsov, V. V. Savkin, and A. I. Lichtenstein (2005). "Continuous-time quantum Monte Carlo method for fermions". Physical Review B 72 (3): 035122. doi:10.1103/PhysRevB.72.035122. 
  7. ^ Vladimir Anisimov, Yuri Izyumov (2010). "DMFT Extensions". Electronic Structure of Strongly Correlated Materials. Springer Series in Solid State Sciences. ISBN 3642048250. http://books.google.com/?id=8YqmOcs5jggC&printsec=frontcover&dq=Electronic+Structure+of+Strongly+Correlated+Materials#v=onepage&q&f=false. 

See also

  • Strongly correlated materials

External links


Wikimedia Foundation. 2010.

Игры ⚽ Поможем сделать НИР

Look at other dictionaries:

  • Mean field theory — (MFT, also known as self consistent field theory) is a method to analyse physical systems with multiple bodies. A many body system with interactions is generally very difficult to solve exactly, except for extremely simple cases (random field… …   Wikipedia

  • Density functional theory — Electronic structure methods Tight binding Nearly free electron model Hartree–Fock method Modern valence bond Generalized valence bond Møller–Plesset perturbation theory …   Wikipedia

  • Many-body theory — The many body theory (or Many body physics) is an area of physics which provides the framework for understanding the collective behavior of vast assemblies of interacting particles. In general terms, the many body theory deals with effects that… …   Wikipedia

  • Theory of relativity — This article is about the scientific concept. For philosophical or sociological theories about relativity, see Relativism. For the silent film, see The Einstein Theory of Relativity. Two dimensional projection of a three dimensional analogy of… …   Wikipedia

  • dynamical time —       specialized timescale used to describe the motion of objects in space.       As a practical matter, time can be defined as that coordinate which can most simply be related to the evolution of closed systems. Proper time is the time measured …   Universalium

  • Number theory — A Lehmer sieve an analog computer once used for finding primes and solving simple diophantine equations. Number theory is a branch of pure mathematics devoted primarily to the study of the integers. Number theorists study prime numbers (the… …   Wikipedia

  • Le Sage's theory of gravitation — is the most common name for the kinetic theory of gravity originally proposed by Nicolas Fatio de Duillier in 1690 and later by Georges Louis Le Sage in 1748. The theory proposed a mechanical explanation for Newton s gravitational force in terms… …   Wikipedia

  • List of physics topics A-E — NOTOC A B C D E F G H I J K L M N O P Q R S T U V W X Y ZList of physics topics F L >>AAb Am* Abbe number * Abbe, Ernst * Absolute zero * Accelerating universe * Acceleration * Acoustic theory * Action (physics) * Active laser medium * Adiabatic… …   Wikipedia

  • Dieter Vollhardt — (* 8. September 1951) ist ein deutscher Physiker und seit 1996 Professor in Augsburg. Inhaltsverzeichnis 1 Leben 2 Werk 3 Schriften 4 …   Deutsch Wikipedia

  • Mott insulator — Mott insulators are a class of materials that should conduct electricity under conventional band theories, but are insulators when measured (particularly at low temperatures). This effect is due to electron electron interactions which are not… …   Wikipedia

Share the article and excerpts

Direct link
Do a right-click on the link above
and select “Copy Link”