Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Polar codes for degradable quantum channels

Definition
QAlgorithms_Papers_Independent_Intakethird100_P10

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

prove2quantum

The objects the formalized statements of Polar codes for degradable quantum channels (arXiv:1109.5346) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • BinaryCQ: PDF p.3: a pair of binary cq letter operators; validity is IsBinaryCQ.
  • IsBinaryCQ: PDF p.3: both letter operators are positive semidefinite and trace one.
  • symmetricInformation: PDF p.3: symmetric Holevo information, with the frozen entropy in bits.
  • polarCoefficient: PDF p.3: coefficient (j,k) of B_N F_2^(tensor n), over GF(2). The address map is the frozen most-significant-bit-first map; B_N reverses j's n-bit address before applying the tensor kernel.
  • polarWord: PDF p.3: row-vector multiplication u -> u G_N over GF(2).
  • combinedOutput: PDF p.3: N independent memoryless uses, after the prescribed row polar transform.
  • SplitRegister: PDF p.3, (1)–(3): prefix classical register and all N quantum outputs.
  • synthesized: PDF p.3, (1)–(3): the synthesized letter density, constructed by uniform averaging over every bit other than i and retaining the preceding bits. The denominator 2^(N-1) includes both prefix and suffix averaging.
  • synthesizedDegrader: PDF p.4, proof of Lemma 3: identity on prefix bits, D^(tensor N) on outputs.
  • polarThreshold: PDF pp.3–4: 2^(-N^beta), for N=2^n; beta is a real exponent.
  • goodSet: PDF p.3, (4): the strict ROOT-fidelity good set. Defined for every real beta; operational uses impose 0<beta<1/2, whereas Theorem 2 only assumes beta<1/2.
  • badSet: PDF p.3, (5): the complement of the strict good set; equality is bad.
  • poorSet: PDF p.4: the strict ROOT-fidelity poor set for Eve.
  • leafProbability: PDF pp.3–6: uniform probability of a set of n-step fair synthesis leaves. N=2^n is always positive, including n=0.
  • Wiretap: PDF pp.4–5: actual joint cq wiretap output. One CPTP map degrades both Bob letters to Eve letters; the joint output need not be a channel complement.
  • Wiretap.bob: PDF p.4, (8): Bob's marginal of the actual joint cq output.
  • Wiretap.eve: PDF p.4, (9): Eve's marginal of the same joint cq output.
  • informationSet: PDF p.5: information indices A=P_E intersect G_B; message count is 2^|A|.
  • frozenSet: PDF p.5: frozen indices B=P_E intersect B_B.
  • assistanceSet: PDF p.5: secret-key/ebit indices X=P_E^c intersect B_B.
  • randomizationSet: PDF p.5: independent randomization indices Y=P_E^c intersect G_B.
  • dilationLetter: PDF pp.6–7: receiver and environment of one fixed qubit-input isometry.
  • dilationBob: PDF pp.6–7: Bob's cq channel induced by the specified orthonormal input basis.
  • dilationEve: PDF pp.6–7: complementary cq channel of the SAME isometry and basis.
  • bellBobOutput: PDF pp.6–7: reference–Bob output from the normalized Bell state in this basis, with E traced out. The input to U is the second Bell half.
  • symmetricCoherentInformation: PDF pp.6–7: symmetric coherent information H(B)-H(AB) for the normalized Bell input.
  • PolarGate: PDF pp.3,7: the elementary gate alphabet used for the polar transform alone. Distinct control and target wires ensure that each CNOT is a reversible quantum gate.
  • PolarGate.apply: PDF pp.3,7, (15): the computational-basis action of a CNOT or SWAP gate.
  • circuitAction: PDF pp.3,7: circuit action in list order. Its elementary-gate count is list length.
  • circuitOperator: PDF p.7: coherent linear extension of a reversible classical gate network.
  • polarOperator: PDF p.7: exact coherent u -> u G_N, with no added phases.
  • pureCQ: PDF p.10, Appendix A: the pure-output cq channel, with validity supplied by unit norms.
  • minusChannel: PDF p.10, Appendix A: the worse channel's even and odd parity mixtures, constructed from the input channel, not independent matrices constrained by a target equality.
  • amplitudeComplement: PDF p.11: the printed amplitude-damping COMPLEMENT at original damping gamma. This is linear on all matrices: the upper-left entry is rho_00+(1-gamma)rho_11. It is damping with parameter 1-gamma, not with parameter gamma.
  • detectedJumpComplement: PDF p.11: the photon-detected-jump complement, extended linearly using trace rather than an affine constant 1 on arbitrary operators.
  • basisColumn: PDF p.11: column x of an arbitrary orthonormal basis, represented by a unitary matrix.
  • erasureData: PDF p.11: embedding into the data sector orthogonal to Option.none, the erasure flag.
  • dephasingComplement: PDF p.11: complement obtained by tracing B out of that isometry.
  • dephasingDiagonal: PDF p.11: both environment outputs when the basis is unbiased to the Pauli eigenbasis.
  • cloningNormalization: PDF p.11, (18)–(19): Delta_m=m(m+1)/2, with m a positive clone count in uses.
  • cloningEnvironment: PDF p.11, (18)–(19): the two displayed complementary cloning outputs. The common environmental orthonormal basis can change with the chosen input basis by the cited covariance. No full cloning channel is reconstructed from these two states.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_QuantumInformation
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_Ensemble
import Definitions.Def_QAlgorithms_Papers_Independent_Support50_P07
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P58
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P67
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P48

/-!
Wilde–Guha, Polar codes for degradable quantum channels, arXiv:1109.5346v3.
All page references are PDF pages. Local definitions are specific to this paper.
Finite nonempty matrix registers; entropy in bits; F is SQUARED fidelity.
Indices in Lean are zero based: i : Fin (2^n) denotes printed index i+1.
Uniform counting over all synthesized leaves is the n-step marginal of the fair
independent plus/minus birth process. No asymptotic operational capacity is
identified with an entropy formula by definition. Deferred nodes are in HARD.md.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P10

open scoped BigOperators Matrix ComplexConjugate
open Filter

-- Accepted definitions, without paper-local aliases.
open QAlgorithms.Papers.Independent.Information100.P58 (basisDensity)
open QAlgorithms.Papers.Independent.Information100.P67 (squaredFidelity)
open QAlgorithms.Papers.Independent.Claimed50.P48 (indexedKrausDilation)

/-- PDF p.3: a pair of binary cq letter operators; validity is `IsBinaryCQ`. -/
abbrev BinaryCQ (H : Type*) := Bool → Matrix H H ℂ

/-- PDF p.3: both letter operators are positive semidefinite and trace one. -/
def IsBinaryCQ {H : Type*} [Fintype H] (W : BinaryCQ H) : Prop :=
  ∀ x, IsDensityMatrix (W x)

/-- PDF p.3: symmetric Holevo information, with the frozen entropy in bits. -/
def symmetricInformation {H : Type*} [Fintype H] [DecidableEq H]
    (W : BinaryCQ H) : ℝ :=
  vnEntropy (((2 : ℂ)⁻¹) • (W false + W true)) -
    (vnEntropy (W false) + vnEntropy (W true)) / 2

/-- PDF p.3: coefficient (j,k) of B_N F_2^(tensor n), over GF(2).
The address map is the frozen most-significant-bit-first map; B_N reverses
j's n-bit address before applying the tensor kernel. -/
def polarCoefficient (n : ℕ) (j k : Fin (2 ^ n)) : Bool :=
  decide (∀ t : Fin n,
    natToBits n k.val t = true → natToBits n j.val (Fin.rev t) = true)

/-- PDF p.3: row-vector multiplication u -> u G_N over GF(2). -/
def polarWord (n : ℕ) (u : Fin (2 ^ n) → Bool) : Fin (2 ^ n) → Bool :=
  fun k => decide ((∑ j : Fin (2 ^ n),
    if u j && polarCoefficient n j k then (1 : ZMod 2) else 0) = 1)

/-- PDF p.3: N independent memoryless uses, after the prescribed row polar transform. -/
def combinedOutput {H : Type*} (W : BinaryCQ H) (n : ℕ)
    (u : Fin (2 ^ n) → Bool) : Matrix (Fin (2 ^ n) → H) (Fin (2 ^ n) → H) ℂ :=
  fun a b => ∏ k, W (polarWord n u k) (a k) (b k)

/-- PDF p.3, (1)–(3): prefix classical register and all N quantum outputs. -/
abbrev SplitRegister (H : Type*) (n : ℕ) (i : Fin (2 ^ n)) :=
  (Fin i.val → Bool) × (Fin (2 ^ n) → H)

/-- PDF p.3, (1)–(3): the synthesized letter density, constructed by uniform
averaging over every bit other than i and retaining the preceding bits.
The denominator 2^(N-1) includes both prefix and suffix averaging. -/
def synthesized {H : Type*} [Fintype H] [DecidableEq H]
    (W : BinaryCQ H) (n : ℕ) (i : Fin (2 ^ n)) : BinaryCQ (SplitRegister H n i) := by
  classical
  exact fun x a b => if a.1 = b.1 then
    ((2 : ℂ) ^ (2 ^ n - 1))⁻¹ *
      ∑ u : Fin (2 ^ n) → Bool,
        if u i = x ∧ ∀ j : Fin i.val,
          u ⟨j.val, Nat.lt_trans j.isLt i.isLt⟩ = a.1 j
        then combinedOutput W n u a.2 b.2 else 0
    else 0

/-- PDF p.4, proof of Lemma 3: identity on prefix bits, D^(tensor N) on outputs. -/
def synthesizedDegrader {H K : Type*} [Fintype H] [DecidableEq H]
    (D : Superop H K) (n : ℕ) (i : Fin (2 ^ n)) :
    Superop (SplitRegister H n i) (SplitRegister K n i) :=
  idTensorSuperop (R := Fin i.val → Bool) (piSuperop (k := 2 ^ n) (fun _ => D))

/-- PDF pp.3–4: 2^(-N^beta), for N=2^n; beta is a real exponent. -/
def polarThreshold (n : ℕ) (beta : ℝ) : ℝ :=
  (2 : ℝ) ^ (-(((2 : ℝ) ^ n) ^ beta))

/-- PDF p.3, (4): the strict ROOT-fidelity good set. Defined for every real beta;
operational uses impose 0<beta<1/2, whereas Theorem 2 only assumes beta<1/2. -/
def goodSet {H : Type*} [Fintype H] [DecidableEq H]
    (W : BinaryCQ H) (n : ℕ) (beta : ℝ) : Finset (Fin (2 ^ n)) := by
  classical
  exact Finset.univ.filter fun i =>
    rootFidelity (synthesized W n i false) (synthesized W n i true) < polarThreshold n beta

/-- PDF p.3, (5): the complement of the strict good set; equality is bad. -/
def badSet {H : Type*} [Fintype H] [DecidableEq H]
    (W : BinaryCQ H) (n : ℕ) (beta : ℝ) : Finset (Fin (2 ^ n)) := by
  classical
  exact Finset.univ \ goodSet W n beta

/-- PDF p.4: the strict ROOT-fidelity poor set for Eve. -/
def poorSet {H : Type*} [Fintype H] [DecidableEq H]
    (W : BinaryCQ H) (n : ℕ) (beta : ℝ) : Finset (Fin (2 ^ n)) := by
  classical
  exact Finset.univ.filter fun i =>
    1 - polarThreshold n beta <
      rootFidelity (synthesized W n i false) (synthesized W n i true)

/-- PDF pp.3–6: uniform probability of a set of n-step fair synthesis leaves.
N=2^n is always positive, including n=0. -/
def leafProbability {n : ℕ} (S : Finset (Fin (2 ^ n))) : ℝ :=
  (S.card : ℝ) / ((2 : ℝ) ^ n)

/-- PDF pp.4–5: actual joint cq wiretap output. One CPTP map degrades both
Bob letters to Eve letters; the joint output need not be a channel complement. -/
structure Wiretap (B E : Type*) [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] where
  output : BinaryCQ (B × E)
  density : IsBinaryCQ output
  degrading : Superop B E
  cptp : IsChannel degrading
  degrades : ∀ x, traceLeft (output x) = degrading (traceRight (output x))
  classicalEve : Commute (traceLeft (output false)) (traceLeft (output true))

/-- PDF p.4, (8): Bob's marginal of the actual joint cq output. -/
def Wiretap.bob {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] (V : Wiretap B E) : BinaryCQ B :=
  fun x => traceRight (V.output x)

/-- PDF p.4, (9): Eve's marginal of the same joint cq output. -/
def Wiretap.eve {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] (V : Wiretap B E) : BinaryCQ E :=
  fun x => traceLeft (V.output x)

/-- PDF p.5: information indices A=P_E intersect G_B; message count is 2^|A|. -/
def informationSet {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] (V : Wiretap B E) (n : ℕ) (beta : ℝ) :=
  poorSet V.eve n beta ∩ goodSet V.bob n beta

/-- PDF p.5: frozen indices B=P_E intersect B_B. -/
def frozenSet {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] (V : Wiretap B E) (n : ℕ) (beta : ℝ) :=
  poorSet V.eve n beta ∩ badSet V.bob n beta

/-- PDF p.5: secret-key/ebit indices X=P_E^c intersect B_B. -/
def assistanceSet {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] (V : Wiretap B E) (n : ℕ) (beta : ℝ) := by
  classical
  exact (Finset.univ \ poorSet V.eve n beta) ∩ badSet V.bob n beta

/-- PDF p.5: independent randomization indices Y=P_E^c intersect G_B. -/
def randomizationSet {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] [DecidableEq E] (V : Wiretap B E) (n : ℕ) (beta : ℝ) := by
  classical
  exact (Finset.univ \ poorSet V.eve n beta) ∩ goodSet V.bob n beta

/-- PDF pp.6–7: receiver and environment of one fixed qubit-input isometry. -/
def dilationLetter {B E : Type*} [Fintype B] [Fintype E]
    (U : Matrix (B × E) Bool ℂ) (basis : Matrix.unitaryGroup Bool ℂ) (x : Bool) :
    Matrix (B × E) (B × E) ℂ :=
  let v : EuclideanSpace ℂ (B × E) := WithLp.toLp 2 fun z =>
    ∑ a : Bool, U z a * (basis : Matrix Bool Bool ℂ) a x
  pureDensity v

/-- PDF pp.6–7: Bob's cq channel induced by the specified orthonormal input basis. -/
def dilationBob {B E : Type*} [Fintype B] [Fintype E]
    (U : Matrix (B × E) Bool ℂ) (basis : Matrix.unitaryGroup Bool ℂ) : BinaryCQ B :=
  fun x => traceRight (dilationLetter U basis x)

/-- PDF pp.6–7: complementary cq channel of the SAME isometry and basis. -/
def dilationEve {B E : Type*} [Fintype B] [Fintype E]
    (U : Matrix (B × E) Bool ℂ) (basis : Matrix.unitaryGroup Bool ℂ) : BinaryCQ E :=
  fun x => traceLeft (dilationLetter U basis x)

/-- PDF pp.6–7: reference–Bob output from the normalized Bell state in this basis,
with E traced out. The input to U is the second Bell half. -/
def bellBobOutput {B E : Type*} [Fintype B] [Fintype E]
    (U : Matrix (B × E) Bool ℂ) (basis : Matrix.unitaryGroup Bool ℂ) :
    Matrix (Bool × B) (Bool × B) ℂ :=
  fun r s => (2 : ℂ)⁻¹ * ∑ e : E,
    (∑ a : Bool, U (r.2, e) a * (basis : Matrix Bool Bool ℂ) a r.1) *
      conj (∑ a : Bool, U (s.2, e) a * (basis : Matrix Bool Bool ℂ) a s.1)

/-- PDF pp.6–7: symmetric coherent information H(B)-H(AB) for the normalized Bell input. -/
def symmetricCoherentInformation {B E : Type*} [Fintype B] [DecidableEq B]
    [Fintype E] (U : Matrix (B × E) Bool ℂ) (basis : Matrix.unitaryGroup Bool ℂ) : ℝ :=
  vnEntropy (traceLeft (bellBobOutput U basis)) - vnEntropy (bellBobOutput U basis)

/-- n=0 has a single polar coefficient, equal to one. -/
example : polarCoefficient 0 0 0 = true := by simp [polarCoefficient]

/-- The uniform probability of no leaves is zero, including n=0. -/
example (n : ℕ) : leafProbability (n := n) ∅ = 0 := by simp [leafProbability]

/-- Averages and independent-use products retain the letter operators at blocklength one. -/
example {H : Type*} (W : BinaryCQ H) (u : Fin 1 → Bool) (a b : Fin 1 → H) :
    combinedOutput W 0 u a b = W (polarWord 0 u 0) (a 0) (b 0) := by
  simp [combinedOutput]

/-- PDF pp.3,7: the elementary gate alphabet used for the polar transform alone.
Distinct control and target wires ensure that each CNOT is a reversible quantum gate. -/
inductive PolarGate (N : ℕ) where
  | cnot (control target : Fin N) (distinct : control ≠ target)
  | swap (first second : Fin N)

/-- PDF pp.3,7, (15): the computational-basis action of a CNOT or SWAP gate. -/
def PolarGate.apply {N : ℕ} (g : PolarGate N) (u : Fin N → Bool) : Fin N → Bool :=
  match g with
  | .cnot c t _ => Function.update u t (xor (u c) (u t))
  | .swap a b => Function.update (Function.update u a (u b)) b (u a)

/-- PDF pp.3,7: circuit action in list order. Its elementary-gate count is list length. -/
def circuitAction {N : ℕ} (C : List (PolarGate N)) (u : Fin N → Bool) : Fin N → Bool :=
  C.foldl (fun v g => g.apply v) u

/-- PDF p.7: coherent linear extension of a reversible classical gate network. -/
def circuitOperator {N : ℕ} (C : List (PolarGate N)) : Matrix (Qubits N) (Qubits N) ℂ :=
  fun v u => if v = circuitAction C u then 1 else 0

/-- PDF p.7: exact coherent u -> u G_N, with no added phases. -/
def polarOperator (n : ℕ) : Matrix (Qubits (2 ^ n)) (Qubits (2 ^ n)) ℂ :=
  fun v u => if v = polarWord n u then 1 else 0

/-- PDF p.10, Appendix A: the pure-output cq channel, with validity supplied by unit norms. -/
def pureCQ {H : Type*} (psi0 psi1 : EuclideanSpace ℂ H) : BinaryCQ H :=
  fun x => pureDensity (if x then psi1 else psi0)

/-- PDF p.10, Appendix A: the worse channel's even and odd parity mixtures,
constructed from the input channel, not independent matrices constrained by a target equality. -/
def minusChannel {H : Type*} (W : BinaryCQ H) : BinaryCQ (H × H) :=
  fun x => if x then ((2 : ℂ)⁻¹) •
    (Matrix.kronecker (W false) (W true) + Matrix.kronecker (W true) (W false))
  else ((2 : ℂ)⁻¹) •
    (Matrix.kronecker (W false) (W false) + Matrix.kronecker (W true) (W true))

/-- PDF p.11: the printed amplitude-damping COMPLEMENT at original damping gamma.
This is linear on all matrices: the upper-left entry is rho_00+(1-gamma)rho_11.
It is damping with parameter 1-gamma, not with parameter gamma. -/
def amplitudeComplement (gamma : ℝ) : Superop Bool Bool :=
  fun rho a b => if a then
    if b then (gamma : ℂ) * rho true true else (Real.sqrt gamma : ℂ) * rho true false
  else if b then (Real.sqrt gamma : ℂ) * rho false true
  else rho false false + ((1 - gamma : ℝ) : ℂ) * rho true true

/-- PDF p.11: the photon-detected-jump complement, extended linearly using trace
rather than an affine constant 1 on arbitrary operators. -/
def detectedJumpComplement (gamma : ℝ) : Superop Bool Bool :=
  fun rho => Matrix.diagonal fun x => if x then (gamma : ℂ) * rho true true
    else rho.trace - (gamma : ℂ) * rho true true

/-- PDF p.11: column x of an arbitrary orthonormal basis, represented by a unitary matrix. -/
def basisColumn {H : Type*} [Fintype H] [DecidableEq H]
    (basis : Matrix.unitaryGroup H ℂ) (x : H) : EuclideanSpace ℂ H :=
  WithLp.toLp 2 fun a => (basis : Matrix H H ℂ) a x

/-- PDF p.11: embedding into the data sector orthogonal to Option.none, the erasure flag. -/
def erasureData (d : ℕ) (rho : Matrix (Fin d) (Fin d) ℂ) :
    Matrix (Option (Fin d)) (Option (Fin d)) ℂ :=
  fun i j => match i, j with
    | some a, some b => rho a b
    | _, _ => 0

/-- PDF p.11: complement obtained by tracing B out of that isometry. -/
def dephasingComplement (p : ℝ) (sigma : Matrix Bool Bool ℂ) : Superop Bool Bool :=
  let V := indexedKrausDilation (fun x : Bool =>
    if x then (Real.sqrt p : ℂ) • sigma
    else (Real.sqrt (1 - p) : ℂ) • (1 : Matrix Bool Bool ℂ))
  fun rho => traceLeft (V * rho * Vᴴ)

/-- PDF p.11: both environment outputs when the basis is unbiased to the Pauli eigenbasis. -/
def dephasingDiagonal (p : ℝ) : Matrix Bool Bool ℂ :=
  Matrix.diagonal fun x => if x then (p : ℂ) else ((1 - p : ℝ) : ℂ)

/-- PDF p.11, (18)–(19): Delta_m=m(m+1)/2, with m a positive clone count in uses. -/
def cloningNormalization (m : ℕ) : ℝ := (m : ℝ) * (m + 1) / 2

/-- PDF p.11, (18)–(19): the two displayed complementary cloning outputs.
The common environmental orthonormal basis can change with the chosen input basis
by the cited covariance. No full cloning channel is reconstructed from these two states. -/
def cloningEnvironment (m : ℕ) (basis : Matrix.unitaryGroup (Fin m) ℂ) (x : Bool) :
    Matrix (Fin m) (Fin m) ℂ :=
  ∑ i : Fin m, (((i.val + 1 : ℝ) / cloningNormalization m : ℝ) : ℂ) •
    pureDensity (basisColumn basis (if x then Fin.rev i else i))

/-- An empty circuit acts as the identity, including the n=0 transform case. -/
example {N : ℕ} (u : Fin N → Bool) : circuitAction [] u = u := rfl

/-- The worse branch is the even-parity mixture, with weight 1/2 on each term. -/
example {H : Type*} (W : BinaryCQ H) : minusChannel W false = ((2 : ℂ)⁻¹) •
    (Matrix.kronecker (W false) (W false) + Matrix.kronecker (W true) (W true)) := rfl

/-- The complementary damping off-diagonal coefficient is sqrt(gamma). -/
example (gamma : ℝ) (rho : Matrix Bool Bool ℂ) :
    amplitudeComplement gamma rho true false = (Real.sqrt gamma : ℂ) * rho true false := rfl

/-- The detected-jump complement is linear at zero, with no affine constant one. -/
example (gamma : ℝ) : detectedJumpComplement gamma 0 = 0 := by
  ext a b
  cases a <;> cases b <;> simp [detectedJumpComplement]

/-- The erasure data embedding has no component in the orthogonal flag sector. -/
example (d : ℕ) (rho : Matrix (Fin d) (Fin d) ℂ) : erasureData d rho none none = 0 := rfl

/-- The no-error branch of the Pauli isometry is sqrt(1-p) times the identity. -/
example (p : ℝ) (sigma : Matrix Bool Bool ℂ) :
    indexedKrausDilation (fun x : Bool =>
      if x then (Real.sqrt p : ℂ) • sigma
      else (Real.sqrt (1 - p) : ℂ) • (1 : Matrix Bool Bool ℂ))
        (false, false) true = 0 := by simp [indexedKrausDilation]

/-- At one clone the normalizing constant is one, not zero. -/
example : cloningNormalization 1 = 1 := by norm_num [cloningNormalization]

end QAlgorithms.Papers.Independent.Intakethird100.P10
Source
Polar codes for degradable quantum channels, arXiv:1109.5346v3, https://arxiv.org/abs/1109.5346v3
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P10.BinaryCQ

For any type HHH, BinaryCQ(H)\mathrm{BinaryCQ}(H)BinaryCQ(H) is the type of functions assigning a complex matrix indexed by H×HH\times HH×H to each Boolean value. Thus an element is an arbitrary pair (W0,W1)(W_0,W_1)(W0​,W1​) of such matrices, with 000 denoting false and 111 denoting true; this abbreviation imposes no finiteness, positivity, or normalization requirements.

QAlgorithms.Papers.Independent.Intakethird100.P10.IsBinaryCQ

For any finite type HHH and any pair W=(W0,W1)W=(W_0,W_1)W=(W0​,W1​) of complex H×HH\times HH×H matrices, IsBinaryCQ(W)\mathrm{IsBinaryCQ}(W)IsBinaryCQ(W) means that both W0W_0W0​ and W1W_1W1​ are positive semidefinite and have complex trace exactly 111. If HHH is empty, this condition cannot hold, since every matrix on HHH has trace 000.

QAlgorithms.Papers.Independent.Intakethird100.P10.PolarGate.apply

For every natural number NNN, every gate description ggg consisting either of a controlled-NOT with distinct indices c,t<Nc,t<Nc,t<N or a swap with arbitrary indices a,b<Na,b<Na,b<N, and every Boolean word uuu of length NNN, this declaration defines the resulting word as follows. A controlled-NOT replaces utu_tut​ by uc⊕utu_c\mathbin{\oplus}u_tuc​⊕ut​ and leaves all other coordinates unchanged, where ⊕\oplus⊕ is Boolean exclusive OR. A swap first replaces coordinate aaa by the original ubu_bub​, then coordinate bbb by the original uau_aua​, leaving other coordinates unchanged. When a=ba=ba=b, this action is the identity. The parameter NNN may be zero, although then no gate description exists.

QAlgorithms.Papers.Independent.Intakethird100.P10.PolarGate

For each natural number NNN, this declaration introduces a type of gate descriptions with exactly two forms: a controlled-NOT description carrying a control index ccc and target index ttt in {0,…,N−1}\{0,\ldots,N-1\}{0,…,N−1}, together with a proof that c≠tc\ne tc=t; and a swap description carrying two indices a,ba,ba,b in the same set, with no requirement that they differ. This declaration specifies the descriptions and their data, without assigning an action to them. For N=0N=0N=0 there are no such descriptions; for N=1N=1N=1, only the swap with both indices 000 can be formed.

QAlgorithms.Papers.Independent.Intakethird100.P10.SplitRegister

For any type HHH, any natural number nnn, and any index i∈{0,…,2n−1}i\in\{0,\ldots,2^n-1\}i∈{0,…,2n−1}, the defined register type is

{0,1} i×H 2n.\{0,1\}^{\,i}\times H^{\,2^n}.{0,1}i×H2n.

Its elements consist of a Boolean word of length iii and an HHH-valued word of length 2n2^n2n. At i=0i=0i=0, the Boolean word is the unique empty word; n=0n=0n=0 forces i=0i=0i=0 and leaves one HHH-valued coordinate. No finiteness or nonemptiness assumption on HHH is required.

QAlgorithms.Papers.Independent.Intakethird100.P10.Wiretap.bob

For finite types B,EB,EB,E with decidable equality and any wiretap datum VVV, this defines the pair of BBB-indexed matrices

WxB(b,b′)=∑e∈Eρx((b,e),(b′,e)),x∈{0,1},W_x^B(b,b')=\sum_{e\in E}\rho_x((b,e),(b',e)),\qquad x\in\{0,1\},WxB​(b,b′)=e∈E∑​ρx​((b,e),(b′,e)),x∈{0,1},

where ρx=V.output(x)\rho_x=V.\mathrm{output}(x)ρx​=V.output(x). Here the datum VVV includes the requirements that both ρx\rho_xρx​ are positive semidefinite of trace 111, that there is a map D(X)=∑j<rKjXKj†D(X)=\sum_{j<r}K_jXK_j^\daggerD(X)=∑j<r​Kj​XKj†​ for some finite family satisfying ∑jKj†Kj=IB\sum_jK_j^\dagger K_j=I_B∑j​Kj†​Kj​=IB​, that its application to each displayed BBB marginal equals the corresponding EEE marginal ρxE(e,e′)=∑bρx((b,e),(b,e′))\rho_x^E(e,e')=\sum_b\rho_x((b,e),(b,e'))ρxE​(e,e′)=∑b​ρx​((b,e),(b,e′)), and that the two EEE marginals commute. The definition itself takes the partial trace over EEE for each Boolean input.

QAlgorithms.Papers.Independent.Intakethird100.P10.Wiretap.eve

For finite types B,EB,EB,E with decidable equality and any wiretap datum VVV, this defines the pair of EEE-indexed matrices

WxE(e,e′)=∑b∈Bρx((b,e),(b,e′)),x∈{0,1},W_x^E(e,e')=\sum_{b\in B}\rho_x((b,e),(b,e')),\qquad x\in\{0,1\},WxE​(e,e′)=b∈B∑​ρx​((b,e),(b,e′)),x∈{0,1},

where ρx=V.output(x)\rho_x=V.\mathrm{output}(x)ρx​=V.output(x). The datum VVV requires both ρx\rho_xρx​ to be positive semidefinite of trace 111; it supplies a map D(X)=∑j<rKjXKj†D(X)=\sum_{j<r}K_jXK_j^\daggerD(X)=∑j<r​Kj​XKj†​ with ∑jKj†Kj=IB\sum_jK_j^\dagger K_j=I_B∑j​Kj†​Kj​=IB​, requires WxE=D(WxB)W_x^E=D(W_x^B)WxE​=D(WxB​), where WxB(b,b′)=∑eρx((b,e),(b′,e))W_x^B(b,b')=\sum_e\rho_x((b,e),(b',e))WxB​(b,b′)=∑e​ρx​((b,e),(b′,e)), and requires W0EW1E=W1EW0EW_0^EW_1^E=W_1^EW_0^EW0E​W1E​=W1E​W0E​. The definition itself takes the partial trace over BBB for each Boolean input.

QAlgorithms.Papers.Independent.Intakethird100.P10.Wiretap

For finite types B,EB,EB,E with decidable equality, a value of this structure consists of two complex matrices ρ0,ρ1\rho_0,\rho_1ρ0​,ρ1​ indexed by B×EB\times EB×E, each positive semidefinite with trace 111, together with a function DDD from complex BBB-indexed square matrices to complex EEE-indexed square matrices that admits a finite Kraus representation: there exist a natural number rrr and matrices KjK_jKj​ with rows indexed by EEE and columns by BBB such that ∑j<rKj†Kj=IB\sum_{j<r}K_j^\dagger K_j=I_B∑j<r​Kj†​Kj​=IB​ and D(X)=∑j<rKjXKj†D(X)=\sum_{j<r}K_jXK_j^\daggerD(X)=∑j<r​Kj​XKj†​ for every matrix XXX. For each x∈{0,1}x\in\{0,1\}x∈{0,1}, the structure requires

ρxE=D(ρxB),ρxB(b,b′)=∑e∈Eρx((b,e),(b′,e)),ρxE(e,e′)=∑b∈Bρx((b,e),(b,e′)).\rho_x^E=D(\rho_x^B),\qquad \rho_x^B(b,b')=\sum_{e\in E}\rho_x((b,e),(b',e)),\qquad \rho_x^E(e,e')=\sum_{b\in B}\rho_x((b,e),(b,e')).ρxE​=D(ρxB​),ρxB​(b,b′)=e∈E∑​ρx​((b,e),(b′,e)),ρxE​(e,e′)=b∈B∑​ρx​((b,e),(b,e′)).

It also requires ρ0Eρ1E=ρ1Eρ0E\rho_0^E\rho_1^E=\rho_1^E\rho_0^Eρ0E​ρ1E​=ρ1E​ρ0E​. No nonemptiness assumptions on B,EB,EB,E are explicitly supplied, but the required trace-one joint matrices cannot exist if either type is empty.

QAlgorithms.Papers.Independent.Intakethird100.P10.amplitudeComplement

For every real number γ\gammaγ, define a function on all complex 2×22\times22×2 matrices, with rows and columns ordered as false, true, by

Aγ(ρ)=(ρ00+(1−γ)ρ11γ ρ01γ ρ10γρ11).\mathcal A_\gamma(\rho)= \begin{pmatrix} \rho_{00}+(1-\gamma)\rho_{11}&\sqrt{\gamma}\,\rho_{01}\\ \sqrt{\gamma}\,\rho_{10}&\gamma\rho_{11} \end{pmatrix}.Aγ​(ρ)=(ρ00​+(1−γ)ρ11​γ​ρ10​​γ​ρ01​γρ11​​).

No restriction on γ\gammaγ or on the input matrix is imposed. The real square root is total and equals 000 for negative γ\gammaγ.

QAlgorithms.Papers.Independent.Intakethird100.P10.assistanceSet

For finite types B,EB,EB,E with decidable equality, a wiretap datum VVV, any natural number nnn, and any real β\betaβ, this defines the set of indices i<2ni<2^ni<2n satisfying

f(Wn,i,0E,Wn,i,1E)≤1−τandτ≤f(Wn,i,0B,Wn,i,1B),τ=2−(2n)β.f(W_{n,i,0}^E,W_{n,i,1}^E)\le1-\tau \quad\text{and}\quad \tau\le f(W_{n,i,0}^B,W_{n,i,1}^B), \qquad \tau=2^{-(2^n)^\beta}.f(Wn,i,0E​,Wn,i,1E​)≤1−τandτ≤f(Wn,i,0B​,Wn,i,1B​),τ=2−(2n)β.

Here VVV supplies joint positive semidefinite trace-one matrices ρx\rho_xρx​, the marginals WxB(b,b′)=∑eρx((b,e),(b′,e))W_x^B(b,b')=\sum_e\rho_x((b,e),(b',e))WxB​(b,b′)=∑e​ρx​((b,e),(b′,e)) and WxE(e,e′)=∑bρx((b,e),(b,e′))W_x^E(e,e')=\sum_b\rho_x((b,e),(b,e'))WxE​(e,e′)=∑b​ρx​((b,e),(b,e′)), and a finite Kraus channel D(X)=∑jKjXKj†D(X)=\sum_jK_jXK_j^\daggerD(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=IB\sum_jK_j^\dagger K_j=I_B∑j​Kj†​Kj​=IB​ and WxE=D(WxB)W_x^E=D(W_x^B)WxE​=D(WxB​); its two EEE marginals commute. For either pair WWW, put N=2nN=2^nN=2n and define the synthesized matrices by

Wn,i,x((v,h),(v′,h′))={2−(N−1)∑u∈{0,1}Nui=x, u<i=v∏k<NWp(u)k(hk,hk′),v=v′,0,v≠v′,W_{n,i,x}((v,h),(v',h'))= \begin{cases} 2^{-(N-1)}\displaystyle\sum_{\substack{u\in\{0,1\}^N\\u_i=x,\ u_{<i}=v}} \prod_{k<N}W_{p(u)_k}(h_k,h'_k),&v=v',\\ 0,&v\ne v', \end{cases}Wn,i,x​((v,h),(v′,h′))=⎩⎨⎧​2−(N−1)u∈{0,1}Nui​=x, u<i​=v​∑​k<N∏​Wp(u)k​​(hk​,hk′​),0,​v=v′,v=v′,​

where vvv has length iii, hhh has length NNN over the marginal's index type, and p(u)k=∑jajkuj(mod2)p(u)_k=\sum_j a_{jk}u_j\pmod2p(u)k​=∑j​ajk​uj​(mod2), with ajk=1a_{jk}=1ajk​=1 precisely when every 111-bit of kkk appears at a reversed 111-bit position of jjj, using length-nnn binary representations. The function f(X,Y)f(X,Y)f(X,Y) is Re⁡tr⁡X YX\operatorname{Re}\operatorname{tr}\sqrt{\sqrt X\,Y\sqrt X}RetrX​YX​​. Both threshold equalities are included, and n=0n=0n=0 is included.

QAlgorithms.Papers.Independent.Intakethird100.P10.badSet

For any finite type HHH with decidable equality, any pair W=(W0,W1)W=(W_0,W_1)W=(W0​,W1​) of complex H×HH\times HH×H matrices, any natural number nnn, and any real β\betaβ, put N=2nN=2^nN=2n. The defined finite set consists exactly of indices i∈{0,…,N−1}i\in\{0,\ldots,N-1\}i∈{0,…,N−1} satisfying

2−Nβ≤Re⁡tr⁡Si,0 Si,1Si,0.2^{-N^\beta} \le \operatorname{Re}\operatorname{tr} \sqrt{\sqrt{S_{i,0}}\,S_{i,1}\sqrt{S_{i,0}}}.2−Nβ≤RetrSi,0​​Si,1​Si,0​​​.

All matrix square roots mean the total real continuous functional calculus of t↦tt\mapsto\sqrt tt↦t​, which returns zero on non-Hermitian matrices. The matrices Si,xS_{i,x}Si,x​, indexed by Boolean prefixes p,qp,qp,q of length iii and words a,b∈HNa,b\in H^Na,b∈HN, are zero when p≠qp\ne qp=q, and otherwise have entries

Si,x((p,a),(p,b))=2−(N−1)∑u∈{0,1}Nui=xuj=pj (j<i)∏k=0N−1(Wvk)ak,bk,vk=∑j=0N−1ujcn(j,k)(mod2),S_{i,x}((p,a),(p,b)) = 2^{-(N-1)} \sum_{\substack{u\in\{0,1\}^{N}\\u_i=x\\u_j=p_j\ (j<i)}} \prod_{k=0}^{N-1}(W_{v_k})_{a_k,b_k}, \qquad v_k=\sum_{j=0}^{N-1}u_jc_n(j,k)\pmod2,Si,x​((p,a),(p,b))=2−(N−1)u∈{0,1}Nui​=xuj​=pj​ (j<i)​∑​k=0∏N−1​(Wvk​​)ak​,bk​​,vk​=j=0∑N−1​uj​cn​(j,k)(mod2),

where cn(j,k)=1c_n(j,k)=1cn​(j,k)=1 exactly when every nonzero bit ttt of kkk's length-nnn, most-significant-first binary representation corresponds to a nonzero bit n−1−tn-1-tn−1−t of jjj's representation. Equality with the threshold is included. No density-matrix assumption on WWW or restriction on β\betaβ is imposed. At n=0n=0n=0, the sole index is included exactly when the root-fidelity expression for W0,W1W_0,W_1W0​,W1​ is at least 1/21/21/2; if HHH is empty, the set is empty.

QAlgorithms.Papers.Independent.Intakethird100.P10.basisColumn

For every finite index type HHH with decidable equality, every unitary complex matrix UUU indexed by HHH, and every x∈Hx\in Hx∈H, define the vector vx∈CHv_x\in\mathbb C^Hvx​∈CH by vx(a)=Ua,xv_x(a)=U_{a,x}vx​(a)=Ua,x​ for every a∈Ha\in Ha∈H; thus vxv_xvx​ is column xxx of UUU, viewed with the Euclidean norm. No separate nonemptiness assumption is imposed, although specifying xxx requires an element of HHH.

QAlgorithms.Papers.Independent.Intakethird100.P10.bellBobOutput

For arbitrary finite sets B,EB,EB,E, any complex matrix UUU with rows indexed by B×EB\times EB×E and Boolean columns, and any unitary 2×22\times22×2 matrix QQQ, this declaration defines a matrix Ω\OmegaΩ indexed by {0,1}×B\{0,1\}\times B{0,1}×B, with entries

Ω((x,b),(y,b′))=12∑e∈E(∑a∈{0,1}U((b,e),a)Q(a,x))(∑a∈{0,1}U((b′,e),a)Q(a,y))‾.\Omega((x,b),(y,b')) =\frac12\sum_{e\in E} \left(\sum_{a\in\{0,1\}}U((b,e),a)Q(a,x)\right) \overline{\left(\sum_{a\in\{0,1\}}U((b',e),a)Q(a,y)\right)}.Ω((x,b),(y,b′))=21​e∈E∑​​a∈{0,1}∑​U((b,e),a)Q(a,x)​​a∈{0,1}∑​U((b′,e),a)Q(a,y)​​.

No isometry assumption on UUU or nonemptiness assumption on B,EB,EB,E is made. An empty EEE makes all entries zero.

QAlgorithms.Papers.Independent.Intakethird100.P10.circuitAction

For every natural number NNN, every finite list C=[g1,…,gℓ]C=[g_1,\ldots,g_\ell]C=[g1​,…,gℓ​] of gate descriptions on indices 0,…,N−10,\ldots,N-10,…,N−1, and every Boolean word uuu of length NNN, this declaration defines the circuit action by starting with uuu and applying the listed gates in order, from g1g_1g1​ through gℓg_\ellgℓ​. Each controlled-NOT description has distinct control ccc and target ttt, and replaces the current target bit by the exclusive OR of the current control and target bits; each swap exchanges its two current bits, allowing equal indices, in which case it leaves the word unchanged. The empty list acts as the identity. This includes N=0N=0N=0, where the list must be empty and the unique empty word is unchanged.

QAlgorithms.Papers.Independent.Intakethird100.P10.circuitOperator

For every natural number NNN, including N=0N=0N=0, and every finite list CCC of gates on NNN Boolean wires, define a complex matrix indexed by bit strings u,v∈{0,1}Nu,v\in\{0,1\}^{N}u,v∈{0,1}N by

(Cop)v,u={1,v=C(u),0,otherwise.(C_{\mathrm{op}})_{v,u}= \begin{cases} 1,&v=C(u),\\ 0,&\text{otherwise}. \end{cases}(Cop​)v,u​={1,0,​v=C(u),otherwise.​

Here C(u)C(u)C(u) is obtained by applying the listed gates in order: a controlled-NOT has distinct control and target wires and replaces the target bit by its XOR with the control bit; a swap exchanges its two specified bits, which may be the same bit. An empty list acts as the identity.

QAlgorithms.Papers.Independent.Intakethird100.P10.cloningEnvironment

For every natural number mmm, every unitary complex m×mm\times mm×m matrix UUU, and every Boolean xxx, let uju_juj​ denote column jjj of UUU. This defines

ρx=∑i=0m−1i+1m(m+1)/2 urx(i)urx(i)†,rfalse(i)=i,rtrue(i)=m−1−i.\rho_x=\sum_{i=0}^{m-1}\frac{i+1}{m(m+1)/2}\, u_{r_x(i)}u_{r_x(i)}^\dagger, \qquad r_{\mathrm{false}}(i)=i,\quad r_{\mathrm{true}}(i)=m-1-i.ρx​=i=0∑m−1​m(m+1)/2i+1​urx​(i)​urx​(i)†​,rfalse​(i)=i,rtrue​(i)=m−1−i.

Each outer product has entries urx(i)(a)urx(i)(b)‾u_{r_x(i)}(a)\overline{u_{r_x(i)}(b)}urx​(i)​(a)urx​(i)​(b)​. The definition includes m=0m=0m=0, when the sum is empty and the resulting matrix is the zero matrix on the empty index set; real division by zero is defined to be zero.

QAlgorithms.Papers.Independent.Intakethird100.P10.cloningNormalization

For every natural number mmm, this defines the real number

cm=m(m+1)2.c_m=\frac{m(m+1)}2.cm​=2m(m+1)​.

In particular, c0=0c_0=0c0​=0.

QAlgorithms.Papers.Independent.Intakethird100.P10.combinedOutput

For any type HHH, any pair W=(W0,W1)W=(W_0,W_1)W=(W0​,W1​) of complex H×HH\times HH×H matrices, any natural number nnn, and any Boolean word uuu of length N=2nN=2^nN=2n, define a complex matrix indexed by functions a,b:{0,…,N−1}→Ha,b:\{0,\ldots,N-1\}\to Ha,b:{0,…,N−1}→H by

CW,n,u(a,b)=∏k=0N−1(Wvk)ak,bk.C_{W,n,u}(a,b)=\prod_{k=0}^{N-1}(W_{v_k})_{a_k,b_k}.CW,n,u​(a,b)=k=0∏N−1​(Wvk​​)ak​,bk​​.

Here vk=∑j=0N−1ujcn(j,k)(mod2)v_k=\sum_{j=0}^{N-1}u_jc_n(j,k)\pmod 2vk​=∑j=0N−1​uj​cn​(j,k)(mod2), where cn(j,k)=1c_n(j,k)=1cn​(j,k)=1 precisely when every nonzero bit ttt of kkk's length-nnn, most-significant-first binary representation corresponds to a nonzero bit n−1−tn-1-tn−1−t of jjj's representation. No finiteness assumption on HHH, and no positivity or normalization assumption on WWW, is imposed. When n=0n=0n=0, this matrix has entries (Wu0)a0,b0(W_{u_0})_{a_0,b_0}(Wu0​​)a0​,b0​​.

QAlgorithms.Papers.Independent.Intakethird100.P10.dephasingComplement

For every real ppp and every complex matrix σ\sigmaσ indexed by the two Boolean values, this defines a map Cp,σ\mathcal C_{p,\sigma}Cp,σ​ on all complex matrices ρ\rhoρ with those indices as follows: set K0=1−p IK_0=\sqrt{1-p}\,IK0​=1−p​I and K1=p σK_1=\sqrt p\,\sigmaK1​=p​σ, with 000 denoting false and 111 denoting true; then

(Cp,σ(ρ))b,c=tr⁡(KbρKc†)(b,c∈{0,1}).\bigl(\mathcal C_{p,\sigma}(\rho)\bigr)_{b,c} =\operatorname{tr}(K_b\rho K_c^\dagger) \qquad(b,c\in\{0,1\}).(Cp,σ​(ρ))b,c​=tr(Kb​ρKc†​)(b,c∈{0,1}).

No restriction on ppp, σ\sigmaσ, or ρ\rhoρ is imposed. The square roots are real square roots, which equal zero for negative arguments.

QAlgorithms.Papers.Independent.Intakethird100.P10.dephasingDiagonal

For every real ppp, this defines the complex matrix indexed in the order false, true by

Dp=(1−p00p).D_p=\begin{pmatrix}1-p&0\\0&p\end{pmatrix}.Dp​=(1−p0​0p​).

There is no assumption that ppp lies between zero and one.

QAlgorithms.Papers.Independent.Intakethird100.P10.detectedJumpComplement

For every real number γ\gammaγ, define a function on all complex 2×22\times22×2 matrices, with rows and columns ordered as false, true, by

Jγ(ρ)=(tr⁡ρ−γρ1100γρ11).\mathcal J_\gamma(\rho)= \begin{pmatrix} \operatorname{tr}\rho-\gamma\rho_{11}&0\\ 0&\gamma\rho_{11} \end{pmatrix}.Jγ​(ρ)=(trρ−γρ11​0​0γρ11​​).

No restriction on γ\gammaγ and no positivity or normalization assumption on ρ\rhoρ is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P10.dilationBob

For arbitrary finite sets B,EB,EB,E, any complex matrix UUU with rows indexed by B×EB\times EB×E and Boolean columns, and any unitary 2×22\times22×2 matrix QQQ, this declaration defines a Boolean-indexed family of matrices on BBB. For x∈{0,1}x\in\{0,1\}x∈{0,1} and b,b′∈Bb,b'\in Bb,b′∈B, its entries are

dilationBob⁡(U,Q)x(b,b′)=∑e∈E(∑a∈{0,1}U((b,e),a)Q(a,x))(∑a∈{0,1}U((b′,e),a)Q(a,x))‾.\operatorname{dilationBob}(U,Q)_x(b,b') =\sum_{e\in E} \left(\sum_{a\in\{0,1\}}U((b,e),a)Q(a,x)\right) \overline{\left(\sum_{a\in\{0,1\}}U((b',e),a)Q(a,x)\right)}.dilationBob(U,Q)x​(b,b′)=e∈E∑​​a∈{0,1}∑​U((b,e),a)Q(a,x)​​a∈{0,1}∑​U((b′,e),a)Q(a,x)​​.

There is no isometry or normalization assumption on UUU, and neither finite set is required to be nonempty. In particular, when EEE is empty the sum is zero.

QAlgorithms.Papers.Independent.Intakethird100.P10.dilationEve

For arbitrary finite sets B,EB,EB,E, any complex matrix UUU with rows indexed by B×EB\times EB×E and Boolean columns, and any unitary 2×22\times22×2 matrix QQQ, this declaration defines a Boolean-indexed family of matrices on EEE. For x∈{0,1}x\in\{0,1\}x∈{0,1} and e,e′∈Ee,e'\in Ee,e′∈E, its entries are

dilationEve⁡(U,Q)x(e,e′)=∑b∈B(∑a∈{0,1}U((b,e),a)Q(a,x))(∑a∈{0,1}U((b,e′),a)Q(a,x))‾.\operatorname{dilationEve}(U,Q)_x(e,e') =\sum_{b\in B} \left(\sum_{a\in\{0,1\}}U((b,e),a)Q(a,x)\right) \overline{\left(\sum_{a\in\{0,1\}}U((b,e'),a)Q(a,x)\right)}.dilationEve(U,Q)x​(e,e′)=b∈B∑​​a∈{0,1}∑​U((b,e),a)Q(a,x)​​a∈{0,1}∑​U((b,e′),a)Q(a,x)​​.

No isometry or normalization assumption on UUU is imposed, and B,EB,EB,E may be empty. In particular, an empty BBB makes this sum zero.

QAlgorithms.Papers.Independent.Intakethird100.P10.dilationLetter

For arbitrary finite sets B,EB,EB,E, any complex matrix UUU with rows indexed by B×EB\times EB×E and columns indexed by {0,1}\{0,1\}{0,1}, any unitary 2×22\times22×2 matrix QQQ, and any Boolean value xxx, define vx(b,e)=∑a∈{0,1}U((b,e),a)Q(a,x)v_x(b,e)=\sum_{a\in\{0,1\}}U((b,e),a)Q(a,x)vx​(b,e)=∑a∈{0,1}​U((b,e),a)Q(a,x). The declaration defines a matrix on B×EB\times EB×E by

dilationLetter⁡(U,Q,x)((b,e),(b′,e′))=vx(b,e) vx(b′,e′)‾.\operatorname{dilationLetter}(U,Q,x)((b,e),(b',e')) =v_x(b,e)\,\overline{v_x(b',e')}.dilationLetter(U,Q,x)((b,e),(b′,e′))=vx​(b,e)vx​(b′,e′)​.

No isometry condition on UUU, normalization of vxv_xvx​, or nonemptiness of B,EB,EB,E is required; if either set is empty, this is a matrix on an empty index set.

QAlgorithms.Papers.Independent.Intakethird100.P10.erasureData

For every natural number ddd and every complex matrix ρ\rhoρ indexed by {0,…,d−1}\{0,\ldots,d-1\}{0,…,d−1}, this defines a matrix ρ~\widetilde{\rho}ρ​ indexed by those same indices together with an additional index ∅\varnothing∅, with entries ρ~a,b=ρa,b\widetilde{\rho}_{a,b}=\rho_{a,b}ρ​a,b​=ρa,b​ for ordinary indices a,ba,ba,b, and zero whenever either index is ∅\varnothing∅. No positivity or normalization is required of ρ\rhoρ. For d=0d=0d=0, the resulting matrix has only the additional index and is zero.

QAlgorithms.Papers.Independent.Intakethird100.P10.frozenSet

For finite types B,EB,EB,E with decidable equality, a wiretap datum VVV, a natural number nnn, and a real β\betaβ, this defines the set of indices i<2ni<2^ni<2n satisfying

1−τ<f(Wn,i,0E,Wn,i,1E)andτ≤f(Wn,i,0B,Wn,i,1B),τ=2−(2n)β.1-\tau<f(W_{n,i,0}^E,W_{n,i,1}^E) \quad\text{and}\quad \tau\le f(W_{n,i,0}^B,W_{n,i,1}^B), \qquad \tau=2^{-(2^n)^\beta}.1−τ<f(Wn,i,0E​,Wn,i,1E​)andτ≤f(Wn,i,0B​,Wn,i,1B​),τ=2−(2n)β.

The datum VVV consists of joint positive semidefinite trace-one matrices ρx\rho_xρx​, their marginals WxB(b,b′)=∑eρx((b,e),(b′,e))W_x^B(b,b')=\sum_e\rho_x((b,e),(b',e))WxB​(b,b′)=∑e​ρx​((b,e),(b′,e)) and WxE(e,e′)=∑bρx((b,e),(b,e′))W_x^E(e,e')=\sum_b\rho_x((b,e),(b,e'))WxE​(e,e′)=∑b​ρx​((b,e),(b,e′)), and a finite Kraus channel D(X)=∑jKjXKj†D(X)=\sum_jK_jXK_j^\daggerD(X)=∑j​Kj​XKj†​, ∑jKj†Kj=IB\sum_jK_j^\dagger K_j=I_B∑j​Kj†​Kj​=IB​, such that WxE=D(WxB)W_x^E=D(W_x^B)WxE​=D(WxB​), with W0E,W1EW_0^E,W_1^EW0E​,W1E​ commuting. For either marginal pair WWW, set N=2nN=2^nN=2n; Wn,i,xW_{n,i,x}Wn,i,x​ is indexed by (v,h)(v,h)(v,h), where vvv has iii Boolean entries and hhh has NNN entries in the marginal's index type, with

Wn,i,x((v,h),(v′,h′))={2−(N−1)∑u∈{0,1}Nui=x, u<i=v∏k<NWp(u)k(hk,hk′),v=v′,0,v≠v′.W_{n,i,x}((v,h),(v',h'))= \begin{cases} 2^{-(N-1)}\displaystyle\sum_{\substack{u\in\{0,1\}^N\\u_i=x,\ u_{<i}=v}} \prod_{k<N}W_{p(u)_k}(h_k,h'_k),&v=v',\\ 0,&v\ne v'. \end{cases}Wn,i,x​((v,h),(v′,h′))=⎩⎨⎧​2−(N−1)u∈{0,1}Nui​=x, u<i​=v​∑​k<N∏​Wp(u)k​​(hk​,hk′​),0,​v=v′,v=v′.​

Here p(u)k=∑jajkuj(mod2)p(u)_k=\sum_j a_{jk}u_j\pmod2p(u)k​=∑j​ajk​uj​(mod2), and ajk=1a_{jk}=1ajk​=1 exactly when each 111-bit of kkk's length-nnn binary representation occurs at a reversed 111-bit position of jjj's representation. The function fff is Re⁡tr⁡X YX\operatorname{Re}\operatorname{tr}\sqrt{\sqrt X\,Y\sqrt X}RetrX​YX​​. Equality at the BBB threshold is included, equality at the EEE threshold is excluded, and n=0n=0n=0 and arbitrary real β\betaβ are allowed.

QAlgorithms.Papers.Independent.Intakethird100.P10.goodSet

For any finite type HHH with decidable equality, any pair W=(W0,W1)W=(W_0,W_1)W=(W0​,W1​) of complex H×HH\times HH×H matrices, any natural number nnn, and any real β\betaβ, put N=2nN=2^nN=2n. The defined finite set consists exactly of indices i∈{0,…,N−1}i\in\{0,\ldots,N-1\}i∈{0,…,N−1} for which

Re⁡tr⁡Si,0 Si,1Si,0<2−Nβ.\operatorname{Re}\operatorname{tr} \sqrt{\sqrt{S_{i,0}}\,S_{i,1}\sqrt{S_{i,0}}} < 2^{-N^\beta}.RetrSi,0​​Si,1​Si,0​​​<2−Nβ.

All matrix square roots here mean the total real continuous functional calculus of t↦tt\mapsto\sqrt tt↦t​, which returns zero on non-Hermitian matrices. The matrices Si,xS_{i,x}Si,x​, indexed by Boolean prefixes p,qp,qp,q of length iii and words a,b∈HNa,b\in H^Na,b∈HN, are zero when p≠qp\ne qp=q, and otherwise have entries

Si,x((p,a),(p,b))=2−(N−1)∑u∈{0,1}Nui=xuj=pj (j<i)∏k=0N−1(Wvk)ak,bk,vk=∑j=0N−1ujcn(j,k)(mod2),S_{i,x}((p,a),(p,b)) = 2^{-(N-1)} \sum_{\substack{u\in\{0,1\}^{N}\\u_i=x\\u_j=p_j\ (j<i)}} \prod_{k=0}^{N-1}(W_{v_k})_{a_k,b_k}, \qquad v_k=\sum_{j=0}^{N-1}u_jc_n(j,k)\pmod2,Si,x​((p,a),(p,b))=2−(N−1)u∈{0,1}Nui​=xuj​=pj​ (j<i)​∑​k=0∏N−1​(Wvk​​)ak​,bk​​,vk​=j=0∑N−1​uj​cn​(j,k)(mod2),

where cn(j,k)=1c_n(j,k)=1cn​(j,k)=1 exactly when every nonzero bit ttt of kkk's length-nnn, most-significant-first binary representation corresponds to a nonzero bit n−1−tn-1-tn−1−t of jjj's representation. No density-matrix assumption on WWW or restriction on β\betaβ is imposed. At n=0n=0n=0, the sole index is included exactly when the same root-fidelity expression for W0,W1W_0,W_1W0​,W1​ is strictly below 1/21/21/2; if HHH is empty, every index is included.

QAlgorithms.Papers.Independent.Intakethird100.P10.informationSet

For finite types B,EB,EB,E with decidable equality, a wiretap datum VVV, a natural number nnn, and a real β\betaβ, this defines the set of indices i<2ni<2^ni<2n satisfying both

1−τ<f(Wn,i,0E,Wn,i,1E)andf(Wn,i,0B,Wn,i,1B)<τ,τ=2−(2n)β.1-\tau<f(W_{n,i,0}^E,W_{n,i,1}^E) \quad\text{and}\quad f(W_{n,i,0}^B,W_{n,i,1}^B)<\tau, \qquad \tau=2^{-(2^n)^\beta}.1−τ<f(Wn,i,0E​,Wn,i,1E​)andf(Wn,i,0B​,Wn,i,1B​)<τ,τ=2−(2n)β.

Here VVV supplies joint positive semidefinite trace-one matrices ρx\rho_xρx​, their marginals WxB(b,b′)=∑eρx((b,e),(b′,e))W_x^B(b,b')=\sum_e\rho_x((b,e),(b',e))WxB​(b,b′)=∑e​ρx​((b,e),(b′,e)) and WxE(e,e′)=∑bρx((b,e),(b,e′))W_x^E(e,e')=\sum_b\rho_x((b,e),(b,e'))WxE​(e,e′)=∑b​ρx​((b,e),(b,e′)), and a finite Kraus channel D(X)=∑jKjXKj†D(X)=\sum_jK_jXK_j^\daggerD(X)=∑j​Kj​XKj†​, ∑jKj†Kj=IB\sum_jK_j^\dagger K_j=I_B∑j​Kj†​Kj​=IB​, satisfying WxE=D(WxB)W_x^E=D(W_x^B)WxE​=D(WxB​); the two EEE marginals commute. For either marginal pair WWW, write N=2nN=2^nN=2n: its synthesized matrix has entries

Wn,i,x((v,h),(v′,h′))={2−(N−1)∑u∈{0,1}Nui=x, u<i=v∏k<NWp(u)k(hk,hk′),v=v′,0,v≠v′,W_{n,i,x}((v,h),(v',h'))= \begin{cases} 2^{-(N-1)}\displaystyle\sum_{\substack{u\in\{0,1\}^N\\u_i=x,\ u_{<i}=v}} \prod_{k<N}W_{p(u)_k}(h_k,h'_k),&v=v',\\ 0,&v\ne v', \end{cases}Wn,i,x​((v,h),(v′,h′))=⎩⎨⎧​2−(N−1)u∈{0,1}Nui​=x, u<i​=v​∑​k<N∏​Wp(u)k​​(hk​,hk′​),0,​v=v′,v=v′,​

where v,v′v,v'v,v′ have length iii, h,h′h,h'h,h′ are length-NNN words over the marginal's index type, and p(u)k=∑jajkuj(mod2)p(u)_k=\sum_j a_{jk}u_j\pmod2p(u)k​=∑j​ajk​uj​(mod2), with ajk=1a_{jk}=1ajk​=1 exactly when every 111-bit of the length-nnn binary representation of kkk is a 111-bit at the reversed position in that of jjj. Finally f(X,Y)=Re⁡tr⁡X YXf(X,Y)=\operatorname{Re}\operatorname{tr}\sqrt{\sqrt X\,Y\sqrt X}f(X,Y)=RetrX​YX​​. Both inequalities are strict; n=0n=0n=0 and every real β\betaβ are included.

QAlgorithms.Papers.Independent.Intakethird100.P10.leafProbability

For every natural number nnn, including n=0n=0n=0, and every finite subset SSS of {0,…,2n−1}\{0,\ldots,2^n-1\}{0,…,2n−1}, this defines the real number ∣S∣/2n|S|/2^n∣S∣/2n. Thus the empty set has value 000, the full set has value 111, and the denominator is always positive; when n=0n=0n=0, the ambient index set has one element.

QAlgorithms.Papers.Independent.Intakethird100.P10.minusChannel

For any type HHH and any Boolean-indexed family WWW of complex matrices indexed by HHH, define a Boolean-indexed family MMM of matrices indexed by H×HH\times HH×H by

M(false)=12(W(false)⊗W(false)+W(true)⊗W(true)),M(true)=12(W(false)⊗W(true)+W(true)⊗W(false)).M(\mathrm{false})=\frac12\bigl(W(\mathrm{false})\otimes W(\mathrm{false}) +W(\mathrm{true})\otimes W(\mathrm{true})\bigr),\qquad M(\mathrm{true})=\frac12\bigl(W(\mathrm{false})\otimes W(\mathrm{true}) +W(\mathrm{true})\otimes W(\mathrm{false})\bigr).M(false)=21​(W(false)⊗W(false)+W(true)⊗W(true)),M(true)=21​(W(false)⊗W(true)+W(true)⊗W(false)).

The tensor product has entries (A⊗B)(a,c),(b,d)=Aa,bBc,d(A\otimes B)_{(a,c),(b,d)}=A_{a,b}B_{c,d}(A⊗B)(a,c),(b,d)​=Aa,b​Bc,d​. No positivity, trace, normalization, finiteness, or nonemptiness hypothesis is required.

QAlgorithms.Papers.Independent.Intakethird100.P10.polarCoefficient

For every natural number nnn, put N=2nN=2^nN=2n, and let j,k∈{0,…,N−1}j,k\in\{0,\ldots,N-1\}j,k∈{0,…,N−1}. Write bt(a)b_t(a)bt​(a) for bit ttt of the length-nnn binary representation of aaa, ordered from most significant to least significant, so that a=∑t=0n−1bt(a)2n−1−ta=\sum_{t=0}^{n-1}b_t(a)2^{n-1-t}a=∑t=0n−1​bt​(a)2n−1−t. The Boolean coefficient cn(j,k)c_n(j,k)cn​(j,k) is true exactly when

∀t∈{0,…,n−1},bt(k)=1 ⟹ bn−1−t(j)=1.\forall t\in\{0,\ldots,n-1\},\qquad b_t(k)=1\ \Longrightarrow\ b_{n-1-t}(j)=1.∀t∈{0,…,n−1},bt​(k)=1 ⟹ bn−1−t​(j)=1.

It is false otherwise. For n=0n=0n=0, the only indices are j=k=0j=k=0j=k=0, and the condition is vacuous, so the coefficient is true.

QAlgorithms.Papers.Independent.Intakethird100.P10.polarOperator

For every natural number nnn, define a complex matrix indexed by Boolean strings of length 2n2^n2n by

(Pn)v,u={1,v=Tn(u),0,otherwise.(P_n)_{v,u}= \begin{cases} 1,&v=T_n(u),\\ 0,&\text{otherwise}. \end{cases}(Pn​)v,u​={1,0,​v=Tn​(u),otherwise.​

For indices 0≤j,k<2n0\le j,k<2^n0≤j,k<2n, let jt,ktj_t,k_tjt​,kt​ denote their nnn-bit binary digits, with position 000 the most significant digit, and set aj,k=1a_{j,k}=1aj,k​=1 exactly when, for every 0≤t<n0\le t<n0≤t<n, kt=1k_t=1kt​=1 implies jn−1−t=1j_{n-1-t}=1jn−1−t​=1. Then Tn(u)kT_n(u)_kTn​(u)k​ is the parity of the number of indices jjj for which uj=1u_j=1uj​=1 and aj,k=1a_{j,k}=1aj,k​=1. When n=0n=0n=0, the digit condition is vacuous, and this transformation is the identity on the single Boolean wire.

QAlgorithms.Papers.Independent.Intakethird100.P10.polarThreshold

For every natural number nnn and every real number β\betaβ, define

τn,β=2−(2n)β,\tau_{n,\beta}=2^{-(2^n)^\beta},τn,β​=2−(2n)β,

using real exponentiation for the powers with real exponents. There is no restriction on β\betaβ. The threshold lies strictly between 000 and 111; for n=0n=0n=0, it equals 1/21/21/2 for every β\betaβ.

QAlgorithms.Papers.Independent.Intakethird100.P10.polarWord

For every natural number nnn, put N=2nN=2^nN=2n, and let u:{0,…,N−1}→{0,1}u:\{0,\ldots,N-1\}\to\{0,1\}u:{0,…,N−1}→{0,1} be a Boolean word. Its defined transformed word vvv has

vk=(∑j=0N−1ujcn(j,k)) mod 2,v_k=\left(\sum_{j=0}^{N-1}u_j c_n(j,k)\right)\bmod 2,vk​=(j=0∑N−1​uj​cn​(j,k))mod2,

where cn(j,k)=1c_n(j,k)=1cn​(j,k)=1 exactly when every nonzero bit ttt of kkk's length-nnn, most-significant-first binary representation corresponds to a nonzero bit n−1−tn-1-tn−1−t of jjj's representation, and cn(j,k)=0c_n(j,k)=0cn​(j,k)=0 otherwise. Equivalently, vkv_kvk​ is true exactly when an odd number of indices jjj have both uju_juj​ true and this bit condition satisfied. For n=0n=0n=0, the word has one coordinate and the transformation is the identity.

QAlgorithms.Papers.Independent.Intakethird100.P10.poorSet

For any finite type HHH with decidable equality, any pair W0,W1W_0,W_1W0​,W1​ of complex matrices indexed by HHH, any natural number nnn, and any real β\betaβ, this defines the finite set

{i∈{0,…,2n−1}:1−2−(2n)β<f(Wn,i,0,Wn,i,1)},\left\{i\in\{0,\ldots,2^n-1\}:1-2^{-(2^n)^\beta}<f(W_{n,i,0},W_{n,i,1})\right\},{i∈{0,…,2n−1}:1−2−(2n)β<f(Wn,i,0​,Wn,i,1​)},

where f(X,Y)=Re⁡tr⁡X YXf(X,Y)=\operatorname{Re}\operatorname{tr}\sqrt{\sqrt X\,Y\sqrt X}f(X,Y)=RetrX​YX​​, with square roots interpreted through the continuous functional calculus using the real square-root function. To specify the synthesized matrices, put N=2nN=2^nN=2n, and let bt(j)b_t(j)bt​(j) be bit ttt of the length-nnn, most-significant-bit-first binary representation of jjj. Define ajk=1a_{jk}=1ajk​=1 exactly when bt(k)=1b_t(k)=1bt​(k)=1 implies bn−1−t(j)=1b_{n-1-t}(j)=1bn−1−t​(j)=1 for every t<nt<nt<n, and define p(u)k=∑jajkuj(mod2)p(u)_k=\sum_j a_{jk}u_j\pmod 2p(u)k​=∑j​ajk​uj​(mod2) for a Boolean word uuu of length NNN. The matrix Wn,i,xW_{n,i,x}Wn,i,x​ is indexed by pairs (v,h)(v,h)(v,h), with vvv a Boolean word of length iii and h∈HNh\in H^Nh∈HN, and its entries are

Wn,i,x((v,h),(v′,h′))={2−(N−1)∑u∈{0,1}Nui=x, u<i=v∏k=0N−1Wp(u)k(hk,hk′),v=v′,0,v≠v′.W_{n,i,x}((v,h),(v',h'))= \begin{cases} 2^{-(N-1)}\displaystyle\sum_{\substack{u\in\{0,1\}^N\\u_i=x,\ u_{<i}=v}} \prod_{k=0}^{N-1}W_{p(u)_k}(h_k,h'_k),&v=v',\\ 0,&v\ne v'. \end{cases}Wn,i,x​((v,h),(v′,h′))=⎩⎨⎧​2−(N−1)u∈{0,1}Nui​=x, u<i​=v​∑​k=0∏N−1​Wp(u)k​​(hk​,hk′​),0,​v=v′,v=v′.​

There is no density-matrix assumption on WWW and no restriction on β\betaβ. The case n=0n=0n=0 is included and has one possible index.

QAlgorithms.Papers.Independent.Intakethird100.P10.pureCQ

For any type HHH and any two complex coordinate vectors ψ0,ψ1\psi_0,\psi_1ψ0​,ψ1​ indexed by HHH, define a Boolean-indexed family of complex matrices by

W(x)a,b=ψx(a) ψx(b)‾,ψfalse=ψ0,ψtrue=ψ1.W(x)_{a,b}=\psi_x(a)\,\overline{\psi_x(b)}, \qquad \psi_{\mathrm{false}}=\psi_0,\quad \psi_{\mathrm{true}}=\psi_1.W(x)a,b​=ψx​(a)ψx​(b)​,ψfalse​=ψ0​,ψtrue​=ψ1​.

No normalization, finiteness, or nonemptiness assumption on HHH is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P10.randomizationSet

For finite types B,EB,EB,E with decidable equality, a wiretap datum VVV, any natural number nnn, and any real β\betaβ, this defines the set of indices i<2ni<2^ni<2n satisfying

f(Wn,i,0E,Wn,i,1E)≤1−τandf(Wn,i,0B,Wn,i,1B)<τ,τ=2−(2n)β.f(W_{n,i,0}^E,W_{n,i,1}^E)\le1-\tau \quad\text{and}\quad f(W_{n,i,0}^B,W_{n,i,1}^B)<\tau, \qquad \tau=2^{-(2^n)^\beta}.f(Wn,i,0E​,Wn,i,1E​)≤1−τandf(Wn,i,0B​,Wn,i,1B​)<τ,τ=2−(2n)β.

The datum VVV includes joint positive semidefinite trace-one matrices ρx\rho_xρx​, marginals WxB(b,b′)=∑eρx((b,e),(b′,e))W_x^B(b,b')=\sum_e\rho_x((b,e),(b',e))WxB​(b,b′)=∑e​ρx​((b,e),(b′,e)) and WxE(e,e′)=∑bρx((b,e),(b,e′))W_x^E(e,e')=\sum_b\rho_x((b,e),(b,e'))WxE​(e,e′)=∑b​ρx​((b,e),(b,e′)), and a finite Kraus channel D(X)=∑jKjXKj†D(X)=\sum_jK_jXK_j^\daggerD(X)=∑j​Kj​XKj†​, ∑jKj†Kj=IB\sum_jK_j^\dagger K_j=I_B∑j​Kj†​Kj​=IB​, taking WxBW_x^BWxB​ to WxEW_x^EWxE​; the two EEE marginals commute. For either marginal pair WWW, with N=2nN=2^nN=2n, its synthesized matrices are

Wn,i,x((v,h),(v′,h′))={2−(N−1)∑u∈{0,1}Nui=x, u<i=v∏k<NWp(u)k(hk,hk′),v=v′,0,v≠v′,W_{n,i,x}((v,h),(v',h'))= \begin{cases} 2^{-(N-1)}\displaystyle\sum_{\substack{u\in\{0,1\}^N\\u_i=x,\ u_{<i}=v}} \prod_{k<N}W_{p(u)_k}(h_k,h'_k),&v=v',\\ 0,&v\ne v', \end{cases}Wn,i,x​((v,h),(v′,h′))=⎩⎨⎧​2−(N−1)u∈{0,1}Nui​=x, u<i​=v​∑​k<N∏​Wp(u)k​​(hk​,hk′​),0,​v=v′,v=v′,​

where vvv is a Boolean word of length iii, hhh is a length-NNN word over the marginal's index type, and p(u)k=∑jajkuj(mod2)p(u)_k=\sum_j a_{jk}u_j\pmod2p(u)k​=∑j​ajk​uj​(mod2); ajk=1a_{jk}=1ajk​=1 exactly when every 111-bit of kkk's length-nnn binary representation is present at the reversed position in that of jjj. Here f(X,Y)=Re⁡tr⁡X YXf(X,Y)=\operatorname{Re}\operatorname{tr}\sqrt{\sqrt X\,Y\sqrt X}f(X,Y)=RetrX​YX​​. The EEE threshold equality is included and the BBB threshold equality is excluded; n=0n=0n=0 is allowed.

QAlgorithms.Papers.Independent.Intakethird100.P10.symmetricCoherentInformation

Let BBB be a finite set with decidable equality and EEE any finite set, let UUU be a complex matrix with rows indexed by B×EB\times EB×E and Boolean columns, and let QQQ be a unitary 2×22\times22×2 matrix. Define vx(b,e)=∑a∈{0,1}U((b,e),a)Q(a,x)v_x(b,e)=\sum_{a\in\{0,1\}}U((b,e),a)Q(a,x)vx​(b,e)=∑a∈{0,1}​U((b,e),a)Q(a,x), and define a matrix on {0,1}×B\{0,1\}\times B{0,1}×B by Ω((x,b),(y,b′))=12∑evx(b,e)vy(b′,e)‾\Omega((x,b),(y,b'))=\tfrac12\sum_e v_x(b,e)\overline{v_y(b',e)}Ω((x,b),(y,b′))=21​∑e​vx​(b,e)vy​(b′,e)​. Define its matrix on BBB obtained by summing the Boolean diagonal blocks by ΩB(b,b′)=∑x∈{0,1}Ω((x,b),(x,b′))\Omega_B(b,b')=\sum_{x\in\{0,1\}}\Omega((x,b),(x,b'))ΩB​(b,b′)=∑x∈{0,1}​Ω((x,b),(x,b′)). Writing S(X)=−Re⁡tr⁡(f(X))S(X)=-\operatorname{Re}\operatorname{tr}(f(X))S(X)=−Retr(f(X)), with f(t)=tlog⁡2tf(t)=t\log_2tf(t)=tlog2​t applied by continuous functional calculus, the declaration defines

symmetricCoherentInformation⁡(U,Q)=S(ΩB)−S(Ω).\operatorname{symmetricCoherentInformation}(U,Q)=S(\Omega_B)-S(\Omega).symmetricCoherentInformation(U,Q)=S(ΩB​)−S(Ω).

Neither nonemptiness of B,EB,EB,E nor an isometry or normalization condition on UUU is required.

QAlgorithms.Papers.Independent.Intakethird100.P10.symmetricInformation

For any finite type HHH with decidable equality and any pair W=(W0,W1)W=(W_0,W_1)W=(W0​,W1​) of complex H×HH\times HH×H matrices, define

Isym(W)=S ⁣(W0+W12)−S(W0)+S(W1)2,S(M)=−Re⁡tr⁡(f[M]),f(t)=tlog⁡2t,I_{\mathrm{sym}}(W) = S\!\left(\frac{W_0+W_1}{2}\right) -\frac{S(W_0)+S(W_1)}{2}, \qquad S(M)=-\operatorname{Re}\operatorname{tr}\bigl(f[M]\bigr), \quad f(t)=t\log_2 t,Isym​(W)=S(2W0​+W1​​)−2S(W0​)+S(W1​)​,S(M)=−Retr(f[M]),f(t)=tlog2​t,

where f[M]f[M]f[M] denotes the total real continuous functional calculus used for matrices. The real logarithm is total, with its value at zero defined as zero, and this functional calculus returns the zero matrix when its Hermitian applicability condition fails. No assumption that either WxW_xWx​ is a density matrix is required. Empty HHH is allowed, in which case the defined value is 000.

QAlgorithms.Papers.Independent.Intakethird100.P10.synthesized

For any finite type HHH with decidable equality, any pair W=(W0,W1)W=(W_0,W_1)W=(W0​,W1​) of complex H×HH\times HH×H matrices, any natural number nnn, and any i∈{0,…,N−1}i\in\{0,\ldots,N-1\}i∈{0,…,N−1}, where N=2nN=2^nN=2n, define a pair of matrices S0,S1S_0,S_1S0​,S1​ indexed by {0,1} i×HN\{0,1\}^{\,i}\times H^N{0,1}i×HN. For x∈{0,1}x\in\{0,1\}x∈{0,1}, prefixes p,q∈{0,1} ip,q\in\{0,1\}^{\,i}p,q∈{0,1}i, and a,b∈HNa,b\in H^Na,b∈HN, their entries are

Sx((p,a),(q,b))={2−(N−1)∑u∈{0,1}Nui=xuj=pj (0≤j<i)∏k=0N−1(Wvk)ak,bk,p=q,0,p≠q.S_x((p,a),(q,b)) = \begin{cases} \displaystyle 2^{-(N-1)} \sum_{\substack{u\in\{0,1\}^{N}\\u_i=x\\u_j=p_j\ (0\le j<i)}} \prod_{k=0}^{N-1}(W_{v_k})_{a_k,b_k}, & p=q,\\[6pt] 0,&p\ne q. \end{cases}Sx​((p,a),(q,b))=⎩⎨⎧​2−(N−1)u∈{0,1}Nui​=xuj​=pj​ (0≤j<i)​∑​k=0∏N−1​(Wvk​​)ak​,bk​​,0,​p=q,p=q.​

Here vk=∑j=0N−1ujcn(j,k)(mod2)v_k=\sum_{j=0}^{N-1}u_jc_n(j,k)\pmod 2vk​=∑j=0N−1​uj​cn​(j,k)(mod2), and cn(j,k)=1c_n(j,k)=1cn​(j,k)=1 exactly when every nonzero bit ttt of kkk's length-nnn, most-significant-first binary representation corresponds to a nonzero bit n−1−tn-1-tn−1−t of jjj's representation. The prefix condition is vacuous at i=0i=0i=0. No density-matrix hypothesis on WWW is imposed. For n=0n=0n=0, N=1N=1N=1, i=0i=0i=0, and Sx(((),a),((),b))=(Wx)a0,b0S_x(((),a),((),b))=(W_x)_{a_0,b_0}Sx​(((),a),((),b))=(Wx​)a0​,b0​​.

QAlgorithms.Papers.Independent.Intakethird100.P10.synthesizedDegrader

For any finite type HHH with decidable equality, any type KKK, any function DDD from complex H×HH\times HH×H matrices to complex K×KK\times KK×K matrices, any natural number nnn, and any i∈{0,…,N−1}i\in\{0,\ldots,N-1\}i∈{0,…,N−1}, where N=2nN=2^nN=2n, define a function D~\widetilde DD from matrices on {0,1} i×HN\{0,1\}^{\,i}\times H^N{0,1}i×HN to matrices on {0,1} i×KN\{0,1\}^{\,i}\times K^N{0,1}i×KN by

(D~X)((p,y),(q,y′))=∑a∈HN∑b∈HNX((p,a),(q,b))∏k=0N−1(D(Eak,bk))yk,yk′,(\widetilde D X)((p,y),(q,y')) = \sum_{a\in H^N}\sum_{b\in H^N} X((p,a),(q,b)) \prod_{k=0}^{N-1} \bigl(D(E_{a_k,b_k})\bigr)_{y_k,y'_k},(DX)((p,y),(q,y′))=a∈HN∑​b∈HN∑​X((p,a),(q,b))k=0∏N−1​(D(Eak​,bk​​))yk​,yk′​​,

where p,qp,qp,q are Boolean prefixes of length iii, y,y′∈KNy,y'\in K^Ny,y′∈KN, and Eh,h′E_{h,h'}Eh,h′​ has entry 111 at (h,h′)(h,h')(h,h′) and zero elsewhere. No linearity, positivity, or channel assumption on DDD, and no finiteness or decidable-equality assumption on KKK, is required. For n=0n=0n=0, the prefixes are empty and the product has one factor.

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