---
title: Projectability-Disentangled Wannier Functions
url: https://www.emergentmind.com/topics/projectability-disentangled-wannier-functions-pdwf
type: topic
---

# Projectability-Disentangled Wannier Functions

Projectability-Disentangled Wannier Functions (PDWFs) are a class of maximally localized Wannier functions generated via an automated protocol that leverages atomic-orbital projectability measures rather than manual inner/outer energy windows. The PDWF methodology constructs Wannier functions that are both highly localized and chemically intuitive, enabling robust and efficient tight-binding models for a wide variety of materials—including metals, insulators, systems with magnetization, and those with strong spin–orbit coupling. PDWFs have been established as a central tool in large-scale electronic-structure workflows and high-throughput computational materials science, offering numerically precise (meV-level) interpolations of band structures with minimal manual intervention [2303.07877][2507.06840].

## 1. Mathematical Foundations and Projectability Measures

The construction of PDWFs is based on the use of a set of “atomic-like” trial orbitals, typically chosen as pseudopotential pseudo-atomic orbitals (PAOs) or hydrogenic atomic orbitals (AOs). For each Bloch state $|\psi_{m\mathbf{k}}\rangle$ (band index $m$ and wavevector $\mathbf{k}$), the core quantity is the projectability:

\[
p_{mn\mathbf{k}} = |\langle g_n | \psi_{m\mathbf{k}} \rangle|^2,
\quad
p_{m\mathbf{k}} = \sum_{n} p_{mn\mathbf{k}}
\]

where $\{|g_n\rangle\}$ are the atomic-like projectors. This total projectability $p_{m\mathbf{k}}$ quantifies the overlap of $|\psi_{m\mathbf{k}}\rangle$ with the projector manifold. For a complete set, $p_{m\mathbf{k}} \to 1$; for free-electron-like states, $p_{m\mathbf{k}} \ll 1$.

The mathematical framework underlying PDWFs closely follows the Marzari–Vanderbilt localization paradigm. The functional minimized is the quadratic spread

\[
\Omega = \sum_{n=1}^J \left[ \langle w_{n0} | \mathbf{r}^2 | w_{n0} \rangle - |\langle w_{n0} | \mathbf{r} | w_{n0} \rangle|^2 \right]
\]

where the Wannier functions $|w_{nR}\rangle$ are built from the Bloch manifold using semi-unitary matrices $U_{mn}(\mathbf{k})$.

## 2. Selection Criteria and Disentanglement Algorithm

PDWF subspace selection is governed by projectability-based decision thresholds:

- If $p_{m\mathbf{k}} < p_{\min}$: discard the state from the Wannier subspace.
- If $p_{m\mathbf{k}} \geq p_{\max}$: freeze the state (include it as-is).
- If $p_{\min} \leq p_{m\mathbf{k}} < p_{\max}$: include the state in the disentanglement optimization (“flexible”).

Typical default values are $p_{\max} \approx 0.95$–$0.99$ and $p_{\min} \approx 0.01$, but these are tunable to target specified interpolation errors. The actual Wannier subspace is constructed by (i) including all frozen bands, (ii) excluding discarded bands, and (iii) variationally building a smooth subspace from the “flexible” bands by minimizing the gauge-invariant part of the spread $\Omega_I$.

This process fully replaces the conventional practice of using inner and outer energy windows in the Souza–Marzari elucidation of entangled manifolds. The two-step minimization, first of $\Omega_I$ (disentanglement) and then the remainder ($\tilde{\Omega} = \Omega - \Omega_I$, for localization), yields compact and well-localized MLWFs [2303.07877][2507.06840].

## 3. Projector Manifold Augmentation and Automated Robustness

In a nontrivial fraction of cases (∼3–12% of materials, especially near $E_F+2$ eV above the Fermi level), the default set of PAOs may fail to fully span the relevant Hilbert space, leading to poor projectability and interpolation discontinuities. To address this, PDWFs incorporate an automated mechanism for extending the projector manifold by generating additional hydrogenic AOs:

- For each chemical element, minimal sets are defined (e.g., $1s$ for Period 1; $ns+np$ for Periods 2–3; $ns+(n-1)d$ for transition metals, etc.).
- Hydrogenic AOs are generated with radial parameter $\alpha$ fitted to PAO shapes and orthogonalized via Löwdin and Gram–Schmidt orthonormalization.
- This augmentation continues iteratively until all target states are sufficiently projectable (i.e., $P_{n\mathbf{k}} \geq P_{\min, \text{threshold}}$ at all relevant $n,\mathbf{k}$).

This procedural step achieves pseudopotential-agnostic, fully automatic Wannierization, and guarantees >99% success without manual tuning across large and diverse materials sets [2507.06840].

## 4. Spin–Orbit Coupling, Magnetization, and Advanced Protocols

The PDWF protocol has been extended to accommodate spin–orbit coupling (SOC) and arbitrary magnetization:

- For SOC, Bloch states are two-component spinors $\psi_{n\mathbf{k}}(r) = (\psi^{\uparrow}, \psi^{\downarrow})$, while projectors become $j$-dependent (total angular momentum eigenstates), typically derived from fully relativistic pseudopotentials (e.g., PseudoDojo or PAW from pslibrary).
- Magnetization is treated as collinear (separate spin-up and spin-down Wannierizations from spin-polarized DFT) or non-collinear (with SOC, via joint spinor analysis and using a dedicated MagneticStructureData plugin in AiiDA).
- Workflow modules (Quantum ESPRESSO, Wannier90, AiiDA-wannier90-workflows) accommodate these extensions, supporting parallel computation and robust automation.

PDWFs thus support accurate tight-binding models even for systems where exchange-correlation spin effects and relativistic interactions are pronounced [2507.06840].

## 5. Workflow Implementation and Automation

PDWFs are implemented as a fully automated protocol combining:

- DFT (Quantum ESPRESSO with suitable pseudopotentials for either non-relativistic or relativistic cases).
- Projector generation scripts that augment PAOs with hydrogenic AOs as needed.
- Automated overlap computation, projectability analysis, subspace selection, and spread functional minimization in the Wannier90 code.
- Post-processing for tight-binding Hamiltonian construction and band-structure interpolation.
- Comprehensive workflow management by AiiDA infrastructure, enabling k-point pool parallelization and end-to-end reproducibility [2303.07877][2507.06840].

Typical workflow steps are:

1. DFT: obtain $\varepsilon_{m\mathbf{k}}, |\psi_{m\mathbf{k}}\rangle$, and projectors $|g_n\rangle$.
2. Projector analysis and augmentation (if needed).
3. Compute all overlaps $A_{\alpha n \mathbf{k}} = \langle \phi_\alpha | \psi_{n\mathbf{k}} \rangle$ and $P_{n\mathbf{k}}$.
4. Apply PDWF selection and perform disentanglement/localization.
5. Construct Wannier functions and real-space Hamiltonian.

## 6. Performance Benchmarks and Comparative Analysis

Comprehensive benchmarks have established the superior accuracy and locality of PDWFs relative to other protocols:

| Method  | Mean $\eta_2$ (meV) | Max $\eta_2$ (meV) | MLWF Mean Spread (Å$^2$) | Hamiltonian Decay $\tau$ (Å) | Interpolation Success (%) |
|---------|---------------------|--------------------|--------------------------|------------------------------|--------------------------|
| SCDM    | 11.2                | 84                 | 3.5                      | 2.66                         | Not reported             |
| PDWF    | 4.2                 | 36.7               | 1.4                      | 2.27                         | $\sim$97                 |
| Robust PDWF (+hydrogenic, no SOC) | 2.00               | $<$15               | Not specified                | Not specified            | 100                      |
| Standard PDWF (SOC, PseudoDojo) | 10.14              | Not specified        | Not specified                | Not specified            | 87.5                     |
| Robust PDWF (+hydrogenic, SOC) | $\approx$2.1           | $<$15               | Not specified                | Not specified            | 100                      |

Key metrics used include:
- Band-distance $\eta_2$ (root mean squared error up to $E_F+2$ eV),
- Wannier function quadratic spread $\Omega$,
- Real-space Hamiltonian decay length $\tau$ via $\|H(R)\| \propto e^{-\|R\|/\tau}$,
- Wannier center positions relative to atoms.

In high-throughput tests across tens of thousands of materials, robust PDWF protocols (with hydrogenic projector augmentation) achieved $\eta_2 < 15$ meV and $\gtrsim$100% success without manual intervention [2507.06840][2303.07877].

## 7. Applications, Case Studies, and Best Practices

PDWFs have been applied to diverse systems and phenomena:

- **Graphene**: Avoids free-electron contamination in $\sigma/\pi$ bands via customized $p_{\max}$/$p_{\min}$, yielding $p_z$, $sp^2$-like Wannier functions with meV-level interpolation.
- **Silicon**: Conduction bands with $3d$ character are accurately Wannierized via either d-type PAO inclusion or joint projectability/disentanglement (PD+ED).
- **Metals (Cu, SrVO$_3$)**: Achieves atomic $d$- and $s$-like MLWFs with $<$5 meV interpolation errors; more localized and less k-point-intensive than SCDM Wannierization.
- **High-throughput materials screening**: Applied to 21,737 nonmagnetic materials, with >97% achieving $\eta_2 < 20$ meV automatically [2303.07877].

**Best practices** include use of high-quality PAOs as projectors, initializing with $p_{\max}=0.95$/$p_{\min}=0.01$ and a frozen energy window $[E_F, E_F+2 {\rm~eV}]$, automating threshold refinement to minimize $\eta_\nu$, and validating output MLWF centers and spreads. For robust operation, the “projector generation” step for hydrogenic AOs is recommended wherever low projectability arises, enabling consistent operation across different pseudopotential sets and aberrant chemical environments [2507.06840].

PDWFs have thus enabled automated, chemically meaningful, and highly accurate tight-binding models suitable for both detailed single-system analyses and large-scale materials informatics.

Source: https://www.emergentmind.com/topics/projectability-disentangled-wannier-functions-pdwf