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 X⁢Z-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 X⁢Z-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,C⁢X,C⁢C⁢X) 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 X⁢Z-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+1⁢∑t=0TUt⁢…⁢U1⁢|ψ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^⟩=|1⟩1⁢…⁢|1⟩t−1⁢|1⟩t⁢|0⟩t+1⁢…⁢|0⟩T.

    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 |t−1^⟩ 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 =12⁢I⊗(|110⟩⁢⟨110|+|100⟩⁢⟨100|)t−1,t,t+1
    −12⁢Ut⊗(|110⟩⁢⟨100|)t−1,t,t+1−12⁢Ut†⊗(|100⟩⁢⟨110|)t−1,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⁢(I⊗It−Ut⊗Xt)⊗(|10⟩⁢⟨10|)t−1,t+1,

which is diagonal when all qubits are placed in the X-basis except for qubits t−1,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⁢(I⊗It−Ut⊗Xt)⊗|-+⟩⁢⟨-+|t−1,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 U⊗I[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 C⁢X, C⁢Z, and SWAP denote the two-qubit CNOT, controlled-Z and swap operations respectively. Let C⁢C⁢X denote the Toffoli gate and C⁢C⁢Z 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=Gm⁢⋯⁢G1, 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 {|0⟩⟨0|1⊗I,|1⟩⟨1|1⊗I}, 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 Lyes∩Lno=∅. Promise problems generalize languages, a notion from classical complexity theory: a language is just a promise problem which additionally satisfies Lyes∪Lno={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

‖(|1⟩⁢⟨1|1⊗I)⁢Vx⁢|ψ⟩⁢|0⁢⋯⁢0⟩anc‖2. (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 x∈Lyes, then there is some proof state for which Vx accepts with probability 1.

  2. 2.

    If x∈Lno, 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,C⁢X,C⁢C⁢X,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 C⁢C⁢X 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 x∈Lyes, 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 x∈Kyes, and A⁢(x)∈Lno if x∈Kno. 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 k≥1 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 n≥1, let ℬn be a set of bases of n-qubit space. Let ℬ=∪n≥1ℬn. 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 Γ⁢(C⁢X) are already included in 𝒢2. Let a denote an ancilla qubit. Γ⁢(C⁢C⁢X) on control qubits i,j,k and target qubit l can be decomposed as

Γ⁢(C⁢C⁢X)i,j,k,l⊗Ia=C⁢C⁢Xi,j,a⁢C⁢C⁢Xk,a,l⁢C⁢C⁢Xi,j,a⁢C⁢C⁢Xk,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,k⊗Ia=Γ⁢(SWAP)i,j,k⁢(H^⊗H^)k,a⁢Γ⁢(SWAP)i,j,k⁢(H^⊗H^)k,a.

Then, since Γ⁢(A⁢B)=Γ⁢(A)⁢Γ⁢(B) for any gates A and B, we may obtain a decomposition of Γ⁢(SWAP) in terms of C⁢C⁢X’s by applying the identity SWAPi,j=C⁢Xi,j⁢C⁢Xj,i⁢C⁢Xi,j. This completes the proof. We illustrate these constructions in Figure 1. ◀

Figure 1: Decomposition of Γ⁢(C⁢C⁢X) 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 C⁢C⁢X, 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

Π=∑z∈SV⁢|z⟩⁢⟨z|⁢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 z⋆∈S, and otherwise returns outcome I−Π. The probability of getting outcome Π is

∑z∈S‖|z⟩⁢⟨z|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,C⁢X,C⁢C⁢X,H^⊗H^}. Define a new gate set

𝒢X⁢Z={X,C⁢Z,C⁢C⁢Z,G^},

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

G^=(H^⊗H^)⁢C⁢Z⁢(H^⊗H^).

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

Lemma 4.1.

𝖰𝖬𝖠⁢1𝒢2=𝖰𝖬𝖠⁢1𝒢X⁢Z.

Proof.

We start by showing that 𝖰𝖬𝖠⁢1𝒢X⁢Z⊆𝖰𝖬𝖠⁢1𝒢2. It suffices to show that each gate in 𝒢X⁢Z 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 U∈S¯ if U can be exactly written as a sequence of gates from S. The following closure property is immediate: if S′⊆S¯ and U∈S′¯, then U∈S¯. In particular, if V∈S¯ and U∈S∪{V}¯, then U∈S¯.

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

C⁢Zi,j⊗Ia=(H^⊗H^)j,a⁢C⁢Xi,j⁢(H^⊗H^)j,a.

This puts C⁢Z∈𝒢2¯. A similar identity

C⁢C⁢Zi,j,k⊗Ia=(H^⊗H^)k,a⁢C⁢C⁢Xi,j,k⁢(H^⊗H^)k,a

puts C⁢C⁢Z∈𝒢2¯. Finally, G^∈𝒢2∪{C⁢Z}¯ by definition, so by the closure property, G^∈𝒢2¯.

We now show that 𝖰𝖬𝖠⁢1𝒢2⊆𝖰𝖬𝖠⁢1𝒢X⁢Z 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^⋅C⁢Z⋅G^ using the circuit

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

where we used that H^j⁢C⁢Zi,j⁢H^j=C⁢Xi,j and SWAPi,j=C⁢Xi,j⁢C⁢Xj,i⁢C⁢Xi,j. This puts F^∈𝒢X⁢Z¯. Next, using that H^⊗H^ and SWAP commute, we have that

C⁢Xi,j⊗Ia=F^j,a⁢C⁢Zi,a⁢F^j,a,

since

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

Thus C⁢X∈𝒢X⁢Z∪{F^}¯, implying that C⁢X∈𝒢X⁢Z¯ by the closure property. The identity

C⁢C⁢Xi,j,k⊗Ia=F^k,a⁢C⁢C⁢Zi,j,a⁢F^k,a

similarly puts C⁢C⁢X∈𝒢X⁢Z¯. Finally, note that H^⊗H^=SWAP⋅F^. Writing SWAP as three C⁢X’s, we have that H^⊗H^∈{F^,C⁢X}¯, so H^⊗H^∈𝒢X⁢Z¯ by the closure property. ◀

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

Lemma 4.2.

Each gate in 𝒢X⁢Z 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⟩=|ψ⟩⁢|0⁢⋯⁢0⟩anc. 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=UT⁢⋯⁢U1, 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^⟩=|1t⁢0T−t⟩ (the definition of clock states will differ in our construction). This reduction has the following properties:

Completeness

If x∈Lyes, then the history state

|hist⟩=1T+1⁢∑t=0TUt⁢⋯⁢U1⁢|init⟩⊗|t^⟩ (2)

is a zero-energy state of H.

Soundness

If x∈Lno, 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=UT⁢⋯⁢U1 (where Ui∈𝒢X⁢Z) 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 𝒢X⁢Z. 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 UT⁢⋯⁢U1 of T=2⁢τ+1 gates which has the desired form. Concretely, for each j∈[τ], if Wj is diagonal in the X basis, set U2⁢j−1=Wj, and otherwise set U2⁢j=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,2−U1†⊗|0+⟩⟨1+|]if ⁢t=1,12[Iℋ⊗(|1+0⟩⟨1+0|t−1,t,t+1+|1-0⟩⟨1-0|t−1,t,t+1)−Ut⊗|1-0⟩⟨1+0|t−1,t,t+1−Ut†⊗|1+0⟩⟨1-0|t−1,t,t+1]if ⁢t⁢ is even,12[Iℋ⊗(|-0+⟩⟨-0+|t−1,t,t+1+|-1+⟩⟨-1+|t−1,t,t+1)−Ut⊗|-1+⟩⟨-0+|t−1,t,t+1−Ut†⊗|-0+⟩⟨-1+|]if t is odd, t≠1, t≠T,12[Iℋ⊗(|-0⟩⟨-0|T−1,T+|-1⟩⟨-1|T−1,T)−UT⊗|-1⟩⟨-0|T−1,T−UT†⊗|-0⟩⟨-1|T−1,T]if ⁢t=T.

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

Lemma 4.4.

For gates Ut∈𝒢X⁢Z, 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 t≠1,T. By ˜4.2, Ut=Ut†, so by factoring out qubits t−1 and t+1 in the clock register, we have

Hprop,t =|-+⟩⁢⟨-+|t−1,t+1⊗12⁢[Iℋ⊗(|0⟩⁢⟨0|t+|1⟩⁢⟨1|t)−Ut⊗(|1⟩⁢⟨0|t+|0⟩⁢⟨1|t)]
=|-+⟩⁢⟨-+|t−1,t+1⊗12⁢[Iℋ⊗It−Ut⊗Xt],

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+⋯+|+-⟩⁢⟨+-|T−3,T−1),
Hformat,Z =Iℋ⊗(|01⟩⁢⟨01|1,3+|01⟩⁢⟨01|3,5+⋯+|01⟩⁢⟨01|T−2,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=∑i∈Qanc|1⟩⁢⟨1|i⊗|0⟩⁢⟨0|1,Hout=|0⟩⁢⟨0|1⊗|1⟩⁢⟨1|T. (3)

Within the subspace of clock states, the local check |0⟩⁢⟨0|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𝒢X⁢Z-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=a1⁢⋯⁢aj+1 and b=b1⁢⋯⁢bj, we use a⋈b=a1⁢b1⁢a2⁢b2⁢⋯⁢aj⁢bj⁢aj+1 to denote the string formed by interleaving a and b. Define the set of strings S={a⋈b:a∈{0,1}τ+1,b∈{+,-}τ}. The states {|s⟩:s∈S} form a basis for ℋclock. We partition S into subsets Sgood, Sfake, and Sbad. First, define

Sformat={1tZ⁢0τ+1−tZ⋈-tX⁢+τ−tX:tZ∈{0,…,τ+1},tX∈{0,…,τ}}.

Let Sgood⊆Sformat contain the strings of the above form which satisfy either tZ=tX or tZ=tX+1. Let Sfake=Sformat∖Sgood and Sbad=S∖Sformat. 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⟩:s∈Sgood}={|0^⟩,|1^⟩,…,|T^⟩}. Note that ℋgood⟂ (i.e. the orthogonal complement of ℋgood in ℋclock) is the span of {|s⟩:s∈Sfake∪Sbad}.

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⁢|φ⟩⁢|t−1^⟩ =|φ⟩⁢|t−1^⟩+Ut⁢|φ⟩⁢|t^⟩, (4)
Hprop,t⁢|φ⟩⁢|t^⟩ =Ut†⁢|φ⟩⁢|t−1^⟩+|φ⟩⁢|t^⟩, (5)
Hprop,t⁢|φ⟩⁢|k^⟩ =0⁢ for ⁢k∈{0,…,T}∖{t−1,t}, (6)

as desired. We can now write

ℋ⊗ℋclock=ℋ⊗(ℋgood⊕ℋgood⟂)=(ℋ⊗ℋ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 z∈W⟂. Denoting the inner product by (⋅,⋅), we have for any w∈W that (w,O†⁢z)=(O⁢w,z)=0, since O⁢w∈W and z∈W⟂. Thus O†⁢z∈W⟂, 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 H↾ℋ⊗ℋclock (i.e. H as an operator on ℋ⊗ℋclock) by taking the union of the eigenvalues of H↾ℋ⊗ℋgood and those of H↾ℋ⊗ℋgood⟂. In particular,

λmin(H↾ℋ⊗ℋclock)=min(λmin(H↾ℋ⊗ℋgood),λmin(H↾ℋ⊗ℋgood⟂)). (8)

We first lower bound λmin(H↾ℋ⊗ℋgood). 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(H↾ℋ⊗ℋgood)≥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⟂=ℋfake⊕ℋbad. Suppose we have a state |ψ⟩ supported only on ℋgood⟂. We write this state as

|ψ⟩=∑i|ψi⟩⊗|i⟩clock, (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|ψ⟩≥∑i∈Sbad‖|ψi⟩‖2. (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⁢(Iℋ⊗It−(Ut)ℋ⊗Xt)⊗|-+⟩⁢⟨-+|t−1,t+1⊗Irestofclock]⁢|ψ⟩ (12)
=12∑i,j⟨ψi|⊗⟨i|clock[(Iℋ⊗It−(Ut)ℋ⊗Xt)⊗|-+⟩⟨-+|t−1,t+1
⊗Irestofclock]|ψj⟩⊗|j⟩clock (13)
=12⁢∑(i,j)∈Et(⟨ψi|ψi⟩+⟨ψj|ψj⟩−⟨ψi|Ut|ψj⟩−⟨ψj|Ut|ψi⟩) (14)
=12⁢∑(i,j)∈Et‖|ψi⟩−Ut⁢|ψj⟩‖2, (15)

where Et consists of all unordered pairs (i,j) of strings in Sfake∪Sbad such that i and j agree at all positions except position t, disagree at position t, and both strings have - in position t−1 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 Ui⁢j: 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‖|ψi⟩−Ui⁢j⁢|ψj⟩‖2. (16)

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

⟨ψ|(Hprop+Hformat)|ψ⟩≥∑(i,j)∈E12⁢‖|ψi⟩−Ui⁢j⁢|ψj⟩‖2+∑i∈Sbad‖|ψi⟩‖2, (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 Sfake∪Sbad 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,j∈Sfake, the sets N⁢(i),N⁢(j) are disjoint, so every string x∈Sbad is contained in at most one N⁢(i) for i∈Sfake. Using this characterization, we write Equation 17 after the bad-to-bad edges have been omitted as follows:

⟨ψ|(Hprop,t+Hformat)|ψ⟩ ≥∑i∈Sfake[∑j∈N⁢(i)(12⁢‖|ψi⟩−Ui⁢j⁢|ψj⟩‖2+‖|ψj⟩‖2)]
+∑i∈Sbad−leftover‖|ψi⟩‖2, (18)

where Sbad−leftover=Sbad∖⋃i∈SfakeN⁢(i).

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

∑i∈Sfake∪Sbad−leftover⟨ψ|Mi|ψ⟩,

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

⟨ψ|Mi|ψ⟩={∑j∈N⁢(i)(12⁢‖|ψi⟩−Ui⁢j⁢|ψj⟩‖2+‖|ψj⟩‖2)if i∈Sfake,‖|ψi⟩‖2if i∈Sbad−leftover. (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 i∈Sbad−leftover, 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 i∈Sfake. 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(H↾ℋ⊗ℋgood⟂)≥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 Sfake∪Sbad 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 Sfake∪Sbad and edge set E=∪i=1TEi, where

  • ■

    (u,v)∈E1 iff u1≠v1, u2=v2=+, and ui=vi for all i≥3.

  • ■

    For t∈{2,4,⋯,T−3,T−1}, (u,v)∈Et iff ut≠vt, ut−1=vt−1=1, ut+1=vt+1=0, and ui=vi for all i≠t−1,t,t+1.

  • ■

    For t∈{3,5,…,T−4,T−2}, (u,v)∈Et iff ut≠vt, ut−1=vt−1=-, ut+1=vt+1=+, and ui=vi for al i≠t−1,t,t+1.

  • ■

    (u,v)∈ET iff uT≠vT, uT−1=vT−1=-, and ui=vi for all i≤T−2.

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 f∈Sfake, there is some b∈Sbad such that (f,b)∈E.

Proof.

Let f∈Sfake. Then f=1tZ⁢0τ+1−tZ⋈-tX⁢+τ−tX for some tZ∈{0,…,τ+1}, tX∈{0,…,τ} satisfying tZ≠tX and tZ≠tX+1. We case on tX. First, suppose tX=0. Then tZ≥2, 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 b∈Sbad, as desired.

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

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

where we have have underlined position 2⁢tX 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 2⁢tX+1 changed to a 1. Then (f,b)∈E2⁢tX+1, and b∈Sbad since b2⁢tX−1=0 while b2⁢tX+1=1.

If instead tZ≥tX+2, then f looks like

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

where we have underlined position 2⁢tX. Define b to be f but with position 2⁢tX+1 changed to a 0. Once again, we have (f,b)∈E2⁢tX+1, and b∈Sbad since b2⁢tX+1=0 while b2⁢tX+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 b∈Sbad, completing the proof. ◀ We now prove an intermediate “locking” lemma.

Lemma 6.3.

Let s∈Sfake∪Sbad.

  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+0⇔1-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 f∈Sfake, f is the only element of Sfake which is in 𝒞⁢(f).

Proof.

Fix distinct f,g∈Sfake. 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 fk≠gk. Suppose k=1 and without loss of generality take f1=0 and g1=1. Since f∈Sfake, we must have fi=0 for all odd i.

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

Applying ˜6.3 to f1⁢f2=0-, every string h∈𝒞⁢(f) satisfies h1⁢h2=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 g∈Sfake, gi=1 for all odd i. Now, note that gT−1=+: if instead gT−1=-, then gi=- for all even i, which puts g=1τ+1⋈-τ∈Sgood, a contradiction. Then, applying ˜6.3 to gT−1⁢gT=+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 gk−1. The case gk−1=+ is immediate: ˜6.3 locks positions k−1 and k of any string in 𝒞⁢(g) to +1, so f∉𝒞⁢(g). In the other case gk−1=-, 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 g∈Sfake, gi=1 for all odd i≤k, and gi=- for all even i≤k−1. 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 i≤k−1. Moreover, since f∈Sfake, fi=0 for all odd i≥k. Together, this implies that fk+1=-; if instead fk+1=+, then fi=+ for all even i≥k+1 and thus f∈Sgood, 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 gk−1=0 then ˜6.3 locks positions k−1 and k of any string in 𝒞⁢(g) to 0-, so f∉𝒞⁢(g). If instead gk−1=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 k≥1 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⁢ψ0⁢∑i=1kψi (23)
≥k⁢ψ02+(1−ψ02)−2⁢ψ0⁢k⋅∑i=1kψi2 (24)
=(k⁢ψ0−1−ψ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⁢ψ0−1−ψ02)2+(1−ψ02) (26)
=(ψ01−ψ02)⋅(k/2−k/2−k/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±1−8⁢k/(3+k)2). (28)

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

Corollary 6.6.

Let k,d≥1 be integers and let |α0⟩,…,|αk⟩ be (not necessarily unit) vectors over ℂd such that ∑i=0k‖|αi⟩‖2=1. Moreover, let U1,…,Uk be d×d unitary matrices. Then

∑i=1k(12⁢‖|α0⟩−Ui⁢|αi⟩‖2+‖|αi⟩‖2)≥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+‖|β⟩‖2−2⁢‖|α⟩‖⋅‖|β⟩‖ (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.