---
title: 'LDA+DMFT: Modeling Correlated Electron Materials'
url: https://www.emergentmind.com/topics/lda-dmft-local-density-approximation-plus-dynamical-mean-field-theory
type: topic
---

# LDA+DMFT: Modeling Correlated Electron Materials

The LDA+DMFT (Local Density Approximation plus Dynamical Mean-Field Theory) approach is a leading first-principles framework for the quantitative study of correlated electron materials. By combining density-functional theory in the local density approximation (LDA) with dynamical mean-field theory (DMFT), LDA+DMFT enables ab initio modeling of complex materials where strong, local electronic interactions and itinerant band-structure effects are both essential. The formalism was developed to overcome the limitations of LDA and related density-functional techniques, which underestimate electronic correlations in $d$- and $f$-electron systems, and to move beyond purely model-based DMFT by supplying realistic, material-specific one-particle Hamiltonians.

## 1. Theoretical Principles of LDA+DMFT

LDA+DMFT starts from a Kohn–Sham (KS) band Hamiltonian $H_{\mathrm{KS}}$ derived from conventional LDA or generalized gradient approximation calculations. This Hamiltonian is represented in a localized basis (often Wannier or linear muffin-tin orbital functions), allowing for the identification of a subset of "correlated" orbitals (typically $3d$, $4f$, or $5f$ shells) on each atomic site $R$ [1110.2606, 1411.6906]. To $H_{\mathrm{KS}}$ is added a local interaction term of the Hubbard–Hund type:
\[
H_{\mathrm{int}} = \frac{1}{2} \sum_{R,\,\xi\xi',\,\sigma\sigma'} U_{\xi\xi'\sigma\sigma'} c^\dagger_{R\xi\sigma} c^\dagger_{R\xi'\sigma'} c_{R\xi'\sigma'} c_{R\xi\sigma}
\]
together with a double-counting correction $H^{DC}_R$, subtracting the LDA-included fraction of the local interaction energy
[1110.2606, 1411.6906, 1403.2474].

The effective many-body problem thus reads:
\[
H = H_{\mathrm{KS}} + \sum_{R} (U_{R} - H^{DC}_R)
\]
In DMFT, the static mean-field $U$ is replaced by a frequency-dependent, orbital-local self-energy $\Sigma_{\mathrm{loc}}(\omega)$ restricted to the correlated subspace. The lattice Green's function is computed as:
\[
G_k(\omega) = [\omega + \mu - H_{\mathrm{KS}}(k) - \Sigma_{\mathrm{loc}}(\omega)]^{-1}
\]
The DMFT loop enforces self-consistency between the lattice and an auxiliary impurity problem, ensuring that the impurity Green's function matches the on-site projected lattice Green's function [1110.2606, 1003.3600, 1412.8405].

## 2. Computational Workflow and Charge Self-Consistency

The LDA+DMFT computational cycle comprises two nested self-consistency loops: (i) the inner DMFT loop finding a self-consistent impurity self-energy, and (ii) the outer loop updating the total electron density to achieve full charge self-consistency (CSC) [1110.2606, 1111.2157, 1411.6906].

- **Initialization**: An initial charge density $\rho(\mathbf{r})$ is constructed and the KS band structure $H_{\mathrm{KS}}(k)$ is generated.
- **Correlated Subspace Construction**: Projectors $P_{R, \xi}(k)$ define the local correlated orbitals from the electronic structure basis.
- **DMFT Loop**: For each inequivalent correlated site, the local Green's function is projected and the Weiss field (noninteracting bath Green's function) $\mathcal{G}_0^{-1}(\omega)$ is obtained. The impurity problem is solved (typically by CT-HYB quantum Monte Carlo, SPTF, or exact diagonalization), yielding a new impurity self-energy $\Sigma_{\mathrm{imp}}(\omega)$.
- **Self-Consistency Condition**: Impose $\Sigma_{\mathrm{loc}}(\omega) = \Sigma_{\mathrm{imp}}(\omega)$ and $G_{\mathrm{loc}}(\omega) = G_{\mathrm{imp}}(\omega)$.
- **Charge-Density Update**: Reconstruct the density matrix in the Bloch basis,
  \[
  N_{nn'}(k) = -\frac{1}{\pi} \int d\omega\, f(\omega) \mathrm{Im} G_{k,nn'}(\omega + i0^+)
  \]
  and recover the new electron density
  \[
  \rho(r) = \sum_{k, n, n'} \psi_{k, n}(r) N_{nn'}(k) \psi_{k, n'}^*(r)
  \]
  [1110.2606, 1111.2157].
- **Total Energy**: The total LDA+DMFT energy is evaluated as
  \[
  E[\rho, \Sigma] = E_{\mathrm{LDA}}[\rho] - \langle H_{\mathrm{KS}} \rangle_{\mathrm{LDA}} + \langle H_{\mathrm{KS}} \rangle_{\mathrm{DMFT}} + \frac{1}{2} \mathrm{Tr}[\Sigma G_{\mathrm{loc}}] - E_{\mathrm{DC}}
  \]
  [1111.2157].
- **Convergence**: The process continues until both the electron density and the impurity self-energy are converged with respect to a tolerance.

Charge self-consistency is crucial to ensure the feedback of the DMFT local self-energy to the entire electronic structure, affecting both total energies and spectral properties [1110.2606, 1111.2157].

## 3. Double-Counting and Projected Subspace Strategies

A central issue in LDA+DMFT is the identification of the optimal double-counting correction, which must subtract those static correlation effects already accounted for in LDA [1403.2474, 1204.2361]. Frequently used forms include the fully-localized limit (FLL):
\[
H_{\mathrm{DC}} = \frac{1}{2} U\, n_{d}(n_{d}-1) - \frac{1}{2} J \sum_\sigma n_{d\sigma}(n_{d\sigma} - 1)
\]
and "around mean field" (AMF), based on a uniform occupancy assumption.

The recent "exact double counting" procedure [1403.2474] recasts the correlation functionals such that the intersection (double-counting) term is given by the LDA correlation evaluated on the local density. For solids, screening is incorporated via a screened interaction $U^\lambda_C(r-r')$, and the double counting is defined with the correlation energy per particle of a screened homogeneous electron gas.

Projector definitions for the correlated subspace can be constructed via maximally-localized Wannier functions, muffin-tin orbitals, or atomic-like pseudo-orbitals, with the projection formalism determining the degree of overlap between DFT and DMFT treatments [1111.2157, 0801.4353]. The flexible choice of projector and basis set allows adaptation to both all-electron and plane-wave pseudopotential codes.

## 4. Materials Applications and Physical Properties

LDA+DMFT has been applied to a wide range of correlated materials—Mott insulators, itinerant magnets, mixed-valence compounds, rare-earths, and actinide oxides—with both spectral and thermodynamic observables accurately computed [1110.2606, 1002.4947, 1411.6906]. Specific examples include:

- **NiO (Charge-Transfer Insulator)**: Full CSC LDA+DMFT corrects the 3d occupancy, enlarges the Mott gap to $\sim$4 eV, and reproduces experimental spin and orbital moments per Ni (e.g., $m_s\approx1.73\,\mu_B$, $m_o\approx0.24\,\mu_B$) [1110.2606].
- **bcc Fe (Itinerant Ferromagnet)**: Spectral functions and mass enhancement are quantitatively described, with spin ($m_s\approx2.14\,\mu_B$) and orbital ($m_o\approx0.06\,\mu_B$) moments in good agreement with experiment [1110.2606].
- **SmCo$_5$ (Permanent Magnet)**: Simultaneous treatment of Sm $4f$ and Co $3d$ electrons (via different impurity solvers) reveals the necessity of CSC for correct multiplet positions and total magnetization ($\sim8.0\,\mu_B$/f.u.) [1110.2606].
- **Fe-pnictides (e.g., LaFePO)**: LDA+DMFT accounts for moderate correlation-induced mass enhancements ($m^*/m \approx 1.9$–2.2) and reproduces ARPES band renormalization [1002.4947].

The DMFT description yields dynamic local self-energies that redistribute spectral weight between quasiparticle bands and Hubbard satellites, capture temperature-driven Mott metal-insulator transitions, and enable quantitative modeling of valence and core-level spectroscopies, including X-ray photoemission and RIXS [1704.01500, 1911.10366, 1812.06432].

## 5. Extensions: Core-Level Spectroscopies and Dynamical Screening

LDA+DMFT has been extended to model spectroscopies beyond one-particle spectral functions. For core-level XPS and RIXS, an Anderson impurity model is constructed with explicit core-shells (e.g., $2p$) and multiplet coupling to the valence shell, using a hybridization bath derived from DMFT [1704.01500, 1911.10366, 2004.01428]. Such methods quantitatively reproduce satellite features, nonlocal screening effects, and fluorescence-like continua in experimental spectra, which are not accessible to traditional cluster models.

Dynamical screening effects—frequency dependence of the Hubbard $U$—are integrated by promoting $U \rightarrow U(\omega)$ and including the retarded interaction either via direct coupling to bosonic modes or within the "Bose factor ansatz" [1412.8405]. This enables the description of plasmon satellites and the reduction of the low-energy mass renormalization by a bandwidth renormalization factor $Z_B$.

## 6. Methodological Challenges and Perspectives

Open methodological questions remain regarding the optimal formulation of double-counting corrections, the unambiguous construction of correlated projectors, and the systematic inclusion of nonlocal correlations [1204.2361, 1411.6906]. Full charge self-consistency is critical for obtaining accurate forces and energetics, and has been demonstrated in all-electron LMTO, LAPW, and plane-wave pseudopotential codes [1110.2606, 1111.2157]. Efficient quantum impurity solvers (CT-HYB, SPTF, ED) are essential for practical calculations, dictating the accessible temperatures and number of correlated orbitals.

Extensions to GW+DMFT and inclusion of cluster or diagrammatic corrections beyond single-site DMFT are being actively developed to capture nonlocal fluctuations and screening. The parameter-free framework defined by the "exact double counting" scheme [1403.2474] addresses one of the longest-standing ambiguities in the field.

LDA+DMFT thus constitutes a predictive and material-specific theory for strongly correlated electrons, capable of addressing both spectroscopy and thermodynamics in real materials, and serves as a foundation for future developments in ab initio many-body electronic structure theory.

Source: https://www.emergentmind.com/topics/lda-dmft-local-density-approximation-plus-dynamical-mean-field-theory