Back to Rankings

Spectrum Estimation is Almost as Hard as Tomography

Marco Fanizza, Ryan O'Donnell, Chirag Wadhwa

Jul 31, 2026arXiv:2607.29680v1
quant-phcs.DS
Share
Scorecard· 16/16
8.0/10 impact

Resolves a long-standing, explicitly-flagged open problem (superlinear spectrum-estimation lower bound) with a novel reusable technique yielding three near-optimal bounds.

Abstract

We study the sample complexity of estimating and testing fundamental unitarily invariant properties of unknown quantum states; namely, the tasks of spectrum estimation, von Neumann entropy estimation, and rank-testing. For dd-dimensional states, and for every γ>0γ>0, we prove a sample complexity lower bound of Ω(d2γ)Ω(d^{2-γ}) for spectrum estimation to constant sorted total-variation error, entropy estimation to constant additive error, and rank-testing to constant trace distance. Our hard instances are constructed from sandwiched products of Haar-random projectors, suitably normalized using a novel technique that lets us derive explicit expressions for high-order tensor moments of the resultant states. These moments can be expressed as symmetric functions of Jucys--Murphy elements of the symmetric group algebra. To show that two such mixtures are indistinguishable, we analyze the log-likelihood ratio and perform moment-matching, i.e., we set its low-order Jucys--Murphy components to zero. Indistinguishability is then obtained by bounding an ff-divergence through the high-order components; the non-zero high-order terms and concentration of functions of Haar-random unitaries also imply separations in typical spectra, entropies, and ranks, proving all our lower bounds.

AI Impact Assessments

(1 models)

Scientific Impact Assessment

Core Contribution

This paper resolves a well-known open problem in quantum learning theory: it proves that estimating the spectrum (eigenvalues) of an unknown dd-dimensional quantum state to constant precision requires Ω(d2γ)\Omega(d^{2-\gamma}) copies for every γ>0\gamma>0 — nearly as many as full state tomography (Θ(d2)\Theta(d^2)). The natural intuition, which underpinned a two-stage paradigm in quantum tomography algorithms, was that estimating dd eigenvalues should be dramatically cheaper than reconstructing the Θ(d2)\Theta(d^2)-parameter density matrix. This work overturns that intuition. The same techniques yield near-quadratic lower bounds for von Neumann entropy estimation (improving from Ω(d/logd)\Omega(d/\log d) to d2o(1)d^{2-o(1)}) and rank-testing (improving from Ω(r)\Omega(r) to d2γd^{2-\gamma}), showing one-sided and two-sided rank testing are almost equally hard.

The methodological core is genuinely novel: (1) a "tilting" technique that decouples the intractable normalization E[Xn/Tr(X)n]\mathbb{E}[X^{\otimes n}/\mathrm{Tr}(X)^n] into E[Xn]/E[Tr(X)n]\mathbb{E}[X^{\otimes n}]/\mathbb{E}[\mathrm{Tr}(X)^n], enabling use of clean closed-form tensor moments; (2) expressing these moments as symmetric functions of Jucys–Murphy elements of the symmetric group algebra, so distinguishability reduces to tractable combinatorics rather than delicate Schur-polynomial cancellations; and (3) reducing moment-matching to the classical Prouhet–Tarry–Escott problem, solved via Thue–Morse sequences. This is a quantum analogue of the classical moment-matching/Poissonization machinery that drove distribution-testing lower bounds.

Methodological Rigor

The paper is a pure theory contribution and the proofs appear careful and complete. The argument structure is sound: they establish statistical indistinguishability (Theorem 5.7) via an ff-divergence (triangular discrimination) bound, using a clean lemma (6.2) splitting the divergence into a "bad set" contribution and a log-likelihood deviation. They then independently establish separation of spectra, entropies, and ranks using free-probability limits and Meckes' concentration for Lipschitz functions of Haar-random unitaries. The combinatorial lemmas (6.5–6.11) on Cayley lengths and orbit counts are self-contained and appear rigorous. A notable strength is the "warmup" Ω(d4/3)\Omega(d^{4/3}) bound (rank-d/2d/2 projectors vs. Haar marginals), which cleanly illustrates the method before the general construction. The result is near-optimal: it matches the O(d2(loglogd/logd)2)O(d^2 (\log\log d/\log d)^2) upper bound up to do(1)d^{o(1)}.

Potential Impact

The impact within quantum learning theory is high. Spectrum estimation is a fundamental primitive; entropy estimation underlies quantum source coding, entanglement quantification, and many-body physics diagnostics. Establishing that these are essentially as hard as tomography reshapes algorithm design — it tells practitioners that the two-stage "eigenvalues first, cheaply" strategy cannot save asymptotically, and it certifies the classical EYD algorithm as near-optimal. More broadly, the Jucys–Murphy + tilting toolkit is likely to be reused: the authors explicitly frame it as a general recipe for mixture-vs-mixture quantum lower bounds, a class of problems for which no superlinear bounds previously existed. This continues a trend of re-deriving representation-theoretic quantum-learning results with lighter machinery. Direct real-world/industrial application is limited — this is foundational complexity theory — but it informs resource estimation for quantum devices.

Timeliness & Relevance

Highly timely. The superlinear lower bound was explicitly flagged as open by Wright's thesis and recent work (PSTW26, PTTW26), with numerical evidence already pointing to d2o(1)d^{2-o(1)}. Multiple 2026 papers cited are contemporaneous, indicating an active competitive subfield. This paper delivers the sought-after proof.

Strengths & Limitations

Strengths: Resolves a clearly-stated, long-standing open problem; introduces a reusable and elegant technical framework; provides three distinct lower bounds from one construction; near-optimal matching to known upper bounds; self-contained proofs that avoid heavy representation theory.

Limitations: (1) The bounds are for the constant-precision regime against fully-entangled measurements; the tight ϵ\epsilon-dependence (conjectured Θ(d2/ϵ2log2d)\Theta(d^2/\epsilon^2 \log^2 d)) and the practically-important unentangled-measurement setting remain open (the latter still has a large d3/2d^{3/2}d3d^3 gap). (2) The γ>0\gamma>0 slack means it is d2o(1)d^{2-o(1)}, not exactly d2d^2. (3) The AI-use disclosure (ChatGPT contributing to the general-degree extension, and a claimed candidate proof of the ϵ\epsilon-dependent conjecture "to be reviewed") is unusual and slightly unsettling from a verification standpoint, though the authors take full responsibility and the presented proofs stand on their own.

Additional observations: The construction connects several deep threads — free multiplicative convolution (Belinschi's atom formula), Selberg integrals, Marchenko–Pastur/Page-curve results, and the number-theoretic Prouhet–Tarry–Escott problem — demonstrating both breadth and craft. The paper is well-organized, with a strong technical overview and a clear separation of the indistinguishability and separation arguments. Reproducibility, in the theoretical sense, is high: the constructions and proofs are fully specified.

Overall, this is a high-impact theoretical contribution that closes a central open question and delivers a transferable method, likely to be cited and built upon across the quantum property-testing subfield.

```json

{

"score": 8.0,

"score_reason": "Resolves a long-standing, explicitly-flagged open problem (superlinear spectrum-estimation lower bound) with a novel reusable technique yielding three near-optimal bounds.",

"significance": 8.0,

"significance_reason": "Overturns the intuition that eigenvalue estimation is far cheaper than tomography, certifies EYD as near-optimal, and provides the first superlinear quantum mixture-vs-mixture lower bound framework.",

"rigor": 8.5,

"rigor_reason": "Complete, careful proofs with a clean warmup, self-contained combinatorial lemmas, and independent indistinguishability and separation arguments matching known upper bounds.",

"novelty": 8.5,

"novelty_reason": "The tilting decoupling of normalization plus Jucys–Murphy moment-matching via the Prouhet–Tarry–Escott problem is a genuinely new and non-obvious toolkit for quantum lower bounds.",

"clarity": 8.0,

"clarity_reason": "Strong technical overview, logical organization separating indistinguishability from separation, and a pedagogical warmup, though the dense combinatorics require careful reading.",

"difficulty": 9.0,

"difficulty_reason": "Requires deep specialist expertise across representation theory, free probability, random matrix concentration, and number-theoretic moment-matching.",

"surprisingness": 7.0,

"surprisingness_reason": "The result contradicts the natural expectation that d eigenvalues are cheaper to learn than d^2 parameters, though numerical evidence had already hinted at it.",

"reproducibility": 9.0,

"reproducibility_reason": "As a theory paper the constructions, lemmas, and proofs are fully specified and self-contained, allowing independent verification.",

"translational_potential": 2.5,

"translational_potential_reason": "Foundational complexity result informing resource estimation for quantum learning, but with no direct near-term industrial or commercial deployment.",

"evidence_strength": 8.5,

"evidence_strength_reason": "Every claim is supported by explicit proofs matching known upper bounds to within d^{o(1)}, with three distinct corollaries from one construction.",

"generalisability": 7.0,

"generalisability_reason": "The method extends across spectrum, entropy, and rank testing and is framed as a general recipe, though results are confined to constant-precision, fully-entangled-measurement settings.",

"interdisciplinarity": 4.0,

"interdisciplinarity_reason": "Primarily serves quantum information theory but draws on and connects free probability, random matrix theory, symmetric-group representation theory, and classical distribution testing.",

"refutation_value": 6.0,

"refutation_value_reason": "It refutes the load-bearing assumption underlying two-stage tomography algorithms that spectrum estimation is asymptotically far easier than full tomography.",

"replication_value": 3.0,

"replication_value_reason": "It converts prior numerical evidence (PTTW26) of a superlinear bound into a proof, confirming a conjectured but unproven claim through a different (analytic) methodology.",

"resource_intensity": 2.0,

"resource_intensity_reason": "Pure theory requiring only expert human effort and no compute or experimental infrastructure to produce or extend.",

"foundationality": 7.0,

"foundationality_reason": "The tilting plus Jucys–Murphy moment-matching framework is explicitly positioned as a reusable primitive for future quantum mixture-vs-mixture lower bounds."

}

```

Rating:8/ 10
Significance 8Rigor 8.5Novelty 8.5Clarity 8

Generated Aug 3, 2026

Comparison History (0)

No comparisons yet.