One-body elements

One-body elements (obe for short) contains all of the one-body data that is typically converted from a density-functional theory (DFT) calculation or Wannier90 calculation. The one-body elements struct contains four components (obe_components) :

local_space defines the local correlated subspace to be solved within DMFT.
band_dispersion contains the band dispersion $\varepsilon(k, \nu, \sigma)$ where $\varepsilon$ is the energy at momentum k, band index $\sigma$, and spin $\sigma$.
downfolding_projector contains the projector which downfolds from Bloch space 𝓑 (indexed by $\nu$) to the correlated space (index by $m$).
ibz_symmetry_ops contains the operations to symmetrize observables over the full Brillouin zone.

To make it easy to build this object, ModEST provides factory functions that convert raw DFT data (from codes like VASP, Wien2k, or Wannier90) into fully prepared OBE objects ready for DMFT.

Embedding

The embedding object is the central controller of the DMFT problem: it defines how the correlated subspace ( $\mathcal{C}$) is partitioned into impurityies and handles the bookkeeping needed to map back and forth between the lattice and the impurities.

Typical use cases covered by embedding:

Equivalent atoms: use the same impurity model (i.e. same $n_{\mathrm{imp}}$) for different blocks $\alpha$.
Magnetic order: just reverse the spin for some atoms, $\tau(\sigma) = 1-\sigma$
Splitting impurity models : Depending on the block structure of the problem, the corresponding impurity model connected to atom $a$ can be split. Different impurity solvers can then be used for each of these split impurity models. The self-energies from both models map back to the a single atom $a$. Or a trivial solver is used to obtain zero self-energy (drop).
cluster DMFT: re-group a subset of atoms $a$ into a larger super atom which will be treated as a cluster DMFT problem.

The embedding object provides:

embedding_factories – convenient ways to construct common embeddings.
embedding_ops – operations to adapt and modify an embedding (e.g., split, drop, etc.).
embedding_methods – methods to map data between impurity and embed spaces.

This abstrction simplifies the DMFT loop, making it easy to construct arbitrary embedding scenarios.

Local Green's function

A central quantity of the DMFT self-consistency condition – the local Green's function defined as $$ [ G_{\mathrm{loc}}^{\sigma} ]_{m m'} = \sum_{\mathbf{k}} P_{m\nu}^{\sigma}(\mathbf{k}) \Big [ (\omega + \mu)\delta_{\nu\nu'} - H^{\sigma}_{\nu\nu'}(\mathbf{k}) - [P_{m\nu}^{\sigma}]^{\dagger}\Sigma_{\mathrm{embed}}P_{m'\nu'}^{\sigma}(\mathbf{k}) \Big ]^{-1} [P_{m'\nu'}^{\sigma}]^{\dagger},$$ where $\omega$ is a frequency (either real- or Matsubra), $\mu$ is the chemical potential, $H(\mathbf{k})$ is the one-body Hamiltonian, $P(\mathbf{k})$ are the projectors from the band to the orbital basis, and $\Sigma_{\mathrm{embed}}$ is the embedded self-energy.

ModEST computes this efficiently:

Woodbury: reduces the cost of matrix inversion from cubic in the number of bands to linear.
Adaptive Brillouin zone integration: for tight-binding models, allows you to specify desired integration accuracy, improving predictions of transport properties.

Hybridization function

The hybridization function $\Delta(\omega)$ is a key input to TRIQS impurity solvers. It describes the coupling of each impurity to the bath:

$$ [\Delta^{\sigma}(\omega)]_{mm'} = \omega - \tilde{\varepsilon} - G_{\mathrm{loc}}^{-1}(\omega) - \Sigma_{\mathrm{imp}}(\omega),$$

where $\tilde{\varepsilon} = \varepsilon - \mu - \Sigma_{\mathrm{DC}}$.

Densities and Chemical Potentials

Calculate the total electronic density from the lattice Green's function and search for the chemical potential $\mu$ that gives a target electorn count. We provide several algorithms optimized for speed and robusness. Due to the modularity you can also leverage exteranl root finding routines to write a custom chemical potential search.

Interaction Hamiltonians

To describe local interactions, ModEST provides tools to construct standard interaction Hamiltonians:

hubb_kan – Hubbard-Kanamori with density-density, spin-flip, and pair-hopping terms.
slater – full rotationally invariant Slater Hamiltonian.
den_den – density-density only Hamiltonian.

All built using TRIQS many-body operators.

This design makes it easier to swap and test different double counting schemes.

Post-processing

Compute spectral functions and related quantities:

$k$-resolved spectral functions $A(k, \omega)$.
orbital- and atom-resolved $k$-summed spectral function $A_{mm'}(\omega)$.

Plot-ready output to help you analyze and publish results.

Utilities

Extra tools that make life easier:

File I/O

No more boilerplated HDF5 code.

The checkpoint object saved DMFT calculation data (Green's functions, self-energies, etc.) in a clean, standardized HDF5 layout – automatically organized by iteration.

Root finders

Green's function symmetries

Analyze Green's functions analyze and symmetrize degenerate Green's function blocks.

Double counting

Double Counting DFT+DMFT requires subtracting a "double counting" term $\Sigma_{\mathrm{DC}}$ to avoid counting static correlations twice. ModEST abstracts this into a dedicated solver that computes $\Sigma_{\mathrm{DC}}$ based on dirrect formulas.

Table of Contents