Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Example 5 / Eq.(67): erasure zero-error cost

Open
QAlgorithms.Papers.Independent.Intakethird100.P65.erasure_zeroError

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

prove2quantum

All systems have finite positive dimension and fixed computational bases. Matrix order is Loewner order and logarithms are base two.

Let the input dimension be d≥1d\ge1d≥1, let 0≤p≤10\le p\le10≤p≤1, and let ∣e⟩|e\rangle∣e⟩ be a new orthogonal erasure flag. Extend Ep\mathcal E_pEp​ linearly by Ep(X)=(1−p)X+pTr⁡(X)∣e⟩⟨e∣\mathcal E_p(X)=(1-p)X+p\operatorname{Tr}(X)|e\rangle\langle e|Ep​(X)=(1−p)X+pTr(X)∣e⟩⟨e∣, with the data block embedded in the output. Then

SNS,0(Ep)=12log⁡2(d2(1−p)+p).S_{\rm NS,0}(\mathcal E_p)=\tfrac12\log_2(d^2(1-p)+p).SNS,0​(Ep​)=21​log2​(d2(1−p)+p).

This regularized zero-error cost includes complete erasure and dimension one. The extra dummy index in the printed witness is not used.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P65
Formal statement
noncomputable section

open Filter QAlgorithms

open scoped BigOperators Matrix ComplexConjugate ComplexOrder Topology

namespace QAlgorithms.Papers.Independent.Intakethird100.P65

attribute [local instance] Classical.propDecidable Classical.decEq

open QAlgorithms.Papers.Independent.Dependency50.P09 (transposeBob)

open QAlgorithms.Papers.Independent.Dependency50.P35 (maxRelative)

open QAlgorithms.Papers.Independent.Independent2.P01 (dampingKraus damping)

open QAlgorithms.Papers.Independent.Gap48.P33 (phaseFlipChannel)

open QAlgorithms.Papers.Independent.Information100.P06 (Below)

open QAlgorithms.Papers.Independent.Information100.P34 (choi)

open QAlgorithms.Papers.Independent.Information100.P62 (Stochastic)

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

open QAlgorithms.Papers.Independent.Information100.P67 (IsSubnormalized purifiedDistance)

open QAlgorithms.Papers.Independent.Information100.P69 (replacer permute)

open QAlgorithms.Papers.Independent.Information100.P71 (diamond channelPower channelRelative)

open QAlgorithms.Papers.Independent.Support50.P39 (quantumMutualInformation)

variable {A B C D R : Type} [Fintype A] [DecidableEq A]
  [Fintype B] [DecidableEq B] [Fintype C] [DecidableEq C]
  [Fintype D] [DecidableEq D] [Fintype R] [DecidableEq R]

variable [Nonempty A] [Nonempty B] [Nonempty C] [Nonempty D] [Nonempty R]

/-- Example 5, Eq.(67), PDF p.14: an orthogonal flag, with the dimension-one
case included. The erroneous extra index in Eq.(66) is not used. -/
theorem erasure_zeroError (p : ℝ) (h0 : 0 ≤ p) (h1 : p ≤ 1) :
    zeroErrorCost (flaggedErasureChannel A p) =
      Real.logb 2 ((Fintype.card A : ℝ)^2*(1-p) + p) / 2 := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P65

end
Source
Quantum Channel Simulation and the Channel's Smooth Max-Information, arXiv:1807.05354v2, PDF p. 14, Example 5 / Eq.(67): erasure zero-error cost
Read-back

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

For every finite type AAA with decidable equality and every real ppp with 0≤p≤10\le p\le10≤p≤1, put d=∣A∣d=|A|d=∣A∣ and define Ep:Mat⁡(A)→Mat⁡(A∪{⊥})\mathcal E_p:\operatorname{Mat}(A)\to\operatorname{Mat}(A\cup\{\bot\})Ep​:Mat(A)→Mat(A∪{⊥}) by Ep(X)a,b=(1−p)Xa,b\mathcal E_p(X)_{a,b}=(1-p)X_{a,b}Ep​(X)a,b​=(1−p)Xa,b​ for a,b∈Aa,b\in Aa,b∈A, Ep(X)⊥,⊥=ptr⁡X\mathcal E_p(X)_{\bot,\bot}=p\operatorname{tr}XEp​(X)⊥,⊥​=ptrX, and zero on the other entries. The assertion is C0(Ep)=12log⁡2(d2(1−p)+p)C_0(\mathcal E_p)=\tfrac12\log_2(d^2(1-p)+p)C0​(Ep​)=21​log2​(d2(1−p)+p). Here C0(T)C_0(T)C0​(T) is the real infimum of all rrr such that, for every s>0s>0s>0, eventually for all positive integers nnn, cns(T⊗n,0)/n≤r+sc_{\mathrm{ns}}(T^{\otimes n},0)/n\le r+scns​(T⊗n,0)/n≤r+s; tensor powers act by products of the individual matrix-unit entries. The quantity cns(H,e)c_{\mathrm{ns}}(H,e)cns​(H,e) is the real infimum of log⁡2m\log_2 mlog2​m over positive integers mmm whose infimum simulation error is at most eee. This error is the infimum of 12∥Wm(P)−H∥⋄\tfrac12\|\mathcal W_m(P)-H\|_\diamond21​∥Wm​(P)−H∥⋄​ over channels P:Mat⁡(I×{0,…,m−1})→Mat⁡({0,…,m−1}×O)P:\operatorname{Mat}(I\times\{0,\ldots,m-1\})\to\operatorname{Mat}(\{0,\ldots,m-1\}\times O)P:Mat(I×{0,…,m−1})→Mat({0,…,m−1}×O) whose output marginals each depend only on the corresponding input marginal, where Wm(P)(X)b,d=∑a,c<mP(X⊗∣a⟩⟨c∣)(a,b),(c,d)\mathcal W_m(P)(X)_{b,d}=\sum_{a,c<m}P(X\otimes|a\rangle\langle c|)_{(a,b),(c,d)}Wm​(P)(X)b,d​=∑a,c<m​P(X⊗∣a⟩⟨c∣)(a,b),(c,d)​. Channels have finite Kraus representations with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I, and ∥F∥⋄=sup⁡∥X∥1≤1∥(idI⊗F)(X)∥1\|F\|_\diamond=\sup_{\|X\|_1\le1}\|(\mathrm{id}_I\otimes F)(X)\|_1∥F∥⋄​=sup∥X∥1​≤1​∥(idI​⊗F)(X)∥1​, over arbitrary matrices on I×II\times II×I, with ∥X∥1=Re⁡tr⁡X†X\|X\|_1=\operatorname{Re}\operatorname{tr}\sqrt{X^\dagger X}∥X∥1​=RetrX†X​. Zero tolerance constrains the infimum error without asserting its attainment. Both endpoints of ppp are included, and no nonemptiness condition on AAA is displayed. Real infima use 000 when empty or unbounded below, and log⁡20=0\log_2 0=0log2​0=0.

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