Abstract 1 Introduction 2 Preliminaries 3 XZ-Quantum-6-Sat is in QMA1 4 Circuit-to-XZ-Hamiltonian reduction 5 XZ-Quantum-6-Sat is QMA1-hard 6 Soundness: technical lemmas References

Two Bases Suffice for QMA1-Completeness

Henry Ma ORCID CSAIL, Massachusetts Institute of Technology, Cambridge, MA, USA Anand Natarajan ORCID CSAIL, Massachusetts Institute of Technology, Cambridge, MA, USA
Abstract

We introduce a basis-restricted variant of the Quantum-k-Sat problem, in which each term in the input Hamiltonian is required to be diagonal in either the standard or Hadamard basis. Our main result is that the Quantum-6-Sat problem with this basis restriction is already 𝖰𝖬𝖠1-complete, defined with respect to a natural gateset. Our construction is based on the Feynman-Kitaev circuit-to-Hamiltonian construction, with a modified clock encoding that interleaves two clocks in the standard and Hadamard bases. In light of the central role played by CSS codes and the uncertainty principle in the proof of the NLTS theorem of Anshu, Breuckmann, and Nirkhe (STOC ’23), we hope that the CSS-like structure of our Hamiltonians will make them useful for progress towards a quantum PCP theorem.

Keywords and phrases:
quantum complexity theory, Hamiltonian complexity, Quantum Merlin Arthur (QMA), QMA1, quantum satisfiability problem
Funding:
Henry Ma: HM was supported by the NSF Graduate Research Fellowship Program.
Anand Natarajan: AN was supported by NSF CAREER grant number 2339948.
Copyright and License:
[Uncaptioned image] © Henry Ma and Anand Natarajan; licensed under Creative Commons License CC-BY 4.0
2012 ACM Subject Classification:
Theory of computation Complexity classes
; Theory of computation Problems, reductions and completeness ; Theory of computation Quantum complexity theory
Related Version:
Full Version: https://arxiv.org/abs/2509.24390
Acknowledgements:
Part of this work was done while both authors were visiting the Simons Institute for the Theory of Computing as part of the 2025 Summer Cluster on Quantum Computing, and the Challenge Institute for Quantum Computation at UC Berkeley. We thank Chinmay Nirkhe for several helpful discussions.
Editor:
Shubhangi Saraf

1 Introduction

Local Hamiltonians play a central role in quantum complexity theory, tying the subject both to applications in condensed matter physics and quantum chemistry, and to the classical theory of 𝖭𝖯-completeness, constraint satisfaction problems, and combinatorial optimization. Much of the complexity theory of local Hamiltonians has developed in analogy to the theory of 𝖭𝖯-completeness, with the class 𝖰𝖬𝖠 playing the role of 𝖭𝖯, and Kitaev’s result showing that the local Hamiltonian problem (estimating the ground energy of a local Hamiltonian) is 𝖰𝖬𝖠-complete playing the role of the Cook-Levin theorem (𝖭𝖯-completeness of 3SAT, and local constraint satisfaction problems more generally). Despite these analogies, there are considerable differences between the classical and quantum settings. One important difference relates to the issue of “perfect completeness.” In the case of classical constraint satisfaction problems, it is typically just as hard to solve the “SAT” problem of deciding whether all constraints are simultaneously satisfiable, as to solve the “decisional MAX-SAT” problem of deciding whether at least a certain fraction of the constraints are satisfiable: both problems are 𝖭𝖯-complete. In contrast, the equivalent SAT and MAX-SAT problems for local Hamiltonians appear to have different complexities: all of our known 𝖰𝖬𝖠-completeness results are MAX-SAT type results, while SAT-type problems appear to capture the complexity class 𝖰𝖬𝖠1 of 𝖰𝖬𝖠 proof systems with perfect completeness (i.e. the verifier accepts good proofs with certainty). By definition, 𝖰𝖬𝖠1𝖰𝖬𝖠, and the two classes are believed to be qualitatively similar in power, but the exact relation between them is murky.

Due to the nature of our tools for showing computational hardness of local Hamiltonian problems, our understanding of the MAX-SAT setting is better than our understanding of the SAT setting: in the MAX-SAT world, we have dichotomy theorems [10] for 2-local Hamiltonians, and a wide variety of families of Hamiltonians are known to be 𝖰𝖬𝖠-complete, including very simple and natural models like the XZ-model [6]. In the SAT setting, while we know that Quantum-3-Sat (the SAT problem for 3-local Hamiltonians) is 𝖰𝖬𝖠1-complete [13], and that various other cases of Quantum-k-Sat are easier (e.g. for k=2 the problem is in 𝖯 [7], and for so-called “stoquastic” Hamiltonians the problem is in 𝖬𝖠 [8]), we don’t have dichotomy theorems, nor do we have hardness results as clean as the XZ-model hardness mentioned above. This is because tools like gadget reductions used to reduce between different types of Hamiltonians do not preserve perfect completeness.

At the same time, the class 𝖰𝖬𝖠1 is scientifically of great interest as a testbed to better understand Hamiltonian complexity. Arguably the single most important open problem in Hamiltonian complexity is the quantum PCP conjecture: that approximating the ground energy of a local Hamiltonian even up to a constant fraction of the total norm of the Hamiltonian remains 𝖰𝖬𝖠-hard. A natural weakening of this conjecture would be to show 𝖰𝖬𝖠1-hardness for this problem, and in fact, the biggest partial milestone we have towards a quantum PCP theorem comes from the SAT setting. Specifically, the NLTS theorem of Anshu, Breuckmann, and Nirkhe [5] shows the existence of local Hamiltonians whose low-energy states satisfy a certain structural property that we expect should hold for quantum PCP Hamiltonians. These Hamiltonians are parent Hamiltonians of quantum stabilizer error correcting codes, and thus are exactly satisfiable (all valid code states are 0-energy ground states) – and this characterization of the ground space is used explicitly in the proof of the NLTS theorem. However, for the same reason, these Hamiltonians contain no computational hardness. The most natural next step beyond the NLTS theorem would be to show that the NLTS structural property holds for a family of Hamiltonians where the local Hamiltonian problem is computationally hard. In order to achieve such a result, can we show a 𝖰𝖬𝖠1-hardness result for a class of local Hamiltonians that is sufficiently structured to use the techniques of Anshu, Breuckmann and Nirkhe’s proof?

In this work, we focus on one structural property of the code Hamiltonians of [5]: every local term in the Hamiltonian is a projector and is diagonal in either the standard basis (the “Pauli Z-basis”), or the Hadamard basis (the “Pauli X-basis”). Stabilizer codes with this property are called CSS codes, and the CSS structure of the code is crucial to the proof of the NLTS theorem in two ways. Firstly, to prove the NLTS property of low-energy states, Anshu et al, following the approach of Eldar and Harrow [12], reduce the problem to proving that low-energy states of their Hamiltonian, when measured in either the standard or Hadamard bases, give rise to “well spread” probability distributions over the Boolean hypercube. They are able to prove this property in turn using a Heisenberg uncertainty principle relating these two bases, together with properties of the codes they consider. Secondly, these code properties in turn are formulated and proved by viewing CSS codes as consisting of a pair of classical error correcting codes, one in each basis. In fact, the view of CSS codes as chain complexes, which is at the foundation of the vast literature on constructing and analyzing such codes, relies on this two-basis structure.

The main result of this work is to show that the SAT problem for two-basis local Hamiltonians is 𝖰𝖬𝖠1-complete. More precisely, our main theorem is the following.

Theorem 1.1.

Let XZ-Quantum-6-Sat be the following problem: given a local Hamiltonian H=iHi on n qubits, where each local term Hi is a projector acting on 6 qubits and is diagonal in either the standard or Hadamard basis, determine whether λmin(H)=0 or λmin(H)b(n) where b(n)=Ω(1/poly(n)), promised that one of the two is the case. Then for an appropriately chosen function b(n), it holds that XZ-Quantum-6-Sat is complete for the class 𝖰𝖬𝖠1𝒢2 of quantum Merlin-Arthur proof systems with perfect completeness where the verifier consists of a circuit made up of gates from the gate set 𝒢2.

There is a technical subtlety in the theorem statement, arising from the dependence of 𝖰𝖬𝖠1 on the choice of gate set for the verifier’s circuit. We show completeness with respect to the gate set 𝒢2 as defined by Rudolph [18], consisting of NOT, controlled-NOT, and Toffoli gates (X,CX,CCX) and the tensor product H^H^ of two Hadamard gates. This gate set is a slight variant of the commonly used Hadamard-plus-Toffoli gate set [1], which is universal for quantum computing.

To better understand the class of Hamiltonians for which we show hardness, it is perhaps useful to draw a classical analogy. A classical linear code is associated with a collection of linear constraints over 𝔽2. By the Gaussian elimination algorithm, the SAT problem for linear constraints can always be solved in polynomial time, but once the constraints are allowed to be nonlinear, the complexity of the problem jumps up to 𝖭𝖯. In the quantum world, the parent Hamiltonian of a CSS code consists of a pair of linear classical constraint satisfaction problems (CSPs), one for each basis: a valid code state is one that, when measured in the Z-basis, yields a satisfying string for the Z-basis linear constraints, and when measured in the X-basis, yields a satisfying string for the X-basis linear constraints. By standard stabilizer techniques, Gaussian elimination is sufficient to solve the SAT problem for these Hamiltonians in polynomial time. The Hamiltonians we consider also consist of a pair of classical constraint satisfaction problems in the two bases, except with nonlinear111To be clear, the constraints are nonlinear in their action on the binary string measurement outcomes viewed as elements of 𝔽2n – the Hamiltonian is still a linear operator over the Hilbert space. constraints now being allowed, and our result shows that, as in the classical case, allowing for nonlinearity causes the complexity to jump from 𝖯 to the maximal possible level of hardness (𝖰𝖬𝖠1 in this case). This classical analogy also helps illustrate the relation between our result and the MAX-SAT case. In the MAX-SAT case, classically, even two-local linear constraints become 𝖭𝖯-hard (this is by the 𝖭𝖯-hardness of the MAX-CUT problem), and similarly, in the quantum case, the result of Biamonte and Love for the XZ-model shows that the MAX-SAT problem for “linear” constraints is 𝖰𝖬𝖠-hard.

Constraints SAT MAX-SAT
Z-basis, linear In 𝖯 (by Gaussian elimination) 𝖭𝖯-hard (by MAX-CUT)
Z-basis, general 𝖭𝖯-hard (Cook-Levin) 𝖭𝖯-hard (by )
Z and X bases, linear In 𝖯 (by Gaussian elimination) 𝖰𝖬𝖠-hard [6]
Z and X bases, general 𝖰𝖬𝖠1-hard (this work) 𝖰𝖬𝖠-hard (by )

The two basis paradigm

Viewing matters more subjectively and at a higher level, it is a striking fact that many interesting phenomena in quantum computing rest on the interplay between the standard and Hadamard bases. Examples of this “two basis paradigm” in action include the BB84 protocol and Wiesner’s quantum money scheme, Simon’s algorithm, the magic square game (and, arguably, the proof that 𝖬𝖨𝖯=𝖱𝖤), the forrelation problem (some versions of which are 𝖡𝖰𝖯-complete), Aaronson and Christiano’s subspace quantum money scheme, and Mahadev’s measurement protocol, besides the previously mentioned examples of CSS codes and Biamonte and Love’s QMA-hardness results. We see our result as another instance of this paradigm.

1.1 Technical overview

Our proof of ˜1.1 consists of two parts. To show 𝖰𝖬𝖠1𝒢2-completeness of our problem, we must show that it is both contained in 𝖰𝖬𝖠1𝒢2 and that it is 𝖰𝖬𝖠1𝒢2-hard. The containment follows along standard lines: we construct a verifier for XZ-Quantum-6-Sat instances that samples a random term from the Hamiltonian and coherently measures it on the witness state. The only nontrivial step in this construction is showing that the verifier’s gate set 𝒢2 can exactly simulate a controlled Hadamard gate, since this gate is needed to perform the measurement of the witness in the appropriate basis.

The bulk of the proof is concerned with the 𝖰𝖬𝖠1𝒢2-hardness. We build on the Kitaev circuit-to-Hamiltonian construction [15], which is a quantum analog of the fundamental construction used by Cook and Levin [9, 16] to show that Sat is 𝖭𝖯-complete. Let us consider a circuit consisting of gates U1,U2,,UT, acting on a Hilbert space . Kitaev’s construction associates this circuit to a Hamiltonian H acting on a space clock, where clock is the Hilbert space of an ancillary “clock” register. The space clock contains at least T+1 orthonormal states |0^,|1^,,|T^, but typically has a much larger dimension – 2T in Kitaev’s original construction. The full Hamiltonian H can be written as a sum of four components:

H=Hprop+Hformat+Hin+Hout,

where each component is a sum of local projectors. Moreover, Kitaev’s analysis shows that, when applied to a 𝖰𝖬𝖠1 circuit U1,,UT that accepts some witness state with certainty, the Hamiltonian H has a 0- energy ground state, and otherwise, if applied to a circuit that rejects all witnesses with high probability, then the minimum energy of H is at least 1/poly(n).

Our approach is to use Kitaev’s construction, but modify the terms Hprop and Hformat, as well as the encoding of the clock states |0^,,|T^, so that the resulting Hamiltonian H satisfies our two-basis constraint. Since we will not modify them significantly, for now let us ignore the terms Hinit and Hout, and study the ground space of the remaining two terms. These are designed in Kitaev’s construction so that they always have a nonempty 0-energy ground space, which consists of history states of the form

|ψhistory=1T+1t=0TUtU1|ψ0|t^.

This is achieved in two ways:

  • The “format” Hamiltonian Hformat forces the state to be supported only on the “good” clock register states |0^,,|T^. This is done in Kitaev’s construction by taking clock to be a space of T qubits, and taking the good clock states to be encodings of 0,,T in unary in the standard basis, so that

    |t^=|11|1t1|1t|0t+1|0T.

    To force the state to be supported only on these states, Hformat consist of projectors that enforce the constraint on each pair of adjacent qubits (t,t+1) that a 0 can never be followed by a 1. These constraints are already diagonal in the Z-basis.

  • The “propagation” Hamiltonian Hprop consists of a sum of T constraints, each of which forces the |t1^ and |t^ components of the state to be related by an application of the unitary Ut. In Kitaev’s construction, these terms take the form

    Hprop,t =12I(|110110|+|100100|)t1,t,t+1
    12Ut(|110100|)t1,t,t+112Ut(|100110|)t1,t,t+1,

    where the second tensor factor acts on the specified qubits of the clock register. It can be checked that for general choices of Ut, these terms are not diagonal in either the X-basis or the Z-basis.

To make progress, we must modify the propagation terms. To do this, we make the following key observation: the terms Hprop,t are “almost” diagonal in the X-basis whenever the unitary Ut itself is diagonal in the X-basis and is Hermitian (i.e. Ut=Ut). Specifically, in this case, Hprop,t can be written as

Hprop,t=12(IItUtXt)(|1010|)t1,t+1,

which is diagonal when all qubits are placed in the X-basis except for qubits t1,t+1 of the clock register, which remain in the Z-basis. Moreover, the same property would hold, with X and Z interchanged, if the good clock states were X-basis states, and Ut were Hermitian and diagonal in the Z-basis.

This suggests the following modifications:

  • Encode good clock states by putting alternating qubits in the X- and Z-bases. This means that valid clock states would look like

    |0^=|0+0+0+,|1^=|1+0+0+,|2^=|1-0+0+,|3^=|1-1+0+,.
  • Choose a universal gate set where each gate is Hermitian and diagonal in either the X- or Z-basis. We construct such a gate set and show that it is computationally equivalent to 𝒢2.

  • Padding the circuit so that each X-basis gate falls on an odd timestep and each Z-basis timestep falls on an even timestep. This makes the propagation terms fully diagonal in the X-basis for odd times, and the Z-basis for even times: e.g. for odd t,

    Hprop,t=12(IItUtXt)|-+-+|t1,t+1.

With these changes, every term in Hprop,t is diagonal in either the X- or Z-basis. However, now that we have changed the encoding of the clock, we must also change Hformat, and it is not hard to see that any Hamiltonian that restricts us only to valid clock states will not be diagonal in either basis. The solution to this is somewhat surprising.

  • We modify Hformat to include checks that only check consistency separately on the X- and Z-basis parts of the clock state. Specifically, our new Hformat checks that the clock register in a state |tZ,tX consisting of an “interleaving” of some valid unary-encoding of an integer tZ in the Z-basis in odd positions, and an integer tX in the X-basis in even positions. It does not check that these two interleaved times are synchronized with each other, meaning that this Hamiltonian accepts invalid sates such as

    |0-0-0+0+0,

    which do not correspond to valid clock states. We call these states “fake” states.

  • Surprisingly, we show that, in fact, Hprop+Hformat together force the ground state to supported only on valid clock states. This is because every fake clock state is coupled by some term in Hprop to a “bad” clock state, which incurs an energy penalty from Hformat. We analyze this quantitatively by relating the Hamiltonian to a graph Laplacian.

The majority of the technical work in our analysis consists in (a) proving universality of our new gate set, and (b) quantitatively bounding the minimum energy of Hformat+Hprop on the subspace of invalid clock strings.

1.2 Open questions

We see several possibilities for future work with a quantum PCP flavor.

  1. 1.

    Our construction, since it is based on the Feynman-Kitaev clock, cannot achieve the NLTS property: indeed, the product state |0|0^ violates only a single propagation term of the Hamiltonian. To get to NLTS while maintaining the two-basis structure, can we apply our ideas to the tensor-network-based circuit-to-Hamiltonian construction of [4]?

  2. 2.

    What do random instances of our two-basis Hamiltonians look like? Can we find a distribution such that Hamiltonians sampled from this distribution have a zero-energy ground state with high probability, but where finding the ground state is computationally hard? It would be especially interesting to relate this to work on the combinatorial NLTS property for tensor network Hamiltonians constructed from random SAT instances [3].

  3. 3.

    One candidate approach to proving the quantum PCP conjecture is to design a “locality-preserving gap amplification” procedure, that operates on instances of local Hamiltonians by increasing the minimum energy of unsatisfiable Hamiltonians (amplifying the “promise gap”), while preserving the locality of the Hamiltonian. To build towards such a procedure, can we apply classical “gap-amplification” or “locality-reduction” transformations separately to one of the two bases and make some kind of meaningful progress towards gap amplification? Moreover, is there an analog of “distance balancing” for CSS codes [20], whereby a code with high distance in one basis but low distance in the other can be converted into a code with decent distance in both bases?

  4. 4.

    In this paper we have taken the point of view that 𝖰𝖬𝖠1𝒢2 is likely to be qualitatively similar to 𝖰𝖬𝖠 in power. But if 𝖰𝖬𝖠1𝒢2 is in fact much weaker, could our completeness result give us a route to showing this by putting the problem XZ-Quantum-6-Sat into a smaller class like 𝖰𝖢𝖬𝖠? One very speculative route to doing this is to give an efficient classical description of ground states of such Hamiltonians, perhaps using tools from additive combinatorics, since Hamiltonian terms in the X-basis can be viewed as additive constraints on the support of the ground state in the standard basis.

2 Preliminaries

2.1 Quantum computation

We briefly set up the basic formalism of quantum computation. For a more detailed exposition, we refer the reader to [17].

An n-qubit quantum state |ψ is a unit vector in a complex Hilbert space (2)n. A k-qubit quantum operation is a unitary operator U acting on a k-qubit space (2)k. We follow the conventions of Dirac notation. We use U to denote the adjoint of an operator. Let [n]={1,,n}. We can extend U to an operation UI[n]S on n-qubit space, where S[n] is the set of k qubits which U acts on, and I[n]S is the identity operator on the remaining qubits.

We now fix some notation for relevant states and operations. As usual, |0 and |1 refer to the single-qubit standard basis states, while |+ and | are the Hadamard basis states. We will often omit tensor products, e.g. |00=|0|0. Let X and Z be the corresponding single-qubit Pauli operations. We will also refer to the standard and Hadamard bases as the Z basis and X basis respectively. Let H^ be the Hadamard gate (we put the hat to avoid confusion with a Hamiltonian H). Let CX, CZ, and SWAP denote the two-qubit CNOT, controlled-Z and swap operations respectively. Let CCX denote the Toffoli gate and CCZ the controlled-controlled-Z gate.

A quantum circuit on n qubits is an n-qubit operation U which can be decomposed into a sequence of gates U=GmG1, where each Gi is a quantum operation from some fixed gate set; typically, this gate set contains operations which each acts nontrivially on just a small number of qubits. In our protocols, we will consider projective measurements of the first qubit in the Z basis, which is given by the orthogonal projectors {|00|1I,|11|1I}, where the subscript indicates the first qubit register, and the identity operation acts on the remaining qubits.

2.2 Quantum complexity theory

We now define the model of quantum verification of interest in this work. A promise problem is a pair (Lyes,Lno) of subsets of bitstrings Lyes,Lno{0,1} which satisfies LyesLno=. Promise problems generalize languages, a notion from classical complexity theory: a language is just a promise problem which additionally satisfies LyesLno={0,1}.

A quantum verifier is a uniform family of quantum circuits {Vx} whose goal is to determine whether a given string x is in Lyes or Lno, promised that one is the case. On an input x of length n, the verifier is given access to a quantum proof state |ψ and a register of ancilla qubits all initialized to |0. The circuit is efficient, in the sense that both the proof and ancilla registers have poly(n) qubits, and there are poly(n) gates in the circuit.

After the circuit is applied, the final step of the verification is to measure the first qubit in the Z basis. We say the verifier accepts if the measurement outcome is 1 and rejects if the outcome is 0. The probability that the verifier Vx accepts is then

(|11|1I)Vx|ψ|00anc2. (1)

We now define 𝖰𝖬𝖠1 and 𝖰𝖬𝖠; we will mainly be interested in 𝖰𝖬𝖠1 in this work.

Definition 2.1.

A promise problem (Lyes,Lno) is in 𝖰𝖬𝖠1 if there is a uniform family of efficient quantum circuits {Vx} such that for any x{0,1}n,

  1. 1.

    If xLyes, then there is some proof state for which Vx accepts with probability 1.

  2. 2.

    If xLno, then for any proof state, Vx accepts with probability 1/3.

We refer to the first condition as completeness and the second condition as soundness. Importantly, the definition of 𝖰𝖬𝖠1 is not known to be independent of the choice of gate set for the circuits. In this work, we choose the gate set

𝒢2={X,CX,CCX,H^H^}

for 𝖰𝖬𝖠1. When needing to explicitly refer to 𝖰𝖬𝖠1 with this gate set, we use the notation 𝖰𝖬𝖠1𝒢2. We now make some observations about the chosen gate set.

Remark 2.2.

First, 𝒢2 is computationally universal, i.e. any quantum circuit can be simulated by a circuit which only uses gates from 𝒢2. This is because gates in 𝒢2 allow us to implement CCX and H^ (which can be implemented by applying H^H^ on the desired qubit and an otherwise unused ancilla qubit); these two gates already form a computationally universal gate set [19, 1].

Second, 𝒢2 is studied in a work of Rudolph [18] which makes progress on the interesting open problem of whether 𝖰𝖬𝖠1 has a universal gate set. They show that the problem of Gapped Clique Homology on weighted graphs (introduced in [14]) is 𝖰𝖬𝖠1𝒢2-complete. This result gives a surprising connection between an important problem in computational topology and 𝖰𝖬𝖠1𝒢2, motivating our study of other properties of 𝖰𝖬𝖠1𝒢2.

𝖰𝖬𝖠 is defined in the same way as 𝖰𝖬𝖠1, except the completeness parameter is 2/3, i.e. for xLyes, the verifier must accept with probability at least 2/3. In contrast to 𝖰𝖬𝖠1, the definition of 𝖰𝖬𝖠 is independent of the choice of gate set, since the Solovay-Kitaev theorem [11] shows that any quantum gate can be efficiently approximated using gates from a universal gate set. This argument does not straightforwardly extend to 𝖰𝖬𝖠1, since approximating a gate may not preserve the perfect completeness required of a 𝖰𝖬𝖠1 protocol.

The classical theory of 𝖭𝖯-completeness generalizes to 𝖰𝖬𝖠1 and 𝖰𝖬𝖠 in a straightforward way. A promise problem K=(Kyes,Kno) has an efficient reduction to a promise problem L=(Lyes,Lno) if there is an efficient deterministic algorithm A which maps a bitstring x to a bitstring A(x) such that A(x)Lyes if xKyes, and A(x)Lno if xKno. A promise problem L is 𝖰𝖬𝖠1-hard if for every K𝖰𝖬𝖠1, there is an efficient reduction from K to L. A promise problem L is 𝖰𝖬𝖠1-complete if L is in 𝖰𝖬𝖠1 and L is 𝖰𝖬𝖠1-hard. The definitions for 𝖰𝖬𝖠 are analogous.

2.3 Hamiltonian complexity

In this section, we define the Quantum-k-Sat problem and introduce its basis-restricted variant. An n-qubit Hamiltonian is a Hermitian operator acting on n qubits. A k-local operator on n qubits is an operator which acts nontrivially on at most k of the qubits. The energy of a state |ψ with respect to Hamiltonian H is ψ|H|ψ. Note that this quantity is always real and non-negative, since H is Hermitian.

Definition 2.3.

Let k1 and let 𝒮 be a set of Hermitian k-local projectors. The problem Quantum-k-Sat is a promise problem whose input is a classical description of a Hamiltonian H=i=1mhi, where each term hi acts on n qubits and belongs to 𝒮. H is a “yes” instance if there is an n-qubit state |ψ such that ψ|H|ψ=0. H is a “no” instance if for all n-qubit states |ψ, ψ|H|ψ1/poly(n).

The 𝖭𝖯-complete problem k-Sat can be viewed as a special case of Quantum-k-Sat with a highly restricted Hamiltonian: each projector term is diagonal in the same, predetermined basis (namely, the Z basis) [2]. Let us refer to this as a classical Hamiltonian. This brings into view one of the guiding questions of this work: if we relax the restrictions placed on a classical Hamiltonian, at what point does the Hamiltonian become “quantum”? This question has been previously studied for the relaxation in which the Hamiltonian is no longer required to be classical, but instead the Hamiltonian terms must pairwise commute. In this work, we introduce a new way of generalizing classical Hamiltonians which we call basis restriction.

Definition 2.4.

For n1, let n be a set of bases of n-qubit space. Let =n1n. Define -Quantum-k-Sat to be the problem Quantum-k-Sat where we choose 𝒮 (the set of allowed projectors) to be the set of Hermitian k-local projectors which are diagonal in some basis in .

As an example, if for all n we let n include just the Z basis on n qubits, then an instance of -Quantum-k-Sat is a Hamiltonian H=ihi for which each hi is diagonal in the Z basis. Then H is a classical Hamiltonian and -Quantum-k-Sat is in 𝖭𝖯.

We are interested in the complexity of basis-restricted Quantum-k-Sat when we allow for more than one “type” of basis. In particular, is there some setting in which the problem already becomes 𝖰𝖬𝖠-hard, even though the number of basis types is small? We answer this question in the following sections by considering the problem XZ-Quantum-k-Sat, which we define as basis-restricted Quantum-k-Sat problem where each hi is diagonal in either the Z or the X basis.

3 XZ-Quantum-6-Sat is in QMA1

In this section, we begin the proof of ˜1.1 by showing that XZ-Quantum-6-Sat is in 𝖰𝖬𝖠𝒢2. We start with a useful property of the gate set 𝒢2. For gate U, let Γ(U) denote the controlled-U gate.

Lemma 3.1.

For every U𝒢2, Γ(U) can be exactly implemented using gates in 𝒢2 (using an ancilla qubit, which can be in any state and will be left unchanged after the simulation).

Proof.

Γ(X) and Γ(CX) are already included in 𝒢2. Let a denote an ancilla qubit. Γ(CCX) on control qubits i,j,k and target qubit l can be decomposed as

Γ(CCX)i,j,k,lIa=CCXi,j,aCCXk,a,lCCXi,j,aCCXk,a,l.

For Γ(H^H^) with control qubit i and target qubits j,k, we first give a decomposition using the Γ(SWAP) gate.

Γ(H^H^)i,j,kIa=Γ(SWAP)i,j,k(H^H^)k,aΓ(SWAP)i,j,k(H^H^)k,a.

Then, since Γ(AB)=Γ(A)Γ(B) for any gates A and B, we may obtain a decomposition of Γ(SWAP) in terms of CCX’s by applying the identity SWAPi,j=CXi,jCXj,iCXi,j. This completes the proof. We illustrate these constructions in Figure 1.

Figure 1: Decomposition of Γ(CCX) and Γ(H^H^) into gates from 𝒢2.
Theorem 3.2.

XZ-Quantum-6-Sat is in 𝖰𝖬𝖠1𝒢2.

Proof.

First, without loss of generality, we can build a 𝖰𝖬𝖠1𝒢2 verification circuit which allows intermediate Z basis measurements and quantum operations conditioned on the measurement results. This is because a conditionally applied gate U can be replaced by the controlled operation Γ(U), removing the intermediate measurement. Lemma 3.1 then says that Γ(U) can then be decomposed back into the gate set 𝒢. Thus, a circuit with intermediate measurements can be rewritten as a circuit with just a single measurement at the end of the computation, as in our definition of 𝖰𝖬𝖠1. Also, classical processing can be performed freely, since our gate set includes CCX, which is universal for classical computation.

We use the verification procedure of Gosset and Nagaj [13], which is correct as long as the following condition holds:

(*)

There is an efficient quantum algorithm using gates in 𝒢 which, given a projector Π𝒮, exactly performs the projective measurement {IΠ,Π} on a given state |ψ.

Recall that 𝒮 is the set of Hermitian 6-local projectors which are diagonal in either the Z basis or the X basis.

We briefly sketch how the 𝖰𝖬𝖠1 protocol works assuming (*) holds. First, choose a random projector term Π in the Hamiltonian H (random bits can be obtained by applying X to a |0 ancilla then measuring). Then, perform the projective measurement in (*) on the proof state |ψ. Accept when the outcome is IΠ, and reject otherwise. The probability of accepting is 1ψ|Π|ψ. When H has a zero-energy state, the protocol accepts this state with probability 1, showing perfect completeness. If instead ψ|H|ψ1 for all |ψ, the protocol will reject with probability at least 1/poly(n), which can be amplified to give the desired soundness; the full analysis is given in [13].

We now prove that (*) holds. A projector Π𝒮 can be specified by a set of six qubits Q on which it acts nontrivially, a change of basis matrix V (which is either I or a tensor product of H^’s), and a set of supported strings S{0,1}6, such that

Π=zSV|zz|V.

The algorithm first applies V on |ψ, by applying H^ on each of the specified qubits (using the H^H^ gate, with the second H^ acting on an otherwise unused ancilla qubit). It then measures the qubits in Q in the Z basis. Let z be the measurement outcome. The algorithm returns outcome Π if zS, and otherwise returns outcome IΠ. The probability of getting outcome Π is

zS|zz|V|ψ2=ψ|Π|ψ,

as desired. This completes the proof.

4 Circuit-to-XZ-Hamiltonian reduction

In this section, we discuss some preliminaries for the 𝖰𝖬𝖠1-hardness proof of Section 5.

4.1 Two gate sets for QMA1

Recall that our gate set for 𝖰𝖬𝖠1 is 𝒢2={X,CX,CCX,H^H^}. Define a new gate set

𝒢XZ={X,CZ,CCZ,G^},

where G^ is defined as the two-qubit operation equivalent to

G^=(H^H^)CZ(H^H^).

We show that 𝒢XZ and 𝒢2 are interchangeable gate sets, in the following sense:

Lemma 4.1.

𝖰𝖬𝖠1𝒢2=𝖰𝖬𝖠1𝒢XZ.

Proof.

We start by showing that 𝖰𝖬𝖠1𝒢XZ𝖰𝖬𝖠1𝒢2. It suffices to show that each gate in 𝒢XZ can be exactly written as a (constant length) sequence of gates from 𝒢2. Such a simulation can freely use the 𝖰𝖬𝖠1 verifier’s ancilla register. In the sequences we construct, the ancilla qubits can be in any state and will be left unchanged after the gate simulation. Also note that exact simulation is critical: merely approximating a gate using another gate set would not necessarily preserve the 𝖰𝖬𝖠1 protocol’s perfect completeness.

To set some notation, for a gate U and a set of gates S, we say that US¯ if U can be exactly written as a sequence of gates from S. The following closure property is immediate: if SS¯ and US¯, then US¯. In particular, if VS¯ and US{V}¯, then US¯.

We now consider the gates in 𝒢XZ. The X gate is already in 𝒢2. We can simulate the CZ gate on qubits i,j with an ancilla qubit a, using the circuit identity

CZi,jIa=(H^H^)j,aCXi,j(H^H^)j,a.

This puts CZ𝒢2¯. A similar identity

CCZi,j,kIa=(H^H^)k,aCCXi,j,k(H^H^)k,a

puts CCZ𝒢2¯. Finally, G^𝒢2{CZ}¯ by definition, so by the closure property, G^𝒢2¯.

We now show that 𝖰𝖬𝖠1𝒢2𝖰𝖬𝖠1𝒢XZ using the same approach. We introduce a useful intermediate gate: define F^ to be the two-qubit operation equivalent to

F^=SWAP(H^H^).

We observe that F^=G^CZG^ using the circuit

[Uncaptioned image]=[Uncaptioned image]=[Uncaptioned image]=[Uncaptioned image],

where we used that H^jCZi,jH^j=CXi,j and SWAPi,j=CXi,jCXj,iCXi,j. This puts F^𝒢XZ¯. Next, using that H^H^ and SWAP commute, we have that

CXi,jIa=F^j,aCZi,aF^j,a,

since

[Uncaptioned image]=[Uncaptioned image]=[Uncaptioned image]=[Uncaptioned image].

Thus CX𝒢XZ{F^}¯, implying that CX𝒢XZ¯ by the closure property. The identity

CCXi,j,kIa=F^k,aCCZi,j,aF^k,a

similarly puts CCX𝒢XZ¯. Finally, note that H^H^=SWAPF^. Writing SWAP as three CX’s, we have that H^H^{F^,CX}¯, so H^H^𝒢XZ¯ by the closure property.

We find it useful to change to the gate set 𝒢XZ since each gate in 𝒢XZ has the following nice properties, which are straightforward to check:

Lemma 4.2.

Each gate in 𝒢XZ is Hermitian, 3-local, and diagonal in either the Z basis or the X basis.

These properties are used in the proof that XZ-Quantum-6-Sat is 𝖰𝖬𝖠1-hard, when we show that the Hamiltonian term Hprop is an instance of XZ-Quantum-6-Sat (˜4.4).

4.2 The Kitaev circuit-to-Hamiltonian construction

We now sketch the original Kitaev construction and note which parts of the construction carry through to the proof of ˜1.1 unchanged. Let (Lyes,Lno). be a promise problem in 𝖰𝖬𝖠1, and let Vx be the corresponding circuit which verifies a length n bitstring x, given a proof state |ψ. Recall that the verification starts from state |init=|ψ|00anc. Let denote the Hilbert space in which the initial state lives. Let Qproof and Qanc denote the indices of the qubits in the proof and ancilla registers. Write the circuit as Vx=UTU1, where each Ui is a gate from the chosen gate set for 𝖰𝖬𝖠1.

The Kitaev construction is an efficient mapping from the circuit Vx to a Hamiltonian

H=Hin+Hout+Hprop+Hformat.

H acts on Hilbert space clock, where clock is an additional T-qubit clock register. The Kitaev clock states |t^clock,t{0,,T} are defined as |t^=|1t0Tt (the definition of clock states will differ in our construction). This reduction has the following properties:

Completeness

If xLyes, then the history state

|hist=1T+1t=0TUtU1|init|t^ (2)

is a zero-energy state of H.

Soundness

If xLno, then any state |φ has at least 1poly(n) energy with respect to H.

Instance of Quantum-Sat

H is a sum of projectors. Moreover, if each gate in the chosen gate is k-local, then H is (k+3)-local.

Each term in the Hamiltonian represents a constraint on the proof state (which is claimed to be the history state); an energy penalty is given if the constraint is not met.

4.3 𝑯𝐩𝐫𝐨𝐩 and 𝑯𝐟𝐨𝐫𝐦𝐚𝐭

Our construction differs from the Kitaev construction in the definition of Hprop and Hformat, which we define in this section. First, we make a simple observation that the verification circuit Vx can be assumed to be in a standard form.

Lemma 4.3.

Without loss of generality, Vx=UTU1 (where Ui𝒢XZ) satisfies the following properties: (i) T is odd, (ii) Ui is diagonal in the X basis when i is odd, and (iii) Ui is diagonal in the Z basis when i is even.

Proof.

Let WτW1 be any 𝖰𝖬𝖠1 verification circuit with gates Wj from the gate set 𝒢XZ. Recall that by ˜4.2, each Wj is diagonal in either the X basis or the Z basis. By padding the circuit with identity operations, we can construct an equivalent circuit UTU1 of T=2τ+1 gates which has the desired form. Concretely, for each j[τ], if Wj is diagonal in the X basis, set U2j1=Wj, and otherwise set U2j=Wj. The gates Ui which remain undefined by this procedure are set to be the identity operation.

We now define the clock states in our construction, which differ from those for the original Kitaev Hamiltonian. Let

|0^ =|0+0+····0,|1^=|1+0+····0,|2^=|1-0+····0,,|T^=|1-1-····1.

To remove ambiguity, we call these states good clock states, and refer to the clock states of the Kitaev Hamiltonian as Kitaev clock states. Unlike a Kitaev clock state, which has Z basis states on all qubits, a good clock state consists of alternating Z and X basis states. On the i-th “tick” of the clock, the |0 in the i-th position changes to a |1 if i is odd and changes from |+ to |- if i is even.

Our new propagation term has the form Hprop=t=1THprop,t, with

Hprop,t={12[I(|0+0+|1,2+|1+1+|1,2)U1|1+0+|1,2U1|0+1+|]if t=1,12[I(|1+01+0|t1,t,t+1+|1-01-0|t1,t,t+1)Ut|1-01+0|t1,t,t+1Ut|1+01-0|t1,t,t+1]if t is even,12[I(|-0+-0+|t1,t,t+1+|-1+-1+|t1,t,t+1)Ut|-1+-0+|t1,t,t+1Ut|-0+-1+|]if t is odd, t1tT,12[I(|-0-0|T1,T+|-1-1|T1,T)UT|-1-0|T1,TUT|-0-1|T1,T]if t=T.

Importantly, when the gates Ut come from 𝒢XZ, Hprop,t satisfies the conditions imposed on instances of XZ-Quantum-6-Sat.

Lemma 4.4.

For gates Ut𝒢XZ, Hprop,t is a 6-local projector which is diagonal in either the Z basis or the X basis.

Proof.

It is straightforward to check that Hprop,t is a projector, i.e. Hprop,t2=Hprop,t. First, consider when t is odd and t1,T. By ˜4.2, Ut=Ut, so by factoring out qubits t1 and t+1 in the clock register, we have

Hprop,t =|-+-+|t1,t+112[I(|00|t+|11|t)Ut(|10|t+|01|t)]
=|-+-+|t1,t+112[IItUtXt],

where It and Xt are the identity and Pauli X operators on clock qubit t. By ˜4.3, Ut is diagonal in the X basis, since t is odd. In this form, it is clear what basis diagonalizes Hprop,t: the change-of-basis matrix applies H^ on all qubits. Thus, Hprop,t is diagonal in the X basis. Moreorer, by ˜4.2, Ut is 3-local, so Hprop,t is 6-local.

By a similar argument, Hprop,t is diagonal in the Z basis for odd t (including t=1 and t=T), noting that |-+|t+|-+|t=Zt. We now define the formatting Hamiltonian term Hformat=Hformat,X+Hformat,Z, where

Hformat,X =I(|+-+-|2,4+|+-+-|4,6++|+-+-|T3,T1),
Hformat,Z =I(|0101|1,3+|0101|3,5++|0101|T2,T).

The Z format terms give an energy penalty unless the odd clock qubits form a Kitaev clock state. Likewise, the X format terms only allow Kitaev clock states on the even clock qubits (under the relabeling |0|+ and |1|-).

It is clear that every good clock state is a zero-energy state of Hformat. However, note that that there are some zero-energy states of Hformat which are not good clock states! We refer to these states as fake clock states. Intuitively, these are states where the X and Z clocks are each valid Kitaev clock states, but they are not properly synchronized with each other. In the following section, we will show the ground space of H has no support on fake clock states, even though they are zero-energy states of Hformat. We also use the term bad clock states to refer to non-zero energy states of Hformat.

4.4 𝑯𝐢𝐧 and 𝑯𝐨𝐮𝐭

The terms Hin and Hout of the Kitaev construction remain unchanged in our construction, and we state them here:

Hin=iQanc|11|i|00|1,Hout=|00|1|11|T. (3)

Within the subspace of clock states, the local check |00|1 projects onto the subspace spanned by |0^. Thus, Hin assigns an energy penalty to states whose |0^ clock component contains ancilla qubits which are not initialized to |0. Likewise, Hout assigns an energy penalty when the |T^ clock component of the state does not have a |1 on the first qubit of (recall that this is qubit measured at the end of the 𝖰𝖬𝖠1 verification). Note that hist|Hin|hist=hist|Hout|hist=0.

5 XZ-Quantum-6-Sat is QMA1-hard

In this section, we show that XZ-Quantum-6-Sat is 𝖰𝖬𝖠1𝒢XZ-hard, which suffices to prove ˜1.1 by ˜4.1. Completeness is immediate: it is straightforward to check that in the “yes” case, the history state given by Equation 2 (where |t^ now refers to a good clock state, not a Kitaev clock state, and |ψ is the accepting proof of the 𝖰𝖬𝖠1 verifier) is a zero-energy state of our Hamiltonian H.

In the remainder of this section, we prove soundness: in the “no” case, any state in clock has energy at least 1/poly(n) with respect to H, i.e. the smallest eigenvalue λmin(H) is at least 1/poly(n). Recall that T is odd by ˜4.3, and let T=2τ+1.

Setup and notation

We first partition clock into subspaces corresponding to good, fake, and bad clock states. For strings a=a1aj+1 and b=b1bj, we use ab=a1b1a2b2ajbjaj+1 to denote the string formed by interleaving a and b. Define the set of strings S={ab:a{0,1}τ+1,b{+,-}τ}. The states {|s:sS} form a basis for clock. We partition S into subsets Sgood, Sfake, and Sbad. First, define

Sformat={1tZ0τ+1tZ-tX+τtX:tZ{0,,τ+1},tX{0,,τ}}.

Let SgoodSformat contain the strings of the above form which satisfy either tZ=tX or tZ=tX+1. Let Sfake=SformatSgood and Sbad=SSformat. It is clear that Sgood and Sfake correspond to good and fake clock states respectively, while Sbad corresponds to non-zero-energy states of Hformat. Let good be the subspace of clock spanned by {|s:sSgood}={|0^,|1^,,|T^}. Note that good (i.e. the orthogonal complement of good in clock) is the span of {|s:sSfakeSbad}.

Factoring out the good subspace

We show that good is an invariant subspace of H. This is clear for Hin, Hout, and Hformat since each term acts as identity or zero on clock. Then, for t[T] and |φ,

Hprop,t|φ|t1^ =|φ|t1^+Ut|φ|t^, (4)
Hprop,t|φ|t^ =Ut|φ|t1^+|φ|t^, (5)
Hprop,t|φ|k^ =0 for k{0,,T}{t1,t}, (6)

as desired. We can now write

clock=(goodgood)=(good)(good). (7)

In fact, good and good are orthogonal, by the orthogonality of good and good. It follows from Equation 7 that good=(good) (the orthogonal complement of good in clock).

We now apply the following linear algebra fact:

Lemma 5.1.

If W be a O-invariant subspace of Hilbert space V, then W is O-invariant.

Proof.

Let zW. Denoting the inner product by (,), we have for any wW that (w,Oz)=(Ow,z)=0, since OwW and zW. Thus OzW, as desired. We conclude that good is H-invariant as well, using ˜5.1 and the fact that H is Hermitian. H can then be block-diagonalized, simplifying our analysis. We get the eigenvalues of Hclock (i.e. H as an operator on clock) by taking the union of the eigenvalues of Hgood and those of Hgood. In particular,

λmin(Hclock)=min(λmin(Hgood),λmin(Hgood)). (8)

We first lower bound λmin(Hgood). Observe that Hformat has zero energy on good states, and the remaining part Hin+Hout+Hprop is identical to the Hamiltonian in Kitaev’s original reduction [15] after doing a change of basis which applies H^ on all even clock qubits. Since the eigenvalues are invariant under change of basis, the lower bound from the analysis in [15] carries over directly: we have that

λmin(Hgood)1/poly(n). (9)

In the next part of this section, we find a lower bound when H acts on good.

Bounding the energy outside the good subspace

First, note that for any |φgood,

φ|(Hin+Hout+Hprop+Hformat)|φφ|(Hprop+Hformat)|φ,

since Hin+Hout is positive definite. It thus suffices to argue that Hprop+Hformat has energy at least 1/poly(n) on good (in fact, we will show that the energy is Ω(1)).

Define fake and bad in terms of Sfake and Sbad analogously to good, so that good=fakebad. Suppose we have a state |ψ supported only on good. We write this state as

|ψ=i|ψi|iclock, (10)

where the vectors |i are clock basis states, and the vectors |ψi are subnormalized. We are now going to bound Hprop+Hformat by decomposing it in terms of the components |ψi. Firstly, Hformat is simple to bound:

ψ|Hformat|ψiSbad|ψi2. (11)

This follows because Hformat assigns an energy penalty of at least 1 to every bad string. Next, let us look at a single propagation term Hprop,t. We will assume that t is odd and 1<t<T for simplicity, and also use the fact that all of our gates are self-adjoint to simplify expressions. We will write this propagation term as a sum of terms that look like graph Laplacians over the graph whose vertices consist of clock basis strings, and whose edges correspond to pairs of strings that are paired by Hprop,t.

ψ|Hprop,t|ψ =ψ|[12(IIt(Ut)Xt)|-+-+|t1,t+1Irestofclock]|ψ (12)
=12i,jψi|i|clock[(IIt(Ut)Xt)|-+-+|t1,t+1
Irestofclock]|ψj|jclock (13)
=12(i,j)Et(ψi|ψi+ψj|ψjψi|Ut|ψjψj|Ut|ψi) (14)
=12(i,j)Et|ψiUt|ψj2, (15)

where Et consists of all unordered pairs (i,j) of strings in SfakeSbad such that i and j agree at all positions except position t, disagree at position t, and both strings have - in position t1 and + in position t+1. (For other values of t other than the ones we have considered here, Et may be defined analogously: we give a full definition encapsulating all edge cases in ˜6.1.) To simplify notation for what will follow, let us replace Ut in the last line with Uij: this is well-defined because it is easy to see that the sets Et are disjoint for different values of t. Thus, we have

ψ|Hprop,t|ψ=12(i,j)Et|ψiUij|ψj2. (16)

Putting these together, the total energy of Hprop+Hformat is

ψ|(Hprop+Hformat)|ψ(i,j)E12|ψiUij|ψj2+iSbad|ψi2, (17)

where E=tEt is the union of all the edges associated with all the terms of Hprop.

To lower-bound this quantity, we need to show that the fake strings are well-connected to the bad strings, thus forcing a large amount of weight onto the second summation. To do this, we will use two facts shown in the next section to study the connected component structure of the graph G whose vertices are SfakeSbad and whose edges are E. By ˜6.2, every fake string is connected to at least one bad string. This means that the connected components of G consist either exclusively bad strings, or a mixture of fake and bad strings: there are no components of G containing only fake strings. Furthermore, by ˜6.4, each connected component contains at most one fake string. Thus, the components either consist entirely of bad strings, or exactly one fake string and at least one bad string.

Moreover, we can freely drop edges from the summation in Equation 17, and only lower the energy. We will try dropping all the “bad to bad” edges, so that the only remaining edges are between fake and bad strings. Let us introduce the notation N(i) to refer to the set of all neighbors in the graph G of a string i. Observe that by the facts about the component structure of G mentioned in the previous paragraph, it holds that for distinct i,jSfake, the sets N(i),N(j) are disjoint, so every string xSbad is contained in at most one N(i) for iSfake. Using this characterization, we write Equation 17 after the bad-to-bad edges have been omitted as follows:

ψ|(Hprop,t+Hformat)|ψ iSfake[jN(i)(12|ψiUij|ψj2+|ψj2)]
+iSbadleftover|ψi2, (18)

where Sbadleftover=SbadiSfakeN(i).

Now, because the neighborhoods N(i) are disjoint for distinct iSfake, we see that the entire RHS of Equation 18 can be written as a sum

iSfakeSbadleftoverψ|Mi|ψ,

where the matrices Mi are Hermitian matrices that act on different factors of the clock space, so that

ψ|Mi|ψ={jN(i)(12|ψiUij|ψj2+|ψj2)if iSfake,|ψi2if iSbadleftover. (19)

From this, it is easy to see that the minimum value of the RHS of Equation 18 is simply the minimum eigenvalue of the matrices Mi, restricted to their associated factors. For iSbadleftover, Mi acts as the identity matrix on its corresponding clock sector, so its minimum eigenvalue is 1. It thus remains to lower bound the minimum eigenvalue of Mi when iSfake. By ˜6.6, this is at least 1/4. So overall, we conclude that for |ψ supported on good,

ψ|(Hprop+Hformat)|ψ1/4,

and hence

λmin(Hgood)1/4. (20)

Combining Equation 8 with Equation 9 and Equation 20, we conclude that

λmin(H)1/poly(n)

as desired, establishing the soundness property.

6 Soundness: technical lemmas

In this section we prove the technical lemmas that were used in the soundness proof in the previous section.

The structure of 𝑯𝐩𝐫𝐨𝐩

We start with some useful facts about how strings in SfakeSbad are coupled together by Hprop. Our first step is to more formally define the graph G capturing this coupling, which was introduced in the previous section.

Definition 6.1.

Define a graph G with vertex set SfakeSbad and edge set E=i=1TEi, where

  • (u,v)E1 iff u1v1, u2=v2=+, and ui=vi for all i3.

  • For t{2,4,,T3,T1}, (u,v)Et iff utvt, ut1=vt1=1, ut+1=vt+1=0, and ui=vi for all it1,t,t+1.

  • For t{3,5,,T4,T2}, (u,v)Et iff utvt, ut1=vt1=-, ut+1=vt+1=+, and ui=vi for al it1,t,t+1.

  • (u,v)ET iff uTvT, uT1=vT1=-, and ui=vi for all iT2.

For vertex v in G, let 𝒞(v) be the connected component of G which contains v.

It is clear by the definition of Hprop that edges in G correspond to states which are connected under the action of Hprop (in the sense of Equation 4). We now prove important lemmas about the connectivity of G (˜6.2, ˜6.4), which allow us to understand Hprop.

Lemma 6.2.

For every fSfake, there is some bSbad such that (f,b)E.

Proof.

Let fSfake. Then f=1tZ0τ+1tZ-tX+τtX for some tZ{0,,τ+1}, tX{0,,τ} satisfying tZtX and tZtX+1. We case on tX. First, suppose tX=0. Then tZ2, so f starts with 1+1. Define b to be the same string as f, but with the first position changed to a 0. Then (f,b)E1, by the definition of E1. Note that there is a 0 followed by a 1 on the odd terms of b, which puts bSbad, as desired.

Suppose instead that tX[τ1]. We can subdivide into two cases: tZ satisfies either tZtX1 or tZtX+2. If tZtX1, then f looks like

* * 0 0 0 0
····           ····,
    - - -¯ + +

where we have have underlined position 2tX and offset the odd and even terms for clarity, with * denoting either a 0 or a 1. Let b be the same as f but with position 2tX+1 changed to a 1. Then (f,b)E2tX+1, and bSbad since b2tX1=0 while b2tX+1=1.

If instead tZtX+2, then f looks like

1 1 1 1 * *
····           ····,
    - -¯ + + +

where we have underlined position 2tX. Define b to be f but with position 2tX+1 changed to a 0. Once again, we have (f,b)E2tX+1, and bSbad since b2tX+1=0 while b2tX+3=1.

Finally, suppose tX=τ. Then tZτ1, so f ends in 0-0. Let b be the same as f but with a 1 in the final position. Then (f,b)ET and bSbad, completing the proof. We now prove an intermediate “locking” lemma.

Lemma 6.3.

Let sSfakeSbad.

  1. 1.

    If s contains +1 as a substring (say, in positions k and k+1), then any string in 𝒞(s) contains +1 in positions k and k+1.

  2. 2.

    If s contains 0- in positions k and k+1, then any string in 𝒞(s) contains 0- in positions k and k+1.

Proof.

It suffices to consider the four (two-way) rewrite rules

[0+[1+1+01-0-0+-1+-0]-1]

on string s, where [ marks the beginning of a string and ] marks the end of a string. We first show (1). By the second rewrite rule, position k+1 must be 0 in order for position k to change under a rewrite starting from s. Likewise, by the remaining rewrite rules, position k must be - in order for position k+1 to change under a rewrite. Together, this implies that positions k and k+1 are invariant under rewrites from s, as desired. The analysis is analogous to show (2).

Lemma 6.4.

For any fSfake, f is the only element of Sfake which is in 𝒞(f).

Proof.

Fix distinct f,gSfake. We want to show that f and g are in different components of G. We case on the the first position k[T] at which fkgk. Suppose k=1 and without loss of generality take f1=0 and g1=1. Since fSfake, we must have fi=0 for all odd i.

We now argue that f2=-. If instead f2=+, then fi=+ for all even i, since fSfake. But then f=0τ+1+τSgood, which is a contradiction.

Applying ˜6.3 to f1f2=0-, every string h𝒞(f) satisfies h1h2=0-. In particular, since g1=1, g1𝒞(f), as desired. We depict the above deductions as

f =0-¯0*0····
g =1****····,

where the underlined positions are those “locked” by ˜6.3, and * indicates a position which is unconstrained (or otherwise unimportant in our argument).

Now suppose k=T. Without loss of generality, take fT=0 and gT=1. We show that

f =····****0
g =····1*1+1¯.

Since gSfake, gi=1 for all odd i. Now, note that gT1=+: if instead gT1=-, then gi=- for all even i, which puts g=1τ+1-τSgood, a contradiction. Then, applying ˜6.3 to gT1gT=+1, every string in 𝒞(g) ends in +1, and thus f𝒞(g) since fT=0.

In the third case, suppose k is odd and 1<k<T. As usual, take fk=0 and gk=1. We consider two subcases corresponding to the possible values of gk1. The case gk1=+ is immediate: ˜6.3 locks positions k1 and k of any string in 𝒞(g) to +1, so f𝒞(g). In the other case gk1=-, we will deduce that

f =1-1-····1-0ˇ-¯0*0····
g =1-1-····1-1****····,

where we have accented the k-th position of f. Since gSfake, gi=1 for all odd ik, and gi=- for all even ik1. By definition, k is the first position in which f and g differ, so fi=1 for all odd i<k, and fi=- for all even ik1. Moreover, since fSfake, fi=0 for all odd ik. Together, this implies that fk+1=-; if instead fk+1=+, then fi=+ for all even ik+1 and thus fSgood, which is a contradiction. ˜6.3 then locks positions k and k+1 of any string in 𝒞(f) to 0-, so g𝒞(f), as desired.

Finally, suppose k is even and take fk=+, gk=-. This case is analogous to the previous one. First, if gk1=0 then ˜6.3 locks positions k1 and k of any string in 𝒞(g) to 0-, so f𝒞(g). If instead gk1=1, then it is easy to check that f and g have the form

f =1-1-····1+ˇ1¯+*+···
g =1-1-····1-****···

using the reasoning of the previous case. Thus g𝒞(f), completing the proof.

Optimizing the energy over one component

We now bound the minimum of the quadratic form appearing in Equation 18 for a component containing a fake string, and show that it is lower-bounded by a constant independent of the number of bad strings in the component. We first show the desired bound in a special case where the computational Hilbert space is one-dimensional, from which the general case will follow as an easy corollary.

Lemma 6.5.

Let k1 be an integer and let ψ=(ψ0,,ψk)k+1 be a unit vector. Then

f(ψ):=i=1k(12|ψ0ψi|2+|ψi|2)1/4. (21)
Proof.

We would like to minimize

f(ψ)=i=1k(12|ψ0ψi|2+|ψi|2) (22)

subject to ψ being a unit vector in k+1. Our first observation is to see that we can without loss of generality take all coordinates of ψ to be real and nonnegative. Next, observe that

i=1k|ψ0ψi|2 =kψ02+(1ψ02)2ψ0i=1kψi (23)
kψ02+(1ψ02)2ψ0ki=1kψi2 (24)
=(kψ01ψ02)2, (25)

where we have used the Cauchy-Schwarz inequality in passing to (24), and then applied the normalization condition. So overall the quantity we want to minimize is

f(ψ)g(ψ0) :=12(kψ01ψ02)2+(1ψ02) (26)
=(ψ01ψ02)(k/2k/2k/23/2)G(ψ01ψ02). (27)

Thus, we have reduced the problem to finding the minimum eigenvalue of the 2×2 matrix G. Through explicit computation, we can find the eigenvalues of this matrix.

λ =3+k4(1±18k/(3+k)2). (28)

It is not hard to show that for all integer k1, λ1/4 (for details, see the full version).

Corollary 6.6.

Let k,d1 be integers and let |α0,,|αk be (not necessarily unit) vectors over d such that i=0k|αi2=1. Moreover, let U1,,Uk be d×d unitary matrices. Then

i=1k(12|α0Ui|αi2+|αi2)1/4. (29)
Proof.

Observe that for any two vectors |α,|β and for any unitary U, it holds that

|αU|β2 =|α2+|β2α|U|ββ|U|α (30)
|α2+|β22|α|β (31)
=||α|β|2. (32)

Thus, if we let ψi=|αi, then the vector (ψ0,,ψk) satisfies the conditions of ˜6.5, and by Equation 32, the quantity we wish to bound is at least f(ψ)1/4 by ˜6.5.

References

  • [1] Dorit Aharonov. A Simple Proof that Toffoli and Hadamard are Quantum Universal, 2003. doi:10.48550/arXiv.quant-ph/0301040.
  • [2] Dorit Aharonov, Itai Arad, and Thomas Vidick. Guest column: The quantum PCP conjecture. SIGACT News, 44(2):47–79, 2013. doi:10.1145/2491533.2491549.
  • [3] Eric R. Anschuetz, David Gamarnik, and Bobak Kiani. Combinatorial NLTS from the overlap gap property. Quantum, 8:1527, 2024. doi:10.22331/Q-2024-11-19-1527.
  • [4] Anurag Anshu, Nikolas P Breuckmann, and Quynh T Nguyen. Circuit-to-hamiltonian from tensor networks and fault tolerance. In Proceedings of the 56th Annual ACM Symposium on Theory of Computing, pages 585–595, 2024. doi:10.1145/3618260.3649690.
  • [5] Anurag Anshu, Nikolas P. Breuckmann, and Chinmay Nirkhe. NLTS Hamiltonians from good quantum codes. In Proceedings of the 55th Annual ACM Symposium on Theory of Computing, STOC 2023, pages 1090–1096, New York, NY, USA, 2023. Association for Computing Machinery. doi:10.1145/3564246.3585114.
  • [6] Jacob D. Biamonte and Peter J. Love. Realizable Hamiltonians for universal adiabatic quantum computers. Physical Review A, 78(1):012352, 2008. doi:10.1103/PhysRevA.78.012352.
  • [7] Sergey Bravyi. Efficient algorithm for a quantum analogue of 2-SAT, 2006. doi:10.48550/arXiv.quant-ph/0602108.
  • [8] Sergey Bravyi and Barbara Terhal. Complexity of stoquastic frustration-free hamiltonians. Siam journal on computing, 39(4):1462–1485, 2010. doi:10.1137/08072689X.
  • [9] Stephen A. Cook. The complexity of theorem-proving procedures. In Proceedings of the Third Annual ACM Symposium on Theory of Computing, STOC ’71, pages 151–158, New York, NY, USA, 1971. Association for Computing Machinery. doi:10.1145/800157.805047.
  • [10] Toby Cubitt and Ashley Montanaro. Complexity Classification of Local Hamiltonian Problems. SIAM Journal on Computing, 45(2):268–316, 2016. doi:10.1137/140998287.
  • [11] Christopher M. Dawson and Michael A. Nielsen. The Solovay-Kitaev algorithm. Quantum Info. Comput., 6(1):81–95, 2006. doi:10.26421/QIC6.1-6.
  • [12] Lior Eldar and Aram W. Harrow. Local Hamiltonians Whose Ground States Are Hard to Approximate. In 2017 IEEE 58th Annual Symposium on Foundations of Computer Science (FOCS), pages 427–438, 2017. doi:10.1109/FOCS.2017.46.
  • [13] David Gosset and Daniel Nagaj. Quantum 3-SAT Is QMA1-Complete. SIAM Journal on Computing, 45(3):1080–1128, 2016. doi:10.1137/140957056.
  • [14] Robbie King and Tamara Kohler. Gapped Clique Homology on Weighted Graphs is QMA1-Hard and Contained in QMA. In 2024 IEEE 65th Annual Symposium on Foundations of Computer Science (FOCS), pages 493–504, 2024. doi:10.1109/FOCS61266.2024.00039.
  • [15] A. Yu. Kitaev, A. Shen, and M. N. Vyalyi. Classical and Quantum Computation. American Mathematical Society, USA, June 2002.
  • [16] L. A. Levin. Universal sequential search problems. Problems of Information Transmission, 9(3):265–266, 1973.
  • [17] Michael A. Nielsen and Isaac L. Chuang. Quantum Computation and Quantum Information: 10th Anniversary Edition, 2010. doi:10.1017/CBO9780511976667.
  • [18] Dorian Rudolph. Towards a universal gateset for QMA1, 2025. doi:10.48550/arXiv.2411.02681.
  • [19] Yaoyun Shi. Both Toffoli and controlled-NOT need little help to do universal quantum computing. Quantum Info. Comput., 3(1):84–92, 2003. doi:10.26421/QIC3.1-7.
  • [20] Adam Wills, Ting-Chun Lin, and Min-Hsiu Hsieh. Local testability of distance-balanced quantum codes. npj Quantum Information, 10(1):120, 2024. doi:10.1038/s41534-024-00908-8.