---
title: Hubbard-Corrected Energy Functionals
url: https://www.emergentmind.com/topics/hubbard-corrected-energy-functionals
type: topic
---

# Hubbard-Corrected Energy Functionals

Hubbard-corrected energy functionals, commonly denoted as LDA+U or, in generalized forms, DFT+U(+J/V), are a class of extended density functionals designed to remedy the failures of local and semilocal density functional theory (DFT) in describing systems with localized electrons—particularly in correlated d and f orbitals of transition metals and rare-earth elements. These methods supplement standard DFT energies with an explicit on-site Coulomb term and an associated double-counting correction, targeting spurious delocalization and self-interaction errors, and have been further extended to more elaborate functionals capturing exchange and intersite effects. Hubbard corrections have evolved into essential computational tools to describe Mott insulators, charge-transfer systems, magnetic order, and beyond, balancing accuracy against computational cost in strongly correlated materials.

## 1. Theoretical Foundation and Formulation

The central feature of Hubbard-corrected functionals is their additive structure. An archetypical formulation is
$$
E_{LDA+U}[\rho] = E_{LDA}[\rho] + E_{\text{Hub}}[{n^{I\sigma}_{mm'}}] - E_{dc}[{n^I}],
$$
with $E_{LDA}$ the standard (local or semilocal) DFT energy, $E_{\text{Hub}}$ the Hubbard on-site (Coulomb/interaction) term, and $E_{dc}$ a double-counting term to subtract the portion of electron–electron interaction already included in $E_{LDA}$ in an averaged manner.

A common "rotationally invariant" form (Dudarev et al.) is:
$$
E_U = \sum_{I, \sigma} \frac{U^I}{2} \operatorname{Tr}\left[\mathbf{n}^{I\sigma}(1 - \mathbf{n}^{I\sigma})\right],
$$
where $\mathbf{n}^{I\sigma}$ is the occupation matrix projected onto a localized basis (typically atomic d or f orbitals) for site $I$ and spin $\sigma$. The occupation matrix elements are obtained by projecting Kohn–Sham orbitals onto localized orbitals,
$$
n^{I\sigma}_{mm'} = \sum_{k,v} f_{kv}^\sigma \langle\psi_{kv}^\sigma|\phi_{m'}^I\rangle \langle\phi_m^I|\psi_{kv}^\sigma\rangle.
$$
Differentiation gives rise to a Hubbard correction in the Kohn–Sham potential, penalizing deviations from integer orbital occupations, and thereby restoring the Mott–Hubbard insulating state where standard DFT fails.

Extensions include fully rotationally invariant functionals:
$$
E_{U+J} = \sum_{I, \sigma} \frac{U^I - J^I}{2} \operatorname{Tr}\left[\mathbf{n}^{I\sigma}(1 - \mathbf{n}^{I\sigma})\right] + \sum_{I,\sigma} \frac{J^I}{2} \operatorname{Tr}\left[\mathbf{n}^{I\sigma}\mathbf{n}^{I,-\sigma}\right],
$$
where $J^I$ encodes Hund's exchange. The double-counting correction is usually handled by either the "fully localized limit" (FLL) or "around mean field" (AMF) schemes, reflecting different assumptions on electronic localization.

The Hubbard $U$ parameter is often determined by a linear-response approach, which gives:
$$
U^I = (\chi_0^{-1} - \chi^{-1})_{II},
$$
with $\chi$ the fully screened, and $\chi_0$ the bare (non-self-consistent) susceptibility matrix measuring the response of $n^I$ to perturbations in the potential on localized sites.

## 2. Applications to Correlated Materials

The LDA+U functional is widely used for systems where conventional DFT underperforms:

- **Transition-metal oxides** (e.g., NiO, MnO, FeO, CuO): introducing $U$ opens insulating band gaps, restores the charge-transfer gap, and improves magnetic state prediction. For example, in NiO, $U$ shifts the O $p$ and Ni $d$ states to correct band alignment.
- **Defective oxides** (e.g., reduced CeO$_2$): LDA+U localizes excess electrons on $f$ or $d$ orbitals.
- **Intermetallics** (e.g., Ni$_2$MnGa): LDA+U modifies the balance of phases by shifting localized $d$ levels.
- **Metal complexes and molecules**: LDA+U and its extensions (incorporating intersite $V$ or Hund’s $J$) improve the description of states with strong metal–ligand charge transfer, as in phosphorescent Ir complexes.
- **High-throughput screening**: Due to modest computational overhead compared to plain DFT, LDA+U has facilitated screening of correlated materials and phase diagrams.

The method is applicable to both localized insulators and systems with variable oxidation states that require discrimination between integer valence configurations.

## 3. Approximations, Ambiguities, and Comparisons

Key approximations and issues include:

- **Mean-field treatment**: The on-site Coulomb correction is factorized at the mean-field level, replacing four-operator expectation values with products of two occupation matrices.
- **Occupation matrix truncation**: Often, off-diagonal terms in $n_{mm'}$ are neglected or only averaged, and Hund’s $J$ is sometimes absorbed into an effective $U_{\text{eff}} = U - J$.
- **Double-counting ambiguity**: $E_{dc}$ is not uniquely defined. FLL is suited for systems with localized electrons, AMF for overdoped/metals.
- **Linear-response U**: The value of $U$ obtained may depend on definition of the localized orbitals, the symmetry, and the chemical environment.
- **Neglected frequency dependence**: Standard DFT+U is "static", not addressing dynamical correlations. Comparisons with hybrid functionals and DFT+DMFT indicate that while LDA+U is cost-effective, it lacks dynamical self-energies and many-body fluctuations.
- **Limitations for metals**: For itinerant or metallic systems, the penalization of fractional occupation can artificially open gaps unless symmetry is broken.

In comparison, hybrid functionals and DMFT include a fraction of explicit exchange or dynamic many–body effects, but at increased computational cost.

## 4. Recent Extensions: LDA+U+V, Hund’s J, and Beyond

Several extensions have been formulated to improve predictive capability:

- **LDA+U+V**: Incorporates intersite Coulomb interactions ($V$), correcting not only on-site self-interaction but also the covalency and hybridization between neighboring atomic sites. Explicitly,
  $$
  E_{UV} = \sum_{I,\sigma} \frac{U^I}{2} \operatorname{Tr}\left[n^{II,\sigma}(1 - n^{II,\sigma})\right] - \sum_{I,J,\sigma}^{*} \frac{V^{IJ}}{2} \operatorname{Tr}\left[n^{IJ,\sigma} n^{JI,\sigma}\right],
  $$
  where the asterisk indicates summation over neighboring $I, J$ pairs.
- **Hund’s coupling ($J$) explicit**: For systems with competing spin and orbital orders (e.g., cubic CuO), a proper $J$ is required for the system to become insulating or to set the correct spin state.
- **Analytic derivatives**: The form of the energy functional allows for straightforward analytical derivation of forces, stress tensors, and even phonon spectra in the presence of Hubbard corrections.
- **Self-consistency of $U$ and $V$**: Linear-response theory is iterated self-consistently over geometry and occupation, ensuring state-dependent and environment-sensitive parameters.

## 5. Open Issues and Outstanding Questions

Key unresolved aspects remain:

- **Calculation and transferability of interaction parameters**: $U$ and $J$ depend sensitively on the choice of localized basis (atomic vs. Wannier), as well as on geometry and magnetic/oxidation state, complicating their unique assignment in practice.
- **Double-counting corrections**: No consensus exists for a universal, first-principles double-counting functional. Different correction schemes produce divergent results in some materials.
- **Appropriate treatment of metals**: The static nature of LDA+U makes its application to metallic phases problematic or ambiguous.
- **Frequency dependence and dynamical screening**: The need to connect with methods such as DFT+DMFT is highlighted, and dynamic corrections remain a target for functional development.
- **Symmetry and basis invariance**: The correction may depend on basis rotations or selections, potentially threatening the method’s universality.
- **Extension to multi-site and multi-orbital interactions**: For highly delocalized or low-symmetry systems, additional interaction channels may need to be considered.

## 6. Methodological Advances and Future Directions

The paper outlines several future priorities:

- **Automated, ab initio determination of $U$, $J$, and $V$** that dynamically follow changes in structure and local environment, enabling adaptive corrections during molecular dynamics and structural optimization.
- **Development of more flexible functionals** with multi-site and multi-orbital terms, to capture both localized and itinerant electronic correlation, especially in low-symmetry or mixed-valence materials.
- **Systematic derivation of double-counting terms** that are transferable across materials classes, bridging strongly localized insulators and metallic, paramagnetic, or even high-entropy systems.
- **Bridging static and dynamic approaches**: Efforts to incorporate dynamical screening, either as an explicit frequency dependence or through a parametric extension, to recover more of the physics addressed in DMFT but at lower computational cost.
- **Application in high-throughput materials discovery**, where fast, transferable corrections are needed for reliable property prediction.

## 7. Summary and Outlook

Hubbard-corrected (LDA+U) energy functionals provide an essential correction to conventional DFT, enabling accurate treatment of correlated electrons by penalizing fractional occupancy in localized subspaces. Their success is notable in transition-metal oxides, rare-earth compounds, intermetallic alloys, magnetic systems, and excited-state calculations, where they enable improved predictions of gaps, charge localization, and relative phase stability at negligible incremental computational expense. Limitations—centered on parameter determination, double-counting, basis dependence, and handling of dynamical effects—persist and motivate ongoing methodological innovation, including self-consistent and adaptive $U$, generalized $J$ and $V$ extensions, improved double-counting corrections, and dynamical variants approaching DMFT-level accuracy. The evolution towards generalized, automated, and more transferable functionals will further enlarge the method’s reach in complex and previously intractable correlated materials [1309.3355].

Source: https://www.emergentmind.com/topics/hubbard-corrected-energy-functionals