Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Exact amplitude-damping phase output matrix

Proved
QAlgorithms.Papers.Independent.Intakethird100.P74.damping_phase_state

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

prove2quantum

For 0≤η≤10\le\eta\le10≤η≤1, let Nη\mathcal N_\etaNη​ be the qubit amplitude-damping channel with Kraus operators 1−η∣0⟩⟨1∣\sqrt{1-\eta}|0\rangle\langle1|1−η​∣0⟩⟨1∣ and ∣0⟩⟨0∣+η∣1⟩⟨1∣|0\rangle\langle0|+\sqrt\eta|1\rangle\langle1|∣0⟩⟨0∣+η​∣1⟩⟨1∣. Let ∣Φ⟩=(∣00⟩+∣11⟩)/2|\Phi\rangle=(|00\rangle+|11\rangle)/\sqrt2∣Φ⟩=(∣00⟩+∣11⟩)/2​, with reference CCC first, and modulate CCC by ZxZ^xZx for x∈{0,1}x\in\{0,1\}x∈{0,1}. In the ordered basis ∣00⟩,∣01⟩,∣10⟩,∣11⟩|00\rangle,|01\rangle,|10\rangle,|11\rangle∣00⟩,∣01⟩,∣10⟩,∣11⟩,

(idC⊗Nη)(ZCx∣Φ⟩⟨Φ∣ZCx)=(1200(−1)xη20000001−η20(−1)xη200η2).(\mathrm{id}_C\otimes\mathcal N_\eta)(Z_C^x|\Phi\rangle\langle\Phi|Z_C^x) =\begin{pmatrix} \frac12&0&0&\frac{(-1)^x\sqrt\eta}{2}\\ 0&0&0&0\\ 0&0&\frac{1-\eta}{2}&0\\ \frac{(-1)^x\sqrt\eta}{2}&0&0&\frac\eta2 \end{pmatrix}.(idC​⊗Nη​)(ZCx​∣Φ⟩⟨Φ∣ZCx​)=​21​002(−1)xη​​​0000​0021−η​0​2(−1)xη​​002η​​​.

This gives an exact symbolic phase-channel state, without inferring a numerical parameter interval from Figure 2.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P74
Formal statement
namespace QAlgorithms.Papers.Independent.Intakethird100.P74

open QAlgorithms

open scoped ComplexConjugate Matrix

noncomputable section

local instance (p : Prop) : Decidable p := Classical.propDecidable p

variable {B E S1 S2 : Type}

variable [Fintype B] [Fintype E] [Fintype S1] [Fintype S2]

variable [DecidableEq B] [DecidableEq E] [DecidableEq S1] [DecidableEq S2]

/-- PDF pp.10–11, Appendix C.C: the normalized symbolic phase output matrix, in
C-first computational order. The undefined gamma in an intermediate line is bypassed
using the printed Kraus operators and the final eta-only formula. No numerical key-rate
interval is inferred from this exact state identity. -/
theorem damping_phase_state (eta : ℝ) (he0 : 0 ≤ eta) (he1 : eta ≤ 1) (x : Bool) :
    dampingPhaseOutput eta x = Matrix.of (fun p q =>
      dampingPhaseMatrix eta x (cbIndex p) (cbIndex q)) := by sorry

end

end QAlgorithms.Papers.Independent.Intakethird100.P74
Source
Polar codes for private classical communication, arXiv:1203.5794v1, PDF p. 11, Appendix C.C, unnumbered: explicit phase-channel output matrices
Read-back

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

For every real η\etaη satisfying 0≤η≤10\leq\eta\leq10≤η≤1, including both endpoints, and every Boolean xxx, put sx=1s_x=1sx​=1 for x=falsex=\mathrm{false}x=false and sx=−1s_x=-1sx​=−1 for x=truex=\mathrm{true}x=true. In the Boolean basis ordered as false,true\mathrm{false},\mathrm{true}false,true, define

vx=∣false,false⟩+sx∣true,true⟩2,K0=(100η),K1=(01−η00),v_x=\frac{|\mathrm{false},\mathrm{false}\rangle+s_x|\mathrm{true},\mathrm{true}\rangle}{\sqrt2}, \qquad K_0=\begin{pmatrix}1&0\\0&\sqrt{\eta}\end{pmatrix}, \qquad K_1=\begin{pmatrix}0&\sqrt{1-\eta}\\0&0\end{pmatrix},vx​=2​∣false,false⟩+sx​∣true,true⟩​,K0​=(10​0η​​),K1​=(00​1−η​0​),

where each ket is the coordinate vector equal to one at its indicated pair and zero elsewhere. The assertion is that applying the map X↦K0XK0†+K1XK1†X\mapsto K_0XK_0^\dagger+K_1XK_1^\daggerX↦K0​XK0†​+K1​XK1†​ to the second Boolean factor of the matrix vxvx†v_xv_x^\daggervx​vx†​ gives, for every pair of Boolean pairs p,qp,qp,q,

[∑r=01(I2⊗Kr) vxvx† (I2⊗Kr)†]p,q=(1200sxη20000001−η20sxη200η2)c(p),c(q),\left[\sum_{r=0}^{1}(I_2\otimes K_r)\,v_xv_x^\dagger\,(I_2\otimes K_r)^\dagger\right]_{p,q} = \begin{pmatrix} \frac12&0&0&\frac{s_x\sqrt{\eta}}2\\ 0&0&0&0\\ 0&0&\frac{1-\eta}2&0\\ \frac{s_x\sqrt{\eta}}2&0&0&\frac{\eta}2 \end{pmatrix}_{c(p),c(q)},[r=0∑1​(I2​⊗Kr​)vx​vx†​(I2​⊗Kr​)†]p,q​=​21​002sx​η​​​0000​0021−η​0​2sx​η​​002η​​​c(p),c(q)​,

where c(false,false)=0c(\mathrm{false},\mathrm{false})=0c(false,false)=0, c(false,true)=1c(\mathrm{false},\mathrm{true})=1c(false,true)=1, c(true,false)=2c(\mathrm{true},\mathrm{false})=2c(true,false)=2, and c(true,true)=3c(\mathrm{true},\mathrm{true})=3c(true,true)=3, and †\dagger† denotes conjugate transpose.

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