alphaFeedback

WorksPhysics and Astronomy

Exact amplitudes of sparse quantum circuits in subexponential time, assuming the new #SAT algorithms

AI contributions
  • Anthropic Claude Fable
  • Anthropic Claude Opus
View details
Ancillary data
3 code links
View

Abstract

Classical simulation lower bounds for quantum circuits have been derived from the Exponential Time Hypothesis (ETH). However, in October 2026 OpenAI and Anthropic claimed deterministic subexponential algorithms for #k-SAT, refuting ETH, and neither manuscript mentions quantum circuits. We record the corollary. Assuming three internal results of the Anthropic manuscript, one deterministic Turing machine computes the amplitude ⟨y∣C∣x⟩\langle y|C|x\rangle of a Clifford+T circuit with GG two-qubit gates in arbitrary connectivity exactly, as an element of Z[ω]/2 k\mathbb{Z}[\omega]/\sqrt2^{\,k} with ω=eiπ/4\omega=e^{i\pi/4}, in time 2O(G/log⁡G)2^{O(G/\sqrt{\log G})} whenever the input length LL satisfies log⁡L≤G0.9\log L\le G^{0.9}. The proof reads the sum over paths, one bit per wire segment, as an Anthropic-type network, with a Kronecker map for Z[ω]\mathbb{Z}[\omega] and a Galois bound. We machine-check this theorem in Lean 4 (amplitude_TM_rate', amplitude_TM_rate, amplitude_TM; axioms propext, Classical.choice, Quot.sound) on top of Anthropic's Lean development (sources at https://github.com/xennygrimmato/eth-hamiltonicity). On paper only, the method gives exact marginals and exact sampling, depth-dd circuits in 2O(dn/log⁡n)2^{O(dn/\sqrt{\log n})}, and, through OpenAI's tensor-contraction theorem, polynomial space when the number of Hadamard-type gates is O(G)O(G). Conversely, constant-depth circuits with O(n)O(n) gates encode #3-SAT with O(n)O(n) clauses, so the rate matches that of the refutations. We check the fine-grained supremacy conjectures we found against their sources; a counting form of the Calabro-Impagliazzo-Paturi width reduction refutes the gap conjecture of Morimae and Tamaki's depth-scaling paper. Dense statements (cubic phases, permanents, random circuits) survive. Nothing here concerns approximate sampling, the constants are galactic, and everything is conditional on the claimed refutations, which are under review.

Introduction

Part of the case for “quantum computational supremacy” rests on fine-grained hypotheses. Morimae and Tamaki [20, 19, 21], Huang, Newman and Szegedy [14, 15] and Dalzell, Harrow, Koh and La Placa [12] derive lower bounds on the classical simulation of quantum circuits from ETH, SETH and relatives. They usually do so by writing #SAT\#\mathrm{SAT}, or the gap gap⁡(f)=∑x(−1)f(x)\operatorname{gap}(f)=\sum_x(-1)^{f(x)} of a Boolean function, as a circuit amplitude. The Exponential Time Hypothesis (ETH) of Impagliazzo and Paturi [16] asserts that there is δ>0\delta>0 such that 3-SAT on nn variables cannot be solved in time O(2δn)O(2^{\delta n}); by the sparsification lemma [17] it is equivalent to the same statement for the number mm of clauses.

However, in October 2026 two manuscripts appeared, by OpenAI [23] and by Anthropic [3], each claiming a deterministic algorithm that counts the satisfying assignments of a kk-CNF in subexponential time, and hence a refutation of ETH. Both come with Lean 4 formalisations of their top-level statements. These are claimed results, still under review by the community. We write [A] and [O] for the two manuscripts. The Main Theorem of [A] counts the models of a kk-CNF with mm clauses on nn variables deterministically in time 2O(nkm/log⁡m, k2nlog⁡n/log⁡n)2^{O(n\sqrt{km/\log m},\,k^2 n\log n/\sqrt{\log n})}; the main theorem of [O] does so in (n+m)O(k)2Ok(nlog⁡log⁡n/log⁡n)(n+m)^{O(k)}2^{O_k(n\log\log n/\sqrt{\log n})} bit operations and polynomial space, and in 2Ok,D(n/log⁡n)2^{O_{k,D}(n/\sqrt{\log n})} when m≤Dnm\le Dn. Neither manuscript mentions quantum circuits (§ 8). So we ask the converse question: what do the new algorithms give for the classical simulation of quantum circuits, and which of the lower bounds above survive? For depth d=O(1)d=O(1) in arbitrary connectivity we know of no previous 2o(n)2^{o(n)} exact method: the interaction graph of such a circuit can be a bounded-degree expander, whose treewidth is Θ(n)\Theta(n).

The problem and a running example

We consider the exact amplitude problem: given a circuit CC and basis states x,yx,y, output ⟨y∣C∣x⟩\langle y|C|x\rangle exactly, as an element of a fixed number ring. We fix the setting as follows. A circuit CC acts on nn qubits and is a sequence of gates from the set

G={H,T,T†,S,S†,X,Z}∪{CNOT,CZ,CS,SWAP,CH}.\mathcal{G}=\{H,T,T^\dagger,S,S^\dagger,X,Z\}\cup\{\mathrm{CNOT},\mathrm{CZ},\mathrm{CS},\mathrm{SWAP},\mathrm{CH}\}.

It has GG two-qubit gates, in arbitrary connectivity, and G1G_1 single-qubit gates. Put ω=eiπ/4\omega=e^{i\pi/4}, so that Z[ω]=Z[x]/(x4+1)\mathbb{Z}[\omega]=\mathbb{Z}[x]/(x^4+1), i=ω2i=\omega^2 and 2=ω−ω3∈Z[ω]\sqrt{2}=\omega-\omega^3\in\mathbb{Z}[\omega]. Every gate U∈GU\in\mathcal{G} is stored as a pair (NU,κU)(N_U,\kappa_U) with NUN_U a matrix over Z[ω]\mathbb{Z}[\omega] and U=NU/2κUU=N_U/\sqrt{2}^{\kappa_U}, where κU=1\kappa_U=1 for HH and CH\mathrm{CH} and κU=0\kappa_U=0 otherwise (for CH\mathrm{CH}, N=diag⁡(2,2)⊕(111−1)N=\operatorname{diag}(\sqrt{2},\sqrt{2})\oplus\left(\begin{smallmatrix}1&1\\1&-1\end{smallmatrix}\right)). Put k(C):=∑UκUk(C):=\sum_U\kappa_U, the number of HH and CH\mathrm{CH} gates. Running times are counted as in [A]: operations on words long enough to hold an address, plus table queries; all logarithms are to base 2. The input has length poly⁡(n+G+G1)\operatorname{poly}(n+G+G_1).

Example 1 (Running example). Let CC act on n=2n=2 qubits: first HH on qubit 1, then CNOT with control 1 and target 2, then CZ. So G=2G=2, G1=1G_1=1 and k(C)=1k(C)=1. The state before CZ is (∣00⟩+∣11⟩)/2(|00\rangle+|11\rangle)/\sqrt{2}, so ⟨11∣C∣00⟩=−1/2\langle11|C|00\rangle=-1/\sqrt{2}, that is, (a0,a1,a2,a3)=(−1,0,0,0)(a_0,a_1,a_2,a_3)=(-1,0,0,0) over 2\sqrt{2}. Each qubit line runs from the first two-qubit gate to the second, so it carries one wire segment between them. The amplitude is therefore a sum over the 222^2 values of the two segment bits; we return to this in Example 2.

The naive routes and the principle

At first one might use a known exact method (Table 1). The Feynman path sum takes time 2O(G)2^{O(G)}. Tensor contraction [18] takes time 2O(tw)2^{O(\mathrm{tw})}, and tw=Θ(n)\mathrm{tw}=\Theta(n) is possible at d=O(1)d=O(1). So neither is 2o(n)2^{o(n)} for shallow circuits on expanders. Instead, we rely on a different principle: a qubit circuit is a sum over binary variables, one per wire segment, coupled by bounded-arity local tables, so it has O(G)O(G) shared bits; that is the kind of input the network engine of [A] consumes, whatever the connectivity. What remains is arithmetic. Amplitudes live in Z[ω]\mathbb{Z}[\omega], which is not a subring of Z\mathbb{Z}, nor of the truncated 2-adic ring of [O]. A Kronecker substitution maps Z[ω]\mathbb{Z}[\omega] to Z/(24b+1)Z\mathbb{Z}/(2^{4b}+1)\mathbb{Z} injectively on a box (Theorem 3.3), and a Galois-conjugate bound shows that the answer lies in the box (Theorem 3.2).

MethodTimeSpaceSource
state vector2npoly⁡2^n\operatorname{poly}2n2^nfolklore
Feynman path sum2O(G)2^{O(G)}poly⁡\operatorname{poly}folklore
recursive (Savitch-type)O(n(2d)n+1)O(n(2d)^{n+1})O(nlog⁡d)O(n\log d)[2 Thm. 4.1]
tensor contraction2O(tw)2^{O(\mathrm{tw})}, tw=Θ(n)\mathrm{tw}=\Theta(n) possible at d=O(1)d=O(1)2O(tw)2^{O(\mathrm{tw})}[18]
geometrically local, dim. ≥2\ge2subexponential in GG–[31]
Theorem 1.12O(G/log⁡G)≤2O(dn/log⁡n)2^{O(G/\sqrt{\log G})} \le2^{O(dn/\sqrt{\log n})}as time under (HA)(\mathrm{H}_{\mathrm{A}})
Theorem 1.2same, if k(C)=O(G)k(C)=O(G)poly under (HO)(\mathrm{H}_{\mathrm{O}})

Table 1. Exact strong simulation of a circuit with nn qubits, GG two-qubit gates and depth dd, arbitrary connectivity. The last two rows are conditional.

What we prove

We use the following hypotheses, which name exactly what we import from the two manuscripts (§ 2).

Hypothesis (H_A). The following statements of [A] hold as stated there: Lemma 8, case (i) (the reduction from general networks to cliques, Theorem 2.3); Theorem A (a),(b) (one algorithm for the (Kℓ,λ)(K_{\ell},\lambda)-cliques of all levels ℓ\ell, Theorem 2.4); and Theorem B (pieces with parameters (1/2,c)(1/\sqrt{2},c), Theorem 2.5).

Hypothesis (H_O). The Boolean tensor-contraction theorem in the appendix of [O] holds as stated there (Theorem 2.6).

Theorem 1.1 (Exact amplitudes; machine-checked). Assume (H_A). There are an absolute constant cc and one deterministic Turing machine that, given a circuit CC on nn qubits with GG two-qubit gates, and x,y∈{0,1}nx,y\in\{0,1\}^n, encoded as a word of length LL with log⁡2(L+2)≤(G+2)0.9\log_2(L+2)\le(G+2)^{0.9}, outputs integers a0,a1,a2,a3a_0,a_1,a_2,a_3 with

⟨y∣C∣x⟩=a0+a1ω+a2ω2+a3ω32k(C).\langle y|C|x\rangle=\frac{a_0+a_1\omega+a_2\omega^2+a_3\omega^3}{\sqrt{2}^{k(C)}}.

within ⌈c⋅2c(G+2)/log⁡2(G+2)⌉\left\lceil c\cdot2^{c(G+2)}/\sqrt{\log_2(G+2)}\right\rceil steps. Without the condition on LL, the same holds within ⌈c⋅2cG/log⁡2(G+2)(L+2)c(G+2)1/15⌉\left\lceil c\cdot2^{cG}/\sqrt{\log_2(G+2)}(L+2)^{c(G+2)^{1/15}}\right\rceil steps, and for every δ>0\delta>0 there is a machine running within ⌈cδ2δG(L+2)cδ⌉\left\lceil c_\delta2^{\delta G}(L+2)^{c_\delta}\right\rceil steps.

Since the input word has length L=3n+1+(G+G1)(2n+4)L=3n+1+(G+G_1)(2n+4), we have L+2≤6n(G+G1+1)L+2\le6n(G+G_1+1) for n≥1n\ge1, so the condition on LL holds, for example, whenever n≥1n\ge1 and 6n(G+G1+1)≤2G0.96n(G+G_1+1)\le2^{G^{0.9}}; on paper, the proof below also gives 2cG/log⁡Gpoly⁡(n+G1)2^{cG}/\sqrt{\log G}\operatorname{poly}(n+G_1) under the weaker condition log⁡(G1+2)≤G0.9\log(G_1+2)\le G^{0.9}, but that variant is not machine-checked and not claimed.

Theorem 1.2 (Polynomial space; proved on paper only). Assume (H_O). Put s:=max⁡{4G,k(C),4}s:=\max\{4G,k(C),4\}. The output of Theorem 1.1 can be computed in poly⁡(n+G+G1)2O(s/log⁡s)\operatorname{poly}(n+G+G_1)2^{O(s/\sqrt{\log s})} bit operations and poly⁡(n+G+G1)\operatorname{poly}(n+G+G_1) bits of space. In particular, if k(C)=O(G)k(C)=O(G) (for example, if the number of HH and CH\mathrm{CH} gates is O(G)O(G)), the time is 2O(G/log⁡G)poly⁡(n+G1)2^{O(G/\sqrt{\log G})}\operatorname{poly}(n+G_1) and the space is polynomial.

Corollary 1.3 (Marginals and exact sampling; proved on paper only). Assume (HA)(\mathrm{H}_{\mathrm{A}}), and let C,xC,x be as in Theorem 1.1, including its condition on LL. For every S⊆[n]S \subseteq[n] and yS∈{0,1}Sy_S \in\{0,1\}^{S}, the marginal probability p(yS)=∑y[n]∖S∣⟨y∣C∣x⟩∣2p(y_S) = \sum_{y_{[n]\setminus S}} \lvert\langle y|C|x\rangle\rvert^2 can be computed exactly, as (r0+r12)/2k(C)(r_0+r_1\sqrt{2})/2^{k(C)} with r0,r1∈Zr_0,r_1 \in\mathbb{Z}, within 2c′G/log⁡Gpoly⁡(n+G1)2^{c'G/\sqrt{\log G}}\operatorname{poly}(n+G^1) operations, for an absolute constant c′c'. Consequently there is a Las Vegas algorithm that outputs y∈{0,1}ny \in\{0,1\}^{n} with probability exactly ∣⟨y∣C∣x⟩∣2\lvert\langle y|C|x\rangle\rvert^2 in expected time n⋅2c′G/log⁡Gpoly⁡(n+G1)n \cdot2^{c'G/\sqrt{\log G}}\operatorname{poly}(n+G_1). Under (HO)(\mathrm{H}_{\mathrm{O}}) the same holds in polynomial space when k(C)=O(G)k(C)=O(G).

Corollary 1.4 (Constant and slowly growing depth; proved on paper only). Assume (HA)(\mathrm{H}_{\mathrm{A}}). Amplitudes, marginals and exact samples of depth-dd circuits over GG on n≥4n \ge4 qubits can be computed within 2O(dn/log⁡n)2^{O(dn/\sqrt{\log n})} operations; the exponent is o(n)o(n) whenever d=o(log⁡n)d=o(\sqrt{\log n}).

Proposition 1.5 (The rate is that of #3-SAT; proved on paper only). There is a polynomial-time map from 3-CNFs φ\varphi with nn variables and mm clauses to circuits CφC_\varphi over GG of depth O(1)O(1) with O(n+m)O(n+m) qubits and O(n+m)O(n+m) gates, and to an integer n′n', such that ⟨0∣Cφ∣0⟩=#φ/2n′\langle0|C_\varphi|0\rangle=\#\varphi/2^{n'}. Hence an algorithm for exact amplitudes of GG-gate circuits in time T(G)T(G) counts the models of 3-CNFs with m≥nm \ge n clauses in time T(O(m))poly⁡(m)T(O(m))\operatorname{poly}(m).

So, assuming the refutations, the exact amplitude problem for circuits with O(n)O(n) gates is exactly as hard as #3-SAT with O(n)O(n) clauses, up to constants in the exponent: a bound 2o(G/log⁡G)2^{o(G/\sqrt{\log G})} in Theorem 1.1 would improve the rate of [A] for m=O(n)m=O(n). Before 2026 the same reduction showed, under the counting version #ETH [13], that 2Ω(n)2^{\Omega(n)} time was needed. Table 1 compares Theorem 1.1 with the classical methods for exact strong simulation.

Against the polynomial-space bound (2d)O(n)=2O(nlog⁡d)(2d)^{O(n)}=2^{O(n\log d)} of Aaronson and Chen [2], Theorem 1.2 gives 2O(dn/log⁡n)2^{O(dn/\sqrt{\log n})} in the worst case G=Θ(dn)G=\Theta(dn), which is asymptotically smaller only when d/log⁡d=O(log⁡n)d/\log d=O(\sqrt{\log n}).

Contributions

Our contributions are:

  • We prove that one deterministic Turing machine outputs the exact amplitudes of Clifford+TT circuits with GG two-qubit gates in time 2O(G/log⁡G)2^{O(G/\sqrt{\log G})}, assuming (HA)(\mathrm{H}_{\mathrm{A}}) (Theorem 1.1, § 3), and we machine-check this theorem in Lean on top of [A]’s development (§ 10.2).

  • We give the route through [O] in polynomial space, and show that its restriction k(C)=O(G)k(C)=O(G) comes from the precision cap of [O]’s ring (Theorem 1.2, § 4).

  • We derive exact marginals, exact sampling and the depth-dd bound, and prove a converse reduction from #3-SAT, so a better rate here would improve the rate of [A] for m=O(n)m=O(n) (Theorems 1.3 to 1.5, § 5).

  • We check the fine-grained supremacy statements we found, 24 rows in all, against their primary sources, and record which are refuted, which lose their premise and which survive; for this we prove a counting form of the width reduction of Calabro, Impagliazzo and Paturi (Theorem 6.1, § 6).

  • We explain why cubic phases resist the clause-count-free bound of [A], and prove the weaker 2-adic truncation that does transfer (Theorem 7.1, § 7).

  • We check the reductions in exact arithmetic on random small circuits (§ 10).

Mathematically, Theorem 1.2 is an immediate corollary of the appendix theorem of [O], which contracts Boolean tensor networks with one-bit wires: anyone who reads that theorem next to Markov and Shi will see it. What we add is the care needed to make it a theorem: the coefficient ring, the precision, which comes from Theorem 3.2, the resulting restriction k(C)=O(G)k(C) = O(G) in the route through [O] and its absence in the route through [A], marginals and exact sampling, the converse reduction, the check of the literature, and the analysis of cubic phases. We claim the careful statement and proof, not the idea.

Status of each claim

1.5. Status of each claim. We separate four kinds of claim. Machine-checked in Lean: Theorem 1.1 only, as EthFollowup.Quantum.amplitude_TM_rate’ (the bound under the condition on LL), EthFollowup.Quantum.amplitude_TM_rate (all circuits) and EthFollowup.Quantum.amplitude_TM (every δ>0\delta> 0), with the axioms propext, Classical.choice and Quot.sound only. The sources are in our public repository [30], https://github.com/xennygrimmato/eth-hamiltonicity, at commit ca17a37 (§ 10.2). The bound under the condition on LL and the bound for all circuits are separate Lean theorems; the Lean proof of the former reuses the machine of the latter, so one machine gives both. In Lean, (HA)(\mathrm{H}_{\mathrm{A}}) is not a hypothesis: [A]’s Lean development proves the forms of the three imported results that our Lean proof uses (§ 2), and the kernel checks them together with ours. Proved on paper only: everything else, in particular Theorem 1.2, Theorems 1.3 and 1.4, Theorem 1.5, Theorem 6.1, Theorem 6.2, Theorem 7.1 and the entries of the list in § 6. We do not claim these as machine-checked. Assumed: on paper, (HA)(\mathrm{H}_{\mathrm{A}}), (HO)(\mathrm{H}_{\mathrm{O}}) or the main theorem of [A] or [O], as each result states; in Lean, nothing beyond the three axioms. Either way the results stand or fall with the claimed refutations (§ 9). The bound after Theorem 7.1 is a sketch and is not claimed. Checked numerically: the reductions of § 3 and § 4 and the marginals of Theorem 1.3, on circuits with n≤6n \le6 and G≤20G \le20, and Theorem 7.1, on cubic phases with n≤7n \le7 and m≤14m \le14, all in exact arithmetic (§ 10). The circuit of Theorem 1.5 is not scripted; only its analogue for IQP circuits is. The running example is checked by hand.

Outline

Outline. § 2 recalls the imported results. § 3 proves Theorem 1.1, § 4 gives the route through [O], and § 5 the corollaries and the converse. § 6 checks the fine-grained supremacy literature, and § 7 treats cubic phases. § 8 to § 10 give related work, limitations, and the numerical checks and the formalisation.

The imported results

In this section we recall the definitions and results of [A] and [O] that we use, in the notation of [A] except where a letter clashes with ours, as in [29]. We state them exactly as we use them, so that (HA)(\mathrm{H}_{\mathrm{A}}) and (HO)(\mathrm{H}_{\mathrm{O}}) can be checked against the sources. Conventions of [A]: log⁡=log⁡2\log= \log_2 and [x=u]∈{0,1}[x = u] \in\{0,1\}.

Definition 2.1 ([A, Definition 1]). A function F ⁣:X×X′→ZF \colon X \times X' \to\mathbb{Z} factors through cc bits if there are f ⁣:X×{0,1}c→Zf \colon X \times\{0,1\}^{c} \to\mathbb{Z} and f′ ⁣:X′×{0,1}c→Zf' \colon X' \times\{0,1\}^{c} \to\mathbb{Z} with F(z,z′)=∑θ∈{0,1}cf(z,θ)f′(z′,θ)F(z,z') = \sum_{\theta\in\{0,1\}^{c}} f(z,\theta)f'(z',\theta). In particular, for any h ⁣:X→{0,1}ch \colon X \to\{0,1\}^{c} and h′ ⁣:X′→{0,1}ch' \colon X' \to\{0,1\}^{c} the function [h(z)=h′(z′)]=∑θ[h(z)=θ][h′(z′)=θ][h(z) = h'(z')] = \sum_{\theta}[h(z) = \theta][h'(z') = \theta] factors through cc bits.

Definition 2.2 ([A, Definition 2]). Let P,B,τ≥1P, B, \tau\ge1 be integers. A network on JJ with parameters P,B,τP,B,\tau consists of a finite set JJ of nodes; for each node aa a finite set SaS_a of settings and a function Wa ⁣:Sa→ZW_a \colon S_a \to\mathbb{Z}; and for each unordered pair {a,b}\{a,b\} of distinct nodes a function Ψa,b ⁣:Sa×Sb→Z\Psi_{a,b} \colon S_a \times S_b \to\mathbb{Z} (with Ψb,a(zb,za)=Ψa,b(za,zb)\Psi_{b,a}(z_b,z_a) = \Psi_{a,b}(z_a,z_b)), an integer ca,b=cb,a≥0c_{a,b} = c_{b,a} \ge0 and a factorisation (fa,b,fb,a)(f_{a,b},f_{b,a}) of Ψa,b\Psi_{a,b} through ca,bc_{a,b} bits. The settings of aa are encoded as bit strings of one fixed length Pa≤PP_a \le P, and (i) ∑b≠aca,b≤P\sum_{b \ne a} c_{a,b} \le P for every a∈Ja \in J; (ii) all values of all WaW_a, Ψa,b\Psi_{a,b}, fa,bf_{a,b} are below 2B2^B in absolute value; (iii) an algorithm receives PP, BB and the numbers PaP_a, ca,bc_{a,b} as input, and can decide membership in SaS_a and compute any value of WaW_a, Ψa,b\Psi_{a,b}, fa,bf_{a,b} with at most τ\tau operations.

The value of the network is Z=∑z∈∏aSa∏a∈JWa(za)∏{a,b}Ψa,b(za,zb)Z = \sum_{z \in\prod_a S_a} \prod_{a \in J} W_a(z_a) \prod_{\{a,b\}} \Psi_{a,b}(z_a,z_b), and its total is U=∑(a,b): a≠bca,bU = \sum_{(a,b):\,a \ne b} c_{a,b}.

In the (K,λ)(K,\lambda)-clique of [A, Definition 3] each of KK nodes holds a λ\lambda-bit string for each other node, cw,w′=λc_{w,w'} = \lambda, Ψw,w′=[zw(w′)=zw′(w)]\Psi_{w,w'} = [z_w(w') = z_{w'}(w)], and the tables ∣Ww∣<2B|W_w| < 2^B are queried at unit cost.

Imported Lemma 2.3 ([A, Lemma 8, case (i)]). Let K≥2K \ge2, g=3K8g=3K^{8} and ζ=1/(2K)\zeta=1/(2K). Take a network on [g][g] with parameters PP, BB, τ\tau and total UU, and put L=⌈(U+gP/K)/K2⌉L=\left\lceil(U+gP/K)/K^{2}\right\rceil and B′=⌈2gP/K⌉+g2BB^{\prime}=\left\lceil2gP/K\right\rceil+g^{2}B. If an algorithm outputs the value of every (K,L)(K,L)-clique with ∣Ww∣<2B′\lvert W_w\rvert< 2^{B^{\prime}} within TT operations, then the value of the network is computed within

(Kg+T⋅22gP/K⋅τ)⋅(g(P+B))C1(1)\left(K^{g}+T\cdot2^{2gP/K}\cdot\tau\right)\cdot\left(g(P+B)\right)^{C_{1}} \tag*{(1)}

operations, by one procedure, for an absolute constant C1C_{1}.

Imported Theorem 2.4 ([A, Theorem A]). Suppose there are pieces with parameters (ρ,c)(\rho,c) in the sense of [A, Definition 7], and let α′=log⁡(1/ρ)\alpha^{\prime}=\log(1/\rho). Fix an integer x0≥18x_{0}\ge18 with 2x0≥c+52^{x_{0}}\ge c+5 and define K0=22x0K_{0}=2^{2x_{0}}, ζℓ=2−ℓ−5\zeta_{\ell}=2^{-\ell-5}, gℓ=23ℓ+x0Kℓg_{\ell}=2^{3\ell+x_{0}}K_{\ell}, Kℓ+1=gℓKℓK_{\ell+1}=g_{\ell}^{K_{\ell}}, e0=1e_{0}=1 and eℓ+1=(ρ+3ζℓ)eℓ+3/Kℓe_{\ell+1}=(\rho+3\zeta_{\ell})e_{\ell}+3/K_{\ell}.

(a) log⁡Kℓ=(3x0+6)2ℓ−6ℓ−2x0−6\log K_{\ell}=(3x_{0}+6)2^{\ell}-6\ell-2x_{0}-6, so Kℓ<Kℓ+1≤Kℓ10K_{\ell}<K_{\ell+1}\le K_{\ell}^{10}; and 12ρℓ≤eℓ≤0.85ρℓ\frac{1}{2}\rho^{\ell}\le e_{\ell}\le0.85\rho^{\ell}, where ρℓ≤((3x0+6)/log⁡Kℓ)α′\rho^{\ell}\le((3x_{0}+6)/\log K_{\ell})^{\alpha^{\prime}}.

(b) One deterministic algorithm, given ℓ≥0\ell\ge0, λ\lambda, B≥1B\ge1 and query access to tables with ∣Ww∣<2B\lvert W_w\rvert<2^{B}, outputs the value of the (Kℓ,λ)(K_{\ell},\lambda)-clique within

2eℓKℓ2λ⋅(Kℓ(λ+B))C0Kℓ2(2)2^{e_{\ell}K_{\ell}^{2}\lambda}\cdot\left(K_{\ell}(\lambda+B)\right)^{C_{0}K_{\ell}^{2}} \tag*{(2)}

operations, with C0C_{0} depending only on cc.

Imported Theorem 2.5 ([A, Theorem B]). There are pieces with parameters (1/2,c)(1/\sqrt{2},c) for an absolute constant cc.

Imported Theorem 2.6 ([O, Appendix, Boolean tensor contraction], up to notation). A Boolean tensor network is a bipartite multigraph in which each edge carries a single bit; each vertex vv has an ordered list δ(v)\delta(v) of incident edges and a tensor TvT_v assigning a ring element to each Boolean assignment to those edges; its contraction is Z=∑x∈{0,1}E∏vTv(x∣δ(v))Z=\sum_{x\in\{0,1\}^{E}}\prod_v T_v(x|_{\delta(v)}), an isolated vertex contributing its scalar entry. Fix positive constants CC, AA, bb. Let s≥2s\ge2 be an integer, and let LL be the total length of the graph and tensor descriptions of a bipartite tensor network with at most CsCs edges and maximum degree at most ⌈s9/10⌉\lceil s^{9/10}\rceil. Suppose that the number of vertices is bounded by a fixed polynomial in ss, and that the entries belong to

BP0,ℓ=(Z/2P0Z)[V]/(Vℓ),\mathcal{B}_{P_{0},\ell}=(\mathbb{Z}/2^{P_{0}}\mathbb{Z})[V]/(V^{\ell}),

1≤P0,ℓ≤Cs1\le P_{0},\ell\le Cs.

Suppose, uniformly for every vertex and Boolean choice of its ports, that its entry can be computed from the description in time (s+L)b2As9/10(s+L)^{b}2^{As^{9/10}} and space (s+L)b(s+L)^{b}. Then the contraction can be computed deterministically in (s+L)O(1)2O(s/log⁡s)(s+L)^{O(1)}2^{O(s/\sqrt{\log s})} bit operations and (s+L)O(1)(s+L)^{O(1)} bits of space, with constants depending only on CC, AA, bb.

[A]’s Lean development proves Theorem 2.3 (in a form weaker by a polynomial factor), Theorem 2.4(b) and Theorem 2.5 for arbitrary networks in its own format, and our Lean proof of Theorem 1.1 applies them directly to the circuit network (§ 10.2). The appendix of [O] is not formalised, so Theorem 1.2 is proved on paper only.

Proof of Theorem 1.1: one bit per wire segment

In this section we write the amplitude as a sum over the bits of the wire segments (§ 3.1), map its coefficients from Z[ω]\mathbb{Z}[\omega] into the integers, read the sum as a network in the sense of [A] (§ 3.3), and bound the running time (§ 3.4).

The sum over paths on wire segments

For a qubit qq let γq,1,…,γq,rq\gamma_{q,1},\ldots,\gamma_{q,r_q} be the two-qubit gates acting on qq, in circuit order; qq is touched if rq≥1r_q\ge1. Every single-qubit gate on a touched line is multiplied into the next two-qubit gate on that line, and the single-qubit gates after γq,rq\gamma_{q,r_q} are multiplied into γq,rq\gamma_{q,r_q} from the left. This is legitimate because a single-qubit gate on line qq commutes with every gate that does not act on qq. The result is a list of GG merged gates N~i=Ni/2κi\widetilde{N}_i=N_i/\sqrt{2}^{\kappa_i} with N~i∈M4(Z[ω])\widetilde{N}_i\in M_4(\mathbb{Z}[\omega]), where κi\kappa_i is the sum of the κ\kappa of the gates merged into it. An untouched line contributes the scalar ⟨yq∣Rq∣xq⟩=uq/2κq\langle y_q|R_q|x_q\rangle=u_q/\sqrt{2}^{\kappa_q}, where RqR_q is the product of its gates. The merging costs poly⁡(n+G+G1)\operatorname{poly}(n+G+G_1) operations, since all coefficients have O(k(C)+G+G1)O(k(C)+G+G_1) bits.

A segment is a piece of line qq between γq,j\gamma_{q,j} and γq,j+1\gamma_{q,j+1}; it carries one bit σe\sigma_e. The input leg of γq,1\gamma_{q,1} is fixed to xqx_q and the output leg of γq,rq\gamma_{q,r_q} to yqy_q. With ntn_t touched lines there are $s_0=\sum_q(r_q-1)=2G-n_t\le 2G segments. For a merged gate ii write σin(i),σout(i)∈{0,1}2\sigma_{\mathrm{in}}(i),\sigma_{\mathrm{out}}(i) \in\{0,1\}^{2} for the values on its input and output legs (a segment bit or a fixed xq,yqx_q,y_q).

Lemma 3.1 (Sum over paths). Put kt:=∑iκ^ik_t := \sum_i \hat{\kappa}_i and A:=∑σ∈{0,1}s0∏i=1GN^i[σout(i),σin(i)]∈Z[ω]A := \sum_{\sigma\in\{0,1\}^{s_0}} \prod_{i=1}^{G} \hat{N}_i[\sigma_{\mathrm{out}}(i),\sigma_{\mathrm{in}}(i)] \in\mathbb{Z}[\omega]. Then

⟨y∣C∣x⟩=(∏q untoucheduq)⋅A2k(C),k(C)=kt+∑q untouchedκq,\langle y|C|x\rangle= \left(\prod_{q\ \mathrm{untouched}} u_q\right)\cdot\frac{A}{\sqrt{2^{k(C)}}}, \qquad k(C)=k_t+\sum_{q\ \mathrm{untouched}}\kappa_q,

and A=2kt⟨yt∣Ct∣xt⟩A=\sqrt{2^{k_t}}\langle y_t|C_t|x_t\rangle, where CtC_t is the circuit CC restricted to the touched lines.

Proof. Insert ∑b∈{0,1}∣b⟩⟨b∣\sum_{b\in\{0,1\}}|b\rangle\langle b| on every segment and expand the matrix product. □

Coefficients: a Galois bound and a Kronecker map

For α=∑jajωj∈Z[ω]\alpha=\sum_j a_j\omega^j\in\mathbb{Z}[\omega] write ∥α∥∞=max⁡j∣aj∣\|\alpha\|_\infty=\max_j|a_j| and ∥α∥1=∑j∣aj∣\|\alpha\|_1=\sum_j|a_j|. Multiplication by ωj\omega^j permutes the coefficients up to sign, so ∥αβ∥1≤∥α∥1∥β∥1\|\alpha\beta\|_1\le\|\alpha\|_1\|\beta\|_1.

Lemma 3.2 (Galois bound). Let MM be a product of gates from G\mathcal{G} (on any number of qubits), let κ\kappa be the sum of their κU\kappa_U, and let NN be the product of their numerators, so M=N/2κM=N/\sqrt{2^\kappa}. Then every entry α\alpha of NN satisfies ∥α∥∞≤2κ/2\|\alpha\|_\infty\le2^{\kappa/2}.

Proof. The Galois group of Q(ω)/Q\mathbb{Q}(\omega)/\mathbb{Q} consists of σt:ω↦ωt\sigma_t:\omega\mapsto\omega^t, t∈{1,3,5,7}t\in\{1,3,5,7\}; it is abelian and contains complex conjugation σ7\sigma_7. Hence, for a unitary UU with entries in Q(ω)\mathbb{Q}(\omega), σt(U)σt(U)∗=σt(UU∗)=I\sigma_t(U)\sigma_t(U)^*=\sigma_t(UU^*)=I, so σt(U)\sigma_t(U) is unitary. Applying σt\sigma_t entrywise to all factors, and using σt(2)=±2\sigma_t(\sqrt{2})=\pm\sqrt{2}, we get σt(N)=(±2)κσt(M)\sigma_t(N)=(\pm\sqrt{2})^\kappa\sigma_t(M) with σt(M)\sigma_t(M) unitary, so ∣σt(α)∣≤2κ/2|\sigma_t(\alpha)|\le2^{\kappa/2}. Finally ∑tωrt=4[r≡0]\sum_t\omega^{rt}=4[r\equiv0] for r∈{−3,…,3}r\in\{-3,\ldots,3\}, which gives aj=14∑tσt(α)ω−jta_j=\frac14\sum_t\sigma_t(\alpha)\omega^{-jt} and ∣aj∣≤2κ/2|a_j|\le2^{\kappa/2}. □

By Theorem 3.1 and Theorem 3.2 (applied to CtC_t), ∥A∥∞≤2kt/2\|A\|_\infty\le2^{k_t/2}, and every entry of every N^i\hat{N}_i has ∥⋅∥∞≤2κ^i/2\|\cdot\|_\infty\le2^{\hat{\kappa}_i/2}, hence ∥⋅∥1≤4⋅2κ^i/2\|\cdot\|_1\le4\cdot2^{\hat{\kappa}_i/2}.

Lemma 3.3 (Kronecker map). Let b≥2b\ge2 and Mb:=24b+1M_b:=2^{4b}+1. The substitution x↦2bx\mapsto2^b defines a ring homomorphism κb:Z[ω]=Z[x]/(x4+1)→Z/MbZ\kappa_b:\mathbb{Z}[\omega]=\mathbb{Z}[x]/(x^4+1)\to\mathbb{Z}/M_b\mathbb{Z}, because 24b≡−1(modMb)2^{4b}\equiv-1\pmod{M_b}. It is injective on {α:∥α∥∞≤2b−1−1}\{\alpha:\|\alpha\|_\infty\le2^{b-1}-1\}, and α\alpha is recovered from κb(α)\kappa_b(\alpha) by taking the representative in (−Mb/2,Mb/2](-M_b/2,M_b/2] and reading its balanced base-2b2^b digits.

Proof. For α\alpha in the box, ∣∑jaj2bj∣≤(2b−1−1)(1+2b+22b+23b)<24b−1<Mb/2|\sum_j a_j2^{bj}|\le(2^{b-1}-1)(1+2^b+2^{2b}+2^{3b})<2^{4b-1}<M_b/2, and balanced digits are unique. □

Put Kron⁡b(α):=∑jaj2bj∈Z\operatorname{Kron}_b(\alpha):=\sum_j a_j2^{bj}\in\mathbb{Z}, an integer representative of κb(α)\kappa_b(\alpha), and fix b:=⌈kt/2⌉+2b:=\lceil k_t/2\rceil+2, so that ∥A∥∞≤2b−1−1\|A\|_\infty\le2^{b-1}-1.

No ring embedding is available (x4=−1x^4=-1 has no solution modulo 44, nor modulo VV); a homomorphism injective on a box containing the answer suffices.

The network

Fix g≥1g\ge1 and a partition of the merged gates into groups V1,…,VgV_1,\ldots,V_g with ∣Va∣≤⌈G/g⌉|V_a|\le\lceil G/g\rceil (some groups may be empty). A segment is internal to aa if both its end gates lie in VaV_a, and crossing otherwise. Let int⁡a\operatorname{int}_a be the number of segments internal to aa and PaP_a the number of crossing segments with an end in VaV_a.

Definition 3.4. The network N(C,x,y,(Va),b)\mathcal{N}(C,x,y,(V_a),b) on [g][g] has settings Sa={0,1}Pa\mathcal{S}_a=\{0,1\}^{P_a}, the bits of the crossing segments at aa;

Wa(za):=Kron⁡b(∑internal bits of a∏i∈VaN^i[σout(i),σin(i)]),W_a(z_a):=\operatorname{Kron}_b\left(\sum_{\text{internal bits of }a}\prod_{i\in V_a}\hat{N}_i[\sigma_{\mathrm{out}}(i),\sigma_{\mathrm{in}}(i)]\right),

the sum taken in Z[ω]\mathbb{Z}[\omega] over the 2int⁡a2^{\operatorname{int}_a} values of the internal segment bits, with the crossing bits read from zaz_a; and Ψa,b(za,zb)=∏e crossing between a and b[za(e)=zb(e)]\Psi_{a,b}(z_a,z_b)=\prod_{e\text{ crossing between }a\text{ and }b}[z_a(e)=z_b(e)] over the segments ee crossing between aa and bb, which factors through ca,b:=#{e crossing between a and b}c_{a,b}:=\#\{e\text{ crossing between }a\text{ and }b\} bits by Theorem 2.1, with fa,b∈{0,1}f_{a,b}\in\{0,1\}.

Example 2 (Running example, continued). In Example 1 both lines are touched, so there are no untouched scalars. The merged gates are Γ^1=CNOT⁡(H⊗I)\hat{\Gamma}_1=\operatorname{CNOT}(H\otimes I), with κ^1=1\hat{\kappa}_1=1, and Γ^2=CZ⁡\hat{\Gamma}_2=\operatorname{CZ}, so s0=2G−nt=2s_0=2G-n_t=2 segments. Since N^1[σ,00]=1\hat{N}_1[\sigma,00]=1 for σ∈{00,11}\sigma\in\{00,11\} and 00 otherwise, and CZ⁡[11,11]=−1\operatorname{CZ}[11,11]=-1, Theorem 3.1 gives A=−1A=-1 and ⟨11∣C∣00⟩=−1/2\langle11|C|00\rangle=-1/\sqrt{2}. Take g=2g=2, V1={1}V_1=\{1\} and V2={2}V_2=\{2\}. Both segments cross, so P1=P2=c1,2=2P_1=P_2=c_{1,2}=2 and U=4≤4GU=4\le4G. Here kt=1k_t=1, so b=3b=3 and Mb=212+1M_b=2^{12}+1. The value of the network is Z=∑σW1(σ)W2(σ)=−1Z=\sum_\sigma W_1(\sigma)W_2(\sigma)=-1, and Theorem 3.3 decodes A=−1A=-1.

Proposition 3.5 (Parameters and value). Let P:=4⌈G/g⌉+1P := 4\lceil G/g\rceil+ 1, B:=3b+2+4⌈G/g⌉+⌈kt/2⌉B := 3b + 2 + 4\lceil G/g\rceil+ \lceil k_t/2\rceil and τ:=22⌈G/g⌉poly⁡(G,k(C))\tau:= 2^{2\lceil G/g\rceil}\operatorname{poly}(G,k(C)). Then N=N(C,x,y,(Va),b)N = \mathcal{N}(C,x,y,(V_a),b) is a network on [g][g] with parameters P,B,τP,B,\tau in the sense of Theorem 2.2, its total is U≤4GU \le4G, and its value ZZ satisfies Z≡κb(A)(modMb)Z \equiv\kappa_b(A) \pmod{M_b}. Hence AA is obtained from ZZ by Theorem 3.3.

Proof. Condition (i). ∑b≠aca,b=Pa\sum_{b\ne a} c_{a,b}=P_a, and every merged gate has four legs, so Pa≤4∣Va∣≤PP_a \le4|V_a| \le P. This is where bounded gate arity enters; nothing about connectivity is used. The total. U=∑aPa=2⋅#{crossing segments}≤2s0≤4GU=\sum_a P_a=2\cdot\#\{\text{crossing segments}\}\le2s_0\le4G. Condition (ii). The values of Ψ\Psi and ff are 0 or 1. For WaW_a, the inner sum αa\alpha_a has ∥αa∥∞≤∥αa∥1≤2int⁡a∏i∈Va4⋅2κ^i/2\|\alpha_a\|_\infty\le\|\alpha_a\|_1\le2^{\operatorname{int}_a}\prod_{i\in V_a}4\cdot2^{\hat{\kappa}_i/2}, and int⁡a≤2∣Va∣\operatorname{int}_a\le2|V_a|, so log⁡∥αa∥∞≤4⌈G/g⌉+kt/2\log\|\alpha_a\|_\infty\le4\lceil G/g\rceil+k_t/2; and ∣Kron⁡b(α)∣<∥α∥∞23b+1|\operatorname{Kron}_b(\alpha)|<\|\alpha\|_\infty^{2^{3b+1}}. So ∣Wa∣<2B|W_a|<2^B. Condition (iii). Membership in SaS_a is trivial, Ψ\Psi and ff compare bits, and Wa(za)W_a(z_a) is a sum of 2int⁡a≤22⌈G/g⌉2^{\operatorname{int}_a}\le2^{2\lceil G/g\rceil} products of ∣Va∣|V_a| elements of Z[ω]\mathbb{Z}[\omega] with O(k(C)+G)O(k(C)+G)-bit coefficients, read from the precomputed merged gates. Value. Kron⁡b(α)≡κb(α)(modMb)\operatorname{Kron}_b(\alpha)\equiv\kappa_b(\alpha)\pmod{M_b} and κb\kappa_b is a ring homomorphism, so Z≡κb(∑z∏aαa(za)∏{a,b}Ψa,b)Z\equiv\kappa_b\left(\sum_z\prod_a\alpha_a(z_a)\prod_{\{a,b\}}\Psi_{a,b}\right). The factors Ψ\Psi restrict the sum to settings that agree on every crossing segment, and the inner sums run over the internal segments, so the argument of κb\kappa_b is AA. □

Running time

The algorithm merges the gates, computes ktk_t, bb and the untouched scalars, chooses a level ℓ\ell (below), puts K:=KℓK:=K_\ell and g:=3K8g:=3K^8, takes any partition into groups of size at most ⌈G/g⌉\lceil G/g\rceil, computes PP, BB, PaP_a, ca,bc_{a,b}, and runs the procedure of Theorem 2.3 on N\mathcal{N} with the algorithm of Theorem 2.4(b) at level ℓ\ell as the clique algorithm. It reduces the value ZZ modulo MbM_b, decodes AA by Theorem 3.3, and multiplies by the untouched scalars.

Hypotheses. Theorem 2.3 needs a network on [g][g] with integer parameters P,B,τ≥1P,B,\tau\ge1 and g=3K8g=3K^8; Theorem 3.5 supplies it. It then needs an algorithm for every (K,L)(K,L)-clique with ∣W∣<2B′|W|<2^{B'}: Theorem 2.4(b) at level ℓ\ell, with λ=L≥1\lambda=L\ge1 and B=B′≥1B=B'\ge1, gives T=2eℓK2L(K(L+B′))C0K2T=2^{e_\ell K^2L(K(L+B'))C_0K^2}. Theorem 2.4 needs pieces, which Theorem 2.5 supplies with ρ=1/2\rho=1/\sqrt{2}, α′=1/2\alpha'=1/2, and c,x0c,x_0 absolute. The circuit reaches the clique algorithm only through the tables that Theorem 2.3 builds from the WaW_a and fa,bf_{a,b}, and Theorem 2.3 itself searches the KgK^g groupings of the gg nodes for one with small pairwise totals, so no balancing beyond the partition is needed.

The exponent. By (1) and (2), with e:=eℓ≤1e:=e_\ell\le1 and K2L≤U+gP/K+K2K^2L\le U+gP/K+K^2, the cost is 2X2^X with

X≤1+max⁡{glog⁡K, e(U+gPK+K2)+C0K2log⁡(K(L+B′))+2gPK+log⁡τ}+C1log⁡(g(P+B)).(3)X\le1+\max\left\{g\log K,\ e\left(U+\frac{gP}{K}+K^2\right)+C_0K^2\log\left(K(L+B')\right)+\frac{2gP}{K}+\log\tau\right\}+C_1\log(g(P+B)). \tag*{(3)}

For G≤G0:=K030G\le G_0:=K_0^{30}, an absolute constant, the algorithm instead sums over the ≤2G\le2G segment bits directly, in O(1)⋅poly⁡(n+G1)O(1)\cdot\operatorname{poly}(n+G_1) operations. For G≥G0G\ge G_0 let ℓ\ell be the largest index with Kℓ30≤GK_\ell^{30}\le G. Then Kℓ+1>G1/30K_{\ell+1}>G^{1/30}, and Kℓ+1≤Kℓ10K_{\ell+1}\le K_\ell^{10} gives G1/300<K≤G1/30G^{1/300}<K\le G^{1/30}, so g=3K8≤3G8/30g=3K^8\le3G^{8/30} and, by Theorem 2.4(a),

e≤0.85(3x0+6log⁡K)1/2≤0.85(300(3x0+6))1/2log⁡G=:celog⁡G.e\le0.85\left(\frac{3x_0+6}{\log K}\right)^{1/2}\le\frac{0.85\left(300(3x_0+6)\right)^{1/2}}{\sqrt{\log G}}=: \frac{c_e}{\sqrt{\log G}}.

Also gP≤g(4G/g+5)=4G+5ggP\le g(4G/g+5)=4G+5g, L≤U+gP+1=O(G)L\le U+gP+1=O(G) and B′≤2gP+1+g2B=O(G1.6(1+k(C)))B'\le2gP+1+g^2B=O(G^{1.6}(1+k(C))). Since k(C)≤G+G1k(C)\le G+G_1 and log⁡(G1+2)≤G0.9\log(G_1+2)\le G^{0.9}, we have log⁡(k(C)+2)=O(G0.9)\log(k(C)+2)=O(G^{0.9}). The terms of (3) are bounded as follows:

eU≤4ceGlog⁡G,eU\le\frac{4c_eG}{\sqrt{\log G}},
(e+2)gPK≤3(4G+5g)G1/300=O(G1−1/300),(e+2)\frac{gP}{K}\le\frac{3(4G+5g)}{G^{1/300}}=O(G^{1-1/300}),
eK2+C0K2log⁡(K(L+B′))=O(G1/15(log⁡G+log⁡(k(C)+2))=O(G0.97),eK^2+C_0K^2\log(K(L+B'))=O(G^{1/15}(\log G+\log(k(C)+2))=O(G^{0.97}),
log⁡τ≤2G/g+O(log⁡G+log⁡(k(C)+2))=O(G1−8/300+G0.9),\log\tau\le2G/g+O(\log G+\log(k(C)+2))=O(G^{1-8/300}+G^{0.9}),
glog⁡K=O(G8/30log⁡G),g\log K=O(G^{8/30}\log G),
C1log⁡(g(P+B))=O(log⁡G+log⁡(k(C)+2))=O(G0.9).C_1\log(g(P+B))=O(\log G+\log(k(C)+2))=O(G^{0.9}).

Every term except the first is o(G/log⁡G)o(G/\sqrt{\log G}), uniformly in the circuit. The value ZZ has O(U+gB)O(U+gB) bits, so reducing and decoding it, and the preparation, cost poly⁡(n+G+G1)\operatorname{poly}(n+G+G_1). Hence the total is 2cG/log⁡Gpoly⁡(n+G1)2^{cG/\sqrt{\log G}}\operatorname{poly}(n+G_1) for an absolute cc (any c>4cec>4c_e works for large GG). This proves Theorem 1.1. □

Remark 3.6 (The single-qubit gate count). G1G_1 enters polynomially (input, merging) and through k(C)≤G+G1k(C)\le G+G_1, which sets the table bit length B=O(G/g+k(C))B=O(G/g+k(C)); BB reaches the exponent only via C0K2log⁡(K(L+B′))C_0K^2\log(K(L+B')), hence the hypothesis log⁡(G1+2)≤G0.9\log(G_1 + 2) \le G^{0.9}. If it fails, the direct sum over 2G2G segment bits is quasi-polynomial in the input length. The route through [A] is exact whatever k(C)k(C) is.

The route through [O]: polynomial space when k(C)=O(G)k(C)=O(G)

In this section we prove Theorem 1.2. We check the hypotheses of Theorem 2.6 one by one, with s=max⁡{4G,k(C),4}s = \max\{4G, k(C), 4\} and C=4C = 4.

  • Bipartite, one-bit edges, at most CsCs edges. The merged gates form one side. Every segment is subdivided by a binary equality tensor, as [O] does in its own CNF encoding; the equality vertices form the other side. There are 2s0≤4G≤s2s_0 \le4G \le s edges. Fixed legs (xq,yq)(x_q, y_q) are absorbed into the gate tensor, and a gate all of whose legs are fixed is an isolated vertex.

  • Maximum degree at most ⌈s9/10⌉\lceil s^{9/10}\rceil. The degree is at most 4, and ⌈s9/10⌉≥4\lceil s^{9/10}\rceil\ge4 for s≥4s \ge4.

  • Number of vertices polynomial in ss. It is G+s0≤3GG + s_0 \le3G.

  • Entries in BP0,ℓB_{P_0,\ell} with 1≤P0,ℓ≤Cs1 \le P_0,\ell\le Cs. Map x↦Vx \mapsto V. The entries of N^i\widehat{N}_i become integer polynomials of degree at most 3, reduced modulo 2P02^{P_0}, and the equality tensors have entries 0, 1. Over Z\mathbb{Z} the contraction is Q(V):=∑i∏iN^i(V)Q(V) := \sum_i \prod_i \widehat{N}_i(V), of degree at most 3G3G, so with ℓ:=3G+1≤s\ell:= 3G + 1 \le s nothing is truncated and the contraction in BP0,ℓB_{P_0,\ell} is Q(V) mod 2P0Q(V) \bmod2^{P_0}. Reducing QQ modulo V4+1V^4 + 1 over Z/2P0\mathbb{Z}/2^{P_0} gives AA modulo 2P02^{P_0} coefficientwise, because Q mod (x4+1)=AQ \bmod(x^4 + 1) = A in Z[ω]\mathbb{Z}[\omega]. By Theorem 3.2, ∣aj∣≤2kt/2<2P0−1|a_j| \le2^{kt/2} < 2^{P_0-1} for P0:=⌈k(C)/2⌉+3P_0 := \lceil k(C)/2\rceil+ 3, so the balanced residues are the aja_j. Here P0≤s/2+4≤2sP_0 \le s/2 + 4 \le2s. This is the only place where k(C)k(C) enters, and it is why s≥k(C)s \ge k(C): [O] works modulo a power of 2 with precision at most CsCs.

  • Entries computable in time (s+L)b24s9/10(s+\mathcal{L})^{b_2} 4^{s^{9/10}} and space (s+L)b(s+\mathcal{L})^b. Each entry is a lookup in a merged gate, precomputed in polynomial time and space; L=poly⁡(G+k(C))\mathcal{L} = \operatorname{poly}(G+k(C)).

The conclusion of Theorem 2.6 gives (s+L)O(1)2O(s/log⁡s)(s+\mathcal{L})^{O(1)}2^{O(s/\sqrt{\log s})} bit operations and (s+L)O(1)(s+\mathcal{L})^{O(1)} space; merging and the untouched scalars cost polynomial time and space. □

When k(C)≫Gk(C) \gg G we see no way around the precision cap inside [O]’s ring: the numerator of a two-qubit Clifford+T unitary can need arbitrarily many bits. A different route avoids the cap at a weaker rate: Bhattacharjee’s factor-system theorem [4 Thm. 5.3] sums products of bounded-scope factors with values in any commutative Q\mathbb{Q}-algebra in poly⁡(n)2O(ρ(n))\operatorname{poly}(n)2^{O(\rho(n))} operations of the algebra and polynomial storage, where ρ(n)=n(log⁡log⁡n)2/log⁡n\rho(n)=n(\log\log n)^2/\sqrt{\log n}. Applied to our path sum with the algebra Q(ω)\mathbb{Q}(\omega) it suggests polynomial space with no restriction on k(C)k(C); we have not checked the bit complexity in Q(ω)\mathbb{Q}(\omega), which that theorem does not state, and we do not claim this variant.

Corollaries, and why the rate is that of #3-SAT

In this section we prove Theorems 1.3 and 1.4 and the converse reduction, Theorem 1.5.

Marginals: the doubled network

Proof of Theorem 1.3, marginals. Untouched lines contribute ∣uq∣2/2κq|u_q|^2/2^{\kappa_q} for q∈Sq \in S and 1 for q∉Sq \notin S. On the touched lines take two copies of the merged gates, the second with every entry conjugated (a0+a1ω+a2ω2+a3ω3=a0−a3ω−a2ω2+a1ω3)(a_0+a_1\omega+a_2\omega^2+a_3\omega^3 = a_0-a_3\omega-a_2\omega^2+a_1\omega^3). Both copies get input xx. For q∈Sq \in S the output legs of both copies are fixed to yqy_q; for q∉Sq \notin S the two output legs are joined into one new trace segment. The resulting sum over paths has 2G2G merged gates and at most 2s0+nt≤4G2s_0+n_t \le4G segments, and its value is R=∑ySˉA(y)A(y)‾=2ktpt(yS)R=\sum_{y_{\bar S}} A(y)\overline{A(y)}=2^{k_t}p_t(y_S), where ptp_t is the marginal of CtC_t. RR is real, so R=r0+r12R=r_0+r_1\sqrt{2} with r0,r1∈Zr_0,r_1\in\mathbb{Z}. For every tt, σt(R)=2kt∑ySˉ∣⟨y∣σt(Ct)∣x⟩∣2∈[0,2kt]\sigma_t(R)=2^{k_t}\sum_{y_{\bar S}}|\langle y|\sigma_t(C_t)|x\rangle|^2\in[0,2^{k_t}] (the proof of Theorem 3.2), so ∥R∥∞≤2kt\|R\|_\infty\le2^{k_t} and b:=kt+2b:=k_t+2 suffices in Theorem 3.3. § 3.3 and § 3.4 now apply verbatim with 2G2G merged gates, U≤8GU\le8G, P=⌈4[G/q]⌉+1P= \lceil4[G/q] \rceil+1 and ktk_t replaced by 2kt2k_t; the exponent changes by a constant factor. Under (HO)(H_O) take s=max⁡{8G,2k(C),4}s=\max\{8G,2k(C),4\}, ℓ=6G+1\ell=6G+1 and P0=k(C)+3P_0=k(C)+3. □

Proof of Theorem 1.3, sampling. Draw u∈[0,1]u\in[0,1] uniformly, revealing its binary digits only when needed, and output the y∈{0,1}ny\in\{0,1\}^n whose interval in the lexicographic partition of [0,1)[0,1) into intervals of lengths ∣⟨y∣C∣x⟩∣2|\langle y|C|x\rangle|^2 contains uu; then Pr⁡[y]=∣⟨y∣C∣x⟩∣2\Pr[y]=|\langle y|C|x\rangle|^2 exactly. The digits y1,y2,…y_1,y_2,\ldots are found one at a time: if the interval of the prefix y<jy_{<j} is [L,L+p(y<j0))[L,L+p(y_{<j}0)), then yj=0y_j=0 iff u<L+p(y<j0)u<L+p(y_{<j}0). This takes nn marginal computations, with endpoints in Z[2]/2k(C)\mathbb{Z}[\sqrt{2}]/2^{k(C)}, compared exactly with dyadic rationals (compare squares of integers). A comparison is undecided after tt digits with probability at most 2−t2^{-t}, so log⁡n+O(1)\log n + O(1) digits are needed in expectation. □

Proof of Theorem 1.4. G≤dn/2G \le dn/2 and G1≤dnG_1 \le dn. If G≥2G \ge2 and log⁡(G1+2)≤G0.9\log(G_1 + 2) \le G^{0.9}, use Theorems 1.1 and 1.3 and G/log⁡G≤dn/log⁡nG/\sqrt{\log G} \le dn/\sqrt{\log n} (n≥4n \ge4). Otherwise G≤(log⁡(dn+2))1.12G \le(\log(dn + 2))^{1.12} and the direct sum over at most 4G4G segment bits suffices. □

The converse reduction

Proof of Theorem 1.5. Make every variable occur at most three times, parsimoniously (Tovey’s copying [28]): replace the d≥2d \ge2 occurrences of a variable vv by fresh copies v1,…,vdv_1,\ldots,v_d and add the 2-clauses (v‾1∨v2),…,(v‾d−1∨vd),(v‾d∨v1)(\overline{v}_1 \lor v_2),\ldots,(\overline{v}_{d-1} \lor v_d),(\overline{v}_d \lor v_1), which force all copies to be equal. The result φ′\varphi' has n′=O(n+m)n'=O(n+m) variables, m′=O(n+m)m'=O(n+m) clauses of width at most 3, each variable in at most three clauses, and #φ′=#φ\#\varphi'=\#\varphi.

Use one qubit per variable, one qubit yjy_j per clause and one ancilla per clause, and put D:=∏jDjD := \prod_j D_j with Dj∣x,y⟩=(−1)yj[clause j is violated by x]∣x,y⟩D_j|x,y\rangle= (-1)^{y_j[\text{clause }j\text{ is violated by }x]}|x,y\rangle. Since [clause jj violated] is the AND of the negated literals, DjD_j is a multi-controlled ZZ on yjy_j and at most three literal qubits, conjugated by XX on negative literals. A controlled-controlled-controlled-ZZ is exact over G\mathcal{G} with one ancilla returned to ∣0⟩|0\rangle: compute the AND of two controls into the ancilla with a Toffoli gate, apply CCZ, and uncompute; Toffoli and CCZ have exact Clifford+TT circuits with seven TT gates and O(1)O(1) CNOT and HH gates [22]. Every gadget touches O(1)O(1) qubits and every variable qubit lies in at most three gadgets, so the gadgets can be coloured greedily into O(1)O(1) layers of disjoint gadgets. Let Cφ=H⊗(n′+m′)DH⊗(n′+m′)C_\varphi= H^{\otimes(n'+m')}DH^{\otimes(n'+m')}, the Hadamards acting on variable and clause qubits only. Then

⟨0∣Cφ∣0⟩=2−(n′+m′)∑x,y∏j(−1)yj[clause j violated]=2−(n′+m′)∑x∏j2[clause j satisfied]=#φ2n.\langle0|C_\varphi|0\rangle= 2^{-(n'+m')}\sum_{x,y}\prod_j(-1)^{y_j[\text{clause }j\text{ violated}]} = 2^{-(n'+m')}\sum_x\prod_j2^{[\text{clause }j\text{ satisfied}]} = \frac{\#\varphi}{2^n}.

CφC_\varphi has depth O(1)O(1) and O(n+m)O(n+m) gates. Reading off #φ\#\varphi needs only the amplitude. □

This construction is not scripted; the analogue for IQP circuits [9] ⟨0∣H⊗nDfH⊗n∣0⟩=gap⁡(f)/2n\langle0|H^{\otimes n}D_fH^{\otimes n}|0\rangle=\operatorname{gap}(f)/2^n for degree-3 ff, with CCZ compiled into Clifford+TT, is checked by qc_check.py (§ 10).

Fine-grained supremacy statements after the refutations

Throughout this section “refuted” means refuted assuming the claimed refutations of ETH. Entries 1, 2, 5, 12 and 13 of the list in § 6 need only the main theorem of either [A] or [O], which is what the labs have formalised in Lean.

Two facts about SETH-type statements

ETH fails if and only if sk=0s_k=0 for every k≥3k\ge3, where sks_k is the infimum of the exponents δ\delta with kk-SAT in time 2δn2^{\delta n} [17]; then s∞=lim⁡ksk=0s_\infty=\lim_k s_k=0, so the Strong Exponential Time Hypothesis s∞=1s_\infty=1 [16] fails too. Several statements below are phrased for CNFs of unbounded width with at most cncn clauses. For the decision version, Calabro, Impagliazzo and Paturi [11] relate them to kk-SAT through Schuler’s width reduction [26]. We need the counting version, which follows from the same reduction because its two branches are disjoint.

Lemma 6.1 (Counting width reduction). Let k≥1k\ge1 and s>0s>0, and suppose a deterministic algorithm counts the models of kk-CNF with n′n' variables and m′m' clauses in time 2sn′poly⁡(n′,m′)2^{sn'}\operatorname{poly}(n',m'). Then the models of an arbitrary CNF with nn variables and mm clauses can be counted deterministically in time 2sn′exp⁡(2−sk(m+n/k))poly⁡(n,m)2^{sn'}\exp(2^{-sk}(m+n/k))\operatorname{poly}(n,m).

Proof. Write #XF\#_X F for the number of models of FF over the variable set XX. If FF has a clause C={l1,…,lk′}C=\{l_1,\ldots,l_{k'}\} with k′>kk'>k, pick C′={l1,…,lk}C'=\{l_1,\ldots,l_k\}, let F1=(F∖{C})∪{C′}F_1=(F\setminus\{C\})\cup\{C'\}, and let F0F_0 be FF with l1,…,lkl_1,\ldots,l_k set to false (satisfied clauses deleted, false literals removed). Since C′⊆CC'\subseteq C, the models of F1F_1 are exactly the models of FF that satisfy C′C'; the models of FF that violate C′C' are the extensions of the models of F0F_0 by l1=⋯=lk=falsel_1=\cdots=l_k=\text{false}. These two sets are disjoint, so #XF=#XF1+#X∖var⁡(C′)F0\#_X F=\#_X F_1+\#_{X\setminus\operatorname{var}(C')}F_0. Recurse until no clause is longer than kk (a leaf whose formula contains the empty clause counts 00), and count each leaf with the given algorithm. Along a root-to-leaf path every left branch removes a clause longer than kk and right branches create none, so there are at most mm left branches, and a leaf reached by rr right branches has n−krn-kr variables, so r≤n/kr\le n/k. The number of leaves reached by exactly rr right branches is at most (m+rr)\binom{m+r}{r} [11 proof of Lemma 5]. Hence the time is at most poly⁡(n,m)∑r≤n/k(m+rr)2s(n−kr)≤poly⁡(n,m)2sn∑r(m+⌈n/k⌉r)2−skr≤poly⁡(n,m)2sn(1+2−sk)m+n/k\operatorname{poly}(n,m)\sum_{r\le n/k}\binom{m+r}{r}2^{s(n-kr)}\le\operatorname{poly}(n,m)2^{s^n}\sum_r\binom{m+\lceil n/k\rceil}{r}2^{-skr}\le\operatorname{poly}(n,m)2^{sn}(1+2^{-sk})^{m+n/k}. □

Corollary 6.2. Assume the main theorem of [A] or of [O]. For all c,ε>0c,\varepsilon>0 there is a deterministic algorithm that counts the models of CNFs (of any width) with nn variables and at most cncn clauses in time O(2εn)O(2^{\varepsilon n}). Hence gap⁡(f)=2n−2#f\operatorname{gap}(f)=2^n-2^{\#f} is computed in the same time.

Proof. Put s=ε/3s=\varepsilon/3 and choose kk with 2−sk(c+1)log⁡e≤ε/32^{-sk}(c+1)\log e\leq\varepsilon/3. For this fixed kk, the main theorems count the models of kk-CNFs in time 2o(n′)poly⁡(n′,m′)≤2εn′poly⁡(n′,m′)2^{o(n')} \operatorname{poly}(n',m')\leq2^{\varepsilon n'}\operatorname{poly}(n',m') for n′n' large; Theorem 6.1 gives 22εn/3poly⁡(n)2^{2\varepsilon n/3}\operatorname{poly}(n). □

Theorem 6.2 does not give CNF-SAT with m=nO(1)m=n^{O(1)} clauses of unbounded width in time 2(1−δ)npoly⁡(m)2^{(1-\delta)n}\operatorname{poly}(m): the term 2−skm2^{-skm} then forces k=Ω(log⁡n)k=\Omega(\log n), where neither main theorem is subexponential. That form of SETH is not known to follow, and some statements below rest on it.

The table

The list below gives the statements we found, with the wording checked against the primary sources (arXiv versions cited). “Witness family” means the circuits built in the proof. Refuted means the statement itself is false under the refutations; premise refuted means a theorem whose hypothesis is now false, while its conclusion is not contradicted by anything here; survives means we know no consequence of the refutations that touches it.

Fine-grained quantum-supremacy statements, assuming the refutations. Each entry gives the statement and its source, its status, and the reason.

  1. [20] Conj. 1: ETH as 2Ω(m)2^{\Omega(m)} for 3-CNF decision, mm clauses. Status: refuted. [A]: #3-SAT in 2O(m/log⁡m)2^{O(m/\sqrt{\log m})}.

  1. [20] Conj. 2 (“SETH”): deterministic CNF-SAT, ≤cn\leq cn clauses, any width, needs 2(1−a)n2^{(1-a)n}. Status: refuted. Theorem 6.2.

  1. [20] Conj. 3 (NSETH, polynomial-size circuits). Status: survives. General circuits.

  1. [20] Conj. 4–6: nondeterministic gap ≠0\ne0 for poly-size, log-depth, reversible circuits needs 2(1−a)n2^{(1-a)n}. Status: survive. poly(n)(n)-size instances; Theorem 1.1 gives nothing below 2n2^n.

  1. [20] Conj. 7: nondeterministic gap⁡(f)≠0\operatorname{gap}(f)\ne0, 3-CNF with mm clauses, needs 2Ω(m)2^{\Omega(m)}. Status: refuted. gap⁡(f)=2n−2#f\operatorname{gap}(f)=2^n-2^{\#f} in deterministic 2O(m/log⁡m)2^{O(m/\sqrt{\log m})}.

  1. [20] Thm. 6, [15] Thm. 1: ETH ⇒\Rightarrow no 2o(t)2^{o(t)} strong simulation of Clifford+T with tt T gates (additive error 2−(3t+14)/72^{-(3t+14)/7}, resp. deciding ⟨0∣C∣0⟩≠0\langle0|C|0\rangle\ne0). Status: premise refuted. Witness circuits have O(t)O(t) gates and are solved exactly in 2O(t/log⁡t)2^{O(t/\sqrt{\log t})}; with poly(n)(n) Cliffords the conclusion is open.

  1. [20] Thms. 7, 8: Conj. 7 ⇒\Rightarrow no 2o(t)2^{o(t)} multiplicative-error sampling of Clifford+T, resp. tt-qubit PBC. Status: premise refuted. Witness circuits (O(t)O(t) gates) sampled exactly in 2O(t/log⁡t)2^{O(t/\sqrt{\log t})}; conclusions open.

  1. [20] Thm. 10: Conj. 2 ⇒\Rightarrow exact probabilities with hh Hadamards not in 2(1−a)h2^{(1-a)h}. Status: premise refuted. Witness circuits have Θ(ch2)\Theta(ch^2) gates; conclusion open.

  1. [20] Thm. 12: ETH ⇒\Rightarrow efficient stabilizer rank of tt-Toffoli resources is 2Ω(t)2^{\Omega(t)}. Status: premise refuted. Conclusion open.

  1. [14] Thm. 3: CNF-SAT needs 2n−o(n)poly⁡(m)⇒2^{n-o(n)}\operatorname{poly}(m)\Rightarrow strong simulation of poly(n)(n)-size circuits to precision 2−n/22^{-n/2} needs 2n−o(n)2^{n-o(n)}. Status: survives. Premise is the CNF-SAT form with poly(m)(m) clauses, not known to follow; unconditional monotone bound unaffected.

  1. [15] Thm. 2: An O(22.2541⋅10−8N)O(2^{2.2541\cdot10^{-8}N}) simulator (NN T gates) gives 3-SAT in O(1.3n)O(1.3^n). Status: moot. The refutations give the consequent outright.

  1. [19] Conj. 1 (“SETH”, ≤cn\leq cn clauses, deterministic decision). Status: refuted. Theorem 6.2; Thm. 1 loses its premise.

  1. [19] Conj. 2: nondeterministic gap⁡(f)≠0\operatorname{gap}(f)\ne0, CNF with ≤cn\leq cn clauses, needs 2(1−a)n2^{(1-a)n} (counting form). Thm. 2 loses its premise; witness circuits have ≥nm\geq nm qubits.

  1. [19] Conj. 3–4 (OV), 5–6 (3-SUM). Status: survive. OV is implied by SETH, not equivalent; 3-SUM unrelated.

  1. [12] poly3-NSETH(a). Status: survives. Worst case has Θ(n3)\Theta(n^3) monomials (§ 7).

  1. [12] per-int-NSETH(b). Status: survives. Dense permanents.

  1. [12] poly3-ave-SBSETH(a′a'). Status: survives. Average case, dense, additive error.

  1. [12] Thm. 2: ETH ⇒\Rightarrow deterministic poly3-NONBALANCED needs 2an2^{an}. Status: premise refuted. Reduction instances have O(n)O(n) monomials; worst case open.

  1. [21] Conj. 1: random cubic, multiplicative error, NTIME[n2][n^2] oracle, 2an2^{an}. Status: survives. Dense, average case, approximate.

  1. [21] Conj. 2: random quadratic plus worst-case log-depth circuit. Status: survives. General log-depth circuits.

  1. [21] Conj. 3: random quadratic plus worst-case 3-CNF with mm clauses, 2am2^{am}. Status: survives (not shown safe). O(n2+m)O(n^2+m) gates; nothing for m=O(n)m=O(n).

  1. [2] QUATH (poly time, random circuits with >n2>n^2 gates) and HOG. Status: unaffected. Average case, approximate.

  1. [10] basic QSETH: no quantum O(2n(1−δ)/2mO(1))O(2^{n(1-\delta)/2}m^{O(1)}) CNF-SAT. Status: survives. CNF-SAT form; the kk-CNF variant of their footnote 8 is refuted.

  1. 2Ω(t)2^{\Omega(t)} TT-count hardness with unboundedly many Cliffords; χ(∣T⟩⊗t)=2Ω(t)\chi(\lvert T\rangle^{\otimes t})=2^{\Omega(t)} [8]. Status: open. ETH support gone; best known exact stabilizer-rank simulation 20.3963tpoly⁡2^{0.3963t}\operatorname{poly} [24], after [8, 6]; compressed circuits have Θ(t2)\Theta(t^2) gates.

Entry 13 was left open in our working notes: there we checked only the decision form of the width–density duality. Theorem 6.1 settles it.

Sampling versus exact amplitudes

By Theorem 1.3, any hardness claim for sampling, with multiplicative or additive error, stated as 2Ω(n)2^{\Omega(n)} for a family with G=o(nlog⁡n)G=o(n\sqrt{\log n}) gates is refuted, since exact samples are approximate ones. Nothing is said about families with G≥nlog⁡nG\ge n\sqrt{\log n}, which include every random-circuit sampling proposal at depth d≫log⁡nd\gg\sqrt{\log n}, nor about polynomial-time hardness or unconditional separations for shallow circuits [27, 7].

The estimates of [12] (208 IQP qubits and 500 gates, assuming poly3-NSETH(1/2) and a cost linear in the number of gates) extrapolate from the dense worst case to a sparse family, which is asymptotically O(n/log⁡n)O(n/\sqrt{\log n}) under the refutations; at n=208n=208 this has no force.

Cubic phases: only a 2-adic truncation transfers

Let f=∑t∈Tmtf=\sum_{t\in T}m_t over F2\mathbb{F}_2, with mtm_t monomials of degree at most 3 in nn variables and ∣T∣=m|T|=m, and gap⁡(f)=∑x(−1)f(x)\operatorname{gap}(f)=\sum_x(-1)^{f(x)}. This is the object of poly3-NSETH in [12], and ⟨0∣H⊗nDfH⊗n∣0⟩=gap⁡(f)/2n\langle0|H^{\otimes n}D_fH^{\otimes n}|0\rangle=\operatorname{gap}(f)/2^n for the diagonal circuit DfD_f built from ZZ, CZ and CCZ. Compiling CCZ into Clifford+T gives G=Θ(m)G=\Theta(m), so Theorem 1.1 computes gap⁡(f)\operatorname{gap}(f) in 2O(m/log⁡m)2^{O(m/\sqrt{\log m})}, which is 2o(n)2^{o(n)} only for m=o(nlog⁡n)m=o(n\sqrt{\log n}). For CNFs, [A] has a second bound that does not depend on mm. Does one hold for cubic phases? As far as we can see, no.

Why the clause-count-free mechanism does not transfer

As we read [A, Appendix D], the mm-free bound writes [no clause is violated at ξ\xi] by inclusion–exclusion over the set V(ξ)V(\xi) of violated clauses, truncates it at a number of marks of order n/log⁡Δn/\log\Delta, and handles the assignments with more violated clauses separately, using that they contribute zero to #SAT. For phases there is no such set: (−1)f(x)=∏t(1−2mt(x))=∑Π⊆V(x)(−2)∣Π∣=(−1)∣V(x)∣(-1)^{f(x)}=\prod_t(1-2m_t(x))=\sum_{\Pi\subseteq V(x)}(-2)^{|\Pi|}=(-1)^{|V(x)|}, where V(x)V(x) is the set of monomials equal to 1 at xx, so every assignment contributes ±1\pm1, and the assignments with many “violated” monomials carry a non-zero total in general. This is checked on random instances (check C4 in § 10). What does transfer is a weaker, 2-adic truncation.

Proposition 7.1 (2-adic truncation; proved on paper). For S⊆TS\subseteq T let V(S)V(S) be the set of variables occurring in the monomials of SS.

(a) gap⁡(f)=∑S⊆T(−2)∣S∣2n−∣V(S)∣\operatorname{gap}(f)=\sum_{S\subseteq T}(-2)^{|S|}2^{n-|V(S)|}.

(b) Every term with ∣S∣≥∣V(S)∣+2|S|\ge|V(S)|+2, in particular every term with ∣S∣≥n+2|S|\ge n+2, is divisible by 2n+22^{n+2}. Since ∣gap⁡(f)∣≤2n|\operatorname{gap}(f)|\le2^n, gap⁡(f)\operatorname{gap}(f) is the representative in (−2n+1,2n+1](-2^{n+1},2^{n+1}] of ∑∣S∣≤n+1(−2)∣S∣2n−∣V(S)∣\sum_{|S|\le n+1}(-2)^{|S|}2^{n-|V(S)|} modulo 2n+22^{n+2}.

Proof. (a) Expand ∏t(1−2mt(x))\prod_t(1-2m_t(x)) and use ∑x∏t∈Smt(x)=2n−∣V(S)∣\sum_x\prod_{t\in S}m_t(x)=2^{n-|V(S)|}. (b) The 2-adic valuation of the term of SS is ∣S∣+n−∣V(S)∣|S|+n-|V(S)|.

The truncation at n+1n+1 marks is sharp up to a constant factor: truncating at ⌊n/2⌋\lfloor n/2\rfloor marks gives wrong answers on some instances (check C3). So the budget is Θ(n)\Theta(n) marks, where [A] pays about n/log⁡Δn/\log\Delta.

Sketch (not proved here). Feeding (b) into the incidence encoding that [O] uses for sparse F2\mathbb{F}_2-polynomial systems (each monomial a test vertex of weight 1−2∏v∈tzt,v1-2\prod_{v\in t}z_{t,v}, precision 2n+22^{n+2}, ℓ=1\ell=1, incidence degree bound 3(n+1)3(n+1), and branching on the variables in more than (3m/n)log⁡n(3m/n)\sqrt{\log n} monomials) should give O(nlog⁡(em/n)+nlog⁡n)O(n\log(em/n)+\sqrt n\log n) wires, hence, under (H0_0), gap⁡(f)\operatorname{gap}(f) in time 2O(nlog⁡(em/n)/log⁡n)poly⁡(m)2^{O(n\log(em/n)/\sqrt{\log n})}\operatorname{poly}(m): modulo 2n+22^{n+2} only sets of at most n+1n+1 marked monomials matter, at about log⁡(em/n)\log(em/n) bits each. We have not checked this application line by line against the degree and entry hypotheses of [O] (the maximum degree should be n3/42O(log⁡n)≤s9/10n^{3/4}2^{O(\sqrt{\log n})}\le s^{9/10} and the entries are products of 1 and −2-2, so we expect no obstacle). If correct, it is 2o(n)2^{o(n)} if and only if m≤n2O(log⁡n)m \le n2^{O(\sqrt{\log n})}, more than Theorem 1.1 gives, and it is worse than brute force for dense ff (m=Θ(n3)m=\Theta(n^3) gives 2O(nlog⁡n)2^{O(n\sqrt{\log n})}). Making cubic phases mm-independent would decide dense poly3-NONBALANCED in 2o(n)2^{o(n)} and refute poly3-NSETH(aa) for every aa; we record it as a target, not a result.

Related work

We organise related work by idea: fine-grained supremacy, exact strong simulation, the counting engines, and permanents.

Fine-grained supremacy from ETH and SETH. Morimae and Tamaki [20, 19, 21], Huang, Newman and Szegedy [14, 15] and Dalzell, Harrow, Koh and La Placa [12] derive lower bounds on classical simulation from ETH, SETH, NSETH and their relatives. Aaronson and Chen [2] base random-circuit sampling on QUATH and HOG, and Buhrman, Patro and Speelman [10] formulate quantum analogues of SETH. In contrast, we go in the converse direction: we assume the refutations and ask which of these statements survive (§ 6).

Exact strong simulation. The state vector and the Feynman path sum are folklore; Aaronson and Chen [2 Thm. 4.1] give a recursive polynomial-space method. Markov and Shi [18] contract tensor networks in time exponential in the treewidth, and Wahl and Strelchuk [31] obtain time subexponential in the number of gates for geometrically local circuits in dimension at least 2. For circuits with tt T gates, stabilizer-rank methods run in time exponential in tt [8, 6, 24]. These methods are unconditional and, unlike ours, some are practical. In contrast, our bound is subexponential in GG in arbitrary connectivity, but it is conditional and galactic (Table 1).

The counting engines. OpenAI [23] and Anthropic [3] claim deterministic subexponential #SAT algorithms, with Lean developments, under review. Our earlier work [29] counts Hamiltonian cycles on top of [A]’s development, and we reuse its conventions. Neither manuscript mentions quantum circuits: in the source of [O] the words “quantum”, “qubit” and “amplitude” do not occur, and [A] mentions tensor networks, including Markov and Shi’s paper on simulating quantum computation by tensor contraction [18], only in its related work. Nor did we find the corollary in the press or discussion threads we searched. Closest to our work is the appendix theorem of [O] (Theorem 2.6), of which Theorem 1.2 is an immediate corollary. Bhattacharjee’s factor-system theorem over commutative Q\mathbb{Q}-algebras [4 Thm. 5.3] is also close: it does not mention quantum circuits, but it suggests a polynomial-space route without the restriction k(C)=O(G)k(C)=O(G), at the weaker rate 2O(G(log⁡log⁡G)2/log⁡G)2^{O(G(\log\log G)^2/\sqrt{\log G})} (§ 4). In contrast, we work out the coefficient ring and the precision, take the route through [A] at the rate 2O(G/log⁡G)2^{O(G/\sqrt{\log G})} without the restriction k(C)=O(G)k(C)=O(G), and machine-check that route.

Permanents and BosonSampling. BosonSampling amplitudes are permanents of n×nn \times n submatrices of Haar-random unitaries [1], and Ryser’s 2npoly⁡(n)2^n\operatorname{poly}(n) formula [25] is still essentially the best exact method (for matrices of poly⁡(n)\operatorname{poly}(n)-bit integers it has been improved only to 2n−Ω(n/log⁡n)2^{n-\Omega(\sqrt{n/\log n})} [5]). Nothing above applies to them. A qubit circuit is a sum over binary variables, one per wire segment, coupled by bounded-arity local tables, so it has O(G)O(G) shared bits. A dense permanent is a sum over permutations; written as an interferometer it has Θ(m2)\Theta(m^2) beam splitters on modes carrying up to nn photons, and any local encoding seems to need about log⁡n\log n bits per row. Informal notes in the project repository [30], not claimed here, argue that “row-local” networks for the permanent in the model of [A] (rows grouped into clusters of at most ss) have total at least nlog⁡(n/s)n\log(n/s), too much for the engine to beat 2n2^n; only sparse matrices go through. In contrast to the qubit case, the conjecture per-int-NSETH of [12] is therefore untouched (§ 6).

Limitations

The results stand or fall with the claimed refutations. We list this and the other limitations plainly.

(a) Everything algorithmic is conditional on (HA_{\mathrm{A}}), (HO_{\mathrm{O}}), or (for Theorem 6.2) the main theorem of either manuscript. Expert review of the manuscripts is pending; if the hypotheses fail, the results lose their proofs, and we offer no independent route. Theorem 1.1 is machine-checked on top of [A]’s Lean development, so it needs that development, built from source, to be sound. Our other results use the imported statements in sharp forms that are not machine-checked in the generality needed here; their proofs are on paper, written with AI assistance, and not yet checked by a human. (b) Approximate sampling; many gates. The results concern exact simulation of circuits with G=o(nlog⁡n)G=o(n\sqrt{\log n}) gates. Random-circuit sampling rests on the non-collapse of the polynomial hierarchy and on average-case conjectures about approximate amplitudes (QUATH, HOG [2]) for circuits with G≥n2G\ge n^2 gates; none of this is touched, nor are cross-entropy benchmarking and polynomial-depth families. For G≥nlog⁡nG\ge n\sqrt{\log n} the bound is no better than 2n2^n.

(c) The constants are galactic. Using [A], Theorem A(a) with the smallest admissible x0=18x_0=18, the leading term of the exponent is eℓUe\ell U with U≈4GU\approx4G when most wire segments cross the partition, and 15/log⁡G≤ℓ≤114/log⁡G15/\sqrt{\log G}\le\ell\le114/\sqrt{\log G} at the level chosen in § 3.4. For depth dd, so G≈dn/2G\approx dn/2, the leading term alone is below nn only if log⁡n\log n exceeds roughly 900d2900d^2 (with the lower estimate) to 5.2⋅104d25.2\cdot10^4d^2 (with the upper one), and lower-order terms only make this worse. [O] gives no explicit constants. Experiments with 50 to 1000 qubits are not affected.

(d) Denominators in the route through [O]. Its 2-adic precision cap forces k(C)=O(s)k(C)=O(s).

(e) Sketches. The bound after Theorem 7.1 is a sketch; the circuit of Theorem 1.5 is not scripted; the reading of [A], Appendix D is ours.

(f) Literature. Every entry of the list in § 6 was checked against the arXiv text of its source, and all citations were checked against the primary sources on 30 Oct 2026, with a further spot-check on 4 Nov 2026; theorem and conjecture numbers refer to the arXiv versions named in the bibliography. We have not searched systematically for prior statements of Theorem 1.1; given how direct it is, it may well have been noticed independently.

(g) Scale of the numerical checks. They cover n≤6n\le6 and G≤20G\le20 and test the reductions (encodings, decodings, parameters, ring emulation), not the engines, which are replaced by brute force and variable elimination.

Numerical checks and formalisation

In this section we describe the numerical checks of the reductions and the Lean development of Theorem 1.1.

Numerical checks

Numerical checks. The scripts are research/scripts/qc_*.py in the project repository [30]; they are not part of this submission. They use exact integer and Z[ω]\mathbb{Z}[\omega] arithmetic only and print PASS or FAIL; qc_lib.py implements the objects of § 3 to § 5.

qc_check.py [seed] [count] draws random circuits (n∈[2,6]n\in[2,6], G∈[1,20]G\in[1,20], G1∈[0,25]G_1\in[0,25], random gates, xx, yy, and balanced or contiguous partitions into g≤4g\le4 groups) and checks: (Q1) Theorem 3.1 against the state vector; (Q2) the integer network by brute force over the crossing bits, by the literal sum of Theorem 2.2 on small cases and by variable elimination, decoded by Theorem 3.3; (Q3) decoding the unreduced Knorrer polynomial; (Q4) condition (i), U≤4GU\le4G and ∣Wa∣≤2B|W_a|\le2^B; (Q5) the contraction in BP0,ℓ\mathcal{B}_{P_0,\ell} with P0=⌈k/2⌉+3P_0=\lceil k/2\rceil+3, ℓ=3G+1\ell=3G+1, reduced modulo V4+1V^4+1 and lifted; (Q6) Theorem 3.2, with ∣σt(A)∣2≤2k|\sigma_t(A)|^2\le2^k checked exactly in Z[2]\mathbb{Z}[\sqrt{2}]; (Q7) marginals via the doubled network, in Z[2]\mathbb{Z}[\sqrt{2}] and equal to 1 for S=∅S=\varnothing. It also checks 20 IQP instances (gap(f)/2n(f)/2^n with CCZ compiled into Clifford+T) and four ablations that must be caught: bb one bit too small, ℓ=4\ell=4, P0P_0 one bit short, and no conjugation in the doubled network. Seed 1 (150 circuits): 529 networks by brute force, 418 also by the literal sum, 65 truncated-ring brute forces, 91 marginals, at most 23 crossing segments, 48 partitions with a node of more than 14 boundary bits skipped. Seed 2: 544, 432, 62, 86, 24, 39. Both PASS, with all ablations caught; the agreement is equality of integers.

qc_cubic_gap.py checks, for random ff with n≤7n\le7, m≤14m\le14: (C1)–(C2) Theorem 7.1; (C3) the valuation bound, and failure of truncation at ⌊n/2⌋\lfloor n/2\rfloor marked as one instance; (C4) a non-zero total of the assignments with more than ⌊n/3⌋\lfloor n/3\rfloor “violated” monomials on some instance. Seeds 1 (120 instances) and 3 (400) PASS.

Formalisation

Formalisation. We machine-check Theorem 1.1 in Lean 4 on top of [A]’s development, as three theorems in the namespace EthFollowup.Quantum: amplitude_TM_rate' (the stated bound under the condition on LL), amplitude_TM_rate (all circuits) and amplitude_TM (every δ>0\delta>0). The sources are in our public repository [30], https://github.com/xennygrimmatto/eth-hamiltonicity, as of commit ca17a37, in the directories lean/EthFollowup/Quantum and lean-bridge-anthropic/bridge/EthFollowup/Quantum; they are not part of this submission. The semantics a reader must check are short: the gate matrices over C\mathbb{C}, the circuit matrix as the ordered product, and the amplitude ⟨y∣C∣x⟩\langle y|C|x\rangle (Circuit.lean, §1), and the input and output encodings (§7). The theorems depend only on the axioms propext, Classical.choice and Quot.sound, and fresh leanchecker replays of their closures (EthFollowup.Quantum.QRate, which contains QTheorem, and EthFollowup.Quantum.Paths) passed. Theorem 1.2, Theorem 1.3 and Theorem 6.2 are proved on paper only and are not claimed as machine-checked; nor are the other results listed as proved on paper in § 1.5.

Acknowledgements

This work was done with AI assistance (Claude Code, Anthropic): the observation, the proofs, the Lean development, the literature checks and the numerical verification scripts were produced in sessions with Claude Code. Theorem 1.1 is checked by the Lean kernel; everything else has been checked numerically or against the primary sources as described, but awaits human expert review. Responsibility for any errors rests with the author.

References

References

  1. [1]S. Aaronson and A. Arkhipov, The computational complexity of linear optics, Theory Comput. 9 (2013), 143–252.
  2. [2]S. Aaronson and L. Chen, Complexity-theoretic foundations of quantum supremacy experiments, Proc. 32nd Computational Complexity Conference (CCC 2017), LIPIcs, vol. 79, 2017, pp. 22:1–22:67; theorem numbers as in arXiv:1612.05903v2.arxiv.org/abs/1612.05903
  3. [3]Claude (Anthropic), A 20ⁿ time framework for 3-SAT, preprint dated 26 October 2026. Paper and Lean formalisation in the repository https://github.com/anthropics/formal-math, folder eth; we use commit 16a5ce41304fcdacae40be7c13e21c361dc26da5.
  4. [4]D. Bhattacharjee, Superexponential weighted counting and exact optimisation, Hexagon preprint hexagon:2610.00594v1, 29 October 2026, https://hexagonmath.org/abs/2610.00594.
  5. [5]A. Björklund, Below all subsets for some permutational counting problems, Proc. 15th Scandinavian Symposium and Workshop on Algorithm Theory (SWAT 2016), LIPIcs, vol. 53, 2016, pp. 17:1–17:11.arxiv.org/abs/1211.0391
  6. [6]S. Bravyi and D. Gosset, Improved classical simulation of quantum circuits dominated by Clifford gates, Phys. Rev. Lett. 116 (2016), no. 25, 250501.DOI
  7. [7]S. Bravyi, D. Gosset and R. König, Quantum advantage with shallow circuits, Science 362 (2018), no. 6412, 308–311.DOI
  8. [8]S. Bravyi, G. Smith and J. A. Smolin, Trading classical and quantum computational resources, Phys. Rev. X 6 (2016), no. 2, 021043.DOI
  9. [9]M. J. Bremner, R. Jozsa and D. J. Shepherd, Classical simulation of commuting quantum computations implies collapse of the polynomial hierarchy, Proc. R. Soc. A 467 (2011), no. 2126, 459–472.DOI
  10. [10]H. Buhrman, S. Patro and F. Speelman, The quantum strong exponential-time hypothesis, arXiv:1911.05686 (v2), 2019. Published as A framework of quantum strong exponential-time hypotheses, Proc. 38th STACS, LIPIcs, vol. 187, 2021, pp. 19:1–19:19.arxiv.org/abs/1911.05686
  11. [11]C. Calabro, R. Impagliazzo and R. Paturi, A duality between clause width and clause density for SAT, Proc. 21st Annual IEEE Conference on Computational Complexity (CCC 2006), 2006, pp. 252–260.DOI
  12. [12]A. M. Dalzell, A. W. Harrow, D. E. Koh and R. L. La Placa, How many qubits are needed for quantum computational supremacy?, Quantum 4 (2020), 264; arXiv:1805.05224.
  13. [13]H. Dell, T. Husefeldt, D. Marx, N. Taslaman and M. Wahlén, Exponential time complexity of the permanent and the Tutte polynomial, ACM Trans. Algorithms 10 (2014), no. 4, Art. 21.
  14. [14]C. Huang, M. Newman and M. Szegedy, Explicit lower bounds on strong quantum simulation, arXiv:1804.10368 (v2), 2018; published in IEEE Trans. Inform. Theory 66 (2020), no. 9, 5585–5600.
  15. [15]C. Huang, M. Newman and M. Szegedy, Explicit lower bounds on strong simulation of quantum circuits in terms of T-gate count, arXiv:1902.04764 (v2), 2019.arxiv.org/abs/1902.04764
  16. [16]R. Impagliazzo and R. Paturi, On the complexity of k-SAT, J. Comput. System Sci. 62 (2001), no. 2, 367–375.DOI
  17. [17]R. Impagliazzo, R. Paturi and F. Zane, Which problems have strongly exponential complexity?, J. Comput. System Sci. 63 (2001), no. 4, 512–530.
  18. [18]I. L. Markov and Y. Shi, Simulating quantum computation by contracting tensor networks, SIAM J. Comput. 38 (2008), no. 3, 963–981.DOI
  19. [19]T. Morimae and S. Tamaki, Depth-scaling fine-grained quantum supremacy based on SETH and qubit-scaling fine-grained quantum supremacy based on Orthogonal Vectors and 3-SUM, arXiv:1902.08382 (v1), 2019. Version 2 (November 2019), with R. Hayakawa as additional author, is retitled Fine-grained quantum supremacy based on Orthogonal Vectors, 3-SUM and All-Pairs Shortest Paths.
  20. [20]T. Morimae and S. Tamaki, Fine-grained quantum computational supremacy, arXiv:1901.01637 (v4), 2019; published in Quantum Inf. Comput. 19 (2019), no. 13&14, 1089–1115.arxiv.org/abs/1901.01637
  21. [21]T. Morimae and S. Tamaki, Additive-error fine-grained quantum supremacy, Quantum 4 (2020), 329; arXiv:1912.06336.DOI
  22. [22]M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2000.
  23. [23]OpenAI, The Exponential Time Hypothesis is false, preprint dated 17 October 2026, revised 25 October 2026. Paper and Lean formalisation in the repository https://github.com/openai/math.
  24. [24]H. Qassim, H. Pashayan and D. Gosset, Improved upper bounds on the stabilizer rank of magic states, Quantum 5 (2021), 606; arXiv:2106.07740.DOI
  25. [25]H. J. Ryser, Combinatorial Mathematics, Carus Mathematical Monographs, no. 14, Mathematical Association of America, 1963.
  26. [26]R. Schuler, An algorithm for the satisfiability problem of formulas in conjunctive normal form, J. Algorithms 54 (2005), no. 1, 40–44.DOI
  27. [27]B. M. Terhal and D. P. DiVincenzo, Adaptive quantum computation, constant depth quantum circuits and Arthur–Merlin games, Quantum Inf. Comput. 4 (2004), no. 2, 134–145.arxiv.org/abs/quant-ph/0205133
  28. [28]C. A. Tovey, A simplified NP-complete satisfiability problem, Discrete Appl. Math. 8 (1984), no. 1, 85–89.DOI
  29. [29]V. Tulsyan, Counting Hamiltonian cycles in sparse graphs in subexponential time, assuming the new #SAT algorithms, Hexagon preprint hexagon:2610.00672v1, October 2026; source and Lean development at https://github.com/xennygrimmato/eth-hamiltonicity.
  30. [30]V. Tulsyan, eth-hamiltonicity: Lean 4 formalisation and research notes, public repository https://github.com/xennygrimmato/eth-hamiltonicity, 2026; Lean sources as of commit ca17a37.
  31. [31]T. B. Wahl and S. Strelchuk, Simulating quantum circuits using efficient tensor network contraction algorithms with subexponential upper bound, Phys. Rev. Lett. 131 (2023), no. 18, 180601; arXiv:2208.01498.DOI

Paper details

vlnc:W_ah43q84e
v1 · latestHistory (1)

published · Fresh — added to the Hub within the last 30 days

added ·

Usage

views · 0

downloads · 0

A view counts once per browser per day, a download each time the PDF is requested.

CC BY 4.0

Original source

  • Physics and Astronomy

Anthropic Claude Fable, Anthropic Claude Opus · View statement →

Looking for related papers…
Contents