---
title: 'QAVG-DMFT: Hybrid Quantum-Classical Scheme'
url: https://www.emergentmind.com/topics/qavg-dmft-scheme
type: topic
---

# QAVG-DMFT: Hybrid Quantum-Classical Scheme

QAVG-DMFT is a quantum–classical hybrid scheme for dynamical mean-field theory at finite temperature in which the impurity one-particle Green’s function is reconstructed from quantum phase estimation data. In the formulation explicitly named in the literature, the acronym refers to “QPE averaged over variable grids”: modified finite-temperature, channel-agnostic QPE circuits extract excitation-energy differences, while a classical reconstruction layer infers a compact spectral representation by fitting QPE histograms collected under multiple grid settings [2605.29681]. Within DMFT, the reconstructed impurity Green’s function is used to obtain the impurity self-energy and close the self-consistency loop. The term is not used uniformly across the broader DMFT literature: several earlier papers explicitly state that “QAVG-DMFT” is not their terminology and only discuss related GW+DMFT, QSGW+DMFT, or Dual Boson constructions [1602.09129] [1411.5180] [1611.07090] [1109.3972] [1703.08446].

## 1. Definition and terminological scope

In the explicit sense introduced for finite-temperature quantum impurity simulation, QAVG-DMFT combines two components: a quantum subroutine based on modified QPE circuits for the impurity Green’s function, and a classical post-processing stage called QAVG, “QPE averaged over variable grids,” that reconstructs the Green’s function from measured probability distributions [2605.29681]. The scheme is designed for DMFT, where the impurity problem is the computational bottleneck and where direct access to real-frequency spectral information is valuable because analytic continuation from imaginary time is ill-posed and unstable [2605.29681].

The term is not standardized across all uses of DMFT-related abbreviations. One paper on double occupancy in DMFT and the Dual Boson approach states that “QAVG-DMFT” does not appear in that work and contrasts a principled self-consistent Dual Boson construction with any naive averaging of impurity- and lattice-based estimates [1602.09129]. Several GW+DMFT papers likewise state that they do not use the term, and instead map a possible “QAVG-DMFT” interpretation to quasiparticle-averaged or quasiparticle self-consistent GW+DMFT variants [1411.5180] [1611.07090] [1109.3972] [1703.08446]. This suggests that the most precise encyclopedia usage is to reserve “QAVG-DMFT” for the channel-agnostic finite-temperature QPE plus variable-grid reconstruction scheme of [2605.29681], while separately noting the distinct quasiparticle-based reinterpretations proposed elsewhere.

## 2. Finite-temperature Green’s functions and channel-agnostic QPE

The finite-temperature formulation starts from standard one-particle Green’s functions. For fermionic operators $c_\alpha$ and $c_\beta$, the retarded and Matsubara Green’s functions are written as
$$
G^R_{\alpha\beta}(\omega) = -i \int_0^{\infty} dt\, e^{i\omega t} \langle\{ c_\alpha(t), c_\beta^\dagger(0) \}\rangle_\beta,
$$
and
$$
G_{\alpha\beta}(i\omega_n) = \int_{-\infty}^{\infty} d\omega \frac{A_{\alpha\beta}(\omega)}{i\omega_n-\omega},
$$
with
$$
A_{\alpha\beta}(\omega) = -\frac{1}{\pi}\,\mathrm{Im}\,G^R_{\alpha\beta}(\omega).
$$
The finite-temperature Lehmann representation is the basis for the QPE encoding, since the relevant spectral information is carried by energy differences between many-body eigenstates [2605.29681].

The central quantum innovation is a channel-agnostic finite-temperature QPE circuit. The data register is prepared in a Gibbs ensemble $\rho_{\mathrm{Gibbs}}=e^{-\beta H}/Z$ or an equivalent purification, and the circuit encodes difference phases $\Delta E = E_\lambda - E_{\lambda_0}$ rather than absolute many-body energies [2605.29681]. The controlled real-time evolution is
$$
U_{\mathrm{RTE}}(s) = \exp\!\left(-i \frac{2\pi H t_0 s}{N_{\mathrm{val}}}\right),
$$
where $N_{\mathrm{val}} = 2^{n_{q\mathrm{val}}}$ sets the energy-register size. Because the first controlled-evolution block contributes the phase of the initial state and the second contributes the phase of the excited state, the energy register records only their difference [2605.29681].

The channel-agnostic property means that the protocol does not require identifying per shot whether a measurement arose from a particle-addition or particle-removal process. For diagonal excitations, one ancilla is sufficient; for off-diagonal excitations, two ancillas and linear combinations of operators $a_m$, $a_{m'}$, and their adjoints generate four outcomes labeled by $(\xi,\sigma)\in\{e,h\}\times\{+,-\}$ [2605.29681]. The measured QPE probability for a bin $j$ is not a delta function but a filtered version of the true spectrum,
$$
P^{(p)}_j(E) = \left| \frac{1}{N_{\mathrm{val}}}\sum_{j'=0}^{N_{\mathrm{val}}-1}
\exp\!\left[i 2\pi j' \frac{(E-E_{\mathrm{orig}}^{(p)})t_0-j}{N_{\mathrm{val}}}\right]\right|^2,
$$
which is the QPE filter kernel for grid setting $p$ [2605.29681].

A plausible implication is that the method trades explicit channel resolution and absolute-energy estimation for a representation more directly aligned with the Lehmann structure of the Green’s function. That trade is central to its finite-temperature applicability.

## 3. QAVG reconstruction over variable grids

QAVG is the classical inference layer that converts QPE histograms into a reconstructed spectral function and hence into impurity Green’s functions [2605.29681]. Rather than relying on a single QPE grid, the method samples multiple grid settings
$$
\mathcal{G}=\{g_p\},
$$
with different $t_0^{(p)}$, $E_{\mathrm{orig}}^{(p)}$, and $N_{\mathrm{val}}^{(p)}$. The stated purpose is to suppress spectral leakage and aliasing tied to any single Dirichlet kernel and to sharpen inference through variable-grid averaging [2605.29681].

For diagonal spectral components, the modeled probabilities are
$$
\mathbb{P}^{(p,m)}_{\xi j}(\Lambda_\xi) = \int dE\, S^{(\xi)}_{\mathrm{rec},mm}(E;\Lambda_\xi)\, P^{(p)}_j(E),
$$
and analogous expressions hold for off-diagonal matrix elements [2605.29681]. The reconstructed spectral matrices are parameterized in the natural-orbital basis. The measured one-body density matrix $\gamma$ is diagonalized to obtain natural-orbital occupancies and a unitary rotation $U_{\mathrm{nat}}$, after which the electron and hole Green’s functions are modeled as sums over fictitious channels with centers $\epsilon_{\xi\ell}$, widths $\Delta E_{\xi\ell}$, and orthonormal transition vectors $v^{(\nu)}_{\xi\ell}$ [2605.29681]. The vectors are constrained by a hyperspherical Householder parametrization, enforcing orthonormality in the natural-orbital index.

The line shape $\rho_{\xi\ell}$ is chosen “quadratic” on a finite interval and normalized to unity [2605.29681]. This is explicitly motivated by finite-temperature peak crowding, where thousands of addition and removal channels can contribute in narrow energy windows. The reconstruction therefore groups many thermally active transitions into a smaller number of fictitious broadened channels rather than attempting line-by-line recovery [2605.29681].

The fit is defined by minimizing a grid-averaged discrepancy between modeled and measured histograms. The basic distance is
$$
D(h,h') = \sum_{j=0}^{N_{\mathrm{val}}-1} g_j \frac{|h_j-h'_j|}{2},
\qquad
g_j \propto \exp(-\tau_{\mathrm{dec}} j/t_0),
$$
and the total per-channel cost is
$$
F_\xi(\Lambda_\xi)=\sum_m w_\xi^{(m)} F_\xi^{(m)}(\Lambda_\xi)
+\sum_{m\neq m'} w_\xi^{(mm')} F_\xi^{(mm')}(\Lambda_\xi),
$$
with the QAVG solution
$$
\Lambda_{\mathrm{opt}}=\arg\min_\Lambda [F_e(\Lambda_e)+F_h(\Lambda_h)].
$$
These weights emphasize low-energy bins and can be chosen from excitation counts or from the measured density matrix $\gamma$ [2605.29681].

Once the spectral model is fixed, the Matsubara Green’s function is recovered from
$$
G(i\omega_n)=\int_{-\infty}^{\infty} d\omega\,\frac{A(\omega)}{i\omega_n-\omega}.
$$
This yields a DMFT-compatible impurity Green’s function without analytic continuation from imaginary time [2605.29681].

## 4. DMFT embedding and self-consistency structure

Within DMFT, the impurity model is coupled to the lattice through the usual self-consistency equations. For a lattice Hamiltonian $H_{\mathrm{KS}}(\mathbf{k})$ and local self-energy $\Sigma_{\mathrm{loc}}(z)$, the lattice and local Green’s functions are
$$
G_{\mathrm{latt}}(\mathbf{k},z)=\left[z+\mu-H_{\mathrm{KS}}(\mathbf{k})-\Sigma_{\mathrm{loc}}(z)\right]^{-1},
$$
and
$$
G_{\mathrm{loc}}(z)=\frac{1}{N_k}\sum_{\mathbf{k}} G_{\mathrm{latt}}(\mathbf{k},z).
$$
Projected onto the correlated subspace,
$$
\mathcal{G}_0^{-1}(z)=P_{\mathrm{corr}}\left[G_{\mathrm{loc}}^{-1}(z)+\Sigma_{\mathrm{loc}}(z)\right]P_{\mathrm{corr}}.
$$
For the impurity model,
$$
G_{\mathrm{imp}}^{-1}(i\omega_n)=i\omega_n+\mu-\Delta(i\omega_n)-\Sigma_{\mathrm{imp}}(i\omega_n),
$$
with $\Sigma_{\mathrm{loc}}=\Sigma_{\mathrm{imp}}$ [2605.29681].

The high-level QAVG-DMFT loop is therefore:

| Step | Operation | Output |
|---|---|---|
| 1 | Build the AIM from current bath parameters | $H_{\mathrm{imp}}$ |
| 2 | Prepare a Gibbs-like input and run channel-agnostic QPE on selected circuits and grid settings | Histograms |
| 3 | Apply QAVG reconstruction | $S^{(\xi)}_{\mathrm{rec}}$, $G_{\mathrm{imp}}$ |
| 4 | Use Dyson’s equation | $\Sigma_{\mathrm{imp}}$ |
| 5 | Update lattice Green’s function, Weiss field, and bath fit | New $\mathcal{G}_0$, $\Delta$ |

The explicit algorithm in the paper initializes from DFT, constructs maximally localized Wannier orbitals, chooses a correlated subspace, sets an initial $\Sigma_{\mathrm{loc}}$, and iterates the impurity reconstruction and bath fitting inside the DMFT loop [2605.29681]. A key point is that the quantum device is used to sample finite-temperature spectral information directly, while the classical side handles the inverse problem and DMFT bath update.

This suggests that QAVG-DMFT is best viewed not as a new lattice approximation, but as a quantum-assisted impurity-solver architecture embedded in standard DMFT self-consistency.

## 5. Validation on SrVO$_3$ and relation to other quantum-assisted DMFT solvers

The explicit materials validation is carried out for cubic SrVO$_3$ [2605.29681]. The Anderson impurity model uses three correlated $t_{2g}$ orbitals per spin and $n_{\mathrm{bath}}=3$ bath sites, giving
$$
n_{\mathrm{sorb}} = 2(n_{\mathrm{corr}}+n_{\mathrm{bath}})=12
$$
spin orbitals. The Kanamori parameters are $U=3.44\ \mathrm{eV}$, $U_0=2.49\ \mathrm{eV}$, and $J=0.46\ \mathrm{eV}$, with temperature $\beta^{-1}=0.025\ \mathrm{eV}$ [2605.29681]. The numerical validation uses exact probabilities convolved with the QPE filter rather than noisy circuit simulations, thereby isolating the reconstruction problem [2605.29681].

For the QPE settings, the energy register has $n_{q\mathrm{val}}=7$, so $N_{\mathrm{val}}=128$, the grid spacing is $1/t_0=1/6\ \mathrm{eV}$, and three origin shifts are used [2605.29681]. In the one-shot reconstruction at DMFT iteration 10, the electron sector uses 6 independent fictitious excitation energies and the hole sector uses 2 [2605.29681]. The reconstructed density of states reproduces the exact-diagonalization benchmark qualitatively, including a central quasiparticle peak and Hubbard bands, while Matsubara $G(i\omega_n)$ agrees well except for modest discrepancies at the lowest frequencies [2605.29681]. In iterative QAVG-DMFT, the reconstruction uses up to 8 independent electron energies and 4 hole energies, each threefold degenerate, with 123 and 51 angle parameters for the electron and hole transition systems, respectively [2605.29681].

A related but distinct quantum-assisted DMFT strategy uses a symmetry-adapted VQE followed by real-time evolution of a discretized Anderson impurity model [2602.03843]. In that approach, for a four-site AIM the relative error in the ground-state energy remains well below $0.01\%$ with $N_p\le 30$, and the single-particle Green’s function is accurately extracted in intermediate to strong interaction regimes, while the weak interaction regime shows noticeable deviations in low-energy spectral features [2602.03843]. The comparison is instructive: both methods target the impurity Green’s function, but the VQE approach relies on real-time evolution from a ground state, whereas QAVG-DMFT uses finite-temperature QPE and a multi-grid statistical reconstruction [2602.03843] [2605.29681].

## 6. Relation to GW+DMFT, QSGW+DMFT, and other nonstandard uses of “QAVG”

Outside the finite-temperature QPE literature, “QAVG-DMFT” is used only as an interpretive label, not as an established scheme name. In one H$_2$ study, the closest match to such a label is identified as QSGW+DMFT with dynamical double-counting, where a nonlocal quasiparticle Hamiltonian is constructed from an averaged nonlocal quasiparticle renormalization matrix $\bar Z$ and DMFT supplies the local dynamical self-energy [1611.07090]. That work explicitly calls this mapping “the closest match to what we would call ‘QAVG-DMFT’,” but it also states that the paper itself studies fully self-consistent GW+DMFT, one-shot $G_0W_0$+DMFT, and QSGW+DMFT variants rather than a named QAVG-DMFT formalism [1611.07090]. In its H$_2$ benchmarks, the dynamical-double-counting variant improves spectra near equilibrium but fails to reach the atomic limit, whereas QSGW+DMFT with static double-counting and GW+DMFT with causal double-counting recover the atomic limit [1611.07090].

A review and methodological paper on QSGW+DMFT for iron pnictides likewise states that the term “QAVG-DMFT” is not used there, and instead develops a QS GW+DMFT scheme in which nonlocal static exchange/correlation is retained in a quasiparticle Hamiltonian and DMFT supplies local dynamical correlations [1411.5180]. The same terminological clarification appears in broader reviews of GW, GW+DMFT, and D$\Gamma$A, where a hypothetical “QAVG-DMFT” is mapped to a quasiparticle-averaged static GW self-energy combined with a local DMFT self-energy and projector-based double-counting subtraction [1109.3972] [1703.08446].

A different ambiguity arises in finite-dimensional DMFT for two-particle observables such as double occupancy. The self-consistent Dual Boson construction resolves discrepancies between impurity- and lattice-based routes by enforcing
$$
X_{\mathrm{loc}}(i\nu_m)=\chi_{\mathrm{imp}}(i\nu_m),
$$
which yields $D=D_{\mathrm{imp}}=D_{\mathrm{lat}}$ by construction [1602.09129]. That paper explicitly rejects any identification of “QAVG-DMFT” with a principled averaging scheme; if the label were interpreted as averaging inconsistent estimators, the paper characterizes that as ad hoc and contrasts it with a frequency-dependent two-particle self-consistency condition [1602.09129].

These uses do not define a common formalism. They instead show that the acronym can be informally attached to at least three conceptually distinct ideas: variable-grid QPE reconstruction for DMFT [2605.29681], quasiparticle-averaged GW+DMFT reinterpretations [1611.07090] [1703.08446], and, negatively, ad hoc averaging procedures that are explicitly disfavored in two-particle DMFT contexts [1602.09129].

## 7. Limitations, error channels, and practical significance

The explicit QAVG-DMFT formulation identifies several error sources. Single-grid QPE suffers from Dirichlet-kernel sidelobes, motivating variable-grid averaging [2605.29681]. Finite-temperature spectra contain large numbers of thermally active transitions within narrow energy windows, which motivates grouping them into fictitious broadened channels rather than attempting exact line resolution [2605.29681]. Hardware imperfections, readout errors, and imperfect thermal-state preparation further perturb the measured histograms; the paper separates optimization, modeling, statistical, and hardware errors through a triangle-inequality bound on the distance between ideal and reconstructed distributions [2605.29681].

The method also assumes the availability of a Gibbs-like input state and moderate-depth controlled time evolutions [2605.29681]. A plausible implication is that the practical frontier is set as much by thermal-state preparation and controlled-evolution quality as by the classical inverse problem. The SrVO$_3$ demonstration therefore validates the reconstruction concept more directly than a full noisy-hardware implementation [2605.29681].

In comparison with real-time VQE impurity solvers, QAVG-DMFT avoids long coherent real-time evolution as the primary spectral estimator and directly targets finite-temperature excitation energies through QPE [2602.03843] [2605.29681]. In comparison with CT-QMC-based DMFT, it aims to bypass analytic continuation [2605.29681]. In comparison with ED, it still uses a finite-bath impurity Hamiltonian but distributes the computational burden between quantum sampling and compact spectral modeling [2605.29681].

Taken in its explicit form, QAVG-DMFT is therefore a finite-temperature quantum-assisted impurity-solver framework for DMFT in which modified QPE produces channel-agnostic spectral data and QAVG converts those data into Green’s functions suitable for self-consistent many-body calculations [2605.29681]. The broader literature shows that the label has been used informally in other contexts, but those usages describe separate quasiparticle or averaging ideas rather than the finite-temperature QPE scheme itself [1411.5180] [1611.07090] [1703.08446].

Source: https://www.emergentmind.com/topics/qavg-dmft-scheme