pith. sign in

arxiv: 2510.01923 · v3 · submitted 2025-10-02 · 🪐 quant-ph

The Constant Geometric Speed Schedule for Adiabatic State Preparation

Pith reviewed 2026-05-18 10:50 UTC · model grok-4.3

classification 🪐 quant-ph
keywords adiabatic state preparationconstant geometric speed scheduleenergy gap scalingadiabatic path lengthquadratic speedupquantum annealingmolecular simulationunstructured search
0
0 comments X

The pith

A constant geometric speed schedule reduces adiabatic evolution time scaling from O(Δ^{-2}) to O(L Δ^{-1}) when path length stays independent of the gap.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

Adiabatic state preparation changes a Hamiltonian slowly from an easy initial form to a hard target, keeping the system in the ground state if the change is slow enough relative to the smallest energy gap Δ. The usual linear schedule forces evolution time T to grow as the square of 1/Δ. The paper replaces that schedule with one that advances at constant speed measured by the geometric distance along the path in parameter space. This change improves the scaling to T proportional to L/Δ, where L is the total path length, whenever L itself does not grow as Δ shrinks. The authors also supply a segmented version that estimates segment lengths from eigenstate overlaps, so only a global lower bound on Δ is required instead of the full gap function. Tests on search problems and small molecules confirm the improved scaling.

Core claim

We introduce the constant geometric speed (CGS) schedule, which traverses the adiabatic path at a uniform rate. We show that this approach reduces the scaling of the evolution time by a factor of Δ^{-1}, provided L remains bounded independently of Δ. We propose a segmented CGS protocol where path segment lengths are computed from eigenstate overlaps on the fly, reducing the prior spectral-knowledge requirement from the full gap function Δ(s) to just a global lower bound on the energy gap.

What carries the argument

Constant geometric speed (CGS) schedule that advances the Hamiltonian at a uniform rate measured by geometric distance along the adiabatic path in parameter space.

If this is right

  • Evolution time T scales as O(L Δ^{-1}) rather than the conventional O(Δ^{-2}).
  • Quadratic speedup over the linear schedule is realized whenever L stays bounded as Δ approaches zero.
  • The segmented protocol requires only a global lower bound on the gap instead of the full function Δ(s).
  • Optimal Δ^{-1} scaling is observed in numerical tests on unstructured search, N₂, and a [2Fe-2S] cluster.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same constant-speed construction could be paired with existing adiabatic optimizations such as counter-diabatic driving to reduce coherence demands further.
  • Similar uniform-rate ideas may apply directly to diabatic or hybrid quantum-classical state-preparation routines that already track path geometry.
  • Hardware implementations on devices with fixed coherence times would see the largest practical gains precisely when the gap is smallest yet the path remains short.

Load-bearing premise

The adiabatic path length L can be bounded independently of the minimum energy gap Δ.

What would settle it

A concrete adiabatic problem in which the path length L grows at least linearly with 1/Δ, so that the CGS schedule loses its claimed scaling advantage over the linear schedule.

Figures

Figures reproduced from arXiv: 2510.01923 by Hyowon Park, Mancheon Han, Sangkook Choi.

Figure 1
Figure 1. Figure 1: Schematic of the constant speed schedule construction. [PITH_FULL_IMAGE:figures/full_fig_p005_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: Application of the constant speed schedule to the adiabatic Grover search. Panels (a–c) show results for [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Application of the constant speed schedule to the nitrogen molecule. (a) Low-lying energy spectra from Density Functional Theory [PITH_FULL_IMAGE:figures/full_fig_p008_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Application of the constant speed schedule to the [2Fe-2S] cluster. (a) Molecular structure. (b) Estimated adiabatic evolution time [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
read the original abstract

The efficiency of adiabatic quantum evolution is governed by the evolution time $T$, which typically scales as $\mathcal{O}(\Delta^{-2})$ with the minimum energy gap $\Delta$. However, the rigorous lower bound is $\mathcal{O}(L\Delta^{-1})$, where $L$ is the adiabatic path length. Although $L$ is formally upper-bounded by $\mathcal{O}(\Delta^{-1})$, such a bound is often too loose in practice, and $L$ can be bounded independently of $\Delta$. This indicates the potential for a quadratic speedup through adiabatic schedule construction. Here, we introduce the constant geometric speed (CGS) schedule, which traverses the adiabatic path at a uniform rate. We show that this approach reduces the scaling of the evolution time by a factor of $\Delta^{-1}$, provided $L$ remains bounded independently of $\Delta$. We propose a segmented CGS protocol where path segment lengths are computed from eigenstate overlaps on the fly, reducing the prior spectral-knowledge requirement from the full gap function $\Delta(s)$ to just a global lower bound on the energy gap. Numerical tests on adiabatic unstructured search, N$_2$, and a [2Fe-2S] cluster demonstrate the optimal $\Delta^{-1}$ scaling, confirming a quadratic speedup over the standard linear schedule.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

1 major / 2 minor

Summary. The manuscript introduces the constant geometric speed (CGS) schedule for adiabatic state preparation. Standard linear schedules yield evolution time T scaling as O(Δ^{-2}) with the minimum gap Δ. The paper notes the rigorous lower bound T = O(L Δ^{-1}) with adiabatic path length L, and claims that CGS achieves the optimal Δ^{-1} scaling (quadratic improvement) by traversing the path at uniform geometric speed, provided L remains bounded independently of Δ. A segmented CGS variant is proposed that determines segment lengths from on-the-fly eigenstate overlaps, requiring only a global lower bound on the gap rather than the full Δ(s) function. Numerical tests on adiabatic unstructured search, the N₂ molecule, and a [2Fe-2S] cluster are reported to confirm the optimal scaling.

Significance. If the condition that L is independent of Δ holds in the relevant regimes, the CGS schedule offers a practical route to quadratic speedups in adiabatic quantum algorithms while reducing the spectral information needed a priori. The segmented protocol is a notable practical contribution. The numerical demonstrations on both unstructured search and molecular systems provide concrete support for improved scaling in those instances.

major comments (1)
  1. [Abstract and Numerical Tests] Abstract and numerical results: The headline claim of a Δ^{-1} scaling improvement (quadratic speedup) is conditional on L remaining bounded independently of Δ. While the abstract correctly notes the formal bound L = O(Δ^{-1}) and asserts that L 'can be bounded independently' in practice, no explicit computation or plot of L versus Δ is provided for the tested systems (unstructured search, N₂, [2Fe-2S]). If L scales as Δ^{-1} in these examples, the total T reverts to O(Δ^{-2}) and the claimed improvement does not materialize. This verification is load-bearing for the central scaling result.
minor comments (2)
  1. [Methods] Clarify in the methods section how the geometric speed is normalized and how the on-the-fly overlap computation is implemented without introducing additional error sources.
  2. [Figures] Ensure all figures reporting scaling include error bars or confidence intervals from multiple runs, and label axes consistently with the symbols defined in the text.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful review of our manuscript on the constant geometric speed schedule for adiabatic state preparation. We address the major comment point by point below.

read point-by-point responses
  1. Referee: [Abstract and Numerical Tests] Abstract and numerical results: The headline claim of a Δ^{-1} scaling improvement (quadratic speedup) is conditional on L remaining bounded independently of Δ. While the abstract correctly notes the formal bound L = O(Δ^{-1}) and asserts that L 'can be bounded independently' in practice, no explicit computation or plot of L versus Δ is provided for the tested systems (unstructured search, N₂, [2Fe-2S]). If L scales as Δ^{-1} in these examples, the total T reverts to O(Δ^{-2}) and the claimed improvement does not materialize. This verification is load-bearing for the central scaling result.

    Authors: We thank the referee for this important observation. The numerical results demonstrate that the CGS schedule achieves T scaling as O(Δ^{-1}) for the tested systems, in contrast to the linear schedule. Given the relation T = O(L Δ^{-1}), the observed scaling implies L remains bounded independently of Δ; otherwise T would revert to O(Δ^{-2}). To make this explicit, we will add plots of the computed path length L versus Δ for the unstructured search, N₂, and [2Fe-2S] systems in the revised manuscript. These will confirm L is independent of Δ in the relevant regimes. revision: yes

Circularity Check

0 steps flagged

No significant circularity; central claim conditional on external bound for L

full rationale

The paper derives the Δ^{-1} scaling improvement for the CGS schedule directly from the known rigorous lower bound T = O(L Δ^{-1}) by constructing a uniform geometric speed traversal that saturates this bound when L is independent of Δ. This independence is explicitly treated as a provided condition rather than derived from the schedule or fitted to data; the abstract and claim acknowledge the formal O(Δ^{-1}) upper bound on L while noting it is often loose in practice. No self-definitional loops, fitted parameters renamed as predictions, or load-bearing self-citations appear in the provided derivation chain. The segmented protocol reduces spectral requirements to a global gap lower bound via on-the-fly overlaps, which is an independent algorithmic choice. Numerical demonstrations on search, N2, and [2Fe-2S] are presented as confirmation rather than inputs to the scaling result. The derivation remains self-contained against the external adiabatic theorem benchmarks.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The main unproven or assumed element is the independence of the path length from the gap, which enables the quadratic speedup claim. No free parameters or new entities are introduced in the abstract.

axioms (1)
  • domain assumption The adiabatic path length L can be bounded independently of the minimum energy gap Δ
    This condition is required for the claimed reduction in evolution time scaling to Δ^{-1}.

pith-pipeline@v0.9.0 · 5760 in / 1275 out tokens · 50162 ms · 2026-05-18T10:50:10.645435+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

44 extracted references · 44 canonical work pages · 2 internal anchors

  1. [1]

    (7) Substituting these estimates into Eq

    can be bounded using deriva- tives of the Hamiltonian and the inverse of the energy gap [ 3]: ∥ ˙P ∥ ≤ ∥ ˙H(τ)∥ ∆( τ) , ∥Q ¨PP ∥ ≤ ∥ ¨H(τ)∥ ∆( τ) + 4 ∥ ˙H(τ)∥2 ∆ 2(τ) . (7) Substituting these estimates into Eq. ( 6) yields the commonly cited scalings for the required evolution time: T = O ( ∥ ˙H∥ ∆ 2 ) + O ( ∥ ¨H∥ ∆ 2 ) + O ( ∥ ˙H∥2 ∆ 3 ) . (8) This expla...

  2. [2]

    In geometrical terms, the natural parametrization of a curve is by its arc length [ 21], defined as l(s) = ∫ s 0 ∥∂s′ |Φ(s′)⟩ ∥ds′

    to geometric properties of the path. In geometrical terms, the natural parametrization of a curve is by its arc length [ 21], defined as l(s) = ∫ s 0 ∥∂s′ |Φ(s′)⟩ ∥ds′. (9) Accordingly, the instantaneous speed along the path is v(τ) = dl dτ = ∥ |˙Φ(τ)⟩ ∥. (10) The quantity ∥ ˙P ∥ in Eq. (6) is exactly this speed v(τ): ∥ ˙P ∥ = ∥ |˙Φ(τ)⟩ ⟨Φ(τ)|+ |Φ(τ)⟩ ⟨˙Φ(...

  3. [3]

    2), the energy gap ∆( s) often becomes small only in a very narrow region of s, so a uniform grid may fail to capture the point where long evolution time is required

    However, as seen in the adiabatic Grover problem (Fig. 2), the energy gap ∆( s) often becomes small only in a very narrow region of s, so a uniform grid may fail to capture the point where long evolution time is required. In contrast, our root-finding proce- dure adaptively places grid points in proportion to the required time, ensuring that such regions a...

  4. [4]

    Furthermore, the root-finding step in Algorithm 1 need not be highly precise

    through a perturbative expansion |Φ(s)⟩ nears = 0 or can be adjusted iteratively until the fidelity f1 exceeds a desired threshold. Furthermore, the root-finding step in Algorithm 1 need not be highly precise. Theorem 1 does not require the segment lengths ∆ lj to be uniform; it is sufficient to identify a point sj+1 that yields an overlap close to the targe...

  5. [5]

    To ensure this projector can resolve the minimum energy gap, ∆ ≡ mins ∆( s), we set β = 2∆ − 1

    The projector is implemented as an integral of time evolutions weighted by a Gaussian with standard deviation β . To ensure this projector can resolve the minimum energy gap, ∆ ≡ mins ∆( s), we set β = 2∆ − 1. Since each step of the algorithm requires two such projection operations to compute the necessary overlap, we estimate the total overhead Tp asTp =...

  6. [6]

    As the initial Hamiltonian Hi, we use the Kohn-Sham Hamiltonian from a converged DFT calculation

    Nitrogen Molecule ( N2) The first system we consider is the nitrogen molecule (N 2), described using a STO-3G basis set [ 30] and restricted to the singlet subspace. As the initial Hamiltonian Hi, we use the Kohn-Sham Hamiltonian from a converged DFT calculation. The system then evolves towards the final Hamiltonian Hf , which is the full electronic Hamilto...

  7. [7]

    [2Fe-2S] cluster Finally, we demonstrate our method on a [2Fe-2S] cluster, a common bioinorganic motif that serves as a representative problem for strongly correlated systems in quantum chem- istry [ 31– 35]. Here, we show that achieving the optimal gap scaling with the constant speed schedule not only provides a speedup but also significantly enhances the...

  8. [8]

    Born and V

    M. Born and V . Fock, Beweis des adiabatensatzes, Zeitschrift f¨ur Physik 51, 165 (1928)

  9. [9]

    Kato, On the adiabatic theorem of quantum mechanics, Jour- nal of the Physical Society of Japan 5, 435 (1950)

    T. Kato, On the adiabatic theorem of quantum mechanics, Jour- nal of the Physical Society of Japan 5, 435 (1950)

  10. [10]

    Jansen, M.-B

    S. Jansen, M.-B. Ruskai, and R. Seiler, Bounds for the adiabatic approximation with applications to quantum computation, Jour- nal of Mathematical Physics 48, 102111 (2007)

  11. [11]

    Quan- tum computation by adiabatic evolution,

    E. Farhi, J. Goldstone, S. Gutmann, and M. Sipser, Quan- tum computation by adiabatic evolution, arXiv preprint quant- ph/0001106 (2000)

  12. [12]

    I. M. Georgescu, S. Ashhab, and F. Nori, Quantum simulation, Reviews of Modern Physics 86, 153 (2014)

  13. [13]

    Albash and D

    T. Albash and D. A. Lidar, Adiabatic quantum computation, Rev. Mod. Phys. 90, 015002 (2018)

  14. [14]

    Aspuru-Guzik, A

    A. Aspuru-Guzik, A. D. Dutoi, P . J. Love, and M. Head- Gordon, Simulated quantum computation of molecular ener- gies, Science 309, 1704 (2005)

  15. [15]

    V eis and J

    L. V eis and J. Pittner, Adiabatic state preparation study of methylene, The Journal of Chemical Physics 140, 214111 (2014)

  16. [16]

    M. A. Nielsen and I. L. Chuang, Quantum computation and quantum information (Cambridge university press, 2010)

  17. [17]

    A. Y . Kitaev, Quantum measurements and the abelian stabilizer problem, arXiv preprint quant-ph/9511026 (1995)

  18. [18]

    D. S. Abrams and S. Lloyd, Quantum algorithm providing ex- ponential speed increase for finding eigenvalues and eigenvec- tors, Physical Review Letters 83, 5162 (1999)

  19. [19]

    Boixo and R

    S. Boixo and R. D. Somma, Necessary condition for the quan- tum adiabatic approximation, Phys. Rev. A 81, 032308 (2010)

  20. [20]

    Roland and N

    J. Roland and N. J. Cerf, Quantum search by local adiabatic evolution, Phys. Rev. A 65, 042308 (2002)

  21. [21]

    A. T. Rezakhani, W.-J. Kuo, A. Hamma, D. A. Lidar, and P . Za- nardi, Quantum adiabatic brachistochrone, Phys. Rev. Lett.103, 080502 (2009)

  22. [22]

    A. T. Rezakhani, D. F. Abasto, D. A. Lidar, and P . Zanardi, Intrinsic geometry of quantum adiabatic evolution and quantum phase transitions, Phys. Rev. A 82, 012321 (2010)

  23. [23]

    Isermann, On the optimal schedule of adiabatic quantum computing, Quantum Information Processing 20, 300 (2021)

    S. Isermann, On the optimal schedule of adiabatic quantum computing, Quantum Information Processing 20, 300 (2021)

  24. [24]

    Matsuura, S

    S. Matsuura, S. Buck, V . Senicourt, and A. Zaribafiyan, V ari- ationally scheduled quantum simulation, Phys. Rev. A 103, 052435 (2021)

  25. [25]

    Shingu and T

    Y . Shingu and T. Hatomura, Geometrical scheduling of adia- batic control without information of energy spectra, Phys. Rev. A 112, 022410 (2025)

  26. [26]

    Quantum adiabatic optimization without heuristics

    M. Jarret, B. Lackey, A. Liu, and K. Wan, Quantum adi- abatic optimization without heuristics, arXiv preprint quant- ph/1810.04686

  27. [27]

    Braida, S

    A. Braida, S. Chakraborty, A. Chaudhuri, J. Cunningham, R. Menavlikar, L. Novo, and J. Roland, Unstructured Adiabatic Quantum Optimization: Optimality with Limitations, Quantum 9, 1790 (2025)

  28. [28]

    O’Neill, Elementary Differential Geometry , revised 2nd ed

    B. O’Neill, Elementary Differential Geometry , revised 2nd ed. (Academic Press, Burlington, 2006)

  29. [29]

    P . M. Alsing and C. Cafaro, From the classical frenet-serret apparatus to the curvature and torsion of quantum-mechanical evolutions. part ii. nonstationary hamiltonians, International Journal of Geometric Methods in Modern Physics 21, 2450151 (2024)

  30. [30]

    Elgart and G

    A. Elgart and G. A. Hagedorn, A note on the switching adi- abatic theorem, Journal of Mathematical Physics 53, 102202 (2012)

  31. [31]

    M. Han, H. Park, and S. Choi, Quantum zeno monte carlo for computing observables, npj Quantum Information 11, 46 (2025)

  32. [32]

    F. N. Fritsch and J. Butland, A method for constructing local monotone piecewise cubic interpolants, SIAM Journal on Sci- entific and Statistical Computing 5, 300 (1984)

  33. [33]

    Virtanen, R

    P . Virtanen, R. Gommers, T. E. Oliphant, M. Haber- land, T. Reddy, D. Cournapeau, E. Burovski, P . Peterson, W. Weckesser, J. Bright, S. J. van der Walt, M. Brett, J. Wil- son, K. J. Millman, N. Mayorov, A. R. J. Nelson, E. Jones, R. Kern, E. Larson, C. J. Carey, ˙ I. Polat, Y . Feng, E. W. Moore, J. V anderPlas, D. Laxalde, J. Perktold, R. Cimrman, I. H...

  34. [34]

    L. K. Grover, Quantum mechanics helps in searching for a nee- dle in a haystack, Phys. Rev. Lett. 79, 325 (1997)

  35. [35]

    A. D. Becke, Density-functional thermochemistry. III. The role of exact exchange, The Journal of Chemical Physics 98, 5648 (1993)

  36. [36]

    C. Lee, W. Y ang, and R. G. Parr, Development of the Colle- Salvetti correlation-energy formula into a functional of the elec- tron density, Physical Review B 37, 785 (1988)

  37. [37]

    R. F. Stewart, Small gaussian expansions of slater-type orbitals, The Journal of Chemical Physics 52, 431 (1970)

  38. [38]

    Beinert, R

    H. Beinert, R. H. Holm, and E. M ¨unck, Iron-sulfur clusters: Nature’s versatile electron-transfer agents, Science 277, 653 (1997)

  39. [39]

    M. K. Johnson and A. D. Smith, Iron-sulfur proteins, En- cyclopedia of Inorganic Chemistry 10.1002/0470862106.ia109 (2005)

  40. [40]

    Noodleman, T

    L. Noodleman, T. Lovell, T. Liu, F. Himo, and R. A. Torres, Functional-density-theory-based calculations of the electronic structure of iron-sulfur clusters, Current Opinion in Chemical Biology 8, 159 (2004)

  41. [41]

    Reiher, N

    M. Reiher, N. Wiebe, K. M. Svore, D. Wecker, and M. Troyer, Elucidating reaction mechanisms on quantum computers, Pro- ceedings of the National Academy of Sciences 114, 7555 (2017)

  42. [42]

    Motta and G

    M. Motta and G. K.-L. Chan, Ab initio calculations of the iron- sulfur cluster Fe 2S2: A benchmark study, Physical Review B 10 97, 235145 (2018)

  43. [43]

    The Constant Speed Schedule for Adiabatic State Preparation: Towards Quadratic Speedup without Prior Spectral Knowledge

    S. Lee, J. Lee, H. Zhai, Y . Tong, A. M. Dalzell, A. Kumar, P . Helms, J. Gray, Z.-H. Cui, W. Liu, et al., Evaluating the evi- dence for exponential quantum advantage in ground-state quan- tum chemistry, Nature Communications 14, 1952 (2023). Supplementary information for “The Constant Speed Schedule for Adiabatic State Preparation: Towards Quadratic Spee...

  44. [44]

    3749, and W− 1 denotes the − 1 branch of the LambertW function

    3802,c1 = 0. 3749, and W− 1 denotes the − 1 branch of the LambertW function. Then there exists a local maximizer E∗ ofg(E) such that |E∗ − E0| ≤ ϵu, (S.7) with ϵu = 1√ 2β a1 − √ a2 1 − 2a0a2 a2 , a0 = √ 2β ∆ e− β 2∆ 2 , a1 =r − (2β 2∆ 2 − 1)e− β 2∆ 2 , a2 =c0r +c1. (S.8) Moreover , w0 ≤ g(E∗ ) ≤ w0 ∆ ∆ − ϵu e− β 2ϵ2 u. (S.9) Proof. Without loss of general...