Alphabeta Math
LemmaStatement: AI-adaptedProof: AI-adaptedPipeline-generated
How statement and proof provenance work

The first chip identifies the source of the statement or construction; the second identifies the source of its local proof or verification.

  • Literature-sourced: the exact statement appears in a cited source; only wording and notation differ.
  • AI-adapted: a semantically identical restatement of literature-sourced material, modulo indexing, notation, and boundary cases adopted by the library.
  • AI-generated: a genuinely novel statement formulated by AI, with no source for the claim itself.

These labels describe origin, not correctness: citations and verification chips remain separate evidence.

CKS stopping-time and success-mass instantiation

Statement

For the CKS dyadic eigenvalue-resolution procedure with m=⌈log⁡2κ⌉+1 stages, a supported eigenvalue λ is routed to one of two adjacent stopping intervals by the source's gapped phase-estimation test; the transition-band and estimator-failure mass must be included in the error budget. If ck are the input spectral amplitudes, the source analysis gives stage costs tj=O~(d2j), maximum time Tmax⁡=O~(dκ), and, up to the allocated approximation error, a success amplitude proportional to 1/(αmax⁡λ) on eigenvalue λ. Consequently Tav2=O~ ⁣(d2∑k∣ck∣2λk2),psucc=Θ~ ⁣(1αmax⁡∑k∣ck∣2λk2). In the explicitly declared efficient-qsvt mode of the local constructor, the proof below establishes these bounds with its actual charged interval budget. The source coefficient estimate αmax⁡=O~(κ) then gives Tav/psucc=O~(dκ). These are the source-efficient stopping quantities; the proof below establishes the corresponding finite estimated-wrapper floor, normalization, error and stopping bounds in this mode, while the retained Fourier-table mode keeps its distinct actual costs.

Facts & Assumptions

Given: The efficient-qsvt mode with m stages, n=m+1, ϕj=2−j, α=2/ϕm, its declared value/block/simulation error budgets at strengthened internal accuracy ϵ∗=min⁡(ϵ,ρ) and actual interval-cost bound H. The state-preparation workspace is clean as specified by the algorithm. In the error estimates below write ε=ϵ∗, so ε≤ϵ and ε≤ρ. Put X=A+, x=∥Xb∥, and let p denote its final good probability.

[F1]

The constructed stage circuit and resource/error budgets are those of Finite CKS sparse QLSP with variable-time inverse estimation.

[F2]

The proved OP1A polynomial and real QSVT extraction have the stated uniform scalar and input-block/phase error bounds (Gapped pseudoinverse polynomial degree and error, Real bounded polynomial extraction in QSVT, QSVT robustness, approximation error, and resources).

[F3]

The variable-time model preserves all stopped branches and defines RMS time from their orthogonal probabilities (Coherent variable-time algorithm model).

[F4]

The generic probability estimator has its promised relative-error, finite cap and confidence bounds (Relative probability estimation with a charged finite cap).

[F5]

The low-target lemma proves surviving probability at least 1/(16n) and multiplicative loss at most 2/(5n) per amplified stage (Low-target amplification has a controlled multiplicative loss).

[F6]

The clean-subspace unitary hybrid controls reference registers, coherent controls, actual inverses and repeated products along the ideal clean trajectory (Unitary hybrids on a preserved clean-work subspace). Actual work is retained coherently rather than reinitialized.

[F7]

The finite-precision value model supplies a Hermitian A~ with ∥A−A~∥op≤ϵval, and the simulator-call-site bounds sum to ηS over one A list and one A′ list (Finite CKS sparse QLSP with variable-time inverse estimation).

Proof

1.1givenF1F2F7algebra

In the finite Hermitian eigenbasis write b=∑kckvk, Avk=λkvk, where 1/κ≤∣λk∣≤1. Thus 1≤x≤κ. The supplied-permutation sparse block circuit uses the certified entry bound h=1 for the decoded Hermitian matrix A~. Hermiticity lets the row enumeration on the swapped system register serve as the column enumeration, and the supplied reverse-location map makes each matching slot unique. Since ∣A~ij∣≤1, its value flag matrix is unitary and has flag-zero amplitude A~ij on the zero input. The row and column neighbor states have coefficient 1/d on each listed slot, so their unique matching pair contributes A~ij/d to the all-zero flag compression; the support-coverage promise ensures every nonzero decoded entry is listed, any listed slot whose decoded value is zero contributes zero, and padded sentinels have no overlap. Hence the actual PUE block is exactly A~/d, within ϵval/d of the ideal A/d. On the ideal block A/d, eigenvalue λ gives normalized eigenvalue λ/d. The odd stage polynomial with gap δj=ϕj/d approximates the target multiplier sgn⁡(λ)ϕj/(2∣λ∣)=ϕj/(2λ) within its declared scalar error for ∣λ∣≥ϕj. Diluting by ϕm/ϕj gives the common ideal multiplier 1/(αλ), where α=2/ϕm. By [F2], the actual QSVT block on A~/d differs from the ideal transform of A/d by at most 4DjϵBE,j plus the declared polynomial/phase errors; these are included in the per-stage compression budget ηI. The real-extraction success projector includes its control and all zero workspace. The finite SVD therefore gives the ideal supported target Xb/α, with actual stage compression within ηI. Also 4κ≤α<8κ (the endpoint inequalities may be replaced by non-strict upper bounds).

1.2F1algebra

In the ideal comparison circuit, QPE at dyadic label size M≥64π/ϕ estimates the eigenphase modulo 2π. The exact Fourier kernel obeys Pr⁡[dcirc(y/M,λ/(2π))≥a]≤1/[2(Ma−1)]: indeed each label probability is at most 1/[4M2dcirc2], since ∣sin⁡(πu)∣≥2∣u∣ for ∣u∣≤1/2, and summing each of the two reciprocal-square tails uses ∑s≥Rs−2≤1/(R−1). Taking a=ϕ/(4π) gives failure at most 1/30. Circular phase error below ϕ/2 places the unwrapped phase magnitude on the correct side of 3ϕ/2 for each of the two promises. Since ∣λ∣≤1<π, this unwrapping is unambiguous on the good event. For r0 independent coherent QPE copies on this eigenvector, the squared amplitudes of label strings factor. If B counts incorrect comparisons, E[4B]≤(11/10)r0, hence Pr⁡[B≥r0/2]≤(11/20)r0. The declared odd r0≥8ln⁡(1/ζ2)+8 makes this at most ζ2. A reversible majority stores this decision without measurement, so the wrong-branch amplitude is at most ζ.

1.3F1F6F7algebra

For a simulator call at signed time t, differentiating e−iA(t−s)e−iA~s and integrating from 0 to t gives ∥e±iAt−e±iA~t∥op≤∣t∣∥A−A~∥op≤∣t∣ϵval. Thus call site q differs from the ideal e±iAtq by at most ξq+∣tq∣ϵval. Taking controlled direct sums or adjoints preserves that norm bound. By [F6] and the budget in [F7], each complete actual base list is within ηS of its ideal-A comparison on clean inputs, including arbitrary references; actual simulator work remains coherent.

2.1step 1.2algebra

For fixed supported λ, let j∗ be the first stage with 2ϕj∗≤∣λ∣. It exists by the chosen m, and 2j∗≤4/∣λ∣. Before stage j∗−1, the low-eigenphase promise holds, so stopping there has probability at most ζ2 per stage. At j∗ the high-eigenphase promise holds, so failing to stop has probability at most ζ2. The intermediate stage is unconstrained and either branch is allowed. Apart from a total probability at most (m+1)ζ2≤2mζ2, all histories stop at j∗−1 or j∗, and satisfy ∣λ∣≥ϕj at their inverse stage. This is the full adjacent-stage ambiguity analysis; it makes no unique-estimate assumption.

3.1F1F3F6step 1.1step 2.1algebra

Correct histories stop by time 4H/∣λ∣. Orthogonality of eigenvectors and histories gives ideal RMS square at most 16H2∑k∣ck∣2/λk2, with wrong-history mass contributing at most 16H2κ2(2mζ2). By step 1.3, the actual controlled simulations have total clean-input isometry difference at most ηS on the base list, with a reference register; the diagonal stopping-time-squared observable on the orthogonal clock histories has norm at most 16H2κ2. For unit vectors at distance at most ηS, the difference of its expectations is at most 2ηS times this norm, by expanding the two quadratic forms and applying the triangle inequality. Consequently Tav2≤16H2x2+16H2κ2(2mζ2+2ηS)≤17H2x2, using [F1]'s budgets and x≥1. Thus Tav≤5Hx, while Tmax⁡=H2m≤4Hκ≤2Hα. The inverse stage occurs only after stopping, so its unsuccessful workspace does not affect any later stopping probability.

3.2F1F2F6step 1.1step 2.1algebra

In the ideal-A comparison circuit let the joint eigenvalue/history projection select first stops with ∣λ∣≥ϕj. Its complement, including a final no-stop history, has mass at most 2mζ2 by step 2.1. On valid histories, compare the actual successful inverse compression with the common scalar X/α; summing the stage errors conservatively gives mηI. On the complement, the actual projected output has norm at most 2mζ by contractivity, while the comparison vector produced by A′Xb/α has norm at most (κ/α)2mζ≤2mζ/2. By step 1.3, simulation hybrids contribute at most ηS for A and at most (κ/α)ηS for the comparison circuit A′. Therefore, with all auxiliary inputs clean and internal F=0, ∥ΠFAB∣0⟩−A′(∣0⟩F,workXb/α)∥≤mηI+3mζ+2ηS≤ε128α. This compares projected output vectors and uses no full-space simulator bound.

4.1step 3.1step 3.2algebra

Put D=ε/(128α). Step 3.2 gives ∣p−x/α∣≤D, hence p≥x2/(4α2) and Tav/p≤10Hα. Applying the actual unitary (A′)† preserves the vector error. For nonzero vectors v,w, the triangle inequality gives ∥v/∥v∥−w/∥w∥∥≤2∥v−w∥/∥w∥; apply it with ∥w∥=x/α≥1/α to obtain conditional output error at most ε/64. The external herald is copied or measured before this uncompute; the comparison internal flag is clean and any actual residual is included in the vector error. The preparation-work register WB was clean initially and is untouched by A′ and (A′)†. Finally project all work, including WB, onto zero as specified in [F1]. Before that test, the pure joint conditional vector is within ε/64 of target times clean work, so the test fails with probability at most ε2/4096≤ρ/4096. Its successful normalized system vector is within ε/32≤ϵ of the desired vector by the same normalization inequality.

5.1F3F4F5step 4.1algebra

Let si be the original good-or-unfinished mass at interval i, defined by the algorithm's projector for 1≤i≤m: the orthogonal union of branches with F=1, T=0, and first-stop clock word C1=⋯=Cj−1=0,Cj=1 for some j≤i (later clock bits remain zero), and branches with F=0, T=0, and C1=⋯=Ci=0. Previously stopped bad branches and the terminal bad branch are excluded; since T is set after stage m, the unfinished term is empty at i=m. It is at least p because all future good mass lies in the unfinished branch; the clean WB factor does not change this mass. If the prior amplified selected mass is ri−1, then the next sampler has qi=ri−1si/si−1. At the pre-stage index set r0=s0=1: preparation has run, every clock bit and F,T is zero, and the entire state is unfinished. In particular q1=s1≥p and r0≥1/(16n), establishing the induction base. Suppose inductively ri−1≥1/(16n); then qi≥p/(16n)≥1/(64nα2)>pmin⁡. On the estimator-good event, q^i>1/(9n) gives the unamplified ri=qi>(100/101)/(9n)>1/(16n), with loss one. If 0<q^i≤1/(9n), [F5] gives the same floor and the claimed multiplicative loss. A zero estimate or exhausted bounded crossing does not occur on this event. The induction proves every promised estimator floor and the final amplifier's lower bound 1/(18n). Fresh production after every classical construction transcript and ordinary amplification preserve the exact conditional good state of the actual sampler.

6.1F1step 3.1step 3.2step 4.1step 5.1algebra∎

Step 3.1 and step 4.1 give precisely the source-scale RMS and success quantities: Tav2=O~(d2∑k∣ck∣2/λk2), p=Θ~(α−1(∑k∣ck∣2/λk2)1/2) and Tav/p=O~(dκ) when the actual quantum interval budget is H=d poly⁡(log⁡(dκ/ϵ),log⁡(1/ρ)) (with the declared sparse/control/implementation overhead). Here the concrete chosen common normalization is α=O(κ); all supplied nonstandard costs remain in H. The original Fourier/table alternate mode retains its distinct charged construction and SELECT costs; its cost is not substituted for this efficient-mode estimate.

Depends on

Used by

Dependency tree · two levels

88 results within two dependency steps of this one, each drawn at its shortest distance from it. An arrow runs from a result to what uses it, so the chart reads left to right and ends at this result, which carries a heavier outline. Every node is a link to that result. Click elsewhere on the chart to enlarge it.

Sources