Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Quantum rate distortion coding with auxiliary resources

Definition
QAlgorithms_Papers_Independent_Intakethird100_P27

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

prove2quantum

The objects the formalized statements of Quantum rate distortion coding with auxiliary resources (arXiv:1212.5316) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • DistortionSource: PDF pp.3--4: a normalized source purification, with reference isomorphic to A, and an arbitrary positive distortion observable on the reference and B. The source density is the actual A marginal, traceLeft (pureDensity psi).
  • IidRegister: PDF p.4: n copies of a register. Only positive n enter operational codes.
  • testOutput: PDF p.4: the single-letter output includes the untouched purifying reference.
  • observableDistortion: PDF p.4: the expectation of delta; no upper bound delta <= I is imposed.
  • letterMarginal: PDF p.4: one-letter reference/reconstruction marginal, tracing all other letters jointly. The kept and forgotten coordinates form a bijective partition.
  • averageDistortion: PDF p.4: average distortion from the actual n-letter reference/output state. The total formula is zero at n=0; every code and block minimum uses n>=1.
  • LossyAssistance: PDF pp.4,6--7: the three distinct allowed auxiliary-resource models.
  • LossyCode: PDF pp.3--4,6--7: actual CPTP encoding and decoding with one quantum message W, one dephased classical message M, and an independent shared resource. The latter two registers are constrained to dimension one when unavailable. Rectangular CPTP maps include arbitrary local scratch systems and discarding.
  • LossyCode.Allowed: PDF pp.4,6: availability constraints, independent of the source and distortion. Unlimited forward classical communication supplies no prior entanglement.
  • LossyCode.initial: PDF pp.4,7: iid purification tensor an independent shared density matrix. The outer register is R^n; Alice owns A^n and the first shared half.
  • LossyCode.output: PDF pp.4,7: actual encoding, forward transmission of MW, and decoding with Bob's resource half; all off-diagonal R^n blocks are preserved.
  • LossyAchievable: PDF p.4, (8): vanishing rate/distortion slack, for all sufficiently large positive blocklengths. The code is selected before any reference indices. No entropic bound occurs in this operational definition.
  • rateDistortion: PDF p.4: infimum over operationally achievable qubit rates. ENNReal keeps infeasibility at infinity and makes zero-weight time-sharing endpoints exact.
  • sourceFidelity: PDF p.4, (7): source-dependent entanglement fidelity, for an arbitrary unit purification rather than only a maximally mixed qubit source.
  • iidChannelOutput: PDF p.6, (16): arbitrary block channel applied to the actual iid purification.
  • formationBlockValue: PDF p.6, (14): normalized block formation infimum over exactly the feasible CPTP maps. For positive k, empty feasible sets have value infinity, never the real sInf default. At k=0 the total expression is zero and is never used.
  • isotropicEntropy: PDF p.8: Bell-error probability vector and its Shannon entropy in bits. Zero probabilities contribute zero, including at D=0.
  • purificationValues: PDF pp.3,10, (3),(22): costs of all finite bipartite purifications, with arbitrary positive ancillary dimensions. The actual AB marginal is fixed.
  • purificationCost: PDF p.3, (3): entanglement of purification, taking a genuine infimum over all finite extensions. The density matrices used in the bound are normalized.
  • purificationBlockValue: PDF pp.4,17, (9),(50): normalized block purification objective, over all feasible block channels rather than just tensor powers of a one-letter channel.
  • regularizedPurification: PDF p.17, (50): the regularized expression; positive blocklengths only. The convexity theorem below asserts actual convergence as well as convexity.
  • isotropicStrategyCost: PDF p.10, (22) and footnote 2: use the accepted isotropic state with local dimension 2 and Bell overlap 1-D, so its eigenvalues are 1-D,D/3,D/3,D/3. Use E_p only up to 3/4; a constant reconstruction costs zero thereafter. Mixture distortions are restricted to [0,1].
  • isotropicConvexUpper: PDF p.10, Remark 5 and footnote 2: all finite time-sharing mixtures whose average distortion is at most D. This is a lower convex envelope, not a sampled plot, and no numerical assertion about evaluating E_p is made.
  • outputWithSide: PDF p.14: apply N to A in R,(A,B), preserving all R and B coherences.
  • conditionalOutputInformation: PDF pp.3,14: I(R;C|B), from the four actual channel-output marginals.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_Ensemble
import Definitions.Def_QAlgorithms_Defs_Measurement
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P19
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P08
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P20
import Definitions.Def_QAlgorithms_Papers_Independent_Bridge50_P37

/-!
Statements from Wilde--Datta--Hsieh--Winter, arXiv:1212.5316v3.
All new definitions are confined to this paper. Registers are finite nonempty
basis types in the results. Entropies and message costs are in bits. Operational
rates are independent infima over actual CPTP codes, taking value infinity below
feasible distortion. Compilation is pending the mechanical batch handoff.
The entropy-duality and redistribution nodes reuse accepted results in items.json.
Seven other nodes are explicitly deferred in HARD.md and GAPS.json.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P27

open Filter
open scoped Matrix ComplexConjugate ComplexOrder ENNReal

attribute [local instance] Classical.propDecidable Classical.decEq

open QAlgorithms.Papers.Independent.Dependency50.P19
  (first marginalRB mutualInformation bipartiteBlock epr)
open QAlgorithms.Papers.Independent.Information100.P08 (formation)
open QAlgorithms.Papers.Independent.Gap48.P20 (redistributionMarginal)
open QAlgorithms.Papers.Independent.Bridge50.P37 (isotropicState normalizedEntangledProjection)

/-- PDF pp.3--4: a normalized source purification, with reference isomorphic to
A, and an arbitrary positive distortion observable on the reference and B.
The source density is the actual A marginal, traceLeft (pureDensity psi). -/
structure DistortionSource (A B : Type) [Fintype A] [Fintype B] where
  psi : EuclideanSpace ℂ (A × A)
  normalized : ‖psi‖ = 1
  delta : Matrix (A × B) (A × B) ℂ
  positive : delta.PosSemidef

/-- PDF p.4: n copies of a register. Only positive n enter operational codes. -/
abbrev IidRegister (A : Type) (n : ℕ) := Fin n → A

/-- PDF p.4: the single-letter output includes the untouched purifying reference. -/
def testOutput {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (N : Superop A B) : Matrix (A × B) (A × B) ℂ :=
  idTensorSuperop N (pureDensity S.psi)

/-- PDF p.4: the expectation of delta; no upper bound delta <= I is imposed. -/
def observableDistortion {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (N : Superop A B) : ℝ :=
  effectProb S.delta (testOutput S N)

/-- PDF p.4: one-letter reference/reconstruction marginal, tracing all other
letters jointly. The kept and forgotten coordinates form a bijective partition. -/
def letterMarginal {A B : Type} {n : ℕ}
    [Fintype A] [Fintype B] (X : Matrix (IidRegister A n × IidRegister B n)
      (IidRegister A n × IidRegister B n) ℂ) (i : Fin n) : Matrix (A × B) (A × B) ℂ :=
  redistributionMarginal X (fun p => (p.1 i, p.2 i))
    (fun p => fun j : {j : Fin n // j ≠ i} => (p.1 j.1, p.2 j.1))

/-- PDF p.4: average distortion from the actual n-letter reference/output state.
The total formula is zero at n=0; every code and block minimum uses n>=1. -/
def averageDistortion {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) {n : ℕ}
    (X : Matrix (IidRegister A n × IidRegister B n) (IidRegister A n × IidRegister B n) ℂ) : ℝ :=
  (∑ i : Fin n, effectProb S.delta (letterMarginal X i)) / (n : ℝ)

/-- PDF pp.4,6--7: the three distinct allowed auxiliary-resource models. -/
inductive LossyAssistance where
  | unassisted
  | entanglement
  | classical
  deriving DecidableEq

/-- PDF pp.3--4,6--7: actual CPTP encoding and decoding with one quantum message
W, one dephased classical message M, and an independent shared resource.
The latter two registers are constrained to dimension one when unavailable.
Rectangular CPTP maps include arbitrary local scratch systems and discarding. -/
structure LossyCode (A B : Type) [Fintype A] [Fintype B] (n : ℕ) where
  quantumDim : ℕ
  classicalDim : ℕ
  sharedDim : ℕ
  quantum_positive : 0 < quantumDim
  classical_positive : 0 < classicalDim
  shared_positive : 0 < sharedDim
  resource : Matrix (Fin sharedDim × Fin sharedDim) (Fin sharedDim × Fin sharedDim) ℂ
  resource_density : IsDensityMatrix resource
  encoder : Superop (IidRegister A n × Fin sharedDim) (Fin classicalDim × Fin quantumDim)
  decoder : Superop ((Fin classicalDim × Fin quantumDim) × Fin sharedDim) (IidRegister B n)
  encoder_channel : IsChannel encoder
  decoder_channel : IsChannel decoder
  classical_message : ∀ X m m' w w', m ≠ m' → encoder X (m, w) (m', w') = 0

/-- PDF pp.4,6: availability constraints, independent of the source and distortion.
Unlimited forward classical communication supplies no prior entanglement. -/
def LossyCode.Allowed {A B : Type} [Fintype A] [Fintype B] {n : ℕ}
    (P : LossyCode A B n) : LossyAssistance → Prop
  | .unassisted => P.classicalDim = 1 ∧ P.sharedDim = 1
  | .entanglement => P.classicalDim = 1
  | .classical => P.sharedDim = 1

/-- PDF pp.4,7: iid purification tensor an independent shared density matrix.
The outer register is R^n; Alice owns A^n and the first shared half. -/
def LossyCode.initial {A B : Type} [Fintype A] [Fintype B] {n : ℕ}
    (P : LossyCode A B n) (S : DistortionSource A B) :
    Matrix (IidRegister A n × ((IidRegister A n × Fin P.sharedDim) × Fin P.sharedDim))
      (IidRegister A n × ((IidRegister A n × Fin P.sharedDim) × Fin P.sharedDim)) ℂ :=
  fun p q =>
    (∏ i, pureDensity S.psi (p.1 i, p.2.1.1 i) (q.1 i, q.2.1.1 i)) *
      P.resource (p.2.1.2, p.2.2) (q.2.1.2, q.2.2)

/-- PDF pp.4,7: actual encoding, forward transmission of MW, and decoding with
Bob's resource half; all off-diagonal R^n blocks are preserved. -/
def LossyCode.output {A B : Type} [Fintype A] [Fintype B] {n : ℕ}
    (P : LossyCode A B n) (S : DistortionSource A B) :
    Matrix (IidRegister A n × IidRegister B n) (IidRegister A n × IidRegister B n) ℂ :=
  idTensorSuperop (P.decoder ∘ first P.encoder) (P.initial S)

/-- PDF p.4, (8): vanishing rate/distortion slack, for all sufficiently large
positive blocklengths. The code is selected before any reference indices.
No entropic bound occurs in this operational definition. -/
def LossyAchievable {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (a : LossyAssistance) (rate D : ℝ) : Prop :=
  0 ≤ rate ∧ 0 ≤ D ∧ ∀ eps > (0 : ℝ), ∀ slack > (0 : ℝ), ∃ n0 : ℕ,
    ∀ n ≥ n0, 1 ≤ n → ∃ P : LossyCode A B n, P.Allowed a ∧
      Real.logb 2 (P.quantumDim : ℝ) ≤ (n : ℝ) * (rate + slack) ∧
      averageDistortion S (P.output S) ≤ D + eps

/-- PDF p.4: infimum over operationally achievable qubit rates. ENNReal keeps
infeasibility at infinity and makes zero-weight time-sharing endpoints exact. -/
def rateDistortion {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (a : LossyAssistance) (D : ℝ) : ℝ≥0∞ :=
  ⨅ (r : ℝ) (_ : LossyAchievable S a r D), ENNReal.ofReal r

/-- PDF p.4, (7): source-dependent entanglement fidelity, for an arbitrary unit
purification rather than only a maximally mixed qubit source. -/
def sourceFidelity {R A : Type} [Fintype R] [Fintype A]
    (v : EuclideanSpace ℂ (R × A)) (N : Superop A A) : ℝ :=
  (∑ i, ∑ j, conj (v i) * idTensorSuperop N (pureDensity v) i j * v j).re

/-- PDF p.6, (16): arbitrary block channel applied to the actual iid purification. -/
def iidChannelOutput {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (k : ℕ) (N : Superop (IidRegister A k) (IidRegister B k)) :
    Matrix (IidRegister A k × IidRegister B k) (IidRegister A k × IidRegister B k) ℂ :=
  idTensorSuperop N (bipartiteBlock k (pureDensity S.psi))

/-- PDF p.6, (14): normalized block formation infimum over exactly the feasible
CPTP maps. For positive k, empty feasible sets have value infinity, never the
real sInf default. At k=0 the total expression is zero and is never used. -/
def formationBlockValue {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (D : ℝ) (k : ℕ) : ℝ≥0∞ :=
  ENNReal.ofReal ((k : ℝ)⁻¹) *
    ⨅ (N : Superop (IidRegister A k) (IidRegister B k))
      (_ : IsChannel N ∧ averageDistortion S (iidChannelOutput S k N) ≤ D),
        ENNReal.ofReal (formation (iidChannelOutput S k N))

/-- PDF p.8: Bell-error probability vector and its Shannon entropy in bits.
Zero probabilities contribute zero, including at D=0. -/
def isotropicEntropy (D : ℝ) : ℝ :=
  let p : Fin 4 → ℝ := ![1 - D, D / 3, D / 3, D / 3];
  -(∑ i, p i * Real.logb 2 (p i))

/-- PDF pp.3,10, (3),(22): costs of all finite bipartite purifications, with
arbitrary positive ancillary dimensions. The actual AB marginal is fixed. -/
def purificationValues {A B : Type} [Fintype A] [Fintype B]
    (rho : Matrix (A × B) (A × B) ℂ) : Set ℝ≥0∞ :=
  {r | ∃ (u v : ℕ) (_ : 0 < u) (_ : 0 < v)
      (psi : EuclideanSpace ℂ ((A × Fin u) × (B × Fin v))),
    ‖psi‖ = 1 ∧
    redistributionMarginal (pureDensity psi) (fun p => (p.1.1, p.2.1))
      (fun p => (p.1.2, p.2.2)) = rho ∧
    r = ENNReal.ofReal (vnEntropy (traceRight (pureDensity psi)))}

/-- PDF p.3, (3): entanglement of purification, taking a genuine infimum over
all finite extensions. The density matrices used in the bound are normalized. -/
def purificationCost {A B : Type} [Fintype A] [Fintype B]
    (rho : Matrix (A × B) (A × B) ℂ) : ℝ≥0∞ :=
  sInf (purificationValues rho)

/-- PDF pp.4,17, (9),(50): normalized block purification objective, over all
feasible block channels rather than just tensor powers of a one-letter channel. -/
def purificationBlockValue {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (D : ℝ) (k : ℕ) : ℝ≥0∞ :=
  ENNReal.ofReal ((k : ℝ)⁻¹) *
    ⨅ (N : Superop (IidRegister A k) (IidRegister B k))
      (_ : IsChannel N ∧ averageDistortion S (iidChannelOutput S k N) ≤ D),
        purificationCost (iidChannelOutput S k N)

/-- PDF p.17, (50): the regularized expression; positive blocklengths only.
The convexity theorem below asserts actual convergence as well as convexity. -/
def regularizedPurification {A B : Type} [Fintype A] [Fintype B]
    (S : DistortionSource A B) (D : ℝ) : ℝ≥0∞ :=
  liminf (fun k => purificationBlockValue S D (k + 1)) atTop

/-- PDF p.10, (22) and footnote 2: use the accepted isotropic state with local
dimension 2 and Bell overlap 1-D, so its eigenvalues are 1-D,D/3,D/3,D/3.
Use E_p only up to 3/4; a constant reconstruction costs zero thereafter.
Mixture distortions are restricted to [0,1]. -/
def isotropicStrategyCost (D : ℝ) : ℝ≥0∞ :=
  if D ≤ 3 / 4 then purificationCost (isotropicState 2 (1 - D)) else 0

/-- PDF p.10, Remark 5 and footnote 2: all finite time-sharing mixtures whose
average distortion is at most D. This is a lower convex envelope, not a sampled
plot, and no numerical assertion about evaluating E_p is made. -/
def isotropicConvexUpper (D : ℝ) : ℝ≥0∞ :=
  sInf {r | ∃ (m : ℕ) (_ : 0 < m) (p d : Fin m → ℝ),
    (∀ i, 0 ≤ p i) ∧ (∑ i, p i) = 1 ∧
    (∀ i, 0 ≤ d i ∧ d i ≤ 1) ∧ (∑ i, p i * d i) ≤ D ∧
    r = ∑ i, ENNReal.ofReal (p i) * isotropicStrategyCost (d i)}

/-- PDF p.14: apply N to A in R,(A,B), preserving all R and B coherences. -/
def outputWithSide {R A B C : Type} (N : Superop A C)
    (phi : EuclideanSpace ℂ (R × (A × B))) :
    Matrix (R × (C × B)) (R × (C × B)) ℂ :=
  idTensorSuperop (first N) (pureDensity phi)

/-- PDF pp.3,14: I(R;C|B), from the four actual channel-output marginals. -/
def conditionalOutputInformation {R A B C : Type}
    [Fintype R] [Fintype A] [Fintype B] [Fintype C]
    (N : Superop A C) (phi : EuclideanSpace ℂ (R × (A × B))) : ℝ :=
  let X := outputWithSide N phi
  vnEntropy (marginalRB X) + vnEntropy (traceLeft X) -
    vnEntropy (traceLeft (traceLeft X)) - vnEntropy X

/-- Small sanity example: the zero-length helper has zero average distortion.
Every operational theorem excludes this blocklength. -/

example {A B : Type} [Fintype A] [Fintype B] (S : DistortionSource A B)
    (X : Matrix (IidRegister A 0 × IidRegister B 0) (IidRegister A 0 × IidRegister B 0) ℂ) :
    averageDistortion S X = 0 := by
  simp [averageDistortion]

example {A : Type} [Fintype A] (S : DistortionSource A A) :
    testOutput S (fun X => X) = pureDensity S.psi := by
  ext p q
  rfl

example {A B : Type} [Fintype A] [Fintype B] {n : ℕ}
    (P : LossyCode A B n) : P.Allowed .classical ↔ P.sharedDim = 1 := Iff.rfl

example {A B : Type} [Fintype A] [Fintype B] {n : ℕ}
    (P : LossyCode A B n) : P.Allowed .unassisted → P.sharedDim = 1 := by
  exact And.right

example : isotropicEntropy 0 = 0 := by
  norm_num [isotropicEntropy, Fin.sum_univ_four]
  change -((0 : ℝ) * Real.logb 2 0) + -(0 * Real.logb 2 0) = 0
  simp

example : isotropicStrategyCost 1 = 0 := by
  norm_num [isotropicStrategyCost]

example : isotropicState 2 1 = normalizedEntangledProjection (Fin 2) := by
  simp [isotropicState]

example {R A B : Type} [Fintype A] (phi : EuclideanSpace ℂ (R × (A × B))) :
    outputWithSide (fun X => X) phi = pureDensity phi := by
  ext p q
  rfl

end QAlgorithms.Papers.Independent.Intakethird100.P27
Source
Quantum rate distortion coding with auxiliary resources, arXiv:1212.5316v3, https://arxiv.org/abs/1212.5316v3
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P27.DistortionSource

For finite types AAA and BBB, a distortion source consists of a complex vector ψ\psiψ indexed by A×AA\times AA×A, satisfying ∥ψ∥=1\|\psi\|=1∥ψ∥=1, and a positive semidefinite complex matrix Δ\DeltaΔ indexed by A×BA\times BA×B. There is no normalization or upper bound imposed on Δ\DeltaΔ. The normalization of ψ\psiψ excludes empty AAA, but the definition does not require BBB to be nonempty.

QAlgorithms.Papers.Independent.Intakethird100.P27.IidRegister

For any type AAA and natural number nnn, the register AnA^nAn is the type of functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to AAA. No finiteness assumption on AAA is required. When n=0n=0n=0, this register has exactly one element, the empty function, including when AAA is empty.

QAlgorithms.Papers.Independent.Intakethird100.P27.LossyAchievable

For finite types A,BA,BA,B, a source consisting of a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B, an assistance label aaa among unassisted\mathrm{unassisted}unassisted, entanglement\mathrm{entanglement}entanglement, and classical\mathrm{classical}classical, and real numbers r,Dr,Dr,D, lossy achievability means r≥0r\ge0r≥0, D≥0D\ge0D≥0, and the following: for every ε>0\varepsilon>0ε>0 and every s>0s>0s>0, there exists a natural number n0n_0n0​ such that, for every natural number n≥n0n\ge n_0n≥n0​ with n≥1n\ge1n≥1, there is a code with positive dimensions q,c,hq,c,hq,c,h, a positive semidefinite trace-one resource τ\tauτ on [h]×[h][h]\times[h][h]×[h], a finite Kraus trace-preserving encoder from matrices on An×[h]A^n\times[h]An×[h] to matrices on [c]×[q][c]\times[q][c]×[q], and a finite Kraus trace-preserving decoder from matrices on ([c]×[q])×[h]([c]\times[q])\times[h]([c]×[q])×[h] to matrices on BnB^nBn. The encoder output must have zero entries between distinct classical indices for every input matrix. The assistance restriction is c=h=1c=h=1c=h=1 for unassisted\mathrm{unassisted}unassisted, c=1c=1c=1 for entanglement\mathrm{entanglement}entanglement, and h=1h=1h=1 for classical\mathrm{classical}classical. Starting with the matrix whose entries are ∏iψ(ri,xi)ψ(ri′,xi′)‾ τ(u,v),(u′,v′)\prod_i\psi(r_i,x_i)\overline{\psi(r'_i,x'_i)}\,\tau_{(u,v),(u',v')}∏i​ψ(ri​,xi​)ψ(ri′​,xi′​)​τ(u,v),(u′,v′)​, apply the encoder to x,ux,ux,u, retain r,vr,vr,v, and then apply the decoder to the encoder output and vvv; call the resulting matrix on An×BnA^n\times B^nAn×Bn XXX. The code must satisfy

log⁡2q≤n(r+s),1n∑i=0n−1Re⁡tr⁡(ΔXi)≤D+ε,\log_2 q\le n(r+s), \qquad \frac1n\sum_{i=0}^{n-1}\operatorname{Re}\operatorname{tr}(\Delta X_i)\le D+\varepsilon,log2​q≤n(r+s),n1​i=0∑n−1​Retr(ΔXi​)≤D+ε,

where XiX_iXi​ is obtained by fixing the two retained indices at position iii and summing entries of XXX over matching row and column indices at every other position. Only the quantum dimension appears in the rate bound; the classical and shared dimensions have no additional size bound.

QAlgorithms.Papers.Independent.Intakethird100.P27.LossyAssistance

The assistance type has exactly three distinct values, named unassisted\mathrm{unassisted}unassisted, entanglement\mathrm{entanglement}entanglement, and classical\mathrm{classical}classical, and equality between its values is decidable. This declaration supplies only these labels and does not itself attach resource restrictions to them.

QAlgorithms.Papers.Independent.Intakethird100.P27.LossyCode.Allowed

For finite types A,BA,BA,B, any natural number nnn, and a lossy code PPP with positive quantum, classical, and shared dimensions q,c,hq,c,hq,c,h, a positive semidefinite trace-one shared resource on two hhh-dimensional registers, and finite Kraus trace-preserving encoder and decoder maps of the dimensions specified by the code, with the encoder output block diagonal in its classical index, the allowed-assistance predicate is defined as follows: unassisted\mathrm{unassisted}unassisted requires c=1c=1c=1 and h=1h=1h=1; entanglement\mathrm{entanglement}entanglement requires c=1c=1c=1; and classical\mathrm{classical}classical requires h=1h=1h=1. Beyond the code’s existing requirements, these cases impose no further restriction on qqq, on the unrestricted dimension in the latter two cases, or on the shared resource.

QAlgorithms.Papers.Independent.Intakethird100.P27.LossyCode.initial

For finite types A,BA,BA,B, any natural number nnn, a lossy code PPP with positive dimensions q,c,hq,c,hq,c,h, shared density matrix τ\tauτ on two hhh-dimensional registers, finite Kraus trace-preserving encoder and decoder, and encoder output block diagonal in its classical index, and a source consisting of a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B, the initial matrix is indexed by An×((An×{0,…,h−1})×{0,…,h−1})A^n\times((A^n\times\{0,\ldots,h-1\})\times\{0,\ldots,h-1\})An×((An×{0,…,h−1})×{0,…,h−1}) and has entries

Ω(r,((x,u),v)),(r′,((x′,u′),v′))=(∏i=0n−1ψ(ri,xi)ψ(ri′,xi′)‾)τ(u,v),(u′,v′).\Omega_{(r,((x,u),v)),(r',((x',u'),v'))} = \left(\prod_{i=0}^{n-1}\psi(r_i,x_i)\overline{\psi(r'_i,x'_i)}\right) \tau_{(u,v),(u',v')}.Ω(r,((x,u),v)),(r′,((x′,u′),v′))​=(i=0∏n−1​ψ(ri​,xi​)ψ(ri′​,xi′​)​)τ(u,v),(u′,v′)​.

Thus it is the product of nnn copies of the matrix [ψ(a,x)ψ(a′,x′)‾][\psi(a,x)\overline{\psi(a',x')}][ψ(a,x)ψ(a′,x′)​] and the shared resource, with the displayed register ordering. When n=0n=0n=0, the product is one and the entries reduce to those of τ\tauτ.

QAlgorithms.Papers.Independent.Intakethird100.P27.LossyCode

For finite types A,BA,BA,B and any natural number nnn, a lossy code consists of positive natural numbers q,c,hq,c,hq,c,h, representing its quantum, classical, and shared dimensions; a positive semidefinite matrix τ\tauτ of trace one on {0,…,h−1}2\{0,\ldots,h-1\}^2{0,…,h−1}2; an encoder from matrices on An×{0,…,h−1}A^n\times\{0,\ldots,h-1\}An×{0,…,h−1} to matrices on {0,…,c−1}×{0,…,q−1}\{0,\ldots,c-1\}\times\{0,\ldots,q-1\}{0,…,c−1}×{0,…,q−1}; and a decoder from matrices on ({0,…,c−1}×{0,…,q−1})×{0,…,h−1}(\{0,\ldots,c-1\}\times\{0,\ldots,q-1\})\times\{0,\ldots,h-1\}({0,…,c−1}×{0,…,q−1})×{0,…,h−1} to matrices on BnB^nBn. Each map must admit a finite Kraus representation: it sends every matrix XXX to ∑jKjXKj†\sum_j K_jXK_j^\dagger∑j​Kj​XKj†​, with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I on its input space. Moreover, for every encoder input matrix XXX, every distinct classical indices m,m′m,m'm,m′, and all quantum indices w,w′w,w'w,w′, the encoder output satisfies Enc⁡(X)(m,w),(m′,w′)=0\operatorname{Enc}(X)_{(m,w),(m',w')}=0Enc(X)(m,w),(m′,w′)​=0. Here AnA^nAn and BnB^nBn denote length-nnn tuples; n=0n=0n=0 is allowed, with each tuple register then a singleton.

QAlgorithms.Papers.Independent.Intakethird100.P27.LossyCode.output

For finite types A,BA,BA,B, any natural number nnn, a lossy code PPP with positive quantum, classical, and shared dimensions q,c,hq,c,hq,c,h, shared density matrix τ\tauτ, finite Kraus trace-preserving encoder Enc⁡:An×[h]→[c]×[q]\operatorname{Enc}:A^n\times[h]\to[c]\times[q]Enc:An×[h]→[c]×[q] and decoder Dec⁡:([c]×[q])×[h]→Bn\operatorname{Dec}:([c]\times[q])\times[h]\to B^nDec:([c]×[q])×[h]→Bn, where these arrows denote maps between matrix spaces and the encoder is block diagonal in its classical output index, and a source (ψ,Δ)(\psi,\Delta)(ψ,Δ) with ∥ψ∥=1\|\psi\|=1∥ψ∥=1 and Δ≥0\Delta\ge0Δ≥0, the output is a matrix on An×BnA^n\times B^nAn×Bn. To define its entries, for each r,r′∈Anr,r'\in A^nr,r′∈An form the matrix Yr,r′Y_{r,r'}Yr,r′​ on An×[h]×[h]A^n\times[h]\times[h]An×[h]×[h] with entries

Yr,r′((x,u),v),((x′,u′),v′)=∏i=0n−1ψ(ri,xi)ψ(ri′,xi′)‾ τ(u,v),(u′,v′).Y_{r,r'}{}_{((x,u),v),((x',u'),v')} = \prod_{i=0}^{n-1}\psi(r_i,x_i)\overline{\psi(r'_i,x'_i)} \,\tau_{(u,v),(u',v')}.Yr,r′​((x,u),v),((x′,u′),v′)​=i=0∏n−1​ψ(ri​,xi​)ψ(ri′​,xi′​)​τ(u,v),(u′,v′)​.

Apply the encoder to its An×[h]A^n\times[h]An×[h] matrix blocks while retaining the second shared register, producing Zr,r′Z_{r,r'}Zr,r′​ with entries Zr,r′((m,w),v),((m′,w′),v′)=Enc⁡([Yr,r′((x,u),v),((x′,u′),v′)](x,u),(x′,u′))(m,w),(m′,w′)Z_{r,r'}{}_{((m,w),v),((m',w'),v')}=\operatorname{Enc}([Y_{r,r'}{}_{((x,u),v),((x',u'),v')}]_{(x,u),(x',u')})_{(m,w),(m',w')}Zr,r′​((m,w),v),((m′,w′),v′)​=Enc([Yr,r′​((x,u),v),((x′,u′),v′)​](x,u),(x′,u′)​)(m,w),(m′,w′)​; the output entry at (r,b),(r′,b′)(r,b),(r',b')(r,b),(r′,b′) is then Dec⁡(Zr,r′)b,b′\operatorname{Dec}(Z_{r,r'})_{b,b'}Dec(Zr,r′​)b,b′​. Here [d]={0,…,d−1}[d]=\{0,\ldots,d-1\}[d]={0,…,d−1}; n=0n=0n=0 is allowed and all products over its positions equal one.

QAlgorithms.Papers.Independent.Intakethird100.P27.averageDistortion

For finite types A,BA,BA,B, a source consisting of a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B, any natural number nnn, and any complex matrix XXX on An×BnA^n\times B^nAn×Bn, the average distortion is

1n∑i=0n−1Re⁡tr⁡(ΔXi),(Xi)(a,b),(a′,b′)=∑ui=a, vi=b, ui′=a′, vi′=b′∀j≠i, (uj,vj)=(uj′,vj′)X(u,v),(u′,v′),\frac1n\sum_{i=0}^{n-1}\operatorname{Re}\operatorname{tr}(\Delta X_i), \qquad (X_i)_{(a,b),(a',b')}= \sum_{\substack{u_i=a,\ v_i=b,\ u'_i=a',\ v'_i=b'\\ \forall j\ne i,\ (u_j,v_j)=(u'_j,v'_j)}} X_{(u,v),(u',v')},n1​i=0∑n−1​Retr(ΔXi​),(Xi​)(a,b),(a′,b′)​=ui​=a, vi​=b, ui′​=a′, vi′​=b′∀j=i, (uj​,vj​)=(uj′​,vj′​)​∑​X(u,v),(u′,v′)​,

where the inner sum ranges over u,u′∈Anu,u'\in A^nu,u′∈An and v,v′∈Bnv,v'\in B^nv,v′∈Bn. No density-matrix assumption is imposed on XXX. For n=0n=0n=0, the numerator is an empty sum and real division by zero is defined to be zero, so the average distortion is zero.

QAlgorithms.Papers.Independent.Intakethird100.P27.conditionalOutputInformation

For finite types R,A,B,CR,A,B,CR,A,B,C, any function NNN from complex matrices on AAA to complex matrices on CCC, and any vector ϕ∈CR×(A×B)\phi\in\mathbb C^{R\times(A\times B)}ϕ∈CR×(A×B), define the matrix XXX on R×(C×B)R\times(C\times B)R×(C×B) by

X(r,(c,b)),(r′,(c′,b′))=N ⁣([ϕ(r,(a,b)) ϕ(r′,(a′,b′))‾]a,a′∈A)c,c′.X_{(r,(c,b)),(r',(c',b'))} = N\!\left( \left[\phi_{(r,(a,b))}\, \overline{\phi_{(r',(a',b'))}}\right]_{a,a'\in A} \right)_{c,c'}.X(r,(c,b)),(r′,(c′,b′))​=N([ϕ(r,(a,b))​ϕ(r′,(a′,b′))​​]a,a′∈A​)c,c′​.

The declaration defines the real number

H(XRB)+H(XCB)−H(XB)−H(X),H(X_{RB})+H(X_{CB})-H(X_B)-H(X),H(XRB​)+H(XCB​)−H(XB​)−H(X),

where

(XRB)(r,b),(r′,b′)=∑cX(r,(c,b)),(r′,(c,b′)),(XCB)(c,b),(c′,b′)=∑rX(r,(c,b)),(r,(c′,b′)),(XB)b,b′=∑c∑rX(r,(c,b)),(r,(c,b′)),(X_{RB})_{(r,b),(r',b')}=\sum_cX_{(r,(c,b)),(r',(c,b'))}, \quad (X_{CB})_{(c,b),(c',b')}=\sum_rX_{(r,(c,b)),(r,(c',b'))}, \quad (X_B)_{b,b'}=\sum_c\sum_rX_{(r,(c,b)),(r,(c,b'))},(XRB​)(r,b),(r′,b′)​=c∑​X(r,(c,b)),(r′,(c,b′))​,(XCB​)(c,b),(c′,b′)​=r∑​X(r,(c,b)),(r,(c′,b′))​,(XB​)b,b′​=c∑​r∑​X(r,(c,b)),(r,(c,b′))​,

and H(Y)=−Re⁡tr⁡(f(Y))H(Y)=-\operatorname{Re}\operatorname{tr}(f(Y))H(Y)=−Retr(f(Y)), with f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x applied by continuous functional calculus. No normalization of ϕ\phiϕ, channel property of NNN, or nonemptiness of the finite types is assumed; the formula also covers empty index types.

QAlgorithms.Papers.Independent.Intakethird100.P27.formationBlockValue

For finite types A,BA,BA,B, a source SSS with unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and positive semidefinite distortion matrix Δ\DeltaΔ on A×BA\times BA×B, any real DDD, and any k∈Nk\in\mathbb Nk∈N, define

max⁡(k−1,0) inf⁡N a channeldˉS(XN)≤Dmax⁡(EF(XN),0)in [0,∞].\max(k^{-1},0)\, \inf_{\substack{N\text{ a channel}\\ \bar d_S(X_N)\le D}} \max(E_F(X_N),0) \quad\text{in }[0,\infty].max(k−1,0)N a channeldˉS​(XN​)≤D​inf​max(EF​(XN​),0)in [0,∞].

Here a channel from matrices on AkA^kAk to matrices on BkB^kBk means a map N(Z)=∑jKjZKj†N(Z)=\sum_jK_jZK_j^\daggerN(Z)=∑j​Kj​ZKj†​ for finitely many matrices KjK_jKj​, with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I; its output XNX_NXN​ is specified by

XN((a,b),(a′,b′))=N ⁣((∏i<kψ(ai,xi)ψ(ai′,yi)‾)x,y)(b,b′).X_N((a,b),(a',b')) = N\!\left( \left(\prod_{i<k}\psi(a_i,x_i)\overline{\psi(a'_i,y_i)}\right)_{x,y} \right)(b,b').XN​((a,b),(a′,b′))=N​(i<k∏​ψ(ai​,xi​)ψ(ai′​,yi​)​)x,y​​(b,b′).

For a matrix XXX on Ak×BkA^k\times B^kAk×Bk, dˉS(X)=k−1∑i<kRe⁡Tr⁡(ΔXi)\bar d_S(X)=k^{-1}\sum_{i<k}\operatorname{Re}\operatorname{Tr}(\Delta X_i)dˉS​(X)=k−1∑i<k​ReTr(ΔXi​), where XiX_iXi​ sums over equal row and column indices outside position iii. The real quantity EF(X)E_F(X)EF​(X) is the infimum over all finite lists of nonnegative weights pjp_jpj​ summing to 111 and unit vectors uj∈CAk×Bku_j\in\mathbb C^{A^k\times B^k}uj​∈CAk×Bk satisfying X=∑jpjujuj†X=\sum_jp_ju_ju_j^\daggerX=∑j​pj​uj​uj†​, of

∑jpjS2 ⁣((∑b∈Bkuj(a,b)uj(a′,b)‾)a,a′),S2(H)=−Re⁡Tr⁡((x↦xlog⁡2x)(H)).\sum_jp_jS_2\!\left( \left(\sum_{b\in B^k}u_j(a,b)\overline{u_j(a',b)}\right)_{a,a'} \right), \qquad S_2(H)=-\operatorname{Re}\operatorname{Tr}\bigl((x\mapsto x\log_2x)(H)\bigr).j∑​pj​S2​​(b∈Bk∑​uj​(a,b)uj​(a′,b)​)a,a′​​,S2​(H)=−ReTr((x↦xlog2​x)(H)).

The function of a matrix here is given by continuous functional calculus, with 0log⁡20=00\log_2 0=00log2​0=0; the real infimum defining EFE_FEF​ is 000 if its defining set is empty. The outer constrained infimum is instead taken in [0,∞][0,\infty][0,∞] and is +∞+\infty+∞ if there is no feasible channel. Inverses at zero are zero, and multiplication in [0,∞][0,\infty][0,∞] has 0⋅∞=00\cdot\infty=00⋅∞=0; consequently the defined block value is 000 when k=0k=0k=0, including when the constraint is infeasible.

QAlgorithms.Papers.Independent.Intakethird100.P27.iidChannelOutput

For finite types A,BA,BA,B, let SSS contain a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B. For any natural number kkk, including zero, and any matrix function NNN from matrices on AkA^kAk to matrices on BkB^kBk, this defines the matrix XXX on Ak×BkA^k\times B^kAk×Bk with entries

X((a,b),(a′,b′))=N ⁣((∏i=0k−1ψ(ai,xi)ψ(ai′,yi)‾)x,y∈Ak)(b,b′).X((a,b),(a',b')) = N\!\left( \left( \prod_{i=0}^{k-1} \psi(a_i,x_i)\overline{\psi(a'_i,y_i)} \right)_{x,y\in A^k} \right)(b,b').X((a,b),(a′,b′))=N​(i=0∏k−1​ψ(ai​,xi​)ψ(ai′​,yi​)​)x,y∈Ak​​(b,b′).

Here AkA^kAk and BkB^kBk mean functions from Fin⁡(k)\operatorname{Fin}(k)Fin(k) to the respective types. There is no channel assumption on NNN, and Δ\DeltaΔ does not enter this definition. For k=0k=0k=0, both function spaces have one element and the empty product is 111, so the output is the one-by-one matrix obtained by applying NNN to the one-by-one matrix [1][1][1].

QAlgorithms.Papers.Independent.Intakethird100.P27.isotropicConvexUpper

For every real DDD, this defines

U(D)=inf⁡m∈N, m>0, pi≥0, ∑i=1mpi=10≤di≤1, ∑i=1mpidi≤D ∑i=1mpiC(di)in [0,∞],U(D)= \inf_{\substack{m\in\mathbb N,\ m>0,\ p_i\ge0,\ \sum_{i=1}^m p_i=1\\ 0\le d_i\le1,\ \sum_{i=1}^m p_id_i\le D}} \ \sum_{i=1}^m p_i C(d_i) \quad\text{in }[0,\infty],U(D)=m∈N, m>0, pi​≥0, ∑i=1m​pi​=10≤di​≤1, ∑i=1m​pi​di​≤D​inf​ i=1∑m​pi​C(di​)in [0,∞],

where C(d)=0C(d)=0C(d)=0 for d>3/4d>3/4d>3/4, and C(d)=Ep((1−d)P+d3(I−P))C(d)=E_{\mathrm p}((1-d)P+\frac d3(I-P))C(d)=Ep​((1−d)P+3d​(I−P)) for d≤3/4d\le3/4d≤3/4. The matrix PPP on two two-dimensional registers has entries P(i,j),(k,l)=1/2P_{(i,j),(k,l)}=1/2P(i,j),(k,l)​=1/2 when i=ji=ji=j and k=lk=lk=l, and zero otherwise. For a matrix ρ\rhoρ on these two registers, Ep(ρ)E_{\mathrm p}(\rho)Ep​(ρ) is the infimum over positive finite ancillary dimensions u,vu,vu,v and unit vectors ξ\xiξ on ({0,1}×{1,…,u})×({0,1}×{1,…,v})(\{0,1\}\times\{1,\ldots,u\})\times(\{0,1\}\times\{1,\ldots,v\})({0,1}×{1,…,u})×({0,1}×{1,…,v}) whose marginal on the original registers is ρ\rhoρ, of the nonnegative part of H(Tr⁡{0,1}×{1,…,v}∣ξ⟩⟨ξ∣)H(\operatorname{Tr}_{\{0,1\}\times\{1,\ldots,v\}}|\xi\rangle\langle\xi|)H(Tr{0,1}×{1,…,v}​∣ξ⟩⟨ξ∣). Here (∣ξ⟩⟨ξ∣)xy=ξxξy‾(|\xi\rangle\langle\xi|)_{xy}=\xi_x\overline{\xi_y}(∣ξ⟩⟨ξ∣)xy​=ξx​ξy​​, and H(X)=−Re⁡tr⁡(f(X))H(X)=-\operatorname{Re}\operatorname{tr}(f(X))H(X)=−Retr(f(X)), with f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x applied by continuous functional calculus. Zero weights and repeated did_idi​ are allowed, and the weighted mean distortion need only be at most DDD. Products and infima use extended nonnegative arithmetic, including 0⋅(+∞)=00\cdot(+\infty)=00⋅(+∞)=0 and inf⁡∅=+∞\inf\varnothing=+\inftyinf∅=+∞; in particular, when D<0D<0D<0, the defining feasible set is empty.

QAlgorithms.Papers.Independent.Intakethird100.P27.isotropicEntropy

For every real DDD, with no restriction to a probability range, define p0=1−Dp_0=1-Dp0​=1−D and p1=p2=p3=D/3p_1=p_2=p_3=D/3p1​=p2​=p3​=D/3. The defined real number is

−∑i=03pilog⁡2pi=−(1−D)log⁡2(1−D)−3D3log⁡2D3.-\sum_{i=0}^{3}p_i\log_2p_i = -(1-D)\log_2(1-D)-3\frac D3\log_2\frac D3.−i=0∑3​pi​log2​pi​=−(1−D)log2​(1−D)−33D​log2​3D​.

The logarithm is the total real function log⁡2x=log⁡x/log⁡2\log_2x=\log x/\log2log2​x=logx/log2, including its value 000 at x=0x=0x=0; thus a zero coefficient contributes zero. For D∉[0,1]D\notin[0,1]D∈/[0,1], the same formula is used without requiring the four coefficients to be nonnegative.

QAlgorithms.Papers.Independent.Intakethird100.P27.isotropicStrategyCost

For every real DDD, this defines an extended nonnegative real number

C(D)={Ep ⁣((1−D)P+D3(I−P)),D≤34,0,D>34,C(D)= \begin{cases} E_{\mathrm p}\!\left((1-D)P+\dfrac D3(I-P)\right),&D\le \dfrac34,\\[2mm] 0,&D>\dfrac34, \end{cases}C(D)=⎩⎨⎧​Ep​((1−D)P+3D​(I−P)),0,​D≤43​,D>43​,​

where PPP is the matrix on {0,1}×{0,1}\{0,1\}\times\{0,1\}{0,1}×{0,1} whose entry P(i,j),(k,l)P_{(i,j),(k,l)}P(i,j),(k,l)​ is 1/21/21/2 when i=ji=ji=j and k=lk=lk=l, and is zero otherwise. Here Ep(ρ)E_{\mathrm p}(\rho)Ep​(ρ) is the infimum in [0,∞][0,\infty][0,∞] over positive natural dimensions u,vu,vu,v and unit vectors ξ\xiξ on ({0,1}×{1,…,u})×({0,1}×{1,…,v})(\{0,1\}\times\{1,\ldots,u\})\times(\{0,1\}\times\{1,\ldots,v\})({0,1}×{1,…,u})×({0,1}×{1,…,v}) whose marginal on the two {0,1}\{0,1\}{0,1} registers equals ρ\rhoρ, of the nonnegative part of H(Tr⁡{0,1}×{1,…,v}∣ξ⟩⟨ξ∣)H(\operatorname{Tr}_{\{0,1\}\times\{1,\ldots,v\}}|\xi\rangle\langle\xi|)H(Tr{0,1}×{1,…,v}​∣ξ⟩⟨ξ∣). The matrix ∣ξ⟩⟨ξ∣|\xi\rangle\langle\xi|∣ξ⟩⟨ξ∣ has entries ξxξy‾\xi_x\overline{\xi_y}ξx​ξy​​, and H(X)=−Re⁡tr⁡(f(X))H(X)=-\operatorname{Re}\operatorname{tr}(f(X))H(X)=−Retr(f(X)), with f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x applied by continuous functional calculus. An empty infimum is +∞+\infty+∞. No restriction on DDD is imposed: the first branch includes D=3/4D=3/4D=3/4 and every negative DDD, and the second includes every D>3/4D>3/4D>3/4.

QAlgorithms.Papers.Independent.Intakethird100.P27.letterMarginal

For finite types A,BA,BA,B, any natural number nnn, any complex matrix XXX indexed by An×BnA^n\times B^nAn×Bn, and any i∈{0,…,n−1}i\in\{0,\ldots,n-1\}i∈{0,…,n−1}, the letter marginal at iii is the matrix on A×BA\times BA×B given by

(Xi)(a,b),(a′,b′)=∑u,u′∈An, v,v′∈Bnui=a, vi=b, ui′=a′, vi′=b′∀j≠i, (uj,vj)=(uj′,vj′)X(u,v),(u′,v′).(X_i)_{(a,b),(a',b')}= \sum_{\substack{u,u'\in A^n,\ v,v'\in B^n\\ u_i=a,\ v_i=b,\ u'_i=a',\ v'_i=b'\\ \forall j\ne i,\ (u_j,v_j)=(u'_j,v'_j)}} X_{(u,v),(u',v')}.(Xi​)(a,b),(a′,b′)​=u,u′∈An, v,v′∈Bnui​=a, vi​=b, ui′​=a′, vi′​=b′∀j=i, (uj​,vj​)=(uj′​,vj′​)​∑​X(u,v),(u′,v′)​.

Here AnA^nAn and BnB^nBn are the sets of length-nnn tuples. No positivity or trace assumption is imposed on XXX. When n=0n=0n=0, there is no index iii, so there is no instance of this definition with such an argument.

QAlgorithms.Papers.Independent.Intakethird100.P27.observableDistortion

For finite types A,BA,BA,B, a source consisting of a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B, and any function NNN from complex matrices on AAA to complex matrices on BBB, the observable distortion is Re⁡tr⁡(ΔT)\operatorname{Re}\operatorname{tr}(\Delta T)Retr(ΔT), where

T(a,b),(a′,b′)=N ⁣([ψ(a,x)ψ(a′,x′)‾]x,x′∈A)b,b′.T_{(a,b),(a',b')}= N\!\left(\left[\psi(a,x)\overline{\psi(a',x')}\right]_{x,x'\in A}\right)_{b,b'}.T(a,b),(a′,b′)​=N([ψ(a,x)ψ(a′,x′)​]x,x′∈A​)b,b′​.

The definition imposes no channel hypothesis on NNN, and Δ\DeltaΔ need not be bounded by the identity.

QAlgorithms.Papers.Independent.Intakethird100.P27.outputWithSide

For arbitrary types R,A,B,CR,A,B,CR,A,B,C, any function NNN from complex matrices on AAA to complex matrices on CCC, and any vector ϕ∈CR×(A×B)\phi\in\mathbb C^{R\times(A\times B)}ϕ∈CR×(A×B) equipped with the stated Euclidean-space structure, this defines a matrix XXX on R×(C×B)R\times(C\times B)R×(C×B) by

X(r,(c,b)),(r′,(c′,b′))=N ⁣([ϕ(r,(a,b)) ϕ(r′,(a′,b′))‾]a,a′∈A)c,c′.X_{(r,(c,b)),(r',(c',b'))} = N\!\left( \left[\phi_{(r,(a,b))}\, \overline{\phi_{(r',(a',b'))}}\right]_{a,a'\in A} \right)_{c,c'}.X(r,(c,b)),(r′,(c′,b′))​=N([ϕ(r,(a,b))​ϕ(r′,(a′,b′))​​]a,a′∈A​)c,c′​.

Thus NNN is applied separately to each matrix block obtained by fixing the two reference indices and the two BBB-indices. There are no finiteness, nonemptiness, normalization, linearity, positivity, or channel hypotheses.

QAlgorithms.Papers.Independent.Intakethird100.P27.purificationBlockValue

For finite types A,BA,BA,B, a source SSS containing a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite distortion matrix Δ\DeltaΔ on A×BA\times BA×B, any real DDD, and any k∈Nk\in\mathbb Nk∈N, define

max⁡(k−1,0) inf⁡N a channeldˉS(XN)≤DP(XN)in [0,∞].\max(k^{-1},0)\, \inf_{\substack{N\text{ a channel}\\\bar d_S(X_N)\le D}}\mathcal P(X_N) \quad\text{in }[0,\infty].max(k−1,0)N a channeldˉS​(XN​)≤D​inf​P(XN​)in [0,∞].

Here NNN ranges over maps from matrices on AkA^kAk to matrices on BkB^kBk admitting finite Kraus representations N(Z)=∑jKjZKj†N(Z)=\sum_jK_jZK_j^\daggerN(Z)=∑j​Kj​ZKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I, and

XN((a,b),(a′,b′))=N ⁣((∏i<kψ(ai,xi)ψ(ai′,yi)‾)x,y)(b,b′).X_N((a,b),(a',b')) = N\!\left( \left(\prod_{i<k}\psi(a_i,x_i)\overline{\psi(a'_i,y_i)}\right)_{x,y} \right)(b,b').XN​((a,b),(a′,b′))=N​(i<k∏​ψ(ai​,xi​)ψ(ai′​,yi​)​)x,y​​(b,b′).

The distortion is dˉS(X)=k−1∑i<kRe⁡Tr⁡(ΔXi)\bar d_S(X)=k^{-1}\sum_{i<k}\operatorname{Re}\operatorname{Tr}(\Delta X_i)dˉS​(X)=k−1∑i<k​ReTr(ΔXi​), where XiX_iXi​ sums over equal row and column indices outside position iii. For a matrix XXX on Ak×BkA^k\times B^kAk×Bk, P(X)\mathcal P(X)P(X) is the infimum in [0,∞][0,\infty][0,∞] of max⁡(S2(Hϕ),0)\max(S_2(H_\phi),0)max(S2​(Hϕ​),0) over positive integers u,vu,vu,v and unit vectors ϕ\phiϕ on (Ak×Fin⁡(u))×(Bk×Fin⁡(v))(A^k\times\operatorname{Fin}(u))\times(B^k\times\operatorname{Fin}(v))(Ak×Fin(u))×(Bk×Fin(v)) whose matrix ϕϕ†\phi\phi^\daggerϕϕ†, after summing over equal auxiliary indices, equals XXX. The matrix HϕH_\phiHϕ​ is obtained by summing ϕϕ†\phi\phi^\daggerϕϕ† over equal Bk×Fin⁡(v)B^k\times\operatorname{Fin}(v)Bk×Fin(v) indices, and S2(H)=−Re⁡Tr⁡((x↦xlog⁡2x)(H))S_2(H)=-\operatorname{Re}\operatorname{Tr}((x\mapsto x\log_2x)(H))S2​(H)=−ReTr((x↦xlog2​x)(H)), with 0log⁡20=00\log_2 0=00log2​0=0. Either infimum is +∞+\infty+∞ if its indexing set is empty. Zero inverses and 0⋅∞0\cdot\infty0⋅∞ are zero, so the defined block value is 000 for k=0k=0k=0, irrespective of feasibility.

QAlgorithms.Papers.Independent.Intakethird100.P27.purificationCost

For finite types A,BA,BA,B and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the infimum in [0,∞][0,\infty][0,∞] of max⁡(S2(Hϕ),0)\max(S_2(H_\phi),0)max(S2​(Hϕ​),0) over all positive integers u,vu,vu,v and unit vectors ϕ∈C(A×Fin⁡(u))×(B×Fin⁡(v))\phi\in\mathbb C^{(A\times\operatorname{Fin}(u))\times(B\times\operatorname{Fin}(v))}ϕ∈C(A×Fin(u))×(B×Fin(v)) such that

ρ((a,b),(a′,b′))=∑α<u, β<vϕ((a,α),(b,β))ϕ((a′,α),(b′,β))‾.\rho((a,b),(a',b')) = \sum_{\alpha<u,\ \beta<v} \phi((a,\alpha),(b,\beta)) \overline{\phi((a',\alpha),(b',\beta))}.ρ((a,b),(a′,b′))=α<u, β<v∑​ϕ((a,α),(b,β))ϕ((a′,α),(b′,β))​.

Here HϕH_\phiHϕ​ is the matrix on A×Fin⁡(u)A\times\operatorname{Fin}(u)A×Fin(u) obtained by summing ϕϕ†\phi\phi^\daggerϕϕ† over equal B×Fin⁡(v)B\times\operatorname{Fin}(v)B×Fin(v) indices, and

S2(Hϕ)=−Re⁡Tr⁡((x↦xlog⁡2x)(Hϕ)),S_2(H_\phi) =-\operatorname{Re}\operatorname{Tr}\bigl((x\mapsto x\log_2x)(H_\phi)\bigr),S2​(Hϕ​)=−ReTr((x↦xlog2​x)(Hϕ​)),

with continuous functional calculus and 0log⁡20=00\log_2 0=00log2​0=0. The auxiliary dimensions have no upper bound. The definition does not assume that ρ\rhoρ is a density matrix and does not assert attainment. If no such purification exists, including when either underlying type is empty, the value is +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P27.purificationValues

For finite types A,BA,BA,B and an arbitrary complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the set of all r∈[0,∞]r\in[0,\infty]r∈[0,∞] for which there exist positive integers u,vu,vu,v and a unit vector ϕ∈C(A×Fin⁡(u))×(B×Fin⁡(v))\phi\in\mathbb C^{(A\times\operatorname{Fin}(u))\times(B\times\operatorname{Fin}(v))}ϕ∈C(A×Fin(u))×(B×Fin(v)) satisfying

ρ((a,b),(a′,b′))=∑α<u∑β<vϕ((a,α),(b,β))ϕ((a′,α),(b′,β))‾,r=max⁡(S2(Hϕ),0),\rho((a,b),(a',b')) = \sum_{\alpha<u}\sum_{\beta<v} \phi((a,\alpha),(b,\beta)) \overline{\phi((a',\alpha),(b',\beta))}, \qquad r=\max(S_2(H_\phi),0),ρ((a,b),(a′,b′))=α<u∑​β<v∑​ϕ((a,α),(b,β))ϕ((a′,α),(b′,β))​,r=max(S2​(Hϕ​),0),

where

Hϕ((a,α),(a′,α′))=∑b∈B∑β<vϕ((a,α),(b,β))ϕ((a′,α′),(b,β))‾,H_\phi((a,\alpha),(a',\alpha')) = \sum_{b\in B}\sum_{\beta<v} \phi((a,\alpha),(b,\beta)) \overline{\phi((a',\alpha'),(b,\beta))},Hϕ​((a,α),(a′,α′))=b∈B∑​β<v∑​ϕ((a,α),(b,β))ϕ((a′,α′),(b,β))​,

and S2(H)=−Re⁡Tr⁡((x↦xlog⁡2x)(H))S_2(H)=-\operatorname{Re}\operatorname{Tr}((x\mapsto x\log_2x)(H))S2​(H)=−ReTr((x↦xlog2​x)(H)), using continuous functional calculus and 0log⁡20=00\log_2 0=00log2​0=0. No density assumption on ρ\rhoρ is a binder; membership requires the displayed normalized purification. There is no upper bound on u,vu,vu,v. If either AAA or BBB is empty, no such unit vector exists and the set is empty.

QAlgorithms.Papers.Independent.Intakethird100.P27.rateDistortion

For finite types A,BA,BA,B, a source SSS consists of a vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A with ∥ψ∥=1\|\psi\|=1∥ψ∥=1 and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B. For an assistance choice a∈{unassisted,entanglement,classical}a\in\{\text{unassisted},\text{entanglement},\text{classical}\}a∈{unassisted,entanglement,classical} and any real DDD, this defines a value in [0,∞][0,\infty][0,∞] by

inf⁡r∈Rr satisfies the conditions belowmax⁡(r,0).\inf_{\substack{r\in\mathbb R\\ r\text{ satisfies the conditions below}}}\max(r,0).r∈Rr satisfies the conditions below​inf​max(r,0).

The conditions on rrr are r≥0r\ge0r≥0, D≥0D\ge0D≥0, and the following: for every ε>0\varepsilon>0ε>0 and s>0s>0s>0, there is n0∈Nn_0\in\mathbb Nn0​∈N such that for every n≥n0n\ge n_0n≥n0​ with n≥1n\ge1n≥1, there is a code with positive integer quantum, classical, and shared dimensions q,c,hq,c,hq,c,h; a positive semidefinite trace-one shared matrix τ\tauτ on Fin⁡(h)×Fin⁡(h)\operatorname{Fin}(h)\times\operatorname{Fin}(h)Fin(h)×Fin(h); an encoder from An×Fin⁡(h)A^n\times\operatorname{Fin}(h)An×Fin(h) to Fin⁡(c)×Fin⁡(q)\operatorname{Fin}(c)\times\operatorname{Fin}(q)Fin(c)×Fin(q); and a decoder from (Fin⁡(c)×Fin⁡(q))×Fin⁡(h)(\operatorname{Fin}(c)\times\operatorname{Fin}(q))\times\operatorname{Fin}(h)(Fin(c)×Fin(q))×Fin(h) to BnB^nBn. Both maps must admit finite Kraus representations Φ(X)=∑jKjXKj†\Phi(X)=\sum_j K_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. Every encoder output must have zero entries between distinct classical labels. The assistance restrictions are c=h=1c=h=1c=h=1 for unassisted codes, c=1c=1c=1 for entanglement assistance, and h=1h=1h=1 for classical assistance. Apply the encoder and decoder to the second register of nnn copies of ψ\psiψ, together with τ\tauτ, leaving the first AnA^nAn register untouched, and call the resulting matrix XXX on An×BnA^n\times B^nAn×Bn. The required bounds are

log⁡2q≤n(r+s),1n∑i=0n−1Re⁡Tr⁡(ΔXi)≤D+ε,\log_2q\le n(r+s), \qquad \frac1n\sum_{i=0}^{n-1}\operatorname{Re}\operatorname{Tr}(\Delta X_i)\le D+\varepsilon,log2​q≤n(r+s),n1​i=0∑n−1​ReTr(ΔXi​)≤D+ε,

where XiX_iXi​ is obtained by summing over equal row and column indices at all positions other than iii. The infimum of an empty family is +∞+\infty+∞; in particular the value is +∞+\infty+∞ when D<0D<0D<0. No positive lower bound on the block length beyond n≥1n\ge1n≥1, and no bound on the classical or shared dimensions beyond the stated assistance restrictions, is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P27.regularizedPurification

For finite types A,BA,BA,B, a source SSS containing a unit vector ψ∈CA×A\psi\in\mathbb C^{A\times A}ψ∈CA×A and a positive semidefinite matrix Δ\DeltaΔ on A×BA\times BA×B, and any real DDD, this defines the extended nonnegative number

lim inf⁡k→∞bk+1=sup⁡K∈Ninf⁡k≥Kbk+1.\liminf_{k\to\infty}b_{k+1} = \sup_{K\in\mathbb N}\inf_{k\ge K}b_{k+1}.k→∞liminf​bk+1​=K∈Nsup​k≥Kinf​bk+1​.

For each positive integer jjj,

bj=1jinf⁡N a channelj−1∑i<jRe⁡Tr⁡(Δ(XN)i)≤DP(XN),XN((a,b),(a′,b′))=N ⁣((∏i<jψ(ai,xi)ψ(ai′,yi)‾)x,y)(b,b′).b_j=\frac1j \inf_{\substack{N\text{ a channel}\\ j^{-1}\sum_{i<j}\operatorname{Re}\operatorname{Tr}(\Delta(X_N)_i)\le D}} \mathcal P(X_N), \qquad X_N((a,b),(a',b')) = N\!\left( \left(\prod_{i<j}\psi(a_i,x_i)\overline{\psi(a'_i,y_i)}\right)_{x,y} \right)(b,b').bj​=j1​N a channelj−1∑i<j​ReTr(Δ(XN​)i​)≤D​inf​P(XN​),XN​((a,b),(a′,b′))=N​(i<j∏​ψ(ai​,xi​)ψ(ai′​,yi​)​)x,y​​(b,b′).

The channel NNN maps matrices on AjA^jAj to matrices on BjB^jBj and must admit a finite representation N(Z)=∑ℓKℓZKℓ†N(Z)=\sum_\ell K_\ell ZK_\ell^\daggerN(Z)=∑ℓ​Kℓ​ZKℓ†​ with ∑ℓKℓ†Kℓ=I\sum_\ell K_\ell^\dagger K_\ell=I∑ℓ​Kℓ†​Kℓ​=I. Each (XN)i(X_N)_i(XN​)i​ sums over equal row and column indices at every position except iii. For any matrix XXX on Aj×BjA^j\times B^jAj×Bj, P(X)\mathcal P(X)P(X) is the infimum in [0,∞][0,\infty][0,∞] of max⁡(S2(Hϕ),0)\max(S_2(H_\phi),0)max(S2​(Hϕ​),0) over positive integers u,vu,vu,v and unit vectors ϕ\phiϕ on (Aj×Fin⁡(u))×(Bj×Fin⁡(v))(A^j\times\operatorname{Fin}(u))\times(B^j\times\operatorname{Fin}(v))(Aj×Fin(u))×(Bj×Fin(v)) whose matrix ϕϕ†\phi\phi^\daggerϕϕ†, after summing over equal auxiliary indices, equals XXX; HϕH_\phiHϕ​ sums that matrix over equal Bj×Fin⁡(v)B^j\times\operatorname{Fin}(v)Bj×Fin(v) indices, and S2(H)=−Re⁡Tr⁡((x↦xlog⁡2x)(H))S_2(H)=-\operatorname{Re}\operatorname{Tr}((x\mapsto x\log_2x)(H))S2​(H)=−ReTr((x↦xlog2​x)(H)), with 0log⁡20=00\log_2 0=00log2​0=0. Empty infima are +∞+\infty+∞. Only positive block lengths enter this limit inferior, and the definition does not assert that the sequence converges or that any infimum is attained.

QAlgorithms.Papers.Independent.Intakethird100.P27.sourceFidelity

For arbitrary finite types R,AR,AR,A, any vector v∈CR×Av\in\mathbb C^{R\times A}v∈CR×A, and any function NNN from complex matrices on AAA to complex matrices on AAA, define Qr,sQ_{r,s}Qr,s​ by

Qr,s(x,y)=v(r,x)v(s,y)‾.Q_{r,s}(x,y)=v(r,x)\overline{v(s,y)}.Qr,s​(x,y)=v(r,x)v(s,y)​.

The defined real number is

Re⁡∑r,s∈R∑a,b∈Av(r,a)‾ N(Qr,s)(a,b) v(s,b).\operatorname{Re}\sum_{r,s\in R}\sum_{a,b\in A} \overline{v(r,a)}\,N(Q_{r,s})(a,b)\,v(s,b).Rer,s∈R∑​a,b∈A∑​v(r,a)​N(Qr,s​)(a,b)v(s,b).

Thus NNN is applied separately to each AAA-matrix block of the matrix vv†vv^\daggervv†, and the resulting matrix is evaluated against vvv. Neither normalization of vvv nor linearity, positivity, or trace preservation of NNN is assumed. Empty finite types are included, in which case the displayed sum is zero.

QAlgorithms.Papers.Independent.Intakethird100.P27.testOutput

For finite types A,BA,BA,B, a source S=(ψ,Δ)S=(\psi,\Delta)S=(ψ,Δ) with ∥ψ∥=1\|\psi\|=1∥ψ∥=1 and Δ≥0\Delta\ge0Δ≥0, and any function NNN from complex A×AA\times AA×A matrices to complex B×BB\times BB×B matrices, the test output is the matrix TTT indexed by A×BA\times BA×B whose entries are

T(a,b),(a′,b′)=N ⁣([ψ(a,x) ψ(a′,x′)‾]x,x′∈A)b,b′.T_{(a,b),(a',b')}= N\!\left(\left[\psi(a,x)\,\overline{\psi(a',x')}\right]_{x,x'\in A}\right)_{b,b'}.T(a,b),(a′,b′)​=N([ψ(a,x)ψ(a′,x′)​]x,x′∈A​)b,b′​.

Thus NNN is applied separately to each block obtained by fixing the first AAA index on each side of the matrix [ψ(a,x)ψ(a′,x′)‾][\psi(a,x)\overline{\psi(a',x')}][ψ(a,x)ψ(a′,x′)​]. No linearity, positivity, or channel hypothesis is imposed on NNN.

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