REVIEW 2 major objections 5 minor 2 cited by
(Sub)Exponential Quantum Speedup for Optimization
T0 review · 2 major / 5 minor · reviewed 2026-08-16 · deepseek-v4-flash
Pith's one-line read The paper proves that two plain quantum optimization algorithms — adiabatic evolution and quantum Hamiltonian descent — can have (sub)exponential query advantage over any classical algorithm.
desk verdict Discrete half is a genuine reduction-based advance; continuous half is a real but conditional result resting on an unproved claim imported from an overlapping preprint—worth refereeing with that flag. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The argument is carried by a chain of perturbative reductions. Starting from a stoquastic sparse adiabatic path with inverse-polynomial spectral gap, edge-subdivision gadgets embed its interaction graph into a hypercube, and then a second-order reduction maps it to a transverse-field diagonal Hamiltonian of the form $-\sum_i a_i X_i + D(t)$; a toy TFI Hamiltonian acting as a unary clock linearizes the path, and a small tilt extends it so that it starts as a pure TFI Hamiltonian and ends diagonal. For the continuous setting, the key object is the double-well operator $\hat X(\lambda)=-\frac{d^2}{d\xi^2}+\lambda^2(\xi^2-1/4)^2$, whose first two eigenstates form symmetric/antisymmetric pairs $|b_0\rangle,|b_1\rangle$ concentrated on one side of the well; the new error analysis shows that a diagonal Hamiltonian $D$ can be represented as a continuous function $\hat D$ over the box with error $\mathrm{poly}(n)\Lambda^{-\Omega(1)}$ rather than the naive $\Theta(2^n)$ bound.
What would settle it
Compute the two lowest eigenvalues of $-\frac{d^2}{d\xi^2}+\lambda^2(\xi^2-1/4)^2$ on $[-1,1]$ with vanishing boundary conditions for large $\lambda$ and check that the gap is $2/\Lambda$ with $\Lambda=e^{\lambda/6(1\pm o(1))}$ and that $\langle b_0|b_0\rangle_{[-1,w]}\le O(\Lambda^{-2/3})$; failure of any of these at increasing $\lambda$ would invalidate the error estimate $\varepsilon\le O(\sqrt{n}\,\Lambda^{-1/6}\|D\|)$ in Lemma 10.5.
Extended reading notes
Core claim
The central claim is Theorem 4.1: for every $n$ there is an $n$-qubit transverse-field Ising Hamiltonian $H_{\mathrm{TFI}}$ and a family of diagonal Hamiltonians $D$ such that the linear path $H(t)=(1-t)H_{\mathrm{TFI}}+tD$ has exponential quantum advantage. Any classical algorithm must use $\exp(n^{\Omega(1)})$ queries to find the diagonal entry $u$ minimizing $\langle u|D|u\rangle$, while the path keeps spectral gap $n^{\Omega(1)}$ and norm $\mathrm{poly}(n)$, and simulating the Schrödinger dynamics for $T=\mathrm{poly}(n)$ yields a state $n^{-\Omega(1)}$-close to $|u\rangle$ with $\mathrm{poly}(n)$ queries and gates. The continuous analogue (Theorem 10.2) transfers the same separation to functions $f$ on a bounded box, where the QHD-type dynamics $H(t)=-\Delta+t f(x)+\nu(t)g(x)$ with $\nu(t)=e^{-\Theta(\sqrt{t})}$ finds the minimizer in polynomial time while classical query algorithms require exponentially many queries.
Load-bearing premise
The continuous part of the construction imports a spectral-gap and concentration bound for the double-well operator from prior work; if that bound has hidden parameter dependencies or errors, the continuous separation does not follow.
Editorial extensions
If this is right
- Any classical query algorithm for the constructed discrete family needs $\exp(n^{\Omega(1)})$ queries even to reach $1/\mathrm{poly}(n)$ accuracy, so no classical black-box method can match the adiabatic algorithm on these instances.
- The plain linear-path adiabatic evolution suffices for the speedup; the construction needs only the standard mixer $D$ plus a fixed TFI starting Hamiltonian, not problem-specific or nonlocal mixers.
- The continuous construction gives the first provable super-polynomial quantum-classical separation for continuous optimization, carried by quantum Hamiltonian descent dynamics.
- The reductions preserve an inverse-polynomial spectral gap and polynomial norms throughout, so the quantum evolutions are implementable by standard digital simulation in polynomial time.
- The modular chain — sparse to hypercube to TFD to linear to Schrödinger — is a reusable template for compiling other adiabatic oracle separations into optimization separations.
Reading between the lines
- Because the separation is oracle-based, converting it into an unconditional complexity-class separation would require additional structure; removing the oracle is a natural next test of the framework.
- The continuous construction suggests that QHD's advantage in this setting is driven by the adiabatic gap and the concentration of the double-well encoding, not by non-adiabatic higher-energy exploration; this could be tested by replacing $\hat X$ with other spectrally similar operators and checking whether the speedup persists.
- The linearization gadget is essentially a quantum tunneling chain through a clock register; this may be reusable to turn other walk-based quantum speedups into optimization speedups.
- The continuous result inherits the cited double-well spectral bounds; an independent re-derivation of Claim 10.6 with explicit $\lambda$-dependencies would make the error analysis self-contained and easier to verify.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper constructs oracle-based families of optimization problems for which it claims provable exponential quantum-classical query separation. For discrete optimization, Theorem 4.1 reduces the GHV stoquastic adiabatic separation to a linear Hamiltonian path H(t)=(1-t)H_TFI+tD, where H_TFI is a transverse-field Ising Hamiltonian and D is diagonal, and shows that plain adiabatic evolution solves the problem in poly(n) time while any classical algorithm needs exp(n^{Ω(1)}) queries. The proof is a sequence of perturbative reductions: sparse stoquastic paths to hypercube paths (Lemma 4.3), hypercube paths to TFD paths (Lemma 4.4), piecewise-linear interpolation (Lemma 4.5), a toy TFI ground-state shifter (Lemma 4.6), linearization via a clock register (Lemma 4.7), and a tilt-and-extension step (Lemma 4.8). For continuous optimization, Theorem 10.2 embeds the discrete construction into Schrödinger operators H(t)=-Δ+ν(t)g(x)+t f(x), with ν(t)=e^{-Θ(√t)}, yielding a Quantum Hamiltonian Descent speedup. The continuous proof depends on Lemma 10.5, whose error analysis in turn depends on Claim 10.6, a set of spectral, symmetry, and concentration properties of the double-well operator bX(λ), imported from the overlapping preprint [ZLLW24].
Significance. If the results are correct, this is a substantial contribution: it provides the first provable exponential quantum-classical separation for optimization in the oracle model using the plain adiabatic algorithm and QHD, with hardness inherited from the GHV separation. The discrete part is a significant technical achievement: the perturbative reductions in Sections 5–9 are explicit, preserve an inverse-polynomial spectral gap, preserve classical hardness, and are stated in reusable form in Propositions B.1 and B.2. The paper also ships detailed proofs rather than relying on numerical experiment. The main caveat is that the continuous result is not self-contained at a load-bearing point: Claim 10.6 is taken from an overlapping-author preprint, and the current manuscript does not prove the quantitative bounds that feed Lemma 10.5 and Theorem 10.2. With that claim supplied or independently verified, the significance is high; as it stands, the continuous theorem is conditional on an external source.
major comments (2)
- [§10.3, Claim 10.6 and Lemma 10.5] The continuous half of the paper rests on Claim 10.6, which is imported from [ZLLW24, Theorem 4.22] rather than proved in this manuscript. The three quantitative properties in the claim—the spectral gap 2/Λ with Λ satisfying (10.16), the reflection identity (10.44), and the concentration bound (10.45)—are exactly the inputs to Claims 10.7 and 10.8, and they control the error ε≤O(√nΛ^{-1/6}∥D∥) in Lemma 10.5 and hence the gap and fidelity statements in Theorem 10.2(ii). If any of these bounds carries a hidden λ- or n-dependence, the estimates δ(λ)≥1/poly(n) and the n^{-Ω(1)} closeness in Theorem 10.2(ii) may fail. Because [ZLLW24] is an overlapping-author preprint, this is not an independent verification. Please include a self-contained proof of Claim 10.6, or state it from a peer-reviewed source with explicit uniform-in-λ constants for λ≥λ0=Θ(log n). As written, Theorem 10.2 is valid only conditionally on that external claim.
- [§10.4, proof of Claim 10.10(i)] The displayed estimate for ⟨bx|bDp|bx⟩ reads "≤⟨x|D|x⟩⟨ex|ex⟩−2∥˚x∥∥D∥"; this inequality is not implied by the preceding line, since the two cross terms can have either sign and the triangle inequality gives an upper bound with a plus sign, not a minus sign. The final O(√nΛ^{-1/3}∥D∥) bound is still recoverable by applying the triangle inequality to obtain both upper and lower bounds, but the proof as printed is mathematically incorrect. Please correct the inequality and explicitly state both bounds.
minor comments (5)
- [§10.3, Eq. (10.16)] The notation d^kΛ/dλ^k = (e^{λ/6})^{1±o_k(1)} for k∈N is not a well-defined definition of Λ; please state precisely what is assumed for each fixed k, whether the o_k(1) terms are uniform in the regime λ≥λ0, and how these derivative bounds are used in the estimates for ν(t), ν̇(t), and ν̈(t).
- [Appendix B] The first sentence contains the typo "impiclitly"; it should read "implicitly".
- [References] Weyl's inequality is cited to a Wikipedia article [Wik25]; please use a standard textbook or published reference for this classical result.
- [Theorem 1.2 vs Theorem 10.2] The informal Theorem 1.2 states T,t_f=poly(n), while the formal Theorem 10.2 states t_f=poly(n,1/η); please align the two statements.
- [§2.6] The text calls the reflection symmetry in (10.44) a "newfound symmetry property," but Claim 10.6 attributes it to [ZLLW24, Theorem 4.22]; please clarify what is new in this paper versus what is imported.
Circularity Check
No circular derivation: discrete reduction anchored to an external GHV separation; continuous extension uses an overlapping-author spectral theorem as an independent premise, not as a restatement of the conclusion.
full rationale
The paper's claimed derivation chain is a sequence of reduction theorems, not a fitted prediction. Theorem 4.1 is proven by starting from the external GHV oracle separation (Theorem 4.2, [GHV21]) and applying in-paper perturbative reductions (Lemmas 4.3 through 4.8). Each reduction explicitly preserves the spectral gap and translates classical queries back to queries of the GHV path, so the classical lower bound is inherited rather than introduced by construction. The continuous result, Theorem 10.2, relies on Claim 10.6, which the paper explicitly imports: "We present a key claim about bX based on [ZLLW24, Theorem 4.22] and its proof." That cited theorem concerns spectral and concentration properties of the double-well operator bX and is stated independently of the target quantum-classical separation; it supplies the parameter Λ, reflection symmetry, and concentration bounds used later in the new error analysis of Claims 10.7-10.11 and Lemma 10.5. While [ZLLW24] has overlapping authors and is a load-bearing premise, it is a prior mathematical theorem with its own stated assumptions, not a rename or a fit of the conclusion, and no equation in the present paper is defined in terms of its own output. The remaining steps, including the construction of bD, the oracle reductions, and the time rescaling, are explicit and self-contained once Claim 10.6 is granted. Thus no circular step can be exhibited; the paper's findings rest on external benchmarks rather than on their own inputs.
Assumptions & free parameters
assumptions (4)
- domain assumption GHV21 oracle separation (Theorem 4.2): there exist stoquastic sparse adiabatic paths with spectral gap Ω(√m), norm O(m), and classical query lower bound exp(n^{1/5-o(1)}).
- domain assumption ZLLW24 double-well operator properties (Claim 10.6, i.e., [ZLLW24, Theorem 4.22]): bX(λ) has spectral gap 2/Λ with Λ=e^{λ/6(1±o(1))}, plus concentration and reflection properties (10.44) and (10.45).
- standard math Perturbative simulation lemmas (Lemmas 3.8 to 3.10) from [BH17]: first- and second-order perturbative reductions preserve low-energy spectra.
- standard math Quantum adiabatic theorem (Theorem A.3) and real-space simulation algorithm [CLLLZ22, Theorem 3] (Lemma A.2).
Cite this review
Pith. "Pith review of (Sub)Exponential Quantum Speedup for Optimization." pith.science (2026). https://pith.science/paper/ZBZOHVDB
@misc{pith2026250414841,
author = {Pith},
title = {Pith review of: (Sub)Exponential Quantum Speedup for Optimization},
year = {2026},
howpublished = {\url{https://pith.science/paper/ZBZOHVDB}},
note = {Machine review of arXiv:2504.14841}
}
read the original abstract
We demonstrate provable (sub)exponential quantum speedups in both discrete and continuous optimization, achieved through simple and natural quantum optimization algorithms, namely the quantum adiabatic algorithm for discrete optimization and quantum Hamiltonian descent for continuous optimization. Our result builds on the Gily\'en--Hastings--Vazirani (sub)exponential oracle separation for adiabatic quantum computing. With a sequence of perturbative reductions, we compile their construction into two standalone objective functions, whose oracles can be directly leveraged by the plain adiabatic evolution and Schr\"odinger operator evolution for discrete and continuous optimization, respectively.
Figures
Figures from the paper (3 more)
Forward citations
Cited by 2 Pith papers
-
Stochastic Quantum Hamiltonian Descent
SQHD is a gate-based quantum algorithm that approximates a Lindblad dynamics blending Hamiltonian descent with stochastic component noise, giving an order-2 weak approximation and an O(1/t + eta sigma*) convergence bo...
-
Quantum Algorithms for Bandits with Knapsacks with Improved Regret and Time Complexities
Quantum algorithms for bandits with knapsacks achieve improved regret and time complexity by replacing classical sampling with quantum Monte Carlo and approximate quantum LP solving.
Reference graph
Works this paper leans on
-
[66]
Time-dependent Hamiltonian simulation with L1-norm scaling
issn: 0370-1573. doi: 10.1016/j.physrep.2024.03.002. [BCSWW20] Dominic W. Berry, Andrew M. Childs, Yuan Su, Xin Wang, and Nathan Wiebe. “Time-dependent Hamiltonian simulation with L1-norm scaling”. In: Quantum 4 (Apr. 2020), p. 254. issn: 2521-327X. doi: 10.22331/q-2020-04-20-254 . [BDL11] Sergey Bravyi, David P. DiVincenzo, and Daniel Loss. “Schrieffer–W...
arXiv 2011
-
[123]
Quantum speed-ups for solving semidefinite programs
isbn: 978-3-939897-89-7. doi: 10.4230/LIPIcs.APPROX-RANDOM.2015.110. [BS17] Fernando G.S.L. Brandao and Krysta M. Svore. “Quantum speed-ups for solving semidefinite programs”. In: 2017 IEEE 58th Annual Symposium on Foundations of Computer Science (FOCS). 2017, pp. 415–426. doi: 10.1109/FOCS.2017.45. 60 [CCD+03] Andrew M. Childs, Richard Cleve, Enrico Deot...
arXiv 2018
-
[2014]
On complexity of the quantum Ising model
arXiv: 1410.0703 [quant-ph]. [BH17] Sergey Bravyi and Matthew Hastings. “On complexity of the quantum Ising model”. In: Communications in Mathematical Physics 349.1 (Jan. 2017), pp. 1–45. issn: 1432-0916. doi: 10.1007/s00220-016-2787-4. [BKL+19] Fernando G. S. L. Brand˜ ao, Amir Kalev, Tongyang Li, et al. “Quantum SDP solvers: large speed-ups, optimality,...
arXiv 2017
-
[2019]
Convex optimization using quantum oracles
arXiv: 1904.03180 [quant-ph]. [AG24] Simon Apers and Sander Gribling. Quantum speedups for linear programming via interior point methods. 2024. arXiv: 2311.03215 [quant-ph]. [AGGW20] Joran van Apeldoorn, Andr´ as Gily´ en, Sander Gribling, and Ronald de Wolf. “Convex optimization using quantum oracles”. In: Quantum 4 (Jan. 2020), p. 220. issn: 2521- 327X....
arXiv 2024
-
[2021]
Quantum algorithm for linear systems of equations
Virtual, Italy: Association for Computing Machinery, 2021, pp. 1357–1369. isbn: 9781450380539. doi: 10.1145/3406325.3451060. [HHL09] Aram W. Harrow, Avinatan Hassidim, and Seth Lloyd. “Quantum algorithm for linear systems of equations”. In: Phys. Rev. Lett. 103 (15 Oct. 2009), p. 150502. doi: 10.1103/PhysRevLett.103.150502. [HSN+21] Matthew P. Harrigan, K...
arXiv 2019
Reviewed August 16, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.