pith. sign in

arxiv: 2510.07439 · v3 · submitted 2025-10-08 · 🪐 quant-ph · cs.DS

Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra

Pith reviewed 2026-05-18 09:02 UTC · model grok-4.3

classification 🪐 quant-ph cs.DS
keywords quantum algorithmeigenvalue multiplicityspectral analysisquantum filteringmany-body systemstopological orderphase transitionsenergy clusters
0
0 comments X

The pith

A quantum algorithm identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physical assumptions.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper presents QFAMES, a quantum method designed to extract fine details about the eigenvalues of many-body Hamiltonians and how many times each value repeats. It targets groups of nearly equal dominant eigenvalues rather than resolving every individual level. This focus on clusters, which often appear in real physical systems, lets the procedure avoid the general computational hardness that blocks exact multiplicity counting on quantum computers. The same filtering step also supports direct estimation of observable averages restricted to states inside a chosen energy window. Tests on the transverse-field Ising model for phase identification and the toric code for degeneracy counting illustrate the approach in practice.

Core claim

QFAMES efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. The method supplies rigorous theoretical guarantees and shows advantages in sample complexity compared with prior subspace techniques.

What carries the argument

QFAMES, a quantum filtering procedure that isolates and counts multiplicities inside targeted energy clusters rather than across the full spectrum.

If this is right

  • Observable expectation values can be estimated inside specific energy clusters to characterize quantum phase transitions.
  • Ground-state degeneracy can be extracted for topologically ordered phases such as the two-dimensional toric code.
  • The procedure generalizes directly to mixed initial states.
  • Sample complexity improves over existing subspace-based spectral methods while retaining rigorous guarantees.
  • Numerical demonstrations confirm utility for both the transverse-field Ising model and topologically ordered systems.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The filtering idea could be combined with existing ground-state preparation routines to focus resources on a chosen energy window.
  • Application to noisy intermediate-scale devices would test whether the cluster assumption survives realistic decoherence.
  • Similar filtering steps might resolve degeneracies in other spectral problems such as those arising in open quantum systems.

Load-bearing premise

The assumption that dominant eigenvalues in physical Hamiltonians form tight clusters is what makes accurate multiplicity extraction feasible.

What would settle it

A concrete Hamiltonian whose dominant eigenvalues are well-separated rather than clustered, yet on which QFAMES still reports incorrect multiplicity counts, would disprove the central claim.

Figures

Figures reproduced from arXiv: 2510.07439 by Lin Lin, Ruizhe Zhang, Yilun Yang, Zhiyan Ding.

Figure 1
Figure 1. Figure 1: Illustration of the QFAMES algorithm. Given the Hamiltonian and two sets of initial states prepared by [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Diagram of the QFAMES algorithm for post-processing the 3-tensor generated from quantum data of the [PITH_FULL_IMAGE:figures/full_fig_p007_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Illustration of the QFAMES algorithm on the illustrative example. The dashed vertical lines are the [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: QFAMES (Algorithm 1) vs QMEGS in the illustrative example as a function of (left): max evolution time Tmax = σT, and (right): total evolution time Ttotal. The dashed line stands for the fitted error proportional to 1/T, as predicted by Eq. (27). A. Illustrative example To illustrate the effectiveness of QFAMES, we consider a simple Hamiltonian with three eigenvalues: λ0 = λ1 = 0 and λ2 = 0.1. There are onl… view at source ↗
Figure 5
Figure 5. Figure 5: Illustration of the QFAMES algorithm on the TFIM model. The upper row shows the Frobenius norm of [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Illustration of Algorithm 2 on the TFIM model. The observable is S z = 1 2L PL i=1 Zi . In the manifold of near-ground-state cluster, the possible ranges of observable expectation values are shown by the gray lines. A transition from ferromagnetic to paramagnetic phase can be observed when g increases from 0.5 to 1.5. where Zi and Xi are the Pauli operators acting on the i-th qubit and g is the coupling co… view at source ↗
Figure 7
Figure 7. Figure 7: 2D Toric codes on a 2-by-4 lattice with (left) torus, and (right) cylinder boundary conditions. Dashed lines stand for periodic boundary conditions in the corresponding directions. The models consist of 16 and 18 qubits, respectively, represented by the small circles. The Toric code ground state is known to exhibit Z2 topological order, which has been realized on a quantum computer based on the prior knowl… view at source ↗
Figure 8
Figure 8. Figure 8: Illustration of GSD estimation for 2D Toric code models with [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Hadamard test circuit. W is either the identity or the phase gate S † . Here INIT |0⟩ = |ϕ⟩. For each n, the Hadamard test circuit guarantees that E[Zn] = ⟨ϕ|e −iHtn |ϕ⟩ = X m∈[M] |⟨ϕ|Em⟩|2 | {z } :=pm ·e −iλmtn . (A1) Thus, the dataset {(tn, Zn)}n∈[N] corresponds to the unbiased noisy samples of the Fourier signal x(t) = X m∈[M] pm exp(−iλmt). (A2) Estimating the dominant eigenvalues {λm}m∈D from these sa… view at source ↗
read the original abstract

Fine-grained spectral properties of quantum Hamiltonians, including both eigenvalues and their multiplicities, provide useful information for characterizing many-body quantum systems as well as for understanding phenomena such as topological order. Extracting such information with small additive error is $\#\textsf{BQP}$-complete in the worst case. In this work, we introduce QFAMES (Quantum Filtering and Analysis of Multiplicities in Eigenvalue Spectra), a quantum algorithm that efficiently identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions, which allows us to bypass worst-case complexity barriers. QFAMES also enables the estimation of observable expectation values within targeted energy clusters, providing a powerful tool for studying quantum phase transitions and other physical properties. We validate the effectiveness of QFAMES through numerical demonstrations, including its applications to characterizing quantum phases in the transverse-field Ising model and estimating the ground-state degeneracy of a topologically ordered phase in the two-dimensional toric code model. We also generalize QFAMES to the setting of mixed initial states. Our approach offers rigorous theoretical guarantees and significant advantages over existing subspace-based quantum spectral analysis methods, particularly in terms of the sample complexity and the ability to resolve degeneracies.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript introduces QFAMES, a quantum algorithm that identifies clusters of closely spaced dominant eigenvalues and determines their multiplicities under physically motivated assumptions on the spectrum. This is claimed to bypass worst-case #BQP-completeness of fine-grained spectral analysis while also enabling estimation of observable expectation values within targeted energy clusters. Numerical validation is provided on the transverse-field Ising model for quantum phases and the 2D toric code for ground-state degeneracy, with a generalization to mixed initial states and claims of rigorous theoretical guarantees plus advantages in sample complexity over subspace methods.

Significance. If the efficiency claims and complexity bypass hold under the stated assumptions, the work would provide a practically relevant tool for quantum many-body spectral analysis, particularly for studying phase transitions and topological order where clusters of dominant eigenvalues are physically natural. The numerical demonstrations on standard models like Ising and toric code add concrete support for applicability, and the mixed-state generalization broadens the scope.

major comments (2)
  1. [Abstract] Abstract: The central claim that physically motivated assumptions on clusters of closely spaced dominant eigenvalues suffice to bypass #BQP-completeness requires explicit quantitative conditions (e.g., lower bounds on inter-cluster gaps relative to 1/poly(n) or intra-cluster spacing thresholds that guarantee polynomial scaling of the filtering subroutine). No such thresholds or gap assumptions are stated, leaving open whether the algorithm remains efficient when gaps are sub-polynomial as noted in the hardness reduction for multiplicity counting.
  2. [Abstract] Abstract and claims of rigorous guarantees: The manuscript asserts 'rigorous theoretical guarantees' and efficiency under the assumptions, but the provided text contains no derivation details, error bounds, or explicit assumption statements for the quantum filtering and multiplicity estimation steps. This makes it impossible to verify that the approach avoids the #BQP-hardness reduction without additional exponential resources.
minor comments (1)
  1. [Abstract] The abstract mentions 'significant advantages over existing subspace-based quantum spectral analysis methods' in sample complexity and degeneracy resolution, but does not cite or compare against specific prior works (e.g., quantum phase estimation variants or subspace diagonalization algorithms); adding these references would improve context.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their detailed and constructive feedback on the abstract and the strength of our claims. We address each major comment below and will revise the manuscript to improve clarity and explicitness while preserving the core contributions.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The central claim that physically motivated assumptions on clusters of closely spaced dominant eigenvalues suffice to bypass #BQP-completeness requires explicit quantitative conditions (e.g., lower bounds on inter-cluster gaps relative to 1/poly(n) or intra-cluster spacing thresholds that guarantee polynomial scaling of the filtering subroutine). No such thresholds or gap assumptions are stated, leaving open whether the algorithm remains efficient when gaps are sub-polynomial as noted in the hardness reduction for multiplicity counting.

    Authors: We agree that the abstract would benefit from explicit quantitative gap conditions to strengthen the complexity claim. The main text (Section 3 and Theorem 1) assumes inter-cluster gaps of at least 1/poly(n) and intra-cluster eigenvalue spacings sufficiently small to ensure the filtering operator achieves polynomial sample complexity via standard quantum signal processing techniques. These conditions are physically motivated for dominant clusters in many-body spectra and suffice to bypass the #BQP-hardness reduction, which requires resolving multiplicities across sub-polynomial gaps. We will revise the abstract to state these thresholds explicitly. revision: yes

  2. Referee: [Abstract] Abstract and claims of rigorous guarantees: The manuscript asserts 'rigorous theoretical guarantees' and efficiency under the assumptions, but the provided text contains no derivation details, error bounds, or explicit assumption statements for the quantum filtering and multiplicity estimation steps. This makes it impossible to verify that the approach avoids the #BQP-hardness reduction without additional exponential resources.

    Authors: The full manuscript contains the derivations, error bounds (additive error O(1/poly(n)) with high probability), and explicit assumption statements in Sections 3–4 and the appendix proofs for the filtering and multiplicity estimation subroutines. The abstract summarizes these results. To address the concern about verifiability from the abstract alone, we will add a concise statement of the key assumptions and polynomial scaling guarantee in the revised abstract, ensuring readers can see that no exponential overhead is introduced beyond the stated conditions. revision: yes

Circularity Check

0 steps flagged

No significant circularity; QFAMES claims rest on external physical assumptions and new algorithm design

full rationale

The paper introduces QFAMES as a quantum algorithm for identifying eigenvalue clusters and multiplicities under physically motivated assumptions on spectra (e.g., clusters of closely spaced dominant eigenvalues). This allows bypassing worst-case #BQP-completeness for fine-grained analysis. The abstract and description tie efficiency and guarantees directly to these external assumptions and the algorithm's filtering steps, without reducing claims to fitted parameters, self-citations, or internal definitions by construction. Numerical demonstrations on Ising and toric code models are presented as validation, not as the source of the core claims. No load-bearing self-citation chains or ansatz smuggling are indicated. The derivation is self-contained against external benchmarks, consistent with a standard algorithm paper.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on unspecified physically motivated assumptions about eigenvalue clustering that are invoked to evade #BQP-completeness; no free parameters or invented entities are mentioned in the abstract.

axioms (1)
  • domain assumption Physically motivated assumptions on the structure of eigenvalue spectra (clusters of closely spaced dominant eigenvalues) hold for the target Hamiltonians.
    Abstract states these assumptions allow bypassing worst-case #BQP-completeness for multiplicity determination.

pith-pipeline@v0.9.0 · 5744 in / 1279 out tokens · 30549 ms · 2026-05-18T09:02:27.909596+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

94 extracted references · 94 canonical work pages · 3 internal anchors

  1. [1]

    there are no tail eigenvalues, i.e.,ptail = 0. To achieveϵ O-accuracy in estimating eigenvalues ofODi, it suffices to choose Tmax = eO(∆−1), N= eΩ(ϵ−2 O ),(42) and the total evolution time Ttotal = eO(N LRTmax) = eO(∆−1ϵ−2 O ).(43) It is instructive to contrast the generalized eigenvalue problem in Eq. (36) with that used in quantum subspace methods for e...

  2. [2]

    exactly allows to measureZl,r,n =⟨ϕ l|e−iHt n |ψr⟩with different left and right initial states|ϕl⟩and|ψ r⟩

  3. [3]

    can produce the data fortin the time interval[0, σT]in a single execution of the algorithm

  4. [4]

    Below we will briefly introduce this algorithm

    can be easily generalized to the cases involving multiple time evolutions, thereby allowing access to the 4-tensor Z O l,r,n,n′ in the observable version of QFAMES. Below we will briefly introduce this algorithm. A key observation is that the main difficulty of measuringZl,r,n = r(tn)eiϕ(tn) lies in the phaseϕ(tn)rather than the absolute valuer(tn), and t...

  5. [5]

    To resolve the degeneracy atλ⋆ 0, it is essential to utilize the off-diagonal entries of the data tensor, which encode important information about the spectral structure

    Moreover, because the overlaps withλ2 are not sufficiently dominant compared to those withλ0, λ1, single-state approaches might not reliably estimate the location ofλ2. To resolve the degeneracy atλ⋆ 0, it is essential to utilize the off-diagonal entries of the data tensor, which encode important information about the spectral structure. We can check that...

  6. [6]

    The first point allows for a more accurate estimation of the location ofλ2, while the second provides insight into the multiplicity ofλ ⋆ 1

    The degeneracy of dominant eigenvalueλ⋆ 0 must be at least two, since the signal corresponding to the frequency λ≈0vanishes in this off-diagonal entry. The first point allows for a more accurate estimation of the location ofλ2, while the second provides insight into the multiplicity ofλ ⋆ 1. Now we perform numerical simulations to demonstrate the main ide...

  7. [7]

    amplified

    Although the diagonal landscapes are the same in Fig. 3(a), if we focus on the off diagonal landscapes, it is evident that there exists one dominant eigenvalue nearλ2 = 0.1. This can be more systematically observed in the plot of the Frobenius norm∥G(θ)∥2 F in Fig. 3(b), which clearly reveals two distinct dominant eigenvalues, including an “amplified” pea...

  8. [8]

    We randomly sample vectors from the Haar distribution and apply imaginary-time evolutionexp(−βH)to enhance their overlap with the ground-state subspace

  9. [9]

    post-Kitaev

    We then apply QFAMES to these boosted vectors to estimate the multiplicity of the cluster corresponding to the smallest eigenvalue. We find that in the first step of the algorithm, it is necessary to generate a sufficient number of initial states such that the ground-state subspace is contained within their linear span. This condition is essential for cor...

  10. [10]

    Generate a sequence of real numbers{tn}n∈[N] as the Hamiltonian evolution times

  11. [11]

    9) to the initial state|ϕ⟩for each timet n, and collect the measurement outcomes to form a dataset{(tn, Zn)}n∈[N], whereZ n ∈ {±1±i}

    Apply Hadamard test circuit (Fig. 9) to the initial state|ϕ⟩for each timet n, and collect the measurement outcomes to form a dataset{(tn, Zn)}n∈[N], whereZ n ∈ {±1±i}

  12. [12]

    |0⟩ H • W H |0⟩ INIT e−itnH Figure 9: Hadamard test circuit.Wis either the identity or the phase gateS†

    Perform classical post-processing on the dataset{(tn, Zn)}n∈[N] to estimate the dominant eigenvalues{λm}m∈D. |0⟩ H • W H |0⟩ INIT e−itnH Figure 9: Hadamard test circuit.Wis either the identity or the phase gateS†. HereINIT|0⟩=|ϕ⟩. For eachn, the Hadamard test circuit guarantees that E[Zn] =⟨ϕ|e −iHt n |ϕ⟩= X m∈[M] |⟨ϕ|Em⟩|2 | {z } :=pm ·e−iλmtn .(A1) Thus...

  13. [13]

    (B2) is applied in Eq

    In eigenvalue estimation, the condition Eq. (B2) is applied in Eq. (C25) to ensure that Φ:,D(πθ) Ψ:,D(πθ) † 2 F and∥G(θ≈λ ⋆ πi)∥F admit sufficiently large lower bounds. This guarantees that the dominant eigenvalues can be detected by identifying the peaks of∥G(θ)∥F

  14. [14]

    This property is then quantitatively exploited in Eq

    In multiplicity estimation, the condition ensures that the overlap matrix Φ:,Dπθ Ψ:,Dπθ † has rank equal to the corresponding degeneracy. This property is then quantitatively exploited in Eq. (C50) to establish a lower bound on the singular values ofG(θ≈θ⋆ πi). Appendix C: Rigorous version of Theorem IV.1 and proof In this section, we introduce and prove ...

  15. [15]

    +E(θ ⋆ 1)∥F ≤ 5ptail 4 , B eθ⋆ 1 +C eθ⋆ 1 +E eθ⋆ 1 F ≤ 5ptail 4 .(C34) Then, we get that W2(θ) = Tr(G†(θ)G(θ)) = Tr(A†(θ)A(θ)) + 2Re Tr(A†(θ)(B(θ) +C(θ) +E(θ))) + Tr((B†(θ) +C †(θ) +E †(θ))(B(θ) +C(θ) +E(θ))). (C35) We first deal with the first term, the last term, and then the second term: •For the first term, we have A(θ) = Φ:,Dπ1 ·diag exp(−(θ−λ ⋆ π1)2...

  16. [16]

    +E(θ ⋆ 1))))≤ 9 4 p2 tail, Tr((B†(eθ⋆

  17. [17]

    (C39) 30 •For the second term, we have 2Re Tr(A†(θ⋆ 1)(B(θ ⋆

    +E(eθ⋆ 1)))≤ 9 4 p2 tail . (C39) 30 •For the second term, we have 2Re Tr(A†(θ⋆ 1)(B(θ ⋆

  18. [18]

    +E(θ ⋆ 1)) −2Re Tr(A†(eθ⋆ 1)(B(eθ⋆

  19. [19]

    +E(eθ⋆ 1))) ≤2 A(θ⋆ 1)−A( eθ⋆ 1) F ∥(B(θ ⋆

  20. [20]

    +E(θ ⋆ 1))∥F + 2 A(eθ⋆ 1) F (B(θ ⋆

  21. [21]

    (C40) Similar to the calculation in Eq

    +E(eθ⋆ 1)) F . (C40) Similar to the calculation in Eq. (C19), we obtain from Eq. (C38) that A(θ⋆ 1)−A( eθ⋆ 1) F ≤T θ⋆ 1 −eθ⋆ 1 Φ:,Dπ1 Ψ:,Dπ1 † F + 2 √ KLRδ·T ≤T θ⋆ 1 −eθ⋆ 1 Φ:,Dπ1 Ψ:,Dπ1 † F + ptail 2 , (C41) where we useδ·T=O(p tail/ √ KLR)in the last inequality. According to Theorem C.6 Eq. (C65) and Theo- rem C.5 Eq. (C58), we have ∥B(θ)−B(θ ′)∥F ≤p ta...

  22. [22]

    +E(θ ⋆ 1))) −2Re Tr(A†(eθ⋆ 1)(B(eθ⋆

  23. [23]

    (C16)) in the second inequality andq=O(p tail/((1 +σ) √ KLR))andδ=O(p tail/(T(1 +σ) √ KLR))in the last inequality

    +E(eθ⋆ 1)))) ≤2 T θ⋆ 1 −eθ⋆ 1 Φ:,Dπ1 Ψ:,Dπ1 † F + ptail 2 · 5 4 ptail + 2· 5 4 Φ:,Dπ1 Ψ:,Dπ1 † F · 2ptailσT θ⋆ 1 −eθ⋆ 1 + p2 tail 4 √ KLR ≤(3 + 5σ) Φ:,Dπ1 Ψ:,Dπ1 † F T θ⋆ 1 −eθ⋆ 1 ptail + 15 8 p2 tail ≤(3 + 5σ) Φ:,Dπ1 Ψ:,Dπ1 † F T θ⋆ 1 −λ ⋆ π1 + q T +δ ptail + 15 8 p2 tail ≤(3 + 5σ) Φ:,Dπ1 Ψ:,Dπ1 † F T θ⋆ 1 −λ ⋆ π1 ptail + 23 8 p2 tail , (C44) where we us...

  24. [24]

    (C48) 32 With probability at least1−η, we have ∥B(θ ⋆

    = Φ :,Dπ1 ·diag exp(−(θ⋆ 1 −λ m)2T 2) m∈Dπ1 · Ψ:,Dπ1 † | {z } :=A(θ⋆ 1) + Φ:,Dc ·diag exp(−(θ⋆ 1 −λ m)2T 2) m∈Dc ·(Ψ :,Dc)† | {z } :=B(θ ⋆ 1) + X i̸=π1 Φ:,Di ·diag exp(−(θ⋆ 1 −λ m)2T 2) m∈Di ·(Ψ :,Di)† | {z } :=C(θ ⋆ 1) +E(θ ⋆ 1). (C48) 32 With probability at least1−η, we have ∥B(θ ⋆

  25. [25]

    +E(θ ⋆ 1)∥2 ≤ ∥B(θ ⋆

  26. [26]

    +E(θ ⋆ 1)∥F ≤2p tail .(C49) Furthermore, for1≤k≤ |D π1 |, Eq. (C5) implies that sk (A(θ⋆ 1)) = exp −O((1 +σ)(p tail/pmin) +δ·T) 2 λk Ψ:,Dπ1 Φ† :,Dπ1 Φ:,Dπ1 Ψ:,Dπ1 † 1/2 ≥exp −O((1 +σ)(p tail/pmin) +δ·T) 2 λmin Φ† :,Dπ1 Φ:,Dπ1 1/2 λk Ψ:,Dπ1 Ψ:,Dπ1 † 1/2 ≥exp −O((1 +σ)(p tail/pmin) +δ·T) 2 λmin Φ† :,Dπ1 Φ:,Dπ1 1/2 λmin Ψ† :,Dπ1 Ψ:,Dπ1 1/2 ≥exp −O((1 +σ)(p t...

  27. [27]

    Let{t n}to be i.i.d sampled from the truncated Gaussiana T (t)defined as Eq.(16)

    Useful lemmas Lemma C.5.Givenβ >0. Let{t n}to be i.i.d sampled from the truncated Gaussiana T (t)defined as Eq.(16). Define E(θ) := 1 N NX n=1 Zn exp(iθtn)−Φ i,: ·diag {exp(−(θ−λ m)2T 2}m∈[M] ·(Ψ †):,j ∈C L×R .(C56) LetΘ :={θ j}J j=1 ∪ {λm}m∈D. Ifσ= Ω log1/2 √ LR/β as in Eq.(16)andN= Ω LR β2 log (J+|D|) LR η , we have Pr max θ∈Θ ∥E(θ)∥ F ≤β ≥1−η .(C57) an...

  28. [28]

    Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev

    X.-G. Wen, Colloquium: Zoo of quantum-topological phases of matter, Rev. Mod. Phys.89, 041004 (2017)

  29. [29]

    Nayak, S

    C. Nayak, S. H. Simon, A. Stern, M. Freedman, and S. Das Sarma, Non-abelian anyons and topological quantum compu- tation, Rev. Mod. Phys.80, 1083 (2008)

  30. [30]

    Brown, S

    B. Brown, S. T. Flammia, and N. Schuch, Computational difficulty of computing the density of states, Phys. Rev. Lett. 107, 040501 (2011)

  31. [31]

    F. G. S. L. Brandao, Entanglement theory and the quantum simulation of many-body physics (2008), arXiv:0810.0026 [quant-ph]

  32. [32]

    Gyurik, C

    C. Gyurik, C. Cade, and V. Dunjko, Towards quantum advantage via topological data analysis, Quantum6, 855 (2022)

  33. [33]

    Y. Yang, A. Christianen, M. C. Bañuls, D. S. Wild, and J. I. Cirac, Phase-sensitive quantum measurement without controlled operations, Phys. Rev. Lett.132, 220601 (2024)

  34. [34]

    Clinton, T

    L. Clinton, T. S. Cubitt, R. Garcia-Patron, A. Montanaro, S. Stanisic, and M. Stroeks, Quantum phase estimation without controlled unitaries (2024), arXiv:2410.21517 [quant-ph]

  35. [35]

    X. Wang, L. Xiong, X. Cai, and X. Yuan, Computingn-time correlation functions without ancilla qubits (2025), arXiv:2504.12975 [quant-ph]

  36. [36]

    B. F. Schiffer, D. S. Wild, N. Maskara, M. D. Lukin, and J. I. Cirac, Hardware-efficient quantum phase estimation via local control (2025), arXiv:2506.18765 [quant-ph]. 37

  37. [37]

    Yi and C

    C. Yi and C. Zhou, Ancilla-free quantum protocol for thermal Green’s functions (2025), arXiv:2509.03288 [quant-ph]

  38. [38]

    P. Zeng, J. Sun, and X. Yuan, Universal quantum algorithmic cooling on a quantum computer (2022), arXiv:2109.15304 [quant-ph]

  39. [39]

    Zhang, G

    R. Zhang, G. Wang, and P. Johnson, Computing ground state properties with early fault-tolerant quantum computers, Quantum6, 761 (2022)

  40. [40]

    J. Sun, P. Zeng, T. Gur, and M. S. Kim, High-precision and low-depth eigenstate property estimation: theory and resource estimation (2024), arXiv:2406.04307 [quant-ph]

  41. [41]

    Sun and P

    J. Sun and P. Zeng, Randomised composite linear-combination-of-unitaries: its role in quantum simulation and observable estimation (2025), arXiv:2506.15658 [quant-ph]

  42. [42]

    A. Y. Kitaev, Quantum measurements and the abelian stabilizer problem (1995), arXiv:quant-ph/9511026 [quant-ph]

  43. [43]

    A. Y. Kitaev, A. H. Shen, and M. N. Vyalyi,Classical and Quantum Computation(American Mathematical Society, USA, 2002)

  44. [44]

    M.Dobšíček, G.Johansson, V.Shumeiko,andG.Wendin,Arbitraryaccuracyiterativequantumphaseestimationalgorithm using a single ancillary qubit: A two-qubit benchmark, Phys. Rev. A76, 030306 (2007)

  45. [45]

    Wiebe and C

    N. Wiebe and C. Granade, Efficient Bayesian phase estimation, Phys. Rev. Lett.117, 010503 (2016)

  46. [46]

    R. D. Somma, Quantum eigenvalue estimation via time series analysis, New J. Phys.21, 123025 (2019)

  47. [47]

    T. E. O’Brien, B. Tarasinski, and B. M. Terhal, Quantum phase estimation of multiple eigenvalues for small-scale (noisy) experiments, New J. Phys.21, 023022 (2019)

  48. [48]

    Ding and L

    Z. Ding and L. Lin, Even shorter quantum circuit for phase estimation on early fault-tolerant quantum computers with applications to ground-state energy estimation, PRX Quantum4, 020331 (2023)

  49. [49]

    Z.DingandL.Lin,Simultaneousestimationofmultipleeigenvalueswithshort-depthquantumcircuitonearlyfault-tolerant quantum computers, Quantum7, 1136 (2023)

  50. [50]

    Z. Ding, H. Li, L. Lin, H. Ni, L. Ying, and R. Zhang, Quantum multiple eigenvalue gaussian filtered search: an efficient and versatile quantum phase estimation method, Quantum8, 1487 (2024)

  51. [51]

    Z. Ding, Y. Dong, Y. Tong, and L. Lin, Robust ground-state energy estimation under depolarizing noise (2024), arXiv:2307.11257 [quant-ph]

  52. [52]

    C. Yi, C. Zhou, and J. Takahashi, Quantum phase estimation by compressed sensing, Quantum8, 1579 (2024)

  53. [53]

    Castaldo and S

    D. Castaldo and S. Corni, Heisenberg limited multiple eigenvalue estimation via off-the-grid compressed sensing (2025), arXiv:2507.12438 [quant-ph]

  54. [54]

    H. Ni, H. Li, and L. Ying, On low-depth algorithms for quantum phase estimation, Quantum7, 1165 (2023)

  55. [55]

    G. Wang, D. S. França, R. Zhang, S. Zhu, and P. D. Johnson, Quantum algorithm for ground state energy estimation using circuit depth with exponentially improved dependence on precision, Quantum7, 1167 (2023)

  56. [56]

    H. Li, H. Ni, and L. Ying, Adaptive low-depth quantum algorithms for robust multiple-phase estimation, Phys. Rev. A 108, 062408 (2023)

  57. [57]

    R. M. Parrish and P. L. McMahon, Quantum filter diagonalization: Quantum eigendecomposition without full quantum phase estimation, preprint 10.48550/arXiv.1909.08925 (2019)

  58. [58]

    W. J. Huggins, J. Lee, U. Baek, B. O’Gorman, and K. B. Whaley, A non-orthogonal variational quantum eigensolver, New J. Phys.22, 073009 (2020)

  59. [59]

    N. H. Stair, R. Huang, and F. A. Evangelista, A multireference quantum Krylov algorithm for strongly correlated electrons, J. Chem. Theory Comput.16, 2236 (2020)

  60. [60]

    Seki and S

    K. Seki and S. Yunoki, Quantum power method by a superposition of time-evolved states, PRX Quantum2, 1 (2021)

  61. [61]

    E. N. Epperly, L. Lin, and Y. Nakatsukasa, a theory of quantum subspace diagonalization, SIMAX43, 1263 (2022)

  62. [62]

    U. Baek, D. Hait, J. Shee, O. Leimkuhler, W. J. Huggins, T. F. Stetina, M. Head-Gordon, and K. B. Whaley, Say NO to optimization: A non-orthogonal quantum eigensolver, PRX Quantum4, 1 (2022)

  63. [63]

    Kirby, Analysis of quantum Krylov algorithms with errors, Quantum8, 1 (2024)

    W. Kirby, Analysis of quantum Krylov algorithms with errors, Quantum8, 1 (2024)

  64. [64]

    Yoshioka, M

    N. Yoshioka, M. Amico, W. Kirby, P. Jurcevic, A. Dutt, B. Fuller, S. Garion, H. Haas, I. Hamamura, A. Ivrii, R. Majumdar, Z. Minev, M. Motta, B. Pokharel, P. Rivero, K. Sharma, C. J. Wood, A. Javadi-Abhari, and A. Mezzacapo, Krylov diagonalization of large many-body Hamiltonians on a quantum processor, Nat. Commun.16, 1 (2025)

  65. [65]

    Peruzzo, J

    A. Peruzzo, J. McClean, P. Shadbolt, M.-H. Yung, X.-Q. Zhou, P. J. Love, A. Aspuru-Guzik, and J. L. O’Brien, A variational eigenvalue solver on a photonic quantum processor, Nat. Commun.5, 4213 (2014)

  66. [66]

    A Quantum Approximate Optimization Algorithm

    E. Farhi, J. Goldstone, and S. Gutmann, A quantum approximate optimization algorithm (2014), arXiv:1411.4028 [quant- ph]

  67. [67]

    Kandala, A

    A. Kandala, A. Mezzacapo, K. Temme, M. Takita, M. Brink, J. M. Chow, and J. M. Gambetta, Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets, Nature549, 242 (2017)

  68. [68]

    Motta, C

    M. Motta, C. Sun, A. T. Tan, M. J. O’Rourke, E. Ye, A. J. Minnich, F. G. Brandão, and G. K. L. Chan, Determining eigenstates and thermal states on a quantum computer using quantum imaginary time evolution, Nat. Phys.16, 205 (2020)

  69. [69]

    S.-H. Lin, R. Dilip, A. G. Green, A. Smith, and F. Pollmann, Real- and imaginary-time evolution with compressed quantum circuits, PRX Quantum2, 010342 (2021)

  70. [70]

    X. Yuan, S. Endo, Q. Zhao, Y. Li, and S. C. Benjamin, Theory of variational quantum simulation, Quantum3, 1 (2019)

  71. [71]

    Slepian, Prolate spheroidal wave functions, fourier analysis, and uncertainty — v: the discrete case, BSTJ57, 1371 (1978)

    D. Slepian, Prolate spheroidal wave functions, fourier analysis, and uncertainty — v: the discrete case, BSTJ57, 1371 (1978)

  72. [72]

    Kaiser and R

    J. Kaiser and R. Schafer, On the use of the i0-sinh window for spectrum analysis, IEEE Trans. Acoust. Speech Signal Process.28, 105 (1980)

  73. [73]

    D. B. Percival and A. T. Walden,Spectral Analysis for Physical ApplicationSs(Cambridge University Press, 1993). 38

  74. [74]

    Y. Dong, L. Lin, and Y. Tong, Ground-state preparation and energy estimation on early fault-tolerant quantum computers via quantum eigenvalue transformation of unitary matrices, PRX Quantum3, 040305 (2022)

  75. [75]

    Kosugi, Y

    T. Kosugi, Y. Nishiya, H. Nishi, and Y.-i. Matsushita, Imaginary-time evolution using forward and backward real-time evolution with a single ancilla: First-quantized eigensolver algorithm for quantum chemistry, Phys. Rev. Res.4, 033121 (2022)

  76. [76]

    Vidal, Efficient simulation of one-dimensional quantum many-body systems, Phys

    G. Vidal, Efficient simulation of one-dimensional quantum many-body systems, Phys. Rev. Lett.93, 040502 (2004)

  77. [77]

    Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics303, 2 (2003)

    A. Kitaev, Fault-tolerant quantum computation by anyons, Annals of Physics303, 2 (2003)

  78. [78]

    K. J. Satzinger, Y. J. Liu, A. Smith, C. Knapp, M. Newman, C. Jones, Z. Chen, C. Quintana, X. Mi, A. Dunsworth, C. Gidney, I. Aleiner, F. Arute, K. Arya, J. Atalaya, R. Babbush, J. C. Bardin, R. Barends, J. Basso, A. Bengtsson, A. Bilmes, M. Broughton, B. B. Buckley, D. A. Buell, B. Burkett, N. Bushnell, B. Chiaro, R. Collins, W. Courtney, S. Demura, A. R...

  79. [79]

    T. G. Kolda and B. W. Bader, Tensor decompositions and applications, SIAM Review51, 455 (2009)

  80. [80]

    explanatory

    R. A. Harshmanet al., Foundations of the parafac procedure: Models and conditions for an “explanatory” multi-modal factor analysis, UCLA working papers in phonetics16, 84 (1970)

Showing first 80 references.