Programs and Methods Based on Wannier Functions

The tools listed below fall into several categories: (i) programs directly interfaced with Wannier90, (ii) programs that read standard Wannier90 tight-binding files such as *_hr.dat, (iii) packages that use Wannier-derived model Hamiltonians after conversion or custom interfacing, and (iv) methods based on the broader Wannier representation but not necessarily on the Wannier90 code itself.

Wannier functions are used far beyond the Wannier90 code itself. The Hamiltonian generated by Wannier90 is not only useful inside Wannier90, but is also widely used as an input for external programs that perform more advanced calculations. In this sense, Wannier90 often serves as a bridge between first-principles electronic-structure calculations and specialized simulation tools.

Common Wannier90 output files that can be generated for downstream applications include seedname_hr.dat, seedname_centres.xyz, and seedname_tb.dat, contain the real-space tight-binding Hamiltonian, Wannier centers, and related matrix elements. These files can be used by many post-processing programs for calculations of electronic transport, topology, electron–phonon coupling, optical properties, many-body effects, and other physical phenomena.

Wannier90 v4 and library workflows

Since Wannier90 v4, most functionality of the standalone wannier90.x executable can also be accessed through a parallel library interface, enabling in-memory integration with electronic-structure and workflow codes. See Wannier90 as a Library and Software Infrastructure for details.

This page has two related goals. First, it provides examples of external programs that interface directly with Wannier90 or use data generated by it. Second, it introduces alternative methods for constructing Wannier functions or related localized orbitals without using Wannier90.

Programs and Methods Using Wannier Representations

The programs in this section either interface directly with Wannier90 or use Wannier90 output files, such as the real-space Hamiltonian seedname_hr.dat.

Wannier90 is normally used to transform Bloch states obtained from a first-principles calculation into a compact localized basis. External programs then use this basis to calculate physical properties that may not be available directly in the original DFT code.

These applications include:

• electronic transport;
• topological properties;
• electron–phonon interactions;
• optical and nonlinear optical response;
• dynamical mean-field theory;
• excitonic calculations;
• model Hamiltonians and many-body simulations.

Topological Materials and Surface States

Several programs use Wannier Hamiltonians to study topological properties of materials:

• WannierTools — surface band structures, surface density of states, Fermi arcs, Weyl points, Dirac points, nodal lines, Berry curvature, Wilson loops, Chern numbers, and $\mathbb{Z}_2$ topological invariants.
• Z2Pack — calculation of $\mathbb{Z}_2$ invariants and topological indices from Wannier-based tight-binding models.
• TBmodels — construction, manipulation, and symmetrization of tight-binding models derived from Wannier90.
• PythTB — simple Python-based tight-binding modeling, including models imported from Wannier90.
• WannierBerri — Berry curvature, anomalous Hall conductivity, orbital magnetization, nonlinear Hall response, and other Berry-related quantities.

Electron–Phonon Coupling and Superconductivity

Wannier functions are also essential for interpolating electron–phonon matrix elements:

• EPW — electron–phonon coupling, phonon-limited carrier mobility, superconductivity, Eliashberg spectral functions, and superconducting critical temperatures.
• Perturbo — electron–phonon interactions, carrier dynamics, ultrafast relaxation, and transport properties.
• EPWpy — a Python-based workflow and automation layer for EPW calculations, providing tools and Jupyter workflows for preparing, running, and analyzing electron–phonon calculations based on Quantum ESPRESSO, Wannier90, and EPW.

Phonons and Topological Lattice Dynamics

• Simphony — a tight-binding package for lattice vibrations and topological phonon analysis. It constructs phonon tight-binding models from real-space dynamical matrices expressed in the Wannier90 *_hr.dat format and supports bulk and slab phonon dispersions, surface spectral functions, Wilson loops/Wannier charge centers, Berry phases and curvature, Weyl-node detection and chirality, and LO–TO splitting in polar materials. Simphony employs a Wannier-like localized representation of lattice-vibration degrees of freedom rather than electronic Bloch states. Its QE2TBDAT utility can generate the required wannier90_hr.dat directly from Quantum ESPRESSO ph.x calculations, so running Wannier90 itself is not required.

Transport Calculations

The following packages either consume Wannier-derived tight-binding models directly or can be coupled to them through conversion or custom workflows.

• BoltzTraP2 — semiclassical Boltzmann transport, electrical conductivity, Seebeck coefficient, thermal conductivity, and thermoelectric properties.
• Kwant — quantum transport in tight-binding systems, nanostructures, and devices.
• NanoDCAL — quantum transport and device simulations based on localized basis representations.
• OpenMX post-processing tools — transport and tight-binding analysis using localized orbitals and Wannier-like Hamiltonians.

Optical, Nonlinear Optical, and Quantum-Geometric Response

Wannier interpolation is very useful for optical calculations because it allows dense $\mathbf{k}$-space sampling at low computational cost:

Berry/Wannier interpolation

• WannierBerri — optical conductivity, Berry curvature dipoles, anomalous Hall effects, nonlinear Hall effects, and related geometric response functions.

Excitonic/BSE optics

• Yambo-based workflows — many-body optical spectra, GW, BSE, and excitonic effects can be combined with Wannier-based interpolation in advanced workflows.
• WanTiBexos, WTB,  WannierBSE, and Xatu use Wannier-based tight-binding representations as an efficient starting point for Bethe–Salpeter equation calculations, enabling studies of exciton energies and wavefunctions, optical spectra, exciton dispersion and transport, and excitonic quantum geometry.
• PyMEX is specialized for excitonic calculations in moiré materials, where Wannier-based models provide an efficient description of large supercells and spatially varying excitonic states.

Nonlinear and strong-field optics

• WannierNLQG — a Julia package for nonlinear optical responses and quantum-geometric quantities from Wannier tight-binding models, including shift current, ordinary/spin/photon-drag injection currents, Berry curvature, quantum metric, geometric multipoles, and degeneracy-aware formulations based on projectors and Wilson loops.

• Semiconductor Wannier Equations (SWE) — a real-time, real-space formulation of ultrafast light–matter dynamics in crystals using a localized Wannier basis. The method provides a gauge-clean alternative to reciprocal-space semiconductor Bloch equations and includes electron–electron interactions at the time-dependent Hartree plus static screened-exchange level, with applications to nonlinear optics, high-harmonic generation, and attosecond spectroscopy. SWE is currently a theoretical formalism rather than a standalone software package.

Magnetic Interactions, Spin Symmetry, and Spin Models

Wannier Hamiltonians can also be used to derive effective magnetic models:

• TB2J — exchange interactions, Dzyaloshinskii–Moriya interactions, magnetic anisotropy, and construction of spin Hamiltonians from Wannier-based tight-binding models.
• WannierBerri — spin Berry curvature, spin Hall conductivity, orbital magnetization, and related spin-orbit effects.
• IRSSG — an open-source package for identifying spin space groups and analyzing magnetic band representations. It can operate directly on Wannier90-format *_hr.dat Hamiltonians and determine spin-space symmetry, little groups, and irreducible corepresentations, making it particularly useful for magnetic materials and altermagnets.

Moiré Materials and Large-Scale Tight-Binding Models

Wannier-based Hamiltonians are especially useful when the system is too large for direct first-principles calculations:

• MoireStudio, TWISTER — construction of twisted bilayers and moiré supercells combined with Wannier tight-binding models.
• PythTB — educational and research-level tight-binding models for large supercells.
• TBmodels — construction and manipulation of large tight-binding Hamiltonians.
• KITE — large-scale tight-binding simulations, density of states, conductivity, and quantum transport in large systems.
• Kwant — quantum transport in finite systems, ribbons, defects, and nanodevices.
• PyMEX - solving the Bethe–Salpeter equation in moiré exciton systems.

Many-Body and Model-Hamiltonian Workflows Using Wannier-Derived Parameters

Wannier Hamiltonians are often used as a starting point for correlated-electron calculations:

• TRIQS/DFTTools — DFT+DMFT calculations using localized correlated subspaces.
• DMFTwDFT3 — DFT+DMFT calculations of strongly correlated materials using orbital-dependent many-body self-energies. 
• ComDMFT — charge self-consistent DFT+DMFT workflows.
• ALPS-based workflows — model Hamiltonian simulations derived from Wannier tight-binding parameters.
• Custom Hubbard-model and exact-diagonalization codes (TRIQS, divERGe, HPhi, PythTB) — construction of effective models using Wannier hopping parameters.

Visualization and Analysis Tools

Wannier90 output can also be used for visualization and post-processing:

• FermiSurfer — visualization of Fermi surfaces and momentum-dependent quantities.
• XCrySDen — visualization of Wannier functions and real-space orbitals.
• VESTA — visualization of Wannier functions, charge densities, and orbital shapes.
• ParaView — visualization of large real-space volumetric data.
• Python/NumPy/SciPy workflows — custom analysis of Wannier Hamiltonians, band structures, density of states, Berry curvature, and real-space hopping amplitudes.

Ab Initio Codes with Wannier90 Interfaces or Wannier Workflows

Many first-principles electronic-structure codes provide direct Wannier90 interfaces, conversion tools, or compatible Wannier workflows:

• Quantum ESPRESSO  
• VASP  
• ABINIT  
• SIESTA  
• WIEN2k  
• FLEUR  
• OpenMX  
• GPAW  
• Elk  
• CASTEP  
• CP2K  
• BigDFT  

These programs provide Bloch wave functions, band energies, and overlap matrices, which Wannier90 then transforms into localized Wannier functions and compact tight-binding Hamiltonians.

In summary, Wannier90 is not only a program for constructing Wannier functions. It is a central component of a much larger computational ecosystem. Once a reliable Wannier Hamiltonian is obtained, it can be used to study band structures, topology, transport, optical response, superconductivity, electron–phonon coupling, magnetism, moiré materials, excitons, and correlated-electron physics.

Automated Wannierization and Meta-Optimization

Recent developments increasingly treat Wannierization itself as an optimization problem, with automated selection of projections, disentanglement windows, and other hyperparameters.

• metaWann — automated meta-optimization of Wannierization parameters using differential evolution and Bayesian optimization. It optimizes projections, disentanglement windows and related hyperparameters in order to obtain high-accuracy Wannier interpolations with minimal user intervention.

Independent Methods for Constructing Wannier Functions

Wannier90 is one of the most widely used implementations of Wannier functions, but it is not the only possible approach. Several electronic-structure packages contain their own algorithms for constructing Wannier functions or related localized orbitals.

These implementations can be divided into three broad groups.

1. Iterative Construction of Maximally Localized Wannier Functions

In this approach, a unitary transformation among the Bloch states is optimized by minimizing the spatial spread of the resulting Wannier functions. For entangled bands, an additional disentanglement procedure is used to identify an optimal low-energy subspace.

This is conceptually closest to the standard method implemented in Wannier90.

Examples include:

• RESPACK — construction of low-energy effective Hamiltonians, screened Coulomb interactions, and Hubbard parameters. RESPACK contains its own independent implementation of disentanglement and spread minimization and does not use Wannier90. The resulting Wannier basis is used to calculate hopping parameters, screened Coulomb interactions, and low-energy effective Hamiltonians.
• JDFTx — contains an internal Wannier module that constructs localized functions and treats both isolated and entangled groups of bands.
• OpenMX — provides an internal implementation of maximally localized Wannier functions, while also offering a separate interface to the external Wannier90 program.

2. Direct, Algebraic, and Non-Variational Wannier Constructions

These methods construct a localized basis directly from projections, singular-value decompositions, density matrices, or selected localized trial orbitals. They may avoid the iterative spread-minimization procedure used for MLWFs.

Such approaches are often faster and more automatic, although the resulting functions are not necessarily maximally localized according to the Marzari–Vanderbilt spread functional.

Examples include:

• VASP LOCPROJ — constructs Wannier orbitals from projections onto localized orbitals using a single-step singular-value decomposition.
• VASP SCDM — uses selected columns of the density matrix to automatically construct a Wannier basis.
• OpenMX Closest Wannier Functions (CWFs) — constructs Wannier functions that remain as close as possible to physically selected pseudo-atomic, hybrid, or molecular guiding orbitals while also allowing disentanglement.
• Discrete Adiabatic Transport — a deterministic, non-variational construction of maximally localized Wannier functions that combines gauge smoothing, projected-position operators, and adiabatic transport across band degeneracies. It provides an alternative to conventional iterative spread-functional minimization.

3. Localization of Occupied or Partially Occupied Orbitals

A third group of methods localizes occupied or selected partially occupied orbitals using criteria such as Boys localization, Berry-phase localization, Pipek–Mezey localization, or related spread measures.

These localized orbitals are often used to analyze chemical bonding, polarization, molecular dynamics, charge transfer, and local electronic structure. They are closely related to Wannier functions, although they are not always intended for band interpolation or the construction of a complete low-energy tight-binding Hamiltonian.

Examples include:

• CP2K — contains internal localization methods for occupied and partially occupied orbitals, including Wannier-type localization and gradient-ascent localization for partially occupied states.
• ASE/GPAW — provides an independent Wannier localization procedure that can construct localized orbitals from GPAW wave functions.

CP2K also provides a separate interface to Wannier90, so its internal orbital-localization methods and its Wannier90 workflow should be regarded as distinct options.

These three groups should not be viewed as completely separate definitions of Wannier functions. All of them construct localized combinations of extended electronic states, but they differ in the localization criterion, treatment of entangled bands, level of automation, and intended physical application.

Wannier90 remains particularly important because it provides a general and transferable framework for constructing maximally localized Wannier functions and exporting the resulting Hamiltonians and matrix elements to a large ecosystem of external simulation programs.

Next: Wannier90 as a Library and Software Infrastructure