---
title: Simple Lehmann Model Overview
url: https://www.emergentmind.com/topics/simple-lehmann-model
type: topic
---

# Simple Lehmann Model Overview

The expression **“Simple Lehmann Model”** is used in several technically distinct but structurally related ways across recent literature. In numerical many-body theory, it commonly denotes a compact **discrete Lehmann representation** of imaginary-time Green’s functions, in which a spectral integral is replaced by a small sum of exponentials in imaginary time and simple poles in Matsubara frequency [2107.13094] [2110.06765] [2404.02334]. In quantum field theory, the same phrase also appears in connection with standard or generalized **Källén–Lehmann** decompositions of two-point functions in local, nonlocal, curved, lattice, and Lorentz-violating settings [2405.14056] [2306.00090]. In nonequilibrium integrable systems, it further refers to explicit **Lehmann sums** for post-quench observables and their Quench Action reductions [2605.20065]. This suggests a common conceptual core: a correlator is represented as a superposition over spectral data, poles, or intermediate states, with the model becoming “simple” when that superposition is reduced to a compact, computationally tractable form.

## 1. Common structure and scope

The usages appearing in the literature share the same formal motif: an interacting object is expressed as a sum or integral over simpler building blocks.

| Context | Object | Representative form |
|---|---|---|
| Imaginary-time numerics | Green’s function | $G(\tau)\approx \sum_{l=1}^r K(\tau,\omega_l)\,g_l$ |
| Spectral QFT | Two-point function | $\int d\mu^2\,\rho(\mu^2)/(p^2-\mu^2+i\epsilon)$ |
| Quench dynamics | One-point function | $\sum_{m,n} c_m^* c_n\,e^{i(E_m-E_n)t}\,\langle m|O|n\rangle$ |

In the DLR setting, the “simple” model is explicitly discrete: the continuous spectral density is replaced by an effective spectral density that is a sum of $\delta$ functions,
$$
A_{\mathrm{eff}}(\omega)=\sum_{k=1}^r w_k\,\delta(\omega-s_k),
$$
so that both the imaginary-time and Matsubara-frequency representations become finite sums [2107.13094]. In standard Källén–Lehmann theory, the same simplifying idea appears as a decomposition over invariant masses with a positive spectral density, often split into a one-particle pole plus continuum. In nonequilibrium quench problems, the simplification is not necessarily a pole model, but rather a reduction of the full double Lehmann sum to a more manageable sampling problem or to a Quench Action single-sum expansion [2605.20065]. A further terminological complication is historical: in low-energy hadronic physics, “Lehmann model” can mean H. Lehmann’s 1972 unitarized chiral approach, which is explicitly distinguished from the Lehmann spectral representation and from LSZ [1002.2519].

## 2. Discrete Lehmann representation in imaginary time

In the DLR formulation, imaginary-time Green’s functions are built from the finite-temperature Lehmann representation. In scaled variables, with $\tau\in[0,1]$ and $\omega\in[-\Lambda,\Lambda]$,
$$
G(\tau) = -\int_{-\Lambda}^{\Lambda} K(\tau,\omega)\,\rho(\omega)\,d\omega,
\qquad
K(\tau,\omega)=\frac{e^{-\tau\omega}}{1+e^{-\omega}}.
$$
For numerical stability at large negative $\omega$, the equivalent form
$$
K(\tau,\omega)=\frac{e^{(1-\tau)\omega}}{1+e^{\omega}}
$$
is also used [2110.06765].

The core approximation is a low-rank decomposition of the kernel,
$$
K(\tau,\omega)\approx \sum_{l=1}^r K(\tau,\omega_l)\,\pi_l(\omega),
$$
which yields the DLR expansion
$$
G(\tau)\approx \sum_{l=1}^r K(\tau,\omega_l)\,g_l.
$$
The same coefficients reconstruct the Matsubara function through
$$
G(i\nu_n)=\sum_{l=1}^r K(i\nu_n,\omega_l)\,g_l
=\sum_{l=1}^r \frac{g_l}{\omega_l+i\nu_n}.
$$
For fermions,
$$
i\nu_n=i(2n+1)\pi/\beta,
$$
and for bosons,
$$
i\nu_n=i(2n)\pi/\beta.
$$
The basis therefore consists of exponentials in imaginary time and simple poles in Matsubara frequency [2110.06765] [2404.02334].

The numerical efficiency of the representation follows from the low numerical rank of the analytic continuation kernel. The observed scaling is
$$
r=O\big(\log\Lambda\,\log(1/\epsilon)\big),
$$
with $\Lambda=\beta\omega_{\max}$ and user-controlled tolerance $\epsilon$ [2107.13094] [2110.06765]. The literature contrasts this with other imaginary-time parameterizations: uniform grid/Fourier series require $O(\Lambda/\epsilon)$ degrees of freedom, orthogonal polynomials such as Legendre and Chebyshev require $O(\sqrt{\Lambda}\,\log(1/\epsilon))$, while IR and DLR both exhibit $O(\log\Lambda\,\log(1/\epsilon))$ scaling [2110.06765]. This is the sense in which the DLR-based simple Lehmann model is both compact and explicit.

## 3. Basis construction, interpolation grids, and software realizations

The DLR basis is constructed by discretizing the kernel on a fine adaptive grid and applying interpolative decomposition through pivoted QR. The procedure first forms
$$
A_{ij}=K(\tau_i^f,\omega_j^f),
$$
then applies a rank-revealing column-pivoted QR to select the DLR frequencies $\{\omega_l\}_{l=1}^r$, and finally a rank-revealing row-pivoted QR to $K(\tau_i^f,\omega_l)$ to select the DLR imaginary-time nodes $\{\tau_k\}_{k=1}^r$ [2110.06765]. The selected time nodes are precisely those from which the coefficients $\{g_l\}$ can be stably recovered by solving a small square linear system. A corresponding DLR Matsubara grid $\{i\nu_{n_k}\}_{k=1}^r$ is built by increasing a cutoff $n_{\max}$, practically set to $n_{\max}\approx\Lambda$, until the selected nodes stabilize [2110.06765].

The implementation ecosystem is unusually explicit. **libdlr** provides `dlr_it_build` to build the basis and imaginary-time grid, `dlr_mf` to construct the Matsubara grid, `dlr_it2cf_init` and `dlr_it2cf` to recover coefficients from imaginary-time grid samples, `dlr_it_fit` for least-squares fitting from noisy data, `dlr_it_eval` and `dlr_mf_eval` for evaluation, `eqpts_rel`, `abs2rel`, and `rel2abs` for relative and absolute time-point conversion, and `dlr_it2itr_init` for the reflection mapping used by the SYK self-energy [2110.06765]. Standard operations such as convolution and integration are available in the explicit exponential basis, and the Dyson equation can be solved either in Matsubara frequency through
$$
G^{-1}(i\nu_n)=G_0^{-1}(i\nu_n)-\Sigma(i\nu_n),
\qquad
G_0(i\nu_n)=(i\nu_n+\mu)^{-1},
$$
or in imaginary time via DLR discretization of convolutions [2110.06765].

The software landscape spans multiple languages. **libdlr** is written in Fortran, provides a C header interface, contains the Python module **pydlr**, and is accompanied by the stand-alone Julia implementation **Lehmann.jl** [2110.06765]. **cppdlr** provides a C++ interface implementing the same DLR logic for functions in imaginary time and Matsubara frequency, with explicit exponentials in time and simple poles in Matsubara frequency, and has been integrated into TRIQS and used in DMFT solvers, Keldysh mixing Green’s function calculations, and diagrammatic evaluations [2404.02334].

The practical workflow is standardized. One chooses $\beta$, estimates $\omega_{\max}$, sets $\Lambda=\beta\omega_{\max}$, chooses $\epsilon$, builds the basis and grids, samples or fits the function, and evaluates the resulting representation in either domain. The reported examples are concrete: for $\Lambda=100$ and $\epsilon=10^{-6}$, there are $r=21$ nodes in each domain; for noisy data with $\beta=1000$, $\Lambda=1000$, noise magnitude $\eta=10^{-2}$, and $\epsilon=10^{-2}$, the DLR rank is $r=13$ and the fit agrees well with $G(\tau)$; and for an SYK solver with $\beta=1000$, $\Lambda=5000$, and $\epsilon=10^{-10}$, there are $r=66$ DLR nodes in each domain [2110.06765].

## 4. Sampled, fitted, and learned Lehmann models in many-body dynamics

Outside equilibrium imaginary-time calculations, the simple Lehmann model becomes a tool for explicit state sums. For a homogeneous quantum quench from an initial pure state $|\Psi_0\rangle=\sum_n c_n |n\rangle$, the time-dependent expectation value of a local operator is
$$
\langle O(t)\rangle
=\sum_{m,n} c_m^* c_n\,e^{i(E_m-E_n)t}\,\langle m|O|n\rangle.
$$
This direct-sum representation is exact but exponentially hard to evaluate in interacting systems because the number of relevant eigenstates grows exponentially with system size, while overlaps and form factors can be exponentially small and strongly fluctuating [2605.20065].

The Quench Action formalism reduces the double sum to a single sum over excitations of a representative saddle-point state,
$$
\langle O(t)\rangle \xrightarrow[L\to\infty]{}
\sum_\mu \Big(\alpha_\mu(t)F_\mu^{O}+\alpha_\mu^*(t)\,[F_\mu^{O^\dagger}]^*\Big),
$$
with
$$
F_\mu^{O}=\frac{\langle\Psi_0|\mu\rangle\langle\mu|O|\lambda_{\rm sp}\rangle}{2\,\langle\Psi_0|\lambda_{\rm sp}\rangle},
\qquad
\alpha_\mu(t)=e^{i(E_\mu-E_{\lambda_{\rm sp}})t}.
$$
The Monte Carlo scheme samples configurations with probability
$$
\mathcal{P}(\boldsymbol{\lambda})=\frac{|F_{\boldsymbol{\lambda}}|}{Z},
$$
using a Metropolis–Hastings acceptance step. The method is benchmarked in the transverse-field Ising chain and the Lieb–Liniger model, and the paper identifies a sign problem for more general dynamical correlators and generic initial states, while integrable initial states with pair structures avoid it [2605.20065].

A different computational development is the Lehmann-representation-based PINN for Anderson impurity models. There the self-energy is parameterized by a nonnegative auxiliary spectral function,
$$
\Sigma(z)=\Sigma_\infty+\int_{-\infty}^{\infty} d\omega\,\frac{\rho_{\rm aux}(\omega)}{z-\omega},
$$
which is implemented as a finite pole sum,
$$
\Sigma(z)=\Sigma_\infty+\sum_{\ell=1}^{L}\frac{w_\ell}{z-p_\ell},
\qquad
w_\ell\ge0.
$$
This hard-wires causality, analyticity, Kramers–Kronig consistency, and the correct $1/\omega_n$ tail into the model, while sparse IR/DLR grids stabilize the representation [2411.18835]. In the reported single-orbital AIM study, the model is trained for $U\in[0.1,4.0]$, $V\in[0.1,1.0]$, and $\beta=100$, and the Lehmann constraints reduce the maximum test error in electron filling by a factor of about $7.8$ relative to a comparable network without the Lehmann projection [2411.18835]. A plausible implication is that the “simple Lehmann model” increasingly functions not only as a compressed representation, but also as an inductive bias in machine-learned many-body surrogates.

## 5. Källén–Lehmann representations in local, nonlocal, curved, lattice, and Lorentz-violating field theory

In local quantum field theory, the canonical scalar Källén–Lehmann representation is
$$
G_F(p^2)=\int_{0}^{\infty} d\mu^2\;\frac{\rho(\mu^2)}{p^2-\mu^2+i\epsilon},
$$
with $\rho(\mu^2)\ge0$, and typically
$$
\rho(\mu^2)=Z\,\delta(\mu^2-m^2)+\theta(\mu^2-\mu_0^2)\,\sigma(\mu^2),
\qquad 0\le Z\le1.
$$
The time-ordered two-point function is correspondingly a superposition of free propagators weighted by the same positive spectral density [2405.14056]. This basic structure extends, with important modifications, to several nonstandard settings.

In **nonlocal quantum gravity** with entire form factors,
$$
\mathcal{F}(\Box)=e^{H(\Box)},
$$
the time-ordered two-point function still admits a generalized Källén–Lehmann representation with the standard momentum dependence
$$
\tilde G_{\rm to}(p^2)=\int_{0}^{\infty} ds\;\frac{\rho(s)}{p^2-s+i\epsilon},
$$
while nonlocality enters through the spectral density. The free spectral density coincides with the local one because $H(0)=0$, the physical spectrum is unchanged, no extra poles or ghosts appear, and the local limit is recovered smoothly as $\ell_*\to0$ [2405.14056].

In **de Sitter spacetime**, the KL decomposition is organized by unitary irreducible representations of $SO(d+1,1)$. For scalars,
$$
G_{\mathcal{O}}(Z)
=
\int_{\mathbb R} d\lambda\; \rho^{\mathcal{P},0}_{\mathcal{O}}(\lambda)\, G_{\lambda,0}(Z)
+
\int_{-\frac d2}^{\frac d2} d\lambda\; \rho^{\mathcal{C},0}_{\mathcal{O}}(\lambda)\, G_{i\lambda,0}(Z),
$$
with nonnegative spectral densities in the Bunch–Davies vacuum; for spinning traceless symmetric operators, the principal-series part sums over spins $\ell=0,\dots,J$ and is supplemented by complementary-series contributions, and in $dS_2$ also by discrete-series terms [2306.00090]. A related development gives compact formulas for spinor–scalar and spinor–spinor loop densities in AdS and dS, derives corresponding Källén–Lehmann decompositions of Wightman products, and writes chain-approximation spectral equations for scalar and spinor one-loop self-energies in Yukawa theory in de Sitter space [2508.07467].

In **Lorentz-violating field theory** with a constant symmetric observer-Lorentz tensor $c^{\mu\nu}$, the exact propagators still admit generalized Källén–Lehmann representations, but the spectral densities depend not only on $p^2$ but also on the observer scalars
$$
I_n\equiv p\cdot c^n\cdot p,\qquad n=1,2,3.
$$
For fermions, the spectral matrix involves the four Dirac structures $1$ and $J_k=\gamma\cdot c^k\cdot p$ for $k=0,1,2,3$, and equal-time commutator and anticommutator sum rules generalize the usual normalization conditions [1112.5739].

In **lattice field theory**, the Umezawa–Kamefuchi–Källén–Lehmann representation follows from Hermiticity, translational invariance, reflection positivity, and polynomial boundedness. In momentum space one obtains
$$
G(p)=\int_{0}^{\infty} d\mu^2\,\frac{\rho(\mu^2)}{\hat p^2+\mu^2},
\qquad
\rho(\mu^2)\ge0,
$$
with $\hat p^2=\sum_\mu 4\sin^2(p_\mu/2)$. The positivity of the spectral density is thus necessary under the stated axioms, and the overlap scalar boson provides a counterexample: its continuum density changes sign, so the model violates reflection positivity [1201.3415].

Across these settings, the common point is not a fixed formula but a conserved architecture: two-point functions are decomposed into simpler free or harmonic building blocks, with positivity, analyticity, and pole structure encoded in the spectral weights.

## 6. Historical usage and terminological ambiguity

A final usage is historically important because it changes the meaning of “Lehmann model” altogether. In the review of low-energy soft hadronic physics, the term refers to **H. Lehmann’s 1972 unitarized chiral approach** to low-energy $\pi\pi$ scattering. There, the one-loop chiral amplitude in the chiral limit is combined with exact elastic unitarity and analyticity through an effective-range expansion for $\cot\delta$, the inverse-amplitude method, or Padé approximants. The relevant formulas are
$$
f(\nu)=\frac{e^{i\delta(\nu)}\sin\delta(\nu)}{\rho(\nu)},
\qquad
g(\nu)\equiv f^{-1}(\nu)=\rho(\nu)\,(-i+\cot\delta(\nu)),
$$
together with the effective-range form
$$
\sqrt{\nu}\,\cot\delta(\nu)=g(0)+\alpha_0\nu+\alpha_1\nu^2+\dots
$$
and the Padé unitarization
$$
f^{[0,1]}(\nu)=\frac{f^0}{1-\frac{f^1(\nu)}{f^0}}.
$$
The review states explicitly that this is **not** the Lehmann spectral representation and **not** the LSZ formalism [1002.2519].

This terminological divergence matters because the same phrase can denote either a **spectral representation of correlators** or a **unitarized low-energy scattering model**. The literature therefore supports a careful, context-dependent reading. In numerical many-body work, “simple Lehmann model” usually means a low-rank discrete spectral ansatz. In field theory, it more often points to Källén–Lehmann spectral decompositions and their generalizations. In older hadronic applications, it denotes a unitarized chiral partial-wave construction. The unifying thread is the use of Lehmann-type structure to compress, organize, or reconstruct nontrivial physics from a minimal spectral or state-sum description.

Source: https://www.emergentmind.com/topics/simple-lehmann-model