Prove2Me
Navigate
DiscoverCollectionsFormalpediaGalaxyBlogsUsersMomentumMy Missions+
Prove2Me
Log in
← Formalpedia

Theorem 4.21 — extremality criterion for unital channels

Open
WatrousTQI.unital_extreme_kraus_criterion

by Mrnaiq · Nov 11, 2026 · Mathlib 0df444a (Lean v4.33.1)

convexitymixed-unitaryquantum-channelsquantum-information

Let d,m>0d,m>0d,m>0, let Φ\PhiΦ be a trace-preserving completely positive map with Φ(Id)=Id\Phi(I_d)=I_dΦ(Id​)=Id​, and let A0,…,Am−1A_0,\ldots,A_{m-1}A0​,…,Am−1​ be a complex-linearly independent Kraus family representing Φ\PhiΦ on every operator. Then Φ\PhiΦ is an extreme point of the real convex set of all unital channels if and only if the m2m^2m2 pairs (Ab†Aa,AaAb†)(A_b^\dagger A_a,A_aA_b^\dagger)(Ab†​Aa​,Aa​Ab†​) are complex-linearly independent. Pairs of matrices encode the block-diagonal operators in equation (4.147) through the natural complex-linear isomorphism; both blocks and the order of their factors are retained.

Preamble
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_WatrousTQI_UnitalChannels
set_option autoImplicit false
noncomputable section
open scoped Matrix ComplexConjugate BigOperators
Formal statement
namespace WatrousTQI

theorem unital_extreme_kraus_criterion (d m : ℕ) (hd : 0 < d) (hm : 0 < m)
    (Φ : QAlgorithms.Superop (Fin d) (Fin d))
    (hΦ : QAlgorithms.IsChannel Φ)
    (hu : Φ (1 : Matrix (Fin d) (Fin d) ℂ) = 1)
    (A : Fin m → Matrix (Fin d) (Fin d) ℂ)
    (hLI : LinearIndependent ℂ A)
    (hrep : ∀ X, Φ X = QAlgorithms.krausApply A X) :
    Φ ∈ (unitalChannels d).extremePoints ℝ ↔
      LinearIndependent ℂ (fun ab : Fin m × Fin m =>
        ((A ab.2)ᴴ * A ab.1, A ab.1 * (A ab.2)ᴴ)) := by sorry

end WatrousTQI
Source
John Watrous, The Theory of Quantum Information, Cambridge University Press, 2018, Chapter 4. https://cs.uwaterloo.ca/~watrous/TQI/TQI.pdf; Theorem 4.21, equations (4.146)–(4.147), printed p.225 (PDF p.233).
Read-back

What the Lean code literally says, in plain math · GPT-6 (Codex)

For every pair of natural numbers d,md,md,m with 0<d0<d0<d and 0<m0<m0<m, every function Φ:Md(C)→Md(C)\Phi:M_d(\mathbb C)\to M_d(\mathbb C)Φ:Md​(C)→Md​(C), and every indexed family A=(Aj)j∈{0,…,m−1}A=(A_j)_{j\in\{0,\ldots,m-1\}}A=(Aj​)j∈{0,…,m−1}​ of complex d×dd\times dd×d matrices, assume that Φ\PhiΦ is a channel in the following precise sense: there exist a natural number nnn and matrices B0,…,Bn−1B_0,\ldots,B_{n-1}B0​,…,Bn−1​ such that ∑rBr∗Br=Id\sum_r B_r^*B_r=I_d∑r​Br∗​Br​=Id​ and Φ(X)=∑rBrXBr∗\Phi(X)=\sum_r B_rXB_r^*Φ(X)=∑r​Br​XBr∗​ for every complex d×dd\times dd×d matrix XXX; assume also that Φ(Id)=Id\Phi(I_d)=I_dΦ(Id​)=Id​, that AAA is complex linearly independent (if ∑jcjAj=0\sum_j c_jA_j=0∑j​cj​Aj​=0 for complex coefficients cjc_jcj​, then every cj=0c_j=0cj​=0), and that Φ(X)=∑jAjXAj∗\Phi(X)=\sum_j A_jXA_j^*Φ(X)=∑j​Aj​XAj∗​ for every XXX. Then Φ\PhiΦ is an extreme point over the real scalars of the set of all channels Ψ:Md(C)→Md(C)\Psi:M_d(\mathbb C)\to M_d(\mathbb C)Ψ:Md​(C)→Md​(C) satisfying Ψ(Id)=Id\Psi(I_d)=I_dΨ(Id​)=Id​, where “channel” has exactly the finite matrix representation and normalization just specified, if and only if the family of ordered pairs of matrices (Ab∗Aa, AaAb∗)(A_b^*A_a,\,A_aA_b^*)(Ab∗​Aa​,Aa​Ab∗​), indexed by every ordered pair (a,b)∈{0,…,m−1}2(a,b)\in\{0,\ldots,m-1\}^2(a,b)∈{0,…,m−1}2, is complex linearly independent in Md(C)×Md(C)M_d(\mathbb C)\times M_d(\mathbb C)Md​(C)×Md​(C). Here ∗^*∗ denotes conjugate transpose; extremality means that Φ\PhiΦ belongs to that set and, for any two members Ψ,Ω\Psi,\OmegaΨ,Ω of the set and any real ttt with 0<t<10<t<10<t<1, the pointwise equality Φ=tΨ+(1−t)Ω\Phi=t\Psi+(1-t)\OmegaΦ=tΨ+(1−t)Ω implies Ψ=Φ=Ω\Psi=\Phi=\OmegaΨ=Φ=Ω; and the linear independence on the right means that, for any complex coefficients cabc_{ab}cab​, the simultaneous equalities ∑a,bcabAb∗Aa=0\sum_{a,b}c_{ab}A_b^*A_a=0∑a,b​cab​Ab∗​Aa​=0 and ∑a,bcabAaAb∗=0\sum_{a,b}c_{ab}A_aA_b^*=0∑a,b​cab​Aa​Ab∗​=0 imply cab=0c_{ab}=0cab​=0 for every ordered pair. Both diagonal pairs a=ba=ba=b and both orders of distinct indices are included. The assumptions exclude d=0d=0d=0 and m=0m=0m=0; the definition of channel itself permits a representation indexed by n=0n=0n=0, with empty sums equal to zero, but such a representation cannot satisfy ∑rBr∗Br=Id\sum_r B_r^*B_r=I_d∑r​Br∗​Br​=Id​ when d>0d>0d>0. The conclusion is conditional on all the stated hypotheses, including the independence of AAA; it asserts nothing for families that fail them.

Human review
  • Endorsed by Shuze Chen · Nov 12, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Mrnaiq · Nov 12, 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 paperProve2Me Team
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me