Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

Quantum Information

126 missions · 62 completed

The study of information encoded in the states of quantum systems, where the qubit, superposition, and entanglement replace the classical bit and measurement is inherently probabilistic. Governed by the constraints of quantum mechanics, it underpins quantum computing and cryptography and the theory of optimal quantum measurements.

Missions

Open64Completed62All126
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 5: Ω(√N) deterministic classical queries for the HSP with N trivially intersecting subgroupsTextbook

Motivation

This mission develops Chapter 5 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

In the hidden subgroup problem (HSP) on a finite group GGG, an algorithm has black-box access to a function f:G→Sf : G \to Sf:G→S, where SSS is a finite set, promised to hide some unknown subgroup H≤GH \le GH≤G:

f(x)=f(y)  ⟺  x−1y∈H(x,y∈G).f(x) = f(y) \iff x^{-1}y \in H \qquad (x, y \in G).f(x)=f(y)⟺x−1y∈H(x,y∈G).

The goal is to determine HHH. A deterministic classical algorithm is a decision tree: it queries group elements g1,g2,…g_1, g_2, \dotsg1​,g2​,… adaptively, each query returning f(gℓ)∈Sf(g_\ell) \in Sf(gℓ​)∈S, and finally outputs a subgroup of GGG. It solves the HSP if, for every subgroup H≤GH \le GH≤G and every fff hiding HHH, its output is HHH. Its query count on input fff is the number of queries it makes on fff.

Theorem. There are an absolute constant c>0c > 0c>0 and a threshold N0N_0N0​ such that the following holds for every N≥N0N \ge N_0N≥N0​. Let GGG be a finite group and SSS a finite set with ∣S∣≥∣G∣|S| \ge |G|∣S∣≥∣G∣, and suppose GGG has a family H\mathcal HH of NNN distinct subgroups any two of which have only the identity in common, i.e. H∩K={1}H \cap K = \{1\}H∩K={1} for all distinct H,K∈HH, K \in \mathcal HH,K∈H. Then every deterministic classical algorithm that solves the HSP on GGG (with values in SSS) makes, on some input fff hiding some subgroup of GGG,

#queries(f)≥cN.\#\text{queries}(f) \ge c\sqrt{N}.#queries(f)≥cN​.

The theorem shows that the HSP is classically hard whenever the group has many trivially intersecting subgroups (for example Simon's problem), in contrast with efficient quantum algorithms for abelian groups.

Formalization Note The algorithm outputs the subgroup itself rather than a generating set, which changes nothing for a query lower bound. The condition ∣S∣≥∣G∣|S| \ge |G|∣S∣≥∣G∣ makes every subgroup of GGG hideable; the proof's adversary labels each new query with a fresh value, and with a very small SSS (say ∣S∣=1|S| = 1∣S∣=1) only H=GH = GH=G can be hidden, so the bound would fail. The constants ccc and N0N_0N0​ are chosen before the group, the set SSS, the family and the algorithm. The threshold N0N_0N0​ is needed because for the trivial group N=1N = 1N=1 and no queries are needed. The phrase "whose only common element is the identity" is read pairwise, as in the section's "trivially intersecting subgroups" and as the proof needs when it counts at most t(t−1)t(t-1)t(t−1) non-identity values gj−1gkg_j^{-1}g_kgj−1​gk​ against NNN subgroups. Read as ⋂H={1}\bigcap \mathcal H = \{1\}⋂H={1}, the theorem is false: in the cyclic group of squarefree order p1⋯pkp_1\cdots p_kp1​⋯pk​, all 2k2^k2k subgroups intersect trivially, yet k+1k+1k+1 queries determine HHH.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Theorem 5.1: Ω(√N) deterministic classical queries for the HSP with N trivially intersecting subgroups.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
1 thm2 active usersReviewed
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 21: Symmetrization: the Hamming-weight average of a multilinear polynomial is a polynomial of no larger degreeTextbook

Motivation

This mission develops Chapter 21 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

This is the symmetrization lemma of Minsky and Papert.

Let p∈R[x1,…,xn]p \in \mathbb R[x_1, \dots, x_n]p∈R[x1​,…,xn​] be multilinear, meaning that every variable has degree at most 111 in ppp. For k∈{0,1,…,n}k \in \{0, 1, \dots, n\}k∈{0,1,…,n}, let

P(k):=E∣x∣=k [p(x)]=(nk)−1∑x∈{0,1}n, ∣x∣=kp(x),P(k) := \mathbb E_{|x| = k}\,[p(x)] = \binom{n}{k}^{-1} \sum_{x \in \{0,1\}^n,\ |x| = k} p(x),P(k):=E∣x∣=k​[p(x)]=(kn​)−1x∈{0,1}n, ∣x∣=k∑​p(x),

where ∣x∣|x|∣x∣ is the Hamming weight of xxx (its number of 1s). Then there is a univariate polynomial Q∈R[k]Q \in \mathbb R[k]Q∈R[k] with deg⁡Q≤deg⁡p\deg Q \le \deg pdegQ≤degp (total degree) such that Q(k)=P(k)Q(k) = P(k)Q(k)=P(k) for every k∈{0,1,…,n}k \in \{0, 1, \dots, n\}k∈{0,1,…,n}.

Symmetrization reduces questions about multivariate polynomials that represent symmetric Boolean functions, such as parity or OR, to questions about univariate polynomials. Those questions can then be settled by counting sign changes or by Markov's inequality.

Formalization Note The statement that "PPP is a polynomial" is formalized as agreement with a univariate polynomial on the points {0,…,n}\{0, \dots, n\}{0,…,n}, the only points where PPP is defined. For these kkk the average is taken over a nonempty set.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Lemma 21.1: Polynomial method: the acceptance probability of a t-query algorithm is a polynomial of degree at most 2t.
  • Lemma 21.2: Symmetrization: the Hamming-weight average of a multilinear polynomial is a polynomial of no larger degree.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
2 thms2 active usersReviewed
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 27: Local d²-edge-colouring of a bipartite graph of maximum degree dTextbook

Motivation

This mission develops Chapter 27 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

Fix integers N,d≥0N,d\ge 0N,d≥0. There is a single colouring rule ccc, chosen before the graph, that assigns to an ordered pair of vertices α,β∈{1,…,N}\alpha,\beta\in\{1,\dots,N\}α,β∈{1,…,N} a colour c((α,s(α),Lα),(β,s(β),Lβ))∈N×Nc\big((\alpha,s(\alpha),L_\alpha),(\beta,s(\beta),L_\beta)\big)\in\mathbb N\times\mathbb Nc((α,s(α),Lα​),(β,s(β),Lβ​))∈N×N that depends only on the two vertices, their sides s(α),s(β)s(\alpha),s(\beta)s(α),s(β) of the bipartition, and their neighbour lists Lα,LβL_\alpha,L_\betaLα​,Lβ​, with the following property.

Let GGG be an undirected simple graph on the vertex set {1,…,N}\{1,\dots,N\}{1,…,N} with a bipartition s:V→{left,right}s:V\to\{\mathrm{left},\mathrm{right}\}s:V→{left,right} (every edge joins the two sides) and maximum degree at most ddd, and let LvL_vLv​ be any duplicate-free list of the neighbours of each vertex vvv. Then the map {α,β}↦c((α,s(α),Lα),(β,s(β),Lβ))\{\alpha,\beta\}\mapsto c\big((\alpha,s(\alpha),L_\alpha),(\beta,s(\beta),L_\beta)\big){α,β}↦c((α,s(α),Lα​),(β,s(β),Lβ​)) is a proper edge colouring of GGG:

  1. it is symmetric on edges, so each edge receives one colour;
  2. two distinct edges sharing a vertex receive distinct colours;
  3. the number of distinct colours used on the edges of GGG satisfies
#{ c(α,β):αβ∈E(G) }≤d2.\#\{\,c(\alpha,\beta) : \alpha\beta\in E(G)\,\}\le d^2 .#{c(α,β):αβ∈E(G)}≤d2.

In the source the colour of an edge αβ\alpha\betaαβ with α\alphaα on the left is (idx(α,β),idx(β,α))(\mathrm{idx}(\alpha,\beta),\mathrm{idx}(\beta,\alpha))(idx(α,β),idx(β,α)), where idx(α,β)\mathrm{idx}(\alpha,\beta)idx(α,β) is the position of β\betaβ in the neighbour list of α\alphaα. The lemma is the colouring step of the simulation of sparse Hamiltonians: the edges of each colour class form a matching, whose Hamiltonian is simulated separately before the pieces are combined by a product formula.

Formalization Note The source's "efficiently computable" has no machine model, so the running time of the rule is not formalized; what is kept is that one rule serves every graph and reads only the two endpoints' local data (two neighbour-list queries per edge). Colours are pairs of natural numbers rather than elements of [d]×[d][d]\times[d][d]×[d] so that the rule is total at d=0d=0d=0; the bound d2d^2d2 is on the colours actually used on edges.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Lemma 27.1: Local d²-edge-colouring of a bipartite graph of maximum degree d.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
1 thm2 active usersReviewed
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 9: Period finding over the reals: an irrational period from two Fourier samplesTextbook

Motivation

This mission develops Chapter 9 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

This is the classical post-processing step of period finding over the reals (Childs, Lemma 9.1).

Write ⌊x⌉\lfloor x \rceil⌊x⌉ for the integer nearest to the real number xxx. Let r>0r > 0r>0 be a real number (the period, possibly irrational), let NNN be a positive integer (the size of the cyclic group ZN\mathbb Z_NZN​ used for Fourier sampling), and let j,j′j, j'j,j′ be relatively prime integers with 1≤j<r1 \le j < r1≤j<r and 1≤j′<r1 \le j' < r1≤j′<r; Fourier sampling returns the integers ⌊jN/r⌉\lfloor jN/r \rceil⌊jN/r⌉ and ⌊j′N/r⌉\lfloor j'N/r \rceil⌊j′N/r⌉. If N≥3r2N \ge 3r^2N≥3r2, then:

  1. the fraction j/j′j/j'j/j′ is one of the convergents of the continued fraction expansion of the rational number
⌊jN/r⌉⌊j′N/r⌉;\frac{\lfloor jN/r \rceil}{\lfloor j'N/r \rceil};⌊j′N/r⌉⌊jN/r⌉​;
  1. the period is recovered to within one unit:
∣r−⌊jN⌊jN/r⌉⌉∣≤1.\left| r - \left\lfloor \frac{jN}{\lfloor jN/r \rceil} \right\rceil \right| \le 1.​r−⌊⌊jN/r⌉jN​⌉​≤1.

Thus, from two samples with coprime multipliers, one learns jjj from the continued fraction expansion and then an approximation of rrr to within 111.

Formalization Note The nearest-integer function is Mathlib's round (ties rounded up; the source leaves ties unspecified and its proof uses only that the rounding error is at most 1/21/21/2). Convergents are Mathlib's Real.convergent, which coincide with the convergents of the regular continued fraction. The ranges 1≤j,j′<r1 \le j, j' < r1≤j,j′<r are those the source's proof uses (they come from ⌊jN/r⌉∈ZN\lfloor jN/r \rceil \in \mathbb Z_N⌊jN/r⌉∈ZN​ and the proof's appeals to j≥1j \ge 1j≥1 and j′<rj' < rj′<r); the proof's "j<j′j < j'j<j′ without loss of generality" is not imposed. Under these hypotheses the denominator ⌊j′N/r⌉\lfloor j'N/r \rceil⌊j′N/r⌉ is at least 333, so no division by zero occurs.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Lemma 9.1: Period finding over the reals: an irrational period from two Fourier samples.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
1 thm2 active usersReviewed
🏆Completed
Captain: Goku

Quantum Computing: Lecture Notes 6: Abelian HSP: the standard algorithm samples uniformly from H^⊥Textbook

Motivation

This mission develops Chapter 6 of Quantum Computing: Lecture Notes by Ronald de Wolf. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

This is the correctness statement of the standard quantum algorithm for the Abelian hidden subgroup problem.

Let G=ZN1×⋯×ZNℓG = \mathbb Z_{N_1}\times\cdots\times\mathbb Z_{N_\ell}G=ZN1​​×⋯×ZNℓ​​ be a finite Abelian group written additively (Ni≥1N_i \ge 1Ni​≥1), with characters χg(h)=∏i=1ℓωNigihi\chi_g(h) = \prod_{i=1}^{\ell} \omega_{N_i}^{g_i h_i}χg​(h)=∏i=1ℓ​ωNi​gi​hi​​, ωN=e2πi/N\omega_N = e^{2\pi i/N}ωN​=e2πi/N, and quantum Fourier transform FG=FN1⊗⋯⊗FNℓF_G = F_{N_1}\otimes\cdots\otimes F_{N_\ell}FG​=FN1​​⊗⋯⊗FNℓ​​, ∣k⟩↦∣χk⟩=∣G∣−1/2∑g∈Gχk(g)∣g⟩|k\rangle \mapsto |\chi_k\rangle = |G|^{-1/2}\sum_{g\in G}\chi_k(g)|g\rangle∣k⟩↦∣χk​⟩=∣G∣−1/2∑g∈G​χk​(g)∣g⟩. Let SSS be a finite set, H≤GH \le GH≤G a subgroup, and f:G→Sf : G \to Sf:G→S a function hiding HHH:

f(g)=f(g′)  ⟺  g+H=g′+H.f(g) = f(g') \iff g + H = g' + H.f(g)=f(g′)⟺g+H=g′+H.

The standard algorithm starts in ∣0⟩∣0⟩|0\rangle|0\rangle∣0⟩∣0⟩ (registers of dimensions ∣G∣|G|∣G∣ and ∣S∣|S|∣S∣), prepares ∣G∣−1/2∑g∣g⟩∣0⟩|G|^{-1/2}\sum_{g}|g\rangle|0\rangle∣G∣−1/2∑g​∣g⟩∣0⟩, computes fff in superposition with one query, ∣g⟩∣0⟩↦∣g⟩∣f(g)⟩|g\rangle|0\rangle\mapsto|g\rangle|f(g)\rangle∣g⟩∣0⟩↦∣g⟩∣f(g)⟩, (measures the second register,) applies FGF_GFG​ to the first register and measures it. For every g∈Gg \in Gg∈G the outcome ggg occurs with probability

Pr⁡[output g]={∣H∣/∣G∣if χg(h)=1 for all h∈H (i.e. χg∈H⊥),0otherwise.\Pr[\text{output } g] = \begin{cases} |H|/|G| & \text{if } \chi_g(h) = 1 \text{ for all } h\in H \ (\text{i.e. } \chi_g \in H^\perp),\\ 0 & \text{otherwise.}\end{cases}Pr[output g]={∣H∣/∣G∣0​if χg​(h)=1 for all h∈H (i.e. χg​∈H⊥),otherwise.​

Since ∣H⊥∣=∣G∣/∣H∣|H^\perp| = |G|/|H|∣H⊥∣=∣G∣/∣H∣, the algorithm samples uniformly from the labels of the elements of H⊥H^\perpH⊥; each sample gives a linear constraint on HHH, which is how the algorithm generalizes Simon's algorithm, period finding and discrete logarithms.

Formalization Note The finite set SSS is given an arbitrary additive group structure so that computing fff is the unitary ∣g,b⟩↦∣g,b+f(g)⟩|g,b\rangle\mapsto|g,b+f(g)\rangle∣g,b⟩↦∣g,b+f(g)⟩, which maps ∣g,0⟩|g,0\rangle∣g,0⟩ to ∣g,f(g)⟩|g,f(g)\rangle∣g,f(g)⟩ as in the source. The intermediate measurement of the second register (step 4) is omitted: it does not change the distribution of the first register, so the output law is the first-register marginal of the state before the final measurement. The source's claim that HHH can then be found efficiently from few samples is asserted without proof and is not part of this statement.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • §6.2.2: Abelian HSP: the standard algorithm samples uniformly from H^⊥.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Ronald de Wolf, Quantum Computing: Lecture Notes, 27 August 2026. https://homepages.cwi.nl/~rdewolf/qcnotes.pdf.
1 thm2 active usersReviewed
🏆Completed
Captain: Goku

Quantum Computing: Lecture Notes 4: Phase estimation returns an exactly n-bit phase with probability 1Textbook

Motivation

This mission develops Chapter 4 of Quantum Computing: Lecture Notes by Ronald de Wolf. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

Let UUU be a unitary matrix acting on a finite-dimensional space, and let ∣ψ⟩|\psi\rangle∣ψ⟩ be a unit eigenvector of UUU with eigenvalue e2πiφe^{2\pi i\varphi}e2πiφ, where φ∈[0,1)\varphi\in[0,1)φ∈[0,1):

U∣ψ⟩=e2πiφ∣ψ⟩.U|\psi\rangle=e^{2\pi i\varphi}|\psi\rangle .U∣ψ⟩=e2πiφ∣ψ⟩.

Fix n≥0n\ge0n≥0, put N=2nN=2^nN=2n, and suppose φ\varphiφ can be written exactly with nnn bits, i.e. 2nφ=k2^n\varphi=k2nφ=k for a natural number kkk (necessarily k<Nk<Nk<N). Run phase estimation: start in ∣0n⟩∣ψ⟩|0^n\rangle|\psi\rangle∣0n⟩∣ψ⟩, apply the quantum Fourier transform FNF_NFN​ to the first nnn qubits, apply the controlled power ∣j⟩∣ψ⟩↦∣j⟩Uj∣ψ⟩|j\rangle|\psi\rangle\mapsto|j\rangle U^j|\psi\rangle∣j⟩∣ψ⟩↦∣j⟩Uj∣ψ⟩, then apply FN−1=FN∗F_N^{-1}=F_N^{*}FN−1​=FN∗​ to the first nnn qubits and measure them in the computational basis (bit strings read most significant bit first). Then

  1. after the controlled power the state is (FN∣k⟩)⊗∣ψ⟩\big(F_N|k\rangle\big)\otimes|\psi\rangle(FN​∣k⟩)⊗∣ψ⟩, i.e. 1N∑j=0N−1e2πiφj∣j⟩∣ψ⟩=FN∣2nφ⟩∣ψ⟩\frac{1}{\sqrt N}\sum_{j=0}^{N-1}e^{2\pi i\varphi j}|j\rangle|\psi\rangle=F_N|2^n\varphi\rangle|\psi\rangleN​1​∑j=0N−1​e2πiφj∣j⟩∣ψ⟩=FN​∣2nφ⟩∣ψ⟩;
  2. the measurement of the first register yields the outcome k=2nφ=φ1…φnk=2^n\varphi=\varphi_1\dots\varphi_nk=2nφ=φ1​…φn​ with probability 111.

This is the exact case of phase estimation, the subroutine underlying Shor's algorithm, quantum counting and many Hamiltonian algorithms.

Formalization Note The final state is the shared phaseEstimationState, (FN∗⊗I) Λ(U) (FN⊗I) ∣0n⟩∣ψ⟩(F_N^{*}\otimes I)\,\Lambda(U)\,(F_N\otimes I)\,|0^n\rangle|\psi\rangle(FN∗​⊗I)Λ(U)(FN​⊗I)∣0n⟩∣ψ⟩ with Λ(U)∣j⟩∣ϕ⟩=∣j⟩Uj∣ϕ⟩\Lambda(U)|j\rangle|\phi\rangle=|j\rangle U^{j}|\phi\rangleΛ(U)∣j⟩∣ϕ⟩=∣j⟩Uj∣ϕ⟩; the probability of outcome kkk on the first register is the marginal ∑b∣⟨k,b∣⋅⟩∣2\sum_b|\langle k,b|\cdot\rangle|^2∑b​∣⟨k,b∣⋅⟩∣2. The source's general-case claim (an nnn-bit approximation with additive error O(2−n)O(2^{-n})O(2−n) 'with high probability') supplies no probability constant and is not stated.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • §4.5: The n-qubit QFT has an exact circuit of at most n^2 Hadamard, controlled-R_s and swap gates.
  • §4.6: Phase estimation returns an exactly n-bit phase with probability 1.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Ronald de Wolf, Quantum Computing: Lecture Notes, 27 August 2026. https://homepages.cwi.nl/~rdewolf/qcnotes.pdf.
2 thms2 active usersReviewed
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 8: Principal ideals of Z[√d] coincide iff the generators differ by a unitTextbook

Motivation

This mission develops Chapter 8 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

Let ddd be a squarefree integer with d≥2d \ge 2d≥2, and let Z[d]={x+yd:x,y∈Z}\mathbb Z[\sqrt d] = \{x + y\sqrt d : x, y \in \mathbb Z\}Z[d​]={x+yd​:x,y∈Z}. For ξ∈Z[d]\xi \in \mathbb Z[\sqrt d]ξ∈Z[d​] write ξZ[d]\xi\mathbb Z[\sqrt d]ξZ[d​] for the principal ideal generated by ξ\xiξ.

For all ξ,ζ∈Z[d]\xi, \zeta \in \mathbb Z[\sqrt d]ξ,ζ∈Z[d​],

ξ Z[d]=ζ Z[d]  ⟺  ξ=ζϵ for some unit ϵ of Z[d].\xi\,\mathbb Z[\sqrt d] = \zeta\,\mathbb Z[\sqrt d] \iff \xi = \zeta\epsilon \text{ for some unit } \epsilon \text{ of } \mathbb Z[\sqrt d].ξZ[d​]=ζZ[d​]⟺ξ=ζϵ for some unit ϵ of Z[d​].

Thus the map ξ↦ξZ[d]\xi \mapsto \xi\mathbb Z[\sqrt d]ξ↦ξZ[d​] is multiplicatively periodic, with the units as its periods; this periodicity, with period the fundamental unit, underlies the periodic function used in Hallgren's algorithm for Pell's equation.

Formalization Note Z[d]\mathbb Z[\sqrt d]Z[d​] is Mathlib's ℤ√d, the principal ideal ξZ[d]\xi\mathbb Z[\sqrt d]ξZ[d​] is Ideal.span {ξ}, and a unit is IsUnit. The statement holds also for ξ=0\xi = 0ξ=0 (both sides then force ζ=0\zeta = 0ζ=0), so no nonzero hypothesis is added.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Proposition 8.1: Units of Z[√d] are exactly the elements of norm ±1.
  • Proposition 8.2: Principal ideals of Z[√d] coincide iff the generators differ by a unit.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
2 thms2 active usersReviewed
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 26: One compressed phase query raises ‖P ψ‖ by at most √(2/2^m)Textbook

Motivation

This mission develops Chapter 26 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

Let n,m≥0n,m\ge0n,m≥0, Σ={0,1}m\Sigma=\{0,1\}^mΣ={0,1}m, and let WWW be a finite workspace. Consider the Hilbert space with computational basis ∣x,z,w,y⟩|x,z,w,y\rangle∣x,z,w,y⟩, where x∈[n]x\in[n]x∈[n] is the query input, z∈Σz\in\Sigmaz∈Σ the query output, w∈Ww\in Ww∈W the workspace and y∈(Σ∪{⊥})ny\in(\Sigma\cup\{\perp\})^ny∈(Σ∪{⊥})n the oracle register. The compressed phase oracle is

Φ=(I⊗H⊗m⊗I⊗C⊗n) Λ (I⊗H⊗m⊗I⊗C⊗n).\Phi=(I\otimes H^{\otimes m}\otimes I\otimes C^{\otimes n})\,\Lambda\,(I\otimes H^{\otimes m}\otimes I\otimes C^{\otimes n}).Φ=(I⊗H⊗m⊗I⊗C⊗n)Λ(I⊗H⊗m⊗I⊗C⊗n).

Here HHH is the Hadamard gate, CCC is the compression unitary of §26.1 (it swaps ∣⊥⟩|\perp\rangle∣⊥⟩ and ∣0^⟩|\hat 0\rangle∣0^⟩ and fixes ∣y^⟩|\hat y\rangle∣y^​⟩ for y≠0y\ne0y=0), and Λ∣x,z,w,y⟩=∣x,z⊕yx,w,y⟩\Lambda|x,z,w,y\rangle=|x,z\oplus y_x,w,y\rangleΛ∣x,z,w,y⟩=∣x,z⊕yx​,w,y⟩ (the identity when yx=⊥y_x=\perpyx​=⊥). Let PPP be the projection onto the basis states whose oracle register contains a zero, i.e. yx∗=0my_{x^*}=0^myx∗​=0m for some x∗x^*x∗. Then for every state ∣ψ⟩|\psi\rangle∣ψ⟩ (unit vector),

∥P Φ∣ψ⟩∥  ≤  ∥P∣ψ⟩∥+22m.\|P\,\Phi|\psi\rangle\|\;\le\;\|P|\psi\rangle\|+\sqrt{\frac{2}{2^m}} .∥PΦ∣ψ⟩∥≤∥P∣ψ⟩∥+2m2​​.

Iterated over ttt queries, this shows that the compressed-oracle probability of having recorded a zero is O(t2/2m)O(t^2/2^m)O(t2/2m). Combined with Lemma 26.1, this gives the Ω(2m/2)\Omega(2^{m/2})Ω(2m/2) query lower bound for unstructured search against a random oracle.

Formalization Note ∣ψ⟩|\psi\rangle∣ψ⟩ is an arbitrary unit vector of the full space; no reachability condition is imposed, and the proof uses none. The register order is ((x,z),w),y((x,z),w),y((x,z),w),y, as in (26.30).

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Lemma 26.1: Standard vs compressed oracle: √p ≤ √p′ + √(k/2^m).
  • Fact 26.2: One compressed phase query raises ‖P ψ‖ by at most √(2/2^m).

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
3 thms2 active usersReviewed
🏆Completed
Captain: Goku

Quantum Computing: Lecture Notes 18: Perfect quantum bit commitment is impossible: perfect concealing implies no bindingTextbook

Motivation

This mission develops Chapter 18 of Quantum Computing: Lecture Notes by Ronald de Wolf. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

This is the impossibility of perfect quantum bit commitment.

In de Wolf's model, a quantum bit-commitment protocol is represented by the two joint pure states ∣ϕ0⟩,∣ϕ1⟩∈CdA⊗CdB|\phi_0\rangle, |\phi_1\rangle \in \mathbb{C}^{d_A} \otimes \mathbb{C}^{d_B}∣ϕ0​⟩,∣ϕ1​⟩∈CdA​⊗CdB​ that Alice and Bob share at the end of the commit phase when Alice honestly commits to the bit b=0b = 0b=0 or b=1b = 1b=1. Here CdA\mathbb{C}^{d_A}CdA​ is Alice's space and CdB\mathbb{C}^{d_B}CdB​ is Bob's. Restricting to pure states loses no generality. The protocol is perfectly concealing if Bob's reduced density matrix does not depend on bbb, where Tr⁡A\operatorname{Tr}_ATrA​ is the partial trace over Alice's register.

Let dA,dB≥0d_A, d_B \ge 0dA​,dB​≥0 and let ∣ϕ0⟩,∣ϕ1⟩|\phi_0\rangle, |\phi_1\rangle∣ϕ0​⟩,∣ϕ1​⟩ be unit vectors in CdA⊗CdB\mathbb{C}^{d_A} \otimes \mathbb{C}^{d_B}CdA​⊗CdB​ with

Tr⁡A(∣ϕ0⟩⟨ϕ0∣)=Tr⁡A(∣ϕ1⟩⟨ϕ1∣).\operatorname{Tr}_A\big(|\phi_0\rangle\langle\phi_0|\big) = \operatorname{Tr}_A\big(|\phi_1\rangle\langle\phi_1|\big).TrA​(∣ϕ0​⟩⟨ϕ0​∣)=TrA​(∣ϕ1​⟩⟨ϕ1​∣).

Then there is a unitary UUU on Alice's space CdA\mathbb{C}^{d_A}CdA​ such that

(U⊗I) ∣ϕ0⟩=∣ϕ1⟩.(U \otimes I)\,|\phi_0\rangle = |\phi_1\rangle.(U⊗I)∣ϕ0​⟩=∣ϕ1​⟩.

So after committing to 000, Alice can turn the joint state into the honest commitment to 111 by acting on her own register only. The protocol is then not binding at all, and no quantum bit-commitment protocol is both perfectly concealing and perfectly binding.

Formalization Note Alice's space is indexed by Fin dA and Bob's by Fin dB. The joint space is indexed by Fin dA × Fin dB with Alice's register first. Tr⁡A(∣ϕ⟩⟨ϕ∣)\operatorname{Tr}_A(|\phi\rangle\langle\phi|)TrA​(∣ϕ⟩⟨ϕ∣) is traceLeft (pureDensity phi), and U⊗IU \otimes IU⊗I is Matrix.kronecker U 1 with UUU in Matrix.unitaryGroup. When dA=0d_A = 0dA​=0 or dB=0d_B = 0dB​=0 no unit vector exists, so the statement is vacuous in those cases, as in the source.

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • §18.4: Perfect quantum bit commitment is impossible: perfect concealing implies no binding.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Ronald de Wolf, Quantum Computing: Lecture Notes, 27 August 2026. https://homepages.cwi.nl/~rdewolf/qcnotes.pdf.
1 thm2 active usersReviewed
🏆Completed
Captain: Goku

Lecture Notes on Quantum Algorithms 24: Effective spectral gap lemmaTextbook

Motivation

This mission develops Chapter 24 of Lecture Notes on Quantum Algorithms by Andrew Childs. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

Let H\mathcal HH be a finite-dimensional complex inner product space and let Π\PiΠ and Δ\DeltaΔ be orthogonal projections on H\mathcal HH. Put U=(2Π−I)(2Δ−I)U = (2\Pi - I)(2\Delta - I)U=(2Π−I)(2Δ−I), a unitary. For ω≥0\omega \ge 0ω≥0, let PωP_\omegaPω​ be the orthogonal projector onto the span of the eigenvectors of UUU whose eigenvalues are eiθe^{i\theta}eiθ with ∣θ∣<ω|\theta| < \omega∣θ∣<ω. If ∣ϕ⟩∈H|\phi\rangle \in \mathcal H∣ϕ⟩∈H is a unit vector with Δ∣ϕ⟩=0\Delta|\phi\rangle = 0Δ∣ϕ⟩=0, then

∥PωΠ∣ϕ⟩∥≤ω2.\|P_\omega \Pi |\phi\rangle\| \le \frac{\omega}{2}.∥Pω​Π∣ϕ⟩∥≤2ω​.

The lemma is the key step in showing that the phase-estimation algorithm built from a dual adversary solution rejects inputs with f(x)=0f(x) = 0f(x)=0: the start state has little weight on eigenvectors of small phase.

Formalization Note Π\PiΠ and Δ\DeltaΔ are given as the orthogonal projections onto submodules KKK and LLL; 2P−I2P - I2P−I is Mathlib's reflection through the subspace. The eigenvalue set is {eiθ:θ∈R,∣θ∣<ω}\{e^{i\theta} : \theta \in \mathbb{R}, |\theta| < \omega\}{eiθ:θ∈R,∣θ∣<ω}, which is the same set as with θ\thetaθ restricted to (−π,π](-\pi,\pi](−π,π].

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • Theorem 24.1: The dual of the general adversary bound.
  • Theorem 24.2: Adversary bound of a composed function is at most the product.
  • Lemma 24.3: Effective spectral gap lemma.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Andrew Childs, Lecture Notes on Quantum Algorithms, 17 April 2025. https://www.cs.umd.edu/~amchilds/qa/qa.pdf.
3 thms2 active usersReviewed
🏆Completed
Captain: Goku

Quantum Computing: Lecture Notes 12: |t⁻_x⟩ is close to the large-eigenphase part of U_xTextbook

Motivation

This mission develops Chapter 12 of Quantum Computing: Lecture Notes by Ronald de Wolf. The chapter's selected results describe the mathematical guarantees below.

Setting and conventions

The linked theorem statements and shared definitions specify the register types, oracle access, normalization, and resource conventions. Their source locators and independent readbacks accompany each target.

Chapter goal

In the setting of Claim 3 (a feasible solution {uxj,vxj}\{u_{xj}, v_{xj}\}{uxj​,vxj​} of the dual SDP for ADV±(f)\mathrm{ADV}^{\pm}(f)ADV±(f) with objective value A>0A > 0A>0, the states ∣tx±⟩|t^\pm_x\rangle∣tx±​⟩, ∣ψy⟩|\psi_y\rangle∣ψy​⟩, the projectors Λ\LambdaΛ, Πx\Pi_xΠx​ and the unitary Ux=(2Πx−I)(2Λ−I)U_x = (2\Pi_x - I)(2\Lambda - I)Ux​=(2Πx​−I)(2Λ−I)), let Θ=1/(1000A)\Theta = 1/(1000A)Θ=1/(1000A). For every x∈Dx \in Dx∈D there is a vector ∣φ⟩|\varphi\rangle∣φ⟩ in the span of the eigenvectors of UxU_xUx​ whose eigenvalues eiθe^{i\theta}eiθ, θ∈(−π,π]\theta \in (-\pi, \pi]θ∈(−π,π], satisfy ∣θ∣>Θ|\theta| > \Theta∣θ∣>Θ, such that

∥∣tx−⟩−∣φ⟩∥≤120+12000A.\big\| |t^-_x\rangle - |\varphi\rangle \big\| \le \frac{1}{20} + \frac{1}{2000A}.​∣tx−​⟩−∣φ⟩​≤201​+2000A1​.

Equivalently, the component of ∣tx−⟩|t^-_x\rangle∣tx−​⟩ in the eigenspaces of UxU_xUx​ with eigenphase ∣θ∣≤Θ|\theta| \le \Theta∣θ∣≤Θ has norm at most 1/20+1/(2000A)1/20 + 1/(2000A)1/20+1/(2000A). Together with Claim 3, this lets phase estimation with precision Θ/2\Theta/2Θ/2 distinguish ∣tx+⟩|t^+_x\rangle∣tx+​⟩ from ∣tx−⟩|t^-_x\rangle∣tx−​⟩.

Formalization Note The eigenphase θ\thetaθ is the principal argument Complex.arg, which takes values in (−π,π](-\pi, \pi](−π,π]; the span of eigenvectors is the supremum of the eigenspaces of UxU_xUx​ over the eigenvalues with ∣arg⁡μ∣>Θ|\arg \mu| > \Theta∣argμ∣>Θ. The informal "close" is pinned to the bound the source's proof establishes, and it is stated for every A>0A > 0A>0 (the proof's own hypothesis; the source adds "we may assume A≥1A \ge 1A≥1" only to call the bound small).

Significance

These are known mathematical results. Their Lean statements have been checked and reviewed against the cited source; formal proofs remain open for contributors.

Supporting targets

  • §12.1, Claim 1: Generalized adversary method: the final progress measure is small.
  • §12.1, Claim 2: Generalized adversary method: one query changes the progress measure by at most 2 max_i ‖Γ_i‖.
  • §12.1, Theorem 2: Generalized (negative-weights) adversary lower bound.
  • §12.3, unnumbered: A dual-adversary solution of value A gives an O(A)-query bounded-error algorithm.
  • §12.3, Claim 3: |t⁺_x⟩ is within 0.01 of an eigenvalue-1 eigenvector of U_x.
  • §12.3, Claim 4: |t⁻_x⟩ is close to the large-eigenphase part of U_x.

Formalization scope

Prove the existing statements under their stated hypotheses and definitions. Any correction to a printed source claim remains explicit in its target description. All listed target proofs are currently sorry; no proof is claimed. The mission goal is the designated chapter capstone; the remaining targets are separate milestones. Parts a, b, c and d are merged into their original chapter.

Selected references

  • Ronald de Wolf, Quantum Computing: Lecture Notes, 27 August 2026. https://homepages.cwi.nl/~rdewolf/qcnotes.pdf.
6 thms2 active usersReviewed
🏆Completed
Condensed Matter PhysicsLinear algebraMathematical Physics+2·Captain: Alien60

Seiberg–Seifnashri–Shao 2024: Non-invertible Kramers–Wannier symmetry on the Ising chainResearch Paper

Motivation

Kramers–Wannier duality was found in 1941 as a map between the high- and low-temperature phases of the two-dimensional Ising model (Kramers–Wannier 1941). In its quantum form it exchanges the transverse field XjX_jXj​ and the nearest-neighbour coupling Zj−1ZjZ_{j-1}Z_jZj−1​Zj​ of the transverse-field Ising chain, and it fixes the critical point. In conformal field theory it is implemented by a topological defect line (Fröhlich–Fuchs–Runkel–Schweigert 2004), the prototypical non-invertible symmetry: its fusion with itself is a sum of operators, not a single invertible one (Shao, TASI lectures 2023).

On a finite periodic chain the duality cannot be a unitary operator, because it maps the Z2\mathbb Z_2Z2​-even sector to itself and annihilates the odd sector. Seiberg–Shao 2023 and Seiberg–Seifnashri–Shao 2024 wrote the exact lattice operator D\mathsf DD on the tensor-product Hilbert space (C2)⊗L(\mathbb C^2)^{\otimes L}(C2)⊗L. They showed that its algebra closes only after the lattice translation is included: D2=(1+η)T−1\mathsf D^2=(1+\eta)T^{-1}D2=(1+η)T−1. This is a lattice relation that a fusion category does not describe, since the continuum fusion rule N2=1+η\mathcal N^2=1+\etaN2=1+η appears only at low energy. The same identity drives the LSM-type constraint proved in that paper.

Timeline. 1941: Kramers and Wannier find the duality of the 2D Ising model. 2004: Fröhlich, Fuchs, Runkel and Schweigert realise it as a conformal defect. 2016: Aasen–Mong–Fendley construct the lattice defects. 2023–2024: Seiberg, Seifnashri and Shao find the exact operator algebra on the tensor-product Hilbert space, including the translation factor.

Setting

Fix L≥1L\ge1L≥1 qubits on a ring, with site indices in Z/L\mathbb Z/LZ/L (Fin L with modular arithmetic). A computational-basis state is labelled by a configuration s=(s0,…,sL−1)∈(Z/2)Ls=(s_0,\dots,s_{L-1})\in(\mathbb Z/2)^Ls=(s0​,…,sL−1​)∈(Z/2)L (Config L). Operators are complex 2L×2L2^L\times2^L2L×2L matrices indexed by configurations (Op L). Write sgn⁡(a)=(−1)a\operatorname{sgn}(a)=(-1)^asgn(a)=(−1)a.

  • Pauli operators. Xj∣s⟩=∣s+ej⟩X_j|s\rangle=|s+e_j\rangleXj​∣s⟩=∣s+ej​⟩ flips spin jjj (pauliX), and Zj∣s⟩=(−1)sj∣s⟩Z_j|s\rangle=(-1)^{s_j}|s\rangleZj​∣s⟩=(−1)sj​∣s⟩ (pauliZ).
  • Z2\mathbb Z_2Z2​ symmetry. η=∏jXj\eta=\prod_j X_jη=∏j​Xj​ flips every spin (eta).
  • Translation. T∣s⟩=∣s′⟩T|s\rangle=|s'\rangleT∣s⟩=∣s′⟩ with sk′=sk−1s'_k=s_{k-1}sk′​=sk−1​ (transl), so that TXjT−1=Xj+1TX_jT^{-1}=X_{j+1}TXj​T−1=Xj+1​.
  • Kramers–Wannier operator (kwD):
⟨t∣ D ∣s⟩  =  2−L/2 (−1)∑ktk (sk+sk+1).\langle t|\,\mathsf D\,|s\rangle \;=\; 2^{-L/2}\,(-1)^{\sum_{k} t_k\,(s_k+s_{k+1})}.⟨t∣D∣s⟩=2−L/2(−1)∑k​tk​(sk​+sk+1​).

This explicit formula is the definition used throughout. Relations (1.5) together with Dη=D\mathsf D\eta=\mathsf DDη=D determine D\mathsf DD up to a scalar, and (1.4) fixes the scalar up to a sign. Every statement below is invariant under D↦−D\mathsf D\mapsto-\mathsf DD↦−D, so this is the operator of Seiberg–Seifnashri–Shao in that paper's normalisation (footnote 1). Its equality with their matrix-product presentation (2.64) is not formalised here.

Formalization targets

Goal: the non-invertible fusion rule (eqs. (1.4), (2.65), (F.18))

D2=(1+η) T−1for every L≥1.\mathsf D^2=(1+\eta)\,T^{-1}\qquad\text{for every }L\ge1 .D2=(1+η)T−1for every L≥1.

Milestones

  1. Invertible algebra (eqs. (1.1)–(1.3), (2.10), (2.65)): η2=1\eta^2=1η2=1, Tη=ηTT\eta=\eta TTη=ηT, TL=1T^L=1TL=1, TXjT−1=Xj+1TX_jT^{-1}=X_{j+1}TXj​T−1=Xj+1​, TZjT−1=Zj+1TZ_jT^{-1}=Z_{j+1}TZj​T−1=Zj+1​, ηXj=Xjη\eta X_j=X_j\etaηXj​=Xj​η, ηZj=−Zjη\eta Z_j=-Z_j\etaηZj​=−Zj​η.
  2. Duality on XXX (eq. (1.5)): DXj=Zj−1Zj D\mathsf D X_j=Z_{j-1}Z_j\,\mathsf DDXj​=Zj−1​Zj​D.
  3. Duality on ZZZZZZ (eq. (1.5)): DZj−1Zj=Xj−1 D\mathsf D Z_{j-1}Z_j=X_{j-1}\,\mathsf DDZj−1​Zj​=Xj−1​D.
  4. Absorbing η\etaη (eqs. (1.4), (2.65)): Dη=ηD=D\mathsf D\eta=\eta\mathsf D=\mathsf DDη=ηD=D.
  5. Translation and adjoint (eq. (2.65), Table 1): TD=DTT\mathsf D=\mathsf DTTD=DT and D†=DT\mathsf D^\dagger=\mathsf DTD†=DT.
  6. Non-invertibility (§1.1): D\mathsf DD is not a unit, and ker⁡D≠0\ker\mathsf D\neq0kerD=0.

Significance

The result. The goal identity is the defining relation of the lattice Kramers–Wannier symmetry. Together with milestones 4–5 it shows that D/2\mathsf D/\sqrt2D/2​ restricted to the η\etaη-even sector squares to a translation, so D\mathsf DD is a square root of T−1T^{-1}T−1 there. The paper uses this to derive its LSM-type constraint, and also to show that the lattice symmetry is not a fusion category, because the continuum rule N2=1+η\mathcal N^2=1+\etaN2=1+η only emerges at low energy. Milestones 2–3 are the lattice Kramers–Wannier map on Z2\mathbb Z_2Z2​-invariant operators, so D\mathsf DD commutes with the critical Ising Hamiltonian H=−∑j(Zj−1Zj+Xj)H=-\sum_j(Z_{j-1}Z_j+X_j)H=−∑j​(Zj−1​Zj​+Xj​).

Formalizing it. The identities are proved in the physics literature by tensor-network manipulations ((F.14)–(F.18)) and checked numerically, but none of this is machine-checked. This mission produces a reusable Lean setup for qubit chains: Pauli strings on (Z/2)L(\mathbb Z/2)^L(Z/2)L, Z2\mathbb Z_2Z2​ and translation symmetries, and an exact non-invertible operator with its full operator algebra for every chain length. No mission on the platform or file in Mathlib currently formalizes a non-invertible lattice symmetry.

Difficulty

Each identity is easy to check for a fixed small LLL by brute force. The content is that it holds for all LLL at once, on a 2L2^L2L-dimensional space. The difficulty is bookkeeping with periodic indices: sums over (Z/2)L(\mathbb Z/2)^L(Z/2)L of characters (−1)⟨t, ⋅⟩(-1)^{\langle t,\,\cdot\rangle}(−1)⟨t,⋅⟩, reindexing k↦k±1k\mapsto k\pm1k↦k±1 in Fin L (including L=1L=1L=1, where k+1=kk+1=kk+1=k, and L=2L=2L=2, where k+1=k−1k+1=k-1k+1=k−1), and the exact normalisation 2−L/22^{-L/2}2−L/2, so that the right-hand side carries no stray factor 2L2^L2L. The goal also involves the matrix inverse T−1T^{-1}T−1, so a proof must establish that TTT is invertible rather than rely on Lean's junk value.

Formalization scope

  • Representation. Operators are Matrix (Fin L → ZMod 2) (Fin L → ZMod 2) ℂ; products are matrix products; T−1T^{-1}T−1 is Mathlib's matrix inverse (transl L)⁻¹; †\dagger† is conjTranspose. Every theorem assumes [NeZero L], i.e. L≥1L\ge1L≥1.
  • Conventions committed to. Zj=(−1)sjZ_j=(-1)^{s_j}Zj​=(−1)sj​ on basis states, translation to the right (sk′=sk−1s'_k=s_{k-1}sk′​=sk−1​), and D\mathsf DD normalised by 2−L/22^{-L/2}2−L/2 as in eqs. (1.4)/(2.65). Other conventions (left translation, unnormalised D\mathsf DD) change the statements to TTT instead of T−1T^{-1}T−1 or add a factor 2L2^L2L; the conventions here reproduce the paper's equations verbatim.
  • No trivialisation. kwD is an explicit formula, not a hypothesis. The goal fixes both sides completely, and milestone 1 pins the translation convention against eq. (2.10).
  • Out of scope. The defect Hamiltonians, the MPO presentation (2.64) itself, the Z2\mathbb Z_2Z2​- and duality-twisted operators of §§2.3.3–2.3.4, and the LSM constraint of §4.
  • Reusable. Pauli operators, translation and η\etaη on qubit chains; character-sum lemmas over (Z/2)L(\mathbb Z/2)^L(Z/2)L. Contributions of general Mathlib lemmas on Walsh–Hadamard-type sums are welcome.

Selected references

  • H. A. Kramers, G. H. Wannier, Statistics of the Two-Dimensional Ferromagnet. Part I, Phys. Rev. 60, 252–262 (1941). doi:10.1103/PhysRev.60.252
  • J. Fröhlich, J. Fuchs, I. Runkel, C. Schweigert, Kramers–Wannier duality from conformal defects, Phys. Rev. Lett. 93, 070601 (2004). arXiv:cond-mat/0404051
  • D. Aasen, R. S. K. Mong, P. Fendley, Topological Defects on the Lattice I: The Ising model, J. Phys. A 49, 354001 (2016). arXiv:1601.07185
  • N. Seiberg, S.-H. Shao, Majorana chain and Ising model – (non-invertible) translations, anomalies, and emanant symmetries, SciPost Phys. 16, 064 (2024). arXiv:2307.02534
  • S.-H. Shao, What's Done Cannot Be Undone: TASI Lectures on Non-Invertible Symmetries (2023). arXiv:2308.00747
  • Source: N. Seiberg, S. Seifnashri, S.-H. Shao, Non-invertible symmetries and LSM-type constraints on a tensor product Hilbert space, SciPost Phys. 16, 154 (2024). arXiv:2401.12281
8 thms2 active usersReviewed
🏆Completed
Mathematical Physics·Captain: Lucas

Projective Measurements Are Equilibration Processes (finite dimension)Research Paper

Motivation

Projective measurement is usually taught together with state update: an observable is measured, a value is obtained, and the conditional state is assigned. These are two different questions. The first is the prediction of the outgoing ensemble from the incoming ensemble and the measured observable alone, before the outcome is known; the second is the incorporation of the outcome once it is known. The technical brief by G. Carcassi for the Assumptions of Physics collaboration (AoP technical briefs) isolates the first question and treats the measurement as a black-box process that only relates incoming and outgoing ensembles, without committing to a mechanism.

Within this view the standard prediction, the nonselective (Lüders) post-measurement ensemble, admits several characterizations that come from different traditions: the algebraic pinching map, the long-time limit of a dephasing Lindblad equation (open quantum systems, see Lindblad 1976), and a closest-point principle, either in Hilbert–Schmidt distance or in quantum relative entropy (Umegaki 1962). The brief's Theorem 1 states that these characterizations select the same process.

Setting

Fix a finite index set nnn; the Hilbert space is Cn\mathbb C^nCn and operators are complex n×nn\times nn×n matrices. A density matrix is a positive semidefinite matrix ρ\rhoρ with tr⁡ρ=1\operatorname{tr}\rho=1trρ=1. An observable is a Hermitian matrix XXX with spectral decomposition

X=∑xx Px,X=\sum_{x}x\,P_x ,X=x∑​xPx​,

where xxx ranges over the distinct eigenvalues of XXX and PxP_xPx​ is the orthogonal projection onto the xxx-eigenspace.

  • The pinching of ρ\rhoρ is PX(ρ)=∑xPxρPx\mathcal P_X(\rho)=\sum_x P_x\rho P_xPX​(ρ)=∑x​Px​ρPx​.
  • The dephasing Lindblad equation with vanishing Hamiltonian and jump operator XXX is dρdt=XρX−12{X2,ρ}\dfrac{d\rho}{dt}=X\rho X-\tfrac12\{X^2,\rho\}dtdρ​=XρX−21​{X2,ρ}, t≥0t\ge0t≥0.
  • The Hilbert–Schmidt distance is ∥ρ−σ∥HS=(tr⁡(ρ−σ)∗(ρ−σ))1/2\|\rho-\sigma\|_{HS}=\big(\operatorname{tr}(\rho-\sigma)^*(\rho-\sigma)\big)^{1/2}∥ρ−σ∥HS​=(tr(ρ−σ)∗(ρ−σ))1/2.
  • The relative entropy is D(ρ∥σ)=tr⁡(ρlog⁡ρ)−tr⁡(ρlog⁡σ)D(\rho\|\sigma)=\operatorname{tr}(\rho\log\rho)-\operatorname{tr}(\rho\log\sigma)D(ρ∥σ)=tr(ρlogρ)−tr(ρlogσ) when supp⁡ρ⊆supp⁡σ\operatorname{supp}\rho\subseteq\operatorname{supp}\sigmasuppρ⊆suppσ and +∞+\infty+∞ otherwise; the von Neumann entropy is S(ρ)=−tr⁡(ρlog⁡ρ)S(\rho)=-\operatorname{tr}(\rho\log\rho)S(ρ)=−tr(ρlogρ).

The spectral equilibration process ΦX\Phi_XΦX​ is an unknown map on density matrices. The brief's four conditions on it are:

  • EQ-MEAS: ΦX(ρ)=∑xPxρPx\Phi_X(\rho)=\sum_xP_x\rho P_xΦX​(ρ)=∑x​Px​ρPx​.
  • EQ-LIND: ΦX(ρ)=lim⁡t→∞ρ(t)\Phi_X(\rho)=\lim_{t\to\infty}\rho(t)ΦX​(ρ)=limt→∞​ρ(t) for the Lindblad evolution with ρ(0)=ρ\rho(0)=\rhoρ(0)=ρ.
  • EQ-HS: ΦX(ρ)=argmin⁡{∥ρ−σ∥HS2: σ a density matrix, [σ,X]=0}\Phi_X(\rho)=\operatorname{argmin}\{\|\rho-\sigma\|_{HS}^2:\ \sigma \text{ a density matrix},\ [\sigma,X]=0\}ΦX​(ρ)=argmin{∥ρ−σ∥HS2​: σ a density matrix, [σ,X]=0}.
  • EQ-KL: ΦX(ρ)=argmin⁡{D(ρ∥σ): σ a density matrix, [σ,X]=0}\Phi_X(\rho)=\operatorname{argmin}\{D(\rho\|\sigma):\ \sigma \text{ a density matrix},\ [\sigma,X]=0\}ΦX​(ρ)=argmin{D(ρ∥σ): σ a density matrix, [σ,X]=0}.

Formalization targets

Goal: Theorem 1 (finite-dimensional form)

For every Hermitian XXX and every map Φ\PhiΦ,

EQ-MEAS  ⟺  EQ-LIND,EQ-MEAS  ⟺  EQ-HS,EQ-MEAS  ⟺  EQ-KL.\text{EQ-MEAS}\iff\text{EQ-LIND},\qquad \text{EQ-MEAS}\iff\text{EQ-HS},\qquad \text{EQ-MEAS}\iff\text{EQ-KL}.EQ-MEAS⟺EQ-LIND,EQ-MEAS⟺EQ-HS,EQ-MEAS⟺EQ-KL.

Milestones

The milestones follow the steps of the brief's proof: the block decomposition ρ=∑x,yPxρPy\rho=\sum_{x,y}P_x\rho P_yρ=∑x,y​Px​ρPy​ (eq. 2); the blockwise solution Pxρ(t)Py=e−t2(x−y)2Pxρ(0)PyP_x\rho(t)P_y=e^{-\frac t2(x-y)^2}P_x\rho(0)P_yPx​ρ(t)Py​=e−2t​(x−y)2Px​ρ(0)Py​ (eqs. 4–5); convergence ρ(t)→PX(ρ)\rho(t)\to\mathcal P_X(\rho)ρ(t)→PX​(ρ) (eqs. 9–10) and EQ-LIND ⇔\Leftrightarrow⇔ EQ-MEAS; the equivalence [σ,X]=0  ⟺  PX(σ)=σ[\sigma,X]=0\iff\mathcal P_X(\sigma)=\sigma[σ,X]=0⟺PX​(σ)=σ; PX\mathcal P_XPX​ preserves density matrices; the Pythagorean identity (eq. 11) and EQ-HS ⇔\Leftrightarrow⇔ EQ-MEAS; the relative-entropy chain rule (eq. 12); D(ρ∥PXρ)=S(PXρ)−S(ρ)D(\rho\|\mathcal P_X\rho)=S(\mathcal P_X\rho)-S(\rho)D(ρ∥PX​ρ)=S(PX​ρ)−S(ρ) (eq. 13); non-negativity and definiteness of DDD; and EQ-KL ⇔\Leftrightarrow⇔ EQ-MEAS.

Significance

The result says that "projective measurement" can be specified without reference to wave-function collapse: any of the three dynamical or variational descriptions pins down the same input–output map. Downstream, this justifies treating repeatable measurements as equilibration processes whose outputs are fixed points, and treating the outcome-update step as classical conditioning on a mixture.

The mathematics in finite dimension is classical; this mission asks for a machine-checked proof of the equivalence in Lean 4 with Mathlib. A complete development produces reusable pieces: spectral projections and the pinching map for Hermitian matrices, explicit solutions and asymptotics of a dephasing Lindblad semigroup, and an extended-real relative entropy with Klein's inequality.

Difficulty

The Hilbert–Schmidt part is linear algebra once block-diagonality is in place. The Lindblad part requires uniqueness of solutions of a linear matrix ODE on the closed half-line, to transfer the explicit blockwise formula to an arbitrary solution. The relative-entropy part is the most delicate: it needs the interaction between the functional-calculus logarithm and block structure (log⁡σ\log\sigmalogσ is block diagonal when σ\sigmaσ commutes with XXX), the support conventions that make DDD take the value +∞+\infty+∞, and Klein's inequality with its equality case for matrices, which the development must supply if the library does not already provide it.

Formalization scope

  • Finite dimension only. The Hilbert space is Cn\mathbb C^nCn for a finite type nnn, so all observables are bounded with pure point spectrum and every density matrix has finite entropy. The brief also covers separable infinite-dimensional spaces and unbounded observables (where Condition [B]FINENT is a genuine restriction); that regime is not part of this mission.
  • Spectral data. PxP_xPx​ is built from Mathlib's eigen-decomposition of a Hermitian matrix and is defined for every real xxx (zero when xxx is not an eigenvalue); the pinching sums over the finite set of distinct eigenvalues.
  • Lindblad solutions. A solution is any function R→Mn(C)\mathbb R\to M_n(\mathbb C)R→Mn​(C) with the given initial value whose derivative within [0,∞)[0,\infty)[0,∞) satisfies the equation for t≥0t\ge0t≥0. EQ-LIND requires existence of a solution and convergence of every solution, so it cannot hold vacuously.
  • Argmin. "Φ(ρ)=argmin⁡\Phi(\rho)=\operatorname{argmin}Φ(ρ)=argmin" is encoded as "Φ(ρ)\Phi(\rho)Φ(ρ) is a feasible minimizer"; uniqueness is a consequence proved by the theorem, not an assumption.
  • Logarithm and relative entropy. log⁡\loglog is the continuous functional calculus applied to Real.log (so log⁡0=0\log0=0log0=0, i.e. the logarithm on the support); DDD takes values in the extended reals and equals +∞+\infty+∞ unless ker⁡σ⊆ker⁡ρ\ker\sigma\subseteq\ker\rhokerσ⊆kerρ.
  • Convergence is in the entrywise topology, which coincides with the trace-norm and Hilbert–Schmidt topologies in finite dimension.

Contributions of general-purpose lemmas (Klein's inequality, uniqueness for linear ODEs on half-lines, properties of pinching) are welcome as separate theorems.

Selected references

  • G. Carcassi (Assumptions of Physics collaboration), Projective measurements are equilibration processes, AoP Technical Brief 017. https://assumptionsofphysics.org/resources/briefs
  • G. Lindblad, On the generators of quantum dynamical semigroups, Comm. Math. Phys. 48 (1976) 119–130. https://doi.org/10.1007/BF01608499
  • H. Umegaki, Conditional expectation in an operator algebra. IV. Entropy and information, Kodai Math. Sem. Rep. 14 (1962) 59–85. https://doi.org/10.2996/kmj/1138844604
15 thms2 active usersReviewed
🏆Completed
Complexity Theory·Captain: marwahaha

Regular trajectories, pruning and quantum parityOpen Problem

Motivation

The parity function tests whether a shallow quantum circuit can combine information from all input bits. Polynomially many ancillary qubits make the measured-output formulation substantially stronger than a restriction to clean final registers. The pinned manuscript supplies the research context.

Setting

Allowed gates are arbitrary one-qubit unitaries and Toffoli gates with finitely many controls. Gates within a layer have disjoint supports. Inputs occupy the first n qubits, with zero ancillas.

Formalization target

The selected goal is OAI.QAC.parity_lower_bound_polynomial_size. Its central assertion is

∃x∈{0,1}n:P(output=parity(x))<2/3.\exists x\in\{0,1\}^n:\quad P(\mathrm{output}=\mathrm{parity}(x))<2/3.∃x∈{0,1}n:P(output=parity(x))<2/3.

The theorem states that, for any positive integers d, k and C, there is a threshold n₀ such that for every n ≥ n₀ and every total qubit count N with n ≤ N ≤ C(n+1)^k, the following holds. Take any circuit given as a list of at most d physical layers on N qubits, where a layer is a list of gates with pairwise disjoint supports and each gate is either a one-qubit unitary on a single qubit or a Toffoli gate with an arbitrary finite set of control qubits and a target outside that set (flipping the target exactly when all controls are 1). The circuit's matrix is the product of the layer matrices, with the first layer acting first. For every choice of output qubit out among the N qubits, there exists an n-bit input x such that successProbability, the total Born probability over all final computational-basis strings y whose out-th bit equals the parity (sum mod 2) of x, is strictly less than 2/3. Here the input state is x placed on the first n qubits with all remaining qubits set to 0, and the remaining qubits are summed over with no requirement that they be clean. So fixed-depth circuits with polynomially many qubits cannot compute parity with worst-case success probability at least 2/3. The theorem is admitted in the source (proof is sorry).

Significance and status

The published goal fixes success threshold 2/3 and positive integer depth, polynomial exponent and coefficient. It does not state the manuscript's full range of positive advantages. The target is currently Open on Prove2Me. The manuscript's mathematical argument and a machine-checked proof of the selected statement are separate deliverables.

Difficulty

The bound covers arbitrary gates, all choices of output qubit and polynomially many total qubits. Discarded registers can retain unrestricted garbage.

Formalization scope

The exact published goal, its hypotheses and its referenced definition blocks specify the requested formalization. The explanatory formula above is a summary; all quantifiers and additional clauses in the linked statement remain required.

The available published material supplies this goal and its necessary definitions. Additional manuscript lemmas are not represented as attached milestones.

Selected references

  • OpenAI, Regular trajectories, pruning and quantum parity, preprint, 2026. Manuscript.
  • OpenAI, accompanying formal statements, commit adc7f1241b42. Selected goal source.
2 thms2 active usersReviewed
🏆Completed
Information Theory·Captain: wurtle

Classical capacity and entropy inequalities for generalized amplitude dampingResearch Paper

Motivation: the classical capacity of a thermal qubit channel

A quantum channel models how a physical system, used to carry information, is corrupted by its environment. The classical capacity of a channel is the largest rate, in bits per use, at which classical messages can be sent with vanishing error when the sender may encode into arbitrary — possibly entangled — states of many uses and the receiver may measure all outputs jointly. The Holevo–Schumacher–Westmoreland coding theorem (Holevo 1998; Schumacher–Westmoreland 1997) expresses this capacity as a regularized limit of the Holevo information over ever larger blocks of uses. Since Holevo information is not additive in general (Hastings 2009), computing the capacity of any specific channel requires a channel-specific additivity argument, and exact formulas are known only for a few families.

Generalized amplitude damping is the standard model of a qubit relaxing toward a thermal state: it loses or gains an excitation with probabilities set by a damping strength γ\gammaγ and a stationary excited population ν\nuν, and its coherences decay. It is one of the basic noise models of quantum information theory.

Background timeline

  • 1996–1998 — Coding theorems for pure and mixed signals establish the regularized capacity formula (Hausladen et al. 1996; Schumacher–Westmoreland 1997; Holevo 1998).
  • 2002 — Additivity for unital qubit channels (King 2002) and for entanglement-breaking channels (Shor 2002) settles some parameter regions of the family; Bennett, Shor, Smolin and Thapliyal describe the two-signal optimization for ordinary amplitude damping (BSST 2002).
  • 2005 — Giovannetti and Fazio derive the optimized one-use expression for ordinary amplitude damping (Giovannetti–Fazio 2005); Berry studies qubit channels achieving capacity with two states (Berry 2005).
  • 2007 — Hou and Fang characterize the one-use Holevo optimum of the generalized channel (Hou–Fang 2007).
  • 2009 — Hastings shows that Holevo information is not additive in general (Hastings 2009).
  • 2018 — Leditzky, Kaur, Datta and Wilde note that the unrestricted classical capacity of amplitude damping had not been determined (LKDW 2018).
  • 2020 — Khatri, Sharma and Wilde collect upper bounds for the generalized family (KSW 2020).
  • 2026 — Tang, Zhu, Bai and Wang prove additivity for qubit channels with a pure output, including ordinary amplitude damping (arXiv:2609.28592); Fang gives a close converse bound for the generalized family (arXiv:2609.14621).
  • 2026 — An OpenAI preprint, Classical capacity and entropy inequalities for generalized amplitude damping (OpenAI Math Release, September 24, 2026), claims the capacity for all (γ,ν)∈[0,1]2(\gamma,\nu)\in[0,1]^2(γ,ν)∈[0,1]2. The preprint has not been peer reviewed, and its results are not formally verified.

Setting

For γ,ν∈[0,1]\gamma,\nu\in[0,1]γ,ν∈[0,1] the generalized amplitude-damping channel acts on a qubit density matrix by

Aγ,ν(1−qzzˉq)=(1−t1−γ z1−γ zˉt),t=(1−γ)q+γν,\mathcal A_{\gamma,\nu}\begin{pmatrix}1-q&z\\\bar z&q\end{pmatrix}=\begin{pmatrix}1-t&\sqrt{1-\gamma}\,z\\\sqrt{1-\gamma}\,\bar z&t\end{pmatrix},\qquad t=(1-\gamma)q+\gamma\nu,Aγ,ν​(1−qzˉ​zq​)=(1−t1−γ​zˉ​1−γ​zt​),t=(1−γ)q+γν,

and nnn independent uses act as Aγ,ν⊗n\mathcal A_{\gamma,\nu}^{\otimes n}Aγ,ν⊗n​ on 2n×2n2^n\times2^n2n×2n matrices. With S(P)=−Tr⁡Pln⁡PS(P)=-\operatorname{Tr}P\ln PS(P)=−TrPlnP, the Holevo information in bits is

χ(N)=1ln⁡2sup⁡{qa,ρa}{S(∑aqaN(ρa))−∑aqaS(N(ρa))},\chi(\mathcal N)=\frac1{\ln2}\sup_{\{q_a,\rho_a\}}\Bigl\{S\Bigl(\sum_aq_a\mathcal N(\rho_a)\Bigr)-\sum_aq_aS(\mathcal N(\rho_a))\Bigr\},χ(N)=ln21​{qa​,ρa​}sup​{S(a∑​qa​N(ρa​))−a∑​qa​S(N(ρa​))},

over all finite ensembles of input states (entangled across uses when N=A⊗n\mathcal N=\mathcal A^{\otimes n}N=A⊗n). The classical capacity C(N)C(\mathcal N)C(N) is the supremum of rates RRR for which there are nnn-use codes — arbitrary input states for each of MnM_nMn​ equiprobable messages and one joint measurement of all outputs — with average error tending to zero and lim inf⁡nn−1log⁡2Mn≥R\liminf_n n^{-1}\log_2M_n\ge Rliminfn​n−1log2​Mn​≥R; there is no shared entanglement or feedback. Write hhh for binary entropy in nats, g(u)=h((1+1−4u)/2)g(u)=h\bigl((1+\sqrt{1-4u})/2\bigr)g(u)=h((1+1−4u​)/2) and vγ,ν(p)=γν(1−ν)+γ(1−γ)(p−ν)2v_{\gamma,\nu}(p)=\gamma\nu(1-\nu)+\gamma(1-\gamma)(p-\nu)^2vγ,ν​(p)=γν(1−ν)+γ(1−γ)(p−ν)2.

Formalization targets

Goal: additivity and capacity formula (Theorem 1.1)

For all γ,ν∈[0,1]\gamma,\nu\in[0,1]γ,ν∈[0,1] and every n≥1n\ge1n≥1,

χ(Aγ,ν⊗n)=n χ(Aγ,ν)=nln⁡2max⁡0≤p≤1{h((1−γ)p+γν)−g(vγ,ν(p))},\chi(\mathcal A_{\gamma,\nu}^{\otimes n})=n\,\chi(\mathcal A_{\gamma,\nu})=\frac n{\ln 2}\max_{0\le p\le1}\bigl\{h((1-\gamma)p+\gamma\nu)-g(v_{\gamma,\nu}(p))\bigr\},χ(Aγ,ν⊗n​)=nχ(Aγ,ν​)=ln2n​0≤p≤1max​{h((1−γ)p+γν)−g(vγ,ν​(p))},

the maximum is attained, C(Aγ,ν)=χ(Aγ,ν)C(\mathcal A_{\gamma,\nu})=\chi(\mathcal A_{\gamma,\nu})C(Aγ,ν​)=χ(Aγ,ν​), and for every maximizer ppp the equiprobable independent products of the two signals ϕ±=1−p ∣0⟩±p ∣1⟩\phi_\pm=\sqrt{1-p}\,|0\rangle\pm\sqrt p\,|1\rangleϕ±​=1−p​∣0⟩±p​∣1⟩ attain χ(Aγ,ν⊗n)\chi(\mathcal A_{\gamma,\nu}^{\otimes n})χ(Aγ,ν⊗n​) at every block length. The goal is published on the platform with status Open.

Significance

The result itself. The theorem gives a closed one-variable formula for the classical capacity of every generalized amplitude-damping channel, removing the regularization over block lengths: entangled input signals across uses give no advantage, and a simple product code with two signals per use is optimal. This settles the ordinary amplitude-damping case listed as open by LKDW 2018 and extends beyond the pure-output criterion of Tang et al. to the thermal regime, where the channel has no pure output. The source further proves additivity of Holevo capacity, minimum output entropy and classical capacity when the channel is used in parallel with any finite-dimensional partner channel (Corollary 1.2); that corollary is not part of this goal.

Formalizing it. The goal ties together an operational capacity (codes and error probabilities), an entropic quantity (a supremum over unrestricted ensembles on 2n2^n2n-dimensional matrices), and a scalar optimization. A machine-checked proof would certify both the matrix entropy inequality at the heart of the argument and the converse direction of the coding theorem. Mathlib does not currently contain quantum channel capacity theory, so the definitions here are new infrastructure.

Difficulty

Optimizing a single use gives a lower bound on capacity; the difficulty is the matching upper bound for all nnn simultaneously, where an input may be entangled across all uses. Holevo additivity is false for general channels, so no abstract argument applies, and for interior thermal occupation 0<ν<10<\nu<10<ν<1 with γ>0\gamma>0γ>0 the channel has no pure output, which rules out the recent pure-output criterion. The proof needs a lower bound on the entropy of the output of an arbitrary entangled block input, i.e. a matrix entropy inequality for 2×22\times22×2 block-structured states that is tight on product inputs. Identifying the capacity with the regularized Holevo quantity also requires the full coding theorem and a Fano-type converse.

Formalization scope

  • Matrices are indexed by Fin n → Fin 2; the nnn-use channel is given by tensor products of the four standard Kraus operators of generalized amplitude damping, which reproduce the displayed action exactly.
  • Entropy is Tr (cfc negMulLog P) (nats, 0ln⁡0=00\ln0=00ln0=0); Holevo values are divided by ln⁡2\ln 2ln2 (bits). holevo γ ν n is the sSup over all finite ensembles of arbitrary 2n2^n2n-dimensional states, with no product or purity restriction.
  • capacity is the sSup of operationally achievable rates: codes with any positive number of messages, arbitrary state encodings, an arbitrary POVM decoder, average error tending to 000, and the liminf rate condition. It is defined operationally, not by regularization, so the identity C=χC=\chiC=χ includes the coding theorem.
  • The scalar objective uses Real.binEntropy (nats); the maximizer is over p∈[0,1]p\in[0,1]p∈[0,1]. The final clause also asserts that the phase-product signals are states.
  • Parameters range over the closed square [0,1]2[0,1]^2[0,1]2, including boundary cases (identity channel at γ=0\gamma=0γ=0, full replacement at γ=1\gamma=1γ=1).
  • Needed infrastructure: matrix entropy and its concavity, CPTP maps and Kraus representations, Holevo's bound, Fano's inequality, the HSW coding theorem for finite ensembles, and Fekete-type superadditivity. These are reusable for any quantum Shannon-theory mission.

Selected references

  • P. Hausladen, R. Jozsa, B. Schumacher, M. Westmoreland, W. K. Wootters, Classical information capacity of a quantum channel, Phys. Rev. A 54 (1996), 1869–1876. https://doi.org/10.1103/PhysRevA.54.1869
  • B. Schumacher and M. D. Westmoreland, Sending classical information via noisy quantum channels, Phys. Rev. A 56 (1997), 131–138. https://doi.org/10.1103/PhysRevA.56.131
  • A. S. Holevo, The capacity of the quantum channel with general signal states, IEEE Trans. Inf. Theory 44 (1998), 269–273. https://doi.org/10.1109/18.651037
  • C. King, Additivity for unital qubit channels, J. Math. Phys. 43 (2002), 4641–4653. https://doi.org/10.1063/1.1500791
  • P. W. Shor, Additivity of the classical capacity of entanglement-breaking quantum channels, J. Math. Phys. 43 (2002), 4334–4340. https://doi.org/10.1063/1.1498000
  • C. H. Bennett, P. W. Shor, J. A. Smolin, A. V. Thapliyal, Entanglement-assisted capacity of a quantum channel and the reverse Shannon theorem, IEEE Trans. Inf. Theory 48 (2002), 2637–2655. https://doi.org/10.1109/TIT.2002.802612
  • V. Giovannetti and R. Fazio, Information-capacity description of spin-chain correlations, Phys. Rev. A 71 (2005), 032314. https://doi.org/10.1103/PhysRevA.71.032314
  • D. W. Berry, Qubit channels that achieve capacity with two states, Phys. Rev. A 71 (2005), 032334. https://doi.org/10.1103/PhysRevA.71.032334
  • L.-Z. Hou and M.-F. Fang, The Holevo capacity of a generalized amplitude-damping channel, Chinese Physics 16 (2007), 1843–1847. https://doi.org/10.1088/1009-1963/16/7/006
  • M. B. Hastings, Superadditivity of communication capacity using entangled inputs, Nature Physics 5 (2009), 255–257. https://doi.org/10.1038/nphys1224
  • F. Leditzky, E. Kaur, N. Datta, M. M. Wilde, Approaches for approximate additivity of the Holevo information of quantum channels, Phys. Rev. A 97 (2018), 012332. https://doi.org/10.1103/PhysRevA.97.012332
  • S. Khatri, K. Sharma, M. M. Wilde, Information-theoretic aspects of the generalized amplitude-damping channel, Phys. Rev. A 102 (2020), 012401. https://doi.org/10.1103/PhysRevA.102.012401
  • Z. Tang, C. Zhu, G. Bai, X. Wang, Classical capacity and entanglement cost of the amplitude damping channel, arXiv:2609.28592 (2026). https://arxiv.org/abs/2609.28592v1
  • K. Fang, A pretty-tight converse on the classical capacity of generalized amplitude-damping channels, arXiv:2609.14621 (2026). https://arxiv.org/abs/2609.14621v1
  • OpenAI, Classical capacity and entropy inequalities for generalized amplitude damping, OpenAI Math Release preprint, September 24, 2026 (Theorem 1.1, p. 3). https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Classical-capacity-and-entropy-inequalities-for-generalized-amplitude-damping-September-24-2026/paper.pdf
37 thms2 active usersReviewed
🏆Completed
Complexity TheoryTheoretical Computer Science·Captain: wurtle

Product-projection localization and the QAC0 parity lower boundResearch Paper

Motivation

Parity, PARITYn(x)=x1⊕⋯⊕xn\mathrm{PARITY}_n(x)=x_1\oplus\cdots\oplus x_nPARITYn​(x)=x1​⊕⋯⊕xn​, is the standard test of whether a shallow circuit can aggregate information from all of its inputs. Classically, constant-depth circuits with unbounded fan-in AND/OR gates (AC0\mathsf{AC}^0AC0) cannot compute it: Furst, Saxe and Sipser (Math. Systems Theory 1984) and Ajtai (APAL 1983) proved superpolynomial lower bounds, and Håstad's switching lemma gave nearly optimal exponential ones (STOC 1986). The quantum analogue QAC0\mathsf{QAC}^0QAC0 allows arbitrary one-qubit unitaries and Toffoli gates of unbounded arity, but gates in one layer must act on disjoint qubits, so a wire cannot control many simultaneous gates. Moore (arXiv:quant-ph/9903046, 1999) showed that coherent parity and quantum fanout are equivalent under Hadamard conjugation, and conjectured that neither is in QAC0\mathsf{QAC}^0QAC0. Because fanout makes constant-depth quantum circuits very powerful (Høyer–Špalek, ToC 2005), whether PARITY∉QAC0\mathrm{PARITY}\notin\mathsf{QAC}^0PARITY∈/QAC0 has been a central question about shallow quantum computation.

Background

  • 1999 — Moore introduces the quantum circuit classes and the parity/fanout conjecture.
  • 2002 — Green, Homer, Moore and Pollett develop quantum counting classes and modular-gate equivalences (QIC 2002).
  • 2006 — Fang, Fenner, Green, Homer and Zhang: depth Ω(log⁡(n/(a+1)))\Omega(\log(n/(a+1)))Ω(log(n/(a+1))) for exact clean parity with aaa ancillas (QIC 2006).
  • 2020–2021 — Padé, Fenner, Grier and Thierauf rule out exact clean parity at entangling depth two (arXiv:2005.12169); Rosenthal bounds approximate parity unitaries (ITCS 2021).
  • 2024–2025 — Pauli-spectrum lower bounds with restricted ancillas (Nadimpalli, Parham, Vasconcelos, Yuen, STOC 2024); no fixed parity advantage with barely superlinear ancillas (Anshu, Dong, Ou, Yao, STOC 2025), extended by Dong, Ou and Yao (arXiv:2510.00593).
  • 2026 — Joshi, Tal, Vasconcelos and Wright: exact parity lower bounds at entangling depth three without size restriction (STOC 2026); Gretta, Gupta and Joshi relate the question to Fourier concentration (arXiv:2604.02793).
  • September 2026 — Two OpenAI preprints claim the lower bound for every fixed depth and every polynomial number of qubits, by independent arguments: product-projection localization (source) and regular trajectories with pruning (source).

This mission asks for a formal proof of the parity lower bound stated as Theorem 1 of an OpenAI preprint dated September 24, 2026 (source). The preprint has not been peer reviewed and its claims have not been independently verified; on the platform the Lean goal is open.

Setting

A circuit acts on NNN qubits, with state space (C2)⊗N(\mathbb C^2)^{\otimes N}(C2)⊗N and computational basis indexed by bit strings y∈{0,1}Ny\in\{0,1\}^Ny∈{0,1}N. The allowed gates are arbitrary one-qubit unitaries and generalized Toffoli gates ∣z1,…,zt,b⟩↦∣z1,…,zt,b⊕z1⋯zt⟩|z_1,\dots,z_t,b\rangle\mapsto|z_1,\dots,z_t,b\oplus z_1\cdots z_t\rangle∣z1​,…,zt​,b⟩↦∣z1​,…,zt​,b⊕z1​⋯zt​⟩ on distinct control and target qubits. A layer is a set of gates with pairwise disjoint supports, and a depth-ddd circuit is a product W=Ld⋯L1W=L_d\cdots L_1W=Ld​⋯L1​ of at most ddd layers. Gates and their complex entries may depend on nnn but not on the input.

On input x∈{0,1}nx\in\{0,1\}^nx∈{0,1}n with n≤Nn\le Nn≤N, the circuit starts in ∣x⟩⊗∣0⟩⊗(N−n)|x\rangle\otimes|0\rangle^{\otimes(N-n)}∣x⟩⊗∣0⟩⊗(N−n). One designated output qubit is measured in the computational basis; all other qubits are discarded in arbitrary states, and the input need not be preserved. The success probability on xxx is

Pr⁡[output=PARITYn(x)]=∑y: yout=PARITYn(x)∣⟨y∣W∣x,0⟩∣2.\Pr[\text{output}=\mathrm{PARITY}_n(x)]=\sum_{y:\ y_{\mathrm{out}}=\mathrm{PARITY}_n(x)}\bigl|\langle y|W|x,0\rangle\bigr|^2 .Pr[output=PARITYn​(x)]=y: yout​=PARITYn​(x)∑​​⟨y∣W∣x,0⟩​2.

A circuit computes parity with worst-case advantage ε\varepsilonε if this probability is at least 1/2+ε1/2+\varepsilon1/2+ε for every xxx.

Formalization targets

Goal: the parity lower bound (Theorem 1, p. 1)

Fix an integer d≥0d\ge0d≥0, a real c≥1c\ge1c≥1 and 0<ε≤1/20<\varepsilon\le1/20<ε≤1/2. Then for all sufficiently large nnn, for every NNN with n≤N≤ncn\le N\le n^cn≤N≤nc, every circuit of depth at most ddd on NNN qubits and every choice of output qubit, there is an input x∈{0,1}nx\in\{0,1\}^nx∈{0,1}n with

Pr⁡[output=PARITYn(x)]<12+ε.\Pr[\text{output}=\mathrm{PARITY}_n(x)]<\tfrac12+\varepsilon .Pr[output=PARITYn​(x)]<21​+ε.

With ε=1/6\varepsilon=1/6ε=1/6 this is the usual bounded-error statement PARITY∉QAC0\mathrm{PARITY}\notin\mathsf{QAC}^0PARITY∈/QAC0 in the measured-output model.

Significance

The result itself. It answers Moore's question for every fixed depth and every polynomial bound on the total number of qubits, with arbitrary final garbage and any fixed positive advantage, whereas earlier results fixed the entangling depth or restricted the number of ancillas. Through Hadamard conjugation it also excludes constant-depth polynomial-size coherent fanout. The preprint derives further consequences (Section 7): state-preparation obstructions (Corollary 8, p. 14), a qualitative quantum LMN-type Fourier tail bound (Corollary 9, p. 16), and, via reductions of Xu and Li, lower bounds for majority and other symmetric functions (Corollaries 10–11, pp. 17–18). The structural Theorem 2 (p. 5), localization of conjugated product projections, is of independent interest.

Formalizing it. The statement is elementary (finite matrices, Born probabilities) while the proof is a two-level induction with polynomial approximation, so a formal proof would give a fully checked circuit lower bound in quantum complexity, an area with few machine-checked results. A companion OpenAI preprint gives an independent proof by regular trajectories and pruning.

Difficulty

Light-cone arguments fail at once: a single Toffoli gate may touch every qubit, so after one layer an output can depend on all inputs. Spectral and Pauli-weight methods previously needed ancilla bounds such as O~(n1+2−d)\tilde O(n^{1+2^{-d}})O~(n1+2−d). The preprint's route controls transitions from a product condition to states with many mismatches, ∥[M≥Nt] U [D=0] U† [M=0]∥≤e−Ns\bigl\|[M\ge N^t]\,U\,[D=0]\,U^\dagger\,[M=0]\bigr\|\le e^{-N^s}​[M≥Nt]U[D=0]U†[M=0]​≤e−Ns for all fixed 0≤s<t0\le s<t0≤s<t (Theorem 2, p. 5). The obstacle is that reflection expansions cost exponentially in the projection's support, and low-degree polynomial surrogates for a product projection are accurate only at low mismatch counts; the uncontrolled high-count error is handled by a second induction on thresholds that exchanges the roles of the two counts.

Formalization scope

  • Word N := Fin N → Fin 2; operators are complex matrices indexed by words. Gate.local q u is a one-qubit unitary tensored with identities; Gate.toffoli controls target is the permutation matrix flipping the target when all controls are 111 (an empty control set gives a NOT gate, itself a one-qubit unitary).
  • PhysicalLayer N is a list of gates whose supports are pairwise disjoint; its matrix is the product of the gate matrices. physicalCircuitMatrix layers applies the layers in list order (rightmost matrix acts first).
  • inputWord x pads xxx with zeros; successProbability W out x sums ∣Wy,(x,0)∣2|W_{y,(x,0)}|^2∣Wy,(x,0)​∣2 over all yyy whose output bit equals the parity of xxx, so final garbage is unrestricted.
  • The output qubit out : Fin N is arbitrary, including an input qubit. The quantifier order is: depth, size bound and advantage fixed first; then a threshold n0n_0n0​; then every n≥n0n\ge n_0n≥n0​, every NNN in range, every circuit and every output qubit.
  • No intermediate measurement, postselection or fanout gate is in the model; the statement is a pure statement about unitary matrices, with no complexity-class machinery.
  • ParityStatement takes the size bound as a real exponent: N≤ncN\le n^cN≤nc with c≥1c\ge1c≥1 real, n≤Nn\le Nn≤N, and depth layers.length ≤ d with d≥0d\ge0d≥0, exactly as in Theorem 1. Its conclusion is the negation of "for every xxx, success ≥1/2+ε\ge 1/2+\varepsilon≥1/2+ε".

Welcome contributions: a reusable Lean model of unitary circuits on qubits, spectral projections of commuting count operators, and Chebyshev-type polynomial approximation with coefficient bounds (Lemma 3, p. 6).

Selected references

  • OpenAI, Product-projection localization and the QAC0 parity lower bound, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Product-projection-localization-and-the-QAC0-parity-lower-bound-September-24-2026/paper.pdf
  • OpenAI, Regular trajectories, pruning and quantum parity, OpenAI Math Release preprint, September 24, 2026. https://github.com/openai/math/blob/adc7f1241b42e322a6451854ab7e4b4c146bf78a/preprints/Regular-trajectories-pruning-and-quantum-parity-September-24-2026/paper.pdf
  • C. Moore, Quantum circuits: fanout, parity, and counting, arXiv:quant-ph/9903046 (1999). https://arxiv.org/abs/quant-ph/9903046v3
  • M. Furst, J. B. Saxe, M. Sipser, Parity, circuits, and the polynomial-time hierarchy, Math. Systems Theory 17 (1984), 13–27. https://doi.org/10.1007/BF01744431
  • M. Ajtai, Σ11\Sigma^1_1Σ11​-formulae on finite structures, Ann. Pure Appl. Logic 24 (1983), 1–48. https://doi.org/10.1016/0168-0072(83)90038-6
  • J. Håstad, Almost optimal lower bounds for small depth circuits, STOC 1986. https://doi.org/10.1145/12130.12132
  • F. Green, S. Homer, C. Moore, C. Pollett, Counting, fanout and the complexity of quantum ACC, QIC 2 (2002). https://doi.org/10.26421/QIC2.1-3
  • P. Høyer, R. Špalek, Quantum fan-out is powerful, Theory of Computing 1 (2005). https://doi.org/10.4086/toc.2005.v001a005
  • M. Fang, S. Fenner, F. Green, S. Homer, Y. Zhang, Quantum lower bounds for fanout, QIC 6 (2006). https://doi.org/10.26421/QIC6.1-3
  • G. Rosenthal, Bounds on the QAC0 complexity of approximating parity, ITCS 2021. https://doi.org/10.4230/LIPIcs.ITCS.2021.32
  • S. Nadimpalli, N. Parham, F. Vasconcelos, H. Yuen, On the Pauli spectrum of QAC0, STOC 2024. https://doi.org/10.1145/3618260.3649662
  • M. R. Joshi, A. Tal, F. Vasconcelos, J. Wright, Improved lower bounds for QAC0, STOC 2026. https://doi.org/10.1145/3798129.3800922
  • L. Gretta, M. Gupta, M. R. Joshi, Parity ∉ QAC0 ⟺ QAC0 is Fourier-concentrated, arXiv:2604.02793 (2026). https://arxiv.org/abs/2604.02793v2
32 thms2 active usersReviewed
🏆Completed
Mathematical Physics·Captain: Lucas

Peres-Terno: the no-communication theorem for commuting Kraus operatorsResearch Paper

Motivation

Two observers, Alice and Bob, share a quantum system. Alice performs some intervention (a measurement, possibly followed by discarding part of her apparatus) and Bob, somewhere else, performs his own. If the statistics of Bob's outcomes could depend on what Alice chose to do, Alice could send Bob a message through the shared system alone, and if the two interventions are spacelike separated this would be a signal faster than light. Quantum mechanics and special relativity coexist peacefully only because this does not happen: the no-communication theorem guarantees that Bob's marginal statistics are independent of Alice's intervention whenever the two interventions are "local" to different parts of the system.

The review article of Peres and Terno, Quantum information and relativity theory (Rev. Mod. Phys. 76, 2004), develops the operational toolkit of quantum information (Kraus matrices, positive-operator-valued measures, completely positive maps) in Sec. II.D and then derives, in Sec. II.E, a sufficient algebraic condition for the absence of instantaneous information transfer: all of Alice's Kraus matrices commute with all of Bob's (their Eq. (9)). This mission formalizes that derivation together with the facts about Kraus matrices and complete positivity that the same pages state.

Setting

A quantum state on a finite-dimensional Hilbert space Cd\mathbb C^dCd is a density matrix ρ\rhoρ: a positive semidefinite d×dd\times dd×d complex matrix with tr⁡ρ=1\operatorname{tr}\rho = 1trρ=1.

An intervention with outcomes μ\muμ is described by Kraus matrices AμmA_{\mu m}Aμm​, where the label mmm ranges over a finite set indexing the subsystems discarded at the end of the interaction. When outcome μ\muμ occurs, the state is updated (Eq. (6)) to the unnormalized matrix

ρμ′=∑mAμm ρ Aμm†,\rho'_\mu = \sum_m A_{\mu m}\,\rho\,A_{\mu m}^\dagger ,ρμ′​=m∑​Aμm​ρAμm†​,

and the probability of outcome μ\muμ is pμ=tr⁡ρμ′p_\mu = \operatorname{tr}\rho'_\mupμ​=trρμ′​. The matrices

Eμ=∑mAμm†AμmE_\mu = \sum_m A_{\mu m}^\dagger A_{\mu m}Eμ​=m∑​Aμm†​Aμm​

(Eq. (8)) are the POVM elements; the intervention is complete when ∑μEμ=1\sum_\mu E_\mu = \mathbb 1∑μ​Eμ​=1.

A map TTT on matrices is positive if it sends positive semidefinite matrices to positive semidefinite matrices, and completely positive if T⊗1T\otimes\mathbb 1T⊗1, acting blockwise on matrices over Cd⊗Cn\mathbb C^d\otimes\mathbb C^nCd⊗Cn, is positive for every ancilla dimension nnn.

For the no-communication theorem, Alice's Kraus matrices AμmA_{\mu m}Aμm​ and Bob's Kraus matrices BνnB_{\nu n}Bνn​ act on the same finite-dimensional space. The probability that Bob obtains ν\nuν, irrespective of Alice's outcome, is (Eq. (10))

pν=∑μtr⁡(∑m,nBνnAμm ρ Aμm†Bνn†).p_\nu = \sum_\mu \operatorname{tr}\Big(\sum_{m,n} B_{\nu n} A_{\mu m}\,\rho\,A_{\mu m}^\dagger B_{\nu n}^\dagger\Big).pν​=μ∑​tr(m,n∑​Bνn​Aμm​ρAμm†​Bνn†​).

Formalization targets

Goal: no-communication (Sec. II.E, Eqs. (9)-(11))

If [Aμm,Bνn]=0[A_{\mu m}, B_{\nu n}] = 0[Aμm​,Bνn​]=0 for all μ,m,ν,n\mu, m, \nu, nμ,m,ν,n, Alice's POVM is complete, and ρ\rhoρ is a density matrix, then for every outcome ν\nuν of Bob

∑μtr⁡(∑m,nBνnAμm ρ Aμm†Bνn†)=tr⁡(∑nBνn ρ Bνn†),\sum_\mu \operatorname{tr}\Big(\sum_{m,n} B_{\nu n} A_{\mu m}\,\rho\,A_{\mu m}^\dagger B_{\nu n}^\dagger\Big) = \operatorname{tr}\Big(\sum_n B_{\nu n}\,\rho\,B_{\nu n}^\dagger\Big),μ∑​tr(m,n∑​Bνn​Aμm​ρAμm†​Bνn†​)=tr(n∑​Bνn​ρBνn†​),

so all of Alice's operators disappear from Bob's statistics.

Milestones (Sec. II.D-II.E)

  1. Eq. (7): ∑mtr⁡(AμmρAμm†)=tr⁡(ρEμ)\sum_m \operatorname{tr}(A_{\mu m}\rho A_{\mu m}^\dagger) = \operatorname{tr}(\rho E_\mu)∑m​tr(Aμm​ρAμm†​)=tr(ρEμ​).
  2. Eq. (8): every EμE_\muEμ​ is positive semidefinite.
  3. Eqs. (7)-(8): for a density matrix and a complete POVM, the pμp_\mupμ​ are nonnegative reals summing to 111.
  4. Eq. (6) defines a completely positive map.
  5. Eq. (6) is the most general completely positive linear map: every completely positive linear map between matrix algebras has a finite Kraus representation.
  6. Time reversal (transposition, i.e. complex conjugation of a Hermitian ρ\rhoρ) is a positive map,
  7. but it is not completely positive.
  8. The exchange step of Sec. II.E: under the commutation hypothesis, tr⁡∑m,nBνnAμmρAμm†Bνn†=tr⁡(Eμ∑nBνnρBνn†)\operatorname{tr}\sum_{m,n} B_{\nu n}A_{\mu m}\rho A_{\mu m}^\dagger B_{\nu n}^\dagger = \operatorname{tr}\big(E_\mu \sum_n B_{\nu n}\rho B_{\nu n}^\dagger\big)tr∑m,n​Bνn​Aμm​ρAμm†​Bνn†​=tr(Eμ​∑n​Bνn​ρBνn†​).

Significance

The no-communication theorem is the consistency check between quantum measurement theory and relativistic causality, and it is the starting point of the relativistic measurement theory developed in Sec. III of the review. The commutation condition (9) is exactly what local quantum field theory provides for spacelike separated regions, so this algebraic form is the one later used for microcausality arguments. The Kraus representation (milestones 4-5) is a basic structural theorem of quantum information theory, and the failure of complete positivity for transposition (milestones 6-7) underlies the partial-transpose entanglement criterion.

The results are classical and their proofs are known; the mission produces machine-checked versions on top of Mathlib's matrix library. The finite Kraus representation theorem (milestone 5) requires Choi-matrix machinery that is not, to the knowledge of this proposal, available in Mathlib in this form; it is reusable well beyond this mission.

Difficulty

The goal and milestones 1-3 and 8 are finite-dimensional trace manipulations; the care needed is in keeping the order of products right, using the adjoint of the commutation relation, and summing over indices in the correct order. Milestone 7 needs an explicit entangled witness on C2⊗Cn\mathbb C^2\otimes\mathbb C^nC2⊗Cn. Milestone 5 is the substantial one: the obvious approach of writing T(ρ)T(\rho)T(ρ) in a basis does not produce Kraus operators directly; a positivity argument (via the Choi matrix of TTT and its spectral decomposition) is needed.

Formalization scope

All Hilbert spaces are finite dimensional: matrices are indexed by arbitrary finite types with complex entries. Outcome sets and discarded-subsystem label sets are finite types; Alice's and Bob's label sets may depend on the outcome. Kraus matrices for the single-intervention milestones may be rectangular (e×de\times de×d), reflecting that the system after the intervention may differ from the original one; in the no-communication theorem both observers' matrices are square on a common space. The order on C\mathbb CC used for "nonnegative probability" is the standard partial order (z≥0z\ge 0z≥0 iff zzz is real and nonnegative). Complete positivity quantifies over ancillas Cn\mathbb C^nCn for every n∈Nn\in\mathbb Nn∈N, with T⊗1T\otimes\mathbb 1T⊗1 acting blockwise. Time reversal is encoded by the linear map ρ↦ρT\rho\mapsto\rho^{T}ρ↦ρT; on Hermitian matrices this coincides with complex conjugation, and only the linear extension gives a meaningful complete-positivity statement.

The goal keeps the hypotheses of the source setting (both POVMs complete, ρ\rhoρ a density matrix); none of them can be dropped silently by a degenerate encoding, and the conclusion is an identity of complex numbers, so no hypothesis is vacuous: for instance, single-outcome trivial interventions with A=B=1A = B = \mathbb 1A=B=1 satisfy all of them.

Note: the sentence on p. 100 of the source claiming that commuting POVM elements are necessarily orthogonal projections is not included; as stated it fails (e.g. E1=E2=121E_1 = E_2 = \tfrac12\mathbb 1E1​=E2​=21​1).

Selected references

  • A. Peres and D. R. Terno, Quantum information and relativity theory, Rev. Mod. Phys. 76, 93-123 (2004). https://doi.org/10.1103/RevModPhys.76.93
  • K. Kraus, States, Effects, and Operations, Lecture Notes in Physics 190, Springer (1983). https://doi.org/10.1007/3-540-12732-1
  • M.-D. Choi, Completely positive linear maps on complex matrices, Linear Algebra Appl. 10, 285-290 (1975). https://doi.org/10.1016/0024-3795(75)90075-0
10 thms2 active usersReviewed
🏆Completed
Number TheoryProbabilityTheoretical Computer Science·Captain: mikedeng1

Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer 4: The Discrete Logarithm Circuit Gives a Good Output with Probability at Least 1/480Research Paper

Motivation

The discrete logarithm problem modulo a prime asks, given a prime ppp, a generator ggg of the multiplicative group modulo ppp, and a nonzero residue xxx, for the exponent rrr with gr≡x(modp)g^r\equiv x \pmod pgr≡x(modp). Its presumed classical hardness underlies Diffie–Hellman key exchange, ElGamal encryption and the Digital Signature Algorithm. The best classical algorithm known when Shor wrote, Gordon's adaptation of the number field sieve, runs in time exp⁡(O((log⁡p)1/3(log⁡log⁡p)2/3))\exp(O((\log p)^{1/3}(\log\log p)^{2/3}))exp(O((logp)1/3(loglogp)2/3)).

In §6 of Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer (SIAM J. Comput. 26(5), 1997; doi:10.1137/S0097539795293172, arXiv:quant-ph/9508027), Shor gave a quantum algorithm that uses two modular exponentiations and two quantum Fourier transforms and outputs, with constant probability, a pair from which rrr can be computed. The quantitative core of that analysis is a single number: the circuit produces a "good" output with probability at least 1/4801/4801/480. This mission formalizes that bound and the three estimates it is assembled from.

Setting

Let ppp be a prime and ggg a generator of (Z/pZ)×(\mathbb Z/p\mathbb Z)^\times(Z/pZ)×, so that 1,g,…,gp−21,g,\dots,g^{p-2}1,g,…,gp−2 are all the nonzero residues. Fix the unknown rrr with 0≤r<p−10\le r<p-10≤r<p−1 and put x=grx=g^rx=gr. Let q=2lq=2^lq=2l be the power of 222 with p<q<2pp<q<2pp<q<2p.

The Fourier matrix AqA_qAq​ is the q×qq\times qq×q matrix with entries (Aq)a,c=q−1/2exp⁡(2πi ac/q)(A_q)_{a,c}=q^{-1/2}\exp(2\pi i\,ac/q)(Aq​)a,c​=q−1/2exp(2πiac/q) for 0≤a,c<q0\le a,c<q0≤a,c<q (§4, eq. (4.1)). Rows index input basis vectors and columns output basis vectors.

The algorithm uses three registers: two holding numbers 0≤a,b<q0\le a,b<q0≤a,b<q and one holding a nonzero residue modulo ppp. It starts from the state

1p−1∑a=0p−2∑b=0p−2∣a,b,gax−b (mod p)⟩(6.1)\frac{1}{p-1}\sum_{a=0}^{p-2}\sum_{b=0}^{p-2}|a,b,g^ax^{-b}\ (\mathrm{mod}\ p)\rangle \qquad (6.1)p−11​a=0∑p−2​b=0∑p−2​∣a,b,gax−b (mod p)⟩(6.1)

(preFourierState), applies AqA_qAq​ to each of the first two registers (finalState), and measures all three registers. The probability of observing ∣c,d,y⟩|c,d,y\rangle∣c,d,y⟩ is the squared modulus of its amplitude (outcomeProb).

For integers zzz and q>0q>0q>0, the symmetric residue {z}q\{z\}_q{z}q​ is the residue of zzz modulo qqq in (−q/2,q/2](-q/2,q/2](−q/2,q/2] (symmRes). Put

T=rc+d−rp−1{c(p−1)}q.T=rc+d-\frac{r}{p-1}\{c(p-1)\}_q .T=rc+d−p−1r​{c(p−1)}q​.

An observed state ∣c,d,y⟩|c,d,y\rangle∣c,d,y⟩ is good (IsGood) when

∣{T}q∣≤12(6.10)and∣{c(p−1)}q∣≤q/12(6.11).|\{T\}_q|\le\tfrac12 \quad (6.10) \qquad\text{and}\qquad |\{c(p-1)\}_q|\le q/12 \quad (6.11).∣{T}q​∣≤21​(6.10)and∣{c(p−1)}q​∣≤q/12(6.11).

Goodness depends only on (c,d)(c,d)(c,d).

Formalization targets

Goal: a good output with probability at least 1/4801/4801/480 (§6, p. 1504)

∑0≤c,d<q(c,d) good ∑y∈(Z/p)×Pr⁡[c,d,y] ≥ 1480.\sum_{\substack{0\le c,d<q\\ (c,d)\ \text{good}}}\ \sum_{y\in(\mathbb Z/p)^\times}\Pr[c,d,y]\ \ge\ \frac1{480}.0≤c,d<q(c,d) good​∑​ y∈(Z/p)×∑​Pr[c,d,y] ≥ 4801​.

The constant is the one the page carries forward. The goal fixes no threshold on ppp: it is stated for every prime ppp that admits a power of two strictly between ppp and 2p2p2p.

Milestones

  1. The output distribution, eq. (6.4). For 0≤k<p−10\le k<p-10≤k<p−1,
Pr⁡[c,d,gk]=∣1(p−1)q∑0≤a,b≤p−2a−rb≡k (p−1)exp⁡(2πiq(ac+bd))∣2.\Pr[c,d,g^k]=\left|\frac{1}{(p-1)q}\sum_{\substack{0\le a,b\le p-2\\ a-rb\equiv k\ (p-1)}}\exp\Bigl(\frac{2\pi i}{q}(ac+bd)\Bigr)\right|^2 .Pr[c,d,gk]=​(p−1)q1​0≤a,b≤p−2a−rb≡k (p−1)​∑​exp(q2πi​(ac+bd))​2.
  1. Each good state is likely, eq. (6.17). If (c,d)(c,d)(c,d) is good, then Pr⁡[c,d,y]≥1/(20q2)\Pr[c,d,y]\ge 1/(20q^2)Pr[c,d,y]≥1/(20q2) for every yyy.
  2. Many good pairs (p. 1504). At least q/12q/12q/12 pairs (c,d)(c,d)(c,d) are good.
  3. Each good ccc is likely (p. 1504). If (c,d)(c,d)(c,d) is good for some ddd, then ∑d′,yPr⁡[c,d′,y]≥(p−1)/(20q2)≥1/(40q)\sum_{d',y}\Pr[c,d',y]\ge(p-1)/(20q^2)\ge1/(40q)∑d′,y​Pr[c,d′,y]≥(p−1)/(20q2)≥1/(40q).

Significance

The result. The bound 1/4801/4801/480 is what turns the circuit into an algorithm. Repeating the circuit O(1)O(1)O(1) times in expectation yields a good output, and from a good pair (c,d)(c,d)(c,d) one reads off an equation that determines rrr modulo divisors of p−1p-1p−1 (§6, eqs. (6.18)–(6.20)). Together with the quantum Fourier transform circuit and reversible modular exponentiation, this places the discrete logarithm modulo a prime in quantum polynomial time. Every later analysis of quantum attacks on discrete-logarithm cryptography starts from this success probability or a sharpened version of it.

Formalizing it. The result has been proved since 1994–1997 and is textbook material; it is not open. As far as is known, no machine-checked proof of Shor's discrete-logarithm analysis exists. The paper's proof of eq. (6.17) replaces a sum by an integral with an error term O(W/(pq))O(W/(pq))O(W/(pq)) whose constant is not given, yet states 1/(20q2)1/(20q^2)1/(20q2) for every prime. A formal proof must therefore either control that error explicitly or find another argument, and so settles a point the paper leaves informal. Numerically, the smallest value of q2Pr⁡[c,d,y]q^2\Pr[c,d,y]q2Pr[c,d,y] over good states is about 0.490.490.49 for all primes p<90p<90p<90, so the unconditional claim is not in doubt for small ppp. The page also contains two small slips, recorded under Formalization scope; a complete development pins down exactly what is true.

Difficulty

The exponential sum (6.4) runs over pairs (a,b)(a,b)(a,b) satisfying a congruence modulo p−1p-1p−1, while the phases are taken modulo qqq. The two moduli are unrelated: qqq is a power of two and p−1p-1p−1 is arbitrary. Eliminating aaa through the congruence introduces a floor function ⌊(br+k)/(p−1)⌋\lfloor(br+k)/(p-1)\rfloor⌊(br+k)/(p−1)⌋, and the resulting phase is not linear in bbb. The obvious estimate treats the sum as a geometric series in bbb and bounds it by its first-order phase; this fails because the floor term perturbs every phase by an amount of size up to ∣{c(p−1)}q∣|\{c(p-1)\}_q|∣{c(p−1)}q​∣. Condition (6.11) only keeps this perturbation within π/6\pi/6π/6 of the main phase; it does not remove it. The per-state bound must survive this perturbation uniformly in ppp, rrr and kkk, including small primes where the paper's integral approximation gives no explicit control.

The count of good pairs needs a separate argument about how often a multiple c(p−1)c(p-1)c(p−1) lies within q/12q/12q/12 of a multiple of qqq when gcd⁡(p−1,q)\gcd(p-1,q)gcd(p−1,q) is large.

Formalization scope

  • States are functions Fin q × Fin q × (ZMod p)ˣ → ℂ. The first two registers range over {0,…,q−1}\{0,\dots,q-1\}{0,…,q−1}; the third over the units modulo ppp.
  • Matrix convention. Following §2, rows are inputs, so the amplitude of ∣c,d,y⟩|c,d,y\rangle∣c,d,y⟩ after the transforms is ∑a,bψ(a,b,y)(Aq)a,c(Aq)b,d\sum_{a,b}\psi(a,b,y)(A_q)_{a,c}(A_q)_{b,d}∑a,b​ψ(a,b,y)(Aq​)a,c​(Aq​)b,d​. finalState is defined this way from (6.1) and AqA_qAq​. It is not typed in as the closed form (6.3) or (6.4). A formalization that defined the final state by (6.4) directly would make milestone 1 trivial, and is ruled out.
  • Probability of a basis state is the squared norm of its amplitude, with no normalization hypothesis.
  • Parameters. ppp is prime (Fact p.Prime). The generator is encoded as orderOf g = p - 1. r<p−1r<p-1r<p−1 is a parameter, with x=grx=g^rx=gr. qqq is given by q = 2 ^ l together with p<q<2pp<q<2pp<q<2p. No large-ppp threshold is added anywhere.
  • Arithmetic. x−bx^{-b}x−b is x⁻¹ ^ b in the unit group. p−1p-1p−1 is computed in Z\mathbb ZZ and R\mathbb RR inside TTT and the congruences, and as natural-number subtraction only where p≥2p\ge2p≥2 makes it exact. TTT is real.
  • Condition (6.10) is stated as "some integer jjj has ∣T−jq∣≤12|T-jq|\le\frac12∣T−jq∣≤21​". Because q≥4q\ge4q≥4, this is equivalent to the page's form with jjj the closest integer to T/qT/qT/q.
  • Not formalized. The preparation of (6.1) by testing and restarting is not formalized; the state (6.1) is taken as displayed. The printed test "whether the number is less than ppp" should read p−1p-1p−1, as the sums in (6.1) show. Also out of scope: the recovery of rrr (eqs. (6.18)–(6.20)), the repetition count "480t480t480t", and all running-time claims.
  • Printed slips.
    • The page asserts that for each ccc there is exactly one ddd satisfying (6.10). At a tie {T}q=±12\{T\}_q=\pm\frac12{T}q​=±21​ there can be two such ddd. Milestone 3 states only the count, which needs at least one.
    • The page's intermediate bound "at least p/(240q)p/(240q)p/(240q)" should be (p−1)/(240q)(p-1)/(240q)(p−1)/(240q). The conclusion 1/4801/4801/480 is unaffected, since qqq and 2p2p2p are both even and so q≤2(p−1)q\le 2(p-1)q≤2(p−1). Only 1/4801/4801/480 is stated.

Needed infrastructure: finite exponential sums and their modulus, the symmetric residue and its basic properties, and counting multiples in residue classes of Z/q\mathbb Z/qZ/q. The exponential-sum estimates of milestones 1 and 2 are reusable in the order-finding analysis of §5 of the same paper. Proofs of any milestone, of the normalization ∑Pr⁡=1\sum\Pr=1∑Pr=1, and of auxiliary lemmas about symmRes are welcome.

Selected references

  • P. W. Shor, Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer, SIAM J. Comput. 26(5):1484–1509, 1997. https://doi.org/10.1137/S0097539795293172 (preprint: https://arxiv.org/abs/quant-ph/9508027)
  • D. M. Gordon, Discrete logarithms in GF(p) using the number field sieve, SIAM J. Discrete Math. 6(1):124–138, 1993. https://doi.org/10.1137/0406010
  • W. Diffie and M. E. Hellman, New directions in cryptography, IEEE Trans. Inform. Theory 22(6):644–654, 1976. https://doi.org/10.1109/TIT.1976.1055638
  • M. A. Nielsen and I. L. Chuang, Quantum Computation and Quantum Information, Cambridge University Press, 2000. https://doi.org/10.1017/CBO9780511976667
11 thms2 active usersReviewed
🏆Completed
Harmonic AnalysisTheoretical Computer Science·Captain: mikedeng1

Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer 1: The Quantum Fourier Transform Circuit Computes A_q up to Bit ReversalResearch Paper

Motivation

Shor's factoring and discrete logarithm algorithms (Shor 1997) reduce both problems to sampling from the output of a quantum Fourier transform: a register holding a superposition with a hidden period is transformed, then measured, and the measured value carries information about the period. The whole speed-up rests on one engineering fact: for q=2lq = 2^lq=2l the q×qq \times qq×q Fourier matrix, which has q2=4lq^2 = 4^lq2=4l entries, can be applied by a quantum circuit of only O(l2)O(l^2)O(l2) elementary gates, each acting on one or two bits.

That circuit was found independently by Coppersmith (IBM RC 19642, 1994) and Deutsch, and Shor's §4 presents it following Ekert and Jozsa (Rev. Mod. Phys. 68, 1996). It is the quantum analogue of the radix-2 fast Fourier transform, and it reappears in phase estimation, in the hidden subgroup algorithms for abelian groups, and in every textbook account of quantum computation. This mission formalizes Shor's statement that the circuit computes the Fourier matrix, up to a reversal of the output bits, together with the two displayed steps of its verification.

Setting

A register of lll bits has one basis state ∣a⟩=∣al−1al−2…a0⟩|a\rangle = |a_{l-1} a_{l-2} \dots a_0\rangle∣a⟩=∣al−1​al−2​…a0​⟩ for every bit string, with a0a_0a0​ the least significant bit; the string encodes the integer

a=∑j=0l−12jaj,0≤a<q=2l.a = \sum_{j=0}^{l-1} 2^j a_j, \qquad 0 \le a < q = 2^l .a=j=0∑l−1​2jaj​,0≤a<q=2l.

A state is a vector of complex amplitudes, one per basis state. A gate is a matrix whose rows are indexed by input basis vectors and whose columns are indexed by output basis vectors (§2, p. 1489). A gate on some of the bits acts on those bits through its matrix and leaves the others alone (§2, p. 1490): the output amplitude at a string bbb is the sum, over the possible input values uuu of the acted-on bits, of the input amplitude at bbb with those bits replaced by uuu, times the matrix entry from uuu to the corresponding bits of bbb.

The Fourier matrix AqA_qAq​ (eq. (4.1)) is the q×qq \times qq×q matrix with (a,c)(a, c)(a,c) entry q−1/2exp⁡(2πi ac/q)q^{-1/2}\exp(2\pi i\,ac/q)q−1/2exp(2πiac/q); it takes ∣a⟩|a\rangle∣a⟩ to q−1/2∑c=0q−1exp⁡(2πi ac/q) ∣c⟩q^{-1/2}\sum_{c=0}^{q-1}\exp(2\pi i\,ac/q)\,|c\rangleq−1/2∑c=0q−1​exp(2πiac/q)∣c⟩.

The circuit uses two gates (eqs. (4.2), (4.3)):

  • RjR_jRj​ acts on bit jjj with matrix 12(111−1)\frac{1}{\sqrt 2}\begin{pmatrix} 1 & 1 \\ 1 & -1\end{pmatrix}2​1​(11​1−1​);
  • Sj,kS_{j,k}Sj,k​, for j<kj < kj<k, acts on bits jjj and kkk with matrix diag(1,1,1,eiθk−j)\mathrm{diag}(1, 1, 1, e^{i\theta_{k-j}})diag(1,1,1,eiθk−j​), where θk−j=π/2k−j\theta_{k-j} = \pi / 2^{k-j}θk−j​=π/2k−j; it multiplies the amplitude of a basis state by eiθk−je^{i\theta_{k-j}}eiθk−j​ when bits jjj and kkk are both 111.

The circuit is the gate sequence (4.4), applied from left to right:

Rl−1 Sl−2,l−1 Rl−2 Sl−3,l−1 Sl−3,l−2 Rl−3⋯R1 S0,l−1 S0,l−2⋯S0,2 S0,1 R0,R_{l-1}\, S_{l-2,l-1}\, R_{l-2}\, S_{l-3,l-1}\, S_{l-3,l-2}\, R_{l-3} \cdots R_1\, S_{0,l-1}\, S_{0,l-2} \cdots S_{0,2}\, S_{0,1}\, R_0 ,Rl−1​Sl−2,l−1​Rl−2​Sl−3,l−1​Sl−3,l−2​Rl−3​⋯R1​S0,l−1​S0,l−2​⋯S0,2​S0,1​R0​,

that is, for j=l−1,…,0j = l-1, \dots, 0j=l−1,…,0 it applies Sj,l−1,…,Sj,j+1S_{j,l-1}, \dots, S_{j,j+1}Sj,l−1​,…,Sj,j+1​ and then RjR_jRj​. On three bits it is R2S1,2R1S0,2S0,1R0R_2 S_{1,2} R_1 S_{0,2} S_{0,1} R_0R2​S1,2​R1​S0,2​S0,1​R0​.

The bit reversal of a string bbb is the string ccc with ck=bl−1−kc_k = b_{l-1-k}ck​=bl−1−k​.

Formalization targets

Goal: the circuit computes AqA_qAq​ up to bit reversal (§4, p. 1496)

For every l≥0l \ge 0l≥0 and every basis state ∣a⟩|a\rangle∣a⟩, with q=2lq = 2^lq=2l,

circuit ∣a⟩=1q1/2∑bexp⁡(2πi ac/q) ∣b⟩,c=bit reversal of b.\text{circuit}\,|a\rangle = \frac{1}{q^{1/2}}\sum_{b}\exp(2\pi i\,ac/q)\,|b\rangle, \qquad c = \text{bit reversal of } b .circuit∣a⟩=q1/21​b∑​exp(2πiac/q)∣b⟩,c=bit reversal of b.

The sum runs over all lll-bit strings bbb. Reading the output register in reverse order therefore yields Aq∣a⟩A_q|a\rangleAq​∣a⟩.

Milestone 1: the amplitude along the circuit (§4, eq. (4.5), p. 1496)

The amplitude of ∣b⟩|b\rangle∣b⟩ in circuit ∣a⟩\text{circuit}\,|a\ranglecircuit∣a⟩ is

2−l/2exp⁡(i(∑0≤j<lπajbj+∑0≤j<k<lπ2k−jajbk)).2^{-l/2}\exp\Big(i\Big(\sum_{0\le j<l}\pi a_jb_j + \sum_{0\le j<k<l}\frac{\pi}{2^{k-j}}a_jb_k\Big)\Big).2−l/2exp(i(0≤j<l∑​πaj​bj​+0≤j<k<l∑​2k−jπ​aj​bk​)).

Milestone 2: the phase identity (§4, eqs. (4.6)–(4.10), pp. 1496–1497)

For bit strings a,ba, ba,b, with ccc the bit reversal of bbb and a,ca, ca,c their values,

exp⁡(i(∑0≤j<lπajbj+∑0≤j<k<lπ2k−jajbk))=exp⁡(2πi ac/q).\exp\Big(i\Big(\sum_{0\le j<l}\pi a_jb_j + \sum_{0\le j<k<l}\frac{\pi}{2^{k-j}}a_jb_k\Big)\Big) = \exp(2\pi i\,ac/q).exp(i(0≤j<l∑​πaj​bj​+0≤j<k<l∑​2k−jπ​aj​bk​))=exp(2πiac/q).

Significance

The result. The goal says that AqA_qAq​, a dense unitary on 2l2^l2l amplitudes, is realized by lll one-bit gates and l(l−1)/2l(l-1)/2l(l−1)/2 two-bit gates, followed by a relabelling of the output. This is what makes the Fourier sampling step of the factoring algorithm (§5) and of the discrete logarithm algorithm (§6) polynomial in the number of bits. Without it, the analyses of those sections describe measurements of states that no efficient circuit is known to prepare. The bit reversal is a real part of the statement: the circuit does not compute AqA_qAq​ itself for l≥2l \ge 2l≥2, and an implementation must either permute the output bits or read them in reverse order.

Formalizing it. The identity is classical and its proof is short on paper; its content lies in bookkeeping that is easy to get wrong: which bit a gate touches, which of the input or output value a bit holds when a phase gate acts, and the order in which the gates are applied. A machine-checked version fixes all of these conventions explicitly and yields a reusable model of gate-level circuits on bit strings. To the drafter's knowledge there is no Lean formalization of this circuit on Prove2Me; the platform's FastFourierTransform rows state the classical Cooley–Tukey recursion for the unnormalized transform with the opposite sign, which is a different object.

Difficulty

Multiplying out the gate matrices is not an argument beyond tiny lll: the difficulty is to control the whole product of l(l+1)/2l(l+1)/2l(l+1)/2 gates symbolically. The paper's argument (§4, p. 1496) is informal at exactly the points a formal proof must make precise: that only one sequence of intermediate basis states from ∣a⟩|a\rangle∣a⟩ to ∣b⟩|b\rangle∣b⟩ carries nonzero amplitude, and that at the moment Sj,kS_{j,k}Sj,k​ acts, bit kkk already holds its output value bkb_kbk​ while bit jjj still holds its input value aja_jaj​. Both facts depend on the order (4.4); applying the same gates in the reverse order gives a different unitary, whose phase pairs bjb_jbj​ with aka_kak​. The phase identity then holds only modulo 2π2\pi2π, not as an equality of real numbers, so it cannot be closed by rearranging sums alone.

Formalization scope

  • States. A bit string on lll bits is Fin l → Fin 2, with bit j equal to aja_jaj​ and value ∑j2jaj\sum_j 2^j a_j∑j​2jaj​ (least significant bit first). A state is a function from bit strings to ℂ. The case l=0l = 0l=0 (q=1q = 1q=1, empty circuit) is included.
  • Gate application follows the row = input convention: a one-bit gate MMM on bit jjj sends ψ\psiψ to b↦∑uψ(b[j↦u]) Mu,bjb \mapsto \sum_u \psi(b[j\mapsto u])\,M_{u, b_j}b↦∑u​ψ(b[j↦u])Mu,bj​​, and a two-bit gate acts analogously on the pair (bit jjj, bit kkk).
  • The circuit is the gate list (4.4) — an explicit list of constructors R j and S j k — run by a left fold, so the leftmost gate acts first. It is not defined as the matrix AqA_qAq​ or by its entries; a definition of that kind would make the goal true by unfolding and is ruled out.
  • Angles. θk−j=π/2k−j\theta_{k-j} = \pi/2^{k-j}θk−j​=π/2k−j uses natural-number subtraction, which is exact because the list only contains Sj,kS_{j,k}Sj,k​ with j<kj < kj<k.
  • Normalization. The prefactor q−1/2q^{-1/2}q−1/2 is written (2l)−1(\sqrt{2^l})^{-1}(2l​)−1.
  • Basis form. The goal is stated on basis states, as on the page; by linearity it determines the circuit on every state.
  • Not stated. The gate count: the paper's sentence "we thus need to use l(l−1)/2l(l-1)/2l(l−1)/2 quantum gates" (p. 1496) counts only the gates Sj,kS_{j,k}Sj,k​; the sequence (4.4) also contains the lll gates RjR_jRj​, l(l+1)/2l(l+1)/2l(l+1)/2 gates in all. The approximate transform of Coppersmith and the one-bit construction of Griffiths and Niu (p. 1497) are out of scope, as is the polynomial-time claim.

Welcome contributions: proofs of the two milestones and the goal; general lemmas about one- and two-bit gate actions on Fin l → Fin 2 states (commutation of gates on disjoint bits, linearity, action on basis states), which are reusable for any gate-level circuit; and a proof that the circuit is unitary.

Selected references

  • P. W. Shor, Polynomial-Time Algorithms for Prime Factorization and Discrete Logarithms on a Quantum Computer, SIAM J. Comput. 26(5):1484–1509, 1997. https://doi.org/10.1137/S0097539795293172
  • D. Coppersmith, An approximate Fourier transform useful in quantum factoring, IBM Research Report RC 19642, 1994. https://arxiv.org/abs/quant-ph/0201067
  • A. Ekert and R. Jozsa, Quantum computation and Shor's factoring algorithm, Rev. Mod. Phys. 68:733–753, 1996. https://doi.org/10.1103/RevModPhys.68.733
  • R. B. Griffiths and C.-S. Niu, Semiclassical Fourier transform for quantum computation, Phys. Rev. Lett. 76:3228–3231, 1996. https://doi.org/10.1103/PhysRevLett.76.3228
5 thms2 active usersReviewed
🏆Completed
Mathematical Physics·Captain: Lucas

Decoherence in Two-Flavour Neutrino Oscillations (Alves 2020)Research Paper

Motivation

Neutrino oscillations are an interference phenomenon, so they are sensitive to any mechanism that degrades quantum coherence during propagation. The open quantum system description models such effects without specifying their microscopic origin: the neutrino is treated as a subsystem whose reduced density matrix evolves under a Markovian master equation. The dissertation of G. F. S. Alves (Instituto de Física, Universidade de São Paulo, 2020) develops this framework for two and three neutrino generations and uses public IceCube atmospheric data to constrain the new parameters. Its two-generation chapter derives closed-form oscillation probabilities under several choices of the dissipative term; this mission formalizes that chapter's mathematical content.

Setting

A density matrix is a positive semidefinite complex matrix of unit trace. For a two-level system (two neutrino mass eigenstates) the reduced state ρ(t)∈C2×2\rho(t)\in\mathbb C^{2\times2}ρ(t)∈C2×2 obeys the Lindblad equation with constant coefficients,

ρ˙=−i[H,ρ]+D[ρ],D[ρ]=12∑i,j=13aij(2σiρσj−{σjσi,ρ}),\dot\rho=-i[H,\rho]+D[\rho],\qquad D[\rho]=\tfrac12\sum_{i,j=1}^3 a_{ij}\bigl(2\sigma_i\rho\sigma_j-\{\sigma_j\sigma_i,\rho\}\bigr),ρ˙​=−i[H,ρ]+D[ρ],D[ρ]=21​i,j=1∑3​aij​(2σi​ρσj​−{σj​σi​,ρ}),

where σ1,σ2,σ3\sigma_1,\sigma_2,\sigma_3σ1​,σ2​,σ3​ are the Pauli matrices and the 3×33\times33×3 Kossakowski matrix a=(aij)a=(a_{ij})a=(aij​) is Hermitian positive semidefinite (the condition for complete positivity). In the mass basis the vacuum Hamiltonian is H=diag⁡(0,Δm2/2E)H=\operatorname{diag}(0,\Delta m^2/2E)H=diag(0,Δm2/2E), and the mixing matrix with Majorana phase α\alphaα is

U=(cos⁡θeiαsin⁡θ−e−iαsin⁡θcos⁡θ).U=\begin{pmatrix}\cos\theta&e^{i\alpha}\sin\theta\\-e^{-i\alpha}\sin\theta&\cos\theta\end{pmatrix}.U=(cosθ−e−iαsinθ​eiαsinθcosθ​).

The muon neutrino is produced in the state ρμ(0)=U†∣ν1⟩⟨ν1∣U\rho_\mu(0)=U^\dagger|\nu_1\rangle\langle\nu_1|Uρμ​(0)=U†∣ν1​⟩⟨ν1​∣U, and the survival probability after a baseline LLL is Pνμ→νμ=Tr⁡{ρ(L)ρμ(0)}P_{\nu_\mu\to\nu_\mu}=\operatorname{Tr}\{\rho(L)\rho_\mu(0)\}Pνμ​→νμ​​=Tr{ρ(L)ρμ​(0)}. The linear entropy Sl(ρ)=1−Tr⁡ρ2S_l(\rho)=1-\operatorname{Tr}\rho^2Sl​(ρ)=1−Trρ2 and the von Neumann entropy S(ρ)=−Tr⁡ρln⁡ρS(\rho)=-\operatorname{Tr}\rho\ln\rhoS(ρ)=−Trρlnρ measure the loss of purity.

Formalization targets

Goal: decoherent survival probability, eq. (3.24)

For the pure-decoherence Kossakowski matrix a=diag⁡(0,0,γ)a=\operatorname{diag}(0,0,\gamma)a=diag(0,0,γ), γ≥0\gamma\ge0γ≥0, every solution with ρ(0)=ρμ(0)\rho(0)=\rho_\mu(0)ρ(0)=ρμ​(0) satisfies, for all L≥0L\ge0L≥0,

Pνμ→νμ=14(3+cos⁡4θ+2e−2γLcos⁡(Δm22EL)sin⁡22θ).P_{\nu_\mu\to\nu_\mu}=\frac14\Bigl(3+\cos4\theta+2e^{-2\gamma L}\cos\Bigl(\frac{\Delta m^2}{2E}L\Bigr)\sin^22\theta\Bigr).Pνμ​→νμ​​=41​(3+cos4θ+2e−2γLcos(2EΔm2​L)sin22θ).

Milestones

  1. Range of the linear entropy, eq. (2.45), and of the von Neumann entropy, eqs. (2.40)–(2.41).
  2. Non-negativity of the relative entropy, eq. (2.43).
  3. Entropy increase whenever the maximally mixed state is a fixed point of DDD (Appendix B.1, two-level case).
  4. Pauli-basis form of the Dissipator, eq. (3.5).
  5. Average-energy conservation holds iff a=diag⁡(0,0,γ)a=\operatorname{diag}(0,0,\gamma)a=diag(0,0,γ), eqs. (3.19)–(3.22).
  6. The evolved state in case (1), eq. (3.23), and the population in case (4), eq. (3.26).

Significance

The goal is the formula used in the dissertation to show how decoherence damps oscillations: it reduces to the standard vacuum formula at γ=0\gamma=0γ=0, it tends to 14(3+cos⁡4θ)\tfrac14(3+\cos4\theta)41​(3+cos4θ) at large LLL instead of oscillating, and it does not depend on α\alphaα. The milestones connect the formula to the structure of the Lindblad generator (energy conservation singles out this Dissipator) and to entropy growth. A formal development yields a checked two-level Lindblad toolkit (Pauli expansion of dissipators, explicit solutions, entropy bounds) reusable in other open-system missions. The results are proved informally in the source; none is known to be machine-checked.

Difficulty

The ODE part requires a uniqueness argument for a linear matrix ODE on [0,∞)[0,\infty)[0,∞) given only one-sided derivatives at t=0t=0t=0, plus exact trigonometric bookkeeping. The entropy milestones need spectral calculus for Hermitian matrices (Klein's inequality, concavity of −xln⁡x-x\ln x−xlnx) and, for Appendix B.1, the monotonicity of relative entropy under completely positive trace-preserving maps together with complete positivity of Lindblad semigroups, neither of which is expected to be available off the shelf.

Formalization scope

Matrices are Matrix (Fin 2) (Fin 2) ℂ; the Kossakowski matrix is Matrix (Fin 3) (Fin 3) ℂ with indices 0,1,20,1,20,1,2 for the thesis' 1,2,31,2,31,2,3. Solutions are curves R→C2×2\mathbb R\to\mathbb C^{2\times2}R→C2×2 satisfying the Lindblad equation entrywise for t≥0t\ge0t≥0 (derivatives within [0,∞)[0,\infty)[0,∞)); the exponential eLte^{\mathcal Lt}eLt is not used. Entropies use the continuous functional calculus. Two printed formulas are corrected in the drafts: the sign of the ρ0 Im(aij)\rho_0\,\mathrm{Im}(a_{ij})ρ0​Im(aij​) terms in eq. (3.5), and the overall sign of the off-diagonal entries in eq. (3.23); both corrections are explained in the corresponding items. The two-generation cases (4) and (5) survival probabilities, eqs. (3.25) and (3.27), and the three-generation and IceCube analyses are out of scope.

Selected references

  • G. F. S. Alves, Decoherence in Neutrino Oscillations in the IceCube Experiment, MSc dissertation, Instituto de Física, Universidade de São Paulo, 2020.
  • G. Lindblad, On the generators of quantum dynamical semigroups, Commun. Math. Phys. 48 (1976) 119–130, https://doi.org/10.1007/BF01608499
  • V. Gorini, A. Kossakowski, E. C. G. Sudarshan, Completely positive dynamical semigroups of N-level systems, J. Math. Phys. 17 (1976) 821–825, https://doi.org/10.1063/1.522979
10 thms2 active usersReviewed
🏆Completed
Quantum Error Correction·Captain: Lucas

Knill-Laflamme-Zurek: Resilient Quantum Computation and the Error ThresholdResearch Paper

Motivation

A quantum computer manipulates superpositions of many degrees of freedom, and every gate it applies is imperfect. Before 1996 it was widely argued that this fragility is fatal: the accuracy demanded of each elementary operation appeared to grow with the length of the computation, so that arbitrarily long computations would require arbitrarily perfect hardware. The discovery of quantum error-correcting codes, of fault tolerant syndrome extraction, and of transversally encoded operations changed that picture, but the first fault tolerant schemes still required the error per operation to shrink as the computation grew.

Knill, Laflamme and Zurek, Resilient Quantum Computation: Error Models and Thresholds (arXiv:quant-ph/9702058), together with the independent work of Kitaev and of Aharonov and Ben-Or, removed that requirement: they exhibited a fixed threshold value such that any physical error rate below it permits arbitrarily accurate encoded computation, at an overhead that is only polylogarithmic in the length of the computation. Their paper is also the one that pushes past independent stochastic noise: its thresholds are proved for quasi-independent error models, which accommodate coherent over-rotation and weak residual couplings between neighbouring qubits.

This mission formalizes the quantitative engine of that paper: the concatenation recursion, the error-model bounds it is applied to, the threshold it produces, and the resulting overhead.

Setting

Fix a physical implementation in which each gate, state preparation and idle memory step carries an error location. A noisy network with error locations indexed by a finite set locs\mathrm{locs}locs is described through an error expansion; each summand fails at a set of locations, and the model constrains how probable simultaneous failures are.

Write μ\muμ for the underlying probability measure and FiF_iFi​ for the event that location iii has failed. Two models from Section I.B of the paper are used here.

  • Quasi-independent stochastic model with error probability ppp: for every set S⊆locsS \subseteq \mathrm{locs}S⊆locs of error locations,
μ(⋂i∈SFi)≤p∣S∣.\mu\Big(\bigcap_{i \in S} F_i\Big) \le p^{|S|}.μ(i∈S⋂​Fi​)≤p∣S∣.
  • Quasi-independent monotonic model with constant CCC and parameter ppp: the same quantity is bounded by C p∣S∣C\,p^{|S|}Cp∣S∣.

The fault tolerant construction of Section I turns physical qubits into encoded qubits: each encoded gate is a transversal operation preceded by fault tolerant recovery networks built from the 777-qubit Steane code. The analysis of Section II declares an encoded gate failed according to an explicit combinatorial rule, and shows that a failure requires two failures among a set of fff minimal pairs of physical error locations. Consequently one level of encoding maps a failure parameter ppp to fp2f p^{2}fp2.

Iterating this reduction is concatenation. Define the level-hhh failure parameter by

E0=p,Eh+1=f Eh2,E_0 = p, \qquad E_{h+1} = f\,E_h^{2},E0​=p,Eh+1​=fEh2​,

written Eh=levelError f p hE_h = \texttt{levelError}\,f\,p\,hEh​=levelErrorfph in the Lean development. Let KKK bound the number of qubits and time steps at the level below that are consumed by one encoded gate, so that one computational gate at level hhh costs KhK^{h}Kh physical resources.

Formalization targets

Goal — threshold theorem with polylogarithmic overhead

For f≥1f \ge 1f≥1, K≥2K \ge 2K≥2, 0≤p<1/f0 \le p < 1/f0≤p<1/f, and every gate count nnn and target failure probability qqq with 0<q<n0 < q < n0<q<n, there exists a level hhh with

n Eh<qandKh  ≤  K2(max⁡(1, log⁡(n/q)log⁡(1/(fp))))log⁡2K.n\,E_h < q \qquad\text{and}\qquad K^{h} \;\le\; K^{2}\Big(\max\Big(1,\ \frac{\log(n/q)}{\log\big(1/(fp)\big)}\Big)\Big)^{\log_2 K}.nEh​<qandKh≤K2(max(1, log(1/(fp))log(n/q)​))log2​K.

The first clause is resilience: a computation of nnn encoded gates can be made to fail with probability below any prescribed qqq. The second is the cost: the overhead per computational gate is polylogarithmic in n/qn/qn/q, with exponent log⁡2K\log_2 Klog2​K. The threshold itself is the hypothesis p<1/fp < 1/fp<1/f; no numerical value is hard-coded, so the statement survives any improvement of the paper's count f≤337195f \le 337195f≤337195 and of the resulting bounds 3.0×10−63.0 \times 10^{-6}3.0×10−6 and 1.3×10−61.3 \times 10^{-6}1.3×10−6.

Supporting targets

The milestone list covers the closed form Eh=f2h−1p2hE_h = f^{2^{h}-1} p^{2^{h}}Eh​=f2h−1p2h of the concatenation recursion and its threshold form Eh=(fp)2h/fE_h = (fp)^{2^{h}}/fEh​=(fp)2h/f; convergence Eh→0E_h \to 0Eh​→0 below threshold; the failure bounds ∣locs∣ p|\mathrm{locs}|\,p∣locs∣p and C ∣locs∣ pC\,|\mathrm{locs}|\,pC∣locs∣p for a network under the two error models; the criterion 2h>log⁡(n/q)/log⁡(1/(fp))2^{h} > \log(n/q)/\log(1/(fp))2h>log(n/q)/log(1/(fp)) under which n(fp)2h<qn (fp)^{2^{h}} < qn(fp)2h<q; and the overhead estimate K⌈log⁡2L⌉≤KLlog⁡2KK^{\lceil \log_2 L\rceil} \le K L^{\log_2 K}K⌈log2​L⌉≤KLlog2​K.

Significance

The threshold theorem is the reason large-scale quantum computation is regarded as a physically meaningful goal rather than an idealization: it converts an engineering requirement that scales with the problem size into a fixed, size-independent accuracy target, and it fixes the currency — threshold value and overhead exponent — in which every later architecture is compared. The quasi-independent models matter separately: they show the conclusion does not depend on noise being stochastic, which is what rules out the objection that coherent over-rotations invalidate the stochastic analysis.

The paper's own argument is a physics-style derivation: a circuit-level construction, a combinatorial count of failure-inducing pairs of error locations, and a recursion. This mission formalizes the recursion, the probabilistic bounds and the resulting threshold and overhead statements. It does not formalize the circuit-level content — the 777-qubit code, the cat-state syndrome extraction networks, transversality, or the count f≤337195f \le 337195f≤337195 itself; these enter the formal statements as the parameters fff and KKK. No machine-checked proof of the full threshold theorem, in the sense of a verified fault tolerant circuit construction, is known to the author of this proposal; formalizing the quantitative layer is a prerequisite for one.

Difficulty

The obvious route to "errors vanish under concatenation" is to iterate the bound p↦fp2p \mapsto f p^{2}p↦fp2 and observe the doubly exponential decay. That step is the easy half, and the milestone list makes it explicit. What the naive argument silently assumes is exactly what Section II must establish: that the level-(h+1)(h+1)(h+1) failure events again obey a quasi-independent bound, with kkk simultaneous encoded failures bounded by (fp2)k(f p^{2})^{k}(fp2)k rather than merely each single failure bounded by fp2f p^{2}fp2. This is why the paper's failure declaration is defined backwards in time and demands disjoint pairs of contributing locations. A formalization that bounds only single-gate failure probabilities cannot close the induction.

The second difficulty is bookkeeping at the edges: the recursion is stated for a failure parameter that is a probability in one model and an operator strength in the other, and the overhead bound involves a real exponent log⁡2K\log_2 Klog2​K and a ceiling, whose interaction with the degenerate cases p=0p = 0p=0 and n≤qn \le qn≤q has to be handled rather than assumed away.

Formalization scope

All parameters are real numbers: f,p,q,K,C∈Rf, p, q, K, C \in \mathbb{R}f,p,q,K,C∈R, and the gate count nnn is a natural number cast to R\mathbb{R}R. levelError is defined by recursion on the level and is total, so it is meaningful for every real f,pf, pf,p; the threshold hypothesis appears explicitly as p<1/fp < 1/fp<1/f where it is needed.

The two error models are predicates on a measure μ\muμ on a measurable space, a Finset of error locations, and a family of failure events, with values in ENNReal; they quantify over all subsets SSS of the location set, matching the paper's "at a given kkk many error locations". The network failure bounds are stated for μ(⋃iFi)\mu\big(\bigcup_{i} F_i\big)μ(⋃i​Fi​), the probability that some location fails.

The goal is not vacuous: the hypotheses f≥1f \ge 1f≥1, K≥2K \ge 2K≥2, 0≤p<1/f0 \le p < 1/f0≤p<1/f, 0<q<n0 < q < n0<q<n are simultaneously satisfiable (for instance f=1f = 1f=1, K=2K = 2K=2, p=10−6p = 10^{-6}p=10−6, n=106n = 10^{6}n=106, q=1/2q = 1/2q=1/2), and the conclusion asserts a strict inequality together with a quantitative overhead bound, so it cannot be discharged by a degenerate reading. Degenerate inputs are inside the statement rather than excluded: p=0p = 0p=0, and the case log⁡(n/q)/log⁡(1/(fp))≤1\log(n/q)/\log(1/(fp)) \le 1log(n/q)/log(1/(fp))≤1, are covered by the max⁡\maxmax and must be handled.

A complete development needs only Mathlib's real analysis (Real.log, Real.logb, Real.rpow, Nat.ceil) and the measure-theoretic union bound; nothing quantum is required by the formal statements. Contributions strengthening the milestone list towards the circuit level — a formal model of error locations in a network, of the induced next-level error model, or of stabilizer-code recovery — are welcome as further targets.

Selected references

  • E. Knill, R. Laflamme, W. H. Zurek, Resilient Quantum Computation: Error Models and Thresholds, arXiv:quant-ph/9702058 (1997), https://arxiv.org/abs/quant-ph/9702058
  • A. M. Steane, Error Correcting Codes in Quantum Theory, Physical Review Letters 77, 793 (1996), https://doi.org/10.1103/PhysRevLett.77.793
  • P. W. Shor, Fault-tolerant quantum computation, Proceedings of the 37th Symposium on Foundations of Computer Science (1996), https://arxiv.org/abs/quant-ph/9605011
  • D. Aharonov, M. Ben-Or, Fault-Tolerant Quantum Computation With Constant Error, arXiv:quant-ph/9611025 (1996), https://arxiv.org/abs/quant-ph/9611025
  • A. Yu. Kitaev, Quantum computations: algorithms and error correction, Russian Mathematical Surveys 52, 1191 (1997), https://doi.org/10.1070/RM1997v052n06ABEH002155
9 thms2 active usersReviewed
🏆Completed
Quantum Error Correction·Captain: Lucas

Aharonov-Ben-Or Threshold: Recursive Sparseness Below the ThresholdResearch Paper

Motivation

Quantum computations are performed on physical devices whose elementary operations are never exact: every gate, every idle qubit, every time step is subject to decoherence, damping and systematic inaccuracy. Without protection these errors accumulate and destroy the computation. In 1997 Dorit Aharonov and Michael Ben-Or proved that this can be overcome: as long as the error rate — the probability, per location per time step, that a fault occurs — is below a fixed positive threshold, an arbitrary quantum circuit can be simulated reliably with only polylogarithmic overhead (arXiv:quant-ph/9906129). This threshold theorem is the reason large-scale quantum computing is regarded as physically possible in principle, and the number it produces is the target every hardware programme is measured against.

The engine of the proof is not quantum mechanical but combinatorial. A circuit is encoded, then the encoded circuit is encoded again, rrr times over. The locations of the resulting circuit are organised into a tree of nested rectangles, and the computation survives precisely when the set of faulty locations is sparse in a recursive sense. The mission formalizes that engine: the recursive sparseness calculus of §7.4–7.6 and §8.2 of the paper, which converts a constant error rate below the threshold into a doubly exponentially small failure probability.

Setting

Fix two natural numbers: AAA, the number of locations in a rectangle, and kkk, the number of faulty sub-rectangles a rectangle tolerates (in the paper k=⌊d/l⌋k=\lfloor d/l\rfloork=⌊d/l⌋, where ddd is the number of errors the code corrects and lll is the spread of a fault). A 000-rectangle is a single location of the circuit, and an (r+1)(r+1)(r+1)-rectangle consists of AAA many rrr-rectangles. A fault pattern of an rrr-rectangle records, for each of its ArA^{r}Ar locations, whether a fault occurred there.

Following Definition 18 of the paper, a fault pattern is (r,k)(r,k)(r,k)-sparse by recursion on rrr: at level 000 it is sparse when the location carries no fault; at level r+1r+1r+1 it is sparse when at most kkk of the AAA constituent rrr-rectangles carry a pattern that is not (r,k)(r,k)(r,k)-sparse. Under independent probabilistic noise of rate η\etaη each location is faulty with probability η\etaη, independently; write P(r)P(r)P(r) for the probability that the pattern of an rrr-rectangle is (r,k)(r,k)(r,k)-sparse, and 1−P(r)1-P(r)1−P(r) for the probability that it is not.

Two thresholds appear. The threshold condition for probabilistic noise (Definition 19) is

(Ak+1) ηk+1<η,\binom{A}{k+1}\,\eta^{k+1} < \eta ,(k+1A​)ηk+1<η,

and the threshold condition for general noise (Definition 21) is

e(Ak+1)(2η)k+1<2η.e\binom{A}{k+1}(2\eta)^{k+1} < 2\eta .e(k+1A​)(2η)k+1<2η.

The corresponding thresholds are ηc=(Ak+1)−1/k\eta_c=\binom{A}{k+1}^{-1/k}ηc​=(k+1A​)−1/k and ηc′=12(e(Ak+1))−1/k\eta_c'=\tfrac12\bigl(e\binom{A}{k+1}\bigr)^{-1/k}ηc′​=21​(e(k+1A​))−1/k: every positive η\etaη below them satisfies the respective condition.

Target

The goal is the quantitative core of Theorem 10. Below the threshold there are constants C,c>0C,c>0C,c>0 such that for every circuit size v≥1v\ge 1v≥1 and every accuracy 0<ε<10<\varepsilon<10<ε<1 some recursion depth rrr satisfies both

v⋅(1−P(r))<εandAr≤C(1+log⁡vε)c.v\cdot\bigl(1-P(r)\bigr) < \varepsilon \qquad\text{and}\qquad A^{r} \le C\Bigl(1+\log\frac{v}{\varepsilon}\Bigr)^{c}.v⋅(1−P(r))<εandAr≤C(1+logεv​)c.

The first inequality says that, over all vvv rectangles of the simulated circuit, the probability that any of them carries a non-sparse fault pattern is below ε\varepsilonε; the second says that the overhead per simulated location, ArA^{r}Ar, is polylogarithmic in v/εv/\varepsilonv/ε — the cost claimed in Theorem 10.

The milestones build up to it: the two threshold conditions hold below the respective thresholds, the exponent gap δ\deltaδ of equation 7.5 exists, the base case P(0)=1−ηP(0)=1-\etaP(0)=1−η, the union-bound recursion 1−P(r+1)≤(Ak+1)(1−P(r))k+11-P(r+1)\le\binom{A}{k+1}(1-P(r))^{k+1}1−P(r+1)≤(k+1A​)(1−P(r))k+1, Lemma 10 itself (P(r)≥1−η(1+δ)rP(r)\ge 1-\eta^{(1+\delta)^r}P(r)≥1−η(1+δ)r), and the analogous recursion 8.8–8.10 governing the bad part in the general-noise analysis.

Significance

Lemma 10 is where a constant error rate becomes a vanishing one. Each level of concatenation only has to improve reliability slightly, because the improvement compounds doubly exponentially in the number of levels; this is exactly why constant-rate noise can be tolerated while Shor's earlier scheme needed a polylogarithmically small rate. The same recursion, with 2η2\eta2η in place of η\etaη and an extra factor of eee, controls the trace norm of the bad part of the fault-path expansion in the general-noise model (§8.2), which is how the result is extended from probabilistic errors to decoherence, damping and systematic inaccuracy.

Formalizing this part of the paper yields a reusable, model-independent statement: nothing in the milestones mentions density matrices, so the sparseness calculus can later be plugged into a formal model of noisy quantum circuits, or reused for classical concatenated fault tolerance. What this mission does not do is formalize the circuit-level statement of Theorem 10 — that requires a full Lean model of quantum circuits with mixed states, fault-tolerant procedures, and the code constructions of §3–§6, and is deliberately out of scope here.

Difficulty

The obvious approach — bound the failure probability by a union bound over all fault sets of size k+1k+1k+1 and iterate — is exactly right, but the iteration is where the work is. The recursion qr+1≤(Ak+1)qrk+1q_{r+1}\le\binom{A}{k+1}q_r^{k+1}qr+1​≤(k+1A​)qrk+1​ only produces the doubly exponential decay η(1+δ)r\eta^{(1+\delta)^r}η(1+δ)r once the threshold condition is converted into the strict exponent gap δ\deltaδ of equation 7.5, and the induction has to keep qr≤ηq_r\le\etaqr​≤η alive to reuse that gap at every level. In the general-noise recursion there is a second trap: the factor (1+br)A−k−1(1+b_r)^{A-k-1}(1+br​)A−k−1 exceeds 111, so the recursion is not contracting termwise, and one needs 2ηA≤12\eta A\le 12ηA≤1 to bound it by eee.

On the probabilistic side, the statements are about a genuine product distribution on an exponentially large sample space, not a recursively defined real number: the union-bound step must be carried out over the actual set of fault patterns, where independence across the AAA sub-rectangles is what makes the product bound valid.

Formalization scope

Fault patterns are a type defined by recursion on the level: FaultPattern(A,0)\mathrm{FaultPattern}(A,0)FaultPattern(A,0) is the booleans and FaultPattern(A,r+1)\mathrm{FaultPattern}(A,r+1)FaultPattern(A,r+1) is the functions from the AAA sub-rectangles to FaultPattern(A,r)\mathrm{FaultPattern}(A,r)FaultPattern(A,r); every level is a finite type with decidable equality. Sparseness is a decidable predicate following Definition 18 verbatim. The noise distribution is given by the weight of a pattern, the product over all locations of η\etaη (faulty) or 1−η1-\eta1−η (clean), and the sparseness probability is the sum of these weights over the sparse patterns — an honest finite product measure, so no hypothesis 0≤η≤10\le\eta\le10≤η≤1 is built into the definitions and each statement carries the range hypotheses it needs.

Two conventions are worth flagging. First, the thresholds are formalized as (Ak+1)−1/k\binom{A}{k+1}^{-1/k}(k+1A​)−1/k and 12(e(Ak+1))−1/k\tfrac12(e\binom{A}{k+1})^{-1/k}21​(e(k+1A​))−1/k; the paper prints the exponent as −k-k−k in equations 7.4 and 8.6, which is inconsistent with its own threshold conditions 7.3 and 8.5, and −1/k-1/k−1/k is the exponent for which the paper's claim "any η<ηc\eta<\eta_cη<ηc​ satisfies the threshold condition" is true. Second, Lemma 10 is stated with ≥\ge≥ rather than the paper's strict >>>, because at r=0r=0r=0 the two sides are equal.

Nothing here is vacuous: all hypotheses are satisfiable — for instance A=7A=7A=7, k=1k=1k=1, η=10−3\eta=10^{-3}η=10−3 satisfies every hypothesis of the goal — and the goal quantifies over all circuit sizes and all accuracies, so it cannot be discharged by a degenerate choice. A complete development needs only real analysis, finite sums and binomial counting from Mathlib. Contributions that generalise the branching factor to a per-rectangle bound Ai≤AA_i\le AAi​≤A, or that connect the calculus to a formal model of noisy circuits, are welcome.

Selected references

  • Dorit Aharonov and Michael Ben-Or, Fault-Tolerant Quantum Computation With Constant Error Rate, arXiv:quant-ph/9906129 (1999); preliminary version in Proc. 29th ACM STOC (1997). https://arxiv.org/abs/quant-ph/9906129
  • Peter W. Shor, Fault-tolerant quantum computation, Proc. 37th FOCS (1996). https://arxiv.org/abs/quant-ph/9605011
  • Dorit Aharonov, Alexei Kitaev and Noam Nisan, Quantum circuits with mixed states, Proc. 30th ACM STOC (1998). https://arxiv.org/abs/quant-ph/9806029
9 thms2 active usersReviewed
🏆Completed
Algebraic TopologyQuantum Error Correction·Captain: Rui Chao

Counterexamples to the Zeng–Pryadko Homological-Distance ConjectureResearch Paper

Background and main question

The tensor product of chain complexes is a basic construction in homological algebra and in the theory of quantum CSS codes. If AAA and BBB are finite based chain complexes over a finite field, their tensor product is graded by total degree,

(A⊗B)j=⨁i=0jAi⊗Bj−i.(A\otimes B)_j=\bigoplus_{i=0}^{j} A_i\otimes B_{j-i}.(A⊗B)j​=i=0⨁j​Ai​⊗Bj−i​.

Each complex carries a basis-dependent homological distance: dj(A)d_j(A)dj​(A) is the least Hamming weight of a degree-jjj cycle that is not a boundary, with dj(A)=∞d_j(A)=\inftydj​(A)=∞ when the degree-jjj homology vanishes. A natural candidate for the distance of the tensor product is therefore

mj(A,B)=min⁡0≤i≤jdi(A)dj−i(B).m_j(A,B)=\min_{0\le i\le j} d_i(A)d_{j-i}(B).mj​(A,B)=0≤i≤jmin​di​(A)dj−i​(B).

Every nonzero pair of homology classes in complementary degrees gives a pure-tensor class of the corresponding product weight, so one always has dj(A⊗B)≤mj(A,B)d_j(A\otimes B)\le m_j(A,B)dj​(A⊗B)≤mj​(A,B). The substantive question is whether this upper bound is always sharp.

In the preprint Higher-Dimensional Quantum Hypergraph-Product Codes with Finite Rates, posted in 2018 and subsequently published in Physical Review Letters, Zeng and Pryadko proved that the bound is sharp when one factor is a binary one-complex. Their exact formula is Eq. (13) of the arXiv version and is the capstone theorem of the Prove2Me mission Higher-Dimensional Quantum Hypergraph-Product Codes with Finite Rates. The present formalization is a direct sequel: it retains the same definitions and degree conventions while examining what happens when the one-complex restriction is removed.

Zeng and Pryadko later considered arbitrary finite chain complexes over arbitrary finite fields in Minimal distances for certain quantum product codes and tensor products of chain complexes, published in 2020. In the corresponding arXiv preprint, Conjecture 18 asserts the unrestricted equality

dj(A⊗B)=mj(A,B).d_j(A\otimes B)=m_j(A,B).dj​(A⊗B)=mj​(A,B).

The conjecture is appealing because it would make tensor-product distance completely compositional: the degreewise distances of the two factors would determine the distance of their product. The obstruction is that a homology class in A⊗BA\otimes BA⊗B need not have a minimum-weight representative supported in a single bidegree. A representative spread across several summands of the total complex can be lighter than every pure-tensor representative. The purpose of this formalization is to turn that observation into a concrete, machine-checked counterexample to Conjecture 18 while preserving the valid one-complex theorem as a sharply delimited special case.

The counterexample mechanism

The common foundation is recorded in the Prove2Me entry Based binary chain complexes and homological distance. In particular, the boundary in degree jjj is a map ∂j:Aj→Aj−1\partial_j:A_j\to A_{j-1}∂j​:Aj​→Aj−1​, and

dj(A)=inf⁡{wt⁡(x):x∈ker⁡∂j, x∉im⁡∂j+1}.d_j(A)=\inf\{\operatorname{wt}(x):x\in\ker\partial_j,\ x\notin\operatorname{im}\partial_{j+1}\}.dj​(A)=inf{wt(x):x∈ker∂j​, x∈/im∂j+1​}.

The one-complex result is separately available as Eq. (13) — Exact distance with a one-complex. The construction below uses exactly the same notion of distance, but both tensor factors are genuine three-term complexes.

Begin with binary CSS check maps

HX:F2n⟶F2rX,HZ:F2n⟶F2rZ,H_X:\mathbb F_2^n\longrightarrow\mathbb F_2^{r_X}, \qquad H_Z:\mathbb F_2^n\longrightarrow\mathbb F_2^{r_Z},HX​:F2n​⟶F2rX​​,HZ​:F2n​⟶F2rZ​​,

assumed surjective and satisfying HXHZT=HZHXT=0H_XH_Z^T=H_ZH_X^T=0HX​HZT​=HZ​HXT​=0. Suppose there are logical vectors x,z∈F2nx,z\in\mathbb F_2^nx,z∈F2n​ such that

HZx=0,HXz=0,x⋅z=1.H_Zx=0,\qquad H_Xz=0,\qquad x\cdot z=1.HZ​x=0,HX​z=0,x⋅z=1.

The check maps determine two dual three-term complexes

A:F2rX←HXF2n←HZTF2rZ,B:F2rZ←HZF2n←HXTF2rX.A:\quad \mathbb F_2^{r_X}\xleftarrow{H_X}\mathbb F_2^n \xleftarrow{H_Z^T}\mathbb F_2^{r_Z}, \qquad B:\quad \mathbb F_2^{r_Z}\xleftarrow{H_Z}\mathbb F_2^n \xleftarrow{H_X^T}\mathbb F_2^{r_X}.A:F2rX​​HX​​F2n​HZT​​F2rZ​​,B:F2rZ​​HZ​​F2n​HXT​​F2rX​​.

Their degree-two tensor space has three bidegree summands, corresponding to (2,0)(2,0)(2,0), (1,1)(1,1)(1,1), and (0,2)(0,2)(0,2). Under the natural matrix identifications, consider the element whose three blocks are

(IrZ,In,IrX).(I_{r_Z},I_n,I_{r_X}).(IrZ​​,In​,IrX​​).

The CSS orthogonality relations make this element a cycle. Its pairing with the chosen logical vectors certifies that it is not a boundary. Its Hamming weight is exactly rZ+n+rXr_Z+n+r_XrZ​+n+rX​, whereas the componentwise candidate in degree two reduces to

m2(A,B)=d1(A)d1(B).m_2(A,B)=d_1(A)d_1(B).m2​(A,B)=d1​(A)d1​(B).

Consequently, any CSS datum satisfying

rZ+n+rX<d1(A)d1(B)r_Z+n+r_X<d_1(A)d_1(B)rZ​+n+rX​<d1​(A)d1​(B)

produces the strict inequality d2(A⊗B)<m2(A,B)d_2(A\otimes B)<m_2(A,B)d2​(A⊗B)<m2​(A,B). For orientation, a binary quantum Golay CSS presentation with parameters [[23,1,7]][[23,1,7]][[23,1,7]] has rX=rZ=11r_X=r_Z=11rX​=rZ​=11, giving the numerical comparison 45<4945<4945<49. The formal proof must supply an explicit CSS instance and verify its algebraic and distance properties, rather than relying on the parameter notation alone.

Formalization objectives

The first milestone proves the general certificate: for every CSS datum satisfying the hypotheses above, the element (IrZ,In,IrX)(I_{r_Z},I_n,I_{r_X})(IrZ​​,In​,IrX​​) is a nontrivial degree-two cycle of weight rZ+n+rXr_Z+n+r_XrZ​+n+rX​, and the componentwise minimum is d1(A)d1(B)d_1(A)d_1(B)d1​(A)d1​(B).

The second milestone constructs and verifies one explicit CSS datum for which rZ+n+rX<d1(A)d1(B)r_Z+n+r_X<d_1(A)d_1(B)rZ​+n+rX​<d1​(A)d1​(B). This is the step that turns the general mechanism into an actual counterexample.

The capstone packages the construction as the direct existential statement

∃ A,Bd2(A⊗B)<min⁡0≤i≤2di(A)d2−i(B).\exists\,A,B\qquad d_2(A\otimes B)<\min_{0\le i\le 2}d_i(A)d_{2-i}(B).∃A,Bd2​(A⊗B)<0≤i≤2min​di​(A)d2−i​(B).

Thus the final theorem is not conditional on the existence of suitable code data: it exhibits finite based binary chain complexes for which the equality proposed in Conjecture 18 fails.

Relation to prior work

The counterexample concerns only the unrestricted passage from a one-complex factor to two arbitrary bounded complexes. It does not conflict with Zeng and Pryadko's Eq. (13), whose one-complex hypothesis rules out the three-bidegree interaction used here.

The broader literature also indicates why additional structure matters. Bravyi and Hastings introduced homological-product codes and analyzed logical representatives in product constructions; Audoux and Couvreur developed tensor products of CSS codes through chain-complex methods. More recently, Akhmechet et al. discussed the Zeng–Pryadko conjecture in the structured setting of complexes derived from Khovanov homology, while Berthusen et al. restated it as Conjecture 5.1 in their study of automorphism gadgets. Golowich and Guruswami obtained strong distance guarantees for iterated homological products under expansion and local-testability hypotheses. These results are compatible with the proposed counterexample: they concern special families or impose hypotheses that are absent from Conjecture 18.

Besides settling the unrestricted statement, the formalization isolates a reusable obstruction. It shows precisely how a low-weight class assembled across several bidegrees can evade a formula based only on the degreewise distances of the factors. This distinction should help guide corrected formulations in which an exact product formula, or a useful lower bound, is recovered from additional geometric, expansion, or local-testability assumptions.

The main formal issues are mathematically substantive rather than presentational: the proof must respect the endpoint conventions for the boundary maps, distinguish a cycle from a non-boundary, compare finite weights with ∞\infty∞-valued distances, and verify a concrete strict-gap instance. Making each of these points explicit is especially important here, because an indexing shift or a merely conditional existence statement would no longer constitute a refutation of the conjecture as stated.

References

  • W. Zeng and L. P. Pryadko, Higher-Dimensional Quantum Hypergraph-Product Codes with Finite Rates, Physical Review Letters 122, 230501 (2019); arXiv:1810.01519 (2018), Eq. (13).
  • W. Zeng and L. P. Pryadko, Minimal distances for certain quantum product codes and tensor products of chain complexes, Physical Review A 102, 062402 (2020); arXiv:2007.12152, Conjecture 18.
  • S. Bravyi and M. B. Hastings, Homological Product Codes, STOC 2014; arXiv:1311.0885.
  • B. Audoux and A. Couvreur, On tensor products of CSS codes, Annales de l'Institut Henri Poincaré D 6 (2019); arXiv:1512.07081.
  • R. Akhmechet et al., Khovanov homology and quantum error-correcting codes, arXiv:2410.11252 (2024).
  • N. Berthusen et al., Automorphism gadgets in homological product codes, arXiv:2508.04794 (2025).
  • L. Golowich and V. Guruswami, Quantum LDPC Codes of Almost Linear Distance via Iterated Homological Products, CCC 2025; full version.
6 thms2 active usersReviewed
🏆Completed
Complexity Theory·Captain: Goku

Shallow Quantum Circuits and Causal ConesTextbook

Motivation

A quantum circuit of depth ddd built from gates of fan-in at most two cannot let an output wire depend on more than 2d2^d2d input wires. The argument is folklore and takes a paragraph on paper: the causal cone of the measured wire grows by at most a factor of two per layer. Formalizing it exposes a subtlety that the paper argument hides, and that is what this mission is about.

The subtlety

Define the backward cone step of a wire set SSS through a layer lll by adjoining the support of every gate of lll that meets SSS. There is a choice here: test each gate against the incoming set SSS, or against the partially accumulated cone. Testing against the accumulator over-approximates, and the doubling bound fails. Testing against the incoming set gives the bound — but is only correct when the gates within a layer act on pairwise disjoint wires.

Without that hypothesis (LayerOk) the semantic statement is false, and the counterexample is small: on three wires, the single layer [cnot 2 1, cnot 1 0] has cone {0,1}\{0,1\}{0,1} around wire 000, yet wire 000 ends up holding x0⊕x1⊕x2x_0 \oplus x_1 \oplus x_2x0​⊕x1​⊕x2​. So the cone under-approximates the true dependence. This mission's development carries LayerOk throughout, and the counterexample is recorded in the source.

What is formalized

Layered circuits over {H,S,T,CNOT}\{H, S, T, \mathrm{CNOT}\}{H,S,T,CNOT} on nnn wires, with states as amplitude functions on bit-strings and no tensor products anywhere. On top of that:

  • the combinatorial half — one layer at most doubles the cone, hence ∣cone∣≤2d ∣S∣|\mathrm{cone}| \le 2^{d}\,|S|∣cone∣≤2d∣S∣;
  • norm preservation, so that acceptProb is a genuine probability in [0,1][0,1][0,1];
  • the semantic half — inputs agreeing on the causal cone of the output wire are accepted with equal probability.

The semantic half is proved in the Heisenberg picture. The measurement observable is conjugated backwards through the circuit and its support tracked: a gate meeting the support enlarges it by that gate's own wires, and a gate missing it commutes with the observable and cancels against its own adjoint. That cancellation is the reason the non-cascading cone step is correct, and it is why unitarity of the gate set is needed at the 2n2^n2n-dimensional level rather than gate by gate. Supporting this is a small reusable algebra of local operators: locality is monotone, closed under adjoint and product, and disjointly supported operators commute.

The frontier

The published depth bound assumes each input wire lies in the syntactic cone of the output. That is weaker than saying the wire matters. Milestone 1 asks for the semantically honest version, stated in terms of genuine functional dependence; the bridge is the semantic cone theorem already in the development.

Beyond that, the natural continuations are the same argument for fan-in-kkk gates (∣cone∣≤kd|\mathrm{cone}| \le k^{d}∣cone∣≤kd), for geometrically local circuits where cone growth is linear rather than exponential, and ultimately the Bravyi–Gosset–König separation QNC0⊄NC0\mathrm{QNC}^{0} \not\subset \mathrm{NC}^{0}QNC0⊂NC0 — which needs machinery (non-local games, magic squares) that this development deliberately does not build.

9 thms2 active usersReviewed
🏆Completed
Complexity Theory·Captain: thomasvidick

MIP* = REResearch Paper

Why the value of a nonlocal game matters

In a two-player one-round nonlocal game, a verifier samples a pair of questions, sends one to each of two players who cannot communicate, and accepts or rejects the pair of answers. The players may share an arbitrary quantum state and measure it; the entangled value of the game is the supremum of their acceptance probability over all such strategies. Nonlocal games are the operational form of Bell's theorem, the model behind device-independent cryptography and self-testing, and the engine of the complexity class MIP∗\mathrm{MIP}^*MIP∗ of languages decidable by a polynomial-time verifier interacting with entangled provers (Cleve, Høyer, Toner, Watrous, 2004, arXiv:quant-ph/0404076).

With classical provers, multiprover interactive proofs decide exactly NEXP\mathrm{NEXP}NEXP (Babai, Fortnow, Lund, 1991, doi:10.1007/BF01200056). For entangled provers the lower bounds grew in steps: NEXP⊆MIP∗\mathrm{NEXP} \subseteq \mathrm{MIP}^*NEXP⊆MIP∗ (Ito, Vidick, 2012, arXiv:1207.0550), NEEXP⊆MIP∗\mathrm{NEEXP} \subseteq \mathrm{MIP}^*NEEXP⊆MIP∗ (Natarajan, Wright, 2019, arXiv:1904.05870), and finally MIP∗=RE\mathrm{MIP}^* = \mathrm{RE}MIP∗=RE (Ji, Natarajan, Vidick, Wright, Yuen, 2020, arXiv:2001.04383): every recursively enumerable language, the halting problem included, has a two-prover one-round interactive proof with entangled provers, and the entangled value of a game cannot be approximated by any algorithm. Through the equivalences of Fritz (2012, arXiv:1008.1168), Junge, Navascués, Palazuelos, Pérez-García, Scholz, Werner (2011, arXiv:1008.1142) and Ozawa (2013, arXiv:1212.1700), the theorem gives a negative answer to Tsirelson's problem and to Connes' embedding problem for von Neumann algebras.

Synchronous games, their value, and game descriptions

Everything below is defined from scratch; the formal development uses exactly these objects.

A synchronous game G=(X,A,μ,D)G = (X, A, \mu, D)G=(X,A,μ,D) has a finite question alphabet XXX and a finite answer alphabet AAA, shared by both players, a probability distribution μ\muμ on X×XX \times XX×X, and a decision predicate D(x,y,a,b)∈{0,1}D(x, y, a, b) \in \{0, 1\}D(x,y,a,b)∈{0,1} which, on equal questions x=yx = yx=y, rejects every pair of unequal answers a≠ba \ne ba=b.

A synchronous strategy of dimension d≥1d \ge 1d≥1 assigns to each question xxx a projective measurement (Max)a∈A(M^x_a)_{a \in A}(Max​)a∈A​ on Cd\mathbb{C}^dCd: positive semidefinite matrices summing to the identity, each idempotent. Both players use the same measurements, and probabilities are taken in the normalized trace τ(M)=Tr⁡(M)/d\tau(M) = \operatorname{Tr}(M)/dτ(M)=Tr(M)/d. The value of the strategy is

ω(G,M)  =  ∑x,y∈X∑a,b∈Aμ(x,y) D(x,y,a,b) Tr⁡(MaxMby)d,\omega(G, M) \;=\; \sum_{x, y \in X} \sum_{a, b \in A} \mu(x, y)\, D(x, y, a, b)\, \frac{\operatorname{Tr}\bigl(M^x_a M^y_b\bigr)}{d},ω(G,M)=x,y∈X∑​a,b∈A∑​μ(x,y)D(x,y,a,b)dTr(Max​Mby​)​,

the acceptance probability when the players measure a maximally entangled state of dimension ddd. The synchronous value ωs(G)\omega_{\mathrm{s}}(G)ωs​(G) is the supremum of ω(G,M)\omega(G, M)ω(G,M) over all synchronous strategies of all dimensions. A synchronous strategy is a tensor-product strategy, so ωs(G)≤ω∗(G)\omega_{\mathrm{s}}(G) \le \omega^*(G)ωs​(G)≤ω∗(G), the entangled value.

A game description ggg is a finite object from which a synchronous game is read off: natural numbers nX,nAn_X, n_AnX​,nA​ fixing X={0,…,nX}X = \{0, \dots, n_X\}X={0,…,nX​} and A={0,…,nA}A = \{0, \dots, n_A\}A={0,…,nA​}, a finite list of weighted question pairs (x,y,w)(x, y, w)(x,y,w) and a finite list of accepted tuples (x,y,a,b)(x, y, a, b)(x,y,a,b). The game GgG_gGg​ has μ(x,y)\mu(x, y)μ(x,y) proportional to the total weight listed for (x,y)(x, y)(x,y) (the point mass at (0,0)(0, 0)(0,0) if the total weight is zero), and DDD accepts exactly the listed tuples, except that unequal answers to equal questions are rejected. Descriptions are coded as natural numbers, so a map from machines to descriptions can be computable.

Turing machines are Gödel numbers c∈Cc \in \mathcal{C}c∈C of partial recursive functions, an equivalent model of computation; "ccc halts on the empty input" means that ccc applied to 000 is defined.

Formalization targets

Goal: halting reduces to the synchronous game value

There is a computable map ggg from machines to game descriptions such that for every machine ccc,

c halts on the empty input  ⟹  ωs(Gg(c))=1,c does not halt  ⟹  ωs(Gg(c))≤12.c \text{ halts on the empty input} \;\Longrightarrow\; \omega_{\mathrm{s}}\bigl(G_{g(c)}\bigr) = 1, \qquad c \text{ does not halt} \;\Longrightarrow\; \omega_{\mathrm{s}}\bigl(G_{g(c)}\bigr) \le \tfrac12 .c halts on the empty input⟹ωs​(Gg(c)​)=1,c does not halt⟹ωs​(Gg(c)​)≤21​.

This is Theorem 12.2 of arXiv:2001.04383 in synchronous form, with the running time of ggg relaxed from polynomial to computable. It leaves the running time and the choice of value unfixed, and it is the weakest statement that already makes the game value uncomputable: no algorithm decides whether a game has synchronous value 111 or at most 12\tfrac1221​.

Stronger: polynomial time, and the entangled value

The paper's statement has ggg computable in polynomial time, and both the synchronous and the entangled value satisfy the two clauses. With polynomial time the goal gives RE⊆MIP∗\mathrm{RE} \subseteq \mathrm{MIP}^*RE⊆MIP∗ for the class MIP1,1/2∗(2,1)\mathrm{MIP}^*_{1, 1/2}(2, 1)MIP1,1/2∗​(2,1) of two-prover one-round proofs with completeness 111 and soundness 12\tfrac1221​, and with the elementary inclusion MIP∗⊆RE\mathrm{MIP}^* \subseteq \mathrm{RE}MIP∗⊆RE this is MIP∗=RE\mathrm{MIP}^* = \mathrm{RE}MIP∗=RE. These statements need a model of polynomial time and a definition of tensor-product strategies, and will be added as milestones once those definitions are published here.

Significance

The result. Approximating ω∗(G)\omega^*(G)ω∗(G) to within a constant is as hard as the halting problem. A hierarchy of semidefinite programs converges from above to the commuting-operator value ωco(G)\omega^{\mathrm{co}}(G)ωco(G), so the two values cannot coincide for every game: a negative answer to Tsirelson's problem and, by the equivalences cited above, to Connes' embedding problem, open since 1976. On the complexity side, MIP∗=RE\mathrm{MIP}^* = \mathrm{RE}MIP∗=RE is the exact characterization of the power of entangled provers.

Formalizing it. The theorem is proved and its proof is machine-checked: HaltingGameValue.halting_reduces_to_gameValue in the Lean 4 development vidick/MIPRE-formalization proves exactly the goal above over Mathlib, with propext, Classical.choice and Quot.sound as its only axioms. The development is about 1,200 modules and 410,000 lines, including vendored components and a separate library for the low individual degree test; its blueprint records every statement, the dependency graph and the correspondence with the paper. This mission makes the result a fixed, auditable target: the definitional layer above is written over Mathlib alone, so the goal can be checked from first principles, and the components of the proof whose developments are self-contained are published as milestones with their proofs, to be imported, reused and re-proved. The first milestone, the recursive compression lemma of Marks, Nezhadi and Yuen (Lemma 5.1 of their manuscript; blueprint node lem:recursive-compression), comes with a complete proof over a dependency tree of 22 definition bundles and 60 theorems: a cost model of computation, an efficient universal machine, the s-m-n and Kleene recursion theorems, succinct descriptions.

Difficulty

The obvious approach fails at once: a verifier runs in time polynomial in its input, while a machine may run for arbitrarily long before halting, so no fixed game can check halting by simulation. What makes RE\mathrm{RE}RE reachable is gap-preserving compression: a polynomial-time transformation of verifiers such that the compressed verifier, at index nnn, has value 111 when the original verifier at index 2n2^n2n does, and value at most 12\tfrac1221​ when the original does. The compressed game is exponentially cheaper than the game it certifies, and a fixed-point argument through Kleene's recursion theorem, in the abstract form of the recursive compression lemma, turns compression into the halting reduction. Compression rests on three components, each with a central ingredient that is a theorem in its own right: question reduction by introspection, in which the players sample their own questions and a Pauli basis self-test built on the quantum low individual degree test certifies that they did; answer reduction, in which the players prove through a probabilistically checkable proof that their answers would have been accepted; and gap amplification by parallel repetition. In each component soundness must hold against strategies of unbounded dimension, where no finite enumeration or compactness argument is available, and the quantitative losses must compose so that the threshold 12\tfrac1221​ is preserved exactly.

Formalization scope

The goal is stated with the synchronous value, projective measurements on Cd\mathbb{C}^dCd with d≥1d \ge 1d≥1, and the normalized trace in place of a state vector; the real part of the trace is taken in the strategy value, and the synchronous value is the supremum, in R\mathbb{R}R, over all strategies of all dimensions (the values lie in [0,1][0, 1][0,1]). Game descriptions are interpreted with the fallback distribution above, so that the interpretation is total. Machines are Mathlib's Nat.Partrec.Code and computability is Mathlib's Computable; polynomial time is not fixed by the goal.

No trivializing formalization is available: the game depends on ccc only through ggg, and the two clauses demand values 111 and at most 12\tfrac1221​ on complementary sets of machines, so no map that ignores its argument, or produces only descriptions of zero total weight, satisfies the statement.

A complete development needs four layers: (i) a model of computation with costs, with polynomial-time maps, an efficient universal machine, the s-m-n and Kleene recursion theorems and succinct descriptions, which is the recursive compression milestone's dependency tree and is reusable in any undecidability argument by compression; (ii) normal form verifiers, conditionally linear samplers, the Pauli basis test and the quantum low individual degree test; (iii) the three transformations, with answer reduction over a PCP decider built from a succinct Cook–Levin theorem and direct parallel repetition in value form, and their soundness; (iv) the halting reduction from compression and the passage from the tensor-product to the synchronous value. Statements and proofs of these components over the shared definitional layer are welcome, as are alternative proofs of the goal from the same definitions, the polynomial-time and entangled-value strengthenings, and the inclusion MIP∗⊆RE\mathrm{MIP}^* \subseteq \mathrm{RE}MIP∗⊆RE.

Selected references

  • Z. Ji, A. Natarajan, T. Vidick, J. Wright, H. Yuen, MIP* = RE, 2020 (revised 2022). arXiv:2001.04383
  • Z. Ji, A. Natarajan, T. Vidick, J. Wright, H. Yuen, Quantum soundness of the classical low individual degree test, 2020. arXiv:2009.12982
  • A. Marks, S. S. Nezhadi, H. Yuen, The recursive compression method for proving undecidability results, manuscript.
  • J. Lin, Tracially embeddable strategies: lifting MIP* tricks to MIP^co, 2023. arXiv:2304.01940
  • A. Natarajan, J. Wright, NEEXP ⊆ MIP*, FOCS 2019. arXiv:1904.05870
  • T. Ito, T. Vidick, A multi-prover interactive proof for NEXP sound against entangled provers, FOCS 2012. arXiv:1207.0550
  • R. Cleve, P. Høyer, B. Toner, J. Watrous, Consequences and limits of nonlocal strategies, CCC 2004. arXiv:quant-ph/0404076
  • L. Babai, L. Fortnow, C. Lund, Non-deterministic exponential time has two-prover interactive protocols, Computational Complexity 1 (1991). doi:10.1007/BF01200056
  • T. Fritz, Tsirelson's problem and Kirchberg's conjecture, Reviews in Mathematical Physics 24 (2012). arXiv:1008.1168
  • M. Junge, M. Navascués, C. Palazuelos, D. Pérez-García, V. B. Scholz, R. F. Werner, Connes' embedding problem and Tsirelson's problem, Journal of Mathematical Physics 52 (2011). arXiv:1008.1142
  • N. Ozawa, About the Connes embedding conjecture: algebraic approaches, Japanese Journal of Mathematics 8 (2013). arXiv:1212.1700
  • The Lean development: vidick/MIPRE-formalization (statement MIPRE/HaltingGameValue.lean, proof MIPRE/MainTheorem.lean) and its blueprint.
106 thms1 active userReviewed
PreviousPage 2 of 3Next

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me