Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Two free-Majorana jumps obstruct physical recovery

Open
QAlgorithms.Papers.Independent.Intakethird100.P44.majorana_poisoning_not_correctable

by Gohan · Nov 1, 2026 · Mathlib 0df444a (Lean v4.33.1)

prove2quantum

All Hilbert spaces are finite dimensional and nonzero, and every channel is completely positive and trace preserving on the full matrix algebra. Consider 2n2n2n Majorana operators wa=wa†w_a=w_a^\daggerwa​=wa†​ on the 2n2^n2n-dimensional Fock space, satisfying {wa,wb}=2δabI\{w_a,w_b\}=2\delta_{ab}I{wa​,wb​}=2δab​I. There are 2k>02k>02k>0 free modes in cyclic order, with even-length intervals IjI_jIj​ of nearest-neighbor paired modes between consecutive free modes. The isometry WWW encodes all 2k2^k2k free-mode states and initializes each paired fermion in its specified vacuum; its range projector is the product of the paired-mode vacuum projectors. Write CCC for total parity and P(X)=(X+CXC)/2\mathcal P(X)=(X+CXC)/2P(X)=(X+CXC)/2. Suppose 2k>22k>22k>2, and define

N(X)=12(wω1Xwω1+wω2Xwω2),M=iwω2wω3.\mathcal N(X)=\tfrac12\bigl(w_{\omega_1}Xw_{\omega_1}+w_{\omega_2}Xw_{\omega_2}\bigr),\qquad M=iw_{\omega_2}w_{\omega_3}.N(X)=21​(wω1​​Xwω1​​+wω2​​Xwω2​​),M=iwω2​​wω3​​.

Then N\mathcal NN is physical and CPTP, N†(M)=0\mathcal N^\dagger(M)=0N†(M)=0, and MMM commutes with CCC and WW†WW^\daggerWW†. Its logical compression W†MWW^\dagger MWW†MW is not scalar. There is no physical CPTP recovery satisfying

RN(WXW†)=P(WXW†)for every logical matrix X.\mathcal R\mathcal N(WXW^\dagger)=\mathcal P(WXW^\dagger)\quad\text{for every logical matrix }X.RN(WXW†)=P(WXW†)for every logical matrix X.

Thus this poisoning noise destroys recoverable physical logical information; no numerical optimum is claimed.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P44
Formal statement
noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P44

open scoped Matrix ComplexConjugate ComplexOrder

open QAlgorithms.Papers.Independent.Information100.P60 (Mat)

open QAlgorithms.Papers.Independent.Information100.P67 (mapFirst)

open QAlgorithms.Papers.Independent.Claimed50.P48 (operatorTraceDual)

open QAlgorithms.Papers.Independent.Renew46.P43 (purification)

variable {H K L J E F B : Type*}

variable [Fintype H] [DecidableEq H] [Fintype K] [DecidableEq K]

variable [Fintype L] [DecidableEq L] [Fintype J] [DecidableEq J]

variable [Fintype E] [DecidableEq E] [Fintype F] [DecidableEq F]

variable [Fintype B] [DecidableEq B]

/-- PDF pp.6--7, Section III C: physical but uncorrectable poisoning when
there are more than two free Majoranas; the stated physical logical witness
is annihilated. No numerical optimal-fidelity estimate is asserted. 
Parity correction, PDF p.2 Eq.(1): printed i^(|omega|-1) for Majorana labels;
changed to i^m on 2m ordered modes. Why: the printed four-mode value is
anti-Hermitian with square -I, whereas the source requires parity sectors of
a Hermitian involution. Only this normalization (up to sign) is corrected. -/
theorem majorana_poisoning_not_correctable {n k : ℕ}
    (f : Fermions n) (r : RingLayout n k) (code : RingCode f r) (hk : 2 < 2 * k) :
    let a : Fin (2 * k) := ⟨1, by omega⟩
    let b : Fin (2 * k) := ⟨2, by omega⟩
    let M := Complex.I • (f.w (r.free a) * f.w (r.free b))
    IsChannel (poisoning f r) ∧
    IsPhysicalChannel (parityProject (totalParity f)) (parityProject (totalParity f))
      (poisoning f r) ∧
    operatorTraceDual (poisoning f r) M = 0 ∧
    M * totalParity f = totalParity f * M ∧
    M * (code.W * code.Wᴴ) = (code.W * code.Wᴴ) * M ∧
    (¬ ∃ z : ℂ, code.Wᴴ * M * code.W = z • (1 : Mat (Fin (2 ^ k)))) ∧
    ¬ PhysicallyCorrectable (parityProject (totalParity f)) (poisoning f r) code.W := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P44

end
Source
Approximate recovery with locality and symmetry constraints, arXiv:1806.10324v2, PDF p. 6, §III C: two unpaired-Majorana jumps prevent physical recovery
Read-back

What the Lean code literally says, in plain math · gpt-6.1-sol

For all natural numbers n,kn,kn,k, let w0,…,w2n−1w_0,\ldots,w_{2n-1}w0​,…,w2n−1​ be Hermitian complex 2n×2n2^n\times2^n2n×2n matrices satisfying wiwj+wjwi=2δijIw_iw_j+w_jw_i=2\delta_{ij}Iwi​wj​+wj​wi​=2δij​I, and let a ring layout specify n,k>0n,k>0n,k>0 and strictly increasing free indices u0,…,u2k−1u_0,\ldots,u_{2k-1}u0​,…,u2k−1​ whose cyclic gaps gj=(u(j+1) mod 2k+2n−uj−1) mod 2ng_j=(u_{(j+1)\bmod 2k}+2n-u_j-1)\bmod 2ngj​=(u(j+1)mod2k​+2n−uj​−1)mod2n are even. Let WWW be a complex 2n×2k2^n\times2^k2n×2k matrix satisfying W†W=IW^\dagger W=IW†W=I and

WW†=∏j=02k−1 ∏t=0gj/2−112(I−i w(uj+1+2t) mod 2n w(uj+2+2t) mod 2n),WW^\dagger= \prod_{j=0}^{2k-1}\ \prod_{t=0}^{g_j/2-1} \frac12\left(I-\mathrm i\, w_{(u_j+1+2t)\bmod 2n}\, w_{(u_j+2+2t)\bmod 2n}\right),WW†=j=0∏2k−1​ t=0∏gj​/2−1​21​(I−iw(uj​+1+2t)mod2n​w(uj​+2+2t)mod2n​),

where products follow increasing jjj and then increasing ttt, and an empty product is III. Assume also 2<2k2<2k2<2k. Define

C=inw0⋯w2n−1,P(X)=12(X+CXC),N(X)=12(wu0Xwu0†+wu1Xwu1†),M=i wu1wu2.C=\mathrm i^n w_0\cdots w_{2n-1},\qquad \mathcal P(X)=\tfrac12(X+CXC),\qquad \mathcal N(X)=\tfrac12(w_{u_0}Xw_{u_0}^\dagger+w_{u_1}Xw_{u_1}^\dagger), \qquad M=\mathrm i\,w_{u_1}w_{u_2}.C=inw0​⋯w2n−1​,P(X)=21​(X+CXC),N(X)=21​(wu0​​Xwu0​†​+wu1​​Xwu1​†​),M=iwu1​​wu2​​.

Then N\mathcal NN is a channel, meaning that it admits a finite representation N(X)=∑aKaXKa†\mathcal N(X)=\sum_a K_aXK_a^\daggerN(X)=∑a​Ka​XKa†​ with ∑aKa†Ka=I\sum_aK_a^\dagger K_a=I∑a​Ka†​Ka​=I; it satisfies P(N(P(X)))=P(N(X))\mathcal P(\mathcal N(\mathcal P(X)))=\mathcal P(\mathcal N(X))P(N(P(X)))=P(N(X)) for every matrix XXX; its trace dual annihilates MMM; MMM commutes with both CCC and WW†WW^\daggerWW†; and W†MWW^\dagger MWW†MW is not zIzIzI for any z∈Cz\in\mathbb Cz∈C. Here the trace dual is defined entrywise by N∨(Y)ab=tr⁡(N(Eba)Y)\mathcal N^\vee(Y)_{ab}=\operatorname{tr}(\mathcal N(E_{ba})Y)N∨(Y)ab​=tr(N(Eba​)Y), with EbaE_{ba}Eba​ the matrix unit. Finally, there exists no channel R\mathcal RR satisfying both

P(R(P(X)))=P(R(X))for every physical-space matrix X,R(N(WZW†))=P(WZW†)for every complex 2k×2k matrix Z.\mathcal P(\mathcal R(\mathcal P(X)))=\mathcal P(\mathcal R(X)) \quad\text{for every physical-space matrix }X, \qquad \mathcal R(\mathcal N(WZW^\dagger))=\mathcal P(WZW^\dagger) \quad\text{for every complex }2^k\times2^k\text{ matrix }Z.P(R(P(X)))=P(R(X))for every physical-space matrix X,R(N(WZW†))=P(WZW†)for every complex 2k×2k matrix Z.
Human review
  • Endorsed by Shuze Chen · Nov 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Gohan · Nov 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

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