Abstract 1 Introduction 2 Preliminaries 3 Results 4 Physical relevance 5 Open problems References

On the Complexity of Unique Quantum Witnesses
and Quantum Approximate Counting

Anurag Anshu ORCID Harvard University, Cambridge, MA, USA Jonas Haferkamp ORCID Saarland University, Saarbrücken, Germany Yeongwoo Hwang ORCID Harvard University, Cambridge, MA, USA Quynh T. Nguyen ORCID Harvard University, Cambridge, MA, USA
Abstract

We study the long-standing open question on the power of unique witnesses in quantum protocols, which asks if 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠, a variant of 𝖰𝖬𝖠 whose accepting witness space is 1-dimensional, contains 𝖰𝖬𝖠 under quantum reductions.

This work rules out any black-box reduction from 𝖰𝖬𝖠 to 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 by showing a quantum oracle separation between 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠. This provides a contrast to the classical case, where the Valiant-Vazirani theorem shows a black-box randomized reduction from 𝖴𝗇𝗂𝗊𝗎𝖾𝖭𝖯 to 𝖭𝖯, and suggests the need for studying the structure of the ground space of local Hamiltonians in distilling a potential unique witness. Via similar techniques, we show, relative to a quantum oracle, that 𝖰𝖬𝖠𝖰𝖬𝖠 cannot decide quantum approximate counting, ruling out a quantum analogue of Stockmeyer’s algorithm in the black-box setting. Our results employ a subspace reflection oracle, previously considered in [4, 3, 51], but we introduce new tools which allow us to exploit the unique witness constraint. We also show a strong “polarization” behavior of 𝖰𝖬𝖠 circuits, which could be of independent interest in studying quantum polynomial hierarchies.

We then ask a natural question; what structural properties of the local Hamiltonian problem can we exploit? We introduce a physically motivated candidate by showing that the ground energy of local Hamiltonians that satisfy a computational variant of the eigenstate thermalization hypothesis (ETH) can be estimated through a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 protocol. Our protocol can be viewed as a quantum expander test in a low energy subspace of the Hamiltonian and verifies a unique entangled state across two copies of the subspace. This allows us to conclude that if 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 is not equivalent to 𝖰𝖬𝖠, then 𝖰𝖬𝖠-hard Hamiltonians must violate ETH under adversarial perturbations (more accurately, further assuming the quantum PCP conjecture if ETH only applies to extensive energy subspaces). Under the same assumption, this also serves as evidence that chaotic local Hamiltonians, such as the SYK model may be computationally simpler than general local Hamiltonians.

Keywords and phrases:
Quantum complexity, approximate counting, Valiant-Vazirani, eigenstate thermalization hypothesis
Funding:
Anurag Anshu: AA acknowledges support through the NSF award QCIS-FF: Quantum Computing & Information Science Faculty Fellow at Harvard University (NSF 2013303).
Jonas Haferkamp: JH acknowledges support through the Harvard Quantum Initiative postdoctoral fellowship.
Yeongwoo Hwang: YH is supported by the National Science Foundation Graduate Research Fellowship under Grant No. 2140743. Any opinion, findings, and conclusions or recommendations expressed in this material are those of the authors(s) and do not necessarily reflect the views of the National Science Foundation.
Quynh T. Nguyen: QTN acknowledges support through the Harvard Quantum Initiative PhD fellowship.
Copyright and License:
[Uncaptioned image] © Anurag Anshu, Jonas Haferkamp, Yeongwoo Hwang, and Quynh T. Nguyen; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Quantum computation theory
Related Version:
Full Version: https://arxiv.org/abs/2410.23811
Acknowledgements:
In an earlier version of this paper, the oracle separation between 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 was given via a quantum oracle. We thank Chinmay Nirkhe and Anand Natarajan for pointing out that our proof did not in fact need quantumness in this case. This version has been updated to provide a classical oracle separation. We thank Soonwon Choi, Sam Garratt, Yingfei Gu, Rahul Jain and Shivaji Sondhi for insightful discussions, and especially thank Soonwon Choi for sharing an example which shows that the computational ETH assumption is incorrect if the Gaussian prescription is applied at high energies (we only applying it at low energies).
Funding:
This work was done in part while the authors were visiting the Simons Institute for the Theory of Computing, supported by DOE QSA grant number FP00010905. AA and QTN acknowledge support through the NSF Award No. 2238836.
Editor:
Shubhangi Saraf

1 Introduction

One of the central aims of quantum complexity theory is to explore the structure of ground, and more generally low-energy, subspaces in quantum many-body systems. A key question in this field is whether these low-energy subspaces of local Hamiltonians are computationally simpler than generic quantum subspaces.111From an information-theoretic point of view, the ground space of a local Hamiltonian is known to be simpler as it has a polynomial-sized tensor network approximation [50]. However this tensor network is not known to be computationally tractable. A long-standing problem which formalizes this inquiry is the relationship between 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 [9]. The complexity class 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 captures verification tasks in which positive instances have a 1-dimensional accepting subspace (unique quantum witness) while 𝖰𝖬𝖠 places no restriction on the dimension of the accepting subspace. The open question is to determine the relationship between these classes: are they equivalent under quantum reductions, or does the uniqueness promise strictly reduce computational power? On one hand, an equivalence would imply that a polynomial-time quantum algorithm can uniquely verify some fixed state in a low-energy subspace of local Hamiltonians, a task that intuitively seems unattainable for generic quantum subspaces. On the other, a separation would formalize the physics intuition that “natural” Hamiltonians with a non-degenerate ground space and inverse polynomial spectral gap are computationally simpler than 𝖰𝖬𝖠-hard Hamiltonians and thus amenable to analysis.

This question is also motivated by the insights produced by studying the class 𝖴𝗇𝗂𝗊𝗎𝖾𝖭𝖯. Classically, the celebrated Valiant-Vazirani isolation lemma [59] was used to show that 𝖴𝗇𝗂𝗊𝗎𝖾𝖭𝖯 contains 𝖭𝖯 under randomized reductions, i.e., 𝖭𝖯𝖱𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖭𝖯 (where 𝖱𝖯 is a one-sided error version of 𝖡𝖯𝖯), suggesting that the uniqueness promise does not significantly reduce the computational power of 𝖭𝖯, at least under randomized reductions. The isolation lemma has underlaid various developments in theoretical computer science, including approximate counting [57], Toda’s theorem [58], parallel computation [43], one-way functions [26], and average-case complexity [15]. One could hope that studying the possibility of a “quantum isolation lemma” would sharpen our understanding in a similarly wide range of quantum problems.

Limitations of unique quantum witnesses.

There has been a series of works attempting to resolve this question. Aharonov et al. [9] (first appearing on arxiv in 2008) initiated the study of 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and possible quantum analogues of the isolation lemma. There, the authors extended the hashing-based arguments in the Valiant-Vazirani isolation lemma to the classes 𝖬𝖠 and 𝖰𝖢𝖬𝖠. However, they noted the failure of a similar strategy for 𝖰𝖬𝖠, citing the fact that two random quantum states have exponentially small overlap. A later work [32] sought a completely different approach to obtain a deterministic isolation procedure in a special case, but also fell short of achieving a unique witness protocol for general 𝖰𝖬𝖠 problems as their strategy worked only for accepting subspaces of polynomial dimension. This led [32] to conjecture that a quantum reduction is required to obtain a unique verifier, suggesting the question of whether the inclusion 𝖰𝖬𝖠𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 holds.222In this work, is implicit that the oracularized complexity classes refers to the promise version (i.e. 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 actually refers to 𝖯𝗋𝗈𝗆𝗂𝗌𝖾𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠). This discussion brings us to our first result.

Theorem 1.

There exists a quantum oracle relative to which 𝖰𝖬𝖠𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠.

This oracle separation implies there is no black-box quantum isolation procedure for quantum witnesses, unlike its classical counterpart [59], and explains the failures of all approaches in [9, 32]. It may also be viewed as some evidence for the inequality between 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠. Furthermore, it suggests the need to exploit the physical structure of quantum low-energy spaces to construct a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 verifier for local Hamiltonians (see our result below).

Our Theorem 1 is based on a quantum oracle reflecting about an unknown subspace, and the task is to determine whether the subspace is non-empty (of some large dimensionality). Although this result in a sense is a continuation of a series of works on quantum oracle separations and unitary property testing [4, 31, 51, 6], the techniques used in our work differ in that they make careful use of the structure of the quantum witness. Previous separations used techniques which were not sensitive to the presence of a witness (e.g. replacing the witness with the maximally mixed state), which would undermine the uniqueness promise in our case. We hope that the techniques presented in this work will help develop similar query lower bounds on other witness-restricted quantum classes, such as 𝖰𝖬𝖠(2) where constructing oracle separations have been challenging.

Along the way to Theorem 1 we also obtain a classical oracle separation between 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 (Corollary 12). In fact, we can even separate 𝖭𝖯 from 𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 under the same classical oracle. However, this separation is perhaps not as interesting as Theorem 1, for one would a priori expect randomized or quantum reductions were needed to reduce 𝖰𝖬𝖠 to 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠, in analogy to the classical case. Regardless, we note that derandomizing the Valiant-Vazirani lemma is a long-standing open question by itself [13, 37], and this intermediate result rules out such a derandomization in the black-box setting. Thus, we believe the relationship between 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 is the proper question through which to study the power of unique quantum witnesses.

On classical vs. quantum oracles.

Both of our separations Theorem 1 and Theorem 2 are by quantum oracles. In our approach in Theorem 1, the quantumness is necessary as otherwise the classical version of our oracle problem is easily shown to be contained in 𝖭𝖯𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠. A similar observation can be made for Theorem 2. Historically, when the class or resource is classical, a classical separation can be viewed as the “gold-standard,” as evidenced by the long series of works aiming to separate 𝖰𝖢𝖬𝖠 and 𝖰𝖬𝖠 with a classical oracle (see the exposition in [60]). However, this decision should be made on a case-by-case basis. In our context, the subspace reflection oracle captures a way in which an algorithm might probe the ground space of a Hamiltonian. Naively, the best algorithmic way to “access” this subspace is to use quantum phase estimation to, say, implement reflections about it. There is a priori no known reason to suspect that this subspace will be well structured for worst-case Hamiltonians and thus we instantiate the ground space with a subspace in a randomly chosen basis. Another consideration when using and interpreting quantum oracle separations is that they only rule out quantum relativizing proof strategies. Indeed, [1, 7] already point out that the proof technique of explicitly representing quantum amplitudes is classically relativizing but quantumly not (and, to our knowledge, this is the only known proof technique with such property). However, this technique has only been used used to show containments in very powerful classical classes like 𝖤𝖷𝖯 and 𝖯𝖲𝖯𝖠𝖢𝖤. Thus, quantum oracle separations are still useful in ruling out black-box reductions at lower complexity regimes like 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠.

Improved quantum approximate counting lowerbounds.

Classically, the question of uniqueness is closely related to approximate counting. In [59], a random hash is used to isolate a single solution to an 𝖭𝖯 circuit. The celebrated Stockmeyer algorithm [57] uses similar tools to give a 𝖡𝖯𝖯𝖭𝖯 algorithm for estimating the number of accepting witnesses to multiplicative precision. A quantum analogue of this is to multiplicatively estimate the dimensionality of the accepting subspace of a 𝖰𝖬𝖠 circuit [18, 53, 16]. The complexity class capturing this quantum approximate counting problem is 𝖰𝖷𝖢. As usual, one should be careful when defining quantum analogues of classical classes. However, this problem is well motivated and shown in [17] to be connected to a number of physical relevant computational problems, such as computing quantum partition functions and estimating local observables of Gibbs states. 𝖰𝖷𝖢 was further studied in [16] and shown to contain the problem of estimating the Kronecker coefficients of the symmetric group. Stockmeyer’s approximate counting algorithm has led to many important applications, including the sampling-counting equivalence [33] and even the theoretical foundations for quantum advantage experiments [28]. Thus, it is natural to ask whether there exists a quantum analogue: Could quantum approximate counting be solved in (quantum) polynomial time with access to an oracle solving 𝖰𝖬𝖠 decision problems? In essence, do we have 𝖰𝖷𝖢𝖡𝖰𝖯𝖰𝖬𝖠?

The phrase “quantum approximate counting” has appeared elsewhere in the literature. In [39, 2], the authors study SBQP, the quantum analogue of SBP (where SBP is the multiplicatively precise analogue of BPP), which captures the task of estimating the acceptance probability of a 𝖡𝖰𝖯 machine to multiplicative precision. This is not known to be equivalent to the definition of QXC (unlike SBP which is equivalent to classical approximate counting). We do not consider 𝖲𝖡𝖰𝖯 in this work. Additionally, in [3], “quantum approximate counting” is used to refer to the task of designing quantum (𝖡𝖰𝖯 or 𝖰𝖬𝖠) algorithms for classical approximate counting, given quantum queries to a Boolean oracle. [51] consider a quantum generalization of the problem in [3], called the quantum approximate dimension problem; the task here is to decide whether an unknown subspace has dimension w or 2w, given access to the corresponding reflection oracle. Using the quantum approximate dimension problem, we show

Theorem 2 (Corollary 22, restated).

There exists a quantum oracle relative to which 𝖰𝖷𝖢𝖡𝖰𝖯𝖰𝖬𝖠, and in fact, 𝖰𝖷𝖢𝖰𝖬𝖠𝖰𝖬𝖠.

Our Theorem 2 rules out the existence of a quantum analogue of Stockmeyer’s theorem in the black-box setting. This is a strengthening of a result of [51], where the authors show a 𝖰𝖬𝖠 query lower bound against this task via a reduction to the classical problem in [3]. Again, the quantum oracle is necessary in our approach, because the classical version of our oracle is the classical approximate counting problem, which is contained in 𝖡𝖯𝖯𝖭𝖯𝖡𝖰𝖯𝖰𝖬𝖠. The polynomial methods used in [3, 51] do not transfer to the relativized setting in the presence of 𝖰𝖬𝖠-oracle calls. Moreover, prior works [2, 6] dealing with quantum versions of the polynomial hierarchy are all based on variants of the switching lemma, which are not known to be compatible with quantum oracles. Therefore, we need to develop new techniques to handle quantum oracle calls in “stacked” complexity classes like 𝖰𝖬𝖠𝖰𝖬𝖠.

Along the way we also prove a classical oracle separation 𝖰𝖷𝖢𝖰𝖬𝖠, in fact, 𝖲𝖡𝖯𝖯𝖰𝖬𝖠, thus giving a much simpler and elementary proof of the 𝖰𝖬𝖠 lower-bound against classical approximate counting of [3] (who showed 𝖲𝖡𝖯𝖰𝖬𝖠). Our separation is perhaps stronger as 𝖯𝖰𝖬𝖠 likely strictly contains 𝖰𝖬𝖠 as it contains 𝖼𝗈𝖰𝖬𝖠.

Unique witness from a well-connected ground space.

Theorem 1 underscores the importance of gaining a deeper understanding of the structure of quantum ground spaces to extract a unique quantum witness. For this, it is crucial to identify classes of Hamiltonians for which we know how to verify a unique low-energy witness. However, the only such classes of Hamiltonian known so far are classical Hamiltonians [59], local Hamiltonians with ground states that can be prepared by a polynomial-sized quantum circuit [9], and local Hamiltonians with a polynomial sized low-energy subspace [32]. Our final result identifies a new class of Hamiltonians for which we can verify a unique low energy state. Our starting point is a simple observation already hinted at in the previous paragraphs: An algorithm that only uses subroutines treating the entire Hamiltonian H as a single operator cannot distinguish between states of the same energy as they are indistinguishable with respect to H. For example, quantum phase estimation, which only uses time evolutions of H, falls into this category. Thus, any quantum algorithm that aims to verify a unique low-energy state must use operators which are more fine-grained than H itself, namely the individual local terms in H.

The action of local operators on Hamiltonian eigenstates is in itself a challenging problem, with few general results known [12]. Despite this, there are families of local Hamiltonians for which useful statements can be made. A particularly interesting such family is the set of local Hamiltonians that satisfy the Eigenstate Thermalization Hypothesis (ETH) from quantum statistical mechanics [56]. This hypothesis roughly says that within an energy window, eigenstates are well-connected in the sense that local operators induce transitions from one eigenstate to many others. We adopt the following informal definition from [19].

ETH (informal). There is a collection of polynomially many local observables A1,Am such that for every eigenstate |ψ with energy roughly E, the vector Ai|ψ roughly stays within a small energy window around E, and has a fairly good overlap with all the eigenstates in this window.

See Section 3.3.1 for the formal definition of ETH (Hypothesis 24) and a discussion on its physical origins. From the computer science perspective, a well-connected low-energy eigenspace should have features of an expanding graph and hence allow one to single out a unique “fixed point”. We solidify this intuition by generalizing the quantum expander test for verifying the maximally entangled state [10] and show the following theorem:

Theorem 3 (Informal restatement of Theorem 25).

There is a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 protocol for estimating the ground energy of a local Hamiltonian satisfying ETH.

The protocol verifies a unique quantum state in two copies of the subspace with energy around E, if ETH holds in that subspace. Assuming E is the ground energy itself, we can extract a unique ground state witness. If E is larger but sub-extensive (that is, nc with c<1, as a function of the number of qubits n), we can still obtain a unique witness for the purposes of 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 protocols due to known reductions concerning the local-Hamiltonian problem.

A likely implication of Theorem 3 is that the complexity of local Hamiltonians satisfying ETH, such as chaotic quantum Hamiltonians including the SYK model [55], is significantly lower than that of general local Hamiltonians. It is unclear if the theorem would help establish that 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 are equivalent under quantum reductions, since we do not know if 𝖰𝖬𝖠-hard local Hamiltonians would satisfy the ETH assumption.

2 Preliminaries

2.1 Complexity classes

Definition 4 (𝖰𝖬𝖠).

A promise problem L=(Lyes,Lno) is in the complexity class Quantum Merlin-Arthur (or 𝖰𝖬𝖠) if there exists a polynomial-time classical algorithm yielding a description of a quantum circuit V which implements a quantum operator Ux: with =RinRaux, such that

  1. 1.

    (Completeness) For all xLyes, there exists a witness |ϕRin such that

    ΠaccUx|ϕRin|0Raux223
  2. 2.

    (Soundness) For all xLno and |ϕRin,

    ΠaccUx|ϕRin|0Raux213

where Πacc can be taken to be the projector into the |0 on the first qubit.

Definition 5 (Unique 𝖰𝖬𝖠).

A promise problem L=(Lyes,Lno) is in the complexity class Unique Quantum Merlin-Arthur (or 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠) if in addition to the above properties of 𝖰𝖬𝖠, we have the further restriction that in the completeness case, there is some accepting state |ϕ and for all states orthogonal to |ϕ, the verifier accepts with probability equal to the soundness case. In particular,

  1. 1.

    (Completeness) For all xLyes, there exists a unique witness |ϕwitness such that

    ΠaccUx|ϕRin|0Raux223

    and for all states |ψ|ϕ,

    ΠaccUx|ψRin|0Raux213

2.2 Quantum oracles

Our first result is a oracle separation between 𝖰𝖬𝖠 and 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠. There have been various types of oracles considered in previous works. Although the “gold standard” for oracle separations is arguably a classical oracle (which can be queried in superposition), quantum separations have required more complex oracles, such as the distributional oracles of [45]. In this work, we work with quantum oracles, as studied in [4]. Thus, in our setting, an oracle 𝒪 refers to some arbitrary unitary operation. We will also allow controlled queries to 𝒪, so that a call to 𝒪 is in fact a call to the unitary,

(𝕀Π)𝕀+Π𝒪,

where Π is a projector representing a control subspace. Given a family of oracles, we can define a “oracle promise problem,”

Definition 6 (Oracle promise problem).

A oracle promise problem is defined by two disjoint sets of oracles (unitary operators), ΩYES and ΩNO. We use Ω=(ΩYES,ΩNO) to refer to the decision problem itself.

Then, any L-query algorithm deciding Ω should be able to receive an oracle 𝒪 either in the YES or NO case and use L calls to 𝒪 to decide which is the case. Formally, a oracular verifier for 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠𝒪 consists of alternating unitaries (consisting of local quantum gates) and calls to 𝒪. The task is to (uniquely) determine whether the unknown oracle corresponds to a YES case or NO case instance.

Definition 7 (𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 Oracle verifier).

Given a n-qubit oracle 𝒪 and mpoly(n), a m-qubit, L-query 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 oracular verifier V𝒪 is defined over a n-qubit oracle register R𝒪, in a naux qubit auxiliary register Raux. Additionally, there is a nin-qubit input register Rin. Thus, m=naux+nin. Given a nin-qubit witness, the verifier then performs a sequence of unitaries on the full m-qubit subspace, interwoven with L queries to an oracle 𝒪 on R𝒪, controlled on the remaining system. Finally, the verifier performs the measurement {|00|,|11|} on the first qubit in the auxiliary register (which we label as the register Rout), rejecting if it obtains |00| the and accepting if obtains |11|. See Figure 1. Thus, the measurement corresponding to the “accept” case can be written as

V (𝕀|11|Rout)UL𝒪UL1U1𝒪U0(𝕀|0naux0naux|Raux),
Πaccept VV. (1)

We say that V uniquely solves the oracle promise problem Ω=(ΩYES,ΩNO) if

  • If 𝒪ΩYES then there exists a witness |ψRin so that

    acc(V,ψ)=V|ψ|0Raux22=tr[Πaccept|ψψ|Rin|00|Raux]1ε,

    and for any witness |ψ where |ψ is orthogonal to |ψ, acc(V,ψ)ε.

  • If 𝒪ΩNO then acc(V)ε; for all witnesses the verifier accepts with probability at most ε.

Figure 1: Circuit diagram for the verifier V. R𝒪 is the oracle register, contained in the auxiliary register Raux, initialized to |0Raux. Rin is the input register. The final measurement M is performed on Rout, which is the first qubit of the auxiliary register.

The standard oracle that we consider in this work are variations of an oracle considered in prior works: the subspace reflection oracle 𝒪𝒮=𝕀2Π𝒮, where Π𝒮 corresponds to a projection onto a k-dimensional subspace 𝒮. In [4], 𝒪𝒮 with k=1 was used to separate 𝖰𝖢𝖬𝖠 from 𝖰𝖬𝖠; in [3], 𝒮 was restricted to be a classical subspace to separate 𝖰𝖬𝖠 from classical approximate counting; and finally [51] quantized the prior work to separate 𝖰𝖬𝖠 from quantum approximate counting. We borrow notation from these last two works to define the oracle promise problems we study.

Definition 8 (Quantum approximate dimension).

In the oracle problem QApxDim(k1,k2,N), the problem is to decide whether an n-qubit unknown oracle corresponds to a reflection about a subspace of dimension k2 (the YES case), or a subspace of dimension k1 (the NO case), promised that one of these two is the case. For technical reasons, we restrict k222n/3.

If we are further promised that the subspace is in the computational basis, then we have the classical problem ApxDim(k1,k2,N).

Variations of this oracle have been used in other works. The case when 𝒮 corresponds to a 1-dimensional subspace is exactly the oracle used by [4] to separate 𝖰𝖢𝖬𝖠 from 𝖰𝖬𝖠. When increasing the dimension and restricting the subspace 𝒮 to be classical, this becomes the classical oracle used to separate 𝖰𝖬𝖠 and classical approximate counting [3]; when the subspace is quantum, this is exactly the QApxDim problem and has been shown by [51] to separate 𝖰𝖬𝖠 and quantum approximate counting. Having said that, our results and proof techniques do not follow from these works in any straightforward way. For example, the [4] oracle would not work in our setting as this problem is easily seen to be in 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠. A better candidate is the oracle 𝒪S corresponding to a large subspace S, as there does not seem to be a natural “unique” witness that the prover could provide. However, when constructing our query lower bounds (which are, as usual, the heart of nearly all oracle separations), new techniques are necessary to account for the uniqueness restriction on the witness. Indeed, prior oracle separations involving 𝖰𝖬𝖠 avoid directly dealing with the witness either by conditioning on a particular classical witness [4], replacing the witness with the maximally mixed state [3, 51], or reducing to a classical circuit lower bound [2, 6]. Instead, our work carefully uses the requirement that the witness be unique.

 Remark 9 (Classical oracles).

We use the term “classical oracle” to refer to phase oracles 𝒪S=𝕀2ΠS, where S is restricted to be in the computational basis. The truly “classical” analogue would be the oracle mapping |x|b|x|b𝟏xS. However, as our separations holds under controlled oracles, these notions are equivalent.

In this work, we also consider oracle separations between relativized classes, e.g. 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠. We instantiate verifiers belong to a relativized class 𝒞1𝒞2 as a verifier for 𝒞1, with access to an black-box oracle 𝒪 for a 𝒞2-complete problem. In our case, 𝒞2 is either 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 or 𝖰𝖬𝖠 and thus 𝒜 is a classical oracle and the corresponding complete problem will be UniqueQCircuitSat or QCircuitSat respectively. However, one difficulty in defining it in this way is that these languages correspond promise problems, but there is not guarantee that an algorithm querying 𝒜 will do so only on valid instances. To deal with this issue, we take an approach used in prior works [2, 31] and imagine that the outer verifier is given access to an oracle 𝒜 for a decision “extension” of the promise problem.

Definition 10 (Promise oracle extension).

For any oracle given by the partial function 𝒪:{0,1}n{0,1,}, an oracle extension is a total function 𝒜:{0,1}n{0,1} such that 𝒪(x)=𝒜(x) for all xDom(𝒪).

Thus, we take the following definition of a verifier accessing a promise oracle 𝒪, as suggested in [2]:

V𝒪(x){1if V𝒜(x)=1 for every 𝒜 extending 𝒪,0if V𝒜(x)=0 for every 𝒜 extending 𝒪,otherwise.

As stated in [2], it is not obvious that this is the “right” way to define oracle access to a promise problem. Nonetheless, we believe this to be a fairly reasonable definition; when querying an oracle for, e.g., QuantumCircuitSat on an instance falling outside of the promise gap, why should the verifier expect a “reasonable” response from oracle?

Lastly, we remark that whenever we have a verifier of the form 𝒞1𝒞2, we assume that the oracle for 𝒞2 acts on instances of size polynomial in the outer problem instance. This is required technically as many of our arguments go through union bounding over possible inputs to the oracle for 𝒞2. Moreover, as noted in [3], there is a distinguisher between (classical) subspaces of dimension k and 2k with a witness of size roughly 𝒪(k). Thus, without a bound on the witness register, this distinguisher would preclude a 𝖯𝖰𝖬𝖠 query lower bound for approximate counting.

3 Results

3.1 Separating 𝗤𝗠𝗔 from 𝗕𝗤𝗣𝗨𝗻𝗶𝗾𝘂𝗲𝗤𝗠𝗔

We consider circuits with access to the unknown subspace reflection oracle 𝒪S, as well as an oracle 𝒜 for UniqueQCircuitSat (more precisely, their “oracularized” versions), the natural complete problems for 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠. Rather than directly working with ApxDim and QApxDim, we consider a slightly different oracle promise problem, QApxDimExact with parameters k1 and k2 so that the YES case corresponds to subspace dimension exactly k2, and the NO case corresponds to dimension exactly k1. Then, the specific problems we consider are,

EmptyNonEmpty ApxDimExact(0,1,N)ApxDimExact(0,k,N) and
QEmptyNonEmpty QApxDimExact(0,1,N)QApxDimExact(0,k,N).

where the union is taking within the NO case and YES case oracles individually (so that the NO case is still the identity oracle, and YES case corresponds to subspaces of dimension 1 and k). Since EmptyNonEmptyApxDim(0,1,N) (and QEmptyNonEmptyQApxDim(0,1,N)), we could’ve considered the harder problem instead, but we use EmptyNonEmpty since this includes all subspace dimensions needed to make the proof go through.

The overall idea for the separation is to use a hybrid argument and separately consider steps of the 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 verifier where a query is made to the circuit sat oracle 𝒜 and steps where the verifier directly queries the oracle 𝒪S. Therefore, our separation requires two ingredients.

  1. 1.

    First, we argue that the verifier’s calls to 𝒜 are either on circuits C which are themselves good verifiers for the oracle distinguishing task or are nearly useless. We show this by establishing a polarization property of oracular quantum circuits.

  2. 2.

    Second, we establish our oracle separation between 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠. This allows us to rule out the first case in the item above, as well as to argue that the verifier’s direct calls to 𝒪S are not too helpful.

Interestingly, for the second step, the separating oracle can actually be made classical and we show a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 query lower bound for 𝒪S restricted to a computation basis subspace.

Theorem 11 (𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 query lower bound, EmptyNonEmpty).

Let V be a L-query, m-qubit verifier which uniquely decides EmptyNonEmpty with error ε2n both on 𝒪No and on at least (12n)-fraction of YES instances. Then, LΩ(k). In particular, if we take k2Ω(n), then V requires exponentially-many queries in n.

Standard diagonalization arguments yield the following.

Corollary 12.

There exists a classical oracle 𝒪 such that 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠𝒪𝖰𝖬𝖠𝒪, and in fact, 𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠𝒪𝖭𝖯𝒪.

In a sense the proof of the above theorem is where the uniqueness property is crucially used. Imagine we have two classical subspaces S and T, where T is a subspace extending S, so that T=SΔ. We will pick the dimension of S and T to be on the order superpoly(n); we will show that as dim(S),dim(T)2n, a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 verifier using poly(n) oracle queries cannot distinguish S and random extensions T. When S is taken from as a YES instance of an oracle distinguishing task with a unique witness ϕS, this observation essentially says that the “same” witness works for T. Flipping the order of quantification, we can equivalently say that for a random T, nearly all of its “roots” S have the same witness. But since dim(T)dim(S), the S’s can be picked so that they are orthogonal. Finally, applying a variation of the hybrid method shows that the same witness must overlap too many orthogonal subspaces, yielding a contradiction.

However, when moving to Item 1, the use of quantum oracle is seemingly necessary. When considering the 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 verifier’s direct calls to 𝒪S, the correlation between the verifier’s state and the oracle is naturally captured by the overlap of some intermediate state and the projector onto the oracle subspace; we can get a lot of mileage by studying quantities of the form “tr[ΠSψ]” and averaging over subspaces S. The calls to the circuit sat oracles 𝒜 require more careful handling. The challenge is, for each oracular circuit C𝒪S on which the outer verifier is called, we would like to show that the answer 𝒜(C𝒪S) cannot help tell apart the yes and no case oracles (with high probability). We establish this via a concentration argument, enabled by our use of quantum oracles. In our case, we use a version of Levy’s lemma for the Grassmannian (the space of k-dimensional subspaces of N), due to [25].

Lemma 13 (Levy’s lemma on the Grassmannian (informal)).

Let Δ be uniform over GN,k and ΠΔ be corresponding projector. Given any quantum state |ψN, with NN, we have that

PrΔGN,k[tr[ΠΔ|ψψ|]2N1/3]exp(poly(N)).

Via Lemma 13 we show that the acceptance probabilities of oracular quantum polynomial time circuits must concentrate when averaged over the random quantum subspace oracles. A similar observation is made in [31] to rule out the possibility of a search-to-decision reduction for 𝖰𝖬𝖠. We use this observation somewhat differently. Whereas [31] used concentration to argue that on across yes instances the circuit behaves similarly, we instead want to argue the same thing across yes and no instances. This requires a different and more intricate argument considering the average acceptance probabilities in the yes and no cases individually. In the specific form used in this work, we call this property “polarization”.

Lemma 14 (Polarization of 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 oracle circuits).

Let V be an polynomial-time oracular verifier and k=Nα, with α2/3. Then for any function f(n)exp(N1/4):

Pr𝒯HaarN,k[𝒜𝒯(V)𝒜(V)]f(n).

The final obstacle when designing oracle separations for relativized (or “stacked”) classes is the fact that the natural complete problem for 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠, quantum circuit sat, is defined as a promise problem. As mentioned above, we address this issue by defining an appropriate “decision-extension” of the promise oracle.

Theorem 15 (Relativized Oracle Separation).

Let V𝒯 be any L-query 𝖡𝖰𝖯 verifier with gates implementing controlled calls to the subspace reflection oracle 𝒪S, as well as an oracle for UniqueQCircuitSat𝒪S. Suppose V𝒯 solves QEmptyNonEmpty, with error at most εexp(n). Then, if kNα for some αΩ(1), then LΩ(N/k).

Corollary 16.

There exists a quantum oracle 𝒪 such that 𝖡𝖰𝖯𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠𝒪𝖰𝖬𝖠𝒪.

3.2 Separating 𝗤𝗫𝗖 from 𝗤𝗠𝗔𝗤𝗠𝗔

In this section, we use the techniques from the previous result to obtain an oracle separation between 𝖰𝖬𝖠 and quantum approximate counting. This definition was introduced by [17].

Definition 17 (Quantum approximate counting, 𝖰𝖷𝖢).

Let V be a m-sized 𝖰𝖬𝖠 verification circuit which takes a n-qubit state |ψ as input, with mpoly(n). Then, given thresholds 0<b<a1 with ab1poly(n), and an error parameter ε1poly(n), the approximate counting problem 𝖰𝖷𝖢 is to compute an estimate D such that (1ε)DaD(1+ε)Db, where Da is defined as the dimension of the witness subspace for which V accepts with probability at least a, respectively for Db.

We formalize the above as oracle promise problem ApxDim(k1,k2,N)=(ΩYes,ΩNo) with k2N and k1k2(11/poly(n)) where,

  • Each 𝒪ΩYes corresponds to a reflection over a subspace SN of dimension k1, with 𝒪=𝕀2ΠS.

  • Each 𝒪ΩNo corresponds to a reflection over a subspace TN of dimension k2, with 𝒪=𝕀2ΠT.

Intuitively, the different parameters can be attributed to the fact that in 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠, the unique witness condition affords structure even within the YES case instances (thus, both subspaces of dimension 2k1n and 2k2n can be treated as YES cases). For 𝖰𝖬𝖠, we need to treat the dimension 2k2n as a NO case in order to carry through similar arguments. As in the case of 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠, we will combine an (classical) oracle separation between 𝖰𝖬𝖠 and 𝖰𝖷𝖢 as well a version of the polarization lemma for 𝖰𝖬𝖠 verifiers.

Theorem 18 (𝖰𝖬𝖠 query lower bound for ApxDim).

Let k1o(N) and k2=2k1. Let V𝒪 be an oracular 𝖰𝖬𝖠 verifier deciding ApxDim(k1,k2,N) with error εexp(n). Then, LΩ(N/k1).

Lemma 19 (Polarization of 𝖰𝖬𝖠 oracle circuits).

Let C be any n=poly(n)-qubit 𝖰𝖬𝖠 circuit which makes Lo(N/k1) queries to a n qubit oracle 𝒪𝒮=𝕀2Π𝒮. Let 𝒮 be a dimension k1 subspace. Let 𝒯 be a k2<N2/3 dimensional extension of 𝒮. Then for some decision-extension 𝒜 of an oracle for QCiruitSAT(η,1η),

Pr𝒯𝒬𝒮,k2[𝒜(C𝒮)𝒜(C𝒯)]2γ,

where γexp(N1/3).

 Remark 20.

A similar lower bound where dimension k1 corresponds to the YES case and dimension k2 corresponds to the NO case is obtained in [51, 3] but through much more intricate methods utilizing a version of the polynomial method for Laurent polynomials. Our proof is much simpler and gets the same quantitative bound, which is shown to be tight in [3, Theorem 5]. However, their techniques are better in that their lower bound is symmetric (it also works when the YES case is the larger subspace, and the NO case is smaller). Intuitively, this is because their proof goes via the inclusion of 𝖰𝖬𝖠 in SBQP, the exponentially precise analogue of 𝖡𝖰𝖯, and thus is not sensitive to a witness.

In both separations, the polarization property puts strong limitations on the way the outer verifiers can interact with the circuit sat oracles. As such, we hope these techniques will generalize beyond our setting to yield lower bounds against verifiers querying quantum oracles, such as those capturing constant stacks of 𝖰𝖬𝖠 (the “𝖰𝖬𝖠 hierarchy”).

Finally, we make a remark on our (classical) oracle separation of 𝖰𝖬𝖠 and 𝖰𝖷𝖢. One might be tempted to leverage an argument similar to [4] to remove the outer witness. This does not work, as the witness in our case is quantum and would require conditioning on a event with probability exp(2poly(n)). Instead, we use a hybridization argument similar to the 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 proof. When extending a (yes case) small subspace S to a (no cases) larger subspace SΔ, the only way the 𝖰𝖬𝖠 algorithm changes its behavior is that the witness overlaps with Δ. Thus, choosing many orthogonal extensions Δ leads to a contradiction. This recovers the 𝖰𝖬𝖠 lower bound for classical approximate counting (with quantum queries) of [3], albeit by a much more straightforward proof without involving the machinery of Laurent polynomials used in [3].

Theorem 21 (Query lower bound for 𝖰𝖬𝖠𝖰𝖬𝖠).

Let V be any m=poly(n)-qubit, L-query oracular 𝖰𝖬𝖠 verifier which purports to solve the n-qubit oracle promise problem QApxDim(k1,k2,N) with error εexp(n) and k1<k2N2/3. Suppose V has gates implementing controlled calls to a n qubit subspace reflection oracle 𝒪𝒮, as well as an non qubit oracle 𝒜 for QCircuitSat𝒪𝒮. Then LΩ(N1/3).

Corollary 22.

There exists a quantum oracle 𝒪 such that 𝖰𝖷𝖢𝒪𝖰𝖬𝖠𝖰𝖬𝖠𝒪.

3.3 Protocol for unique verification

The oracle separations in the previous section suggest that we need to use the structure of 𝖰𝖬𝖠 problems if we want to have a chance to design a unique verifier. In this section, we describe a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 verification protocol for local Hamiltonians that are sufficiently chaotic in the sense of the eigenstate thermalization hypothesis (ETH) – an influential hypothesis in quantum statistical mechanics.

3.3.1 Eigenstate thermalization hypothesis

Thermodynamics expects the evolution of a generic many-body system to be ergodic, meaning that the system evolves to eventually explore all possible configurations and thermalize. However, the evolution of a closed quantum system is governed by the Schrodinger’s equation which is reversible and also leaves eigenstates invariant. Why do individual energy eigenstates of a large, interacting quantum system look thermal when probed by “simple” observables (like local operators), even though the full wavefunction evolves unitarily without randomness? A widespread explanation to this seeming paradox is the eigenstate thermalization hypothesis (ETH) [23, 56], which has produced many predictions in agreements with thermalization on many physical systems [21]. ETH has many slight variants, but they are typically a set of assumptions about generic Hamiltonians that, informally, predicts that the eigenstates locally resemble a Gibbs state, and furthermore, local observables obey certain random matrix behaviors in the Hamiltonian energy eigenbasis inspired by quantum chaos and random matrix theory.

The first mathematically clean and consistent description of ETH was proposed by Srednicki [56], which is in turn inspired by Wigner’s pioneering random matrix theory (RMT) framework of analyzing quantum systems using statistical methods. Srednicki’s ETH proposes a “chaotic” behavior of local operators in the eigenbasis of a local Hamiltonian. It is hard to characterize precisely how a local operator A may act on an energy eigenstate |νi, so ETH leverages randomness to state that the matrix elements of a local operator A “looks like” νj|A|νi=ciδi,j+eS/2Rij, where ci is some quantity capturing νi|A|νi and |νi,|νj are eigenstates with nearby eigen-energies (a local operator cannot connect far away eigenstates - see [12]) and S is the entropy within this window. Crucially, Rij are chosen to be i.i.d random Gaussians. More formally, we consider the following formulation closely following [19, Section IV], which is arguably more mathematically concrete than [56] to work with at a rigorous level.

Let us first introduce some notations. Fix an energy e0>n0.99 and Δ=1/poly(n) throughout this section333ETH is usually defined for “extensive” energies, those that scale linearly with the system size. Our threshold n0.99 is just a way to capture this without introducing another parameter to set the extensive threshold.. Denote by Δ the energy window of width Δ/2 around e0, Δ=[e0Δ2,e0+Δ2]. Denote by ΠΔ the projector onto the subspace of eigenstates with energy eΔ. We will always think of the eigenstates as ordered by non-decreasing energies.

Hypothesis 23 (Srednicki’s ETH ansatzes).

Consider an n-qubit local Hamiltonian H and a narrow energy window of width Δ around an energy e0. The corresponding “ETH ansatz ensemble”, denoted 𝒢, is a distribution over Hermitian operators A𝒢 with parameters {μα}α and {fα,β}α,β, whose entries within the said window have the form

να|A|νβ=μαδα,β+fα,βTr(ΠΔ)gα,β, (2)

whenever the energies of the eigenstates |να, |νβ are in the interval Δ. In the off-diagonal entries, gα,β=gβ,α are complex i.i.d. Gaussians with mean 0 and variance 1. The real coefficients fα,β represent flexibility in the variance of the Gaussians, subject to fα,β=fβ,α due to the Hermitian condition and ffα,β1, where f is a constant. The diagonal entries μα=να|A|να are non-random, but smoothly vary within the window Δ (as a function of the eigenenergy). The normalization factor Tr(ΠΔ) is simply the number of eigenstates within this window, as ΠΔ denotes the projector onto them. All these prescriptions are subject to the (global) constraints A=A and A=1. 444This can be guaranteed by the countering fluctuations in the entries outside of the window Δ.

We say that H and a collection of local 555The locality of Ai is typically allowed up to o(n), but Ai should be simple (like Pauli, or efficiently implementable). For this work, it suffices to think of Ai as constant-local Pauli operators. observables A1,,Am satisfy ETH within the Δ window around energy e0 if these operators “look like” randomly drawn from the ETH ensemble.

A quick objection to the ETH ansatz is - where is the randomness coming from? And what does “look like” mean? Clearly the Hamiltonian is fixed and the statement is being made for all eigenstates in the energy window. The view is that we make a choice of gαβ by drawing from iid Gaussian and then fix them for the entire calculation (and any predictions are taken “with high probability”). Given that H is a complicated object, the belief is that we would be unable to tell a difference (with high probability over the ETH ensemble) if we did calculations using the precise entries of the A’s in the energy eigenbasis. Such behaviors have been experimentally and numerically verified in many physical quantum systems [48], as well as generic, random Hamiltonians, such as the famous SYK model, where ETH is postulated to hold, see [21, Section 4.3].

In this work, we take this computational indistinguishability viewpoint to design a unique verification protocol for ETH Hamiltonians, but we hope the formulation below could be useful for other applications of ETH in complexity theory.

Hypothesis 24 (Computational ETH).

Setup:

(The Hamiltonian)

Consider an n-qubit local Hamiltonian H and a fixed energy window Δ around e0. Suppose that H has non-clustering of eigenstates in the window Δ.

(The Experiment)

Consider a polynomial-time algorithm 𝒩 that takes in an input state |ξ promised to be inside Δ and runs a polynomial-sized circuit (possibly with ancillas) that ends with a binary POVM {𝖠𝖼𝖼,I𝖠𝖼𝖼}. Furthermore, there is a set of m=poly(n) local Pauli operators Alg={A1,,Am}, such that 𝒩 consists of gates that depend directly as a black-box on {A1,,Am} and the individual local terms in H (which for examples could be time evolutions of any of these operators, or measurements, etc.)

(The ETH Ensemble)

Let 𝒩 be an algorithm obtained by the following procedure: (1) iid sample a set of operators Alg={A1,,Am} from the operator ensemble defined in Hypothesis 23, (2) replace the components related to Alg in 𝒩 by their Alg versions 666These A may require exponential circuit complexity, but we think of them as a black box, and these replacements only appears in our calculations.. More precisely, within the window Δ, each Ai has the form

να|Ai|νβ=μαiδα,β+fα,βTr(ΠΔ)gα,βi, (3)

where the coefficients and random variables are as described in Hypothesis 23. Note that the diagonal entries μαi are non-random and thus depend on Ai, so that distinct Ai is sampled from a distinct distribution. The gα,βi are normal complex Gaussian variables sampled iid across different i.

Then, we say H and 𝒩 satisfy the computational ETH within the energy window Δ if it holds, with high probability over the randomness of the ETH ensemble, that

Tr((𝒩(ξ)𝒩(ξ))𝖠𝖼𝖼)<εETH=negl(n), states ξ of energy Δ. (4)

The above computational version is an attempt to capture in the computer science language the common physicists’ interpretations of ETH. Above, we allow the “experiment” to be beyond 𝖡𝖰𝖯: it can be given a proof (untrusted, except for the fact that its energy is within Δ which can be efficiently verified). This allows us to work in the context of 𝖰𝖬𝖠. The εETH is similar to the distinguishability advantage encountered in cryptography. Here we take it to be negl(n), but we actually only need εETH to be a sufficiently small constant to guarantee the performance of our 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 verifier in this section.

3.3.2 Verifying a unique state in an energy window satisfying ETH

Given a local Hamiltonian H on n qubits, in this section, we give a protocol that uniquely accepts a specific state in a known energy window around e0 satisfying certain properties that will be instantiated by ETH.

The intuition behind the protocol comes from tests of entanglement via quantum expanders [10]. More precisely, the EPR state |Ω=1Di=1D|ii can be tested locally using unitaries of the form UU¯. We briefly explain the idea: Clearly, UU¯|Ω=|Ω for any unitary USU(D) so any averaged operator over the actions UiU¯i for unitaries UiSU(D),i=1,,M satisfies 1Mi=1MUiU¯i|Ω=|Ω. A quantum expander is roughly a collection of unitaries Ui, such that the second highest eigenvalue of 1Mi=1MUiU¯i is bounded away from 1. Such quantum expander families exist with surprisingly small M and constant gap between second highest eigenvalue and 1 (see e.g. [27, 29]). Another common approach, used e.g. for self-testing, is to pick the local Pauli operators, which results in an inverse polynomial gap between the second highest eigenvalue and 1 (see e.g. [46]). The expectation value 1Mi=1Mψ|UiU¯i|ψ can be efficiently estimated in many settings and suggests a simple test: If the above expectation is close to 1, then |ψ is close to |Ω. Otherwise, |ψ and |Ω are far from each other. In particular, |Ω is uniquely accepted by this text. A plausible strategy is therefore to use this principle to uniquely verify the maximally entangled state on two copies of the space spanned by the eigenstates in an energy window. It will turn out that we cannot actually guarantee that the uniquely accepted state is the maximally entangled state on these subspaces. But for our purposes it will suffice that the Perron-Frobenius theorem guarantees the existence of a non-degenerate highest eigenvalue for an eigenstate with non-negative entries. Then, Merlin can simply send this state and Arthur can uniquely verify the ground state energy via a quantum expander.

Two obstacles prevent us from applying this naively: 1) A general unitary will move not respect the subspace spanned by eigenstates with energies in a window and 2) even if we can apply such unitaries we need to bound the gap of the averaging operator to guarantee uniqueness. We solve 1) by considering unitaries of the form eιεA for sufficiently small, but constant ε>0, which approximately preserve subspaces. We then use phase estimation to prepare measurements of the energy window and effectively apply the projector onto the respective subspaces. In order to approach 2) we invoke the ETH: By choosing the Ai to be local (think local Pauli matrices), their matrix entries in the energy window are indistinguishable from the Gaussians in Hypothesis 23. We then use probabilistic techniques to show with high probability over the Gaussians Ai that the operators 𝔼ieιεAieιεAi are sufficiently gapped. In particular, for the Gaussian Ai we obtain a gap and, consequently, unique verification. The ETH then guarantees that any quantum algorithm using the correct Ai will have the same behavior.

Finally, we assume there is a known energy e0 and interval Δ such that the above properties holds.

Assumption 1 (Existence of “good” Energy).

There exists a known energy e0(0,1) and value Δ such that Hypothesis 24 holds. Furthermore, e0 can be represented by O(logn) bits.

Given the existence of this energy window, we use quantum phase estimation to approximately implement the projector ΠΔ onto this subspace. In the algorithm, we use Π to indicate this projector, with subscripts indicating the registers on which the projector acts.

Algorithm 1 Energy subspace test Alg.
Input: state |ψ on the registers Rsys(1),Rsys(2).

We show that the above algorithm indeed succeeds in verifying the existence of a unique ground state.

Theorem 25 (Energy Subspace Test).

Fix a parameter e0 and consider a Hamiltonian H. There exists an algorithm Alg=Alg(H,e0,L) (where L=poly(n) is a precision parameter) such that

  1. 1.

    If λ0(H)>e0+3L then Alg accepts with probability at most 13.

  2. 2.

    If λ0(H)e0 and e0 satisfies Assumption 1 w.r.t. H, then there exists a unique state |Ψ such that Alg accepts on input |Ψ with probability at least 23, and on states orthogonal to |Ψ, Alg accepts with probability at most 13.

4 Physical relevance

On the surface, our first result can be viewed as evidence that 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 is not equal to 𝖰𝖬𝖠. An alternative interpretation is via a characterization due to [9]; in that paper, the authors show that the class of local Hamiltonians with an inverse polynomial spectral gap is 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 complete. Thus, our oracle separations give evidence that the local Hamiltonian problem is easier for inverse polynomially gapped instances777A more recent work [22] also studied the role of the spectral gap in small promise gap regimes.. This is interesting from a physics perspective as many classes of physically relevant Hamiltonians are known to have an inverse polynomial spectral gap (e.g., ferromagnetic Heisenberg [38], 1D AKLT [5], Movassagh-Shor [42]). In particular, the presence of a spectral gap in these models immediately makes studying various properties more tractable [47, 40]. Thus our work provides some evidence for the physical intuition that inverse polynomially gapped Hamiltonians are “easy” by separating such Hamiltonians from 𝖰𝖬𝖠, under an oracle.

Our third result has implications for the Eigenstate Thermalization Hypothesis. If 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 are not equivalent (under quantum reductions), then hard instances of the local Hamiltonian problem cannot satisfy ETH in its low energy window (near-extensive). This violation of ETH is robust against adversarial perturbations of near-extensive strength. Under the quantum PCP conjecture [8], which states that there are local Hamiltonians such that ground energy estimation is 𝖰𝖬𝖠-hard even up to extensively scaling accuracy, the latter statement can be extended to extensive energy scales. In other words, if 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 are not equivalent, then quantum PCP Hamiltonians need to robustly violate ETH. Note that the same assumption implies 𝖰𝖬𝖠-hard Hamiltonians are not many-body localized [44] either. Thus, our results suggest that these Hamiltonians may correspond to an interesting phase of matter that is robust against adversarial perturbations.

We expect our unique verification protocol to work for families of random quantum Hamiltonians such as random spin systems [24] or the SYK model [49, 34, 35, 30]. If random quantum Hamiltonians satisfy a sufficiently strong variant of ETH, our result implies that the average-case is considerably easier than the worst-case, under the assumption 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 and 𝖰𝖬𝖠 are not equivalent. Indeed, randomized versions of 𝖭𝖯-hard problems such as optimizing the Sherrington-Kirkpatrick model [52] have been proven to be average-case easy [41].

5 Open problems

There are multiple avenues to continue this work. Here, we list a couple:

  • Mulmuley, Vazirani, and Vazirani [43] introduced a variant of the isolation lemma in [59], by introducing random weights to the edges of a graph to isolate a unique maximally weighted clique. The quantum analogue of this approach is to consider the quantum variant of the Clique problem, which is 𝖰𝖬𝖠 hard [14]. It seems to suffer with issues analogous to those raised in [9]: introducing poly(n) amount of randomness does not suffice to single out a unique quantum witness. Can the ideas in [43] still be useful in some quantum problems?

  • Theorem 2 demonstrates a quantum oracle to which 𝖰𝖷𝖢𝒪𝖰𝖬𝖠𝖰𝖬𝖠𝒪. We expect our tools can be extended to showing separations for constants stacks of 𝖰𝖬𝖠 (i.e. for the 𝖰𝖬𝖠 hierarchy [2])? Extending our argument to, e.g., 𝖰𝖬𝖠𝖰𝖬𝖠𝖰𝖬𝖠 would require 1) showing polarization for 𝖰𝖬𝖠𝖰𝖬𝖠 circuits, and 2) a two-sided query lower bound, along the lines of [51], but for relativized classes.

  • Theorem 3 constructs a quantum algorithm that verifies some low energy state under ETH. Can we verify a more natural low energy state assuming ETH, such as the purification of Gibbs state (called the thermofield double state)? Note that, due to the Feynman-Kitaev circuit-to-Hamiltonian mapping [36], uniquely verifying a state is equivalent to the state being a near ground state of a local Hamiltonian with inverse polynomial spectral gap. Such local Hamiltonians for the thermofield double state have been constructed under physics-motivated assumptions in [20].

  • Our results hint at a separation between worst-case and average-case local Hamiltonians. Random spin Hamiltonians [24] or the SYK model [49, 34, 35] are expected to satisfy the ETH and are therefore solvable by a 𝖴𝗇𝗂𝗊𝗎𝖾𝖰𝖬𝖠 machine. It is interesting to find more evidence for the easiness of average-case Hamiltonians. Recent work studies the complexity of approximating the maximal energy of the SYK model. Ref. [30] provides a quantum algorithm to prepare a state with energy cn – a constant fraction of the expected optimum. Ref. [11] finds further evidence that estimating the energy of the SYK might even be quantumly easy. More concretely, they show that the annealed and quenched free energies agree at inverse polynomial temperatures. This is in contrast to classical random spin systems where these quantities disagree when the system is in a “glassy” phase that is algorithmically hard.

  • It is interesting to find more connections between quantum many-body physics assumptions such as ETH and quantum computer science. The works [19, 54] show that ETH implies fast preparation of Gibbs quantum states, if there is no bottleneck throughout the eigenvalue distribution and Gibbs distribution. However, the latter assumption cannot be justified on general grounds. On the other hand, our results work perfectly fine without this assumption.

References

  • [1] Scott Aaronson. On perfect completeness for qma. arXiv preprint arXiv:0806.0450, 2008.
  • [2] Scott Aaronson, DeVon Ingram, and William Kretschmer. The acrobatics of bqp. In 37th Computational Complexity Conference, 2022.
  • [3] Scott Aaronson, Robin Kothari, William Kretschmer, and Justin Thaler. Quantum lower bounds for approximate counting via laurent polynomials. In Shubhangi Saraf, editor, 35th Computational Complexity Conference, CCC 2020, Saarbrücken, Germany (Virtual Conference), July 28-31, 2020, volume 169 of LIPIcs, pages 7:1–7:47. Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2020. doi:10.4230/LIPIcs.CCC.2020.7.
  • [4] Scott Aaronson and Greg Kuperberg. Quantum versus classical proofs and advice. Theory Comput., 3(1):129–157, 2007. doi:10.4086/TOC.2007.V003A007.
  • [5] Ian Affleck, Tom Kennedy, Elliott H Lieb, and Hal Tasaki. Rigorous results on valence-bond ground states in antiferromagnets. Condensed Matter Physics and Exactly Soluble Models: Selecta of Elliott H. Lieb, pages 249–252, 2004.
  • [6] Avantika Agarwal and Shalev Ben-David. Oracle separations for the quantum-classical polynomial hierarchy. arXiv preprint arXiv:2410.19062, 2024. doi:10.48550/arXiv.2410.19062.
  • [7] Avantika Agarwal and Srijita Kundu. A cautionary note on quantum oracles. arXiv preprint arXiv:2504.19470, 2025. doi:10.48550/arXiv.2504.19470.
  • [8] Dorit Aharonov, Itai Arad, and Thomas Vidick. Guest column: the quantum pcp conjecture. Acm sigact news, 44(2):47–79, 2013. doi:10.1145/2491533.2491549.
  • [9] Dorit Aharonov, Michael Ben-Or, Fernando G. S. L. Brandão, and Or Sattath. The pursuit of uniqueness: Extending valiant-vazirani theorem to the probabilistic and quantum settings. Quantum, 6:668, 2022. doi:10.22331/Q-2022-03-17-668.
  • [10] Dorit Aharonov, Aram W. Harrow, Zeph Landau, Daniel Nagaj, Mario Szegedy, and Umesh V. Vazirani. Local tests of global entanglement and a counterexample to the generalized area law. In 55th IEEE Annual Symposium on Foundations of Computer Science, FOCS 2014, Philadelphia, PA, USA, October 18-21, 2014, pages 246–255. IEEE Computer Society, 2014. doi:10.1109/FOCS.2014.34.
  • [11] Eric R Anshuetz, Chi-Fang Chen, Bobak T Kiani, and Robbie King. Strongly interacting fermions are non-trivial yet non-glassy. arXiv preprint arXiv:2408.15699, 2024.
  • [12] Itai Arad, Tomotaka Kuwahara, and Zeph Landau. Connecting global and local energy distributions in quantum spin models on a lattice. Journal of Statistical Mechanics: Theory and Experiment, 2016(3):033301, 2016. doi:10.1088/1742-5468/2016/03/033301.
  • [13] Richard Beigel, Harry Buhrman, and Lance Fortnow. Np might not be as easy as detecting unique solutions. In Proceedings of the thirtieth annual ACM symposium on Theory of computing, pages 203–208, 1998. doi:10.1145/276698.276737.
  • [14] Salman Beigi and Peter W. Shor. On the complexity of computing zero-error and holevo capacity of quantum channels, 2008. arXiv:0709.2090.
  • [15] Shai Ben-David, Benny Chor, and Oded Goldreich. On the theory of average case complexity. In Proceedings of the twenty-first annual ACM symposium on Theory of computing, pages 204–216, 1989.
  • [16] Sergey Bravyi, Anirban Chowdhury, David Gosset, Vojtech Havlicek, and Guanyu Zhu. Quantum complexity of the kronecker coefficients. PRX Quantum, 5(1):010329, 2024.
  • [17] Sergey Bravyi, Anirban Chowdhury, David Gosset, and Pawel Wocjan. On the complexity of quantum partition functions. CoRR, abs/2110.15466, 2021. arXiv:2110.15466.
  • [18] Brielin Brown, Steven T Flammia, and Norbert Schuch. Computational difficulty of computing the density of states. Physical review letters, 107(4):040501, 2011.
  • [19] Chi-Fang Chen and Fernando G. S. L. Brandão. Fast Thermalization from the Eigenstate Thermalization Hypothesis. arXiv preprint arXiv:2112.07646, 2023.
  • [20] William Cottrell, Ben Freivogel, Diego M. Hofman, and Sagar F. Lokhande. How to build the thermofield double state. Journal of High Energy Physics, 2019(2):58, February 2019. doi:10.1007/JHEP02(2019)058.
  • [21] Luca D’Alessio, Yariv Kafri, Anatoli Polkovnikov, and Marcos Rigol. From quantum chaos and eigenstate thermalization to statistical mechanics and thermodynamics. Advances in Physics, 65(3):239–362, 2016.
  • [22] Abhinav Deshpande, Alexey V Gorshkov, and Bill Fefferman. Importance of the spectral gap in estimating ground-state energies. PRX Quantum, 3(4):040327, 2022.
  • [23] Josh M Deutsch. Quantum statistical mechanics in a closed system. Physical review a, 43(4):2046, 1991.
  • [24] László Erdős and Dominik Schröder. Phase transition in the density of states of quantum spin glasses. Mathematical Physics, Analysis and Geometry, 17(3):441–464, 2014.
  • [25] Friedrich Götze and Holger Sambale. Higher order concentration on stiefel and grassmann manifolds. Electronic Journal of Probability, 28:1–30, 2023.
  • [26] Joachim Grollmann and Alan L Selman. Complexity measures for public-key cryptosystems. SIAM Journal on Computing, 17(2):309–335, 1988. doi:10.1137/0217018.
  • [27] David Gross and Jens Eisert. Quantum margulis expanders. arXiv preprint arXiv:0710.0651, 2007.
  • [28] Dominik Hangleiter and Jens Eisert. Computational advantage of quantum random sampling. Reviews of Modern Physics, 95(3):035001, 2023.
  • [29] Matthew B Hastings. Random unitaries give quantum expanders. Physical Review A—Atomic, Molecular, and Optical Physics, 76(3):032315, 2007.
  • [30] Matthew B Hastings and Ryan O’Donnell. Optimizing strongly interacting fermionic hamiltonians. In Proceedings of the 54th Annual ACM SIGACT Symposium on Theory of Computing, pages 776–789, 2022. doi:10.1145/3519935.3519960.
  • [31] Sandy Irani, Anand Natarajan, Chinmay Nirkhe, Sujit Rao, and Henry Yuen. Quantum search-to-decision reductions and the state synthesis problem. arXiv preprint arXiv:2111.02999, 2021. arXiv:2111.02999.
  • [32] Rahul Jain, Iordanis Kerenidis, Greg Kuperberg, Miklos Santha, Or Sattath, and Shengyu Zhang. On the power of a unique quantum witness. Theory Comput., 8(1):375–400, 2012. doi:10.4086/TOC.2012.V008A017.
  • [33] Mark R Jerrum, Leslie G Valiant, and Vijay V Vazirani. Random generation of combinatorial structures from a uniform distribution. Theoretical computer science, 43:169–188, 1986. doi:10.1016/0304-3975(86)90174-X.
  • [34] Alexei Kitaev. A simple model of quantum holography (part 1). Entanglement in strongly-correlated quantum matter, page 38, 2015.
  • [35] Alexei Kitaev. A simple model of quantum holography (part 2). Entanglement in strongly-correlated quantum matter, page 38, 2015.
  • [36] Alexei Y. Kitaev, A. H. Shen, and Mikhail N. Vyalyi. Classical and Quantum Computation, volume 47 of Graduate studies in mathematics. American Mathematical Society, 2002. URL: https://bookstore.ams.org/gsm-47/.
  • [37] Adam R Klivans and Dieter Van Melkebeek. Graph nonisomorphism has subexponential size proofs unless the polynomial-time hierarchy collapses. In Proceedings of the Thirty-First Annual ACM Symposium on Theory of Computing, pages 659–667, 1999. doi:10.1145/301250.301428.
  • [38] Tohru Koma and Bruno Nachtergaele. The spectral gap of the ferromagnetic xxz-chain. Letters in Mathematical Physics, 40(1):1–16, 1997.
  • [39] Greg Kuperberg. How hard is it to approximate the jones polynomial? arXiv preprint arXiv:0908.0512, 2009. arXiv:0908.0512.
  • [40] Zeph Landau, Umesh Vazirani, and Thomas Vidick. A polynomial time algorithm for the ground state of one-dimensional gapped local hamiltonians. Nature Physics, 11(7):566–569, 2015.
  • [41] Andrea Montanari. Optimization of the sherrington-kirkpatrick hamiltonian. SIAM J. Comput., 54(4):S19–1, 2025. doi:10.1137/20M132016X.
  • [42] Ramis Movassagh and Peter W Shor. Power law violation of the area law in quantum spin chains. arXiv preprint arXiv:1408.1657, 2014.
  • [43] Ketan Mulmuley, Umesh V. Vazirani, and Vijay V. Vazirani. Matching is as easy as matrix inversion. Combinatorica, 7(1):105–113, March 1987. doi:10.1007/BF02579206.
  • [44] Rahul Nandkishore and David A Huse. Many-body localization and thermalization in quantum statistical mechanics. Annu. Rev. Condens. Matter Phys., 6(1):15–38, 2015.
  • [45] Anand Natarajan and Chinmay Nirkhe. A distribution testing oracle separation between qma and qcma. Quantum, 8:1377, 2024. doi:10.22331/Q-2024-06-17-1377.
  • [46] Anand Natarajan and Thomas Vidick. Robust self-testing of many-qubit states. arXiv preprint arXiv:1610.03574, 2016. arXiv:1610.03574.
  • [47] Román Orús. A practical introduction to tensor networks: Matrix product states and projected entangled pair states. Annals of physics, 349:117–158, 2014.
  • [48] Marcos Rigol, Vanja Dunjko, and Maxim Olshanii. Thermalization and its mechanism for generic isolated quantum systems. Nature, 452(7189):854–858, 2008.
  • [49] Subir Sachdev and Jinwu Ye. Gapless spin-fluid ground state in a random quantum heisenberg magnet. Physical review letters, 70(21):3339, 1993.
  • [50] Norbert Schuch, Michael M. Wolf, Frank Verstraete, and J. Ignacio Cirac. Computational complexity of projected entangled pair states. Phys. Rev. Lett., 98:140506, April 2007. doi:10.1103/PhysRevLett.98.140506.
  • [51] Adrian She and Henry Yuen. Unitary property testing lower bounds by polynomials. In Leibniz International Proceedings in Informatics (LIPIcs): 14th Innovations in Theoretical Computer Science Conference (ITCS 2023). Schloss Dagstuhl – Leibniz-Zentrum für Informatik, 2023. doi:10.4230/LIPIcs.ITCS.2023.96.
  • [52] David Sherrington and Scott Kirkpatrick. Solvable model of a spin-glass. Physical review letters, 35(26):1792, 1975.
  • [53] Yaoyun Shi and Shengyu Zhang. Note on quantum counting classes, 2009. URL: http://www.cse.cuhk.edu.hk/˜syzhang/papers/SharpBQP.pdf.
  • [54] Oles Shtanko and Ramis Movassagh. Preparing thermal states on noiseless and noisy programmable quantum processors, 2023. arXiv:2112.14688.
  • [55] Julian Sonner and Manuel Vielma. Eigenstate thermalization in the sachdev-ye-kitaev model. Journal of High Energy Physics, 2017(11):149, November 2017. doi:10.1007/JHEP11(2017)149.
  • [56] Mark Srednicki. The approach to thermal equilibrium in quantized chaotic systems. Journal of Physics A: Mathematical and General, 32(7):1163, 1999.
  • [57] Larry Stockmeyer. The complexity of approximate counting. In Proceedings of the fifteenth annual ACM symposium on Theory of computing, pages 118–126, 1983.
  • [58] Seinosuke Toda. Pp is as hard as the polynomial-time hierarchy. SIAM Journal on Computing, 20(5):865–877, 1991. doi:10.1137/0220053.
  • [59] Leslie G. Valiant and Vijay V. Vazirani. NP is as easy as detecting unique solutions. Theor. Comput. Sci., 47(3):85–93, 1986. doi:10.1016/0304-3975(86)90135-0.
  • [60] Mark Zhandry. Toward Separating QMA from QCMA with a Classical Oracle, 2024. arXiv:2411.01718 [quant-ph]. doi:10.48550/arXiv.2411.01718.