---
title: Parallelized Amplitude Estimation (PAE)
url: https://www.emergentmind.com/topics/parallelized-amplitude-estimation-pae
type: topic
---

# Parallelized Amplitude Estimation (PAE)

Searching arXiv for recent and foundational papers on parallelized amplitude estimation.
arxiv_search(query="parallelized amplitude estimation quantum", max_results=10)
arxiv_search(query="Parallelized Amplitude Estimation", max_results=10)
Parallelized amplitude estimation (PAE) is a family of quantum amplitude-estimation methods in which information about an unknown amplitude is gathered from multiple independent quantum circuits—typically at different Grover depths, phase-estimation windows, or low-depth signal-processing schedules—and then fused by classical post-processing. In contrast to canonical quantum amplitude estimation based on quantum phase estimation (QPE), which uses controlled powers of the Grover operator and an inverse quantum Fourier transform in a single long coherent circuit, PAE reorganizes the task into batchable or modular subproblems that are amenable to parallel execution across shots, circuit families, phase windows, or distributed branches [1904.10246][2204.01337][2405.14697][2508.06121]. The literature includes nonadaptive multi-depth maximum-likelihood methods, explicit parallel QPE/QAE constructions, signal-processing-based noniterative schemes, modular windowed variants, and low-depth batchable algorithms based on QSP, Bayesian inference, or statistical eigengap estimation [2006.16223][2207.08628][2508.05805][2412.04394][2603.05475].

## 1. Formal setting and signal model

The underlying estimation problem is the standard amplitude-estimation task. A unitary $\mathcal A$ prepares a state of the form
\[
\ket{\Psi}=\mathcal A\ket{0}_{n+1}
=\sqrt{p}\,\ket{\Psi_1}\ket{1}+\sqrt{1-p}\,\ket{\Psi_0}\ket{0},
\]
where the last qubit indicates whether a basis state is “good” or “bad,” and the target quantity is the success probability $p$ or, equivalently, the good-state amplitude $\sqrt p$. The conventional angular parameterization is
\[
\theta=\arcsin\sqrt p,\qquad p=\sin^2\theta,\qquad \theta\in[0,\pi/2].
\]
With Grover operator
\[
\mathcal Q=-\mathcal A\,\mathcal S_0\,\mathcal A^{-1}\,\mathcal S_\chi,
\]
one has the usual amplitude-amplification rotation
\[
\mathcal Q^m\ket{\Psi}
=\sin((2m+1)\theta)\ket{\Psi_1}\ket{1}
+\cos((2m+1)\theta)\ket{\Psi_0}\ket{0},
\]
so that the probability of measuring “good” after depth $m$ is $\sin^2((2m+1)\theta)$ [2206.08449].

This sinusoidal dependence on Grover depth is the common signal model exploited across PAE variants. Some works parameterize the unknown directly as a success probability, while others parameterize the amplitude itself. In particular, the signal-processing-based csAE formulation writes
\[
U\lvert 0^l\rangle=\cos\theta\,\lvert x,0\rangle+\sin\theta\,\lvert x',1\rangle
\]
and sets $a=\sin\theta$, so that the depth-$0$ success probability is $a^2=\sin^2\theta$ rather than $a$ [2405.14697]. Despite such notational differences, the core estimation object remains the same: an angle or phase governing oscillatory statistics across amplified experiments.

From the standpoint of architecture, the central departure from canonical QAE is that the oscillatory information is extracted from many independent executions rather than from one monolithic phase-estimation circuit. In PAE-like schemes, the computational burden shifts from long coherent control to classical inference over measurements gathered at multiple depths, windows, or branches.

## 2. Emergence from non-QPE multi-depth estimation

A major precursor to PAE is the multi-depth maximum-likelihood approach of Suzuki et al., which estimates an unknown amplitude by running a family of circuits $\mathbf Q^{m_k}\ket{\Psi}$ at predetermined depths $\{m_0,m_1,\dots,m_M\}$, collecting counts $h_k$, and maximizing the joint likelihood over $\theta_a\in[0,\pi/2]$ [1904.10246]. The depth schedules studied explicitly are the linearly incremental sequence
\[
m_k=k
\]
and the exponentially incremental sequence
\[
m_0=0,\quad m_1=2^0,\quad m_2=2^1,\quad \dots,\quad m_M=2^{M-1}.
\]
Because these schedules are fixed in advance, every depth circuit can be executed independently. The paper explicitly contrasts this with adaptive simplified-counting approaches and states that the method “can be run in parallel on multiple quantum devices” [1904.10246].

The statistical mechanism is classical maximum likelihood. If depth $m_k$ is repeated $N_k$ times and produces $h_k$ good outcomes, then the likelihood contribution is
\[
L_k(h_k;\theta_a)
=
\left[\sin^2((2m_k+1)\theta_a)\right]^{h_k}
\left[\cos^2((2m_k+1)\theta_a)\right]^{N_k-h_k},
\]
and the estimator maximizes the product over depths. The asymptotic analysis distinguishes two regimes: LIS yields $\hat\epsilon\sim N_{\mathrm q}^{-3/4}$, whereas EIS yields $\hat\epsilon\sim N_{\mathrm q}^{-1}$, matching Heisenberg-like scaling in total query count [1904.10246]. This fixed-schedule, multi-depth, classically fused architecture is one of the clearest antecedents of later PAE terminology.

The noise-aware extension formulates the same multi-depth architecture as a two-parameter estimation problem for amplitude and noise. Under a depolarizing-noise model with parameter $\kappa=-\ln p$, the success probability becomes
\[
P(m;a,\kappa)=\frac{1}{2}-\frac{1}{2}e^{-\kappa m}\cos\!\big(2(2m+1)\theta_a\big),
\]
and the estimator jointly maximizes the corresponding likelihood over $(a,\kappa)$ [2006.16223]. This work shows that noisy PAE can retain quantum query advantage for a range of budgets, but also that the estimation error saturates and that the Fisher information matrix can become nearly singular for anomalous target amplitudes such as $a=\sin^2(\pi/8)\approx 0.146$ [2006.16223].

## 3. Principal parallelized architectures

The literature now contains several distinct architectural interpretations of PAE. Some parallelize the phase-kickback mechanism itself, some parallelize predetermined Grover-depth experiments, and some modularize QPE into independently executable windows.

| Approach | Parallel element | Characteristic claim |
|---|---|---|
| Parallel QPE/QAE [2204.01337] | Parallel kickbacks, duplicated work registers, entangled controls, or reinitialization | Parallelization can reduce depth and improve resilience for specific error models |
| csAE [2405.14697] | Predetermined Grover-depth set $D$ with independent execution and ESPRIT post-processing | Reported $\sim 4.3/\varepsilon$ sequential and $\sim 0.26/\varepsilon$ parallel query complexity at $95\%$ confidence |
| AWQAE [2508.05805] | Independent phase windows/chunks with blockwise QPE-like circuits | Precision $\mathcal O(1/2^n)$ with per-block coherent qubit cost $m_{\rm pr}+m_i+n_T+1$ |
| GHZ/QSP PAE [2508.06121] | Global GHZ state followed by separated low-depth branches | $N=\mathcal O(1/\varepsilon+P\log P)$ and depth $\mathcal O(1/(P\varepsilon)+\log P)$ |

The explicit parallel QPE/QAE construction of Imhof et al. reorganizes amplitude estimation by parallelizing the phase-kickback generation that underlies QPE. The work describes three styles: simple parallelization with many copies of the work register, indirect parallelization using entangled control qubits so that many single kickbacks reproduce higher powers, and reinitialization-based recycling when width is limited [2204.01337]. For amplitude estimation, the key caveat is that each branch must be initialized in a single Grover eigenstate rather than in the usual superposition used by canonical QAE.

The csAE method recasts amplitude estimation as a classical direction-of-arrival problem. The amplified state
\[
\phi_n=G^nU\lvert 0^l\rangle
=\cos((2n+1)\theta)\lvert x,0\rangle+\sin((2n+1)\theta)\lvert x',1\rangle
\]
induces the signal
\[
y_n=e^{i(4\theta n+2\theta)},
\]
so $\theta$ appears as a frequency that can be extracted by ESPRIT from measurement data obtained at sparse, predetermined depths [2405.14697]. The method is explicitly nonadaptive, phase-estimation-free, and designed so that all depth experiments can be run independently in parallel.

AWQAE preserves the QPE/QAE viewpoint but partitions the phase register into windows of sizes $[m_1,m_2,\dots,m_B]$ and reconstructs the full phase classically from small QPE-like blocks [2508.05805]. Each block uses a 2-qubit phase-resolution register, an $m_i$-qubit counting register, the target register, and one ancilla. The architecture is described as modular, scalable, and adaptive, with blocks “independent” and “amenable to parallel processing,” although the final estimate still depends on an ambiguity-aware LSB-to-MSB correction routine [2508.05805].

The most explicit recent use of the term PAE is the GHZ/QSP framework of Yang et al. It prepares a $P$-qubit GHZ state,
\[
\ket{\mathrm{GHZ}_P}=\frac{\ket{0}^{\otimes P}+\ket{1}^{\otimes P}}{\sqrt2},
\]
then applies separated QSP-constructed phase-shift circuits $V_{\varphi,T}$ on each branch, so that the effective signal amplification is proportional to $M=PT$ [2508.06121]. This yields a tunable depth-width-query tradeoff and, in the fully parallel regime $P=\Theta(1/\varepsilon)$, a logarithmic-depth algorithm with near-Heisenberg-limited total query scaling.

## 4. Statistical post-processing paradigms

Although the quantum front ends differ, PAE methods are equally distinguished by their classical inference layers. Maximum-likelihood estimation is the foundational paradigm: all observed counts from different depths are combined into a single one-parameter or two-parameter likelihood and optimized globally [1904.10246][2006.16223]. In this formulation, the nonuniqueness of $\sin^2((2m+1)\theta)$ at a single depth is resolved by cross-consistency across several depths.

QSP-based amplitude-estimation frameworks replace Grover-only oscillations by designed polynomial response functions. The sampling theorem of the QSP approach states that one can sample a bit $\mathbf b_P\in\{0,1\}$ with
\[
\Pr[\mathbf b_P=1]=|P(a)|^2
\]
for QSP-implementable polynomials $P$, using $O(\deg P)$ oracle queries [2207.08628]. Its hybrid theorem makes the depth-for-repetitions tradeoff explicit:
\[
\mathbf D\in \tilde O\!\left((a^{-1}+\varepsilon^{-(1+\beta)})\log(\delta^{-1})\right),\qquad
\mathbf d\in \tilde O\!\left(a^{-1/(1-\beta)}+\varepsilon^{-(1-\beta)}\right),
\]
where $\mathbf D$ is total degree and $\mathbf d$ is maximum degree used in any one circuit [2207.08628]. The framework is adaptive across batches, but all repetitions within a batch are independent and parallelizable.

Bayesian amplitude estimation supplies a different inference backbone. BAE represents the posterior over $\theta$ by sequential Monte Carlo particles, chooses the next Grover depth by expected utility, and updates after batches of repeated shots [2412.04394]. The paper emphasizes that reweighting and resampling are parallelizable over particles, that utility evaluations are parallelizable over candidate controls, and that $Nshots>1$ shots at a chosen depth can be run in parallel on multiple devices [2412.04394]. The annealed variant replaces variance-minimizing utility with an ESS-based design criterion, but the outer loop remains adaptive.

A further statistical reformulation treats amplitude estimation as eigengap estimation of an effective Hamiltonian generated by Grover dynamics. In this picture,
\[
a=\sin^2(\lambda),\qquad \Delta_{\rm eff}=4\lambda,
\]
and the estimation target is extracted from expectation values such as $\mathbb E[Z_m\mid m]=\cos(2\lambda m)$ under nonadaptive Gaussian schedules [2603.05475]. GLSAE minimizes a least-squares objective
\[
\tilde{\mathcal L}(\theta)=\frac{1}{N}\sum_{m\sim \tilde p_T}\big(Z_m-\cos(2\theta m)\big)^2,
\]
while GDMAE, in the flag-qubit setting, maximizes
\[
\tilde{\mathcal M}(\theta)=\frac{1}{N}\sum_{m\sim \tilde q_T}\Big(Z_m\cos(2\theta m)+X_m\sin(2\theta m)\Big).
\]
Both algorithms are noniterative and highly parallelizable across sampled schedule points [2603.05475].

## 5. Resource tradeoffs, performance, and noise

PAE does not remove the fundamental query cost of amplitude estimation; rather, it redistributes that cost across depth, width, and classical computation. This is explicit in the GHZ/QSP PAE theorem,
\[
N=\mathcal O\!\left(\frac{1}{\varepsilon}+P\log P\right),\qquad
\text{depth}=\mathcal O\!\left(\frac{1}{P\varepsilon}+\log P\right),\qquad
\text{width}=P(n+1),
\]
which reduces to
\[
N=\mathcal O\!\left(\frac{\log(1/\varepsilon)}{\varepsilon}\right),\qquad
\text{depth}=\mathcal O(\log(1/\varepsilon))
\]
for $P=\Theta(1/\varepsilon)$ [2508.06121]. Here the

Source: https://www.emergentmind.com/topics/parallelized-amplitude-estimation-pae