Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for A triangle of dualities: reversibly decomposable quantum channels, source-channel duality, and time reversal

Definition
QAlgorithms_Papers_Independent_Intakethird100_P64

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

prove2quantum

The objects the formalized statements of A triangle of dualities: reversibly decomposable quantum channels, source-channel duality, and time reversal (arXiv:quant-ph/0505138) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • feedbackOutputVector: PDF p.2, (6), and p.3: the output vector is computed from U and phi.
  • feedbackEntropyA: PDF p.2, (8): H(A) for the actual pure tripartite vector.
  • mutualRA: PDF p.2: I(R;A)=H(R)+H(A)-H(RA), without using pure-state identities.
  • mutualAB: PDF p.2: I(A;B)=H(A)+H(B)-H(AB).
  • mutualRB: PDF p.2: I(R;B)=H(R)+H(B)-H(RB).
  • localOperation: PDF pp.2–4: a local operation from a retained isometric dilation. W is local work, discarded only when evaluating the reduced output. V†V=1 is carried by codes. Arbitrary work dimension permits typical-subspace failure flags on the full input.
  • FeedbackSimulationCode: PDF p.2, (5): q transmitted qubits and k supplied Bell pairs; A stays at Alice, B goes to Bob. Local work starts pure through the isometry description. No entropy inequalities or desired performance properties are fields of this code.
  • appendBell: PDF p.2: append exactly k normalized, independently supplied Bell pairs.
  • FeedbackSimulationCode.channel: PDF p.2, (5): the actual simulated channel, before extending by any reference.
  • isometryPower: PDF p.2: U applied independently n times, grouping the A and B outputs.
  • RelativeSimulationAchievable: PDF pp.1–2: FQRS simulation of the relative channel at input rho. For every error and rate slack, a strictly positive marginal tolerance is chosen. All large blocklengths have codes fixed before every reference/input test. Fin d ranges over all finite reference dimensions. The promise constrains the input marginal only; the error tests the entire correlated output. Rate zero allows sublinear resources.
  • SlepianWolfCode: PDF pp.3–4, (12), Fig.4: Alice receives A^n and Bob B^n from the Source. The q-qubit message is sent forward. Alice keeps l Bell halves; Bob outputs the other l halves and S^n, the copy of the original source. Local work remains in the dilations, so no impossible full-input compression isometry is assumed.
  • sourceDensity: PDF pp.3–4: the normalized source purification in grouped copy registers.
  • distributedSource: PDF p.3: apply only the specified source isometry; Alice and Bob cannot encode or choose the S input. The reference is untouched.
  • sourceTarget: PDF pp.3–4: source/reference and generated Bell pairs are tested jointly. The ideal S^n output is a copy at Bob, not a free matrix unrelated to phi.
  • SlepianWolfCode.error: PDF pp.3–4: full joint trace-norm error of the physical FQSW code.
  • SourceTransferAchievable: PDF pp.1,3–4: source transfer with Q forward qubits consumed and E ebits generated per source copy. No initial Bell pairs or other nonlocal assistance. The integer output is rounded down and clipped at zero only when E-delta<0.
  • sourceToAlice: PDF p.3: the identity from the Source to Alice, with Bob's initial share one-dimensional. S, A-hat, A-prime and B-hat are copies of the same basis type.
  • RelativeIdentityTest: PDF pp.2–3: the promise supplied by a relative n-use identity channel. It guarantees a joint reference test for every input with eta-close marginal.
  • RoutingCode: PDF p.3, Fig.3: local encoding, n uses of the supplied relative identity, and local decoding. The middle channel is supplied, not chosen by the code.
  • SourceRoutingAchievable: PDF p.3: direct source/identity concatenation, with the actual encoder marginal tested against the supplied channel's promise. The quantified F represents any physical implementation satisfying that promise, not unrestricted access. The code uses exactly n supplied identity calls for n source copies; no nonlocal ancillas or extra communication are available. delta records allowed rate slack.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_Ensemble
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P20
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P48

/-!
Devetak, arXiv:quant-ph/0505138v1, PDF pp.1–4.
Only the five results with complete task-specific operational semantics are stated.
Other nodes, including the printed capacity region, are explicitly deferred in HARD.md.
All register bases are finite and nonempty; Unit represents an absent register.
Entropy is in bits. Error is full trace norm, an equivalent vanishing-error convention
to half trace norm. Maps are selected before all test inputs. One forward quantum
message is counted by its actual number of qubits; Bell counters are actual pairs.
Source inputs are supplied by the Source and are never selected by a code.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P64

open QAlgorithms
open QAlgorithms.Papers.Independent.Gap48.P20
  (Copies redistributionMarginal subsystemEntropy forwardApply)
open QAlgorithms.Papers.Independent.Dependency50.P48 (densityPower)
open scoped Matrix ComplexConjugate

attribute [local instance] Classical.propDecidable Classical.decEq

/-- PDF p.2, (6), and p.3: the output vector is computed from U and phi. -/
def feedbackOutputVector {R S A B : Type} [Fintype S]
    (U : Matrix (A × B) S ℂ) (phi : EuclideanSpace ℂ (R × S)) :
    EuclideanSpace ℂ (R × (A × B)) :=
  WithLp.toLp 2 fun p => ∑ s, U p.2 s * phi (p.1, s)

section Entropies

variable {R A B : Type} [Fintype R] [Fintype A] [Fintype B]

/-- PDF p.2, (8): H(A) for the actual pure tripartite vector. -/
def feedbackEntropyA (psi : EuclideanSpace ℂ (R × (A × B))) : ℝ :=
  subsystemEntropy (pureDensity psi) (fun p => p.2.1) (fun p => (p.1, p.2.2))

/-- PDF p.2: I(R;A)=H(R)+H(A)-H(RA), without using pure-state identities. -/
def mutualRA (psi : EuclideanSpace ℂ (R × (A × B))) : ℝ :=
  subsystemEntropy (pureDensity psi) Prod.fst Prod.snd + feedbackEntropyA psi -
    subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.1)) (fun p => p.2.2)

/-- PDF p.2: I(A;B)=H(A)+H(B)-H(AB). -/
def mutualAB (psi : EuclideanSpace ℂ (R × (A × B))) : ℝ :=
  feedbackEntropyA psi +
    subsystemEntropy (pureDensity psi) (fun p => p.2.2) (fun p => (p.1, p.2.1)) -
    subsystemEntropy (pureDensity psi) Prod.snd Prod.fst

/-- PDF p.2: I(R;B)=H(R)+H(B)-H(RB). -/
def mutualRB (psi : EuclideanSpace ℂ (R × (A × B))) : ℝ :=
  subsystemEntropy (pureDensity psi) Prod.fst Prod.snd +
    subsystemEntropy (pureDensity psi) (fun p => p.2.2) (fun p => (p.1, p.2.1)) -
    subsystemEntropy (pureDensity psi) (fun p => (p.1, p.2.2)) (fun p => p.2.1)

end Entropies

/-- PDF pp.2–4: a local operation from a retained isometric dilation. W is local
work, discarded only when evaluating the reduced output. V†V=1 is carried by codes.
Arbitrary work dimension permits typical-subspace failure flags on the full input. -/
def localOperation {X Y W : Type} [Fintype X] [Fintype W]
    (V : Matrix (Y × W) X ℂ) : Superop X Y :=
  fun rho y y' => ∑ w, (conjSuperop V rho) (y, w) (y', w)

/-- PDF p.2, (5): q transmitted qubits and k supplied Bell pairs; A stays at
Alice, B goes to Bob. Local work starts pure through the isometry description.
No entropy inequalities or desired performance properties are fields of this code. -/
structure FeedbackSimulationCode (S A B : Type) [Fintype S] [Fintype A] [Fintype B] (n : ℕ) where
  q : ℕ
  k : ℕ
  aliceWork : ℕ
  bobWork : ℕ
  enc : Matrix (((Copies A n × Qubits q) × Fin (aliceWork + 1)))
    (Copies S n × Qubits k) ℂ
  dec : Matrix (Copies B n × Fin (bobWork + 1)) (Qubits q × Qubits k) ℂ
  enc_isometry : encᴴ * enc = 1
  dec_isometry : decᴴ * dec = 1

/-- PDF p.2: append exactly k normalized, independently supplied Bell pairs. -/
def appendBell {X : Type} (k : ℕ) (rho : Matrix X X ℂ) :
    Matrix ((X × Qubits k) × Qubits k) ((X × Qubits k) × Qubits k) ℂ :=
  fun p p' => rho p.1.1 p'.1.1 *
    pureDensity (maxEntangled k) (p.1.2, p.2) (p'.1.2, p'.2)

/-- PDF p.2, (5): the actual simulated channel, before extending by any reference. -/
def FeedbackSimulationCode.channel {S A B : Type} [Fintype S] [Fintype A] [Fintype B]
    {n : ℕ} (P : FeedbackSimulationCode S A B n) :
    Superop (Copies S n) (Copies A n × Copies B n) :=
  fun rho => forwardApply (localOperation P.enc) (localOperation P.dec) (appendBell P.k rho)

/-- PDF p.2: U applied independently n times, grouping the A and B outputs. -/
def isometryPower {S A B : Type} (U : Matrix (A × B) S ℂ) (n : ℕ) :
    Matrix (Copies A n × Copies B n) (Copies S n) ℂ :=
  fun p s => ∏ i, U (p.1 i, p.2 i) (s i)

/-- PDF pp.1–2: FQRS simulation of the relative channel at input rho. For every
error and rate slack, a strictly positive marginal tolerance is chosen. All large
blocklengths have codes fixed before every reference/input test. Fin d ranges over
all finite reference dimensions. The promise constrains the input marginal only;
the error tests the entire correlated output. Rate zero allows sublinear resources. -/
def RelativeSimulationAchievable {S A B : Type} [Fintype S] [Fintype A] [Fintype B]
    (U : Matrix (A × B) S ℂ) (rho : Matrix S S ℂ) (Q E : ℝ) : Prop :=
  ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ eta > (0 : ℝ), ∃ N : ℕ,
    ∀ n ≥ N, 1 ≤ n → ∃ P : FeedbackSimulationCode S A B n,
      P.q ≤ ⌊(n : ℝ) * (Q + delta)⌋₊ ∧
      P.k ≤ ⌊(n : ℝ) * (E + delta)⌋₊ ∧
      ∀ d : ℕ, 0 < d →
        ∀ sigma : Matrix (Fin d × Copies S n) (Fin d × Copies S n) ℂ,
          IsDensityMatrix sigma →
          traceNorm (traceLeft sigma - densityPower rho n) ≤ eta →
          traceNorm (idTensorSuperop P.channel sigma -
            idTensorSuperop (conjSuperop (isometryPower U n)) sigma) ≤ eps

/-- PDF pp.3–4, (12), Fig.4: Alice receives A^n and Bob B^n from the Source.
The q-qubit message is sent forward. Alice keeps l Bell halves; Bob outputs the
other l halves and S^n, the copy of the original source. Local work remains in
the dilations, so no impossible full-input compression isometry is assumed. -/
structure SlepianWolfCode (S A B : Type) [Fintype S] [Fintype A] [Fintype B] (n : ℕ) where
  q : ℕ
  l : ℕ
  aliceWork : ℕ
  bobWork : ℕ
  enc : Matrix ((Qubits l × Qubits q) × Fin (aliceWork + 1)) (Copies A n) ℂ
  dec : Matrix ((Qubits l × Copies S n) × Fin (bobWork + 1))
    (Qubits q × Copies B n) ℂ
  enc_isometry : encᴴ * enc = 1
  dec_isometry : decᴴ * dec = 1

/-- PDF pp.3–4: the normalized source purification in grouped copy registers. -/
def sourceDensity {R X : Type} (v : EuclideanSpace ℂ (R × X)) (n : ℕ) :
    Matrix (Copies R n × Copies X n) (Copies R n × Copies X n) ℂ :=
  pureDensity (WithLp.toLp 2 fun p => ∏ i, v (p.1 i, p.2 i))

/-- PDF p.3: apply only the specified source isometry; Alice and Bob cannot
encode or choose the S input. The reference is untouched. -/
def distributedSource {R S A B : Type} [Fintype S]
    (U : Matrix (A × B) S ℂ) (phi : EuclideanSpace ℂ (R × S)) (n : ℕ) :
    Matrix (Copies R n × (Copies A n × Copies B n))
      (Copies R n × (Copies A n × Copies B n)) ℂ :=
  pureDensity (WithLp.toLp 2 fun p =>
    ∏ i, feedbackOutputVector U phi (p.1 i, (p.2.1 i, p.2.2 i)))

/-- PDF pp.3–4: source/reference and generated Bell pairs are tested jointly.
The ideal S^n output is a copy at Bob, not a free matrix unrelated to phi. -/
def sourceTarget {R S : Type} (phi : EuclideanSpace ℂ (R × S)) (n l : ℕ) :
    Matrix (Copies R n × (Qubits l × (Qubits l × Copies S n)))
      (Copies R n × (Qubits l × (Qubits l × Copies S n))) ℂ :=
  pureDensity (WithLp.toLp 2 fun p =>
    (∏ i, phi (p.1 i, p.2.2.2 i)) * maxEntangled l (p.2.1, p.2.2.1))

/-- PDF pp.3–4: full joint trace-norm error of the physical FQSW code. -/
def SlepianWolfCode.error {R S A B : Type} [Fintype R] [Fintype S] [Fintype A] [Fintype B]
    {n : ℕ} (P : SlepianWolfCode S A B n) (U : Matrix (A × B) S ℂ)
    (phi : EuclideanSpace ℂ (R × S)) : ℝ :=
  traceNorm (idTensorSuperop
    (forwardApply (localOperation P.enc) (localOperation P.dec)) (distributedSource U phi n) -
      sourceTarget phi n P.l)

/-- PDF pp.1,3–4: source transfer with Q forward qubits consumed and E ebits
generated per source copy. No initial Bell pairs or other nonlocal assistance.
The integer output is rounded down and clipped at zero only when E-delta<0. -/
def SourceTransferAchievable {R S A B : Type}
    [Fintype R] [Fintype S] [Fintype A] [Fintype B]
    (U : Matrix (A × B) S ℂ) (phi : EuclideanSpace ℂ (R × S)) (Q E : ℝ) : Prop :=
  ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
    ∀ n ≥ N, 1 ≤ n → ∃ P : SlepianWolfCode S A B n,
      P.q ≤ ⌊(n : ℝ) * (Q + delta)⌋₊ ∧
      ⌊(n : ℝ) * (E - delta)⌋₊ ≤ P.l ∧ P.error U phi ≤ eps

/-- PDF p.3: the identity from the Source to Alice, with Bob's initial share
one-dimensional. S, A-hat, A-prime and B-hat are copies of the same basis type. -/
def sourceToAlice (S : Type) : Matrix (S × Unit) S ℂ :=
  fun p s => if p.1 = s then 1 else 0

/-- PDF pp.2–3: the promise supplied by a relative n-use identity channel.
It guarantees a joint reference test for every input with eta-close marginal. -/
def RelativeIdentityTest {S : Type} [Fintype S] (n : ℕ) (rho : Matrix S S ℂ)
    (F : Superop (Copies S n) (Copies S n)) (eta eps : ℝ) : Prop :=
  ∀ d : ℕ, 0 < d →
    ∀ sigma : Matrix (Fin d × Copies S n) (Fin d × Copies S n) ℂ,
      IsDensityMatrix sigma → traceNorm (traceLeft sigma - densityPower rho n) ≤ eta →
      traceNorm (idTensorSuperop F sigma - sigma) ≤ eps

/-- PDF p.3, Fig.3: local encoding, n uses of the supplied relative identity,
and local decoding. The middle channel is supplied, not chosen by the code. -/
structure RoutingCode (S : Type) [Fintype S] (n : ℕ) where
  enc : Superop (Copies S n) (Copies S n)
  dec : Superop (Copies S n) (Copies S n)
  enc_channel : IsChannel enc
  dec_channel : IsChannel dec

/-- PDF p.3: direct source/identity concatenation, with the actual encoder
marginal tested against the supplied channel's promise. The quantified F represents
any physical implementation satisfying that promise, not unrestricted access.
The code uses exactly n supplied identity calls for n source copies; no nonlocal
ancillas or extra communication are available. delta records allowed rate slack. -/
def SourceRoutingAchievable {R S : Type} [Fintype R] [Fintype S]
    (phi : EuclideanSpace ℂ (R × S)) : Prop :=
  let rho := traceLeft (pureDensity phi)
  ∀ eps > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
    ∀ n ≥ N, 1 ≤ n → ∃ P : RoutingCode S n,
      (n : ℝ) ≤ (n : ℝ) * (1 + delta) ∧
      ∀ eta > (0 : ℝ),
        traceNorm (traceLeft (idTensorSuperop P.enc (sourceDensity phi n)) -
          densityPower rho n) ≤ eta ∧
        ∀ F : Superop (Copies S n) (Copies S n), IsChannel F →
          RelativeIdentityTest n rho F eta eps →
          traceNorm (idTensorSuperop (P.dec ∘ F ∘ P.enc) (sourceDensity phi n) -
            sourceDensity phi n) ≤ eps

/-! Sanity examples for the new definitions; none assumes an entropy identity
or a resource inequality. These are definition checks, not the paper's proofs. -/

example {X : Type} (rho : Matrix X X ℂ) (x y : Copies X 0) :
    densityPower rho 0 x y = 1 := by
  simp [densityPower]

example {R S : Type} (phi : EuclideanSpace ℂ (R × S))
    (r : Copies R 0) (s : Copies S 0) :
    sourceDensity phi 0 (r, s) (r, s) = 1 := by
  simp [sourceDensity, pureDensity]

example {X : Type} [Fintype X] (rho : Matrix X X ℂ) :
    redistributionMarginal rho id (fun _ => ()) = rho := by
  ext a b
  unfold redistributionMarginal
  rw [Finset.sum_eq_single a]
  · simp
  · intro x _ hx
    simp [hx]
  · simp

example {X Y : Type} [Fintype X] (V : Matrix (Y × Unit) X ℂ)
    (rho : Matrix X X ℂ) (y y' : Y) :
    localOperation V rho y y' = conjSuperop V rho (y, ()) (y', ()) := by
  simp [localOperation]

example (s t : Bool) : sourceToAlice Bool (s, ()) t = if s = t then 1 else 0 := by
  by_cases h : s = t <;> simp [sourceToAlice, h]

example {X : Type} (rho : Matrix X X ℂ) (x y : X) :
    appendBell 0 rho ((x, fun _ => false), fun _ => false)
      ((y, fun _ => false), fun _ => false) = rho x y := by
  simp [appendBell, pureDensity, maxEntangled]

example {R S : Type} [Fintype S] [DecidableEq S] (phi : EuclideanSpace ℂ (R × S)) :
    feedbackOutputVector (sourceToAlice S) phi =
      WithLp.toLp 2 (fun p : R × (S × Unit) => phi (p.1, p.2.1)) := by
  ext p
  simp [feedbackOutputVector, sourceToAlice]

example {S : Type} [Fintype S] (n : ℕ) (rho : Matrix S S ℂ) (eta eps : ℝ)
    (heps : 0 ≤ eps) :
    RelativeIdentityTest n rho (fun X => X) eta eps := by
  intro d hd sigma hs hm
  have h : idTensorSuperop (fun X : Matrix (Copies S n) (Copies S n) ℂ => X) sigma = sigma := by
    ext p q
    rfl
  rw [h, sub_self]
  simpa [traceNorm, cfc_zero] using heps

end QAlgorithms.Papers.Independent.Intakethird100.P64
Source
A triangle of dualities: reversibly decomposable quantum channels, source-channel duality, and time reversal, arXiv:quant-ph/0505138v1, https://arxiv.org/abs/quant-ph/0505138v1
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P64.FeedbackSimulationCode.channel

For finite types S,A,BS,A,BS,A,B, a natural number nnn, and a feedback simulation code PPP, write Bj\mathcal B_jBj​ for the length-jjj Boolean strings, V=P.encV=P.\mathrm{enc}V=P.enc, W=P.decW=P.\mathrm{dec}W=P.dec, q=P.qq=P.qq=P.q, k=P.kk=P.kk=P.k, wA=P.aliceWorkw_A=P.\mathrm{aliceWork}wA​=P.aliceWork, and wB=P.bobWorkw_B=P.\mathrm{bobWork}wB​=P.bobWork; these matrices satisfy V†V=IV^\dagger V=IV†V=I and W†W=IW^\dagger W=IW†W=I. This defines a map NP\mathcal N_PNP​ from matrices on SnS^nSn to matrices on An×BnA^n\times B^nAn×Bn whose entries, for every input matrix ρ\rhoρ, are

NP(ρ)((a,b),(a′,b′))=2−k∑s,s′∈Sn∑h,h′∈Bk∑m,m′∈Bq∑u=0wA∑v=0wBρ(s,s′)⋅V(((a,m),u),(s,h)) V(((a′,m′),u),(s′,h′))‾⋅W((b,v),(m,h)) W((b′,v),(m′,h′))‾.\begin{aligned} \mathcal N_P(\rho)((a,b),(a',b')) ={}&2^{-k} \sum_{s,s'\in S^n}\sum_{h,h'\in\mathcal B_k} \sum_{m,m'\in\mathcal B_q} \sum_{u=0}^{w_A}\sum_{v=0}^{w_B} \rho(s,s')\\ &\quad{}\cdot V(((a,m),u),(s,h)) \,\overline{V(((a',m'),u),(s',h'))}\\ &\quad{}\cdot W((b,v),(m,h)) \,\overline{W((b',v),(m',h'))}. \end{aligned}NP​(ρ)((a,b),(a′,b′))=​2−ks,s′∈Sn∑​h,h′∈Bk​∑​m,m′∈Bq​∑​u=0∑wA​​v=0∑wB​​ρ(s,s′)⋅V(((a,m),u),(s,h))V(((a′,m′),u),(s′,h′))​⋅W((b,v),(m,h))W((b′,v),(m′,h′))​.​

Here TnT^nTn denotes the functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to TTT. This formula expands the appended vector 2−k/2∑h∣h,h⟩2^{-k/2}\sum_h|h,h\rangle2−k/2∑h​∣h,h⟩, the two matrix actions, the summation over their work coordinates, and the intermediate message coordinates m,m′m,m'm,m′. It applies to arbitrary ρ\rhoρ, including matrices that are not density matrices, and includes zero values of n,q,k,wA,wBn,q,k,w_A,w_Bn,q,k,wA​,wB​.

QAlgorithms.Papers.Independent.Intakethird100.P64.FeedbackSimulationCode

For finite types S,A,BS,A,BS,A,B and any natural number nnn, this structure consists of natural numbers q,k,wA,wBq,k,w_A,w_Bq,k,wA​,wB​, a complex matrix

V:C Sn×Bk⟶C (An×Bq)×{0,…,wA},V:\mathbb C^{\,S^n\times\mathcal B_k} \longrightarrow \mathbb C^{\,(A^n\times\mathcal B_q)\times\{0,\ldots,w_A\}},V:CSn×Bk​⟶C(An×Bq​)×{0,…,wA​},

and a complex matrix

W:C Bq×Bk⟶C Bn×{0,…,wB},W:\mathbb C^{\,\mathcal B_q\times\mathcal B_k} \longrightarrow \mathbb C^{\,B^n\times\{0,\ldots,w_B\}},W:CBq​×Bk​⟶CBn×{0,…,wB​},

together with the identities V†V=IV^\dagger V=IV†V=I and W†W=IW^\dagger W=IW†W=I on their respective input spaces. Here TnT^nTn means the functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to TTT, and Bj\mathcal B_jBj​ is the set of Boolean strings of length jjj. The fields wA,wBw_A,w_BwA​,wB​ are aliceWork and bobWork, so their associated work spaces have dimensions wA+1,wB+1w_A+1,w_B+1wA​+1,wB​+1. All four natural-number fields may be zero, n=0n=0n=0 is allowed, and no nonemptiness assumption is imposed on S,A,BS,A,BS,A,B.

QAlgorithms.Papers.Independent.Intakethird100.P64.RelativeIdentityTest

For a type SSS equipped with a finite enumeration, n∈Nn\in\mathbb Nn∈N, an arbitrary complex matrix ρ\rhoρ indexed by SSS, an arbitrary function FFF from complex matrices indexed by SnS^nSn to matrices with the same indices, and real numbers η,ε\eta,\varepsilonη,ε, this predicate requires the following: for every positive integer ddd and every positive semidefinite matrix σ\sigmaσ indexed by {0,…,d−1}×Sn\{0,\ldots,d-1\}\times S^n{0,…,d−1}×Sn with tr⁡σ=1\operatorname{tr}\sigma=1trσ=1,

∥tr⁡{0,…,d−1}σ−ρ⊗n∥1≤η⟹∥(id⁡d⊗F)(σ)−σ∥1≤ε.\left\|\operatorname{tr}_{\{0,\ldots,d-1\}}\sigma-\rho^{\otimes n}\right\|_1\le\eta \quad\Longrightarrow\quad \|(\operatorname{id}_d\otimes F)(\sigma)-\sigma\|_1\le\varepsilon.​tr{0,…,d−1}​σ−ρ⊗n​1​≤η⟹∥(idd​⊗F)(σ)−σ∥1​≤ε.

Explicitly, (ρ⊗n)x,y=∏i=0n−1ρxi,yi(\rho^{\otimes n})_{x,y}=\prod_{i=0}^{n-1}\rho_{x_i,y_i}(ρ⊗n)x,y​=∏i=0n−1​ρxi​,yi​​, the partial trace has entry ∑a=0d−1σ(a,x),(a,y)\sum_{a=0}^{d-1}\sigma_{(a,x),(a,y)}∑a=0d−1​σ(a,x),(a,y)​, and (id⁡d⊗F)(σ)(a,x),(b,y)=F(σ(a,⋅),(b,⋅))x,y(\operatorname{id}_d\otimes F)(\sigma)_{(a,x),(b,y)}=F(\sigma_{(a,\cdot),(b,\cdot)})_{x,y}(idd​⊗F)(σ)(a,x),(b,y)​=F(σ(a,⋅),(b,⋅)​)x,y​; the trace norm is ∥M∥1=Re⁡tr⁡M†M\|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M}∥M∥1​=RetrM†M​. There is no density-matrix assumption on ρ\rhoρ, channel assumption on FFF, or sign restriction on η,ε\eta,\varepsilonη,ε. A negative η\etaη makes the implication vacuous. If SSS is empty and n>0n>0n>0, no such trace-one σ\sigmaσ exists, so the predicate is vacuous; n=0n=0n=0 is also included, and its tensor-product matrix is [1][1][1].

QAlgorithms.Papers.Independent.Intakethird100.P64.RelativeSimulationAchievable

For finite types S,A,BS,A,BS,A,B, an arbitrary complex matrix U:(A×B)×S→CU:(A\times B)\times S\to\mathbb CU:(A×B)×S→C, an arbitrary complex matrix ρ\rhoρ on SSS, and real numbers Q,EQ,EQ,E, this proposition says that for every ε>0\varepsilon>0ε>0 and every δ>0\delta>0δ>0, there exist η>0\eta>0η>0 and a natural number NNN such that, for every natural number n≥Nn\ge Nn≥N with n≥1n\ge1n≥1, there is a code PPP with natural numbers q,k,wA,wBq,k,w_A,w_Bq,k,wA​,wB​ and matrices

V:CSn×Bk→C(An×Bq)×{0,…,wA},W:CBq×Bk→CBn×{0,…,wB},V†V=W†W=I,V:\mathbb C^{S^n\times\mathcal B_k}\to \mathbb C^{(A^n\times\mathcal B_q)\times\{0,\ldots,w_A\}}, \qquad W:\mathbb C^{\mathcal B_q\times\mathcal B_k}\to \mathbb C^{B^n\times\{0,\ldots,w_B\}}, \quad V^\dagger V=W^\dagger W=I,V:CSn×Bk​→C(An×Bq​)×{0,…,wA​},W:CBq​×Bk​→CBn×{0,…,wB​},V†V=W†W=I,

for which q≤⌊n(Q+δ)⌋+q\le\lfloor n(Q+\delta)\rfloor_+q≤⌊n(Q+δ)⌋+​ and k≤⌊n(E+δ)⌋+k\le\lfloor n(E+\delta)\rfloor_+k≤⌊n(E+δ)⌋+​, and the following test holds for every positive natural number ddd and every positive semidefinite matrix σ\sigmaσ of trace 111 on {0,…,d−1}×Sn\{0,\ldots,d-1\}\times S^n{0,…,d−1}×Sn:

∥Tr⁡{0,…,d−1}σ−ρ[n]∥1≤η⟹∥(idd⊗NP)(σ)−(idd⊗Un)(σ)∥1≤ε.\left\|\operatorname{Tr}_{\{0,\ldots,d-1\}}\sigma-\rho^{[n]}\right\|_1\le\eta \quad\Longrightarrow\quad \left\|(\mathrm{id}_d\otimes\mathcal N_P)(\sigma) -(\mathrm{id}_d\otimes\mathcal U_n)(\sigma)\right\|_1\le\varepsilon.​Tr{0,…,d−1}​σ−ρ[n]​1​≤η⟹∥(idd​⊗NP​)(σ)−(idd​⊗Un​)(σ)∥1​≤ε.

Here TnT^nTn denotes length-nnn strings over TTT, Bj={false,true}j\mathcal B_j=\{\mathrm{false},\mathrm{true}\}^jBj​={false,true}j, ρ[n](s,s′)=∏iρ(si,si′)\rho^{[n]}(s,s')=\prod_i\rho(s_i,s'_i)ρ[n](s,s′)=∏i​ρ(si​,si′​), U[n]((a,b),s)=∏iU((ai,bi),si)U^{[n]}((a,b),s)=\prod_iU((a_i,b_i),s_i)U[n]((a,b),s)=∏i​U((ai​,bi​),si​), and Un(X)=U[n]X(U[n])†\mathcal U_n(X)=U^{[n]}X(U^{[n]})^\daggerUn​(X)=U[n]X(U[n])†. The code map NP\mathcal N_PNP​ appends the outer product of 2−k/2∑h∣h,h⟩2^{-k/2}\sum_h|h,h\rangle2−k/2∑h​∣h,h⟩, applies VVV to SnS^nSn and the first appended register, applies WWW to the resulting Bq\mathcal B_qBq​ register and the second appended register, and sums over equal work coordinates in the resulting matrix. Explicitly, its matrix entries are the sums of

2−kX(s,s′) V(((a,m),u),(s,h))V(((a′,m′),u),(s′,h′))‾ W((b,v),(m,h))W((b′,v),(m′,h′))‾2^{-k}X(s,s')\, V(((a,m),u),(s,h)) \overline{V(((a',m'),u),(s',h'))}\, W((b,v),(m,h)) \overline{W((b',v),(m',h'))}2−kX(s,s′)V(((a,m),u),(s,h))V(((a′,m′),u),(s′,h′))​W((b,v),(m,h))W((b′,v),(m′,h′))​

over s,s′,h,h′,m,m′,u,vs,s',h,h',m,m',u,vs,s′,h,h′,m,m′,u,v in their respective indexing sets. For any matrix map Φ\PhiΦ, its reference extension acts on each reference block by (idd⊗Φ)(σ)((r,y),(r′,y′))=Φ(σr,r′)(y,y′)(\mathrm{id}_d\otimes\Phi)(\sigma)((r,y),(r',y'))=\Phi(\sigma_{r,r'})(y,y')(idd​⊗Φ)(σ)((r,y),(r′,y′))=Φ(σr,r′​)(y,y′), where σr,r′(s,s′)=σ((r,s),(r′,s′))\sigma_{r,r'}(s,s')=\sigma((r,s),(r',s'))σr,r′​(s,s′)=σ((r,s),(r′,s′)); the partial trace sums the diagonal reference blocks, and ∥M∥1=Re⁡tr⁡M†M\|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M}∥M∥1​=RetrM†M​, with the square root defined by continuous functional calculus. The natural-number floor ⌊x⌋+\lfloor x\rfloor_+⌊x⌋+​ is zero for negative xxx. Neither U†U=IU^\dagger U=IU†U=I, nor positivity or normalization of ρ\rhoρ, nor nonnegativity of Q,EQ,EQ,E is assumed; if no density matrix satisfies the stated proximity condition, that conditional test is vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P64.RoutingCode

For a type SSS equipped with a finite enumeration and n∈Nn\in\mathbb Nn∈N, this structure consists of two functions E,D\mathcal E,\mathcal DE,D from complex matrices indexed by SnS^nSn to matrices with the same indices, together with a finite Kraus representation for each: there exist natural numbers a,ba,ba,b, matrices E0,…,Ea−1E_0,\ldots,E_{a-1}E0​,…,Ea−1​ and D0,…,Db−1D_0,\ldots,D_{b-1}D0​,…,Db−1​, such that

∑j=0a−1Ej†Ej=I,E(X)=∑j=0a−1EjXEj†for every X,\sum_{j=0}^{a-1}E_j^\dagger E_j=I,\qquad \mathcal E(X)=\sum_{j=0}^{a-1}E_jXE_j^\dagger\quad\text{for every }X,j=0∑a−1​Ej†​Ej​=I,E(X)=j=0∑a−1​Ej​XEj†​for every X,

and

∑j=0b−1Dj†Dj=I,D(X)=∑j=0b−1DjXDj†for every X.\sum_{j=0}^{b-1}D_j^\dagger D_j=I,\qquad \mathcal D(X)=\sum_{j=0}^{b-1}D_jXD_j^\dagger\quad\text{for every }X.j=0∑b−1​Dj†​Dj​=I,D(X)=j=0∑b−1​Dj​XDj†​for every X.

No further relationship between the two functions is required. The type SSS may be empty; n=0n=0n=0 is allowed and gives a singleton tuple index set. The Kraus-family sizes may be zero, in which case their sums are zero and the displayed completeness equality still must hold.

QAlgorithms.Papers.Independent.Intakethird100.P64.SlepianWolfCode.error

Let R,S,A,BR,S,A,BR,S,A,B be types equipped with finite enumerations, possibly empty, and let n∈Nn\in\mathbb Nn∈N. A code PPP supplies natural numbers q,l,wA,wBq,l,w_A,w_Bq,l,wA​,wB​ and complex matrices

V:CAn⟶C(Bl×Bq)×{0,…,wA},W:CBq×Bn⟶C(Bl×Sn)×{0,…,wB},V:\mathbb C^{A^n}\longrightarrow \mathbb C^{(\mathcal B_l\times\mathcal B_q)\times\{0,\ldots,w_A\}}, \qquad W:\mathbb C^{\mathcal B_q\times B^n}\longrightarrow \mathbb C^{(\mathcal B_l\times S^n)\times\{0,\ldots,w_B\}},V:CAn⟶C(Bl​×Bq​)×{0,…,wA​},W:CBq​×Bn⟶C(Bl​×Sn)×{0,…,wB​},

satisfying V†V=IV^\dagger V=IV†V=I and W†W=IW^\dagger W=IW†W=I, where Bj={0,1}j\mathcal B_j=\{0,1\}^jBj​={0,1}j. For any matrix U:CS→CA×BU:\mathbb C^S\to\mathbb C^{A\times B}U:CS→CA×B and vector ϕ∈CR×S\phi\in\mathbb C^{R\times S}ϕ∈CR×S, define

vn(r,a,b)=∏i=0n−1∑s∈SU(ai,bi),sϕ(ri,s),tn(r,u,v,s)={2−l/2∏i=0n−1ϕ(ri,si),u=v,0,u≠v.v_n(r,a,b)=\prod_{i=0}^{n-1}\sum_{s\in S}U_{(a_i,b_i),s}\phi(r_i,s), \qquad t_n(r,u,v,s)= \begin{cases} \displaystyle 2^{-l/2}\prod_{i=0}^{n-1}\phi(r_i,s_i),&u=v,\\ 0,&u\ne v. \end{cases}vn​(r,a,b)=i=0∏n−1​s∈S∑​U(ai​,bi​),s​ϕ(ri​,s),tn​(r,u,v,s)=⎩⎨⎧​2−l/2i=0∏n−1​ϕ(ri​,si​),0,​u=v,u=v.​

For each work-index pair α∈{0,…,wA}\alpha\in\{0,\ldots,w_A\}α∈{0,…,wA​}, β∈{0,…,wB}\beta\in\{0,\ldots,w_B\}β∈{0,…,wB​}, define a matrix from An×BnA^n\times B^nAn×Bn to Bl×(Bl×Sn)\mathcal B_l\times(\mathcal B_l\times S^n)Bl​×(Bl​×Sn) by

Kαβ[(u,(v,s)),(a,b)]=∑m∈BqV((u,m),α),a W((v,s),β),(m,b).K_{\alpha\beta}[(u,(v,s)),(a,b)] = \sum_{m\in\mathcal B_q}V_{((u,m),\alpha),a}\, W_{((v,s),\beta),(m,b)}.Kαβ​[(u,(v,s)),(a,b)]=m∈Bq​∑​V((u,m),α),a​W((v,s),β),(m,b)​.

The declared error is

∥∑α,β(IRn⊗Kαβ)vnvn†(IRn⊗Kαβ)†−tntn†∥1,∥M∥1=Re⁡tr⁡M†M,\left\| \sum_{\alpha,\beta}(I_{R^n}\otimes K_{\alpha\beta})v_nv_n^\dagger (I_{R^n}\otimes K_{\alpha\beta})^\dagger -t_nt_n^\dagger \right\|_1, \qquad \|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M},​α,β∑​(IRn​⊗Kαβ​)vn​vn†​(IRn​⊗Kαβ​)†−tn​tn†​​1​,∥M∥1​=RetrM†M​,

with the square root defined by continuous functional calculus. No isometry assumption on UUU or normalization assumption on ϕ\phiϕ is made. All four code parameters may be zero, and n=0n=0n=0 is included, with empty products equal to 111.

QAlgorithms.Papers.Independent.Intakethird100.P64.SlepianWolfCode

For finite types S,A,BS,A,BS,A,B and any natural number nnn, this structure consists of natural numbers q,l,wA,wBq,l,w_A,w_Bq,l,wA​,wB​, a complex matrix

V:CAn⟶C(Bl×Bq)×{0,…,wA},V:\mathbb C^{A^n}\longrightarrow \mathbb C^{(\mathcal B_l\times\mathcal B_q)\times\{0,\ldots,w_A\}},V:CAn⟶C(Bl​×Bq​)×{0,…,wA​},

and a complex matrix

W:CBq×Bn⟶C(Bl×Sn)×{0,…,wB},W:\mathbb C^{\mathcal B_q\times B^n}\longrightarrow \mathbb C^{(\mathcal B_l\times S^n)\times\{0,\ldots,w_B\}},W:CBq​×Bn⟶C(Bl​×Sn)×{0,…,wB​},

together with the identities V†V=IV^\dagger V=IV†V=I and W†W=IW^\dagger W=IW†W=I on their respective input spaces. Here TnT^nTn is the set of functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to TTT, Bj\mathcal B_jBj​ is the set of length-jjj Boolean strings, and wA,wBw_A,w_BwA​,wB​ are the fields aliceWork and bobWork, whose work spaces have dimensions wA+1,wB+1w_A+1,w_B+1wA​+1,wB​+1. Zero values of all natural-number fields and n=0n=0n=0 are included, and the finite indexing types are not assumed nonempty.

QAlgorithms.Papers.Independent.Intakethird100.P64.SourceRoutingAchievable

For types R,SR,SR,S equipped with finite enumerations, possibly empty, and an arbitrary vector ϕ∈CR×S\phi\in\mathbb C^{R\times S}ϕ∈CR×S, put

ρs,t=∑r∈Rϕ(r,s)ϕ(r,t)‾,Dn[(r,s),(r′,t)]=∏i=0n−1ϕ(ri,si) ∏i=0n−1ϕ(ri′,ti)‾,(ρ⊗n)s,t=∏i=0n−1ρsi,ti.\rho_{s,t}=\sum_{r\in R}\phi(r,s)\overline{\phi(r,t)},\qquad D_n[(r,s),(r',t)] =\prod_{i=0}^{n-1}\phi(r_i,s_i)\, \overline{\prod_{i=0}^{n-1}\phi(r'_i,t_i)},\qquad (\rho^{\otimes n})_{s,t}=\prod_{i=0}^{n-1}\rho_{s_i,t_i}.ρs,t​=r∈R∑​ϕ(r,s)ϕ(r,t)​,Dn​[(r,s),(r′,t)]=i=0∏n−1​ϕ(ri​,si​)i=0∏n−1​ϕ(ri′​,ti​)​,(ρ⊗n)s,t​=i=0∏n−1​ρsi​,ti​​.

This predicate means that for every ε>0,δ>0\varepsilon>0,\delta>0ε>0,δ>0, there is N∈NN\in\mathbb NN∈N such that for every n≥Nn\ge Nn≥N with n≥1n\ge1n≥1, there are channels E,D\mathcal E,\mathcal DE,D on matrices indexed by SnS^nSn satisfying n≤n(1+δ)n\le n(1+\delta)n≤n(1+δ) and, for every η>0\eta>0η>0,

∥tr⁡Rn(id⁡Rn⊗E)(Dn)−ρ⊗n∥1≤η,\left\|\operatorname{tr}_{R^n}(\operatorname{id}_{R^n}\otimes\mathcal E)(D_n) -\rho^{\otimes n}\right\|_1\le\eta,​trRn​(idRn​⊗E)(Dn​)−ρ⊗n​1​≤η,

and the following implication for every channel FFF on that same matrix space: if for every positive integer ddd and every positive semidefinite trace-one matrix σ\sigmaσ on {0,…,d−1}×Sn\{0,\ldots,d-1\}\times S^n{0,…,d−1}×Sn,

∥tr⁡{0,…,d−1}σ−ρ⊗n∥1≤η ⟹ ∥(id⁡d⊗F)(σ)−σ∥1≤ε,\left\|\operatorname{tr}_{\{0,\ldots,d-1\}}\sigma-\rho^{\otimes n}\right\|_1\le\eta \ \Longrightarrow\ \|(\operatorname{id}_d\otimes F)(\sigma)-\sigma\|_1\le\varepsilon,​tr{0,…,d−1}​σ−ρ⊗n​1​≤η ⟹ ∥(idd​⊗F)(σ)−σ∥1​≤ε,

then

∥(id⁡Rn⊗(D∘F∘E))(Dn)−Dn∥1≤ε.\|(\operatorname{id}_{R^n}\otimes(\mathcal D\circ F\circ\mathcal E))(D_n)-D_n\|_1 \le\varepsilon.∥(idRn​⊗(D∘F∘E))(Dn​)−Dn​∥1​≤ε.

Here a channel means a function T(X)=∑j=0a−1KjXKj†T(X)=\sum_{j=0}^{a-1}K_jXK_j^\daggerT(X)=∑j=0a−1​Kj​XKj†​ for some natural number aaa and matrices satisfying ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I; partial traces sum diagonal blocks in the traced index, and (id⁡Z⊗T)(X)(z,s),(z′,t)=T(X(z,⋅),(z′,⋅))s,t(\operatorname{id}_Z\otimes T)(X)_{(z,s),(z',t)}=T(X_{(z,\cdot),(z',\cdot)})_{s,t}(idZ​⊗T)(X)(z,s),(z′,t)​=T(X(z,⋅),(z′,⋅)​)s,t​. The norm is ∥M∥1=Re⁡tr⁡M†M\|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M}∥M∥1​=RetrM†M​. No normalization of ϕ\phiϕ is assumed. The chosen channels precede the quantifier over η\etaη, so the first norm bound must hold for every positive η\etaη, forcing that difference to be zero; the numerical inequality involving δ\deltaδ is also explicitly required.

QAlgorithms.Papers.Independent.Intakethird100.P64.SourceTransferAchievable

For types R,S,A,BR,S,A,BR,S,A,B equipped with finite enumerations, possibly empty, any complex matrix U:CS→CA×BU:\mathbb C^S\to\mathbb C^{A\times B}U:CS→CA×B, vector ϕ∈CR×S\phi\in\mathbb C^{R\times S}ϕ∈CR×S, and real numbers Q,EQ,EQ,E, this predicate means that for every ε>0\varepsilon>0ε>0 and δ>0\delta>0δ>0, there is N∈NN\in\mathbb NN∈N such that for every n≥Nn\ge Nn≥N with n≥1n\ge1n≥1, there exist natural numbers q,l,wA,wBq,l,w_A,w_Bq,l,wA​,wB​ and matrices

V:CAn→C(Bl×Bq)×{0,…,wA},W:CBq×Bn→C(Bl×Sn)×{0,…,wB},V:\mathbb C^{A^n}\to \mathbb C^{(\mathcal B_l\times\mathcal B_q)\times\{0,\ldots,w_A\}}, \qquad W:\mathbb C^{\mathcal B_q\times B^n}\to \mathbb C^{(\mathcal B_l\times S^n)\times\{0,\ldots,w_B\}},V:CAn→C(Bl​×Bq​)×{0,…,wA​},W:CBq​×Bn→C(Bl​×Sn)×{0,…,wB​},

where Bj={0,1}j\mathcal B_j=\{0,1\}^jBj​={0,1}j, such that V†V=IV^\dagger V=IV†V=I, W†W=IW^\dagger W=IW†W=I, and

q≤⌊n(Q+δ)⌋+,⌊n(E−δ)⌋+≤l,∥∑α,β(I⊗Kαβ)vnvn†(I⊗Kαβ)†−tntn†∥1≤ε.q\le \lfloor n(Q+\delta)\rfloor_+,\qquad \lfloor n(E-\delta)\rfloor_+\le l,\qquad \left\|\sum_{\alpha,\beta}(I\otimes K_{\alpha\beta})v_nv_n^\dagger (I\otimes K_{\alpha\beta})^\dagger-t_nt_n^\dagger\right\|_1\le\varepsilon.q≤⌊n(Q+δ)⌋+​,⌊n(E−δ)⌋+​≤l,​α,β∑​(I⊗Kαβ​)vn​vn†​(I⊗Kαβ​)†−tn​tn†​​1​≤ε.

Here ⌊x⌋+=max⁡(0,⌊x⌋)\lfloor x\rfloor_+=\max(0,\lfloor x\rfloor)⌊x⌋+​=max(0,⌊x⌋), ∥M∥1=Re⁡tr⁡M†M\|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M}∥M∥1​=RetrM†M​, the work indices run through α=0,…,wA\alpha=0,\ldots,w_Aα=0,…,wA​, β=0,…,wB\beta=0,\ldots,w_Bβ=0,…,wB​, and

vn(r,a,b)=∏i=0n−1∑s∈SU(ai,bi),sϕ(ri,s),tn(r,u,v,s)=2−l/21{u=v}∏i=0n−1ϕ(ri,si),Kαβ[(u,(v,s)),(a,b)]=∑m∈BqV((u,m),α),aW((v,s),β),(m,b).\begin{aligned} v_n(r,a,b)&=\prod_{i=0}^{n-1}\sum_{s\in S}U_{(a_i,b_i),s}\phi(r_i,s),\\ t_n(r,u,v,s)&=2^{-l/2}\mathbf1_{\{u=v\}}\prod_{i=0}^{n-1}\phi(r_i,s_i),\\ K_{\alpha\beta}[(u,(v,s)),(a,b)] &=\sum_{m\in\mathcal B_q}V_{((u,m),\alpha),a} W_{((v,s),\beta),(m,b)}. \end{aligned}vn​(r,a,b)tn​(r,u,v,s)Kαβ​[(u,(v,s)),(a,b)]​=i=0∏n−1​s∈S∑​U(ai​,bi​),s​ϕ(ri​,s),=2−l/21{u=v}​i=0∏n−1​ϕ(ri​,si​),=m∈Bq​∑​V((u,m),α),a​W((v,s),β),(m,b)​.​

No normalization of ϕ\phiϕ, isometry of UUU, or nonnegativity of Q,EQ,EQ,E is required. In particular, a nonpositive argument of either natural-number floor produces zero; the entanglement condition is only the displayed lower bound on lll. The code parameters may be zero, while the quantified block lengths are positive.

QAlgorithms.Papers.Independent.Intakethird100.P64.appendBell

For any type XXX, any natural number kkk, and any complex matrix ρ\rhoρ on XXX, this defines a matrix on (X×Bk)×Bk(X\times\mathcal B_k)\times\mathcal B_k(X×Bk​)×Bk​, where Bk\mathcal B_kBk​ is the set of length-kkk Boolean strings, by

Ωρ((x,u,v),(x′,u′,v′))={2−kρ(x,x′),u=v and u′=v′,0,otherwise.\Omega_\rho((x,u,v),(x',u',v')) = \begin{cases} 2^{-k}\rho(x,x'),&u=v\ \text{and}\ u'=v',\\ 0,&\text{otherwise}. \end{cases}Ωρ​((x,u,v),(x′,u′,v′))={2−kρ(x,x′),0,​u=v and u′=v′,otherwise.​

Equivalently, it multiplies each entry of ρ\rhoρ by the corresponding entry of βkβk†\beta_k\beta_k^\daggerβk​βk†​, where βk(u,v)=2−k/2\beta_k(u,v)=2^{-k/2}βk​(u,v)=2−k/2 when u=vu=vu=v, and zero otherwise. No positivity or trace condition on ρ\rhoρ is required. For k=0k=0k=0, each Boolean-string register has one value and the added factor is 111.

QAlgorithms.Papers.Independent.Intakethird100.P64.distributedSource

For arbitrary types R,S,A,BR,S,A,BR,S,A,B, with SSS equipped with a finite enumeration, any complex matrix UUU with rows indexed by A×BA\times BA×B and columns by SSS, any vector ϕ\phiϕ in the complex Euclidean space indexed by R×SR\times SR×S, and n∈Nn\in\mathbb Nn∈N, define

w(r,a,b)=∑s∈SU(a,b),sϕ(r,s),vn(r,a,b)=∏i=0n−1w(ri,ai,bi).w(r,a,b)=\sum_{s\in S}U_{(a,b),s}\phi(r,s), \qquad v_n(r,a,b)=\prod_{i=0}^{n-1}w(r_i,a_i,b_i).w(r,a,b)=s∈S∑​U(a,b),s​ϕ(r,s),vn​(r,a,b)=i=0∏n−1​w(ri​,ai​,bi​).

The declaration defines the matrix indexed by Rn×(An×Bn)R^n\times(A^n\times B^n)Rn×(An×Bn) whose (p,p′)(p,p')(p,p′)-entry is vn(p)vn(p′)‾v_n(p)\overline{v_n(p')}vn​(p)vn​(p′)​. Neither U†U=IU^\dagger U=IU†U=I nor normalization of ϕ\phiϕ is assumed. Empty types are allowed; in particular, if SSS is empty, each sum defining www is zero. When n=0n=0n=0, the product is 111 and the matrix is [1][1][1].

QAlgorithms.Papers.Independent.Intakethird100.P64.feedbackEntropyA

For finite types R,A,BR,A,BR,A,B and any complex vector ψ\psiψ indexed by R×(A×B)R\times(A\times B)R×(A×B), this defines the real number H(ρA)H(\rho_A)H(ρA​), where

ρA(a,a′)=∑r∈R∑b∈Bψ(r,(a,b)) ψ(r,(a′,b))‾,H(M)=−Re⁡tr⁡ ⁣[f(M)],f(x)=xlog⁡2x.\rho_A(a,a')=\sum_{r\in R}\sum_{b\in B} \psi(r,(a,b))\,\overline{\psi(r,(a',b))}, \qquad H(M)=-\operatorname{Re}\operatorname{tr}\!\left[f(M)\right], \quad f(x)=x\log_2 x.ρA​(a,a′)=r∈R∑​b∈B∑​ψ(r,(a,b))ψ(r,(a′,b))​,H(M)=−Retr[f(M)],f(x)=xlog2​x.

Here f(M)f(M)f(M) is evaluated using continuous functional calculus. Thus the matrix whose entropy is taken is obtained from the outer product ψψ†\psi\psi^\daggerψψ† by summing over equal RRR- and BBB-coordinates. The vector need not have norm one, and the finite types may be empty; the definition also applies to the zero vector.

QAlgorithms.Papers.Independent.Intakethird100.P64.feedbackOutputVector

For arbitrary types R,S,A,BR,S,A,BR,S,A,B, with SSS finite, a complex matrix UUU with rows indexed by A×BA\times BA×B and columns indexed by SSS, and a complex coordinate vector ϕ\phiϕ indexed by R×SR\times SR×S, this defines the vector ψ\psiψ indexed by R×(A×B)R\times(A\times B)R×(A×B) through

ψ(r,(a,b))=∑s∈SU((a,b),s) ϕ(r,s).\psi(r,(a,b))=\sum_{s\in S}U((a,b),s)\,\phi(r,s).ψ(r,(a,b))=s∈S∑​U((a,b),s)ϕ(r,s).

There is no normalization assumption on ϕ\phiϕ and no isometry assumption on UUU. If SSS is empty, the sum is empty and the output is zero.

QAlgorithms.Papers.Independent.Intakethird100.P64.isometryPower

For arbitrary types S,A,BS,A,BS,A,B, any complex matrix UUU with rows indexed by A×BA\times BA×B and columns by SSS, and any natural number nnn, this defines the matrix U[n]U^{[n]}U[n] with rows indexed by An×BnA^n\times B^nAn×Bn and columns by SnS^nSn through

U[n]((a,b),s)=∏i=0n−1U((ai,bi),si),U^{[n]}((a,b),s)=\prod_{i=0}^{n-1}U((a_i,b_i),s_i),U[n]((a,b),s)=i=0∏n−1​U((ai​,bi​),si​),

where TnT^nTn means the functions from {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} to TTT. There is no finiteness or isometry hypothesis on UUU or its indexing types. When n=0n=0n=0, the product is 111, and each indexing function space has one element, so the result is the one-by-one matrix with entry 111.

QAlgorithms.Papers.Independent.Intakethird100.P64.localOperation

For arbitrary types X,Y,WX,Y,WX,Y,W, with XXX and WWW finite, and any complex matrix VVV whose rows are indexed by Y×WY\times WY×W and columns by XXX, this defines a map LV\mathcal L_VLV​ from complex matrices on XXX to complex matrices on YYY by

LV(ρ)(y,y′)=∑w∈W(VρV†)((y,w),(y′,w))=∑w∈W∑x,x′∈XV((y,w),x)ρ(x,x′)V((y′,w),x′)‾.\mathcal L_V(\rho)(y,y') =\sum_{w\in W}(V\rho V^\dagger)((y,w),(y',w)) =\sum_{w\in W}\sum_{x,x'\in X} V((y,w),x)\rho(x,x')\overline{V((y',w),x')}.LV​(ρ)(y,y′)=w∈W∑​(VρV†)((y,w),(y′,w))=w∈W∑​x,x′∈X∑​V((y,w),x)ρ(x,x′)V((y′,w),x′)​.

The input ρ\rhoρ is arbitrary, and VVV is not required to satisfy an isometry identity. If XXX or WWW is empty, the resulting map is zero.

QAlgorithms.Papers.Independent.Intakethird100.P64.mutualAB

For finite types R,A,BR,A,BR,A,B and any complex vector ψ\psiψ indexed by R×(A×B)R\times(A\times B)R×(A×B), this defines

H(ρA)+H(ρB)−H(ρAB),H(\rho_A)+H(\rho_B)-H(\rho_{AB}),H(ρA​)+H(ρB​)−H(ρAB​),

where

ρA(a,a′)=∑r,bψ(r,(a,b)) ψ(r,(a′,b))‾,ρB(b,b′)=∑r,aψ(r,(a,b)) ψ(r,(a,b′))‾,ρAB((a,b),(a′,b′))=∑rψ(r,(a,b)) ψ(r,(a′,b′))‾.\begin{aligned} \rho_A(a,a')&=\sum_{r,b}\psi(r,(a,b))\,\overline{\psi(r,(a',b))},\\ \rho_B(b,b')&=\sum_{r,a}\psi(r,(a,b))\,\overline{\psi(r,(a,b'))},\\ \rho_{AB}((a,b),(a',b'))&=\sum_r\psi(r,(a,b))\,\overline{\psi(r,(a',b'))}. \end{aligned}ρA​(a,a′)ρB​(b,b′)ρAB​((a,b),(a′,b′))​=r,b∑​ψ(r,(a,b))ψ(r,(a′,b))​,=r,a∑​ψ(r,(a,b))ψ(r,(a,b′))​,=r∑​ψ(r,(a,b))ψ(r,(a′,b′))​.​

The real-valued function used here is H(M)=−Re⁡tr⁡[f(M)]H(M)=-\operatorname{Re}\operatorname{tr}[f(M)]H(M)=−Retr[f(M)], where f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x is applied by continuous functional calculus. The definition includes zero and unnormalized vectors and permits empty finite indexing types.

QAlgorithms.Papers.Independent.Intakethird100.P64.mutualRA

For finite types R,A,BR,A,BR,A,B and any complex vector ψ\psiψ indexed by R×(A×B)R\times(A\times B)R×(A×B), this defines

H(ρR)+H(ρA)−H(ρRA),H(\rho_R)+H(\rho_A)-H(\rho_{RA}),H(ρR​)+H(ρA​)−H(ρRA​),

where

ρR(r,r′)=∑a,bψ(r,(a,b)) ψ(r′,(a,b))‾,ρA(a,a′)=∑r,bψ(r,(a,b)) ψ(r,(a′,b))‾,ρRA((r,a),(r′,a′))=∑bψ(r,(a,b)) ψ(r′,(a′,b))‾.\begin{aligned} \rho_R(r,r')&=\sum_{a,b}\psi(r,(a,b))\,\overline{\psi(r',(a,b))},\\ \rho_A(a,a')&=\sum_{r,b}\psi(r,(a,b))\,\overline{\psi(r,(a',b))},\\ \rho_{RA}((r,a),(r',a'))&=\sum_b\psi(r,(a,b))\,\overline{\psi(r',(a',b))}. \end{aligned}ρR​(r,r′)ρA​(a,a′)ρRA​((r,a),(r′,a′))​=a,b∑​ψ(r,(a,b))ψ(r′,(a,b))​,=r,b∑​ψ(r,(a,b))ψ(r,(a′,b))​,=b∑​ψ(r,(a,b))ψ(r′,(a′,b))​.​

For each matrix MMM, H(M)=−Re⁡tr⁡[f(M)]H(M)=-\operatorname{Re}\operatorname{tr}[f(M)]H(M)=−Retr[f(M)], with f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x evaluated by continuous functional calculus. No normalization of ψ\psiψ or nonemptiness of the finite types is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P64.mutualRB

For finite types R,A,BR,A,BR,A,B and any complex vector ψ\psiψ indexed by R×(A×B)R\times(A\times B)R×(A×B), this defines

H(ρR)+H(ρB)−H(ρRB),H(\rho_R)+H(\rho_B)-H(\rho_{RB}),H(ρR​)+H(ρB​)−H(ρRB​),

where

ρR(r,r′)=∑a,bψ(r,(a,b)) ψ(r′,(a,b))‾,ρB(b,b′)=∑r,aψ(r,(a,b)) ψ(r,(a,b′))‾,ρRB((r,b),(r′,b′))=∑aψ(r,(a,b)) ψ(r′,(a,b′))‾.\begin{aligned} \rho_R(r,r')&=\sum_{a,b}\psi(r,(a,b))\,\overline{\psi(r',(a,b))},\\ \rho_B(b,b')&=\sum_{r,a}\psi(r,(a,b))\,\overline{\psi(r,(a,b'))},\\ \rho_{RB}((r,b),(r',b'))&=\sum_a\psi(r,(a,b))\,\overline{\psi(r',(a,b'))}. \end{aligned}ρR​(r,r′)ρB​(b,b′)ρRB​((r,b),(r′,b′))​=a,b∑​ψ(r,(a,b))ψ(r′,(a,b))​,=r,a∑​ψ(r,(a,b))ψ(r,(a,b′))​,=a∑​ψ(r,(a,b))ψ(r′,(a,b′))​.​

Here H(M)=−Re⁡tr⁡[f(M)]H(M)=-\operatorname{Re}\operatorname{tr}[f(M)]H(M)=−Retr[f(M)], with f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x interpreted through continuous functional calculus. There is no unit-norm hypothesis on ψ\psiψ, and no nonemptiness hypothesis on the finite types.

QAlgorithms.Papers.Independent.Intakethird100.P64.sourceDensity

For arbitrary types R,XR,XR,X, a vector v∈CR×Xv\in\mathbb C^{R\times X}v∈CR×X in the specified Euclidean space, and n∈Nn\in\mathbb Nn∈N, write RnR^nRn and XnX^nXn for tuples indexed by {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}. This defines the matrix with entries

Dn(v)(r,x),(r′,x′)=(∏i=0n−1v(ri,xi))(∏i=0n−1v(ri′,xi′))‾.D_n(v)_{(r,x),(r',x')} = \left(\prod_{i=0}^{n-1}v(r_i,x_i)\right) \overline{\left(\prod_{i=0}^{n-1}v(r'_i,x'_i)\right)}.Dn​(v)(r,x),(r′,x′)​=(i=0∏n−1​v(ri​,xi​))(i=0∏n−1​v(ri′​,xi′​))​.

There is no normalization assumption on vvv and no finiteness or nonemptiness assumption on R,XR,XR,X. For n=0n=0n=0, both tuple sets are singletons, the products are 111, and the resulting matrix is the one-by-one matrix [1][1][1].

QAlgorithms.Papers.Independent.Intakethird100.P64.sourceTarget

For arbitrary types R,SR,SR,S, any vector ϕ\phiϕ in the complex Euclidean space indexed by R×SR\times SR×S, and natural numbers n,ln,ln,l, let Bl={0,1}l\mathcal B_l=\{0,1\}^lBl​={0,1}l. This defines the outer-product matrix tt†tt^\daggertt†, indexed by Rn×(Bl×(Bl×Sn))R^n\times(\mathcal B_l\times(\mathcal B_l\times S^n))Rn×(Bl​×(Bl​×Sn)), where

t(r,u,v,s)={12l∏i=0n−1ϕ(ri,si),u=v,0,u≠v.t(r,u,v,s)= \begin{cases} \displaystyle \frac{1}{\sqrt{2^l}}\prod_{i=0}^{n-1}\phi(r_i,s_i),&u=v,\\ 0,&u\ne v. \end{cases}t(r,u,v,s)=⎩⎨⎧​2l​1​i=0∏n−1​ϕ(ri​,si​),0,​u=v,u=v.​

Thus its entry at p,p′p,p'p,p′ is t(p)t(p′)‾t(p)\overline{t(p')}t(p)t(p′)​. No normalization, finiteness, or nonemptiness assumptions are imposed on ϕ,R,S\phi,R,Sϕ,R,S. For n=0n=0n=0, the product is 111; for l=0l=0l=0, both bit-string registers are singletons and their equality condition always holds.

QAlgorithms.Papers.Independent.Intakethird100.P64.sourceToAlice

For every type SSS, including an empty or infinite type, this defines the complex matrix with rows indexed by S×{∗}S\times\{\ast\}S×{∗} and columns indexed by SSS, whose entries are

J(a,∗),s={1,a=s,0,a≠s.J_{(a,\ast),s}= \begin{cases} 1,&a=s,\\ 0,&a\ne s. \end{cases}J(a,∗),s​={1,0,​a=s,a=s.​

No finiteness assumption on SSS is imposed.

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