---
title: Fermi-Löwdin-Orbital Self-Interaction Correction
url: https://www.emergentmind.com/topics/fermi-lowdin-orbital-self-interaction-correction-flosic
type: topic
---

# Fermi-Löwdin-Orbital Self-Interaction Correction

Fermi-Löwdin-Orbital Self-Interaction Correction (FLOSIC) is a formally unitary-invariant, size-extensive approach to correcting the unphysical self-interaction errors (SIE) present in common density functional approximations (DFAs). It implements the Perdew–Zunger self-interaction correction (PZ-SIC) using a localized orbital basis built from Fermi orbitals parameterized by real-space descriptors, known as Fermi-Orbital Descriptors (FODs), which are subsequently orthonormalized via the Löwdin procedure. This methodology reduces the scaling bottleneck of variational SIC (scaling as $O(N^2)$ to $O(N^3)$ in orbital rotations) to a much more tractable $O(N)$ real-space parameter optimization, ensuring computational viability for molecular and periodic systems. FLOSIC has facilitated major advances in the accuracy of calculated ionization potentials, dissociation energies, reaction barriers, charge localizations, and spectroscopic properties across atoms, molecules, and solids.

## 1. Formalism: SIC Functional and Fermi–Löwdin Construction

The PZ-SIC energy functional augments the standard Kohn–Sham total energy by subtracting for each occupied orbital the Hartree and exchange–correlation self-energies:
\[
E_{\rm SIC}[\{\phi_{i\sigma}\}] =
E_{\rm DFA}[n_\uparrow, n_\downarrow]
-\sum_{i\sigma}
\left\{
U[n_{i\sigma}]
+E_{\rm xc}^{\rm DFA}[n_{i\sigma},0]
\right\}
\]
where $n_{i\sigma}(r) = |\phi_{i\sigma}(r)|^2$, $U[n] = \frac12\int\int n(r)n(r')/|r-r'|\,dr\,dr'$, and $E_{\rm xc}^{\rm DFA}[n,0]$ is the XC functional of a one-electron density.

In FLOSIC, the $\{\phi_{i\sigma}\}$ are not arbitrary localized orbitals: instead, for each spin $\sigma$, the Fermi orbitals are defined using FODs $\{\mathbf{a}_{i\sigma}\}$ as
\[
F_{i\sigma}(r) = \frac{\sum_{j=1}^{N_\sigma} \psi_{j\sigma}^*(a_{i\sigma}) \psi_{j\sigma}(r)}{\sqrt{\sum_{j=1}^{N_\sigma} |\psi_{j\sigma}(a_{i\sigma})|^2}}
\]
where $\{\psi_{j\sigma}\}$ are the canonical KS orbitals. The set $\{F_{i\sigma}\}$ is Löwdin-orthonormalized to yield the final Fermi–Löwdin orbitals (FLOs) $\{\phi_{i\sigma}\}$ used in $E_{\rm SIC}$.

This procedure yields a size-extensive, unitarily invariant SIC functional dependent only on the $3N$ FOD coordinates, rather than the $O(N^2)$ possible occupied–occupied orbital rotations [1703.10742][1412.3101][1910.10677].

## 2. FOD Optimization and Self-Consistency

The total FLOSIC energy is minimized with respect to both:
- The electronic density (standard SCF procedure)
- The $3N$ FODs via their analytic gradients:
\[
\frac{\partial E_{\rm SIC}}{\partial a_{i\sigma}} = \text{analytic expression in terms of FOs, overlap derivatives, and SIC potentials}
\]
Algorithms alternate between optimizing the FODs at fixed density (often using gradient-based optimizers such as L-BFGS or conjugate-gradient) and SCF updates of the orbitals and density (including the SIC potentials from the current set of FODs) [1703.10742][1905.02631].

The frozen-density loop—a practical acceleration—fixes the density while optimizing FODs to convergence before recomputing the density. This approach can yield significant speedups, as demonstrated in converging electronic structures for transition metal macrocycles [2104.10327].

## 3. Advantages and Theoretical Implications

FLOSIC corrects major known deficiencies of DFAs:
- **Localized spin/charge densities:** Restoration of correct metal $d$-character and quantitation of spin populations in transition-metal complexes, as validated by comparison with multireference wavefunction calculations and EPR experiment [2104.10327].
- **Band gaps and orbital eigenvalues:** FLOSIC systematically opens underestimated HOMO–LUMO (Kohn–Sham) gaps and improves photoemission spectra, extending even to periodic systems if combined with a Wannier–Fermi–Löwdin (WFL) localization or effective potential approaches [2203.08953][2107.13631].
- **Reaction and dissociation energetics:** Barrier heights and dissociation curves—especially those suffering from stretched-bond SIE—are dramatically improved; the majority of the SIC correction to a reaction barrier is attributed quantitatively to the stretched-bond orbital at the transition state [2601.11454].

Size-extensivity and unitarity invariance make FLOSIC robust for both finite and extended systems [1703.10742][2203.08953].

## 4. Symmetry Breaking, Multiple Minima, and Ensemble Correction

The PZ-SIC functional and its FLOSIC implementation can admit multiple symmetry-broken minima due to overlocalization of electrons. In aromatic molecules, e.g., benzene, canonical FLOSIC produces localized Lewis/Kekulé structures, breaking $D_{6h}$ symmetry, or spin-symmetry broken Linnett double-quartet states (with incorrect $\langle S^2\rangle$). Which minimum is found depends on the FOD initialization (LT1, LT2, LDQ) and artifactual dipole moments or symmetry breaking are frequently observed [2109.08199][2405.18394].

The ensemble generalization, E-FLO-SIC, restores symmetry by averaging the SIC energy over all chemically equivalent localized FOD configurations
\[
E_{\rm SIC}^{\rm ens} = -\sum_{\sigma} \frac1{N_{\rm conf}^\sigma} \sum_{P=1}^{N_{\rm conf}^\sigma}\sum_{i=1}^{N_\sigma} \left\{ U[\rho_{i,\sigma}^P] + E_{\rm xc}[\rho_{i,\sigma}^P, 0] \right\}
\]
which yields electron densities and energies invariant under equivalent resonance transformations, as demonstrated for benzene [2405.18394]. This approach can be extended to scaled and local SIC variants to eliminate symmetry breaking across a broad class of systems.

## 5. Scaling Strategies and Computational Efficiency

Practical application to large or periodic systems, or for balancing overcorrection, has prompted the development of:
- **Orbitalwise (OSIC) and selective (SOSIC/vSOSIC) scaling:** Reduces SIC on core or many-electron-like orbitals, while applying full SIC on valence or frontier orbitals—restoring HOMO asymptotics and mitigating overcorrection in total energies and equilibrium properties [2005.09292][2308.04664].
- **Locally scaled SIC (LSIC):** Applies a local scaling function based on kinetic energy-density indicators (such as the iso-orbital indicator $z_\sigma(r)$), yielding restoration of correct uniform-gas limits and improved response properties (e.g. polarizabilities) [2103.05063].
- **WFL-SIC in periodic solids:** SIC is applied to compact Wannier-localized orbitals, minimizing the self-interaction energy with respect to the Wannier charge centers for efficient and accurate structural and band-gap predictions at $O(N)$ scaling [2203.08953].

Other computational efficiency innovations include real-space implementations, generalized Slater averages, and exploitation of analytic FOD gradients [2108.01186][1905.02631][1412.3101].

## 6. Representative Applications

Significant application domains and validation benchmarks include:
- **Transition-metal complexes:** Recovery of correct d/p character, opening of the HOMO–LUMO gap, and proper magnetic exchange coupling [2104.10327][2308.04664].
- **Water clusters and hydrogen-bonded systems:** Correction of overbinding in DFA and alignment with CCSD(T) reference energetics [2012.13469].
- **Polyacenes and conjugated π-systems:** Improved static polarizabilities and IPs, especially with local scaling (LSIC), eliminating both overestimation (DFA) and overcorrection (unscaled SIC) [2103.05063].
- **Reaction barriers:** Dramatic improvement in the accuracy of stretched-bond barrier heights, with the largest SIC correction traced universally to the most delocalized (stretched) orbital; combination with SCAN meta-GGA achieves state-of-the-art semi-local results [2601.11454].
- **Photoemission and orbital energies:** When FLOSIC is coupled with DCEP or OEP/KLI schemes, it yields physically meaningful orbital eigenvalues and improved fundamental gap predictions [2108.00606][2107.13631].

## 7. Limitations, Ongoing Methodological Refinements, and Future Directions

Key limitations and current directions include:
- **Overcorrection and loss of equilibrium accuracy:** Full SIC often overcorrects properties well-described by DFAs. Orbitalwise/local scaling approaches address this by limiting the SIC to where SIE is dominant [1909.13175][2005.09292][2103.05063].
- **Symmetry-breaking minima:** The prevalence of multiple local minima and symmetry breaking in delocalized or aromatic systems is addressed by ensemble corrections and FOD-initialization strategies [2405.18394][2109.08199].
- **Scaling and parallelization:** FLOSIC’s computational tractability is enhanced through the reduction from $O(N^2)$ to $O(N)$ parameters, nested SCF plus FOD loops, frozen-density acceleration, and efficient real-space and periodic implementations [1703.10742][2203.08953][2108.01186].
- **Extensions:** Ongoing developments include further tuning of local indicators for scaling (e.g. ELF, $z_\sigma$, $\beta$), hybrid-FLOSIC methods, excited-state and vibrational analyses via analytic forces, and robust FOD-optimization algorithms [2103.05063][1910.10677].

In sum, FLOSIC—by combining the formal rigor of unitary-invariant, size-extensive SIC with algorithmic efficiency and chemical interpretability—has established itself as a critical tool for overcoming key failures of standard DFAs, with broad impact across quantum chemistry and computational materials science [1703.10742][1910.10677][2104.10327][2203.08953].

---

**Key References:**  
[1703.10742] | [2104.10327] | [2308.04664] | [2203.08953] | [2103.05063] | [2012.13469] | [2109.08199] | [2107.13631] | [1909.13175] | [2005.09292] | [2405.18394] | [2601.11454] | [1412.3101] | [1905.02631]

Source: https://www.emergentmind.com/topics/fermi-lowdin-orbital-self-interaction-correction-flosic