Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for On the distributed compression of quantum information

Definition
QAlgorithms_Papers_Independent_Intakethird100_P11

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

prove2quantum

The objects the formalized statements of On the distributed compression of quantum information (arXiv:quant-ph/0403042) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • DistributedBlockCode: PDF pp.2–3, (2)–(3): an actual blind distributed code. Only the two compressed quantum registers are transmitted; there is no other resource register or message.
  • DistributedBlockCode.output: PDF p.2: independent encoders followed by the joint receiver's decoder. The accepted both applies each local map to all matrix blocks, including coherences.
  • sourceWordVector: PDF p.2: tensor word with Alice's and Bob's factors grouped by owner. The empty word is the unit vector of the singleton zero-copy space.
  • letterFidelity: PDF p.2, (2): pure-input squared fidelity, including possible mixed output.
  • finiteFidelity: PDF p.2: the i.i.d. finite source average in (2); repeated labels remain labels.
  • DistributedAchievable: PDF p.3, (3): every positive accuracy and rate slack, then every sufficiently large positive integer block. Dimension one permits zero communication.
  • FiniteAchievable: PDF pp.2–3: achievable rates for the given finite pure ensemble. Validity of p and the pure vectors is supplied in theorem binders, or by explicit constructions.
  • sourceMean: PDF p.2: average density, using the frozen finite-ensemble definition after enumerating the abstract finite label type.
  • EnsembleIrreducible: PDF p.5: no nontrivial orthogonal partition of the positive-probability support. A singleton support is irreducible; a zero-weight letter cannot connect components.
  • bellSourceVectors: PDF p.8, (13), in the order phi+, phi-, psi+, psi-. Both registers are qubits.
  • erasureBasis: PDF p.5, (8): four orthonormal erasure-code words. Fin 4 denotes the fixed most-significant-bit-first two-qubit basis 00,01,10,11 on each party.
  • erasureVector: PDF p.5: the isometric image of the first column of U(4), sampling all pure superpositions in the code subspace, rather than only the four code words.
  • erasureFidelity: PDF p.5 with p.2, (2): Haar-uniform pure ensemble in the whole code subspace, with independent U(4) samples at the n source positions.
  • hiddenAlice: PDF p.11, (23)–(25), in letter order 1,2,3.
  • hiddenBob: PDF p.11, (26)–(28): Bob's original qutrit vectors.
  • hiddenBobAfter: PDF p.11: Bob's vectors after projecting onto |2> and its complement and resetting the |2> branch to |0>. The output on the third letter is pure |0>.
  • hiddenJoint: PDF p.11: the original equiprobable product source, before Bob's local reset.
  • hiddenRateA: PDF p.11: Alice's actual Schumacher rate, from the average marginal density.
  • hiddenRateB: PDF p.11: Bob's actual Schumacher rate after the project/reset operation.
  • hiddenComparison: PDF p.11, (30): the expression from III.3, evaluated on the original source.
  • hybridVector: PDF p.12, (31)–(32): the displayed opposite-parity family. False/true denote 0/1 for the source label and both computational qubits. This is a restriction of scope: the printed universal local-unitary reduction to this family is false.
  • hybridQ: PDF p.13, (33): Alice's measurement-outcome probabilities.
  • hybridConditional: PDF p.13, (33): conditional density on Bob. A zero-probability outcome is assigned |0><0|; it contributes zero to every averaged quantity.
  • hybridChi: PDF p.13: E' Holevo information; H(E') in (36)–(38) is instead H(q).
  • hybridPiggyback: PDF p.13, (35): chi(E'')=1+S(E_B)-S(E_AB), since Alice is a qubit.
  • hybridRateA: PDF p.13, (37): exact hybrid Alice rate; its denominator is always at least one for normalized p and coefficient pairs, including zero probabilities.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_Ensemble
import Definitions.Def_QAlgorithms_Defs_SchurUnitaryEnsembles
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P19
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P27

/-!
Ahn, Doherty, Hayden and Winter, On the distributed compression of quantum information,
quant-ph/0403042v3. Codes reconstruct the supplied pure ensemble in average squared
fidelity, not with a reference. All rates and entropies are in bits/qubits per signal.
Only pinned accepted modules are imported. The universal local normal form on p.12
is false and is deferred in HARD.md; the hybrid result retains its displayed-family scope.
Compilation is pending the mechanical batch handoff.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P11

open MeasureTheory Filter
open scoped BigOperators ComplexConjugate Matrix
open QAlgorithms.Papers.Independent.Dependency50.P19 (both)
open QAlgorithms.Papers.Independent.Information100.P27 (shannon)

attribute [local instance] Classical.propDecidable Classical.decEq

/-- PDF pp.2–3, (2)–(3): an actual blind distributed code. Only the two compressed
quantum registers are transmitted; there is no other resource register or message. -/
structure DistributedBlockCode (A B : Type*) [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (n da db : ℕ) where
  length_pos : 0 < n
  alice_dim_pos : 0 < da
  bob_dim_pos : 0 < db
  encodeA : Superop (Fin n → A) (Fin da)
  encodeB : Superop (Fin n → B) (Fin db)
  decode : Superop (Fin da × Fin db) ((Fin n → A) × (Fin n → B))
  alice_channel : IsChannel encodeA
  bob_channel : IsChannel encodeB
  decoder_channel : IsChannel decode

/-- PDF p.2: independent encoders followed by the joint receiver's decoder. The
accepted `both` applies each local map to all matrix blocks, including coherences. -/
def DistributedBlockCode.output {A B : Type*} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] {n da db : ℕ}
    (c : DistributedBlockCode A B n da db) :
    Superop ((Fin n → A) × (Fin n → B)) ((Fin n → A) × (Fin n → B)) :=
  c.decode ∘ both c.encodeA c.encodeB

/-- PDF p.2: tensor word with Alice's and Bob's factors grouped by owner. The
empty word is the unit vector of the singleton zero-copy space. -/
def sourceWordVector {I A B : Type*} (v : I → EuclideanSpace ℂ (A × B))
    {n : ℕ} (w : Fin n → I) : EuclideanSpace ℂ ((Fin n → A) × (Fin n → B)) :=
  WithLp.toLp 2 fun x => ∏ t, v (w t) (x.1 t, x.2 t)

/-- PDF p.2, (2): pure-input squared fidelity, including possible mixed output. -/
def letterFidelity {I A B : Type*} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (v : I → EuclideanSpace ℂ (A × B))
    {n da db : ℕ} (c : DistributedBlockCode A B n da db) (w : Fin n → I) : ℝ :=
  effectProb (pureDensity (sourceWordVector v w)) (c.output (pureDensity (sourceWordVector v w)))

/-- PDF p.2: the i.i.d. finite source average in (2); repeated labels remain labels. -/
def finiteFidelity {I A B : Type*} [Fintype I] [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (p : I → ℝ) (v : I → EuclideanSpace ℂ (A × B))
    {n da db : ℕ} (c : DistributedBlockCode A B n da db) : ℝ :=
  ∑ w : Fin n → I, (∏ t, p (w t)) * letterFidelity v c w

/-- PDF p.3, (3): every positive accuracy and rate slack, then every sufficiently
large positive integer block. Dimension one permits zero communication. -/
def DistributedAchievable {A B : Type*} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B]
    (fidelity : {n da db : ℕ} → DistributedBlockCode A B n da db → ℝ) (ra rb : ℝ) : Prop :=
  ∀ delta > (0 : ℝ), ∀ eps > (0 : ℝ), ∃ N : ℕ, ∀ n : ℕ, N < n →
    ∃ da db : ℕ, ∃ c : DistributedBlockCode A B n da db,
      1 - eps ≤ fidelity c ∧
      Real.logb 2 da / n ≤ ra + delta ∧ Real.logb 2 db / n ≤ rb + delta

/-- PDF pp.2–3: achievable rates for the given finite pure ensemble. Validity of
p and the pure vectors is supplied in theorem binders, or by explicit constructions. -/
def FiniteAchievable {I A B : Type*} [Fintype I] [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (p : I → ℝ) (v : I → EuclideanSpace ℂ (A × B))
    (ra rb : ℝ) : Prop :=
  DistributedAchievable (fun c => finiteFidelity p v c) ra rb

/-- PDF p.2: average density, using the frozen finite-ensemble definition after
enumerating the abstract finite label type. -/
def sourceMean {I X : Type*} [Fintype I] (p : I → ℝ)
    (v : I → EuclideanSpace ℂ X) : Matrix X X ℂ :=
  ensembleDensity (fun j => p ((Fintype.equivFin I).symm j))
    (fun j => v ((Fintype.equivFin I).symm j))

/-- PDF p.5: no nontrivial orthogonal partition of the positive-probability support.
A singleton support is irreducible; a zero-weight letter cannot connect components. -/
def EnsembleIrreducible {I X : Type*} [Fintype X] (p : I → ℝ)
    (v : I → EuclideanSpace ℂ X) : Prop :=
  (∃ i, 0 < p i) ∧ ∀ S : Set I,
    (∃ i, 0 < p i ∧ i ∈ S) → (∃ j, 0 < p j ∧ j ∉ S) →
    ∃ i j, 0 < p i ∧ 0 < p j ∧ i ∈ S ∧ j ∉ S ∧ inner ℂ (v i) (v j) ≠ 0

/-- PDF p.8, (13), in the order phi+, phi-, psi+, psi-. Both registers are qubits. -/
def bellSourceVectors : Fin 4 → EuclideanSpace ℂ (Bool × Bool) :=
  let s : ℂ := (Real.sqrt 2 : ℝ)⁻¹
  ![s • (ket (false, false) + ket (true, true)),
    s • (ket (false, false) - ket (true, true)),
    s • (ket (false, true) + ket (true, false)),
    s • (ket (false, true) - ket (true, false))]

/-- PDF p.5, (8): four orthonormal erasure-code words. Fin 4 denotes the fixed
most-significant-bit-first two-qubit basis 00,01,10,11 on each party. -/
def erasureBasis : Fin 4 → EuclideanSpace ℂ (Fin 4 × Fin 4) :=
  let s : ℂ := (Real.sqrt 2 : ℝ)⁻¹
  ![s • (ket (0, 0) + ket (3, 3)), s • (ket (0, 3) + ket (3, 0)),
    s • (ket (1, 1) + ket (2, 2)), s • (ket (2, 1) + ket (1, 2))]

/-- PDF p.5: the isometric image of the first column of U(4), sampling all pure
superpositions in the code subspace, rather than only the four code words. -/
def erasureVector (U : Matrix.unitaryGroup (Fin 4) ℂ) :
    EuclideanSpace ℂ (Fin 4 × Fin 4) :=
  ∑ j : Fin 4, (U : Matrix (Fin 4) (Fin 4) ℂ) j 0 • erasureBasis j

/-- PDF p.5 with p.2, (2): Haar-uniform pure ensemble in the whole code subspace,
with independent U(4) samples at the n source positions. -/
def erasureFidelity {n da db : ℕ} (c : DistributedBlockCode (Fin 4) (Fin 4) n da db) : ℝ :=
  ∫ w : Fin n → Matrix.unitaryGroup (Fin 4) ℂ,
    letterFidelity erasureVector c w ∂Measure.pi (fun _ : Fin n => Schur.unitaryHaar 4)

/-- PDF p.11, (23)–(25), in letter order 1,2,3. -/
def hiddenAlice (alpha : ℝ) : Fin 3 → EuclideanSpace ℂ (Fin 2) :=
  ![ket 0, (Real.sqrt alpha : ℂ) • ket 0 + (Real.sqrt (1 - alpha) : ℂ) • ket 1, ket 1]

/-- PDF p.11, (26)–(28): Bob's original qutrit vectors. -/
def hiddenBob (beta : ℝ) : Fin 3 → EuclideanSpace ℂ (Fin 3) :=
  ![(Real.sqrt (1 - beta) : ℂ) • ket 0 + (Real.sqrt beta : ℂ) • ket 1,
    ket 1, (Real.sqrt (1 - beta) : ℂ) • ket 2 + (Real.sqrt beta : ℂ) • ket 0]

/-- PDF p.11: Bob's vectors after projecting onto |2> and its complement and
resetting the |2> branch to |0>. The output on the third letter is pure |0>. -/
def hiddenBobAfter (beta : ℝ) : Fin 3 → EuclideanSpace ℂ (Fin 3) :=
  ![hiddenBob beta 0, hiddenBob beta 1, ket 0]

/-- PDF p.11: the original equiprobable product source, before Bob's local reset. -/
def hiddenJoint (alpha beta : ℝ) (i : Fin 3) :=
  tensorVec (hiddenAlice alpha i) (hiddenBob beta i)

/-- PDF p.11: Alice's actual Schumacher rate, from the average marginal density. -/
def hiddenRateA (alpha : ℝ) : ℝ :=
  vnEntropy (sourceMean (fun _ : Fin 3 => (1 / 3 : ℝ)) (hiddenAlice alpha))

/-- PDF p.11: Bob's actual Schumacher rate after the project/reset operation. -/
def hiddenRateB (beta : ℝ) : ℝ :=
  vnEntropy (sourceMean (fun _ : Fin 3 => (1 / 3 : ℝ)) (hiddenBobAfter beta))

/-- PDF p.11, (30): the expression from III.3, evaluated on the original source. -/
def hiddenComparison (alpha beta : ℝ) : ℝ :=
  (hiddenRateA alpha +
    vnEntropy (sourceMean (fun _ : Fin 3 => (1 / 3 : ℝ)) (hiddenBob beta)) +
    vnEntropy (sourceMean (fun _ : Fin 3 => (1 / 3 : ℝ)) (hiddenJoint alpha beta))) / 2

/-- PDF p.12, (31)–(32): the displayed opposite-parity family. False/true denote
0/1 for the source label and both computational qubits. This is a restriction of
scope: the printed universal local-unitary reduction to this family is false. -/
def hybridVector (alpha beta : Bool → ℂ) (i : Bool) : EuclideanSpace ℂ (Bool × Bool) :=
  WithLp.toLp 2 fun x => if x.2 = Bool.xor x.1 i then
    (if x.1 then beta i else alpha i) else 0

/-- PDF p.13, (33): Alice's measurement-outcome probabilities. -/
def hybridQ (p : Bool → ℝ) (alpha beta : Bool → ℂ) (j : Bool) : ℝ :=
  ∑ i, p i * (if j then Complex.normSq (beta i) else Complex.normSq (alpha i))

/-- PDF p.13, (33): conditional density on Bob. A zero-probability outcome is
assigned |0><0|; it contributes zero to every averaged quantity. -/
def hybridConditional (p : Bool → ℝ) (alpha beta : Bool → ℂ) (j : Bool) :
    Matrix Bool Bool ℂ :=
  if hybridQ p alpha beta j = 0 then pureDensity (ket false) else
    ((hybridQ p alpha beta j)⁻¹ : ℂ) •
      ∑ i, ((p i * (if j then Complex.normSq (beta i) else Complex.normSq (alpha i)) : ℝ) : ℂ) •
        pureDensity (ket (Bool.xor j i))

/-- PDF p.13: E' Holevo information; H(E') in (36)–(38) is instead H(q). -/
def hybridChi (p : Bool → ℝ) (alpha beta : Bool → ℂ) : ℝ :=
  vnEntropy (∑ j, (hybridQ p alpha beta j : ℂ) • hybridConditional p alpha beta j) -
    ∑ j, hybridQ p alpha beta j * vnEntropy (hybridConditional p alpha beta j)

/-- PDF p.13, (35): chi(E'')=1+S(E_B)-S(E_AB), since Alice is a qubit. -/
def hybridPiggyback (p : Bool → ℝ) (alpha beta : Bool → ℂ) : ℝ :=
  let rho := sourceMean p (hybridVector alpha beta)
  1 + vnEntropy (traceLeft rho) - vnEntropy rho

/-- PDF p.13, (37): exact hybrid Alice rate; its denominator is always at least one
for normalized p and coefficient pairs, including zero probabilities. -/
def hybridRateA (p : Bool → ℝ) (alpha beta : Bool → ℂ) : ℝ :=
  let a := shannon (hybridQ p alpha beta) - hybridChi p alpha beta
  a / (a + hybridPiggyback p alpha beta)

section SanityExamples

example {I A B : Type*} (v : I → EuclideanSpace ℂ (A × B))
    (w : Fin 0 → I) (x : (Fin 0 → A) × (Fin 0 → B)) : sourceWordVector v w x = 1 := by
  simp [sourceWordVector]

example {I A B : Type*} (v : I → EuclideanSpace ℂ (A × B))
    (w : Fin 1 → I) (x : (Fin 1 → A) × (Fin 1 → B)) :
    sourceWordVector v w x = v (w 0) (x.1 0, x.2 0) := by
  simp [sourceWordVector]

example {A B : Type*} [Fintype A] [Fintype B] [DecidableEq A] [DecidableEq B]
    {n da db : ℕ} (c : DistributedBlockCode A B n da db)
    (X : Matrix ((Fin n → A) × (Fin n → B)) ((Fin n → A) × (Fin n → B)) ℂ) :
    c.output X = c.decode (both c.encodeA c.encodeB X) := rfl

example : hiddenBobAfter 0 2 = ket (0 : Fin 3) := by
  simp [hiddenBobAfter]

example (alpha beta : Bool → ℂ) : hybridVector alpha beta false (false, false) = alpha false := by
  simp [hybridVector]

example (alpha beta : Bool → ℂ) : hybridVector alpha beta true (true, false) = beta true := by
  simp [hybridVector]

example (p : Bool → ℝ) (alpha beta : Bool → ℂ) (j : Bool)
    (hq : hybridQ p alpha beta j = 0) :
    hybridConditional p alpha beta j = pureDensity (ket false) := by
  simp [hybridConditional, hq]

end SanityExamples

end QAlgorithms.Papers.Independent.Intakethird100.P11
Source
On the distributed compression of quantum information, arXiv:quant-ph/0403042v3, https://arxiv.org/abs/quant-ph/0403042v3
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P11.DistributedAchievable

For finite types A,BA,BA,B with decidable equality, an arbitrary real-valued function FFF assigning a value to every distributed block code at every choice of natural parameters n,dA,dBn,d_A,d_Bn,dA​,dB​, and real numbers rA,rBr_A,r_BrA​,rB​, distributed achievability means that for every δ>0\delta>0δ>0 and every ε>0\varepsilon>0ε>0, there exists a natural number NNN such that for every natural number n>Nn>Nn>N, there exist positive natural numbers dA,dBd_A,d_BdA​,dB​ and a code ccc of positive length nnn satisfying

1−ε≤F(c),log⁡2dAn≤rA+δ,log⁡2dBn≤rB+δ.1-\varepsilon\le F(c), \qquad \frac{\log_2d_A}{n}\le r_A+\delta, \qquad \frac{\log_2d_B}{n}\le r_B+\delta.1−ε≤F(c),nlog2​dA​​≤rA​+δ,nlog2​dB​​≤rB​+δ.

Such a code consists of encoders from matrices on AnA^nAn and BnB^nBn to matrices of dimensions dA,dBd_A,d_BdA​,dB​, and a decoder from their joint compressed space back to matrices on An×BnA^n\times B^nAn×Bn; every map must have a finite Kraus representation Φ(X)=∑jKjXKj†\Phi(X)=\sum_jK_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. The function FFF has no assumed bounds or other properties, and the rates have no separately imposed nonnegativity condition. The threshold uses the strict inequality N<nN<nN<n, so length zero is never required.

QAlgorithms.Papers.Independent.Intakethird100.P11.DistributedBlockCode

For finite types A,BA,BA,B with decidable equality and natural numbers n,dA,dBn,d_A,d_Bn,dA​,dB​, a distributed block code consists of proofs that n,dA,dB>0n,d_A,d_B>0n,dA​,dB​>0, two maps EA\mathcal E_AEA​ and EB\mathcal E_BEB​ from complex matrices indexed respectively by AnA^nAn and BnB^nBn to matrices indexed by {0,…,dA−1}\{0,\ldots,d_A-1\}{0,…,dA​−1} and {0,…,dB−1}\{0,\ldots,d_B-1\}{0,…,dB​−1}, and a map D\mathcal DD from matrices indexed by {0,…,dA−1}×{0,…,dB−1}\{0,\ldots,d_A-1\}\times\{0,\ldots,d_B-1\}{0,…,dA​−1}×{0,…,dB​−1} to matrices indexed by An×BnA^n\times B^nAn×Bn. Each of these three maps must admit a finite Kraus representation: for the corresponding input and output index sets, there are a natural number mmm and rectangular complex matrices K0,…,Km−1K_0,\ldots,K_{m-1}K0​,…,Km−1​ such that

∑j=0m−1Kj†Kj=I,Φ(X)=∑j=0m−1KjXKj†for every input matrix X.\sum_{j=0}^{m-1}K_j^\dagger K_j=I, \qquad \Phi(X)=\sum_{j=0}^{m-1}K_jXK_j^\dagger \quad\text{for every input matrix }X.j=0∑m−1​Kj†​Kj​=I,Φ(X)=j=0∑m−1​Kj​XKj†​for every input matrix X.

Here AnA^nAn means functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to AAA, and similarly for BnB^nBn. The parameters n,dA,dBn,d_A,d_Bn,dA​,dB​ may initially be any natural numbers, but no code exists when any of them is zero. No nonemptiness assumptions are imposed on A,BA,BA,B; the channel requirements also prevent such a code from existing if either type is empty.

QAlgorithms.Papers.Independent.Intakethird100.P11.DistributedBlockCode.output

For finite types A,BA,BA,B with decidable equality, natural numbers n,dA,dBn,d_A,d_Bn,dA​,dB​, and a distributed block code ccc, whose length and compressed dimensions are positive and whose encoders EA,EB\mathcal E_A,\mathcal E_BEA​,EB​ and decoder D\mathcal DD each have a finite Kraus representation Φ(X)=∑jKjXKj†\Phi(X)=\sum_jK_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I, the output map sends every matrix XXX indexed by An×BnA^n\times B^nAn×Bn to D(Z)\mathcal D(Z)D(Z). Explicitly, for compressed indices u,u′<dAu,u'<d_Au,u′<dA​ and v,v′<dBv,v'<d_Bv,v′<dB​,

Z(u,v),(u′,v′)=[EB ⁣([[EA ⁣([X(a,b),(a′,b′)]a,a′∈An)]u,u′]b,b′∈Bn)]v,v′.Z_{(u,v),(u',v')} = \left[ \mathcal E_B\!\left( \left[ \left[ \mathcal E_A\!\left( [X_{(a,b),(a',b')}]_{a,a'\in A^n} \right) \right]_{u,u'} \right]_{b,b'\in B^n} \right) \right]_{v,v'}.Z(u,v),(u′,v′)​=[EB​([[EA​([X(a,b),(a′,b′)​]a,a′∈An​)]u,u′​]b,b′∈Bn​)]v,v′​.

Thus the definition first applies Alice’s encoder to each matrix block with Bob’s two indices fixed, then applies Bob’s encoder to each block with Alice’s compressed indices fixed, and finally applies the decoder. It is defined on every complex input matrix, without requiring that matrix to be a density matrix.

QAlgorithms.Papers.Independent.Intakethird100.P11.EnsembleIrreducible

For an arbitrary label type III, a finite coordinate type XXX, arbitrary real weights pip_ipi​, and arbitrary vectors vi∈CXv_i\in\mathbb C^Xvi​∈CX, ensemble irreducibility means that at least one label has strictly positive weight, and every subset S⊆IS\subseteq IS⊆I containing a positive-weight label and excluding another positive-weight label has positive-weight labels i∈Si\in Si∈S and j∉Sj\notin Sj∈/S such that

⟨vi,vj⟩=∑x∈Xvi(x)‾vj(x)≠0.\langle v_i,v_j\rangle = \sum_{x\in X}\overline{v_i(x)}v_j(x)\ne0.⟨vi​,vj​⟩=x∈X∑​vi​(x)​vj​(x)=0.

The label type need not be finite, and no normalization of the vectors, nonnegativity of all weights, or condition on their sum is required. With exactly one positive-weight label, the subset condition is vacuous; with no positive-weight labels, the predicate is false.

QAlgorithms.Papers.Independent.Intakethird100.P11.FiniteAchievable

For finite types I,A,BI,A,BI,A,B, with decidable equality on A,BA,BA,B, arbitrary real weights pip_ipi​, arbitrary vectors vi∈CA×Bv_i\in\mathbb C^{A\times B}vi​∈CA×B, and real rates rA,rBr_A,r_BrA​,rB​, finite achievability means that for every δ,ε>0\delta,\varepsilon>0δ,ε>0, there is a natural number NNN such that for every natural number n>Nn>Nn>N, there exist positive dimensions dA,dBd_A,d_BdA​,dB​ and encoders EA,EB\mathcal E_A,\mathcal E_BEA​,EB​ and decoder D\mathcal DD satisfying

1−ε≤∑w∈In(∏t=0n−1pw(t))Re⁡tr⁡ ⁣(ρw D((EA⊗EB)(ρw))),log⁡2dAn≤rA+δ,log⁡2dBn≤rB+δ.1-\varepsilon\le \sum_{w\in I^n} \left(\prod_{t=0}^{n-1}p_{w(t)}\right) \operatorname{Re}\operatorname{tr}\!\left( \rho_w\,\mathcal D((\mathcal E_A\otimes\mathcal E_B)(\rho_w)) \right), \qquad \frac{\log_2d_A}{n}\le r_A+\delta, \qquad \frac{\log_2d_B}{n}\le r_B+\delta.1−ε≤w∈In∑​(t=0∏n−1​pw(t)​)Retr(ρw​D((EA​⊗EB​)(ρw​))),nlog2​dA​​≤rA​+δ,nlog2​dB​​≤rB​+δ.

Here Vw(a,b)=∏t=0n−1vw(t)(a(t),b(t))V_w(a,b)=\prod_{t=0}^{n-1}v_{w(t)}(a(t),b(t))Vw​(a,b)=∏t=0n−1​vw(t)​(a(t),b(t)), (ρw)x,y=Vw(x)Vw(y)‾(\rho_w)_{x,y}=V_w(x)\overline{V_w(y)}(ρw​)x,y​=Vw​(x)Vw​(y)​, the encoders map matrices on An,BnA^n,B^nAn,Bn to dimensions dA,dBd_A,d_BdA​,dB​, the tensor action applies them separately to those two sets of indices, and the decoder returns a matrix on An×BnA^n\times B^nAn×Bn. Each map must admit a finite Kraus representation Φ(X)=∑jKjXKj†\Phi(X)=\sum_jK_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. The definition imposes no probability or normalization conditions on ppp or vvv, and no separate sign restrictions on the rates.

QAlgorithms.Papers.Independent.Intakethird100.P11.bellSourceVectors

The Bell source vectors are the following four vectors in C{false,true}×{false,true}\mathbb C^{\{\mathrm{false},\mathrm{true}\}\times\{\mathrm{false},\mathrm{true}\}}C{false,true}×{false,true}, indexed in the displayed order by 0,1,2,30,1,2,30,1,2,3:

v0=e(false,false)+e(true,true)2,v1=e(false,false)−e(true,true)2,v2=e(false,true)+e(true,false)2,v3=e(false,true)−e(true,false)2.\begin{aligned} v_0&=\frac{e_{(\mathrm{false},\mathrm{false})}+e_{(\mathrm{true},\mathrm{true})}}{\sqrt2}, & v_1&=\frac{e_{(\mathrm{false},\mathrm{false})}-e_{(\mathrm{true},\mathrm{true})}}{\sqrt2},\\ v_2&=\frac{e_{(\mathrm{false},\mathrm{true})}+e_{(\mathrm{true},\mathrm{false})}}{\sqrt2}, & v_3&=\frac{e_{(\mathrm{false},\mathrm{true})}-e_{(\mathrm{true},\mathrm{false})}}{\sqrt2}. \end{aligned}v0​v2​​=2​e(false,false)​+e(true,true)​​,=2​e(false,true)​+e(true,false)​​,​v1​v3​​=2​e(false,false)​−e(true,true)​​,=2​e(false,true)​−e(true,false)​​.​

Here e(a,b)e_{(a,b)}e(a,b)​ has coordinate 111 at (a,b)(a,b)(a,b) and 000 at every other coordinate.

QAlgorithms.Papers.Independent.Intakethird100.P11.erasureBasis

This defines four vectors b0,b1,b2,b3b_0,b_1,b_2,b_3b0​,b1​,b2​,b3​ in C{0,1,2,3}×{0,1,2,3}\mathbb C^{\{0,1,2,3\}\times\{0,1,2,3\}}C{0,1,2,3}×{0,1,2,3} by

b0=e00+e332,b1=e03+e302,b2=e11+e222,b3=e21+e122,b_0=\frac{e_{00}+e_{33}}{\sqrt2},\qquad b_1=\frac{e_{03}+e_{30}}{\sqrt2},\qquad b_2=\frac{e_{11}+e_{22}}{\sqrt2},\qquad b_3=\frac{e_{21}+e_{12}}{\sqrt2},b0​=2​e00​+e33​​,b1​=2​e03​+e30​​,b2​=2​e11​+e22​​,b3​=2​e21​+e12​​,

where eabe_{ab}eab​ is the vector whose coordinate at (a,b)(a,b)(a,b) is 111 and whose other coordinates are 000.

QAlgorithms.Papers.Independent.Intakethird100.P11.erasureFidelity

Let n,dA,dBn,d_A,d_Bn,dA​,dB​ be natural numbers and let ccc be a distributed block code for two four-dimensional letter registers. Such a code includes n,dA,dB>0n,d_A,d_B>0n,dA​,dB​>0, separate encoding channels EA,EBE_A,E_BEA​,EB​ to dimensions dA,dBd_A,d_BdA​,dB​, and a decoding channel DDD back to the two length-nnn registers; each channel has a finite Kraus representation X↦∑jKjXKj†X\mapsto\sum_jK_jXK_j^\daggerX↦∑j​Kj​XKj†​ with ∑jKj†Kj=1\sum_jK_j^\dagger K_j=1∑j​Kj†​Kj​=1. Put C=D∘(EA⊗EB)\mathcal C=D\circ(E_A\otimes E_B)C=D∘(EA​⊗EB​), with the encoders acting separately. For a 4×44\times44×4 unitary UUU, set

v(U)=U00(e00+e33)+U10(e03+e30)+U20(e11+e22)+U30(e21+e12)2,v(U)=\frac{ U_{00}(e_{00}+e_{33})+U_{10}(e_{03}+e_{30}) +U_{20}(e_{11}+e_{22})+U_{30}(e_{21}+e_{12}) }{\sqrt2},v(U)=2​U00​(e00​+e33​)+U10​(e03​+e30​)+U20​(e11​+e22​)+U30​(e21​+e12​)​,

where the eabe_{ab}eab​ are coordinate unit vectors. For w=(U0,…,Un−1)w=(U_0,\ldots,U_{n-1})w=(U0​,…,Un−1​), define Vw(a,b)=∏t=0n−1v(Ut)(at,bt)V_w(a,b)=\prod_{t=0}^{n-1}v(U_t)(a_t,b_t)Vw​(a,b)=∏t=0n−1​v(Ut​)(at​,bt​) and Pw=VwVw†P_w=V_wV_w^\daggerPw​=Vw​Vw†​. The defined real number is

∫U(4)nRe⁡tr⁡ ⁣(Pw C(Pw)) dH⊗n(w),\int_{U(4)^n} \operatorname{Re}\operatorname{tr}\!\left(P_w\,\mathcal C(P_w)\right) \,dH^{\otimes n}(w),∫U(4)n​Retr(Pw​C(Pw​))dH⊗n(w),

where HHH is Haar measure on U(4)U(4)U(4), normalized to have total mass 111.

QAlgorithms.Papers.Independent.Intakethird100.P11.erasureVector

For every complex 4×44\times44×4 unitary matrix UUU, this defines a vector in C{0,1,2,3}×{0,1,2,3}\mathbb C^{\{0,1,2,3\}\times\{0,1,2,3\}}C{0,1,2,3}×{0,1,2,3} using its column indexed by 000:

v(U)=12[U00(e00+e33)+U10(e03+e30)+U20(e11+e22)+U30(e21+e12)],v(U)=\frac1{\sqrt2}\Bigl[ U_{00}(e_{00}+e_{33})+ U_{10}(e_{03}+e_{30})+ U_{20}(e_{11}+e_{22})+ U_{30}(e_{21}+e_{12}) \Bigr],v(U)=2​1​[U00​(e00​+e33​)+U10​(e03​+e30​)+U20​(e11​+e22​)+U30​(e21​+e12​)],

where eabe_{ab}eab​ has coordinate 111 at (a,b)(a,b)(a,b) and 000 elsewhere.

QAlgorithms.Papers.Independent.Intakethird100.P11.finiteFidelity

For finite types I,A,BI,A,BI,A,B, with decidable equality on A,BA,BA,B, arbitrary real weights pip_ipi​, arbitrary vectors vi∈CA×Bv_i\in\mathbb C^{A\times B}vi​∈CA×B, natural numbers n,dA,dBn,d_A,d_Bn,dA​,dB​, and a distributed block code ccc, the finite fidelity is

∑w∈In(∏t=0n−1pw(t))Re⁡tr⁡ ⁣(ρw C(ρw)),Vw(a,b)=∏t=0n−1vw(t)(a(t),b(t)),(ρw)x,y=Vw(x)Vw(y)‾.\sum_{w\in I^n} \left(\prod_{t=0}^{n-1}p_{w(t)}\right) \operatorname{Re}\operatorname{tr}\!\left(\rho_w\,\mathcal C(\rho_w)\right), \quad V_w(a,b)=\prod_{t=0}^{n-1}v_{w(t)}(a(t),b(t)), \quad (\rho_w)_{x,y}=V_w(x)\overline{V_w(y)}.w∈In∑​(t=0∏n−1​pw(t)​)Retr(ρw​C(ρw​)),Vw​(a,b)=t=0∏n−1​vw(t)​(a(t),b(t)),(ρw​)x,y​=Vw​(x)Vw​(y)​.

Here C=D∘(EA⊗EB)\mathcal C=\mathcal D\circ(\mathcal E_A\otimes\mathcal E_B)C=D∘(EA​⊗EB​) applies the two encoders separately to the Alice and Bob indices and then the decoder. The code requires n,dA,dB>0n,d_A,d_B>0n,dA​,dB​>0, and each of its three maps has a finite Kraus representation Φ(X)=∑jKjXKj†\Phi(X)=\sum_jK_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. Neither nonnegative or normalized weights nor normalized vectors are required. In particular, if III is empty, the sum is empty because n>0n>0n>0, and its value is zero.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenAlice

For every real number α\alphaα, this defines three vectors in C2\mathbb C^2C2:

a0=e0,a1=α e0+1−α e1,a2=e1,a_0=e_0,\qquad a_1=\sqrt\alpha\,e_0+\sqrt{1-\alpha}\,e_1,\qquad a_2=e_1,a0​=e0​,a1​=α​e0​+1−α​e1​,a2​=e1​,

where e0,e1e_0,e_1e0​,e1​ are coordinate unit vectors and the real square roots are regarded as complex scalars. There is no restriction α∈[0,1]\alpha\in[0,1]α∈[0,1]; the real square-root function used here is 000 on negative inputs, and no normalization condition is part of this definition.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenBob

For every real number β\betaβ, this defines three vectors in C3\mathbb C^3C3:

b0=1−β e0+β e1,b1=e1,b2=1−β e2+β e0,b_0=\sqrt{1-\beta}\,e_0+\sqrt\beta\,e_1,\qquad b_1=e_1,\qquad b_2=\sqrt{1-\beta}\,e_2+\sqrt\beta\,e_0,b0​=1−β​e0​+β​e1​,b1​=e1​,b2​=1−β​e2​+β​e0​,

where e0,e1,e2e_0,e_1,e_2e0​,e1​,e2​ are coordinate unit vectors and the real square roots are regarded as complex scalars. The parameter is not restricted to [0,1][0,1][0,1]; square roots of negative real inputs are 000, and the definition imposes no normalization requirement.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenBobAfter

For every real number β\betaβ, this defines three vectors in C3\mathbb C^3C3:

b~0=1−β e0+β e1,b~1=e1,b~2=e0,\widetilde b_0=\sqrt{1-\beta}\,e_0+\sqrt\beta\,e_1,\qquad \widetilde b_1=e_1,\qquad \widetilde b_2=e_0,b0​=1−β​e0​+β​e1​,b1​=e1​,b2​=e0​,

where e0,e1,e2e_0,e_1,e_2e0​,e1​,e2​ are coordinate unit vectors. The square roots are real and then embedded into C\mathbb CC, with value 000 on negative inputs. No restriction on β\betaβ or normalization requirement is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenComparison

For every pair of real numbers α,β\alpha,\betaα,β, define vectors

(a0,a1,a2)=(e0,α e0+1−α e1,e1)in C2,(a_0,a_1,a_2)= \left(e_0,\sqrt\alpha\,e_0+\sqrt{1-\alpha}\,e_1,e_1\right) \quad\text{in }\mathbb C^2,(a0​,a1​,a2​)=(e0​,α​e0​+1−α​e1​,e1​)in C2, (b0,b1,b2)=(1−β e0+β e1, e1, 1−β e2+β e0)in C3.(b_0,b_1,b_2)= \left(\sqrt{1-\beta}\,e_0+\sqrt\beta\,e_1,\ e_1,\ \sqrt{1-\beta}\,e_2+\sqrt\beta\,e_0\right) \quad\text{in }\mathbb C^3.(b0​,b1​,b2​)=(1−β​e0​+β​e1​, e1​, 1−β​e2​+β​e0​)in C3.

With eke_kek​ denoting coordinate unit vectors, this defines the real number

12[H ⁣(13∑i=02aiai†)+H ⁣(13∑i=02bibi†)+H ⁣(13∑i=02(ai⊗bi)(ai⊗bi)†)].\frac12\left[ H\!\left(\frac13\sum_{i=0}^2a_ia_i^\dagger\right) + H\!\left(\frac13\sum_{i=0}^2b_ib_i^\dagger\right) + H\!\left(\frac13\sum_{i=0}^2 (a_i\otimes b_i)(a_i\otimes b_i)^\dagger\right) \right].21​[H(31​i=0∑2​ai​ai†​)+H(31​i=0∑2​bi​bi†​)+H(31​i=0∑2​(ai​⊗bi​)(ai​⊗bi​)†)].

Here (ai⊗bi)(a,b)=ai(a)bi(b)(a_i\otimes b_i)(a,b)=a_i(a)b_i(b)(ai​⊗bi​)(a,b)=ai​(a)bi​(b), an outer product xx†xx^\daggerxx† has entries xrxs‾x_r\overline{x_s}xr​xs​​, and H(M)=−Re⁡tr⁡(f(M))H(M)=-\operatorname{Re}\operatorname{tr}(f(M))H(M)=−Retr(f(M)), with f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x applied through continuous functional calculus and f(0)=0f(0)=0f(0)=0. Both parameters are unrestricted, and negative real inputs to the square-root function give 000; no trace-one conditions are imposed on these averaged matrices.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenJoint

For every pair of real numbers α,β\alpha,\betaα,β and every i∈{0,1,2}i\in\{0,1,2\}i∈{0,1,2}, this defines a vector ji∈C{0,1}×{0,1,2}j_i\in\mathbb C^{\{0,1\}\times\{0,1,2\}}ji​∈C{0,1}×{0,1,2} by ji(a,b)=ai(a)bi(b)j_i(a,b)=a_i(a)b_i(b)ji​(a,b)=ai​(a)bi​(b), where

(a0,a1,a2)=(e0,α e0+1−α e1,e1)(a_0,a_1,a_2)= \left(e_0,\sqrt\alpha\,e_0+\sqrt{1-\alpha}\,e_1,e_1\right)(a0​,a1​,a2​)=(e0​,α​e0​+1−α​e1​,e1​)

in C2\mathbb C^2C2, and

(b0,b1,b2)=(1−β e0+β e1, e1, 1−β e2+β e0)(b_0,b_1,b_2)= \left(\sqrt{1-\beta}\,e_0+\sqrt\beta\,e_1,\ e_1,\ \sqrt{1-\beta}\,e_2+\sqrt\beta\,e_0\right)(b0​,b1​,b2​)=(1−β​e0​+β​e1​, e1​, 1−β​e2​+β​e0​)

in C3\mathbb C^3C3. Thus ji=ai⊗bij_i=a_i\otimes b_iji​=ai​⊗bi​, with the tensor product specified by the coordinate multiplication above. All eke_kek​ are coordinate unit vectors; both parameters are unrestricted, and real square roots of negative inputs are 000.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenRateA

For every real number α\alphaα, let u=α e0+1−α e1∈C2u=\sqrt\alpha\,e_0+\sqrt{1-\alpha}\,e_1\in\mathbb C^2u=α​e0​+1−α​e1​∈C2. This defines

H ⁣(e0e0†+uu†+e1e1†3),H(M)=−Re⁡tr⁡ ⁣(f(M)),f(x)=xlog⁡2x,H\!\left(\frac{e_0e_0^\dagger+uu^\dagger+e_1e_1^\dagger}{3}\right), \qquad H(M)=-\operatorname{Re}\operatorname{tr}\!\bigl(f(M)\bigr), \quad f(x)=x\log_2x,H(3e0​e0†​+uu†+e1​e1†​​),H(M)=−Retr(f(M)),f(x)=xlog2​x,

where f(M)f(M)f(M) is formed by continuous functional calculus, e0,e1e_0,e_1e0​,e1​ are coordinate unit vectors, and uu†uu^\daggeruu† has entries uaub‾u_a\overline{u_b}ua​ub​​. The scalar convention gives f(0)=0f(0)=0f(0)=0. The parameter α\alphaα is unrestricted, real square roots of negative inputs are 000, and the matrix inside HHH is not required by this definition to have trace 111.

QAlgorithms.Papers.Independent.Intakethird100.P11.hiddenRateB

For every real number β\betaβ, let u=1−β e0+β e1∈C3u=\sqrt{1-\beta}\,e_0+\sqrt\beta\,e_1\in\mathbb C^3u=1−β​e0​+β​e1​∈C3. This defines

H ⁣(uu†+e1e1†+e0e0†3),H(M)=−Re⁡tr⁡ ⁣(f(M)),f(x)=xlog⁡2x,H\!\left(\frac{uu^\dagger+e_1e_1^\dagger+e_0e_0^\dagger}{3}\right), \qquad H(M)=-\operatorname{Re}\operatorname{tr}\!\bigl(f(M)\bigr), \quad f(x)=x\log_2x,H(3uu†+e1​e1†​+e0​e0†​​),H(M)=−Retr(f(M)),f(x)=xlog2​x,

where f(M)f(M)f(M) is defined by continuous functional calculus, e0,e1,e2e_0,e_1,e_2e0​,e1​,e2​ are coordinate unit vectors, and uu†uu^\daggeruu† has entries uaub‾u_a\overline{u_b}ua​ub​​. Here f(0)=0f(0)=0f(0)=0. There is no restriction on β\betaβ; square roots of negative real inputs are 000, and the definition does not require the matrix inside HHH to have trace 111.

QAlgorithms.Papers.Independent.Intakethird100.P11.hybridChi

For arbitrary Boolean-indexed real weights pip_ipi​ and complex coefficients αi,βi\alpha_i,\beta_iαi​,βi​, define wji=pi∣αi∣2w_{ji}=p_i|\alpha_i|^2wji​=pi​∣αi​∣2 for j=falsej=\mathrm{false}j=false and wji=pi∣βi∣2w_{ji}=p_i|\beta_i|^2wji​=pi​∣βi​∣2 for j=truej=\mathrm{true}j=true, and set qj=∑iwjiq_j=\sum_iw_{ji}qj​=∑i​wji​. Let Cj=∣false⟩⟨false∣C_j=|\mathrm{false}\rangle\langle\mathrm{false}|Cj​=∣false⟩⟨false∣ when qj=0q_j=0qj​=0, and otherwise let Cj=qj−1∑iwji∣j⊕i⟩⟨j⊕i∣C_j=q_j^{-1}\sum_iw_{ji}|j\mathbin{\oplus}i\rangle\langle j\mathbin{\oplus}i|Cj​=qj−1​∑i​wji​∣j⊕i⟩⟨j⊕i∣, with ∣b⟩|b\rangle∣b⟩ the coordinate vector at bbb. This defines

χ=S ⁣(∑jqjCj)−∑jqjS(Cj),S(X)=−Re⁡Tr⁡ ⁣[cfc⁡(t↦tlog⁡2t)(X)],\chi=S\!\left(\sum_jq_jC_j\right)-\sum_jq_jS(C_j), \qquad S(X)=-\operatorname{Re}\operatorname{Tr}\!\left[ \operatorname{cfc}\bigl(t\mapsto t\log_2t\bigr)(X)\right],χ=S(j∑​qj​Cj​)−j∑​qj​S(Cj​),S(X)=−ReTr[cfc(t↦tlog2​t)(X)],

where cfc⁡\operatorname{cfc}cfc denotes continuous functional calculus and the real logarithm is the total real logarithm, with value 000 at 000. No distribution, positivity, or normalization hypotheses are required.

QAlgorithms.Papers.Independent.Intakethird100.P11.hybridConditional

For arbitrary real weights pip_ipi​, arbitrary complex coefficients αi,βi\alpha_i,\beta_iαi​,βi​, all indexed by Booleans, and a Boolean jjj, put wji=pi∣αi∣2w_{ji}=p_i|\alpha_i|^2wji​=pi​∣αi​∣2 if j=falsej=\mathrm{false}j=false, put wji=pi∣βi∣2w_{ji}=p_i|\beta_i|^2wji​=pi​∣βi​∣2 if j=truej=\mathrm{true}j=true, and let qj=∑iwjiq_j=\sum_iw_{ji}qj​=∑i​wji​. This defines the Boolean-indexed complex matrix

Cj={∣false⟩⟨false∣,qj=0,qj−1∑iwji∣j⊕i⟩⟨j⊕i∣,qj≠0,C_j= \begin{cases} |\mathrm{false}\rangle\langle\mathrm{false}|,&q_j=0,\\[1mm] \displaystyle q_j^{-1}\sum_iw_{ji}|j\mathbin{\oplus}i\rangle \langle j\mathbin{\oplus}i|,&q_j\ne0, \end{cases}Cj​=⎩⎨⎧​∣false⟩⟨false∣,qj−1​i∑​wji​∣j⊕i⟩⟨j⊕i∣,​qj​=0,qj​=0,​

where ∣b⟩⟨b∣|b\rangle\langle b|∣b⟩⟨b∣ has its only nonzero entry, equal to 111, at (b,b)(b,b)(b,b). The zero case uses this fixed matrix even when qj=0q_j=0qj​=0 results from cancellation of nonzero signed weights. No positivity or normalization assumptions are imposed on the inputs.

QAlgorithms.Papers.Independent.Intakethird100.P11.hybridPiggyback

For arbitrary Boolean-indexed real weights pip_ipi​ and complex coefficient functions αi,βi\alpha_i,\beta_iαi​,βi​, let vi=αi∣false,i⟩+βi∣true,¬i⟩v_i=\alpha_i|\mathrm{false},i\rangle+\beta_i|\mathrm{true},\neg i\ranglevi​=αi​∣false,i⟩+βi​∣true,¬i⟩, where the kets are coordinate vectors, and let ρ=∑ipi∣vi⟩⟨vi∣\rho=\sum_i p_i|v_i\rangle\langle v_i|ρ=∑i​pi​∣vi​⟩⟨vi​∣, so ρab=∑ipivi(a)vi(b)‾\rho_{ab}= \sum_i p_i v_i(a)\overline{v_i(b)}ρab​=∑i​pi​vi​(a)vi​(b)​. This defines the real number

B=1+S(Tr⁡firstρ)−S(ρ),(Tr⁡firstρ)b,b′=∑a∈{false,true}ρ(a,b),(a,b′),B=1+S(\operatorname{Tr}_{\mathrm{first}}\rho)-S(\rho), \qquad (\operatorname{Tr}_{\mathrm{first}}\rho)_{b,b'} =\sum_{a\in\{\mathrm{false},\mathrm{true}\}}\rho_{(a,b),(a,b')},B=1+S(Trfirst​ρ)−S(ρ),(Trfirst​ρ)b,b′​=a∈{false,true}∑​ρ(a,b),(a,b′)​,

where S(X)=−Re⁡Tr⁡[cfc⁡(t↦tlog⁡2t)(X)]S(X)=-\operatorname{Re}\operatorname{Tr}[\operatorname{cfc}(t\mapsto t\log_2t)(X)]S(X)=−ReTr[cfc(t↦tlog2​t)(X)], using continuous functional calculus and the total real logarithm with ln⁡0=0\ln0=0ln0=0. The inputs need not be probabilities or normalized coefficients, and ρ\rhoρ is not assumed to be a density matrix.

QAlgorithms.Papers.Independent.Intakethird100.P11.hybridQ

For arbitrary functions p:{false,true}→Rp:\{\mathrm{false},\mathrm{true}\}\to\mathbb Rp:{false,true}→R and α,β:{false,true}→C\alpha,\beta:\{\mathrm{false},\mathrm{true}\}\to\mathbb Cα,β:{false,true}→C, and any Boolean jjj, this defines the real number

qj={∑i∈{false,true}pi∣αi∣2,j=false,∑i∈{false,true}pi∣βi∣2,j=true.q_j= \begin{cases} \displaystyle\sum_{i\in\{\mathrm{false},\mathrm{true}\}}p_i|\alpha_i|^2, &j=\mathrm{false},\\[2mm] \displaystyle\sum_{i\in\{\mathrm{false},\mathrm{true}\}}p_i|\beta_i|^2, &j=\mathrm{true}. \end{cases}qj​=⎩⎨⎧​i∈{false,true}∑​pi​∣αi​∣2,i∈{false,true}∑​pi​∣βi​∣2,​j=false,j=true.​

The weights pip_ipi​ need not be nonnegative or sum to 111, and the coefficients need not be normalized, so these numbers are not required to form a probability distribution.

QAlgorithms.Papers.Independent.Intakethird100.P11.hybridRateA

For arbitrary Boolean-indexed real weights pip_ipi​ and complex coefficients αi,βi\alpha_i,\beta_iαi​,βi​, set vi=αi∣false,i⟩+βi∣true,¬i⟩v_i=\alpha_i|\mathrm{false},i\rangle+\beta_i|\mathrm{true},\neg i\ranglevi​=αi​∣false,i⟩+βi​∣true,¬i⟩ and ρ=∑ipi∣vi⟩⟨vi∣\rho=\sum_i p_i|v_i\rangle\langle v_i|ρ=∑i​pi​∣vi​⟩⟨vi​∣, with kets denoting coordinate vectors. Put wji=pi∣αi∣2w_{ji}=p_i|\alpha_i|^2wji​=pi​∣αi​∣2 for j=falsej=\mathrm{false}j=false, wji=pi∣βi∣2w_{ji}=p_i|\beta_i|^2wji​=pi​∣βi​∣2 for j=truej=\mathrm{true}j=true, qj=∑iwjiq_j=\sum_iw_{ji}qj​=∑i​wji​, and

Cj={∣false⟩⟨false∣,qj=0,qj−1∑iwji∣j⊕i⟩⟨j⊕i∣,qj≠0.C_j= \begin{cases} |\mathrm{false}\rangle\langle\mathrm{false}|,&q_j=0,\\ q_j^{-1}\sum_iw_{ji}|j\mathbin{\oplus}i\rangle\langle j\mathbin{\oplus}i|,&q_j\ne0. \end{cases}Cj​={∣false⟩⟨false∣,qj−1​∑i​wji​∣j⊕i⟩⟨j⊕i∣,​qj​=0,qj​=0.​

With S(X)=−Re⁡Tr⁡[cfc⁡(t↦tlog⁡2t)(X)]S(X)=-\operatorname{Re}\operatorname{Tr}[\operatorname{cfc}(t\mapsto t\log_2t)(X)]S(X)=−ReTr[cfc(t↦tlog2​t)(X)], define H(q)=−∑jqjlog⁡2qjH(q)=-\sum_jq_j\log_2q_jH(q)=−∑j​qj​log2​qj​, χ=S(∑jqjCj)−∑jqjS(Cj)\chi=S(\sum_jq_jC_j)-\sum_jq_jS(C_j)χ=S(∑j​qj​Cj​)−∑j​qj​S(Cj​), and B=1+S(Tr⁡firstρ)−S(ρ)B=1+S(\operatorname{Tr}_{\mathrm{first}}\rho)-S(\rho)B=1+S(Trfirst​ρ)−S(ρ), where (Tr⁡firstρ)b,b′=∑aρ(a,b),(a,b′)(\operatorname{Tr}_{\mathrm{first}}\rho)_{b,b'}=\sum_a\rho_{(a,b),(a,b')}(Trfirst​ρ)b,b′​=∑a​ρ(a,b),(a,b′)​. The definition is

RA=H(q)−χH(q)−χ+B.R_A=\frac{H(q)-\chi}{H(q)-\chi+B}.RA​=H(q)−χ+BH(q)−χ​.

No probability or normalization assumptions are imposed. Logarithms use the total real logarithm, including ln⁡0=0\ln0=0ln0=0, and division by zero has value 000, so RA=0R_A=0RA​=0 whenever its denominator is zero.

QAlgorithms.Papers.Independent.Intakethird100.P11.hybridVector

For arbitrary functions α,β:{false,true}→C\alpha,\beta:\{\mathrm{false},\mathrm{true}\}\to\mathbb Cα,β:{false,true}→C and a Boolean label iii, this defines a vector viv_ivi​ in the complex Euclidean space indexed by pairs of Booleans by

vi=αi ∣false,i⟩+βi ∣true,¬i⟩.v_i=\alpha_i\,|\mathrm{false},i\rangle+ \beta_i\,|\mathrm{true},\neg i\rangle.vi​=αi​∣false,i⟩+βi​∣true,¬i⟩.

Here each displayed coordinate vector has entry 111 at its indicated pair and 000 elsewhere; equivalently, the amplitude at (a,b)(a,b)(a,b) is zero unless b=a⊕ib=a\mathbin{\oplus}ib=a⊕i, and in that case it is αi\alpha_iαi​ when a=falsea=\mathrm{false}a=false and βi\beta_iβi​ when a=truea=\mathrm{true}a=true. No normalization or other restriction is imposed on either coefficient function.

QAlgorithms.Papers.Independent.Intakethird100.P11.letterFidelity

For an arbitrary label type III, finite types A,BA,BA,B with decidable equality, any family of vectors vi∈CA×Bv_i\in\mathbb C^{A\times B}vi​∈CA×B, natural numbers n,dA,dBn,d_A,d_Bn,dA​,dB​, a distributed block code ccc, and a word w∈Inw\in I^nw∈In, define Vw(a,b)=∏t=0n−1vw(t)(a(t),b(t))V_w(a,b)=\prod_{t=0}^{n-1}v_{w(t)}(a(t),b(t))Vw​(a,b)=∏t=0n−1​vw(t)​(a(t),b(t)) and the matrix ρw\rho_wρw​ by (ρw)x,y=Vw(x)Vw(y)‾(\rho_w)_{x,y}=V_w(x)\overline{V_w(y)}(ρw​)x,y​=Vw​(x)Vw​(y)​. The letter fidelity is the real number

Re⁡tr⁡ ⁣(ρw C(ρw)),C=D∘(EA⊗EB),\operatorname{Re}\operatorname{tr}\!\left(\rho_w\,\mathcal C(\rho_w)\right), \qquad \mathcal C=\mathcal D\circ(\mathcal E_A\otimes\mathcal E_B),Retr(ρw​C(ρw​)),C=D∘(EA​⊗EB​),

where ccc has positive length nnn, positive compressed dimensions dA,dBd_A,d_BdA​,dB​, and encoders EA,EB\mathcal E_A,\mathcal E_BEA​,EB​ and decoder D\mathcal DD, each admitting a finite representation Φ(X)=∑jKjXKj†\Phi(X)=\sum_jK_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. The tensor action applies Alice’s encoder to the Alice matrix indices and Bob’s encoder to the Bob matrix indices. The vectors viv_ivi​ are not assumed normalized, so the definition does not itself require ρw\rho_wρw​ to be a density matrix or the resulting number to lie in [0,1][0,1][0,1].

QAlgorithms.Papers.Independent.Intakethird100.P11.sourceMean

For a finite label type III, an arbitrary type XXX, arbitrary real weights pip_ipi​, and arbitrary vectors vi∈CXv_i\in\mathbb C^Xvi​∈CX, the source mean is the complex matrix indexed by XXX whose entries are

(sourceMean⁡(p,v))x,y=∑i∈Ipi vi(x)vi(y)‾.(\operatorname{sourceMean}(p,v))_{x,y} = \sum_{i\in I}p_i\,v_i(x)\overline{v_i(y)}.(sourceMean(p,v))x,y​=i∈I∑​pi​vi​(x)vi​(y)​.

The definition obtains this sum by enumerating III with the integers from 000 to ∣I∣−1|I|-1∣I∣−1 and summing the weighted matrices vivi†v_iv_i^\daggervi​vi†​. No finiteness assumption on XXX, probability conditions on the weights, or normalization conditions on the vectors are imposed. If III is empty, this is the zero matrix.

QAlgorithms.Papers.Independent.Intakethird100.P11.sourceWordVector

For arbitrary types I,A,BI,A,BI,A,B, any family vi∈CA×Bv_i\in\mathbb C^{A\times B}vi​∈CA×B of Euclidean-space vectors indexed by i∈Ii\in Ii∈I, any natural number nnn, and any word w:{0,…,n−1}→Iw:\{0,\ldots,n-1\}\to Iw:{0,…,n−1}→I, the source word vector is the vector indexed by An×BnA^n\times B^nAn×Bn whose coordinates are

Vw(a,b)=∏t=0n−1vw(t)(a(t),b(t)).V_w(a,b)=\prod_{t=0}^{n-1}v_{w(t)}(a(t),b(t)).Vw​(a,b)=t=0∏n−1​vw(t)​(a(t),b(t)).

No finiteness, decidable-equality, or normalization assumptions are imposed here. When n=0n=0n=0, the empty product is 111, so the vector has value 111 at the unique pair of empty functions.

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