Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Exact constrained recovery with positive conditional states

Open
QAlgorithms.Papers.Independent.Intakethird100.P44.exact_constrained_recovery

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. Let P1,…,PrP_1,\ldots,P_rP1​,…,Pr​ be nonzero orthogonal Hermitian sector projectors summing to III, let C=∑iciPiC=\sum_i c_iP_iC=∑i​ci​Pi​ have pairwise distinct real eigenvalues, and put P(X)=∑iPiXPi\mathcal P(X)=\sum_iP_iXP_iP(X)=∑i​Pi​XPi​. Let N(X)=∑nEnXEn†\mathcal N(X)=\sum_nE_nXE_n^\daggerN(X)=∑n​En​XEn†​ have a finite normalized Kraus family and satisfy PNP=PN\mathcal P\mathcal N\mathcal P=\mathcal P\mathcal NPNP=PN. Let WWW be an isometry from a nonzero logical space. Correctability means existence of a CPTP R\mathcal RR with PRP=PR\mathcal P\mathcal R\mathcal P=\mathcal P\mathcal RPRP=PR and RN(WXW†)=P(WXW†)\mathcal R\mathcal N(WXW^\dagger)=\mathcal P(WXW^\dagger)RN(WXW†)=P(WXW†) for every logical matrix XXX. The code is correctable if and only if there exist environment matrices σj∣i\sigma_{j|i}σj∣i​ such that

σj∣i⪰0,∑jTr⁡σj∣i=1for each i,\sigma_{j|i}\succeq0,\qquad\sum_j\operatorname{Tr}\sigma_{j|i}=1\quad\text{for each }i,σj∣i​⪰0,j∑​Trσj∣i​=1for each i, W†En†PjEmW=∑i⟨n∣σj∣i∣m⟩W†PiWfor all j,n,m.W^\dagger E_n^\dagger P_jE_mW=\sum_i\langle n|\sigma_{j|i}|m\rangle W^\dagger P_iW\quad\text{for all }j,n,m.W†En†​Pj​Em​W=i∑​⟨n∣σj∣i​∣m⟩W†Pi​Wfor all j,n,m.

This is a positive, normalized classical–quantum conditional-output criterion; it does not normalize each block separately.

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 p.5, Eqs.(18),(20): necessary AND sufficient, including positivity
and total conditional normalization. -/
theorem exact_constrained_recovery [Nonempty H] [Nonempty L] {r : ℕ}
    (s : Sectors H r) (A : J → Mat H) (hA : IsKrausFamily A)
    (hphys : IsPhysicalChannel (sectorProject s) (sectorProject s) (krausApply A))
    (W : Matrix H L ℂ) (hW : Wᴴ * W = 1) :
    PhysicallyCorrectable (sectorProject s) (krausApply A) W ↔
      ChargeStateCondition s A W := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P44

end
Source
Approximate recovery with locality and symmetry constraints, arXiv:1806.10324v2, PDF p. 5, §III B, Eqs.(18),(20): exact constrained recovery criterion
Read-back

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

Let H,LH,LH,L be nonempty finite index types with decidable equality, JJJ a finite index type, and rrr a natural number. Suppose Pi∈Mat⁡(H)P_i\in\operatorname{Mat}(H)Pi​∈Mat(H), i<ri<ri<r, are nonzero Hermitian projections, mutually orthogonal and summing to III, with pairwise distinct real values viv_ivi​. Put P(X)=∑iPiXPi\mathcal P(X)=\sum_iP_iXP_iP(X)=∑i​Pi​XPi​. Let An∈Mat⁡(H)A_n\in\operatorname{Mat}(H)An​∈Mat(H), n∈Jn\in Jn∈J, satisfy ∑nAn†An=I\sum_nA_n^\dagger A_n=I∑n​An†​An​=I, put N(X)=∑nAnXAn†\mathcal N(X)=\sum_nA_nXA_n^\daggerN(X)=∑n​An​XAn†​, and assume 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 XXX. Let WWW be an H×LH\times LH×L matrix with W†W=IW^\dagger W=IW†W=I. Then there exists a channel RRR, meaning R(Y)=∑aVaYVa†R(Y)=\sum_aV_aYV_a^\daggerR(Y)=∑a​Va​YVa†​ for some finite family with ∑aVa†Va=I\sum_aV_a^\dagger V_a=I∑a​Va†​Va​=I, satisfying

P(R(P(Y)))=P(R(Y))for every Y,R(N(WXW†))=P(WXW†)for every X∈Mat⁡(L),\mathcal P(R(\mathcal P(Y)))=\mathcal P(R(Y))\quad\text{for every }Y, \qquad R(\mathcal N(WXW^\dagger))=\mathcal P(WXW^\dagger)\quad\text{for every }X\in\operatorname{Mat}(L),P(R(P(Y)))=P(R(Y))for every Y,R(N(WXW†))=P(WXW†)for every X∈Mat(L),

if and only if there exist positive semidefinite matrices σij∈Mat⁡(J)\sigma_{ij}\in\operatorname{Mat}(J)σij​∈Mat(J) satisfying ∑jtr⁡(σij)=1\sum_j\operatorname{tr}(\sigma_{ij})=1∑j​tr(σij​)=1 for every i<ri<ri<r and W†An†PjAmW=∑i(σij)nmW†PiWW^\dagger A_n^\dagger P_jA_mW=\sum_i(\sigma_{ij})_{nm}W^\dagger P_iWW†An†​Pj​Am​W=∑i​(σij​)nm​W†Pi​W for every j<rj<rj<r and n,m∈Jn,m\in Jn,m∈J. The assumptions exclude r=0r=0r=0 and empty JJJ.

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