Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for The resource theory of steering

Definition
QAlgorithms_Papers_Independent_Intakethird100_P36

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

prove2quantum

The objects the formalized statements of The resource theory of steering (arXiv:1409.5804) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • AssemblageBlocks: PDF p.2, (2): subnormalized density blocks indexed by output and input.
  • bobTotal: PDF p.2, (4): the reduced matrix. Averaging over inputs avoids an arbitrary chosen input. On a positive-size no-signaling assemblage this is each input's sum.
  • assemblageMass: PDF p.2, (5): common trace of an assemblage.
  • IsAssemblage: PDF p.2, (1)–(5): positive finite dimensions, positive blocks, no signaling, and a common trace in [0,1]. Empty physical systems and alphabets are excluded.
  • IsNormalized: PDF p.2, (5): physical input assemblages have trace one.
  • HasFiniteLHSModel: PDF p.2, (6): a normalized LHS model multiplied by the assemblage's assemblageMass. Finite hidden alphabets suffice: responses can be expanded into deterministic strings of outputs and hidden positive matrices grouped by those strings. In particular this definition includes the zero assemblage at positive d,r,s.
  • SteeringProtocol: PDF pp.3–4, (7),(10), Appendix A: Bob's subinstrument followed by his classical message and Alice's two normalized conditional wirings.
  • SteeringProtocol.Complete: PDF p.4, (12): completeness means Bob never aborts.
  • bobStep: PDF p.9, Step 1: Bob acts before Alice chooses her input.
  • protocolBranch: PDF p.9, Steps 2–3: sample x using (w,xf), query Alice's device, then sample af using (w,a,x,xf), retaining Bob's branch w.
  • protocolAction: PDF p.3, (9): the operational map with the outcome forgotten.
  • oneWayFormula: PDF p.10, (A6), equivalently (9)–(10), written as a closed block formula.
  • steeringBranchProbability: PDF p.3, (8): Bob's branch success probability, independent of Alice's input.
  • normalizeAssemblage: PDF p.4: conditional normalization. All physical uses require positive assemblageMass; at zero this raw array is zero and is never called a normalized assemblage.
  • SteeringConverts: PDF p.7 and Appendices E–F: exact single-copy conversion with strictly positive input-independent success. All allowed wirings are quantified.
  • assemblageMixture: PDF p.5, (14): a binary affine combination, also used with nu>1 in (21).
  • assemblageFamilyMixture: PDF p.11, (C2): a finite weighted assemblageMixture with the same index set throughout.
  • SteeringMeasure: PDF p.5, Definition 2: a measure for every finite dimension and alphabet.
  • measureBranch: PDF p.5, (13): zero-probability outcomes contribute zero, without assigning a physical conditional assemblage to an impossible event.
  • StrongMonotonicity: PDF p.5, Definition 2(ii): average monotonicity for complete instruments, with the outcome available to Alice before the average is taken.
  • MeasureConvexity: PDF p.5, Definition 2(iii): convexity, including both endpoints.
  • IsConvexSteeringMonotone: PDF p.5, Definition 2: nonnegative, vanishes on LHS, average monotone, and convex. There is no converse assertion in the vanishing clause.
  • assemblageCQMatrix: PDF p.2, (3): Alice's classical outcome flag, with Bob's actual block.
  • densityD: PDF p.5, (15): the pinned Umegaki divergence, with natural logarithms and +infinity on support failure. The non-density branch is a totalization outside the physical domain; it is never a finite substitute for divergence.
  • steeringClassicalDivergence: PDF p.5, (16): KL divergence of two normalized finite laws, computed as the Umegaki divergence of their diagonal density matrices. This retains negative individual log-ratio contributions before forming the nonnegative total.
  • EntropyStrategy: PDF p.6, (18)–(19): arbitrary finite-outcome complete Bob instrument and an input law chosen after its outcome. The number of outcomes is unrestricted.
  • strategyProbability: PDF p.6, (19): the two input assemblages give separately computed laws.
  • filteredCQ: PDF p.6, numerator in (18): Bob's filtered flagged state at fixed input.
  • conditionalD: PDF p.6, (18): p D(A/p || B/q), zero when p=0 and infinite when p>0=q. This homogeneous convention avoids a fictitious zero-over-zero density matrix.
  • entropyScore: PDF p.6, (18): both the classical outcome KL term and the conditional quantum term weighted by P(x|g) P_g, with separate denominators P_g and P'_g.
  • assemblageEntropy: PDF p.6, Definition 3: complete-order optimization, allowing +infinity. The source's maximum attainment is asserted in relativeEntropy_convexMonotone; it is not assumed or hidden by replacing the range with a bounded real set.
  • EntropyMaximum: PDF p.6, (18): the explicit maximum-attainment obligation.
  • steeringEntropy: PDF p.6, (20): optimize over normalized LHS assemblages on the same space and alphabets. Minimum attainment and finiteness are conclusions below.
  • relativeEntropyMeasure: PDF p.6, Theorem 2: real-valued relative entropy of steering. Its finiteness on every normalized finite assemblage is an explicit theorem obligation.
  • entropyBranch: PDF p.11, (C1): first-input probability weights, with each output separately normalized. If the second branch is impossible and the first is not, the term is infinite; no support-inclusion hypothesis has been added.
  • WeightCandidate: PDF p.6, Definition 4: all raw nu>=0 candidates are retained, including nu>1. Feasible noise and free assemblages are normalized and no-signaling.
  • steerableWeight: PDF p.6, Definition 4: the real infimum; the theorem explicitly asserts that it is an attained minimum, which rules out junk empty/unbounded values.
  • RobustnessCandidate: PDF p.6, Definition 5: arbitrary assemblage noise, not necessarily LHS.
  • steeringRobustness: PDF p.6, Definition 5: the attained minimum is asserted in Theorem 3.
  • conditionalKetProjector: PDF p.6, (23): pure states with no chosen global phase.
  • IsPureAssemblage: PDF p.6, (23): every conditional state has a unit-vector representative, including zero-probability events; purity does not require a pure Bob marginal.
  • IsPureOrthogonal: PDF p.6, following (23): an orthonormal family at each input.
  • angleKet: PDF p.15, (F5): the two qubit bases used for the infinite-class witnesses.
  • angleAssemblage: PDF p.15, (F5): each rank-one block has probability one half.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_ExtendedQuantumEntropy
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P56

/-!
# The resource theory of steering (arXiv:1409.5804v3)

All locators are PDF pages. Steering-specific definitions below are local to this
paper; IsLaw is imported from the pinned accepted Information100.P56 module. The
finite positive dimensions are unrestricted. Raw block arrays let the same
operation act on normalized and subnormalized assemblages, including zero.
Theorem 4 and the qualitative security motivation are deferred in HARD.md.
Compilation is pending the separate mechanical handoff.
-/

noncomputable section

open scoped BigOperators ENNReal ComplexConjugate ComplexOrder

namespace QAlgorithms.Papers.Independent.Intakethird100.P36

open QAlgorithms.Papers.Independent.Information100.P56 (IsLaw)

attribute [local instance] Classical.propDecidable

/-- PDF p.2, (2): subnormalized density blocks indexed by output and input. -/
abbrev AssemblageBlocks (d r s : ℕ) := Fin r → Fin s → Matrix (Fin d) (Fin d) ℂ

variable {d r s df rf sf m : ℕ}

/-- PDF p.2, (4): the reduced matrix. Averaging over inputs avoids an arbitrary
chosen input. On a positive-size no-signaling assemblage this is each input's sum. -/
def bobTotal (rho : AssemblageBlocks d r s) : Matrix (Fin d) (Fin d) ℂ :=
  ((s : ℂ)⁻¹) • ∑ x, ∑ a, rho a x

/-- PDF p.2, (5): common trace of an assemblage. -/
def assemblageMass (rho : AssemblageBlocks d r s) : ℝ := (bobTotal rho).trace.re

/-- PDF p.2, (1)–(5): positive finite dimensions, positive blocks, no signaling,
and a common trace in [0,1]. Empty physical systems and alphabets are excluded. -/
def IsAssemblage (rho : AssemblageBlocks d r s) : Prop :=
  0 < d ∧ 0 < r ∧ 0 < s ∧
  (∀ a x, (rho a x).PosSemidef) ∧
  (∀ x y, ∑ a, rho a x = ∑ a, rho a y) ∧ 0 ≤ assemblageMass rho ∧ assemblageMass rho ≤ 1

/-- PDF p.2, (5): physical input assemblages have trace one. -/
def IsNormalized (rho : AssemblageBlocks d r s) : Prop := IsAssemblage rho ∧ assemblageMass rho = 1

/-- PDF p.2, (6): a normalized LHS model multiplied by the assemblage's assemblageMass.
Finite hidden alphabets suffice: responses can be expanded into deterministic
strings of outputs and hidden positive matrices grouped by those strings.
In particular this definition includes the zero assemblage at positive d,r,s. -/
def HasFiniteLHSModel (rho : AssemblageBlocks d r s) : Prop :=
  IsAssemblage rho ∧ ∃ (n : ℕ) (p : Fin n → ℝ)
    (xi : Fin n → Matrix (Fin d) (Fin d) ℂ)
    (response : Fin n → Fin s → Fin r → ℝ),
    IsLaw p ∧ (∀ l, QAlgorithms.IsDensityMatrix (xi l)) ∧
    (∀ l x, IsLaw (response l x)) ∧
    ∀ a x, rho a x = ∑ l, ((assemblageMass rho * p l * response l x a : ℝ) : ℂ) • xi l

/-- PDF pp.3–4, (7),(10), Appendix A: Bob's subinstrument followed by
his classical message and Alice's two normalized conditional wirings. -/
structure SteeringProtocol (d r s df rf sf m : ℕ) where
  dimensions : 0 < d ∧ 0 < r ∧ 0 < s ∧ 0 < df ∧ 0 < rf ∧ 0 < sf
  kraus : Fin m → Matrix (Fin df) (Fin d) ℂ
  trace_decreasing : (1 - ∑ w, (kraus w).conjTranspose * kraus w).PosSemidef
  input : Fin m → Fin sf → Fin s → ℝ
  output : Fin m → Fin rf → Fin r → Fin s → Fin sf → ℝ
  input_law : ∀ w xf, IsLaw (input w xf)
  output_law : ∀ w a x xf, IsLaw (fun af => output w af a x xf)

/-- PDF p.4, (12): completeness means Bob never aborts. -/
def SteeringProtocol.Complete (P : SteeringProtocol d r s df rf sf m) : Prop :=
  ∑ w, (P.kraus w).conjTranspose * P.kraus w = 1

/-- PDF p.9, Step 1: Bob acts before Alice chooses her input. -/
def bobStep (P : SteeringProtocol d r s df rf sf m) (rho : AssemblageBlocks d r s)
    (w : Fin m) : AssemblageBlocks df r s :=
  fun a x => P.kraus w * rho a x * (P.kraus w).conjTranspose

/-- PDF p.9, Steps 2–3: sample x using (w,xf), query Alice's device,
then sample af using (w,a,x,xf), retaining Bob's branch w. -/
def protocolBranch (P : SteeringProtocol d r s df rf sf m) (rho : AssemblageBlocks d r s)
    (w : Fin m) : AssemblageBlocks df rf sf :=
  fun af xf => ∑ x, ∑ a,
    (P.output w af a x xf : ℂ) • ((P.input w xf x : ℂ) • bobStep P rho w a x)

/-- PDF p.3, (9): the operational map with the outcome forgotten. -/
def protocolAction (P : SteeringProtocol d r s df rf sf m) (rho : AssemblageBlocks d r s) :
    AssemblageBlocks df rf sf := fun af xf => ∑ w, protocolBranch P rho w af xf

/-- PDF p.10, (A6), equivalently (9)–(10), written as a closed block formula. -/
def oneWayFormula (P : SteeringProtocol d r s df rf sf m) (rho : AssemblageBlocks d r s) :
    AssemblageBlocks df rf sf :=
  fun af xf => ∑ w, ∑ a, ∑ x,
    ((P.input w xf x * P.output w af a x xf : ℝ) : ℂ) •
      (P.kraus w * rho a x * (P.kraus w).conjTranspose)

/-- PDF p.3, (8): Bob's branch success probability, independent of Alice's input. -/
def steeringBranchProbability (P : SteeringProtocol d r s df rf sf m) (rho : AssemblageBlocks d r s)
    (w : Fin m) : ℝ :=
  (P.kraus w * bobTotal rho * (P.kraus w).conjTranspose).trace.re

/-- PDF p.4: conditional normalization. All physical uses require positive assemblageMass;
at zero this raw array is zero and is never called a normalized assemblage. -/
def normalizeAssemblage (rho : AssemblageBlocks d r s) : AssemblageBlocks d r s :=
  fun a x => ((assemblageMass rho : ℂ)⁻¹) • rho a x

/-- PDF p.7 and Appendices E–F: exact single-copy conversion with strictly
positive input-independent success. All allowed wirings are quantified. -/
def SteeringConverts (rho : AssemblageBlocks d r s) (sigma : AssemblageBlocks df rf sf) : Prop :=
  IsNormalized rho ∧ IsNormalized sigma ∧
  ∃ (m : ℕ) (P : SteeringProtocol d r s df rf sf m) (p : ℝ),
    0 < p ∧ p ≤ 1 ∧ ∀ af xf, protocolAction P rho af xf = (p : ℂ) • sigma af xf

/-- PDF p.5, (14): a binary affine combination, also used with nu>1 in (21). -/
def assemblageMixture (rho sigma : AssemblageBlocks d r s) (t : ℝ) : AssemblageBlocks d r s :=
  fun a x => (t : ℂ) • rho a x + ((1 - t : ℝ) : ℂ) • sigma a x

/-- PDF p.11, (C2): a finite weighted assemblageMixture with the same index set throughout. -/
def assemblageFamilyMixture {n : ℕ} (p : Fin n → ℝ) (rho : Fin n → AssemblageBlocks d r s) :
    AssemblageBlocks d r s := fun a x => ∑ j, (p j : ℂ) • rho j a x

/-- PDF p.5, Definition 2: a measure for every finite dimension and alphabet. -/
abbrev SteeringMeasure := (d r s : ℕ) → AssemblageBlocks d r s → ℝ

/-- PDF p.5, (13): zero-probability outcomes contribute zero, without assigning
a physical conditional assemblage to an impossible event. -/
def measureBranch (S : SteeringMeasure) (P : SteeringProtocol d r s df rf sf m)
    (rho : AssemblageBlocks d r s) (w : Fin m) : ℝ :=
  let p := steeringBranchProbability P rho w
  if 0 < p then p * S df rf sf (normalizeAssemblage (protocolBranch P rho w)) else 0

/-- PDF p.5, Definition 2(ii): average monotonicity for complete instruments,
with the outcome available to Alice before the average is taken. -/
def StrongMonotonicity (S : SteeringMeasure) : Prop :=
  ∀ (d r s df rf sf m : ℕ) (P : SteeringProtocol d r s df rf sf m)
    (rho : AssemblageBlocks d r s), IsNormalized rho → P.Complete →
    ∑ w, measureBranch S P rho w ≤ S d r s rho

/-- PDF p.5, Definition 2(iii): convexity, including both endpoints. -/
def MeasureConvexity (S : SteeringMeasure) : Prop :=
  ∀ (d r s : ℕ) (rho sigma : AssemblageBlocks d r s) (t : ℝ),
    IsNormalized rho → IsNormalized sigma → 0 ≤ t → t ≤ 1 →
    S d r s (assemblageMixture rho sigma t) ≤ t * S d r s rho + (1 - t) * S d r s sigma

/-- PDF p.5, Definition 2: nonnegative, vanishes on LHS, average monotone,
and convex. There is no converse assertion in the vanishing clause. -/
def IsConvexSteeringMonotone (S : SteeringMeasure) : Prop :=
  (∀ d r s (rho : AssemblageBlocks d r s), IsNormalized rho → 0 ≤ S d r s rho) ∧
  (∀ d r s (rho : AssemblageBlocks d r s), IsNormalized rho → HasFiniteLHSModel rho → S d r s rho = 0) ∧
  StrongMonotonicity S ∧ MeasureConvexity S

/-- PDF p.2, (3): Alice's classical outcome flag, with Bob's actual block. -/
def assemblageCQMatrix (rho : AssemblageBlocks d r s) (x : Fin s) :
    Matrix (Fin r × Fin d) (Fin r × Fin d) ℂ :=
  fun i j => if i.1 = j.1 then rho i.1 x i.2 j.2 else 0

/-- PDF p.5, (15): the pinned Umegaki divergence, with natural logarithms
and +infinity on support failure. The non-density branch is a totalization
outside the physical domain; it is never a finite substitute for divergence. -/
def densityD {I : Type*} [Fintype I] [DecidableEq I]
    (rho sigma : Matrix I I ℂ) : ℝ≥0∞ :=
  if h : QAlgorithms.IsDensityMatrix rho ∧ QAlgorithms.IsDensityMatrix sigma then
    QAlgorithms.QuantumLoss.umegakiEntropy rho sigma h.1 h.2
  else ⊤

/-- PDF p.5, (16): KL divergence of two normalized finite laws, computed as
the Umegaki divergence of their diagonal density matrices. This retains negative
individual log-ratio contributions before forming the nonnegative total. -/
def steeringClassicalDivergence {n : ℕ} (p q : Fin n → ℝ) : ℝ≥0∞ :=
  densityD (Matrix.diagonal (fun i => (p i : ℂ)))
    (Matrix.diagonal (fun i => (q i : ℂ)))

/-- PDF p.6, (18)–(19): arbitrary finite-outcome complete Bob instrument and
an input law chosen after its outcome. The number of outcomes is unrestricted. -/
structure EntropyStrategy (d s : ℕ) where
  outcomes : ℕ
  kraus : Fin outcomes → Matrix (Fin d) (Fin d) ℂ
  complete : ∑ g, (kraus g).conjTranspose * kraus g = 1
  input : Fin outcomes → Fin s → ℝ
  input_law : ∀ g, IsLaw (input g)

/-- PDF p.6, (19): the two input assemblages give separately computed laws. -/
def strategyProbability (T : EntropyStrategy d s) (rho : AssemblageBlocks d r s)
    (g : Fin T.outcomes) : ℝ :=
  (T.kraus g * bobTotal rho * (T.kraus g).conjTranspose).trace.re

/-- PDF p.6, numerator in (18): Bob's filtered flagged state at fixed input. -/
def filteredCQ (T : EntropyStrategy d s) (rho : AssemblageBlocks d r s)
    (g : Fin T.outcomes) (x : Fin s) : Matrix (Fin r × Fin d) (Fin r × Fin d) ℂ :=
  assemblageCQMatrix (fun a y => T.kraus g * rho a y * (T.kraus g).conjTranspose) x

/-- PDF p.6, (18): p D(A/p || B/q), zero when p=0 and infinite when p>0=q.
This homogeneous convention avoids a fictitious zero-over-zero density matrix. -/
def conditionalD {I : Type*} [Fintype I] [DecidableEq I]
    (p q : ℝ) (A B : Matrix I I ℂ) : ℝ≥0∞ :=
  if 0 < p then
    if 0 < q then ENNReal.ofReal p * densityD (((p : ℂ)⁻¹) • A) (((q : ℂ)⁻¹) • B)
    else ⊤
  else 0

/-- PDF p.6, (18): both the classical outcome KL term and the conditional
quantum term weighted by P(x|g) P_g, with separate denominators P_g and P'_g. -/
def entropyScore (T : EntropyStrategy d s) (rho sigma : AssemblageBlocks d r s) : ℝ≥0∞ :=
  steeringClassicalDivergence (strategyProbability T rho) (strategyProbability T sigma) +
    ∑ g, ∑ x, ENNReal.ofReal (T.input g x) *
      conditionalD (strategyProbability T rho g) (strategyProbability T sigma g)
        (filteredCQ T rho g x) (filteredCQ T sigma g x)

/-- PDF p.6, Definition 3: complete-order optimization, allowing +infinity.
The source's maximum attainment is asserted in relativeEntropy_convexMonotone;
it is not assumed or hidden by replacing the range with a bounded real set. -/
def assemblageEntropy (rho sigma : AssemblageBlocks d r s) : ℝ≥0∞ :=
  sSup {v | ∃ T : EntropyStrategy d s, entropyScore T rho sigma = v}

/-- PDF p.6, (18): the explicit maximum-attainment obligation. -/
def EntropyMaximum (rho sigma : AssemblageBlocks d r s) : Prop :=
  ∃ T : EntropyStrategy d s, entropyScore T rho sigma = assemblageEntropy rho sigma

/-- PDF p.6, (20): optimize over normalized LHS assemblages on the same space
and alphabets. Minimum attainment and finiteness are conclusions below. -/
def steeringEntropy (rho : AssemblageBlocks d r s) : ℝ≥0∞ :=
  sInf {v | ∃ sigma : AssemblageBlocks d r s,
    IsNormalized sigma ∧ HasFiniteLHSModel sigma ∧ assemblageEntropy rho sigma = v}

/-- PDF p.6, Theorem 2: real-valued relative entropy of steering. Its finiteness
on every normalized finite assemblage is an explicit theorem obligation. -/
def relativeEntropyMeasure : SteeringMeasure :=
  fun _ _ _ rho => (steeringEntropy rho).toReal

/-- PDF p.11, (C1): first-input probability weights, with each output separately
normalized. If the second branch is impossible and the first is not, the term
is infinite; no support-inclusion hypothesis has been added. -/
def entropyBranch (P : SteeringProtocol d r s df rf sf m) (rho sigma : AssemblageBlocks d r s)
    (w : Fin m) : ℝ≥0∞ :=
  let p := steeringBranchProbability P rho w
  let q := steeringBranchProbability P sigma w
  if 0 < p then
    if 0 < q then ENNReal.ofReal p *
      assemblageEntropy (normalizeAssemblage (protocolBranch P rho w))
        (normalizeAssemblage (protocolBranch P sigma w))
    else ⊤
  else 0

/-- PDF p.6, Definition 4: all raw nu>=0 candidates are retained, including
nu>1. Feasible noise and free assemblages are normalized and no-signaling. -/
def WeightCandidate (rho : AssemblageBlocks d r s) (nu : ℝ) : Prop :=
  0 ≤ nu ∧ ∃ noise sigma : AssemblageBlocks d r s,
    IsNormalized noise ∧ IsNormalized sigma ∧ HasFiniteLHSModel sigma ∧ rho = assemblageMixture noise sigma nu

/-- PDF p.6, Definition 4: the real infimum; the theorem explicitly asserts
that it is an attained minimum, which rules out junk empty/unbounded values. -/
def steerableWeight (rho : AssemblageBlocks d r s) : ℝ := sInf {nu | WeightCandidate rho nu}

/-- PDF p.6, Definition 5: arbitrary assemblage noise, not necessarily LHS. -/
def RobustnessCandidate (rho : AssemblageBlocks d r s) (nu : ℝ) : Prop :=
  0 ≤ nu ∧ ∃ noise : AssemblageBlocks d r s, IsNormalized noise ∧
    HasFiniteLHSModel (fun a x => (((1 + nu : ℝ) : ℂ)⁻¹) • rho a x +
      ((nu / (1 + nu) : ℝ) : ℂ) • noise a x)

/-- PDF p.6, Definition 5: the attained minimum is asserted in Theorem 3. -/
def steeringRobustness (rho : AssemblageBlocks d r s) : ℝ :=
  sInf {nu | RobustnessCandidate rho nu}

/-- PDF p.6, (23): pure states with no chosen global phase. -/
def conditionalKetProjector {d : ℕ} (v : Fin d → ℂ) : Matrix (Fin d) (Fin d) ℂ :=
  fun i j => v i * star (v j)

/-- PDF p.6, (23): every conditional state has a unit-vector representative,
including zero-probability events; purity does not require a pure Bob marginal. -/
def IsPureAssemblage (rho : AssemblageBlocks d r s) : Prop :=
  ∃ (p : Fin s → Fin r → ℝ) (v : Fin r → Fin s → Fin d → ℂ),
    (∀ x, IsLaw (p x)) ∧ (∀ a x, ∑ i, Complex.normSq (v a x i) = 1) ∧
    ∀ a x, rho a x = (p x a : ℂ) • conditionalKetProjector (v a x)

/-- PDF p.6, following (23): an orthonormal family at each input. -/
def IsPureOrthogonal (rho : AssemblageBlocks d r s) : Prop :=
  ∃ (p : Fin s → Fin r → ℝ) (v : Fin r → Fin s → Fin d → ℂ),
    (∀ x, IsLaw (p x)) ∧
    (∀ a b x, ∑ i, star (v a x i) * v b x i = if a = b then 1 else 0) ∧
    ∀ a x, rho a x = (p x a : ℂ) • conditionalKetProjector (v a x)

/-- PDF p.15, (F5): the two qubit bases used for the infinite-class witnesses. -/
def angleKet (theta : ℝ) (a x : Fin 2) : Fin 2 → ℂ :=
  if x = 0 then (fun i => if i = a then 1 else 0)
  else if a = 0 then (fun i => if i = 0 then (Real.cos theta : ℂ) else (Real.sin theta : ℂ))
  else (fun i => if i = 0 then -(Real.sin theta : ℂ) else (Real.cos theta : ℂ))

/-- PDF p.15, (F5): each rank-one block has probability one half. -/
def angleAssemblage (theta : ℝ) : AssemblageBlocks 2 2 2 :=
  fun a x => (1 / 2 : ℂ) • conditionalKetProjector (angleKet theta a x)

/-! Sanity checks for the new local interfaces. These are definition checks,
not proofs of the paper's theorems. -/

example : IsLaw (fun _ : Fin 1 => (1 : ℝ)) := by
  simp [IsLaw]

example (rho sigma : AssemblageBlocks d r s) : assemblageMixture rho sigma 0 = sigma := by
  funext a x
  simp [assemblageMixture]

example (rho sigma : AssemblageBlocks d r s) : assemblageMixture rho sigma 1 = rho := by
  funext a x
  simp [assemblageMixture]

example {I : Type*} [Fintype I] [DecidableEq I]
    (A B : Matrix I I ℂ) (q : ℝ) : conditionalD 0 q A B = 0 := by
  simp [conditionalD]

example (P : SteeringProtocol d r s df rf sf 0) (rho : AssemblageBlocks d r s) :
    protocolAction P rho = 0 := by
  funext a x
  simp [protocolAction]

example (theta : ℝ) : angleKet theta 0 0 0 = 1 := by
  simp [angleKet]

example (theta : ℝ) : angleKet theta 0 0 1 = 0 := by
  simp [angleKet]

end QAlgorithms.Papers.Independent.Intakethird100.P36
Source
The resource theory of steering, arXiv:1409.5804v3, https://arxiv.org/abs/1409.5804v3
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P36.AssemblageBlocks

For natural numbers d,r,sd,r,sd,r,s, this denotes an arbitrary family ρ=(ρa,x)\rho=(\rho_{a,x})ρ=(ρa,x​) of complex d×dd\times dd×d matrices, indexed by 0≤a<r0\le a<r0≤a<r and 0≤x<s0\le x<s0≤x<s. No positivity, normalization, or compatibility conditions are imposed. All three parameters may be zero; if r=0r=0r=0 or s=0s=0s=0, there are no blocks, and if d=0d=0d=0, each block is an empty matrix.

QAlgorithms.Papers.Independent.Intakethird100.P36.EntropyMaximum

For arbitrary natural numbers d,r,sd,r,sd,r,s and arbitrary complex d×dd\times dd×d matrix families ρ,σ\rho,\sigmaρ,σ indexed by a<r,x<sa<r,x<sa<r,x<s, this predicate says that there exists a strategy whose score equals the supremum of all strategy scores. A strategy consists of a natural count hhh, matrices KgK_gKg​ satisfying ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​, and real weights ug,x≥0u_{g,x}\ge0ug,x​≥0 satisfying ∑xug,x=1\sum_xu_{g,x}=1∑x​ug,x​=1. For each such strategy, put Bτ=s−1∑x,aτa,xB_\tau=s^{-1}\sum_{x,a}\tau_{a,x}Bτ​=s−1∑x,a​τa,x​, pg=Re⁡Tr⁡(KgBρKg†)p_g=\operatorname{Re}\operatorname{Tr}(K_gB_\rho K_g^\dagger)pg​=ReTr(Kg​Bρ​Kg†​), qg=Re⁡Tr⁡(KgBσKg†)q_g=\operatorname{Re}\operatorname{Tr}(K_gB_\sigma K_g^\dagger)qg​=ReTr(Kg​Bσ​Kg†​), and Fg,xτ=⨁aKgτa,xKg†F^\tau_{g,x}=\bigoplus_aK_g\tau_{a,x}K_g^\daggerFg,xτ​=⨁a​Kg​τa,x​Kg†​; its score is

D(diag⁡p,diag⁡q)+∑g,xug,x C(pg,qg;Fg,xρ,Fg,xσ).\mathcal D(\operatorname{diag}p,\operatorname{diag}q) +\sum_{g,x}u_{g,x}\,C(p_g,q_g;F^\rho_{g,x},F^\sigma_{g,x}).D(diagp,diagq)+g,x∑​ug,x​C(pg​,qg​;Fg,xρ​,Fg,xσ​).

The divergence D(R,S)\mathcal D(R,S)D(R,S) is infinite unless R,SR,SR,S are positive semidefinite with trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; when these conditions hold it is max⁡{0,Re⁡Tr⁡(f(R)−Rlog⁡S)}\max\{0,\operatorname{Re}\operatorname{Tr}(f(R)-R\log S)\}max{0,ReTr(f(R)−RlogS)}, with f(t)=tlog⁡tf(t)=t\log tf(t)=tlogt and continuous functional calculus. The conditional term CCC is 000 for p≤0p\le0p≤0, +∞+\infty+∞ for p>0,q≤0p>0,q\le0p>0,q≤0, and pD(A/p,B/q)p\mathcal D(A/p,B/q)pD(A/p,B/q) for p,q>0p,q>0p,q>0. No assemblage assumptions are imposed, and the attained supremum may be infinite. If there is no strategy, the supremum is 000, but this existence predicate is false. Zero dimensions and empty sums are included with total reciprocal and extended arithmetic conventions.

QAlgorithms.Papers.Independent.Intakethird100.P36.EntropyStrategy

For arbitrary natural numbers d,sd,sd,s, an object of this structure consists of a natural number hhh of outcomes, complex d×dd\times dd×d matrices KgK_gKg​ for 0≤g<h0\le g<h0≤g<h, and real numbers ug,xu_{g,x}ug,x​ for 0≤g<h0\le g<h0≤g<h and 0≤x<s0\le x<s0≤x<s, together with the requirements

∑g=0h−1Kg†Kg=Id,ug,x≥0,∑x=0s−1ug,x=1for every g.\sum_{g=0}^{h-1}K_g^\dagger K_g=I_d, \qquad u_{g,x}\ge0, \qquad \sum_{x=0}^{s-1}u_{g,x}=1\quad\text{for every }g.g=0∑h−1​Kg†​Kg​=Id​,ug,x​≥0,x=0∑s−1​ug,x​=1for every g.

Here †\dagger† denotes conjugate transpose. None of d,s,hd,s,hd,s,h is explicitly required to be positive. If h=0h=0h=0, completeness is possible only when d=0d=0d=0; if s=0s=0s=0, the input laws are impossible for any existing outcome, so hhh must be 000. In particular, d=s=h=0d=s=h=0d=s=h=0 is permitted.

QAlgorithms.Papers.Independent.Intakethird100.P36.HasFiniteLHSModel

For natural numbers d,r,sd,r,sd,r,s and a family ρa,x\rho_{a,x}ρa,x​ of complex d×dd\times dd×d matrices, this predicate first requires d,r,s>0d,r,s>0d,r,s>0, positive semidefinite blocks, setting-independent sums ∑aρa,x\sum_a\rho_{a,x}∑a​ρa,x​, and total mass mρ=Re⁡tr⁡(s−1∑x∑aρa,x)m_\rho=\operatorname{Re}\operatorname{tr}(s^{-1}\sum_x\sum_a\rho_{a,x})mρ​=Retr(s−1∑x​∑a​ρa,x​) in [0,1][0,1][0,1]; it additionally requires the existence of a natural number nnn, real numbers pℓp_\ellpℓ​, complex d×dd\times dd×d matrices ξℓ\xi_\ellξℓ​, and real numbers Rℓ,x,aR_{\ell,x,a}Rℓ,x,a​, indexed by 0≤ℓ<n0\le\ell<n0≤ℓ<n, 0≤x<s0\le x<s0≤x<s, 0≤a<r0\le a<r0≤a<r, such that

pℓ≥0,∑ℓpℓ=1,ξℓ⪰0,tr⁡ξℓ=1,Rℓ,x,a≥0,∑aRℓ,x,a=1,ρa,x=∑ℓ=0n−1mρpℓRℓ,x,a ξℓ.p_\ell\ge0,\quad \sum_\ell p_\ell=1,\qquad \xi_\ell\succeq0,\quad\operatorname{tr}\xi_\ell=1,\qquad R_{\ell,x,a}\ge0,\quad\sum_aR_{\ell,x,a}=1, \qquad \rho_{a,x}=\sum_{\ell=0}^{n-1}m_\rho p_\ell R_{\ell,x,a}\,\xi_\ell.pℓ​≥0,ℓ∑​pℓ​=1,ξℓ​⪰0,trξℓ​=1,Rℓ,x,a​≥0,a∑​Rℓ,x,a​=1,ρa,x​=ℓ=0∑n−1​mρ​pℓ​Rℓ,x,a​ξℓ​.

Although n=0n=0n=0 is quantified over, it cannot satisfy ∑ℓpℓ=1\sum_\ell p_\ell=1∑ℓ​pℓ​=1. Mass zero is allowed and makes every block in the displayed representation zero; mass one is not required.

QAlgorithms.Papers.Independent.Intakethird100.P36.IsAssemblage

For natural numbers d,r,sd,r,sd,r,s and a family ρa,x\rho_{a,x}ρa,x​ of complex d×dd\times dd×d matrices indexed by 0≤a<r0\le a<r0≤a<r and 0≤x<s0\le x<s0≤x<s, this predicate requires d,r,s>0d,r,s>0d,r,s>0, every ρa,x\rho_{a,x}ρa,x​ to be positive semidefinite, and

∑a=0r−1ρa,x=∑a=0r−1ρa,yfor every 0≤x,y<s,0≤Re⁡tr⁡(1s∑x=0s−1∑a=0r−1ρa,x)≤1.\sum_{a=0}^{r-1}\rho_{a,x}=\sum_{a=0}^{r-1}\rho_{a,y} \quad\text{for every }0\le x,y<s, \qquad 0\le \operatorname{Re}\operatorname{tr}\left( \frac1s\sum_{x=0}^{s-1}\sum_{a=0}^{r-1}\rho_{a,x}\right)\le1.a=0∑r−1​ρa,x​=a=0∑r−1​ρa,y​for every 0≤x,y<s,0≤Retr(s1​x=0∑s−1​a=0∑r−1​ρa,x​)≤1.

Thus zero total mass is allowed, while a zero value of any dimension makes the predicate false.

QAlgorithms.Papers.Independent.Intakethird100.P36.IsConvexSteeringMonotone

For an arbitrary function SSS assigning a real number to every array ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ of complex d×dd\times dd×d matrices for every d,r,s∈Nd,r,s\in\mathbb Nd,r,s∈N, this predicate is the conjunction of four requirements. In these requirements, an array is normalized when d,r,s>0d,r,s>0d,r,s>0, every ρax\rho_{ax}ρax​ is positive semidefinite, ∑aρax\sum_a\rho_{ax}∑a​ρax​ is independent of xxx, and its mass Mρ=Re⁡tr⁡(s−1∑x∑aρax)M_\rho=\operatorname{Re}\operatorname{tr}(s^{-1}\sum_x\sum_a\rho_{ax})Mρ​=Retr(s−1∑x​∑a​ρax​) lies in [0,1][0,1][0,1] and equals 111. First, S(d,r,s,ρ)≥0S(d,r,s,\rho)\ge0S(d,r,s,ρ)≥0 for every normalized array. Second, S(d,r,s,ρ)=0S(d,r,s,\rho)=0S(d,r,s,ρ)=0 for every normalized array that admits some finite representation

ρax=∑ℓ<nMρ pℓ hℓ,x,a ξℓ,\rho_{ax}=\sum_{\ell<n}M_\rho\,p_\ell\,h_{\ell,x,a}\,\xi_\ell,ρax​=ℓ<n∑​Mρ​pℓ​hℓ,x,a​ξℓ​,

where n∈Nn\in\mathbb Nn∈N, pℓ≥0p_\ell\ge0pℓ​≥0, ∑ℓpℓ=1\sum_\ell p_\ell=1∑ℓ​pℓ​=1, every ξℓ\xi_\ellξℓ​ is a positive semidefinite complex d×dd\times dd×d matrix with trace 111, and hℓ,x,a≥0h_{\ell,x,a}\ge0hℓ,x,a​≥0, ∑a<rhℓ,x,a=1\sum_{a<r}h_{\ell,x,a}=1∑a<r​hℓ,x,a​=1 for every ℓ,x\ell,xℓ,x; the finite-model condition also requires positive dimensions, positive semidefinite blocks, setting-independent outcome sums, and mass in [0,1][0,1][0,1]. Third, for every normalized ρ\rhoρ, every positive output dimensions df,rf,sfd_f,r_f,s_fdf​,rf​,sf​, every natural mmm, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d with ∑w<mKw†Kw=Id\sum_{w<m}K_w^\dagger K_w=I_d∑w<m​Kw†​Kw​=Id​ and Id−∑wKw†Kw⪰0I_d-\sum_wK_w^\dagger K_w\succeq0Id​−∑w​Kw†​Kw​⪰0, and nonnegative arrays uw,xf,x,vw,af,a,x,xfu_{w,x_f,x},v_{w,a_f,a,x,x_f}uw,xf​,x​,vw,af​,a,x,xf​​ summing to 111 over x<sx<sx<s and af<rfa_f<r_faf​<rf​, respectively, put

pw=Re⁡tr⁡ ⁣(Kw(1s∑x∑aρax)Kw†),τafxf(w)=∑x∑avw,af,a,x,xfuw,xf,xKwρaxKw†,Mw=Re⁡tr⁡ ⁣(1sf∑xf∑afτafxf(w));p_w=\operatorname{Re}\operatorname{tr}\!\left(K_w\left(\frac1s\sum_x\sum_a\rho_{ax}\right)K_w^\dagger\right),\qquad \tau^{(w)}_{a_fx_f}=\sum_x\sum_a v_{w,a_f,a,x,x_f}u_{w,x_f,x}K_w\rho_{ax}K_w^\dagger,\qquad M_w=\operatorname{Re}\operatorname{tr}\!\left(\frac1{s_f}\sum_{x_f}\sum_{a_f}\tau^{(w)}_{a_fx_f}\right);pw​=Retr(Kw​(s1​x∑​a∑​ρax​)Kw†​),τaf​xf​(w)​=x∑​a∑​vw,af​,a,x,xf​​uw,xf​,x​Kw​ρax​Kw†​,Mw​=Retr​sf​1​xf​∑​af​∑​τaf​xf​(w)​​;

then the sum over w<mw<mw<m of pwS(df,rf,sf,Mw−1τ(w))p_wS(d_f,r_f,s_f,M_w^{-1}\tau^{(w)})pw​S(df​,rf​,sf​,Mw−1​τ(w)) when pw>0p_w>0pw​>0, and 000 otherwise, must be at most S(d,r,s,ρ)S(d,r,s,\rho)S(d,r,s,ρ), with entrywise rescaling and 0−1=00^{-1}=00−1=0. Fourth, for every pair of normalized arrays ρ,σ\rho,\sigmaρ,σ of the same dimensions and every t∈[0,1]t\in[0,1]t∈[0,1],

S(d,r,s,tρ+(1−t)σ)≤tS(d,r,s,ρ)+(1−t)S(d,r,s,σ).S(d,r,s,t\rho+(1-t)\sigma) \le tS(d,r,s,\rho)+(1-t)S(d,r,s,\sigma).S(d,r,s,tρ+(1−t)σ)≤tS(d,r,s,ρ)+(1−t)S(d,r,s,σ).

All natural dimensions are quantified over, but normalization and protocol hypotheses exclude zero dimensions from these obligations; there is no restriction on the values of SSS on other arrays.

QAlgorithms.Papers.Independent.Intakethird100.P36.IsNormalized

For natural numbers d,r,sd,r,sd,r,s and a family ρa,x\rho_{a,x}ρa,x​ of complex d×dd\times dd×d matrices indexed by 0≤a<r0\le a<r0≤a<r and 0≤x<s0\le x<s0≤x<s, this predicate requires d,r,s>0d,r,s>0d,r,s>0, positive semidefiniteness of every block, equality of ∑aρa,x\sum_a\rho_{a,x}∑a​ρa,x​ for every pair of settings x,yx,yx,y, and, writing

mρ=Re⁡tr⁡(1s∑x=0s−1∑a=0r−1ρa,x),m_\rho=\operatorname{Re}\operatorname{tr}\left( \frac1s\sum_{x=0}^{s-1}\sum_{a=0}^{r-1}\rho_{a,x}\right),mρ​=Retr(s1​x=0∑s−1​a=0∑r−1​ρa,x​),

both 0≤mρ≤10\le m_\rho\le10≤mρ​≤1 and mρ=1m_\rho=1mρ​=1. It cannot hold when any dimension is zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.IsPureAssemblage

For natural numbers d,r,sd,r,sd,r,s and a family ρ=(ρa,x)a<r, x<s\rho=(\rho_{a,x})_{a<r,\ x<s}ρ=(ρa,x​)a<r, x<s​ of complex d×dd\times dd×d matrices, this predicate means that there exist real numbers px(a)p_x(a)px​(a) and complex vectors va,x∈Cdv_{a,x}\in\mathbb C^dva,x​∈Cd such that, for every setting x<sx<sx<s, px(a)≥0p_x(a)\geq0px​(a)≥0 for every a<ra<ra<r and ∑apx(a)=1\sum_a p_x(a)=1∑a​px​(a)=1; for every a,xa,xa,x, ∑i<d∣va,x,i∣2=1\sum_{i<d}|v_{a,x,i}|^2=1∑i<d​∣va,x,i​∣2=1; and, for every a,x,i,ja,x,i,ja,x,i,j,

(ρa,x)ij=px(a)va,x,iva,x,j‾.(\rho_{a,x})_{ij}=p_x(a)v_{a,x,i}\overline{v_{a,x,j}}.(ρa,x​)ij​=px​(a)va,x,i​va,x,j​​.

The vector normalization is required even when px(a)=0p_x(a)=0px​(a)=0. No condition equating the outcome sums for different settings is included. If s=0s=0s=0, all the indexed requirements are vacuous; if s>0s>0s>0 and r=0r=0r=0, the probability requirement is impossible, and if s,r>0s,r>0s,r>0 and d=0d=0d=0, vector normalization is impossible.

QAlgorithms.Papers.Independent.Intakethird100.P36.IsPureOrthogonal

For natural numbers d,r,sd,r,sd,r,s and a family ρ=(ρa,x)a<r, x<s\rho=(\rho_{a,x})_{a<r,\ x<s}ρ=(ρa,x​)a<r, x<s​ of complex d×dd\times dd×d matrices, this predicate means that there exist real probabilities px(a)≥0p_x(a)\geq0px​(a)≥0 with ∑a<rpx(a)=1\sum_{a<r}p_x(a)=1∑a<r​px​(a)=1 for every x<sx<sx<s, and complex vectors va,x∈Cdv_{a,x}\in\mathbb C^dva,x​∈Cd, such that for every a,b<ra,b<ra,b<r and x<sx<sx<s,

∑i<dva,x,i‾vb,x,i={1,a=b,0,a≠b,(ρa,x)ij=px(a)va,x,iva,x,j‾(i,j<d).\sum_{i<d}\overline{v_{a,x,i}}v_{b,x,i} =\begin{cases}1,&a=b,\\0,&a\ne b,\end{cases} \qquad (\rho_{a,x})_{ij}=p_x(a)v_{a,x,i}\overline{v_{a,x,j}} \quad(i,j<d).i<d∑​va,x,i​​vb,x,i​={1,0,​a=b,a=b,​(ρa,x​)ij​=px​(a)va,x,i​va,x,j​​(i,j<d).

Thus all the outcome vectors at each setting are required to be orthonormal, including vectors assigned zero probability. No condition equating the outcome sums between settings is included. For s=0s=0s=0 the requirements are vacuous; for s>0s>0s>0, they are impossible if r=0r=0r=0, and otherwise require an orthonormal family of rrr vectors in Cd\mathbb C^dCd.

QAlgorithms.Papers.Independent.Intakethird100.P36.MeasureConvexity

For an arbitrary real-valued function SSS on arrays of complex matrices of all natural-number dimensions, this predicate requires that, for every d,r,s∈Nd,r,s\in\mathbb Nd,r,s∈N, every two arrays ρ,σ\rho,\sigmaρ,σ of complex d×dd\times dd×d matrices indexed by a<r,x<sa<r,x<sa<r,x<s, and every real ttt, if both arrays are normalized and 0≤t≤10\le t\le10≤t≤1, then

S ⁣(d,r,s,(tρax+(1−t)σax)a,x)≤t S(d,r,s,ρ)+(1−t) S(d,r,s,σ).S\!\left(d,r,s,(t\rho_{ax}+(1-t)\sigma_{ax})_{a,x}\right) \le t\,S(d,r,s,\rho)+(1-t)\,S(d,r,s,\sigma).S(d,r,s,(tρax​+(1−t)σax​)a,x​)≤tS(d,r,s,ρ)+(1−t)S(d,r,s,σ).

For each array, normalized means d,r,s>0d,r,s>0d,r,s>0, all its matrices are positive semidefinite, its sums over aaa are independent of xxx, and the real trace of s−1∑x∑aρaxs^{-1}\sum_x\sum_a\rho_{ax}s−1∑x​∑a​ρax​, or of the corresponding expression for σ\sigmaσ, lies in [0,1][0,1][0,1] and equals 111. The endpoints t=0,1t=0,1t=0,1 are included. No condition is imposed by this predicate when either array is unnormalized, when a dimension is zero, or when t∉[0,1]t\notin[0,1]t∈/[0,1].

QAlgorithms.Papers.Independent.Intakethird100.P36.RobustnessCandidate

For natural numbers d,r,sd,r,sd,r,s, an arbitrary family ρ=(ρa,x)a<r, x<s\rho=(\rho_{a,x})_{a<r,\ x<s}ρ=(ρa,x​)a<r, x<s​ of complex d×dd\times dd×d matrices, and ν∈R\nu\in\mathbb Rν∈R, this predicate holds precisely when ν≥0\nu\geq0ν≥0 and there exists a normalized family η\etaη such that the family

Ya,x=ρa,x1+ν+ν1+νηa,xY_{a,x}=\frac{\rho_{a,x}}{1+\nu}+\frac{\nu}{1+\nu}\eta_{a,x}Ya,x​=1+νρa,x​​+1+νν​ηa,x​

has a finite hidden-state model. Explicitly, normalization of η\etaη requires d,r,s>0d,r,s>0d,r,s>0, positive semidefinite blocks, setting-independent sums ∑aηa,x\sum_a\eta_{a,x}∑a​ηa,x​, and mass 111, where the mass of any family FFF is μ(F)=Re⁡tr⁡(s−1∑x∑aFa,x)\mu(F)=\operatorname{Re}\operatorname{tr}(s^{-1}\sum_x\sum_aF_{a,x})μ(F)=Retr(s−1∑x​∑a​Fa,x​). The condition on YYY requires positive semidefinite blocks, setting-independent outcome sums, 0≤μ(Y)≤10\leq\mu(Y)\leq10≤μ(Y)≤1, and the existence of a natural number hhh, probabilities pℓp_\ellpℓ​ summing to 111, positive semidefinite matrices ξℓ\xi_\ellξℓ​ of trace 111, and probability distributions qℓ,xq_{\ell,x}qℓ,x​ on the rrr outcomes such that Ya,x=μ(Y)∑ℓ<hpℓqℓ,x(a)ξℓY_{a,x}=\mu(Y)\sum_{\ell<h}p_\ell q_{\ell,x}(a)\xi_\ellYa,x​=μ(Y)∑ℓ<h​pℓ​qℓ,x​(a)ξℓ​ for all a,xa,xa,x. No normalization or assemblage condition is separately imposed on ρ\rhoρ, and ν\nuν has no upper bound. The denominator 1+ν1+\nu1+ν is positive for a candidate; if any of d,r,sd,r,sd,r,s is zero, the predicate cannot hold.

QAlgorithms.Papers.Independent.Intakethird100.P36.SteeringConverts

For natural numbers d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, and arrays ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ of complex d×dd\times dd×d matrices and σ=(σafxf)af<rf,xf<sf\sigma=(\sigma_{a_fx_f})_{a_f<r_f,x_f<s_f}σ=(σaf​xf​​)af​<rf​,xf​<sf​​ of complex df×dfd_f\times d_fdf​×df​ matrices, this predicate holds precisely when both arrays are normalized and the following protocol exists. Here normalized means that the three dimensions of the array are positive, every matrix entry of the array is positive semidefinite, its sums over outcomes are independent of the setting, and its mass—defined as the real trace of the average of those sums—is in [0,1][0,1][0,1] and equals 111. There must exist a natural number mmm, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d for w<mw<mw<m satisfying Id−∑wKw†Kw⪰0I_d-\sum_wK_w^\dagger K_w\succeq0Id​−∑w​Kw†​Kw​⪰0, nonnegative real numbers uw,xf,xu_{w,x_f,x}uw,xf​,x​ summing to 111 over x<sx<sx<s for every w,xfw,x_fw,xf​, nonnegative real numbers vw,af,a,x,xfv_{w,a_f,a,x,x_f}vw,af​,a,x,xf​​ summing to 111 over af<rfa_f<r_faf​<rf​ for every w,a,x,xfw,a,x,x_fw,a,x,xf​, and a real number ppp with 0<p≤10<p\le10<p≤1, such that

∀af<rf, ∀xf<sf,∑w<m∑x<s∑a<rvw,af,a,x,xf uw,xf,x KwρaxKw†=p σafxf.\forall a_f<r_f,\ \forall x_f<s_f,\qquad \sum_{w<m}\sum_{x<s}\sum_{a<r} v_{w,a_f,a,x,x_f}\,u_{w,x_f,x}\, K_w\rho_{ax}K_w^\dagger =p\,\sigma_{a_fx_f}.∀af​<rf​, ∀xf​<sf​,w<m∑​x<s∑​a<r∑​vw,af​,a,x,xf​​uw,xf​,x​Kw​ρax​Kw†​=pσaf​xf​​.

The protocol is not required to satisfy ∑wKw†Kw=Id\sum_wK_w^\dagger K_w=I_d∑w​Kw†​Kw​=Id​. Although m=0m=0m=0 is included in the existential quantifier, its zero action cannot satisfy the displayed equality for positive ppp and normalized σ\sigmaσ.

QAlgorithms.Papers.Independent.Intakethird100.P36.SteeringMeasure

This is an abbreviation for the type of all functions SSS that take three natural numbers d,r,sd,r,sd,r,s and an arbitrary array (ρax)a<r,x<s(\rho_{ax})_{a<r,x<s}(ρax​)a<r,x<s​ of complex d×dd\times dd×d matrices and return a real number S(d,r,s,ρ)S(d,r,s,\rho)S(d,r,s,ρ). All natural-number dimensions, including zero, and all such arrays are allowed. The type imposes no positivity, normalization, invariance, convexity, or monotonicity requirement on SSS.

QAlgorithms.Papers.Independent.Intakethird100.P36.SteeringProtocol.Complete

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, let PPP consist of positive dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d satisfying Id−∑wKw†Kw⪰0I_d-\sum_wK_w^\dagger K_w\succeq0Id​−∑w​Kw†​Kw​⪰0, nonnegative input weights Iw,xf,xI_{w,x_f,x}Iw,xf​,x​ summing to one over 0≤x<s0\le x<s0≤x<s for every w,xfw,x_fw,xf​, and nonnegative output weights Ow,af,a,x,xfO_{w,a_f,a,x,x_f}Ow,af​,a,x,xf​​ summing to one over 0≤af<rf0\le a_f<r_f0≤af​<rf​ for every w,a,x,xfw,a,x,x_fw,a,x,xf​. This predicate requires the additional equality

∑w=0m−1Kw†Kw=Id.\sum_{w=0}^{m-1}K_w^\dagger K_w=I_d.w=0∑m−1​Kw†​Kw​=Id​.

It cannot hold for m=0m=0m=0, since d>0d>0d>0 and the empty sum is zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.SteeringProtocol

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, this structure consists of proofs that d,r,s,df,rf,sf>0d,r,s,d_f,r_f,s_f>0d,r,s,df​,rf​,sf​>0, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d for 0≤w<m0\le w<m0≤w<m, and real arrays Iw,xf,xI_{w,x_f,x}Iw,xf​,x​ and Ow,af,a,x,xfO_{w,a_f,a,x,x_f}Ow,af​,a,x,xf​​, where 0≤xf<sf0\le x_f<s_f0≤xf​<sf​, 0≤x<s0\le x<s0≤x<s, 0≤af<rf0\le a_f<r_f0≤af​<rf​, and 0≤a<r0\le a<r0≤a<r, satisfying

Id−∑w=0m−1Kw†Kw⪰0,Iw,xf,x≥0,∑x=0s−1Iw,xf,x=1,Ow,af,a,x,xf≥0,∑af=0rf−1Ow,af,a,x,xf=1.I_d-\sum_{w=0}^{m-1}K_w^\dagger K_w\succeq0,\qquad I_{w,x_f,x}\ge0,\quad\sum_{x=0}^{s-1}I_{w,x_f,x}=1,\qquad O_{w,a_f,a,x,x_f}\ge0,\quad \sum_{a_f=0}^{r_f-1}O_{w,a_f,a,x,x_f}=1.Id​−w=0∑m−1​Kw†​Kw​⪰0,Iw,xf​,x​≥0,x=0∑s−1​Iw,xf​,x​=1,Ow,af​,a,x,xf​​≥0,af​=0∑rf​−1​Ow,af​,a,x,xf​​=1.

The input normalization holds for every w,xfw,x_fw,xf​, and the output normalization for every w,a,x,xfw,a,x,x_fw,a,x,xf​. The parameter mmm may be zero: then these arrays have no entries, their normalization requirements are vacuous, and the matrix inequality reduces to Id⪰0I_d\succeq0Id​⪰0.

QAlgorithms.Papers.Independent.Intakethird100.P36.StrongMonotonicity

For an arbitrary real-valued function SSS on arrays of complex matrices of all natural-number dimensions, this predicate requires the following inequality for every choice of natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, every protocol PPP, and every normalized array ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​. Normalization means d,r,s>0d,r,s>0d,r,s>0, ρax⪰0\rho_{ax}\succeq0ρax​⪰0, ∑aρax\sum_a\rho_{ax}∑a​ρax​ is independent of xxx, and Mρ=Re⁡tr⁡(s−1∑x∑aρax)M_\rho=\operatorname{Re}\operatorname{tr}(s^{-1}\sum_x\sum_a\rho_{ax})Mρ​=Retr(s−1∑x​∑a​ρax​) lies in [0,1][0,1][0,1] and equals 111. The protocol has all six dimensions positive, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d, nonnegative coefficients uw,xf,xu_{w,x_f,x}uw,xf​,x​ summing to 111 over xxx, and nonnegative coefficients vw,af,a,x,xfv_{w,a_f,a,x,x_f}vw,af​,a,x,xf​​ summing to 111 over afa_faf​; it satisfies Id−∑wKw†Kw⪰0I_d-\sum_wK_w^\dagger K_w\succeq0Id​−∑w​Kw†​Kw​⪰0, and the hypothesis additionally requires completeness, ∑wKw†Kw=Id\sum_wK_w^\dagger K_w=I_d∑w​Kw†​Kw​=Id​. Define

pw=Re⁡tr⁡ ⁣(Kw(1s∑x∑aρax)Kw†),τafxf(w)=∑x∑avw,af,a,x,xfuw,xf,xKwρaxKw†,Mw=Re⁡tr⁡ ⁣(1sf∑xf∑afτafxf(w)).p_w=\operatorname{Re}\operatorname{tr}\!\left(K_w\left(\frac1s\sum_x\sum_a\rho_{ax}\right)K_w^\dagger\right),\qquad \tau^{(w)}_{a_fx_f}=\sum_x\sum_a v_{w,a_f,a,x,x_f}u_{w,x_f,x}K_w\rho_{ax}K_w^\dagger,\qquad M_w=\operatorname{Re}\operatorname{tr}\!\left(\frac1{s_f}\sum_{x_f}\sum_{a_f}\tau^{(w)}_{a_fx_f}\right).pw​=Retr(Kw​(s1​x∑​a∑​ρax​)Kw†​),τaf​xf​(w)​=x∑​a∑​vw,af​,a,x,xf​​uw,xf​,x​Kw​ρax​Kw†​,Mw​=Retr​sf​1​xf​∑​af​∑​τaf​xf​(w)​​.

The required inequality is

∑w<m{pw S(df,rf,sf,Mw−1τ(w)),pw>0,0,pw≤0 ≤ S(d,r,s,ρ),\sum_{w<m} \begin{cases} p_w\,S(d_f,r_f,s_f,M_w^{-1}\tau^{(w)}),&p_w>0,\\ 0,&p_w\le0 \end{cases} \ \le\ S(d,r,s,\rho),w<m∑​{pw​S(df​,rf​,sf​,Mw−1​τ(w)),0,​pw​>0,pw​≤0​ ≤ S(d,r,s,ρ),

where rescaling is entrywise and 0−1=00^{-1}=00−1=0. There is no requirement here for unnormalized input arrays or incomplete protocols; zero dimensions cannot satisfy the stated hypotheses, and m=0m=0m=0 cannot satisfy completeness when d>0d>0d>0.

QAlgorithms.Papers.Independent.Intakethird100.P36.WeightCandidate

For arbitrary natural numbers d,r,sd,r,sd,r,s, an arbitrary family ρ=(ρa,x)a<r,x<s\rho=(\rho_{a,x})_{a<r,x<s}ρ=(ρa,x​)a<r,x<s​ of complex d×dd\times dd×d matrices, and a real number ν\nuν, this predicate requires ν≥0\nu\ge0ν≥0 and the existence of two families η,σ\eta,\sigmaη,σ of the same shape such that

ρa,x=ν ηa,x+(1−ν)σa,xfor every a<r, x<s.\rho_{a,x}=\nu\,\eta_{a,x}+(1-\nu)\sigma_{a,x} \quad\text{for every }a<r,\ x<s.ρa,x​=νηa,x​+(1−ν)σa,x​for every a<r, x<s.

Both η\etaη and σ\sigmaσ must be normalized in the following precise sense: d,r,s>0d,r,s>0d,r,s>0, every block is positive semidefinite, the sums over aaa are identical for all inputs xxx, and the mass Re⁡Tr⁡(s−1∑x,aτa,x)\operatorname{Re}\operatorname{Tr}(s^{-1}\sum_{x,a}\tau_{a,x})ReTr(s−1∑x,a​τa,x​) equals 111 for each family τ\tauτ. In addition, σ\sigmaσ must admit some finite representation

σa,x=∑ℓ<npℓtℓ,x,aξℓ,\sigma_{a,x}=\sum_{\ell<n}p_\ell t_{\ell,x,a}\xi_\ell,σa,x​=ℓ<n∑​pℓ​tℓ,x,a​ξℓ​,

where nnn is a natural number, pℓ≥0p_\ell\ge0pℓ​≥0 and ∑ℓpℓ=1\sum_\ell p_\ell=1∑ℓ​pℓ​=1, every ξℓ\xi_\ellξℓ​ is positive semidefinite with trace 111, and tℓ,x,a≥0t_{\ell,x,a}\ge0tℓ,x,a​≥0 with ∑atℓ,x,a=1\sum_a t_{\ell,x,a}=1∑a​tℓ,x,a​=1 for every ℓ,x\ell,xℓ,x; the mass factor in the model definition is 111 because σ\sigmaσ is normalized. There is no upper bound ν≤1\nu\le1ν≤1, and no separate assemblage assumption is imposed on ρ\rhoρ. The predicate is impossible if any of d,r,sd,r,sd,r,s is zero; the probability law on the hidden indices also rules out n=0n=0n=0.

QAlgorithms.Papers.Independent.Intakethird100.P36.angleAssemblage

For every real θ\thetaθ, this definition gives a family of complex 2×22\times22×2 matrices indexed by outcome a∈{0,1}a\in\{0,1\}a∈{0,1} and setting x∈{0,1}x\in\{0,1\}x∈{0,1}:

ρ0,0=12(1000),ρ1,0=12(0001),ρ0,1=12(cos⁡2θcos⁡θsin⁡θcos⁡θsin⁡θsin⁡2θ),ρ1,1=12(sin⁡2θ−sin⁡θcos⁡θ−sin⁡θcos⁡θcos⁡2θ).\rho_{0,0}=\frac12\begin{pmatrix}1&0\\0&0\end{pmatrix}, \quad \rho_{1,0}=\frac12\begin{pmatrix}0&0\\0&1\end{pmatrix}, \quad \rho_{0,1}=\frac12\begin{pmatrix}\cos^2\theta&\cos\theta\sin\theta\\ \cos\theta\sin\theta&\sin^2\theta\end{pmatrix}, \quad \rho_{1,1}=\frac12\begin{pmatrix}\sin^2\theta&-\sin\theta\cos\theta\\ -\sin\theta\cos\theta&\cos^2\theta\end{pmatrix}.ρ0,0​=21​(10​00​),ρ1,0​=21​(00​01​),ρ0,1​=21​(cos2θcosθsinθ​cosθsinθsin2θ​),ρ1,1​=21​(sin2θ−sinθcosθ​−sinθcosθcos2θ​).

Equivalently, each block is one half of vv†vv^\daggervv†, using respectively v=(1,0),(0,1),(cos⁡θ,sin⁡θ),(−sin⁡θ,cos⁡θ)v=(1,0),(0,1),(\cos\theta,\sin\theta),(-\sin\theta,\cos\theta)v=(1,0),(0,1),(cosθ,sinθ),(−sinθ,cosθ). There is no restriction on θ\thetaθ.

QAlgorithms.Papers.Independent.Intakethird100.P36.angleKet

For every real θ\thetaθ and indices a,x∈{0,1}a,x\in\{0,1\}a,x∈{0,1}, this definition gives a vector in C2\mathbb C^2C2, with coordinates ordered as 0,10,10,1, by

vθ,0,0=(1,0),vθ,1,0=(0,1),vθ,0,1=(cos⁡θ,sin⁡θ),vθ,1,1=(−sin⁡θ,cos⁡θ).v_{\theta,0,0}=(1,0),\qquad v_{\theta,1,0}=(0,1),\qquad v_{\theta,0,1}=(\cos\theta,\sin\theta),\qquad v_{\theta,1,1}=(-\sin\theta,\cos\theta).vθ,0,0​=(1,0),vθ,1,0​=(0,1),vθ,0,1​=(cosθ,sinθ),vθ,1,1​=(−sinθ,cosθ).

The real coordinates are regarded as complex numbers, and there is no restriction on θ\thetaθ.

QAlgorithms.Papers.Independent.Intakethird100.P36.assemblageCQMatrix

For any natural numbers d,r,sd,r,sd,r,s, any array ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ of complex d×dd\times dd×d matrices, and any x<sx<sx<s, the defined matrix is indexed by pairs (a,i)(a,i)(a,i) with a<ra<ra<r and i<di<di<d, and has entries

C(a,i),(b,j)={(ρax)ij,a=b,0,a≠b.C_{(a,i),(b,j)} = \begin{cases} (\rho_{ax})_{ij},&a=b,\\ 0,&a\ne b. \end{cases}C(a,i),(b,j)​={(ρax​)ij​,0,​a=b,a=b.​

Thus its diagonal blocks are the matrices ρax\rho_{ax}ρax​ at the chosen setting xxx. No positivity, normalization, or setting-independence assumption is imposed. If r=0r=0r=0 or d=0d=0d=0, the matrix has an empty index set; if s=0s=0s=0, there is no possible argument xxx.

QAlgorithms.Papers.Independent.Intakethird100.P36.assemblageEntropy

For arbitrary natural numbers d,r,sd,r,sd,r,s and arbitrary matrix families ρ,σ=(ρa,x),(σa,x)\rho,\sigma=(\rho_{a,x}),(\sigma_{a,x})ρ,σ=(ρa,x​),(σa,x​), with a<r,x<sa<r,x<sa<r,x<s and complex d×dd\times dd×d blocks, this defines the supremum in [0,+∞][0,+\infty][0,+∞] of the following scores over all natural outcome counts hhh, complex d×dd\times dd×d matrices KgK_gKg​, and real weights ug,xu_{g,x}ug,x​ satisfying ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​, ug,x≥0u_{g,x}\ge0ug,x​≥0, and ∑xug,x=1\sum_xu_{g,x}=1∑x​ug,x​=1. Put Bτ=s−1∑x,aτa,xB_\tau=s^{-1}\sum_{x,a}\tau_{a,x}Bτ​=s−1∑x,a​τa,x​, pg=Re⁡Tr⁡(KgBρKg†)p_g=\operatorname{Re}\operatorname{Tr}(K_gB_\rho K_g^\dagger)pg​=ReTr(Kg​Bρ​Kg†​), qg=Re⁡Tr⁡(KgBσKg†)q_g=\operatorname{Re}\operatorname{Tr}(K_gB_\sigma K_g^\dagger)qg​=ReTr(Kg​Bσ​Kg†​), and Fg,xτ=⨁a<rKgτa,xKg†F^\tau_{g,x}=\bigoplus_{a<r}K_g\tau_{a,x}K_g^\daggerFg,xτ​=⨁a<r​Kg​τa,x​Kg†​. The score being maximized is

D(diag⁡p,diag⁡q)+∑g,xug,x C(pg,qg;Fg,xρ,Fg,xσ).\mathcal D(\operatorname{diag}p,\operatorname{diag}q) +\sum_{g,x}u_{g,x}\,C(p_g,q_g;F^\rho_{g,x},F^\sigma_{g,x}).D(diagp,diagq)+g,x∑​ug,x​C(pg​,qg​;Fg,xρ​,Fg,xσ​).

Here D(R,S)\mathcal D(R,S)D(R,S) is +∞+\infty+∞ unless R,SR,SR,S are positive semidefinite of trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; otherwise it is max⁡{0,Re⁡Tr⁡(f(R)−Rlog⁡S)}\max\{0,\operatorname{Re}\operatorname{Tr}(f(R)-R\log S)\}max{0,ReTr(f(R)−RlogS)}, where f(t)=tlog⁡tf(t)=t\log tf(t)=tlogt and matrix functions use continuous functional calculus. Also, C(p,q;A,B)C(p,q;A,B)C(p,q;A,B) is zero for p≤0p\le0p≤0, infinite for p>0,q≤0p>0,q\le0p>0,q≤0, and pD(A/p,B/q)p\mathcal D(A/p,B/q)pD(A/p,B/q) for p,q>0p,q>0p,q>0. Neither family is required to satisfy positivity, normalization, or input consistency. The supremum need not be attained; an empty set of scores has supremum 000. Zero dimensions and outcome counts are allowed subject to the strategy equations, and all reciprocal and extended arithmetic conventions apply.

QAlgorithms.Papers.Independent.Intakethird100.P36.assemblageFamilyMixture

For any natural numbers n,d,r,sn,d,r,sn,d,r,s, any real coefficients pjp_jpj​ indexed by j<nj<nj<n, and any family ρ(j)=(ρax(j))a<r,x<s\rho^{(j)}=(\rho^{(j)}_{ax})_{a<r,x<s}ρ(j)=(ρax(j)​)a<r,x<s​ of arrays of complex d×dd\times dd×d matrices, the defined array has entries

∑j<npj ρax(j),\sum_{j<n}p_j\,\rho^{(j)}_{ax},j<n∑​pj​ρax(j)​,

where each coefficient is regarded as a complex scalar. No nonnegativity or sum-to-one condition is imposed on the coefficients, and no positivity or normalization condition is imposed on the arrays. In particular, n=0n=0n=0 gives the zero array; the other dimensions may also be zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.assemblageMass

For arbitrary natural numbers d,r,sd,r,sd,r,s and an arbitrary family of complex d×dd\times dd×d matrices ρa,x\rho_{a,x}ρa,x​, indexed by 0≤a<r0\le a<r0≤a<r and 0≤x<s0\le x<s0≤x<s, this defines the real number

mρ=Re⁡tr⁡((s:C)−1∑x=0s−1∑a=0r−1ρa,x).m_\rho=\operatorname{Re}\operatorname{tr}\left((s:\mathbb C)^{-1} \sum_{x=0}^{s-1}\sum_{a=0}^{r-1}\rho_{a,x}\right).mρ​=Retr((s:C)−1x=0∑s−1​a=0∑r−1​ρa,x​).

No positivity or normalization is required, so this number is not required to be nonnegative. Empty sums and empty traces are zero, and 0−1=00^{-1}=00−1=0; consequently mρ=0m_\rho=0mρ​=0 if any of d,r,sd,r,sd,r,s is zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.assemblageMixture

For any natural numbers d,r,sd,r,sd,r,s, any two arrays ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ and σ=(σax)a<r,x<s\sigma=(\sigma_{ax})_{a<r,x<s}σ=(σax​)a<r,x<s​ of complex d×dd\times dd×d matrices, and any real number ttt, the defined array has entries

t ρax+(1−t) σax.t\,\rho_{ax}+(1-t)\,\sigma_{ax}.tρax​+(1−t)σax​.

Both real coefficients are regarded as complex scalars. There is no assumption that 0≤t≤10\le t\le10≤t≤1, or that either input array is positive, normalized, or has setting-independent outcome sums. Zero dimensions and empty arrays are included.

QAlgorithms.Papers.Independent.Intakethird100.P36.bobStep

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, take a protocol consisting of positive dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d with ∑wKw†Kw⪯Id\sum_wK_w^\dagger K_w\preceq I_d∑w​Kw†​Kw​⪯Id​, nonnegative input weights Iw,xf,xI_{w,x_f,x}Iw,xf​,x​ normalized over x<sx<sx<s, and nonnegative output weights Ow,af,a,x,xfO_{w,a_f,a,x,x_f}Ow,af​,a,x,xf​​ normalized over af<rfa_f<r_faf​<rf​, for every choice of their other indices. Given any family ρa,x∈Cd×d\rho_{a,x}\in\mathbb C^{d\times d}ρa,x​∈Cd×d, a<r,x<sa<r,x<sa<r,x<s, and a branch index w<mw<mw<m, this defines a family of complex df×dfd_f\times d_fdf​×df​ matrices by

(ρ(w))a,x=Kwρa,xKw†.(\rho^{(w)})_{a,x}=K_w\rho_{a,x}K_w^\dagger.(ρ(w))a,x​=Kw​ρa,x​Kw†​.

No positivity, compatibility, or normalization is assumed for ρ\rhoρ. If m=0m=0m=0, no branch index www exists.

QAlgorithms.Papers.Independent.Intakethird100.P36.bobTotal

For arbitrary natural numbers d,r,sd,r,sd,r,s and an arbitrary family of complex d×dd\times dd×d matrices ρa,x\rho_{a,x}ρa,x​, indexed by 0≤a<r0\le a<r0≤a<r and 0≤x<s0\le x<s0≤x<s, this defines the complex d×dd\times dd×d matrix

Bρ=(s:C)−1∑x=0s−1∑a=0r−1ρa,x.B_\rho=(s:\mathbb C)^{-1}\sum_{x=0}^{s-1}\sum_{a=0}^{r-1}\rho_{a,x}.Bρ​=(s:C)−1x=0∑s−1​a=0∑r−1​ρa,x​.

There are no positivity or normalization assumptions. Empty sums are zero, and the inverse of zero is defined to be zero, so Bρ=0B_\rho=0Bρ​=0 whenever s=0s=0s=0 or r=0r=0r=0.

QAlgorithms.Papers.Independent.Intakethird100.P36.conditionalD

For any finite type III equipped with decidable equality, arbitrary real numbers p,qp,qp,q, and arbitrary complex matrices A,BA,BA,B indexed by I×II\times II×I, this defines

C(p,q;A,B)={p D(A/p,B/q),p>0, q>0,+∞,p>0, q≤0,0,p≤0.C(p,q;A,B)= \begin{cases} p\,\mathcal D(A/p,B/q),&p>0,\ q>0,\\ +\infty,&p>0,\ q\le0,\\ 0,&p\le0. \end{cases}C(p,q;A,B)=⎩⎨⎧​pD(A/p,B/q),+∞,0,​p>0, q>0,p>0, q≤0,p≤0.​

The values lie in [0,+∞][0,+\infty][0,+∞], and D(R,S)\mathcal D(R,S)D(R,S) means +∞+\infty+∞ unless R,SR,SR,S are positive semidefinite, both have trace 111, and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; when all these conditions hold, it means

D(R,S)=max⁡ ⁣{0,Re⁡Tr⁡(f(R)−Rlog⁡S)},f(t)=tlog⁡t,\mathcal D(R,S)= \max\!\left\{0,\operatorname{Re}\operatorname{Tr}\bigl(f(R)-R\log S\bigr)\right\}, \qquad f(t)=t\log t,D(R,S)=max{0,ReTr(f(R)−RlogS)},f(t)=tlogt,

using continuous functional calculus and the total scalar convention log⁡0=0\log0=0log0=0. The displayed divisions mean scalar multiplication by the complex reciprocals. No sign, trace, or positivity conditions are hypotheses. In particular, p≤0p\le0p≤0 gives zero regardless of the other arguments, while for an empty III and p,q>0p,q>0p,q>0, the density trace conditions fail and the value is +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P36.conditionalKetProjector

For every natural number ddd and every complex vector v=(vi)i<dv=(v_i)_{i<d}v=(vi​)i<d​, this definition gives the complex d×dd\times dd×d matrix whose (i,j)(i,j)(i,j) entry is vivj‾v_i\overline{v_j}vi​vj​​ for all i,j<di,j<di,j<d. No normalization or nonzero condition is imposed on vvv; in particular, v=0v=0v=0 gives the zero matrix, and d=0d=0d=0 gives the empty matrix.

QAlgorithms.Papers.Independent.Intakethird100.P36.densityD

For any finite index type III with decidable equality, including the empty type, and any complex matrices ρ,σ\rho,\sigmaρ,σ indexed by I×II\times II×I, this definition returns a nonnegative extended real number as follows. If either matrix fails to be positive semidefinite with trace exactly 111, the value is +∞+\infty+∞. If both satisfy those conditions but ker⁡σ⊈ker⁡ρ\ker\sigma\not\subseteq\ker\rhokerσ⊆kerρ, where the kernels are those of their matrix actions on complex coordinate vectors, the value is also +∞+\infty+∞. Otherwise the value is

max⁡ ⁣{0, Re⁡tr⁡ ⁣(F(ρ)−ρ L(σ))},\max\!\left\{0,\, \operatorname{Re}\operatorname{tr}\!\left(F(\rho)-\rho\,L(\sigma)\right) \right\},max{0,Retr(F(ρ)−ρL(σ))},

regarded as a nonnegative extended real number, where FFF and LLL are applied by continuous functional calculus and are given by F(t)=tlog⁡tF(t)=t\log tF(t)=tlogt and L(t)=log⁡tL(t)=\log tL(t)=logt, with the natural logarithm and the total-function convention log⁡0=0\log0=0log0=0. The maximum with zero is the conversion from a real number to a nonnegative extended real number. For empty III, every matrix has trace zero, so the value is +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P36.entropyBranch

Let d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m be natural numbers, and let PPP consist of positive dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, complex df×dd_f\times ddf​×d matrices KwK_wKw​ for w<mw<mw<m satisfying Id−∑wKw†Kw≥0I_d-\sum_wK_w^\dagger K_w\ge0Id​−∑w​Kw†​Kw​≥0, and real coefficients uw,y,x≥0u_{w,y,x}\ge0uw,y,x​≥0, vw,b,a,x,y≥0v_{w,b,a,x,y}\ge0vw,b,a,x,y​≥0 satisfying ∑x<suw,y,x=1\sum_{x<s}u_{w,y,x}=1∑x<s​uw,y,x​=1 and ∑b<rfvw,b,a,x,y=1\sum_{b<r_f}v_{w,b,a,x,y}=1∑b<rf​​vw,b,a,x,y​=1. For arbitrary input matrix families ρ,σ\rho,\sigmaρ,σ indexed by a<r,x<sa<r,x<sa<r,x<s, and an outcome w<mw<mw<m, put

Bτ=1s∑x,aτa,x,p=Re⁡Tr⁡(KwBρKw†),q=Re⁡Tr⁡(KwBσKw†),Γb,yτ=∑x,avw,b,a,x,yuw,y,xKwτa,xKw†,B_\tau=\frac1s\sum_{x,a}\tau_{a,x},\quad p=\operatorname{Re}\operatorname{Tr}(K_wB_\rho K_w^\dagger),\quad q=\operatorname{Re}\operatorname{Tr}(K_wB_\sigma K_w^\dagger),\quad \Gamma^\tau_{b,y}=\sum_{x,a}v_{w,b,a,x,y}u_{w,y,x}K_w\tau_{a,x}K_w^\dagger,Bτ​=s1​x,a∑​τa,x​,p=ReTr(Kw​Bρ​Kw†​),q=ReTr(Kw​Bσ​Kw†​),Γb,yτ​=x,a∑​vw,b,a,x,y​uw,y,x​Kw​τa,x​Kw†​,

and Mτ=Re⁡Tr⁡(sf−1∑y,bΓb,yτ)M_\tau=\operatorname{Re}\operatorname{Tr}(s_f^{-1}\sum_{y,b}\Gamma^\tau_{b,y})Mτ​=ReTr(sf−1​∑y,b​Γb,yτ​). This defines the value to be 000 when p≤0p\le0p≤0, +∞+\infty+∞ when p>0,q≤0p>0,q\le0p>0,q≤0, and p H(Γρ/Mρ,Γσ/Mσ)p\,H(\Gamma^\rho/M_\rho,\Gamma^\sigma/M_\sigma)pH(Γρ/Mρ​,Γσ/Mσ​) when p,q>0p,q>0p,q>0, with division applied to every block. For families α,β\alpha,\betaα,β of df×dfd_f\times d_fdf​×df​ blocks indexed by b<rf,y<sfb<r_f,y<s_fb<rf​,y<sf​, H(α,β)H(\alpha,\beta)H(α,β) is the supremum over finite matrices LgL_gLg​ satisfying ∑gLg†Lg=Idf\sum_gL_g^\dagger L_g=I_{d_f}∑g​Lg†​Lg​=Idf​​ and weights zg,y≥0z_{g,y}\ge0zg,y​≥0, ∑yzg,y=1\sum_yz_{g,y}=1∑y​zg,y​=1, of D(diag⁡a,diag⁡b)+∑g,yzg,yC(ag,bg;Fg,yα,Fg,yβ)\mathcal D(\operatorname{diag}a,\operatorname{diag}b)+\sum_{g,y}z_{g,y}C(a_g,b_g;F^\alpha_{g,y},F^\beta_{g,y})D(diaga,diagb)+∑g,y​zg,y​C(ag​,bg​;Fg,yα​,Fg,yβ​), where ag=Re⁡Tr⁡(Lg(sf−1∑y,bαb,y)Lg†)a_g=\operatorname{Re}\operatorname{Tr}(L_g(s_f^{-1}\sum_{y,b}\alpha_{b,y})L_g^\dagger)ag​=ReTr(Lg​(sf−1​∑y,b​αb,y​)Lg†​), bgb_gbg​ is the same expression for β\betaβ, and Fg,yτ=⨁bLgτb,yLg†F^\tau_{g,y}=\bigoplus_bL_g\tau_{b,y}L_g^\daggerFg,yτ​=⨁b​Lg​τb,y​Lg†​. Here D(R,S)\mathcal D(R,S)D(R,S) is infinite unless both matrices are positive semidefinite with trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; otherwise it is max⁡{0,Re⁡Tr⁡(f(R)−Rlog⁡S)}\max\{0,\operatorname{Re}\operatorname{Tr}(f(R)-R\log S)\}max{0,ReTr(f(R)−RlogS)}, with f(t)=tlog⁡tf(t)=t\log tf(t)=tlogt and continuous functional calculus. The term C(a,b;A,B)C(a,b;A,B)C(a,b;A,B) is zero for a≤0a\le0a≤0, infinite for a>0,b≤0a>0,b\le0a>0,b≤0, and aD(A/a,B/b)a\mathcal D(A/a,B/b)aD(A/a,B/b) for a,b>0a,b>0a,b>0. Neither ρ\rhoρ nor σ\sigmaσ is assumed normalized or positive, the protocol need not satisfy ∑wKw†Kw=Id\sum_wK_w^\dagger K_w=I_d∑w​Kw†​Kw​=Id​, and normalization uses each branch’s own mass MτM_\tauMτ​, with 0−1=00^{-1}=00−1=0. The index www requires m>0m>0m>0; extended arithmetic includes 0⋅(+∞)=00\cdot(+\infty)=00⋅(+∞)=0.

QAlgorithms.Papers.Independent.Intakethird100.P36.entropyScore

For arbitrary natural numbers d,r,sd,r,sd,r,s, arbitrary families ρ,σ\rho,\sigmaρ,σ of complex d×dd\times dd×d matrices indexed by a<r,x<sa<r,x<sa<r,x<s, and a strategy TTT consisting of hhh outcomes, matrices KgK_gKg​ with ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​, and weights ug,x≥0u_{g,x}\ge0ug,x​≥0 with ∑xug,x=1\sum_xu_{g,x}=1∑x​ug,x​=1, define Bτ=s−1∑x,aτa,xB_\tau=s^{-1}\sum_{x,a}\tau_{a,x}Bτ​=s−1∑x,a​τa,x​, pg=Re⁡Tr⁡(KgBρKg†)p_g=\operatorname{Re}\operatorname{Tr}(K_gB_\rho K_g^\dagger)pg​=ReTr(Kg​Bρ​Kg†​), and qg=Re⁡Tr⁡(KgBσKg†)q_g=\operatorname{Re}\operatorname{Tr}(K_gB_\sigma K_g^\dagger)qg​=ReTr(Kg​Bσ​Kg†​). Let Fg,xτF^\tau_{g,x}Fg,xτ​ be the block diagonal matrix with blocks Kgτa,xKg†K_g\tau_{a,x}K_g^\daggerKg​τa,x​Kg†​ indexed by a<ra<ra<r. The score is

D(diag⁡p,diag⁡q)+∑g<h∑x<smax⁡{ug,x,0} C(pg,qg;Fg,xρ,Fg,xσ).\mathcal D(\operatorname{diag}p,\operatorname{diag}q) +\sum_{g<h}\sum_{x<s}\max\{u_{g,x},0\}\, C(p_g,q_g;F^\rho_{g,x},F^\sigma_{g,x}).D(diagp,diagq)+g<h∑​x<s∑​max{ug,x​,0}C(pg​,qg​;Fg,xρ​,Fg,xσ​).

Here D(R,S)=+∞\mathcal D(R,S)=+\inftyD(R,S)=+∞ unless R,SR,SR,S are positive semidefinite with trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; otherwise it is max⁡{0,Re⁡Tr⁡(f(R)−Rlog⁡S)}\max\{0,\operatorname{Re}\operatorname{Tr}(f(R)-R\log S)\}max{0,ReTr(f(R)−RlogS)}, with f(t)=tlog⁡tf(t)=t\log tf(t)=tlogt and continuous functional calculus. The quantity C(p,q;A,B)C(p,q;A,B)C(p,q;A,B) is 000 if p≤0p\le0p≤0, +∞+\infty+∞ if p>0,q≤0p>0,q\le0p>0,q≤0, and pD(A/p,B/q)p\mathcal D(A/p,B/q)pD(A/p,B/q) if p,q>0p,q>0p,q>0. All arithmetic in the score is extended nonnegative arithmetic, including 0⋅(+∞)=00\cdot(+\infty)=00⋅(+∞)=0. No assemblage or normalization assumptions are imposed on ρ,σ\rho,\sigmaρ,σ; zero dimensions, empty sums, and the total reciprocal 0−1=00^{-1}=00−1=0 are included.

QAlgorithms.Papers.Independent.Intakethird100.P36.filteredCQ

For natural numbers d,r,sd,r,sd,r,s, let TTT have hhh outcomes, complex d×dd\times dd×d matrices KgK_gKg​ with ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​, and real input weights ug,x≥0u_{g,x}\ge0ug,x​≥0 with ∑x<sug,x=1\sum_{x<s}u_{g,x}=1∑x<s​ug,x​=1 for each outcome. Given an arbitrary family ρ=(ρa,y)a<r,y<s\rho=(\rho_{a,y})_{a<r,y<s}ρ=(ρa,y​)a<r,y<s​ of complex d×dd\times dd×d matrices, an outcome g<hg<hg<h, and an input x<sx<sx<s, this defines a complex matrix indexed by pairs (a,i)∈{0,…,r−1}×{0,…,d−1}(a,i)\in\{0,\ldots,r-1\}\times\{0,\ldots,d-1\}(a,i)∈{0,…,r−1}×{0,…,d−1}, whose entries are

F(a,i),(b,j)={(Kgρa,xKg†)i,j,a=b,0,a≠b.F_{(a,i),(b,j)} = \begin{cases} (K_g\rho_{a,x}K_g^\dagger)_{i,j},&a=b,\\ 0,&a\ne b. \end{cases}F(a,i),(b,j)​={(Kg​ρa,x​Kg†​)i,j​,0,​a=b,a=b.​

The input weights do not enter this formula. No positivity, normalization, or consistency condition is imposed on ρ\rhoρ. The binders ggg and xxx require existing outcome and input indices; r=0r=0r=0 or d=0d=0d=0 still permits an empty resulting matrix.

QAlgorithms.Papers.Independent.Intakethird100.P36.measureBranch

Let SSS be any real-valued function on arrays of complex matrices of every natural-number dimension. For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, let PPP be a protocol with all six dimensions positive, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d satisfying Id−∑w<mKw†Kw⪰0I_d-\sum_{w<m}K_w^\dagger K_w\succeq0Id​−∑w<m​Kw†​Kw​⪰0, and real arrays uw,xf,x,vw,af,a,x,xfu_{w,x_f,x},v_{w,a_f,a,x,x_f}uw,xf​,x​,vw,af​,a,x,xf​​ that are nonnegative and sum to 111 over x<sx<sx<s and af<rfa_f<r_faf​<rf​, respectively, for every choice of their other indices. Given any array ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ of complex d×dd\times dd×d matrices and any w<mw<mw<m, put

bρ=1s∑x∑aρax,pw=Re⁡tr⁡(KwbρKw†),τafxf(w)=∑x∑avw,af,a,x,xf uw,xf,x KwρaxKw†,Mw=Re⁡tr⁡ ⁣(1sf∑xf∑afτafxf(w)).b_\rho=\frac1s\sum_x\sum_a\rho_{ax},\qquad p_w=\operatorname{Re}\operatorname{tr}(K_wb_\rho K_w^\dagger),\qquad \tau^{(w)}_{a_fx_f}=\sum_x\sum_a v_{w,a_f,a,x,x_f}\,u_{w,x_f,x}\,K_w\rho_{ax}K_w^\dagger, \qquad M_w=\operatorname{Re}\operatorname{tr}\!\left(\frac1{s_f}\sum_{x_f}\sum_{a_f}\tau^{(w)}_{a_fx_f}\right).bρ​=s1​x∑​a∑​ρax​,pw​=Retr(Kw​bρ​Kw†​),τaf​xf​(w)​=x∑​a∑​vw,af​,a,x,xf​​uw,xf​,x​Kw​ρax​Kw†​,Mw​=Retr​sf​1​xf​∑​af​∑​τaf​xf​(w)​​.

The defined value is pw S(df,rf,sf,(Mw−1τafxf(w))af,xf)p_w\,S(d_f,r_f,s_f,(M_w^{-1}\tau^{(w)}_{a_fx_f})_{a_f,x_f})pw​S(df​,rf​,sf​,(Mw−1​τaf​xf​(w)​)af​,xf​​) if pw>0p_w>0pw​>0, and is 000 otherwise. No validity or normalization condition is imposed on ρ\rhoρ, and SSS is unrestricted. In particular, zero and negative pwp_wpw​ both select the zero value, and any zero mass in the rescaling uses 0−1=00^{-1}=00−1=0.

QAlgorithms.Papers.Independent.Intakethird100.P36.normalizeAssemblage

For any natural numbers d,r,sd,r,sd,r,s and any array ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ of complex d×dd\times dd×d matrices, define

bρ=1s∑x<s∑a<rρax,Mρ=Re⁡tr⁡(bρ).b_\rho=\frac1s\sum_{x<s}\sum_{a<r}\rho_{ax}, \qquad M_\rho=\operatorname{Re}\operatorname{tr}(b_\rho).bρ​=s1​x<s∑​a<r∑​ρax​,Mρ​=Retr(bρ​).

The resulting array has entries (Mρ−1)ρax(M_\rho^{-1})\rho_{ax}(Mρ−1​)ρax​, with the real scalar regarded as complex. No positivity, normalization, or nonzero-mass assumption is imposed. In particular, inversion at zero is defined to be zero, so Mρ=0M_\rho=0Mρ​=0 produces the zero array. Zero dimensions and empty index sets are included; when s=0s=0s=0, the scalar 1/s1/s1/s is also zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.oneWayFormula

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, take a protocol with positive dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d satisfying ∑wKw†Kw⪯Id\sum_wK_w^\dagger K_w\preceq I_d∑w​Kw†​Kw​⪯Id​, nonnegative input weights Iw,xf,xI_{w,x_f,x}Iw,xf​,x​ summing to one over x<sx<sx<s, and nonnegative output weights Ow,af,a,x,xfO_{w,a_f,a,x,x_f}Ow,af​,a,x,xf​​ summing to one over af<rfa_f<r_faf​<rf​, for every choice of their other indices. For any family ρa,x∈Cd×d\rho_{a,x}\in\mathbb C^{d\times d}ρa,x​∈Cd×d, a<r,x<sa<r,x<sa<r,x<s, this defines a family of complex df×dfd_f\times d_fdf​×df​ matrices, indexed by af<rf,xf<sfa_f<r_f,x_f<s_faf​<rf​,xf​<sf​, through

F(ρ)af,xf=∑w=0m−1∑a=0r−1∑x=0s−1(Iw,xf,xOw,af,a,x,xf)Kwρa,xKw†,F(\rho)_{a_f,x_f} =\sum_{w=0}^{m-1}\sum_{a=0}^{r-1}\sum_{x=0}^{s-1} \left(I_{w,x_f,x}O_{w,a_f,a,x,x_f}\right) K_w\rho_{a,x}K_w^\dagger,F(ρ)af​,xf​​=w=0∑m−1​a=0∑r−1​x=0∑s−1​(Iw,xf​,x​Ow,af​,a,x,xf​​)Kw​ρa,x​Kw†​,

where the real product in parentheses is interpreted as a complex scalar. The input family is arbitrary, and protocol completeness is not assumed. If m=0m=0m=0, the sum gives the zero family.

QAlgorithms.Papers.Independent.Intakethird100.P36.protocolAction

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, take a protocol with positive dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d satisfying ∑wKw†Kw⪯Id\sum_wK_w^\dagger K_w\preceq I_d∑w​Kw†​Kw​⪯Id​, nonnegative input weights Iw,xf,xI_{w,x_f,x}Iw,xf​,x​ normalized over x<sx<sx<s, and nonnegative output weights Ow,af,a,x,xfO_{w,a_f,a,x,x_f}Ow,af​,a,x,xf​​ normalized over af<rfa_f<r_faf​<rf​, for every choice of their other indices. For any family ρa,x∈Cd×d\rho_{a,x}\in\mathbb C^{d\times d}ρa,x​∈Cd×d, a<r,x<sa<r,x<sa<r,x<s, this defines the complex df×dfd_f\times d_fdf​×df​ blocks

P(ρ)af,xf=∑w=0m−1(∑x=0s−1∑a=0r−1Ow,af,a,x,xf(Iw,xf,x(Kwρa,xKw†))),af<rf, xf<sf.\mathcal P(\rho)_{a_f,x_f} =\sum_{w=0}^{m-1}\left( \sum_{x=0}^{s-1}\sum_{a=0}^{r-1} O_{w,a_f,a,x,x_f} \left(I_{w,x_f,x}\left(K_w\rho_{a,x}K_w^\dagger\right)\right)\right), \quad a_f<r_f,\ x_f<s_f.P(ρ)af​,xf​​=w=0∑m−1​(x=0∑s−1​a=0∑r−1​Ow,af​,a,x,xf​​(Iw,xf​,x​(Kw​ρa,x​Kw†​))),af​<rf​, xf​<sf​.

No assemblage or normalization condition is imposed on ρ\rhoρ, and no equality ∑wKw†Kw=Id\sum_wK_w^\dagger K_w=I_d∑w​Kw†​Kw​=Id​ is required. For m=0m=0m=0, every resulting block is zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.protocolBranch

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, take a protocol with positive dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​, matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d satisfying ∑wKw†Kw⪯Id\sum_wK_w^\dagger K_w\preceq I_d∑w​Kw†​Kw​⪯Id​, nonnegative weights Iw,xf,xI_{w,x_f,x}Iw,xf​,x​ summing to one over x<sx<sx<s, and nonnegative weights Ow,af,a,x,xfO_{w,a_f,a,x,x_f}Ow,af​,a,x,xf​​ summing to one over af<rfa_f<r_faf​<rf​, for every choice of the remaining indices. For any family ρa,x∈Cd×d\rho_{a,x}\in\mathbb C^{d\times d}ρa,x​∈Cd×d, a<r,x<sa<r,x<sa<r,x<s, and any w<mw<mw<m, this defines the family indexed by af<rf,xf<sfa_f<r_f,x_f<s_faf​<rf​,xf​<sf​ whose blocks are

Baf,xf(w)=∑x=0s−1∑a=0r−1Ow,af,a,x,xf(Iw,xf,x(Kwρa,xKw†)),B^{(w)}_{a_f,x_f} =\sum_{x=0}^{s-1}\sum_{a=0}^{r-1} O_{w,a_f,a,x,x_f}\left(I_{w,x_f,x} \left(K_w\rho_{a,x}K_w^\dagger\right)\right),Baf​,xf​(w)​=x=0∑s−1​a=0∑r−1​Ow,af​,a,x,xf​​(Iw,xf​,x​(Kw​ρa,x​Kw†​)),

with real weights acting as complex scalars. The family ρ\rhoρ is arbitrary and need not satisfy assemblage or normalization conditions; when m=0m=0m=0, there is no admissible www.

QAlgorithms.Papers.Independent.Intakethird100.P36.relativeEntropyMeasure

This defines a real-valued function of every triple of natural numbers d,r,sd,r,sd,r,s and every family ρ=(ρa,x)a<r,x<s\rho=(\rho_{a,x})_{a<r,x<s}ρ=(ρa,x​)a<r,x<s​ of complex d×dd\times dd×d matrices by taking the real conversion of the following extended nonnegative infimum. Minimize H(ρ,σ)H(\rho,\sigma)H(ρ,σ) over families σ\sigmaσ with d,r,s>0d,r,s>0d,r,s>0, positive semidefinite blocks, input-independent sums ∑aσa,x\sum_a\sigma_{a,x}∑a​σa,x​, mass Re⁡Tr⁡(s−1∑x,aσa,x)=1\operatorname{Re}\operatorname{Tr}(s^{-1}\sum_{x,a}\sigma_{a,x})=1ReTr(s−1∑x,a​σa,x​)=1, and a finite representation σa,x=∑ℓ<nvℓtℓ,x,aξℓ\sigma_{a,x}=\sum_{\ell<n}v_\ell t_{\ell,x,a}\xi_\ellσa,x​=∑ℓ<n​vℓ​tℓ,x,a​ξℓ​, where vℓ≥0v_\ell\ge0vℓ​≥0, ∑ℓvℓ=1\sum_\ell v_\ell=1∑ℓ​vℓ​=1, each ξℓ\xi_\ellξℓ​ is positive semidefinite with trace 111, and tℓ,x,a≥0t_{\ell,x,a}\ge0tℓ,x,a​≥0, ∑atℓ,x,a=1\sum_a t_{\ell,x,a}=1∑a​tℓ,x,a​=1. Here H(ρ,σ)H(\rho,\sigma)H(ρ,σ) is the supremum, over finite families KgK_gKg​ with ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​ and weights ug,x≥0u_{g,x}\ge0ug,x​≥0, ∑xug,x=1\sum_xu_{g,x}=1∑x​ug,x​=1, of

D(diag⁡p,diag⁡q)+∑g,xug,x C(pg,qg;Fg,xρ,Fg,xσ),\mathcal D(\operatorname{diag}p,\operatorname{diag}q) +\sum_{g,x}u_{g,x}\,C(p_g,q_g;F^\rho_{g,x},F^\sigma_{g,x}),D(diagp,diagq)+g,x∑​ug,x​C(pg​,qg​;Fg,xρ​,Fg,xσ​),

where Bτ=s−1∑x,aτa,xB_\tau=s^{-1}\sum_{x,a}\tau_{a,x}Bτ​=s−1∑x,a​τa,x​, pg=Re⁡Tr⁡(KgBρKg†)p_g=\operatorname{Re}\operatorname{Tr}(K_gB_\rho K_g^\dagger)pg​=ReTr(Kg​Bρ​Kg†​), qg=Re⁡Tr⁡(KgBσKg†)q_g=\operatorname{Re}\operatorname{Tr}(K_gB_\sigma K_g^\dagger)qg​=ReTr(Kg​Bσ​Kg†​), and Fg,xτ=⨁aKgτa,xKg†F^\tau_{g,x}=\bigoplus_aK_g\tau_{a,x}K_g^\daggerFg,xτ​=⨁a​Kg​τa,x​Kg†​. The divergence D(R,S)\mathcal D(R,S)D(R,S) is +∞+\infty+∞ unless both matrices are positive semidefinite of trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; otherwise it is max⁡{0,Re⁡Tr⁡(f(R)−Rlog⁡S)}\max\{0,\operatorname{Re}\operatorname{Tr}(f(R)-R\log S)\}max{0,ReTr(f(R)−RlogS)}, with f(t)=tlog⁡tf(t)=t\log tf(t)=tlogt and continuous functional calculus. The conditional term CCC is zero when p≤0p\le0p≤0, infinite when p>0,q≤0p>0,q\le0p>0,q≤0, and pD(A/p,B/q)p\mathcal D(A/p,B/q)pD(A/p,B/q) when p,q>0p,q>0p,q>0. Real conversion preserves finite values and sends +∞+\infty+∞ to 000. Thus an empty minimizing set, including any case with d=0d=0d=0, r=0r=0r=0, or s=0s=0s=0, yields the real value 000. No assemblage assumption on ρ\rhoρ, monotonicity property, or attainment condition is part of this definition.

QAlgorithms.Papers.Independent.Intakethird100.P36.steerableWeight

For arbitrary natural numbers d,r,sd,r,sd,r,s and a family ρ=(ρa,x)a<r, x<s\rho=(\rho_{a,x})_{a<r,\ x<s}ρ=(ρa,x​)a<r, x<s​ of complex d×dd\times dd×d matrices, the steerable weight is the real infimum of all ν≥0\nu\geq0ν≥0 for which there exist normalized families η,σ\eta,\sigmaη,σ of the same dimensions satisfying ρa,x=νηa,x+(1−ν)σa,x\rho_{a,x}=\nu\eta_{a,x}+(1-\nu)\sigma_{a,x}ρa,x​=νηa,x​+(1−ν)σa,x​ for every a,xa,xa,x, with σ\sigmaσ admitting a finite hidden-state representation. Here a normalized family has d,r,s>0d,r,s>0d,r,s>0, positive semidefinite blocks, sums ∑aρa,x\sum_a\rho_{a,x}∑a​ρa,x​ independent of xxx, and Re⁡tr⁡B(ρ)=1\operatorname{Re}\operatorname{tr}B(\rho)=1RetrB(ρ)=1, where B(ρ)=s−1∑x∑aρa,xB(\rho)=s^{-1}\sum_x\sum_a\rho_{a,x}B(ρ)=s−1∑x​∑a​ρa,x​. The required representation of the normalized σ\sigmaσ consists of a natural number hhh, probabilities pℓ≥0p_\ell\geq0pℓ​≥0 summing to 111, positive semidefinite matrices ξℓ\xi_\ellξℓ​ of trace 111, and response probabilities qℓ,x(a)≥0q_{\ell,x}(a)\geq0qℓ,x​(a)≥0 summing to 111 over aaa, such that σa,x=∑ℓ<hpℓqℓ,x(a)ξℓ\sigma_{a,x}=\sum_{\ell<h}p_\ell q_{\ell,x}(a)\xi_\ellσa,x​=∑ℓ<h​pℓ​qℓ,x​(a)ξℓ​. There is no explicit restriction ν≤1\nu\leq1ν≤1 or assumption that ρ\rhoρ is an assemblage. If the candidate set is empty, its real infimum is defined to be 000; this occurs in particular if any of d,r,sd,r,sd,r,s is zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.steeringBranchProbability

For natural numbers d,r,s,df,rf,sf,md,r,s,d_f,r_f,s_f,md,r,s,df​,rf​,sf​,m, let ρ=(ρax)a<r,x<s\rho=(\rho_{ax})_{a<r,x<s}ρ=(ρax​)a<r,x<s​ be any array of complex d×dd\times dd×d matrices, and let PPP be a protocol with all six dimensions d,r,s,df,rf,sfd,r,s,d_f,r_f,s_fd,r,s,df​,rf​,sf​ positive. Its data are matrices Kw∈Cdf×dK_w\in\mathbb C^{d_f\times d}Kw​∈Cdf​×d, indexed by w<mw<mw<m, satisfying Id−∑wKw†Kw⪰0I_d-\sum_wK_w^\dagger K_w\succeq0Id​−∑w​Kw†​Kw​⪰0, and real arrays uw,xf,xu_{w,x_f,x}uw,xf​,x​ and vw,af,a,x,xfv_{w,a_f,a,x,x_f}vw,af​,a,x,xf​​, where for each w,xfw,x_fw,xf​, the numbers uw,xf,xu_{w,x_f,x}uw,xf​,x​ are nonnegative and sum to 111 over x<sx<sx<s, and for each w,a,x,xfw,a,x,x_fw,a,x,xf​, the numbers vw,af,a,x,xfv_{w,a_f,a,x,x_f}vw,af​,a,x,xf​​ are nonnegative and sum to 111 over af<rfa_f<r_faf​<rf​. For any w<mw<mw<m, the defined real number is

Re⁡tr⁡ ⁣(Kw(1s∑x<s∑a<rρax)Kw†).\operatorname{Re}\operatorname{tr}\!\left( K_w\left(\frac1s\sum_{x<s}\sum_{a<r}\rho_{ax}\right)K_w^\dagger \right).Retr(Kw​(s1​x<s∑​a<r∑​ρax​)Kw†​).

No positivity, normalization, or equality of the sums ∑aρax\sum_a\rho_{ax}∑a​ρax​ is required of ρ\rhoρ; consequently, the definition does not itself require this real number to lie in [0,1][0,1][0,1]. The protocol’s arrays u,vu,vu,v do not enter this expression.

QAlgorithms.Papers.Independent.Intakethird100.P36.steeringClassicalDivergence

For every natural number nnn, including 000, and arbitrary real functions p,q:{0,…,n−1}→Rp,q:\{0,\ldots,n-1\}\to\mathbb Rp,q:{0,…,n−1}→R, this defines an extended nonnegative real number as follows. Put P=diag⁡(pi)P=\operatorname{diag}(p_i)P=diag(pi​) and Q=diag⁡(qi)Q=\operatorname{diag}(q_i)Q=diag(qi​), regarded as complex matrices. If both matrices are positive semidefinite with trace 111, and ker⁡Q⊆ker⁡P\ker Q\subseteq\ker PkerQ⊆kerP, the value is

max⁡ ⁣{0, Re⁡Tr⁡(f(P)−Plog⁡Q)},f(t)=tlog⁡t,\max\!\left\{0,\ \operatorname{Re}\operatorname{Tr}\bigl(f(P)-P\log Q\bigr)\right\}, \qquad f(t)=t\log t,max{0, ReTr(f(P)−PlogQ)},f(t)=tlogt,

where functions of matrices use continuous functional calculus and the scalar logarithm has its total convention log⁡0=0\log 0=0log0=0. Otherwise the value is +∞+\infty+∞. Thus the density conditions require pi,qi≥0p_i,q_i\ge0pi​,qi​≥0 and ∑ipi=∑iqi=1\sum_i p_i=\sum_i q_i=1∑i​pi​=∑i​qi​=1, and the kernel condition requires qi=0⇒pi=0q_i=0\Rightarrow p_i=0qi​=0⇒pi​=0. No such conditions are assumed on the arguments beforehand. When n=0n=0n=0, the trace condition fails and the value is +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P36.steeringEntropy

For arbitrary natural numbers d,r,sd,r,sd,r,s and any family ρ=(ρa,x)a<r,x<s\rho=(\rho_{a,x})_{a<r,x<s}ρ=(ρa,x​)a<r,x<s​ of complex d×dd\times dd×d matrices, this defines the infimum in [0,+∞][0,+\infty][0,+∞], over families σ\sigmaσ described below, of a quantity H(ρ,σ)H(\rho,\sigma)H(ρ,σ). An admissible σ\sigmaσ requires d,r,s>0d,r,s>0d,r,s>0, positive semidefinite blocks, equal sums ∑aσa,x\sum_a\sigma_{a,x}∑a​σa,x​ for all inputs xxx, and mass Re⁡Tr⁡(s−1∑x,aσa,x)=1\operatorname{Re}\operatorname{Tr}(s^{-1}\sum_{x,a}\sigma_{a,x})=1ReTr(s−1∑x,a​σa,x​)=1; it must also admit a finite model: for some natural nnn, numbers vℓ≥0v_\ell\ge0vℓ​≥0 summing to 111, positive semidefinite matrices ξℓ\xi_\ellξℓ​ of trace 111, and response numbers tℓ,x,a≥0t_{\ell,x,a}\ge0tℓ,x,a​≥0 summing to 111 over aaa, one has σa,x=∑ℓ<nvℓtℓ,x,aξℓ\sigma_{a,x}=\sum_{\ell<n}v_\ell t_{\ell,x,a}\xi_\ellσa,x​=∑ℓ<n​vℓ​tℓ,x,a​ξℓ​. To define H(ρ,σ)H(\rho,\sigma)H(ρ,σ), take the supremum over natural hhh, matrices KgK_gKg​ with ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​, and weights ug,x≥0u_{g,x}\ge0ug,x​≥0 with ∑xug,x=1\sum_xu_{g,x}=1∑x​ug,x​=1, of

D(diag⁡p,diag⁡q)+∑g,xug,x C(pg,qg;Fg,xρ,Fg,xσ),\mathcal D(\operatorname{diag}p,\operatorname{diag}q) +\sum_{g,x}u_{g,x}\,C(p_g,q_g;F^\rho_{g,x},F^\sigma_{g,x}),D(diagp,diagq)+g,x∑​ug,x​C(pg​,qg​;Fg,xρ​,Fg,xσ​),

where Bτ=s−1∑x,aτa,xB_\tau=s^{-1}\sum_{x,a}\tau_{a,x}Bτ​=s−1∑x,a​τa,x​, pg=Re⁡Tr⁡(KgBρKg†)p_g=\operatorname{Re}\operatorname{Tr}(K_gB_\rho K_g^\dagger)pg​=ReTr(Kg​Bρ​Kg†​), qg=Re⁡Tr⁡(KgBσKg†)q_g=\operatorname{Re}\operatorname{Tr}(K_gB_\sigma K_g^\dagger)qg​=ReTr(Kg​Bσ​Kg†​), and Fg,xτ=⨁aKgτa,xKg†F^\tau_{g,x}=\bigoplus_aK_g\tau_{a,x}K_g^\daggerFg,xτ​=⨁a​Kg​τa,x​Kg†​. Here D(R,S)\mathcal D(R,S)D(R,S) is infinite unless R,SR,SR,S are positive semidefinite of trace 111 and ker⁡S⊆ker⁡R\ker S\subseteq\ker RkerS⊆kerR; otherwise it is max⁡{0,Re⁡Tr⁡(f(R)−Rlog⁡S)}\max\{0,\operatorname{Re}\operatorname{Tr}(f(R)-R\log S)\}max{0,ReTr(f(R)−RlogS)}, with f(t)=tlog⁡tf(t)=t\log tf(t)=tlogt and continuous functional calculus. The term C(p,q;A,B)C(p,q;A,B)C(p,q;A,B) is zero for p≤0p\le0p≤0, infinite for p>0,q≤0p>0,q\le0p>0,q≤0, and pD(A/p,B/q)p\mathcal D(A/p,B/q)pD(A/p,B/q) for p,q>0p,q>0p,q>0. There is no assumption on ρ\rhoρ, and no attainment is required for either extremum. An empty admissible set has infimum +∞+\infty+∞, in particular whenever any of d,r,sd,r,sd,r,s is zero; an empty score set has supremum 000, and extended arithmetic uses 0⋅(+∞)=00\cdot(+\infty)=00⋅(+∞)=0.

QAlgorithms.Papers.Independent.Intakethird100.P36.steeringRobustness

For arbitrary natural numbers d,r,sd,r,sd,r,s and a family ρ=(ρa,x)a<r, x<s\rho=(\rho_{a,x})_{a<r,\ x<s}ρ=(ρa,x​)a<r, x<s​ of complex d×dd\times dd×d matrices, the steering robustness is the real infimum of the numbers ν≥0\nu\geq0ν≥0 for which there exists a normalized family η\etaη such that Ya,x=(ρa,x+νηa,x)/(1+ν)Y_{a,x}=(\rho_{a,x}+\nu\eta_{a,x})/(1+\nu)Ya,x​=(ρa,x​+νηa,x​)/(1+ν) has a finite hidden-state model. For a family FFF, write μ(F)=Re⁡tr⁡(s−1∑x∑aFa,x)\mu(F)=\operatorname{Re}\operatorname{tr}(s^{-1}\sum_x\sum_aF_{a,x})μ(F)=Retr(s−1∑x​∑a​Fa,x​). Normalization of η\etaη means that d,r,s>0d,r,s>0d,r,s>0, its blocks are positive semidefinite, its outcome sums are independent of the setting, and μ(η)=1\mu(\eta)=1μ(η)=1. The hidden-state condition on YYY means that its blocks are positive semidefinite, its outcome sums are setting-independent, 0≤μ(Y)≤10\leq\mu(Y)\leq10≤μ(Y)≤1, and there are finitely many density matrices ξℓ\xi_\ellξℓ​, probabilities pℓp_\ellpℓ​ summing to 111, and response probabilities qℓ,x(a)q_{\ell,x}(a)qℓ,x​(a) summing to 111 over aaa, with Ya,x=μ(Y)∑ℓpℓqℓ,x(a)ξℓY_{a,x}=\mu(Y)\sum_\ell p_\ell q_{\ell,x}(a)\xi_\ellYa,x​=μ(Y)∑ℓ​pℓ​qℓ,x​(a)ξℓ​. Each density matrix is positive semidefinite with trace 111. The definition assumes no such conditions on ρ\rhoρ itself and places no upper bound on ν\nuν. An empty candidate set has real infimum 000, including when any dimension d,r,sd,r,sd,r,s is zero.

QAlgorithms.Papers.Independent.Intakethird100.P36.strategyProbability

For natural numbers d,r,sd,r,sd,r,s, let TTT consist of a natural outcome count hhh, complex d×dd\times dd×d matrices KgK_gKg​ satisfying ∑gKg†Kg=Id\sum_gK_g^\dagger K_g=I_d∑g​Kg†​Kg​=Id​, and real input weights ug,x≥0u_{g,x}\ge0ug,x​≥0 satisfying ∑x<sug,x=1\sum_{x<s}u_{g,x}=1∑x<s​ug,x​=1 for each g<hg<hg<h. For any family ρ=(ρa,x)a<r,x<s\rho=(\rho_{a,x})_{a<r,x<s}ρ=(ρa,x​)a<r,x<s​ of complex d×dd\times dd×d matrices and any outcome g<hg<hg<h, this defines the real number

Re⁡Tr⁡ ⁣(Kg(1s∑x=0s−1∑a=0r−1ρa,x)Kg†).\operatorname{Re}\operatorname{Tr}\!\left( K_g\left(\frac1s\sum_{x=0}^{s-1}\sum_{a=0}^{r-1}\rho_{a,x}\right)K_g^\dagger \right).ReTr(Kg​(s1​x=0∑s−1​a=0∑r−1​ρa,x​)Kg†​).

The input weights are part of TTT but do not enter this formula. There is no assumption that ρ\rhoρ is positive, normalized, or has equal outcome sums across inputs; consequently no probability bounds are asserted. Empty sums are zero, and the reciprocal of 000 is interpreted as 000.

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