---
title: Quantum Likelihood Estimation
url: https://www.emergentmind.com/topics/quantum-likelihood-estimation-qle
type: topic
---

# Quantum Likelihood Estimation

Quantum Likelihood Estimation (QLE) is a hybrid quantum-classical method for Hamiltonian learning in which measurement data generated by a quantum system evolving under an unknown Hamiltonian are converted into likelihoods over a candidate set of Hamiltonians, and those likelihoods are incorporated through Bayesian inference into an iteratively refined classical weight distribution. In the formulation developed most explicitly in "Optimal Quantum Likelihood Estimation" [2509.00825], QLE is a prepare–evolve–measure–update loop whose efficiency depends strongly on experiment design at each iteration; the paper makes that design information-theoretically optimal, iteration by iteration, by choosing controls that maximize the mutual information between the hidden Hamiltonian label and the measurement outcome. Earlier work had already realized QLE experimentally as a component of quantum Hamiltonian learning, using a photonic quantum simulator to estimate likelihoods for an NV-center Hamiltonian and to drive sequential Bayesian updates [1703.05402].

## 1. Definition and historical placement

In its quantum-information usage, QLE addresses the problem of identifying which Hamiltonian in a candidate set \(\{H_1,\dots,H_N\}\) generated observed dynamics under an unknown true Hamiltonian \(H_{\text{true}}\). The method maintains a classical probability distribution over the candidate Hamiltonians and refines that distribution as additional quantum measurement data are acquired [2509.00825].

This places QLE squarely in the NISQ-era class of hybrid quantum-classical algorithms: a quantum device performs the expensive data-generating step, while a classical routine updates beliefs and chooses subsequent experiments. The original experimental realization in quantum Hamiltonian learning interfaced two quantum systems through a classical computer: a negatively charged nitrogen-vacancy center in diamond as the target system and a programmable silicon photonic quantum device as the likelihood engine. In that experiment, the likelihood for a candidate Hamiltonian parameter \(\vec x\) was the return probability
\[
\Pr(D|\vec x)=\left|\langle D|e^{-i\hat H(\vec x)t}|\psi\rangle\right|^2,
\]
with \(|D\rangle=|\psi\rangle=|m_s=0\rangle\), and Bayesian inference was implemented approximately by Sequential Monte Carlo (SMC) [1703.05402].

That experimental program established several points that remain central to QLE. Using **50 steps** and **20 SMC particles**, the protocol learned an NV-center Rabi frequency with final estimate
\[
f_{\text{QLE}}=(6.93\pm0.09)\,\text{MHz},
\]
consistent with the independently fitted value \(f_0=6.90\,\text{MHz}\), and the paper reports a final quadratic loss of approximately \(10^{-5}\). The same work also showed that saturation of posterior variance can diagnose model mismatch: in its simplest model the variance saturated at about \(\sigma^2(\omega)\simeq 4.2\times 10^{-5}\), motivating an improved chirped Hamiltonian model with Bayes factor \(K=560\). An interactive variant, IQLE, achieved
\[
f_{\text{IQLE}}=(6.92\pm0.08)\,\text{MHz}
\]
with quadratic loss about \(10^{-7}\) [1703.05402]. These results established QLE not only as a parameter-estimation procedure but also as a model-validation method.

## 2. Core algorithmic structure

At its core, QLE is a prepare–evolve–measure–update protocol. In each iteration, the algorithm begins with an initial state \(\ket{\psi_1}\), evolves it under the unknown unitary
\[
U(t)=e^{-iH_{\text{true}}t},
\]
applies a final unitary \(W\) selecting the measurement basis, and measures in the computational basis. The outcome \(Y\in\{0,1\}^n\) is then used to compute likelihoods for each candidate Hamiltonian \(H_j\), after which Bayesian inference updates the weights \(w_j^{(k)}\) [2509.00825].

The Bayesian update rule is
\[
w_j^{(k+1)}=\frac{w_j^{(k)}\cdot \mathcal{L}_j}{\sum_{i=1}^{N} w_i^{(k)}\cdot \mathcal{L}_i},
\]
where \(\mathcal{L}_j=p_j^{(Y)}\) is the likelihood of observing outcome \(Y\) under hypothesis \(H_j\). The posterior from one round becomes the prior for the next. In this sense, QLE is not a one-shot estimator but an adaptive sequential inference scheme.

For the optimized single-qubit setting analyzed in [2509.00825], the controls are the initial-state angles \((\alpha,\beta)\),
\[
\ket{\psi_1}=\cos(\alpha)\ket{0}+e^{i\beta}\sin(\alpha)\ket{1},
\]
the measurement-basis rotation angles \((\theta,\phi)\),
\[
W=
\begin{bmatrix}
\cos(\theta/2) & e^{-i\phi}\sin(\theta/2)\\
e^{i\phi}\sin(\theta/2) & -\cos(\theta/2)
\end{bmatrix},
\]
and the evolution time \(t\), chosen in the range
\[
t\in[0,2\pi/\Delta_{\min}],
\]
where \(\Delta_{\min}\) is the minimum energy gap among the candidate Hamiltonians.

The optimized algorithm, denoted **QLE\_optimization**, dynamically selects the five-parameter control vector \((\alpha,\beta,\theta,\phi,t)\) at every iteration. This is a substantive departure from earlier QLE practice, in which the initial state is often fixed and \(t\) is chosen by the Particle Guess Heuristic (PGH). The point of the revised formulation is that discrimination power depends jointly on state preparation, measurement basis, and evolution time, so optimizing only one control leaves information on the table [2509.00825].

## 3. Information-theoretic formulation and adaptive control

The defining conceptual move in optimal QLE is to recast each iteration as a single-query oracle-identification problem. The hidden variable is the Hamiltonian label \(F\), distributed according to the current candidate probabilities \(p_f\), and the observable is the measurement outcome \(Y\). The joint classical-quantum state is written as
\[
\rho_{FY}=\sum_{f\in\mathscr{F}} p_f\,\ket{f}\bra{f}\otimes WU_f\ket{\psi_1}\bra{\psi_1}U_f^\dagger W^\dagger.
\]
The experiment-design problem is then to choose the controls so as to maximize the mutual information \(I(F;Y)\) [2509.00825].

The paper expresses the mutual information as
\[
I(F;Y)=S(\rho_Y)-D_Y(\rho_{FY};Z^{\otimes n}),
\]
with the basis-dependent discord-like quantity
\[
D_Y(\rho_{FY};Z^{\otimes n})=S(\rho_Y)-S(\rho_{FY})+S(\rho_{FY}\mid Z^{\otimes n}).
\]
Because the candidate probabilities \(p_f\) are fixed within a given iteration, \(S(\rho_{FY})=H(F)\) is constant for that step. Maximizing \(I(F;Y)\) is therefore equivalent to minimizing the conditional von Neumann entropy, which becomes the cost function actually optimized in the algorithm. Operationally, the selected controls are those that make the post-measurement state maximally informative about the true Hamiltonian [2509.00825].

The search over controls is performed by **simulated annealing**. Starting from a randomly initialized parameter vector
\[
\mathbf{x}=[t,\theta,\phi,\alpha,\beta],
\]
the algorithm generates neighbors according to
\[
\mathbf{x}_{i+1}=\mathbf{x}_i+T\cdot \mathrm{rand}([-1,1]^5)\circ (\mathrm{range}/2),
\]
where \(\circ\) denotes the Hadamard product. Candidates with lower cost are always accepted; higher-cost candidates are accepted with probability
\[
\exp\!\left(\frac{\mathrm{cost}_{\text{old}}-\mathrm{cost}_{\text{new}}}{T}\right).
\]
The temperature is then reduced geometrically,
\[
T\leftarrow \alpha T,\qquad \alpha\in(0,1),
\]
with \(\alpha=0.9\) in the reported implementation. The role of annealing is to explore broadly at high temperature and concentrate on low-entropy regions later, mitigating local-minimum trapping in a highly nonconvex control landscape [2509.00825].

A common simplification is to describe QLE as merely Bayesian Hamiltonian identification. That is accurate but incomplete. In the optimal formulation, the distinguishing feature is not Bayesian updating alone; it is adaptive quantum experiment design driven by an information-theoretic objective.

## 4. Convergence behavior and reported benchmarks

The most direct evidence for the revised formulation comes from the benchmark problems reported in [2509.00825]. In a four-Hamiltonian test set,
\[
H_1=\sigma_x,\qquad H_2=2\sigma_x,\qquad H_3=\sigma_z,\qquad H_4=2\sigma_z,
\]
the best static baseline required about **144 iterations on average** to make all four cases converge, whereas the optimized QLE required only **9 iterations**. When the success threshold was tightened from \(0.99\) to \(0.9999\), the improvement became even more pronounced.

The same study also tested a control variant in which a discrete grid search was repeated at every iteration. That version converged in about **10 iterations**, very close to the **9 iterations** of the continuous simulated-annealing version. The explicit interpretation given in the paper is that the speedup comes mainly from the optimization strategy itself, not merely from the availability of a larger continuous parameter space [2509.00825].

A harder six-Hamiltonian benchmark, containing mixed Pauli terms and asymmetric shifts, sharpened the comparison. The original QLE failed to converge reliably, whereas the optimized version succeeded for all six Hamiltonians and required only about **4–5 iterations on average** per Hamiltonian [2509.00825]. The paper therefore frames the principal gain not as asymptotic complexity theory but as practical reduction in the number of adaptive experimental rounds needed to drive the posterior toward a single Hamiltonian hypothesis.

Two points follow from these benchmarks. First, QLE is highly parameter-sensitive if controls are chosen heuristically. Second, the large empirical improvement is consistent with the information-theoretic interpretation: the algorithm is faster when each query is selected to maximize expected information gain about the hidden Hamiltonian label.

## 5. Relation to other likelihood-based quantum inference methods

QLE belongs to a broader family of likelihood-based quantum inference procedures, but its target object, data model, and adaptation mechanism are specific. Related methods differ by whether they infer Hamiltonian parameters, amplitudes, density matrices, or latent states in continuously monitored or hidden-state models.

| Setting | Likelihood object | Representative papers |
|---|---|---|
| Hamiltonian learning from thermal or projector data | Negative log-likelihood over Gibbs states | [2212.13718] |
| Amplitude estimation | Binomial likelihood over amplification depths | [2006.16223], [2209.03321] |
| Density-matrix inference and tomography | Quantum relative entropy or Gaussian likelihood with physical-state constraints | [1904.05903], [1604.04691] |
| Continuous monitoring | \(L_t=\mathrm{Tr}(\tilde\rho_t)\) from an unnormalized SME | [1212.5700] |
| Quantum walks | Exact displacement/return likelihood \(p^{(k)}(d|\theta)\) | [2202.11846] |
| Quantum hidden-state models | Sequence likelihood from Kraus-operator evolution | [2507.21337] |

In maximum-likelihood Hamiltonian learning for many-body systems, the candidate model is updated directly in Hamiltonian space rather than by maintaining a discrete hypothesis distribution. The update
\[
\hat H_{k+1}=\hat H_k-\gamma \hat R_k
\]
is constructed so that the negative log-likelihood decreases, at least to linear order in the step size, while the state remains Gibbsian throughout the optimization [2212.13718]. This is close in spirit to QLE, but it is not the same algorithmic object: the likelihood is defined from measurement frequencies on thermal or ground states, and the optimization variable is the Hamiltonian itself rather than a control policy over sequential experiments.

In amplitude estimation, likelihood-based quantum inference takes the form of maximum-likelihood fitting of measurement counts produced by repeated amplitude amplification. In the noisy formulation of [2006.16223], the model jointly estimates the target amplitude \(a\) and a depolarization parameter \(\kappa\); in the improved MLQAE analysis of [2209.03321], the likelihood is built from counts gathered at different Grover depths, and exceptional amplitude values together with a depth-jittering fix are analyzed explicitly. These methods share QLE’s hybrid structure—quantum sampling plus classical likelihood optimization—but they estimate an amplitude parameter rather than a Hamiltonian label.

Density-matrix inference provides another neighboring lineage. "Inferring the quantum density matrix with machine learning" formulates **Quantum Maximum Likelihood (QML)** and **Quantum Variational Inference (QVI)** in terms of the quantum relative entropy, using a spectral variational family and optional **quantum flows** to learn a thermal density matrix and thereby recover eigenstates and eigenvalues [1904.05903]. In experimental NMR tomography, maximum-likelihood estimation is used to enforce physicality by parameterizing the state as
\[
\rho=T^\dagger T,\qquad \mathrm{Tr}(T^\dagger T)=1,
\]
thereby guaranteeing positivity and normalization throughout reconstruction [1604.04691]. These approaches are likelihood-based quantum estimation, but their target is the density operator, not an unknown Hamiltonian chosen from a candidate set.

Likelihood generation from continuously monitored systems proceeds differently again. There, the stochastic master equation conditioned on the measurement record is simultaneously an equation for the conditional quantum state and for the likelihood of the record under candidate parameters. The unnormalized conditional state \(\tilde\rho_t\) yields
\[
L(D|\theta)\propto \mathrm{Tr}\,\tilde\rho_t,
\]
which can be used directly in Bayesian updating or MCMC posterior sampling [1212.5700]. The shared structural feature with QLE is the centrality of likelihoods derived from quantum dynamics, but the data are full time-resolved records rather than discrete adaptive interrogation rounds.

## 6. Generalizations, scope, and nomenclature

The 2025 optimal-QLE formulation explicitly generalizes beyond Hamiltonian identification. The paper describes a class of quantum learning problems consisting of an initial state \(\rho_i(\mathbf{c}_i,\mathbf{h}_i)\), a CPTP map \(\mathcal{P}(\mathbf{c}_p,\mathbf{h}_p)\), an attribute \(a(\mathcal{P}\{\rho_i\},\mathbf{c}_a,\mathbf{h}_a)\), and hidden parameters \(\mathbf{h}\) drawn from a prior distribution. In that setting, the controlled parameters \((\mathbf{c}_i,\mathbf{c}_p,\mathbf{c}_a)\) are chosen to maximize
\[
I\!\left(a(\cdot); f(\mathbf{h}_i,\mathbf{h}_p,\mathbf{h}_a)\right),
\]
with Bayesian or likelihood-based inference updating the hidden-parameter distribution after each query [2509.00825]. This makes QLE a specific instance of a more general prepare–measure–update–optimize paradigm.

The same paper notes that the framework can be extended to noisy scenarios and to multi-qubit systems, while also observing that the parameter search becomes more difficult there and may benefit from quantum annealing in place of classical simulated annealing [2509.00825]. A plausible implication is that the principal bottleneck in scaling QLE is not Bayesian updating itself but adaptive control optimization in larger experimental design spaces.

The term **QLE** is, however, acronymically overloaded. In mathematical physics and probability theory, it also denotes **quantum Loewner evolution**, especially **QLE \((8/3,0)\)**. In that usage, QLE is a growth process on a doubly marked \(\sqrt{8/3}\)-Liouville quantum gravity sphere, constructed by repeatedly growing small pieces of an \(\mathrm{SLE}_6\)-type exploration, resampling the tip uniformly along the quantum boundary, and taking a subsequential limit as \(\delta\to 0\). That process is used to define a metric on a dense subset of the surface and to prove symmetry and strict triangle inequality for the resulting distance on sampled points [1507.00719]. This meaning is unrelated to Hamiltonian learning except for the shared initials.

There are also broader statistical uses of the language of quantum likelihood. One line of work develops a quantum framework for likelihood ratios using Hilbert-space encodings of contingency tables and overlap structure interpreted through non-orthonormal inner products [1508.00936]. Another studies maximum-likelihood estimation in quantum probability on orthomodular posets and Greechie diagrams, where likelihoods
\[
L(p)=\prod_{x\in A} p(x)^{n(x)}
\]
are optimized subject to contextual normalization constraints on operations [2401.11253]. These usages widen the semantic field of “quantum likelihood estimation,” but they are conceptually distinct from the Hamiltonian-learning algorithm now commonly denoted QLE in quantum computing.

Taken together, the literature supports a narrow and a broad sense of the term. In the narrow and technically dominant quantum-computing sense, QLE is an adaptive Bayesian Hamiltonian-learning protocol accelerated by quantum-generated likelihood information and, in its current optimal form, by mutual-information-based experiment design [2509.00825]. In the broader sense, it belongs to a wider class of likelihood-based quantum inference procedures that use quantum dynamics, quantum probability, or quantum-hidden-state representations to make latent physical parameters statistically identifiable.

Source: https://www.emergentmind.com/topics/quantum-likelihood-estimation-qle