Definitions for Channel simulation with quantum side information
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P42The objects the formalized statements of Channel simulation with quantum side information (arXiv:quant-ph/0611008) are stated in. Each was checked against the paper by a blind read-back and a moderator.
SimulationBasis: PDF p.4: classical blocks X^n,Y_A^n,Y_B^n and the retained quantum B^n.CQFeedbackCode: PDF p.4: one forward classical message, uniform independent shared bits, private stochastic encoding, a coin/message-controlled POVM, and one deterministic output function used by both parties. The POVM is implemented by its canonical square-root instrument, as required by the gentle measurement paragraph on p.12.feedbackMeasurementBranch: PDF pp.4,12, Eq.(6): reuse the accepted square-root instrument on a finite basis. The equivalence only relabels B as Fin (card B) and back; no coordinate or measurement outcome is discarded, and the branch remains unnormalized.simulationTarget: PDF pp.3--4, Eqs.(2),(4),(5): the memoryless feedback target, including X and B.simulationOutput: PDF p.4, Eq.(6): actual joint output, averaged over the independent uniform coin and Alice's stochastic record. Bob's record is measured, not supplied by Alice.feedbackTraceError: PDF p.4, Eq.(6): full trace norm, without the trace-distance factor 1/2.SimulationAchievable: PDF pp.4--5: operational achievability. Integer bit budgets replace nR,nC; rounding is covered by every positive rate slack. Codes precede the sampled input.SourceSimulation: PDF p.14, Eq.(36): source-relative resource simulation, with identities initially distributing X to Alice and B to Bob and a feedback W at the output. This is the same physical Eq.(6) task, not universal arbitrary-input simulation.bookProbability: PDF pp.6--9: probability of an event of a mutually independent q^n codebook. The codebook is fixed before a uniform index is sampled to define its empirical law.empiricalLaw: PDF p.6: the conditional law of G(W) for a fixed codebook, with W uniform on [M].dilutionError: PDF p.6: full classical l1 error; there is no factor 1/2.CoverEncoder: PDF p.7: a fully normalized stochastic encoder, including zero-marginal rows.coveringMarginal: PDF p.7: p(x) is derived from q(y)P(x|y), not a free marginal parameter.coveringError: PDF pp.7,9: full l1 comparison on [M] times X, not merely an output marginal.productConditional: PDF p.9: the conditional product channel P^n(x^n|y^n).RDDatum: PDF pp.15--16: finite-block variational datum, over all finite positive classical auxiliary alphabets and all message-controlled reconstruction POVMs.wordDistortion: PDF p.15: additive per-letter distortion; theorems only use positive n.rdExpected: PDF p.16: expected distortion includes the decoder's actual Born probability.blockRD: PDF p.16: finite-block variational value in extended nonnegative reals. Empty feasible sets have value infinity. Finite attainment is asserted below, not built into a hypothesis or lost by use of real sInf's empty-set default.ExactDistortionCode: PDF pp.15--16: a deterministic encoder and a message-controlled reconstruction POVM. Outcomes have already been coarsened to the reproduced word. No common randomness, entanglement, or additional communication is supplied.exactDistortionExpected: PDF p.16: literal finite-block expected distortion, including the actual Born probabilities of Bob's reconstruction; the source is not replaced by its mean.ExactDistortionAchievable: PDF p.16: the printed operational convention has rate slack but exact expected distortion at each sufficiently large block. It is kept here for the p.17 converse; no equality with Theorem 4.2's variational limit is asserted.exactDistortionRate: PDF p.16: infimum of literal achievable nonnegative rates. The empty set has value infinity, so an impossible distortion target cannot become a zero rate.QuantumRDChoice: PDF p.17: a finite-dimensional CQ channel X to C and a joint C B reconstruction POVM. The positive dimension is part of the datum and has NO upper bound. Every channel output is a density, even at source labels of weight zero.quantumRDExpected: PDF p.17: expected one-letter distortion. Given X=x the quantum systems C and B are in the PRODUCT state W(x) tensor rho_x; the decoder is joint on CB.quantumRDCost: PDF p.17: I(X;C)-I(C;B) in bits. The accepted label information and conditional-product average compute both information terms from the same source.quantumRDValue: PDF p.17: optimize over ALL finite positive auxiliary dimensions and all joint reconstruction POVMs. Nonnegativity follows from the conditional-product source; ofReal does not clip valid costs. Finite attainment is a theorem conclusion.redistributionReceiverInformation: PDF p.18: I(Bhat;B), computed in the same pure four-register source.senderMutual: PDF p.18: I(Bhat;Ahat), computed by actual partial traces.fqrsQuantum: PDF p.18: FQRS communication I(Bhat;RB)/2, without dropping Bob's register.printedQuantumCost: PDF p.18, Proposition 5.1: printed I(Bhat;R|Ahat)/2. For a pure source this equals the accepted quantumCost I(Bhat;R|B)/2 by complementary entropy duality; retaining the printed conditioning avoids silently changing the result.diagonalCQSource: PDF pp.17,19: classical side information as diagonal conditional densities. The channel Z(.|x) is normalized even on zero-probability source labels.CQFeedbackCode.IsClassical: PDF p.19: diagonal POVMs are classical stochastic decoding of the retained Z^n. Their square-root instruments leave each classical basis label intact.ClassicalSimulationAchievable: PDF p.19: fully classical version of the Eq.(6) task. Diagonal instruments preserve the original Z^n; X^n and both Y^n outputs are still in the error test.
import Mathlib
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P40
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P07
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P20
/-!
Luo--Devetak, arXiv:quant-ph/0611008v1. Entropies and resource rates are in bits.
The accepted paper modules imported here are frozen, reviewed library interfaces.
Distinct source-specific definitions are paper-local; the accepted Luders branch
is reused through finite-basis relabelling. The simulation relation is computed from an
encoder, a square-root measurement instrument and two feedback outputs; it is
not defined by its proposed entropy region. The two source-ambiguity deferrals and the
exact Theorem 4.1 reuse are recorded in the unit artifacts. Compilation pending.
-/
noncomputable section
namespace QAlgorithms.Papers.Independent.Intakethird100.P42
open scoped BigOperators ComplexConjugate ENNReal
open Filter
open QAlgorithms.Papers.Independent.Dependency50.P07
open QAlgorithms.Papers.Independent.Gap48.P20
open QAlgorithms.Papers.Independent.Information100.P27 (IsProbability shannon)
attribute [local instance] Classical.decEq Classical.propDecidable
/-- PDF p.4: classical blocks X^n,Y_A^n,Y_B^n and the retained quantum B^n. -/
abbrev SimulationBasis (X B : Type) (y n : ℕ) :=
Copies X n × (Copies (Fin y) n × (Copies (Fin y) n × Copies B n))
/-- PDF p.4: one forward classical message, uniform independent shared bits,
private stochastic encoding, a coin/message-controlled POVM, and one deterministic
output function used by both parties. The POVM is implemented by its canonical
square-root instrument, as required by the gentle measurement paragraph on p.12. -/
structure CQFeedbackCode (X B : Type) [Fintype B] (y n : ℕ) where
coinBits : ℕ
messageBits : ℕ
privateBits : ℕ
encode : Copies X n → Qubits coinBits → Qubits messageBits → Qubits privateBits → ℝ
encode_nonnegative : ∀ x l m s, 0 ≤ encode x l m s
encode_normalized : ∀ x l, (∑ m, ∑ s, encode x l m s) = 1
measure : Qubits coinBits → Qubits messageBits → Qubits privateBits →
Matrix (Copies B n) (Copies B n) ℂ
measure_valid : ∀ l m, IsPOVM (measure l m)
decode : Qubits coinBits → Qubits messageBits → Qubits privateBits → Copies (Fin y) n
/-- PDF pp.4,12, Eq.(6): reuse the accepted square-root instrument on a finite
basis. The equivalence only relabels B as Fin (card B) and back; no coordinate or
measurement outcome is discarded, and the branch remains unnormalized. -/
def feedbackMeasurementBranch {B : Type} [Fintype B]
(effect rho : Matrix B B ℂ) : Matrix B B ℂ :=
let e := Fintype.equivFin B
(QAlgorithms.Papers.Independent.Claimed50.P40.ludersBranch
(fun _ : Unit => effect.submatrix e.symm e.symm)
(rho.submatrix e.symm e.symm) ()).submatrix e e
/-- PDF pp.3--4, Eqs.(2),(4),(5): the memoryless feedback target, including X and B. -/
def simulationTarget {X B : Type} (E : CQData X B) (V : Auxiliary X) (n : ℕ) :
Matrix (SimulationBasis X B V.size n) (SimulationBasis X B V.size n) ℂ :=
fun a b =>
if a.1 = b.1 ∧ a.2.1 = b.2.1 ∧ a.2.2.1 = b.2.2.1 ∧ a.2.1 = a.2.2.1 then
((wordProb E.prob n a.1 * ∏ i, V.channel (a.1 i) (a.2.1 i)) : ℂ) *
wordState E n a.1 a.2.2.2 b.2.2.2
else 0
/-- PDF p.4, Eq.(6): actual joint output, averaged over the independent uniform
coin and Alice's stochastic record. Bob's record is measured, not supplied by Alice. -/
def simulationOutput {X B : Type} [Fintype B] (E : CQData X B) {y n : ℕ}
(P : CQFeedbackCode X B y n) :
Matrix (SimulationBasis X B y n) (SimulationBasis X B y n) ℂ :=
fun a b => if a.1 = b.1 ∧ a.2.1 = b.2.1 ∧ a.2.2.1 = b.2.2.1 then
∑ l : Qubits P.coinBits, ∑ m : Qubits P.messageBits,
∑ s : Qubits P.privateBits, ∑ t : Qubits P.privateBits,
if P.decode l m s = a.2.1 ∧ P.decode l m t = a.2.2.1 then
((wordProb E.prob n a.1 * ((2 : ℝ) ^ P.coinBits)⁻¹ *
P.encode a.1 l m s) : ℂ) *
feedbackMeasurementBranch (P.measure l m t) (wordState E n a.1) a.2.2.2 b.2.2.2
else 0
else 0
/-- PDF p.4, Eq.(6): full trace norm, without the trace-distance factor 1/2. -/
def feedbackTraceError {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (V : Auxiliary X) {n : ℕ} (P : CQFeedbackCode X B V.size n) : ℝ :=
traceNorm (simulationTarget E V n - simulationOutput E P)
/-- PDF pp.4--5: operational achievability. Integer bit budgets replace nR,nC;
rounding is covered by every positive rate slack. Codes precede the sampled input. -/
def SimulationAchievable {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (V : Auxiliary X) (R C : ℝ) : Prop :=
0 ≤ R ∧ 0 ≤ C ∧ ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
∀ n ≥ N, 1 ≤ n → ∃ P : CQFeedbackCode X B V.size n,
feedbackTraceError E V P ≤ eps ∧ (P.messageBits : ℝ) ≤ n * (R + delta) ∧
(P.coinBits : ℝ) ≤ n * (C + delta)
/-- PDF p.14, Eq.(36): source-relative resource simulation, with identities
initially distributing X to Alice and B to Bob and a feedback W at the output.
This is the same physical Eq.(6) task, not universal arbitrary-input simulation. -/
def SourceSimulation {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (V : Auxiliary X) (R C : ℝ) : Prop :=
SimulationAchievable E V R C
example {B : Type} [Fintype B] (rho : Matrix B B ℂ) :
feedbackMeasurementBranch 0 rho = 0 := by
simp [feedbackMeasurementBranch,
QAlgorithms.Papers.Independent.Claimed50.P40.ludersBranch, cfc_zero]
example {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (V : Auxiliary X) (R C : ℝ) :
SourceSimulation E V R C ↔ SimulationAchievable E V R C := Iff.rfl
/-- PDF pp.6--9: probability of an event of a mutually independent q^n codebook.
The codebook is fixed before a uniform index is sampled to define its empirical law. -/
def bookProbability {Y : Type} [Fintype Y] (q : Y → ℝ) (n M : ℕ)
(event : (Fin M → Copies Y n) → Prop) : ℝ :=
∑ book : Fin M → Copies Y n,
if event book then ∏ i, wordProb q n (book i) else 0
/-- PDF p.6: the conditional law of G(W) for a fixed codebook, with W uniform on [M]. -/
def empiricalLaw {Y : Type} [Fintype Y] {M : ℕ} (book : Fin M → Y) (y : Y) : ℝ :=
(∑ i, if book i = y then (1 : ℝ) else 0) / M
/-- PDF p.6: full classical l1 error; there is no factor 1/2. -/
def dilutionError {Y : Type} [Fintype Y] (q : Y → ℝ) (n M : ℕ)
(book : Fin M → Copies Y n) : ℝ :=
∑ y, |wordProb q n y - empiricalLaw book y|
/-- PDF p.7: a fully normalized stochastic encoder, including zero-marginal rows. -/
abbrev CoverEncoder (X : Type) [Fintype X] (M : ℕ) :=
{e : X → Fin M → ℝ // ∀ x, IsProbability (e x)}
/-- PDF p.7: p(x) is derived from q(y)P(x|y), not a free marginal parameter. -/
def coveringMarginal {X Y : Type} [Fintype Y] (q : Y → ℝ) (P : Y → X → ℝ) (x : X) : ℝ :=
∑ y, q y * P y x
/-- PDF pp.7,9: full l1 comparison on [M] times X, not merely an output marginal. -/
def coveringError {X Y : Type} [Fintype X] [Fintype Y]
(q : Y → ℝ) (P : Y → X → ℝ) (M : ℕ) (book : Fin M → Y)
(e : CoverEncoder X M) : ℝ :=
∑ i, ∑ x, |P (book i) x / M - coveringMarginal q P x * e.val x i|
/-- PDF p.9: the conditional product channel P^n(x^n|y^n). -/
def productConditional {X Y : Type} (P : Y → X → ℝ) (n : ℕ)
(y : Copies Y n) (x : Copies X n) : ℝ := ∏ i, P (y i) (x i)
example {Y : Type} [Fintype Y] (q : Y → ℝ) (n M : ℕ) :
bookProbability q n M (fun _ => False) = 0 := by
simp [bookProbability]
example {X Y : Type} (P : Y → X → ℝ) (x : Copies X 0) (y : Copies Y 0) :
productConditional P 0 y x = 1 := by
simp [productConditional]
example {Y : Type} [Fintype Y] [Subsingleton Y] [Nonempty Y]
(M : ℕ) (hM : 0 < M) (book : Fin M → Y) (y : Y) :
empiricalLaw book y = 1 := by
have hbook : ∀ i, book i = y := fun _ => Subsingleton.elim _ _
simp [empiricalLaw, hbook, Nat.ne_of_gt hM]
example {X : Type} [Fintype X] : Nonempty (CoverEncoder X 1) := by
refine ⟨⟨fun _ _ => 1, ?_⟩⟩
intro x
exact ⟨fun _ => by norm_num, by simp⟩
/-- PDF pp.15--16: finite-block variational datum, over all finite positive
classical auxiliary alphabets and all message-controlled reconstruction POVMs. -/
structure RDDatum (X B : Type) [Fintype X] [Fintype B] (n : ℕ) where
auxiliary : Auxiliary (Copies X n)
decoder : Fin auxiliary.size → Copies X n → Matrix (Copies B n) (Copies B n) ℂ
decoder_valid : ∀ y, IsPOVM (decoder y)
/-- PDF p.15: additive per-letter distortion; theorems only use positive n. -/
def wordDistortion {X : Type} (d : X → X → ℝ) (n : ℕ)
(x z : Copies X n) : ℝ := (∑ i, d (x i) (z i)) / n
/-- PDF p.16: expected distortion includes the decoder's actual Born probability. -/
def rdExpected {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
(d : X → X → ℝ) {n : ℕ} (P : RDDatum X B n) : ℝ :=
∑ x, ∑ y, ∑ z, wordProb E.prob n x * P.auxiliary.channel x y *
effectProb (P.decoder y z) (wordState E n x) * wordDistortion d n x z
/-- PDF p.16: finite-block variational value in extended nonnegative reals.
Empty feasible sets have value infinity. Finite attainment is asserted below,
not built into a hypothesis or lost by use of real sInf's empty-set default. -/
def blockRD {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
(d : X → X → ℝ) (D : ℝ) (n : ℕ) : ℝ≥0∞ :=
sInf {r | ∃ P : RDDatum X B n, rdExpected E d P ≤ D ∧
r = ENNReal.ofReal (auxCost (iidCQ E n) P.auxiliary / n)}
example {X : Type} (d : X → X → ℝ) (x z : Copies X 0) :
wordDistortion d 0 x z = 0 := by simp [wordDistortion]
/-- PDF pp.15--16: a deterministic encoder and a message-controlled reconstruction
POVM. Outcomes have already been coarsened to the reproduced word. No common
randomness, entanglement, or additional communication is supplied. -/
structure ExactDistortionCode (X B : Type) [Fintype X] [Fintype B] (n : ℕ) where
messageBits : ℕ
encode : Copies X n → Qubits messageBits
decoder : Qubits messageBits → Copies X n → Matrix (Copies B n) (Copies B n) ℂ
decoder_valid : ∀ m, IsPOVM (decoder m)
/-- PDF p.16: literal finite-block expected distortion, including the actual
Born probabilities of Bob's reconstruction; the source is not replaced by its mean. -/
def exactDistortionExpected {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (d : X → X → ℝ) {n : ℕ} (P : ExactDistortionCode X B n) : ℝ :=
∑ x, ∑ z, wordProb E.prob n x *
effectProb (P.decoder (P.encode x) z) (wordState E n x) * wordDistortion d n x z
/-- PDF p.16: the printed operational convention has rate slack but exact
expected distortion at each sufficiently large block. It is kept here for the
p.17 converse; no equality with Theorem 4.2's variational limit is asserted. -/
def ExactDistortionAchievable {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (d : X → X → ℝ) (D R : ℝ) : Prop :=
0 ≤ R ∧ ∀ delta > (0 : ℝ), ∃ N : ℕ, ∀ n ≥ N, 1 ≤ n →
∃ P : ExactDistortionCode X B n,
exactDistortionExpected E d P ≤ D ∧ (P.messageBits : ℝ) ≤ n * (R + delta)
/-- PDF p.16: infimum of literal achievable nonnegative rates. The empty set
has value infinity, so an impossible distortion target cannot become a zero rate. -/
def exactDistortionRate {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (d : X → X → ℝ) (D : ℝ) : ℝ≥0∞ :=
sInf {r | ∃ R : ℝ, ExactDistortionAchievable E d D R ∧ r = ENNReal.ofReal R}
/-- PDF p.17: a finite-dimensional CQ channel X to C and a joint C B
reconstruction POVM. The positive dimension is part of the datum and has NO
upper bound. Every channel output is a density, even at source labels of weight zero. -/
structure QuantumRDChoice (X B : Type) [Fintype X] [Fintype B] where
dimension : ℕ
positive : 0 < dimension
output : X → Matrix (Fin dimension) (Fin dimension) ℂ
output_valid : ∀ x, IsDensityMatrix (output x)
decoder : X → Matrix (Fin dimension × B) (Fin dimension × B) ℂ
decoder_valid : IsPOVM decoder
/-- PDF p.17: expected one-letter distortion. Given X=x the quantum systems
C and B are in the PRODUCT state W(x) tensor rho_x; the decoder is joint on CB. -/
def quantumRDExpected {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (d : X → X → ℝ) (P : QuantumRDChoice X B) : ℝ :=
∑ x, ∑ z, E.prob x *
effectProb (P.decoder z) (Matrix.kronecker (P.output x) (E.state x)) * d x z
/-- PDF p.17: I(X;C)-I(C;B) in bits. The accepted label information and
conditional-product average compute both information terms from the same source. -/
def quantumRDCost {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (P : QuantumRDChoice X B) : ℝ :=
QAlgorithms.Papers.Independent.Information100.P27.labelInformation E.prob P.output -
(vnEntropy (QAlgorithms.Papers.Independent.Information100.P27.average E.prob P.output) +
vnEntropy (ensembleMean E) -
vnEntropy (QAlgorithms.Papers.Independent.Information100.P27.memoryEnvironment
E.prob E.state P.output))
/-- PDF p.17: optimize over ALL finite positive auxiliary dimensions and all
joint reconstruction POVMs. Nonnegativity follows from the conditional-product
source; ofReal does not clip valid costs. Finite attainment is a theorem conclusion. -/
def quantumRDValue {X B : Type} [Fintype X] [Fintype B]
(E : CQData X B) (d : X → X → ℝ) (D : ℝ) : ℝ≥0∞ :=
sInf {r | ∃ P : QuantumRDChoice X B, quantumRDExpected E d P ≤ D ∧
r = ENNReal.ofReal (quantumRDCost E P)}
example {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
(P : ExactDistortionCode X B 0) (d : X → X → ℝ) :
exactDistortionExpected E d P = 0 := by
simp [exactDistortionExpected, wordDistortion]
example {X B : Type} [Fintype X] [Fintype B] (E : CQData X B)
(P : QuantumRDChoice X B) : quantumRDExpected E (fun _ _ => 0) P = 0 := by
simp [quantumRDExpected]
/-- PDF p.18: I(Bhat;B), computed in the same pure four-register source. -/
def redistributionReceiverInformation {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
(psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
(fun p => (p.1, p.2.1, p.2.2.2)) +
subsystemEntropy (pureDensity psi) (fun p => p.2.2.2)
(fun p => (p.1, p.2.1, p.2.2.1)) -
subsystemEntropy (pureDensity psi) (fun p => p.2.2) (fun p => (p.1, p.2.1))
/-- PDF p.18: I(Bhat;Ahat), computed by actual partial traces. -/
def senderMutual {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
(psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
(fun p => (p.1, p.2.1, p.2.2.2)) +
subsystemEntropy (pureDensity psi) (fun p => p.2.1) (fun p => (p.1, p.2.2)) -
subsystemEntropy (pureDensity psi) (fun p => (p.2.1, p.2.2.1))
(fun p => (p.1, p.2.2.2))
/-- PDF p.18: FQRS communication I(Bhat;RB)/2, without dropping Bob's register. -/
def fqrsQuantum {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
(psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
(subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
(fun p => (p.1, p.2.1, p.2.2.2)) +
subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2.2))
(fun p => (p.2.1, p.2.2.1)) -
subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2))
(fun p => p.2.1)) / 2
/-- PDF p.18, Proposition 5.1: printed I(Bhat;R|Ahat)/2. For a pure source
this equals the accepted quantumCost I(Bhat;R|B)/2 by complementary entropy
duality; retaining the printed conditioning avoids silently changing the result. -/
def printedQuantumCost {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
(psi : EuclideanSpace ℂ (RedistributionSource A C B R)) : ℝ :=
(subsystemEntropy (pureDensity psi) (fun p => (p.2.1, p.2.2.1))
(fun p => (p.1, p.2.2.2)) +
subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.1))
(fun p => p.2.2) -
subsystemEntropy (pureDensity psi) (fun p => p.2.1)
(fun p => (p.1, p.2.2)) -
subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.1, p.2.2.1))
(fun p => p.2.2.2)) / 2
example {A C B R : Type} [Fintype A] [Fintype C] [Fintype B] [Fintype R]
(psi : EuclideanSpace ℂ (RedistributionSource A C B R)) :
fqrsQuantum psi = (subsystemEntropy (pureDensity psi) (fun p => p.2.2.1)
(fun p => (p.1, p.2.1, p.2.2.2)) +
subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2.2))
(fun p => (p.2.1, p.2.2.1)) -
subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2))
(fun p => p.2.1)) / 2 := rfl
/-- PDF pp.17,19: classical side information as diagonal conditional densities.
The channel Z(.|x) is normalized even on zero-probability source labels. -/
def diagonalCQSource {X : Type} (p : X → ℝ) (Z : Auxiliary X) : CQData X (Fin Z.size) where
prob := p
state := fun x => Matrix.diagonal (fun z => (Z.channel x z : ℂ))
/-- PDF p.19: diagonal POVMs are classical stochastic decoding of the retained
Z^n. Their square-root instruments leave each classical basis label intact. -/
def CQFeedbackCode.IsClassical {X B : Type} [Fintype B] {y n : ℕ}
(P : CQFeedbackCode X B y n) : Prop :=
∀ l m s u v, u ≠ v → P.measure l m s u v = 0
/-- PDF p.19: fully classical version of the Eq.(6) task. Diagonal instruments
preserve the original Z^n; X^n and both Y^n outputs are still in the error test. -/
def ClassicalSimulationAchievable {X : Type} [Fintype X] (p : X → ℝ)
(Z V : Auxiliary X) (R C : ℝ) : Prop :=
0 ≤ R ∧ 0 ≤ C ∧ ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
∀ n ≥ N, 1 ≤ n → ∃ P : CQFeedbackCode X (Fin Z.size) V.size n,
P.IsClassical ∧ feedbackTraceError (diagonalCQSource p Z) V P ≤ eps ∧
(P.messageBits : ℝ) ≤ n * (R + delta) ∧ (P.coinBits : ℝ) ≤ n * (C + delta)
example {X : Type} (p : X → ℝ) (Z : Auxiliary X) (x : X) (z : Fin Z.size) :
(diagonalCQSource p Z).state x z z = (Z.channel x z : ℂ) := by
simp [diagonalCQSource]
end QAlgorithms.Papers.Independent.Intakethird100.P42
Read-back
What the Lean code literally says, in plain math · gpt-6.1-sol
QAlgorithms.Papers.Independent.Intakethird100.P42.CQFeedbackCode.IsClassical
Let be any type, let be a finite type, and let be arbitrary natural numbers, including zero. Let consist of natural numbers , a function for , , , and , satisfying for every ; matrices on , positive semidefinite and satisfying for every ; and a decoding function . The predicate defined here is
Thus it requires every one of the specified measurement matrices to be diagonal. When , has one element and this off-diagonal condition is vacuous. When , no such exists because its decoding function would have a nonempty domain and an empty codomain.
QAlgorithms.Papers.Independent.Intakethird100.P42.CQFeedbackCode
For an arbitrary type , a finite type with a specified finite enumeration, and natural numbers , a code consists of natural numbers , called its coin, message, and private bit counts; a real function , where , , , and , satisfying
matrices on , such that for every , every is positive semidefinite and ; and a function . Here , and words are functions on . The type need not be finite, and may be empty. Zero bit counts are allowed and give singleton bitstring sets. Length is allowed and gives singleton word sets; if , no such code exists because the required decoder would map a nonempty bitstring domain into an empty word type.
QAlgorithms.Papers.Independent.Intakethird100.P42.ClassicalSimulationAchievable
Let be a finite type, possibly empty, let be arbitrary, and let be auxiliary data with positive finite output sizes , respectively, and conditional probabilities , , each nonnegative and summing to over its output for every . For real numbers , this predicate means , , and that for every there is a natural such that every with admits the following data: natural numbers ; encoder weights , for , , , , satisfying ; diagonal positive semidefinite matrices on , satisfying for every ; and a decoding function . To specify the required error, put , , and . On the basis indexed by , let and be diagonal matrices with entries
The required bounds are , , and , where . No bound is imposed on , and may be zero. No sign or normalization condition is imposed on ; the displayed matrices and error are used even when is not a probability law.
QAlgorithms.Papers.Independent.Intakethird100.P42.CoverEncoder
For a finite type and a natural number , a cover encoder is a function satisfying, for every , for every and . No nonemptiness of or positivity of is required in the definition: if is empty these conditions are vacuous, whereas if is nonempty and , no cover encoder exists.
QAlgorithms.Papers.Independent.Intakethird100.P42.ExactDistortionAchievable
For finite, finitely enumerated types , arbitrary data with and complex matrices on , any real-valued function , and real numbers , this property means that and that, for every real , there exists a natural number such that every natural number with admits a natural number , an encoder , and matrices on , positive semidefinite and summing to the identity over for each message , satisfying
where . The code may depend on and ; the distortion bound is exactly . There are no validity assumptions on , no sign assumptions on , and no nonemptiness assumptions on . The block length is excluded from the code requirement, while zero message bits are allowed.
QAlgorithms.Papers.Independent.Intakethird100.P42.ExactDistortionCode
For finite types and , each equipped with a finite enumeration, and any natural number , an exact distortion code consists of a natural number of message bits, a function , and matrices on for every message and reconstruction , where and mean functions from into and , respectively. For every message , the matrices satisfy for all and ; this requirement includes messages outside the image of . Neither nor must be positive, and neither type must be nonempty. At , each word type has one element, including when its alphabet is empty.
QAlgorithms.Papers.Independent.Intakethird100.P42.QuantumRDChoice
For finite types , each equipped with a finite enumeration, a quantum rate-distortion choice consists of a positive natural number , matrices on for every , and matrices on for every , satisfying
Neither nor is required to be nonempty. If is empty, the output-state conditions are vacuous, but the decoder condition still requires its empty sum to equal the identity; consequently, when is nonempty, no such choice exists for empty .
QAlgorithms.Papers.Independent.Intakethird100.P42.RDDatum
For finite types and a natural number , a rate-distortion datum consists of a positive natural number , a real function satisfying and for every word , and matrices on , indexed by and , such that for each , every is positive semidefinite and . Here and are the functions from into and , respectively. The definition allows and empty types , subject to these normalization and matrix conditions; for example, when is empty and is nonempty, the required empty sum cannot equal the identity, so no such datum exists.
QAlgorithms.Papers.Independent.Intakethird100.P42.SimulationAchievable
For finite types with specified finite enumerations, arbitrary data , an auxiliary system of positive size with and , and real numbers , this property means that and, for every , there is such that every with admits a code with coin, message, and private bit counts satisfying
The code has nonnegative weights summing to over for every , decoded words , and matrices on summing to over for every ; are bitstrings of lengths . With and , the matrices on tuples have entries
Here relabels onto a finite ordinal, forms , and relabels back; . The private bit count is unrestricted. No validity assumption on the ensemble data or nonemptiness assumption on either finite type is part of this definition; only positive block lengths are tested.
QAlgorithms.Papers.Independent.Intakethird100.P42.SimulationBasis
For arbitrary types and natural numbers , this is the type
where and denotes functions from to . Thus an element consists of an -word, two separate -words, and a -word, all of length . No finiteness or nonemptiness assumptions on or are imposed. When , every word type here is a singleton, even if its alphabet is empty; when , an empty alphabet makes the corresponding word type empty.
QAlgorithms.Papers.Independent.Intakethird100.P42.SourceSimulation
For finite types with specified finite enumerations, arbitrary real weights , complex matrices on , an auxiliary system of positive size with probabilities , and real numbers , this is the following property: , and for every , all sufficiently large positive integers admit natural bit counts , encoding weights summing to over , a decoder , and matrices on summing to over , such that , , and . Here , the bitstrings have lengths , , , and
The matrix indices range over ; relabels onto a finite ordinal, forms , and relabels back, while . There is no bound on , no ensemble-validity hypothesis, and no nonemptiness hypothesis on or .
QAlgorithms.Papers.Independent.Intakethird100.P42.blockRD
For finite types , arbitrary data with and complex matrices on , an arbitrary real function , a real threshold , and a natural number , this definition takes an infimum in over the following quantities. Choose any , nonnegative conditional weights with for every , and positive semidefinite matrices on with , subject to
where and . For each such choice set , , , and . The quantity entered into the infimum is
embedded as a finite extended nonnegative real, where and for , evaluated by continuous functional calculus. No probability or density assumptions on , or sign assumptions on , are imposed. All divisions are total, so if , and the displayed quotient is when . The infimum of an empty feasible set is ; at , every feasible datum contributes .
QAlgorithms.Papers.Independent.Intakethird100.P42.bookProbability
For a finite type with a specified finite enumeration, an arbitrary real function , natural numbers , and any predicate on functions , this defines the real number
There is no assumption that is nonnegative or sums to . Both and are included: in either case there is exactly one book and its product weight is , so the result is or according to the predicate on that book. If is empty and both are positive, there are no books and the sum is .
QAlgorithms.Papers.Independent.Intakethird100.P42.coveringError
For finite types , arbitrary real functions and , a natural number , an indexed book , and a function with for every , the covering error is
No probability assumptions on or are imposed. Division by has value , and the outer sum is empty when , giving error whenever the supplied encoder exists; such an encoder cannot exist for and nonempty .
QAlgorithms.Papers.Independent.Intakethird100.P42.coveringMarginal
For arbitrary types with finite, functions and , and , the covering marginal at is . Neither function is required to be nonnegative or normalized, and need not be finite. If is empty, the value is .
QAlgorithms.Papers.Independent.Intakethird100.P42.diagonalCQSource
For any type , any function , and auxiliary data consisting of a positive natural number and numbers , indexed by and , with for every , this definition constructs data with label type , matrix index type , weight function , and matrices
viewing the real entries as complex numbers. There is no finiteness or nonemptiness assumption on , and no nonnegativity or normalization assumption on .
QAlgorithms.Papers.Independent.Intakethird100.P42.dilutionError
For any finite type , any function , natural numbers , and list of indexed words , where denotes the functions , the dilution error is
No nonnegativity or normalization of is assumed. Empty products are , empty sums are , and division by is defined as ; thus the empirical term is when . When , there is exactly one empty word, and the error is for and for .
QAlgorithms.Papers.Independent.Intakethird100.P42.empiricalLaw
For a finite type with a specified finite enumeration, any natural number , any function , and any , this defines
with numerator and denominator regarded as real numbers. The definition includes : the numerator is then zero, and real division by zero is defined to be zero, so the value is for every . No positive-size assumption is imposed.
QAlgorithms.Papers.Independent.Intakethird100.P42.exactDistortionExpected
For finite, finitely enumerated types , arbitrary data consisting of a real function and complex matrices on , an arbitrary function , a natural number , and an exact distortion code , the defined real number is
Here are length- word types, has some number of message bits, , and, for every message , and on . No normalization or positivity of or , and no sign condition on , is assumed. Empty alphabets and are allowed; at , empty products equal , and the distortion factor is , so the defined value is .
QAlgorithms.Papers.Independent.Intakethird100.P42.exactDistortionRate
For finite, finitely enumerated types , arbitrary data with real weights and complex matrices on , any , and any real , the defined value is the infimum in of all nonnegative real rates having the following property: for every , all sufficiently large positive integers admit an encoder , for some integer , and matrices on with and for every message , such that
More precisely, each such enters the infimum through its nonnegative extended-real image, which equals because . If no such rate exists, the value is ; attainment is not part of this definition. Empty alphabets are permitted, and no distribution, density-matrix, or nonnegativity assumptions on are imposed.
QAlgorithms.Papers.Independent.Intakethird100.P42.feedbackMeasurementBranch
For a finite type with a specified finite enumeration and arbitrary complex matrices indexed by , choose the enumeration equivalence , and write and . The resulting matrix has entries
where is the continuous functional calculus and the square root is the real square-root function, which is zero on negative inputs. No positivity, Hermiticity, trace, or normalization hypothesis is imposed on either matrix; the functional-calculus expression uses its total definition even outside its usual hypotheses. The empty finite type is also allowed.
QAlgorithms.Papers.Independent.Intakethird100.P42.feedbackTraceError
For finite types with specified finite enumerations, arbitrary ensemble data , an auxiliary system with positive size , probabilities summing to over , a natural number , and a code, this is
Here both matrices are indexed by . Write and . The target entry is zero unless , and otherwise is . To define , the code supplies bit counts , nonnegative encoding weights with , decoded words , and positive semidefinite matrices on with . Its entry is zero unless , and otherwise equals
where are bitstrings of lengths , and is obtained by relabeling onto a finite ordinal, applying the functional-calculus expression , and relabeling back. No probability or density-matrix assumptions on , and no nonemptiness assumptions on , are imposed. The definition includes , when products are empty products and word sets are singleton.
QAlgorithms.Papers.Independent.Intakethird100.P42.fqrsQuantum
For arbitrary finite types , possibly empty, and any vector , form . For a register collection , obtain by summing over equal complementary labels, so that , with coordinates in their original order. Define , where is applied by continuous functional calculus and . The assigned real number is exactly
No normalization or nonzero condition is imposed on , and the definition also applies when a register type is empty.
QAlgorithms.Papers.Independent.Intakethird100.P42.printedQuantumCost
For arbitrary finite types , possibly empty, and any vector , let . For any collection of registers, let , summing over the complementary labels with coordinates restored to their original order, and put , where is applied by continuous functional calculus, with . The definition assigns
This formula requires neither nor , and includes empty register types.
QAlgorithms.Papers.Independent.Intakethird100.P42.productConditional
For arbitrary types , any real function , any natural number , and words and , the product conditional is . No finiteness, nonemptiness, nonnegativity, or normalization assumptions are imposed on . For , the product is .
QAlgorithms.Papers.Independent.Intakethird100.P42.quantumRDCost
For finite, finitely enumerated types , arbitrary data with and complex matrices on , and a quantum rate-distortion choice , let be its output dimension and let be its positive semidefinite, trace-one output matrices on . Write
where is evaluated by continuous functional calculus. The defined cost is
The choice also includes matrices on with , although these matrices do not occur in the cost. No validity assumptions on or nonemptiness assumptions on are made, and no sign condition on the resulting cost is included.
QAlgorithms.Papers.Independent.Intakethird100.P42.quantumRDExpected
For finite, finitely enumerated types , arbitrary data consisting of and complex matrices on , any function , and a quantum rate-distortion choice , the defined real number is
Here supplies an integer , positive semidefinite matrices on with trace , and positive semidefinite matrices on satisfying . The tensor product has entries . No probability or density-matrix conditions are imposed on , and may take negative values. Empty or is permitted whenever such a choice exists; if is empty, the displayed sum is .
QAlgorithms.Papers.Independent.Intakethird100.P42.quantumRDValue
For finite, finitely enumerated types , arbitrary data with real weights and complex matrices on , any , and any real , consider all choices of an integer , matrices on with , and matrices on with , for which
For each choice, set , , , and
using continuous functional calculus. The defined value is the infimum in of the numbers from these choices. Thus negative costs are mapped to , and an empty set of admissible choices gives . Neither attainment, validity of , nonnegativity of , nor nonemptiness of is assumed.
QAlgorithms.Papers.Independent.Intakethird100.P42.rdExpected
For finite types , data consisting of an arbitrary real function and arbitrary complex matrices on , an arbitrary function , a natural number , and a rate-distortion datum consisting of , nonnegative normalized conditional weights for , and positive semidefinite matrices on satisfying , the expected-distortion expression is
No validity assumption on is imposed: need not be a probability function and the matrices need not be density matrices. The function may take arbitrary real values. Empty products are , empty sums are , and division by is , so this expression is at .
QAlgorithms.Papers.Independent.Intakethird100.P42.redistributionReceiverInformation
For arbitrary finite types , possibly empty, and any vector , form the matrix . For a collection of these registers, let be the matrix obtained by summing 's entries over equal labels in all complementary registers; equivalently, , with the coordinates placed in their original order. Write , where is applied by continuous functional calculus, with . The definition assigns the real number
There is no assumption that ; the same formula applies to the zero vector and to empty register types.
QAlgorithms.Papers.Independent.Intakethird100.P42.senderMutual
For arbitrary finite types , possibly empty, and any vector , set . For each collection of registers, define its reduced matrix by , summing over all complementary register labels and restoring the original coordinate order. Let , with applied by continuous functional calculus and . The definition assigns
The vector need not be normalized, and zero vectors and empty register types are included.
QAlgorithms.Papers.Independent.Intakethird100.P42.simulationOutput
For an arbitrary type , a finite type with a specified finite enumeration, arbitrary real weights and complex matrices on , natural numbers , and a code, let its bit counts be , its encoding numbers be with , its decoder be , and its matrices on satisfy and for every ; the indices range respectively over bitstrings of lengths . Put and . For matrices on , let mean: relabel both matrices along the specified equivalence , form in those coordinates, and relabel back. The defined output matrix is
Thus supplies the encoding weight and the first decoded word, while the independently summed supplies the measurement matrix and the second decoded word; equality of those two words is not required. No validity hypothesis is imposed on the ensemble data. Zero lengths and zero bit counts are permitted, with empty products equal to and zero-length bitstring sets singleton.
QAlgorithms.Papers.Independent.Intakethird100.P42.simulationTarget
For arbitrary types , data consisting of a real function and complex matrices on , an auxiliary system consisting of a positive natural number and numbers with for every , and any natural number , this defines a matrix indexed by tuples . Its entries are
Real coefficients in this expression are regarded as complex numbers. The functions are not assumed to form valid ensemble data: neither probability normalization nor density-matrix conditions are required. Neither nor is assumed finite or nonempty. Empty products equal , so at this is the one-by-one matrix with entry .
QAlgorithms.Papers.Independent.Intakethird100.P42.wordDistortion
For any type , any function , any natural number , and words , the word distortion is . No sign, symmetry, diagonal, or metric conditions are imposed on , and need not be finite. When , both the empty sum and division by yield .
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.