---
title: Gaussian-Filtered Quantum Phase Estimation
url: https://www.emergentmind.com/topics/gaussian-filtered-quantum-phase-estimation
type: topic
---

# Gaussian-Filtered Quantum Phase Estimation

Gaussian-filtered quantum phase estimation denotes a family of phase-estimation constructions in which a Gaussian profile is used to localize spectral or phase information before, during, or instead of textbook quantum phase estimation. In the formulation most directly identified with the term, a Gaussian spectral filter
$$
f(H)=\exp\!\left[-\frac{(H-\mu)^2}{2\sigma^2}\right]
$$
is applied to an input state $|\phi_0\rangle=\sum_i \gamma_i |E_i\rangle$ so that the desired eigencomponent is amplified prior to a final high-precision QPE stage [2510.04294]. Related uses of Gaussian filtering include truncated Gaussian time-sampling for multiple-eigenvalue search [2402.01013], Gaussian spin states whose measurement update acts as a local Gaussian phase filter [2010.04001], and linearized continuous-time quantum filtering for small interferometric phase shifts [1601.04374].

## 1. Spectral-filter formulation

In filtered QPE, the starting point is a normalized Hamiltonian $H$ with spectrum in $[-1,1]$ and a reference state
$$
|\phi_0\rangle=\sum_i \gamma_i |E_i\rangle.
$$
The Gaussian filter is centered at $\mu$, ideally close to the target energy $E_0$, and has bandwidth set by $\sigma$:
$$
f(H)=\exp\!\left[-\frac{(H-\mu)^2}{2\sigma^2}\right].
$$
Acting on the reference state gives the unnormalized filtered state
$$
|\tilde\phi_f\rangle=f(H)|\phi_0\rangle=\sum_i \gamma_i f(E_i)|E_i\rangle,
$$
and the normalized post-selected state is
$$
|\phi_f\rangle=\frac{|\tilde\phi_f\rangle}{\|\tilde\phi_f\|}.
$$
The corresponding block-encoding requirement is a unitary $U_f$ on $n+m$ ancilla qubits such that
$$
(\langle 0^m|\otimes I)\,U_f\,(|0^m\rangle\otimes I)\approx f(H),
$$
that is, a $(1,m,\epsilon_f)$ block-encoding of $f(H)$ [2510.04294].

The central quantity is the overlap amplification achieved after filtering. The success probability of post-selection is
$$
p_f=\langle \phi_0|f(H)^\dagger f(H)|\phi_0\rangle=\sum_i |\gamma_i f(E_i)|^2,
$$
while the amplified ground-state overlap is
$$
\gamma_{f0}=\langle E_0|\phi_f\rangle=\frac{\gamma_0 f(E_0)}{\sqrt{p_f}},
$$
so that
$$
|\gamma_{f0}|^2=\frac{1}{1+R_f},
\qquad
R_f=\sum_{i>0}\left|\frac{\gamma_i f(E_i)}{\gamma_0 f(E_0)}\right|^2.
$$
If the Gaussian is peaked at $E_0$ with $f(E_0)=1$ and satisfies $|f(E)|\le \epsilon_g$ outside the main lobe $|E-E_0|\ge \Delta$, then
$$
R_f\le (|\gamma_0|^{-2}-1)\epsilon_g^2,
$$
hence
$$
|\gamma_{f0}|^2\ge \frac{1}{1+(|\gamma_0|^{-2}-1)\epsilon_g^2}.
$$
This makes explicit how excited-state suppression improves the effective input to the subsequent QPE stage [2510.04294].

The same formalism also makes explicit the trade-off between overlap amplification and preparation cost. The filtered success probability obeys
$$
p_f\le \frac{|\gamma_0|^2}{|\gamma_{f0}|^2}.
$$
Accordingly, sharper filtering can improve the final overlap but may reduce post-selection probability, so the advantage of Gaussian-filtered QPE is governed by the balance between overlap amplification and the query cost of implementing the filter.

## 2. Implementations and the two-stage protocol

The Gaussian filter can be realized in several ways. In the trigonometric realization, one uses
$$
\exp\!\left(-\frac{x^2}{2\sigma^2}\right)\approx \sum_{k=-N/2}^{N/2} c_k e^{i\pi k x},
$$
with $N=O((1/\Delta)\log \epsilon_f^{-1})$, where $\Delta$ is the desired spectral resolution, approximately $\sigma$. This is implemented via generalized quantum signal processing, with depth $\sim N$ calls to $e^{\pm i\pi H}$. In the polynomial realization, the Gaussian is expanded in Chebyshev polynomials,
$$
\exp\!\left(-\frac{x^2}{2\sigma^2}\right)\approx \sum_{k=0}^N d_k T_k(x),
$$
with $N=O(\Delta^{-1}\log \epsilon_f^{-1})$, and is implemented by QSVT or qubitization using $O(N)$ calls to the qubitization oracle. A third route uses Krylov-subspace filters: choosing a basis $\{b_k(H)\}$, forming
$$
S_{kl}=\langle \phi_0|b_k(H)^\dagger b_l(H)|\phi_0\rangle,
\qquad
H_{kl}=\langle \phi_0|b_k(H)^\dagger H b_l(H)|\phi_0\rangle,
$$
solving
$$
Hc=Sc\,E^{(N)},
$$
and taking the minimizer $c=c^{(0)}$ to define
$$
f(H)=\sum_k c_k b_k(H).
$$
This yields an adaptive filter concentrated on the low-lying spectrum [2510.04294].

The filtered-QPE protocol is organized in three stages. Stage 1 is optional coarse estimation: one runs standard low-precision QPE of depth $O((\epsilon' \Delta E_0)^{-1})$ and shots $O(|\gamma_0|^{-2}\log \delta_1^{-1})$ to obtain $\tilde E_0$ and $\tilde E_1$ within error $\epsilon' \Delta E_0$. Stage 2 prepares the filtered state by fixing
$$
f(H)=\exp\!\left[-\frac{(H-\tilde E_0)^2}{2\sigma^2}\right],
$$
with
$$
\sigma \simeq \frac{(1-\epsilon')(\tilde E_1-\tilde E_0)}{\sqrt{2\log \epsilon_g^{-1}}},
$$
implementing a block-encoding of order
$$
N_{\rm sp}=O[(\Delta E_0)^{-1}\log \epsilon_g^{-1}],
$$
and repeating the post-selection procedure until success, requiring approximately $M_{\rm sp}\simeq p_f^{-1}$ attempts. Stage 3 then runs standard QPE on $|\phi_f\rangle$ with depth $D\simeq \epsilon^{-1}$ and number of shots
$$
M_{\rm QPE}\simeq |\gamma_{f0}|^{-2}\log \delta_2^{-1}.
$$
The overall success probability satisfies
$$
(1-\delta_1)(1-\delta_2)\simeq 1-\delta
$$
[2510.04294].

A notable aspect of this construction is that the Gaussian filter is not itself the final estimator. It is a state-preparation primitive inserted before the standard phase-readout stage, and its value lies in improving the overlap properties of the state on which high-precision QPE is executed.

## 3. Complexity and empirical performance

For filtered QPE, the total expected cost in queries to $e^{\pm i\pi H}$ or to qubitization oracles is
$$
C_{\rm tot}\simeq M_{\rm sp}N_{\rm sp}+M_{\rm QPE}D
      \simeq |\gamma_0|^{-2}N_{\rm sp}+|\gamma_{f0}|^{-2}\epsilon^{-1},
$$
using the identity
$$
p_f^{-1}|\gamma_{f0}|^{-2}=|\gamma_0|^{-2}.
$$
By contrast, standard QPE has cost
$$
C_{\rm std}=O(|\gamma_0|^{-2}\epsilon^{-1}\log \delta^{-1}).
$$
The point of the Gaussian filter is therefore to replace the direct $|\gamma_0|^{-2}\epsilon^{-1}$ dependence by a two-term cost in which the overlap penalty is paid primarily at the filter-preparation stage [2510.04294].

A parameter choice is obtained by balancing filter preparation and final QPE depth through $N_{\rm sp}\simeq D=\epsilon^{-1}$. This gives
$$
\epsilon_g=
\sqrt{\frac{5\epsilon}{4\pi \Delta E_0}\,
W\!\left(\frac{4\pi \Delta E_0}{5\epsilon}\right)},
$$
where $W(\cdot)$ is the Lambert $W$-function. Under the conditions
$$
8.3\times 10^{-11}<\epsilon/\Delta E_0\le 5.4\times 10^{-2},
\qquad
\epsilon'\le 0.2,
$$
Theorem 1 shows
$$
C_{\rm tot}\le
C_{\rm std}\times
2e
\left[
\frac{5\epsilon}{4\pi \Delta E_0}
W\!\left(\frac{4\pi \Delta E_0}{5\epsilon}\right)
\right]^{1-\epsilon'}
=
O(|\gamma_0|^{-2}\epsilon^{-1})\cdot
(\epsilon/\Delta E_0\log(\Delta E_0/\epsilon))^{1-\epsilon'}.
$$
The corresponding corollary gives
$$
C_{\rm tot}=\tilde O\!\left(\epsilon^{-1}+|\gamma_0|^{-2}\Delta E_0^{-1}\right),
$$
so that in the high-precision regime $\epsilon\ll \Delta E_0$, the $\epsilon^{-1}$ term dominates and the dependence on $|\gamma_0|^{-2}$ is reduced [2510.04294].

Numerical experiments were carried out on 1D $(6,7)$-site and 2D $(2\times 3)$-site Fermi-Hubbard Hamiltonians with Bravyi-Kitaev tapering and spectrum normalized to $[-1,1]$, using reference Néel product states. The reported parameter range includes initial overlaps
$$
|\gamma_0|^2\simeq 1.5\times 10^{-2}\ldots 2.7\times 10^{-3}
$$
and spectral gaps
$$
\Delta E_0\simeq 3\times 10^{-3}\ldots 3\times 10^{-2}.
$$
For Gaussian FQPE, both trigonometric and polynomial filters yielded overlap amplification factors
$$
|\gamma_{f0}|^2/|\gamma_0|^2\gtrsim 100\times,
$$
best cost-reduction factors of approximately $1/300\ldots 1/350$ for $\epsilon/\Delta E_0=10^{-5}$, and still $\gtrsim 1/30$ for $\epsilon/\Delta E_0=10^{-3}$. The worst case over prior errors up to $\epsilon'=0.2$ remained below unity in the high-precision regime. Krylov filters with dimension $N\sim 60$ gave similar overlap amplification but very small success probabilities unless modified; with a small penalty $\Lambda$, the modified KSD procedure recovered
$$
p_f\sim 10^{-3}\ldots 10^{-2},
\qquad
|\gamma_{f0}|^2\sim 0.1,
$$
and a cost reduction of about $1/312$ [2510.04294].

## 4. Related Gaussian-filtered phase-estimation paradigms

The broader literature uses Gaussian filtering in several technically distinct ways [2510.04294] [2402.01013] [2010.04001] [1601.04374].

| Method | Gaussian object | Salient feature |
|---|---|---|
| FQPE | $f(H)=\exp[-(H-\mu)^2/(2\sigma^2)]$ | Filtered-state preparation followed by standard QPE |
| QMEGS | Truncated Gaussian time sampling and $F_T(x)\approx e^{-T^2x^2/2}$ | Hadamard-test filter-and-search for multiple eigenvalues |
| GSS cascade | Local Gaussian phase filter of width $\sigma_k\simeq s_k/\sqrt{N_k}$ | Cascade QPE with Gaussian spin states |
| Interferometric filtering | Linearized Gaussian quantum filter | Continuous nondemolition phase tracking |

In Quantum Multiple Eigenvalue Gaussian filtered Search, the Gaussian enters through the truncated time-sampling distribution
$$
a_T(t)=
\left(1-\int_{|s|\le \sigma T}(2\pi T^2)^{-1/2}e^{-s^2/(2T^2)}ds\right)\delta_0(t)
+(2\pi T^2)^{-1/2}e^{-t^2/(2T^2)}1_{|t|\le \sigma T}(t),
$$
whose Fourier transform is the filter kernel
$$
F_T(x)=\int_{-\sigma T}^{\sigma T}a_T(t)e^{ixt}dt\approx e^{-T^2x^2/2}.
$$
The algorithm assumes access to controlled-$U=\text{controlled-}e^{-iHt}$ and one ancilla qubit, uses Hadamard-test subroutines to obtain unbiased estimates of $\langle \psi|e^{-iHt_n}|\psi\rangle$, and then defines
$$
G(\theta)=\left|\frac1N\sum_{n=1}^N Z_n e^{i\theta t_n}\right|.
$$
On average,
$$
G(\theta)\approx \sum_{m=1}^M p_m F_T(\theta-\lambda_m),
$$
so peaks of $G(\theta)$ locate eigenvalues. The paper states that QMEGS is the first algorithm to simultaneously satisfy two properties: Heisenberg-limited scaling without any spectral gap assumption, and, under a positive energy gap and additional assumptions on the initial state, estimation of all dominant eigenvalues to $\epsilon$ accuracy with significantly reduced circuit depth compared to standard QPE. In the most favorable scenario, the maximal runtime can be as low as $\log(1/\epsilon)$ [2402.01013].

In the Gaussian-spin-state cascade, the Gaussian is built into the ancilla states themselves. A Gaussian spin state is
$$
|\psi(N,s)\rangle=
\frac1{\sqrt{\mathcal N}}
\sum_{\mu=-N/2}^{N/2}
\exp[-\mu^2/(s^2N)]\,|\mu\rangle_y,
$$
with squeezing parameter $s<1$. At step $k$, the outcome distribution can be interpreted through the filter function
$$
W_k(\mu;s_k)=|c_k(\mu)|^2=\exp[-\mu^2/(s_k^2N_k)]/\mathcal N_k,
$$
which in the phase variable corresponds to a local Gaussian filter of width
$$
\sigma_k\simeq s_k/\sqrt{N_k}.
$$
The cascade yields
$$
\Delta \theta_{\rm est},K \simeq \alpha_K N_T^{-(1-1/(2\cdot 3^K))},
\qquad
\alpha_{K\to\infty}=4,
$$
so that for $K\simeq 3$, $\Delta \theta_{\rm est}\approx O(N_T^{-0.98})$, and in the full-cascade limit $\Delta \theta_{\rm est}\to 4/N_T$, with total time $T\approx N_T t_U$ [2010.04001].

The interferometric work of Gough uses “Gaussian” in a different sense: a continuous-time quantum filter becomes exactly linear and Gaussian after linearization around a small phase offset. For homodyne detection and $\phi_0=\pi/2$,
$$
S_{11}(\phi)\approx \frac{1-i}{2}-\frac12\,\delta\phi,
$$
and the conditional state estimate and variance,
$$
x_t=\mathrm{Tr}[\rho_t\,\delta\phi],
\qquad
V_t=\mathrm{Tr}[\rho_t(\delta\phi)^2]-x_t^2,
$$
satisfy
$$
dx_t=-K_t\,dI(t),
\qquad
dV_t=-K_t^2\,dt,
$$
with
$$
K_t=\tfrac12\,\mathrm{Re}\,\beta
$$
for constant coherent amplitude. This is exactly the Kalman-Bucy form, and the mean-square estimation error tends to zero under continuous homodyne probing [1601.04374].

## 5. Conceptual distinctions and common misconceptions

The coexistence of these constructions suggests that “Gaussian-filtered quantum phase estimation” is not a single algorithmic template but a family of methods in which the Gaussian appears at different points of the estimation pipeline [2510.04294] [2402.01013] [2010.04001] [1601.04374]. In FQPE, the Gaussian acts directly on the Hamiltonian spectrum through $f(H)$. In QMEGS, the Gaussian is a time-domain sampling distribution whose Fourier transform yields a spectral search kernel. In the GSS cascade, the Gaussian is encoded in the spin-state amplitudes and induces a local Gaussian phase filter after measurement. In the interferometric setting, the Gaussian structure arises from a linearized conditional dynamics rather than from a discrete spectral filter.

A recurrent misconception is to equate overlap amplification with a free improvement in QPE complexity. The filtered-state formalism makes the opposite point explicit: overlap amplification and success probability are linked by
$$
p_f\le |\gamma_0|^2/|\gamma_{f0}|^2.
$$
The advantage of FQPE comes from a favorable redistribution of cost between state preparation and final phase readout, not from eliminating the preparation burden altogether [2510.04294].

Another misconception is that Gaussian filtering necessarily requires a large QPE register or a quantum Fourier transform. QMEGS explicitly avoids both: it uses one ancilla, Hadamard-test circuits, and a classical grid search, with no Fourier transform or root-finding beyond the search over $\theta_j=-\pi+j\cdot(q/T)$ [2402.01013]. Conversely, filtered QPE in the sense of FQPE does use standard QPE as its final high-precision stage; the Gaussian filter is a preconditioning step, not a replacement for phase readout.

A further distinction concerns target structure. FQPE is organized around improving the overlap with a single target eigenstate, especially the ground state. QMEGS is formulated as a multiple-eigenvalue estimator for the set of dominant eigenvalues. The GSS cascade estimates an unknown phase $\theta\in[-\pi,\pi)$ encoded by a controlled-$U$ gate through successive residual-phase updates. The interferometric Kalman-Bucy filter tracks a small continuous phase fluctuation $\delta\phi_t$ under continuous nondemolition measurement. These are related estimation problems, but they are not interchangeable.

## 6. Limitations and open research directions

Gaussian-filtered QPE remains parameter sensitive. In the FQPE formulation, the filter center must satisfy $\mu\simeq$ prior $E_0$, and the width
$$
\sigma \simeq \frac{(1-\epsilon')(E_1-E_0)}{\sqrt{2\log \epsilon_g^{-1}}}
$$
must be chosen so that the Gaussian covers $E_0$ while rejecting $E_1$. The implementation error, including Trotter or imperfect qubitization error, must satisfy
$$
\le O(\Delta E_0\,\epsilon_f)
$$
to preserve the filter profile accurately. The reported advantage is therefore tied to the regime $\epsilon\ll \Delta E_0$ and prior accuracy $\epsilon'\lesssim 0.2$ [2510.04294].

QMEGS has a different set of open problems. The classical search over $O(T/\epsilon)$ grid points has cost $O(\epsilon^{-1})$, and multiscale or binary-search strategies are proposed as a possible route to polylogarithmic cost in $1/\epsilon$. Adaptive tuning of $\sigma$, $T$, $N$, and $q$ without prior spectral-gap or overlap information remains open. The analysis assumes ideal controlled-$U$ access, and although the incorporation of realistic Hamiltonian-simulation errors or noise models is described as straightforward in many cases, detailed bounds are left for further study. The paper also identifies Gaussian-derivative, Heaviside-derivative, and Kaiser windows as possible alternative kernels [2402.01013].

The GSS cascade is robust against several noise sources, but its asymptotics change under strong depolarization. Full depolarization in the symmetric subspace imposes a minimum squeezing
$$
s_{\min}^2=\mathcal E N/3.
$$
As long as $s_k\ge s_{\min}$, the ideal Heisenberg scaling can be recovered; beyond that point, one can no longer increase $k$ and instead repeats the best-squeezed step, yielding
$$
\Delta \theta_{\rm est}\simeq \beta_{\mathcal E} N_T^{-1/2},
\qquad
\beta_{\mathcal E}\sim \mathcal E^{1/4}
$$
[2010.04001].

The continuous-time interferometric filter is restricted to the linear Gaussian regime. Its exact Kalman-Bucy form requires $|\delta\phi_t|\ll 1$ and neglects higher-order phase curvature. A nonzero coherent intensity $\beta$ is required; with vacuum input, the innovations vanish. Detector inefficiency, optical loss, and finite detector bandwidth are omitted, although the paper notes that such effects can be incorporated by standard quantum filtering extensions. Possible extensions include squeezed-state probes and adaptive feedback of the form
$$
\beta\mapsto \beta+\text{(gain)}\times I(t)
$$
[1601.04374].

Taken together, these limitations clarify the present status of Gaussian filtering in phase estimation. The Gaussian profile is a versatile device for suppressing unwanted spectral or phase contributions, but its practical benefit depends on how precisely the filter can be centered, how cheaply it can be implemented, and whether the underlying estimation problem is discrete, multiple-eigenvalue, collective-spin, or continuous-time interferometric.

Source: https://www.emergentmind.com/topics/gaussian-filtered-quantum-phase-estimation