Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for On Mutual Information in Multipartite Quantum States and Equality in Strong Subadditivity of Entropy

Definition
QAlgorithms_Papers_Independent_Intakethird100_P63

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

prove2quantum

The objects the formalized statements of On Mutual Information in Multipartite Quantum States and Equality in Strong Subadditivity of Entropy (arXiv:quant-ph/0311193) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • SpectralExpansion: PDF pp.2,5: a countable spectral resolution of a normalized density operator. Only positive eigenvalues require unit, mutually orthogonal eigenvectors. Zero padding permits finite rank and does not require an infinite orthonormal set in a finite space.
  • spectralKernel: PDF p.2: the operator itself, computed from its spectral resolution. The sum converges in trace norm, since each active projector has trace norm one.
  • herbutEntropy: PDF pp.2,6: intrinsic von Neumann entropy, in nats, including infinity. The defining set is a singleton value on densities, even with spectral degeneracies. The formula -x log x is zero at x=0 and x=1.
  • frameCompression: PDF p.6: finite operator compression in an arbitrary orthonormal frame. Each nested coefficient series is absolutely convergent on density kernels.
  • traceClassDivergence: PDF p.6: full trace-class Umegaki divergence in nats, via all finite orthogonal compressions, retaining infinite values. This is the same complete projection formula as the accepted Dependency50.P13.relativeEntropy, expressed in countable coordinates; its finiteRelativeEntropy definition is imported directly. The rank-zero frame contributes zero. No arbitrary eigenbasis or finite cutoff is fixed.
  • DiscreteDensityEnsemble: PDF pp.5–7: an arbitrary discrete ensemble. The actual component set is the positive-weight support; a zero label carries no normalization or orthogonality premise.
  • herbutMixture: PDF pp.5–7, (11),(14a),(17a): the computed weighted density operator. The norm of each active density is at most its trace norm one; the weighted series converges in trace norm. Inactive kernels contribute zero, regardless of their values.
  • Cluster: PDF p.2: finite sets of parties, numbered from zero in Lean.
  • Coordinates: PDF p.2: the orthonormal coordinates of a completed cluster tensor product. The empty cluster has the unique empty coordinate, hence the scalar tensor factor.
  • restrictCoordinate: PDF p.2: restrict a full tensor coordinate to a cluster.
  • joinCoordinate: PDF p.2: combine retained and discarded coordinates, without changing party order.
  • herbutReduction: PDF p.2: the actual partial trace, computed from the single full density kernel. Positivity and trace one guarantee absolute convergence of every coefficient series.
  • herbutMarginalEntropy: PDF pp.2–3: real marginal entropy, used only under finite-entropy hypotheses.
  • FiniteSubsystemEntropies: PDF pp.2–8: finite-entropy arithmetic domain. Finiteness of all one-party entropies is equivalent to this condition by subadditivity. Empty-cluster entropy is zero. This imposes no bound on Hilbert-space dimension or the number of mixture terms.
  • clusterInformation: PDF p.3, (4b): within-cluster correlation information; singletons contribute zero.
  • clusterMutualInformation: PDF pp.3–4: mutual information of disjoint clusters, using their computed marginals.
  • IsPartition: PDF pp.2–3, (3): unordered partitions into nonempty, pairwise disjoint blocks.
  • amongCorrelation: PDF p.3, (5b): among-cluster correlation information.
  • tensorMarginals: PDF p.2, Lemma 1: the tensor product of the state's own one-party marginals.
  • productAcross: PDF p.4, (8), and Corollary 2: independence of the two complementary blocks. Both blocks may have arbitrary internal correlations.
  • kernelProduct: PDF pp.5–7: ordinary operator multiplication in countable orthonormal coordinates.
  • Biorthogonal: PDF p.5, (10b): biorthogonality on two specified disjoint clusters. The two zero-product conditions are imposed for every pair of distinct active labels.
  • Monoorthogonal: PDF p.7, (17b): only the designated cluster's marginals are orthogonal.
  • BinaryTree: PDF p.3, Corollary 1: syntax for arbitrary successive binary splits.
  • BinaryTree.parties: PDF p.3: the parties in a binary subtree.
  • BinaryTree.Valid: PDF p.3: every split is into disjoint nonempty child clusters. Nonemptiness is inherited from singleton leaves; validity prohibits duplicated parties.
  • BinaryTree.information: PDF p.3: sum of the binary mutual information at every internal node, including the last split into two singletons.
  • attachThird: PDF p.7, Remark 2: product of a two-party density and an arbitrary third factor.
  • productMixture: PDF p.7, Remark 2: the explicit computed mixture of the product components. Third factors may coincide or differ, and zero labels are ignored.
  • pureExpansion: A one-term spectral expansion exists on any unit ray, including finite spaces.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_RealLabelRegisters
import Definitions.Def_QAlgorithms_Defs_SDP
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P33
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P13
import Definitions.Def_QAlgorithms_Papers_Independent_Renew46_P08

/-!
# Herbut: multipartite correlation information and equality in SSA

arXiv:quant-ph/0311193v1. All page numbers below are PDF page numbers.
Countable orthonormal coordinates model completed complex Hilbert tensor products.
Finite dimensional factors embed with zero coefficients; finite Hilbert dimension is
never a theorem hypothesis. Countably many densities share separable local supports.
The accepted Dependency50.P33.IsDensity predicate supplies nuclear trace class,
positivity and trace one. Spectral expansions compute intrinsic entropy on these
operator coefficient matrices; they are not a replacement density predicate. Entropy is the
nonnegative extended spectral sum with natural logarithms. Every expansion of the same
kernel has the same spectrum with multiplicities; the extended infimum below therefore
selects its intrinsic entropy. Empty expansion sets give infinity, and are never used
as states. Real entropy differences are used only on the stated finite-entropy domain.
Lemma 3 is stated with its attribution to Herbut [8], Lemma 9, retained in the metadata.
Compilation and fresh independent readbacks are pending the batch handoff.
-/

noncomputable section

open scoped ENNReal BigOperators

namespace QAlgorithms.Papers.Independent.Intakethird100.P63

attribute [local instance] Classical.propDecidable Classical.decEq

-- Reuse the accepted nuclear trace-class, positive, trace-one predicate verbatim.
open QAlgorithms.Papers.Independent.Dependency50.P33 (IsDensity)
open QAlgorithms.Papers.Independent.Dependency50.P13 (finiteRelativeEntropy)

/-- PDF pp.2,5: a countable spectral resolution of a normalized density operator.
Only positive eigenvalues require unit, mutually orthogonal eigenvectors. Zero padding
permits finite rank and does not require an infinite orthonormal set in a finite space. -/
structure SpectralExpansion (I : Type*) where
  weight : ℕ → ℝ
  vector : ℕ → QAlgorithms.RealRegisters.Hilbert I
  nonnegative : ∀ k, 0 ≤ weight k
  summable : Summable weight
  mass : (∑' k, weight k) = 1
  unit : ∀ k, 0 < weight k → ‖vector k‖ = 1
  orthogonal : ∀ k l, k ≠ l → 0 < weight k → 0 < weight l →
    inner ℂ (vector k) (vector l) = 0

/-- PDF p.2: the operator itself, computed from its spectral resolution.
The sum converges in trace norm, since each active projector has trace norm one. -/
def spectralKernel {I : Type*} (e : SpectralExpansion I) : (Matrix I I ℂ) :=
  fun i j => ∑' k, (e.weight k : ℂ) * e.vector k i * star (e.vector k j)

/-- PDF pp.2,6: intrinsic von Neumann entropy, in nats, including infinity.
The defining set is a singleton value on densities, even with spectral degeneracies.
The formula -x log x is zero at x=0 and x=1. -/
def herbutEntropy {I : Type*} (rho : (Matrix I I ℂ)) : ℝ≥0∞ :=
  sInf {s | ∃ e : SpectralExpansion I, spectralKernel e = rho ∧
    s = ∑' k, ENNReal.ofReal (-e.weight k * Real.log (e.weight k))}

/-- PDF p.6: finite operator compression in an arbitrary orthonormal frame.
Each nested coefficient series is absolutely convergent on density kernels. -/
def frameCompression {I : Type*} {m : ℕ}
    (rho : (Matrix I I ℂ)) (v : Fin m → QAlgorithms.RealRegisters.Hilbert I) :
    Matrix (Fin m) (Fin m) ℂ :=
  fun a b => ∑' i, star (v a i) * (∑' j, rho i j * v b j)

/-- PDF p.6: full trace-class Umegaki divergence in nats, via all finite orthogonal
compressions, retaining infinite values. This is the same complete projection formula
as the accepted Dependency50.P13.relativeEntropy, expressed in countable coordinates;
its finiteRelativeEntropy definition is imported directly.
The rank-zero frame contributes zero. No arbitrary eigenbasis or finite cutoff is fixed. -/
def traceClassDivergence {I : Type*} (rho sigma : (Matrix I I ℂ)) : ℝ≥0∞ :=
  ⨆ m : ℕ, ⨆ v : {v : Fin m → QAlgorithms.RealRegisters.Hilbert I // Orthonormal ℂ v},
    finiteRelativeEntropy (frameCompression rho v.val) (frameCompression sigma v.val)

/-- PDF pp.5–7: an arbitrary discrete ensemble. The actual component set is the
positive-weight support; a zero label carries no normalization or orthogonality premise. -/
structure DiscreteDensityEnsemble (I : Type) where
  weight : ℕ → ℝ
  component : ℕ → (Matrix I I ℂ)
  nonnegative : ∀ k, 0 ≤ weight k
  summable : Summable weight
  mass : (∑' k, weight k) = 1
  density : ∀ k, 0 < weight k → IsDensity (component k)

/-- PDF pp.5–7, (11),(14a),(17a): the computed weighted density operator.
The norm of each active density is at most its trace norm one; the weighted series
converges in trace norm. Inactive kernels contribute zero, regardless of their values. -/
def herbutMixture {I : Type} (e : DiscreteDensityEnsemble I) : (Matrix I I ℂ) :=
  fun i j => ∑' k, (e.weight k : ℂ) * e.component k i j

/-- PDF p.2: finite sets of parties, numbered from zero in Lean. -/
abbrev Cluster (n : ℕ) := Finset (Fin n)

/-- PDF p.2: the orthonormal coordinates of a completed cluster tensor product.
The empty cluster has the unique empty coordinate, hence the scalar tensor factor. -/
abbrev Coordinates {n : ℕ} (C : Cluster n) := C → ℕ

/-- PDF p.2: restrict a full tensor coordinate to a cluster. -/
def restrictCoordinate {n : ℕ} (C : Cluster n) (x : Fin n → ℕ) : Coordinates C :=
  fun i => x i.val

/-- PDF p.2: combine retained and discarded coordinates, without changing party order. -/
def joinCoordinate {n : ℕ} (C : Cluster n)
    (x : Coordinates C) (z : Coordinates Cᶜ) : Fin n → ℕ :=
  fun i => if h : i ∈ C then x ⟨i, h⟩ else z ⟨i, by simpa using h⟩

/-- PDF p.2: the actual partial trace, computed from the single full density kernel.
Positivity and trace one guarantee absolute convergence of every coefficient series. -/
def herbutReduction {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (C : Cluster n) :
    (Matrix (Coordinates C) (Coordinates C) ℂ) :=
  fun x y => ∑' z : Coordinates Cᶜ, rho (joinCoordinate C x z) (joinCoordinate C y z)

/-- PDF pp.2–3: real marginal entropy, used only under finite-entropy hypotheses. -/
def herbutMarginalEntropy {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (C : Cluster n) : ℝ :=
  (herbutEntropy (herbutReduction rho C)).toReal

/-- PDF pp.2–8: finite-entropy arithmetic domain. Finiteness of all one-party
entropies is equivalent to this condition by subadditivity. Empty-cluster entropy is zero.
This imposes no bound on Hilbert-space dimension or the number of mixture terms. -/
def FiniteSubsystemEntropies {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) : Prop :=
  ∀ C : Cluster n, herbutEntropy (herbutReduction rho C) ≠ ⊤

/-- PDF p.3, (4b): within-cluster correlation information; singletons contribute zero. -/
def clusterInformation {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (C : Cluster n) : ℝ :=
  (∑ i ∈ C, herbutMarginalEntropy rho {i}) - herbutMarginalEntropy rho C

/-- PDF pp.3–4: mutual information of disjoint clusters, using their computed marginals. -/
def clusterMutualInformation {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (A B : Cluster n) : ℝ :=
  herbutMarginalEntropy rho A + herbutMarginalEntropy rho B - herbutMarginalEntropy rho (A ∪ B)

/-- PDF pp.2–3, (3): unordered partitions into nonempty, pairwise disjoint blocks. -/
def IsPartition {n : ℕ} (P : Finset (Cluster n)) : Prop :=
  (∀ C ∈ P, C.Nonempty) ∧
    (∀ C ∈ P, ∀ D ∈ P, C ≠ D → Disjoint C D) ∧
    ∀ i : Fin n, ∃ C ∈ P, i ∈ C

/-- PDF p.3, (5b): among-cluster correlation information. -/
def amongCorrelation {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ))
    (P : Finset (Cluster n)) : ℝ :=
  (∑ C ∈ P, herbutMarginalEntropy rho C) - herbutMarginalEntropy rho Finset.univ

/-- PDF p.2, Lemma 1: the tensor product of the state's own one-party marginals. -/
def tensorMarginals {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) :=
  fun x y => ∏ i : Fin n, herbutReduction rho {i} (fun _ => x i) (fun _ => y i)

/-- PDF p.4, (8), and Corollary 2: independence of the two complementary blocks.
Both blocks may have arbitrary internal correlations. -/
def productAcross {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (C : Cluster n) :
    (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) :=
  fun x y => herbutReduction rho C (restrictCoordinate C x) (restrictCoordinate C y) *
    herbutReduction rho Cᶜ (restrictCoordinate Cᶜ x) (restrictCoordinate Cᶜ y)

/-- PDF pp.5–7: ordinary operator multiplication in countable orthonormal coordinates. -/
def kernelProduct {I : Type*} (rho sigma : (Matrix I I ℂ)) : (Matrix I I ℂ) :=
  fun i j => ∑' z, rho i z * sigma z j

/-- PDF p.5, (10b): biorthogonality on two specified disjoint clusters.
The two zero-product conditions are imposed for every pair of distinct active labels. -/
def Biorthogonal {n : ℕ} (e : DiscreteDensityEnsemble (Fin n → ℕ)) (A B : Cluster n) : Prop :=
  ∀ k l, k ≠ l → 0 < e.weight k → 0 < e.weight l →
    kernelProduct (herbutReduction (e.component k) A) (herbutReduction (e.component l) A) = 0 ∧
    kernelProduct (herbutReduction (e.component k) B) (herbutReduction (e.component l) B) = 0

/-- PDF p.7, (17b): only the designated cluster's marginals are orthogonal. -/
def Monoorthogonal {n : ℕ} (e : DiscreteDensityEnsemble (Fin n → ℕ)) (B : Cluster n) : Prop :=
  ∀ k l, k ≠ l → 0 < e.weight k → 0 < e.weight l →
    kernelProduct (herbutReduction (e.component k) B) (herbutReduction (e.component l) B) = 0

/-- PDF p.3, Corollary 1: syntax for arbitrary successive binary splits. -/
inductive BinaryTree (n : ℕ) where
  | leaf (i : Fin n)
  | split (left right : BinaryTree n)

/-- PDF p.3: the parties in a binary subtree. -/
def BinaryTree.parties {n : ℕ} : BinaryTree n → Cluster n
  | .leaf i => {i}
  | .split l r => l.parties ∪ r.parties

/-- PDF p.3: every split is into disjoint nonempty child clusters. Nonemptiness is
inherited from singleton leaves; validity prohibits duplicated parties. -/
def BinaryTree.Valid {n : ℕ} : BinaryTree n → Prop
  | .leaf _ => True
  | .split l r => l.Valid ∧ r.Valid ∧ Disjoint l.parties r.parties

/-- PDF p.3: sum of the binary mutual information at every internal node,
including the last split into two singletons. -/
def BinaryTree.information {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) : BinaryTree n → ℝ
  | .leaf _ => 0
  | .split l r => clusterMutualInformation rho l.parties r.parties + l.information rho + r.information rho

/-- PDF p.7, Remark 2: product of a two-party density and an arbitrary third factor. -/
def attachThird (sigma : (Matrix (Fin 2 → ℕ) (Fin 2 → ℕ) ℂ)) (tau : (Matrix ℕ ℕ ℂ)) : (Matrix (Fin 3 → ℕ) (Fin 3 → ℕ) ℂ) :=
  fun x y => sigma (fun i => x i.castSucc) (fun i => y i.castSucc) * tau (x 2) (y 2)

/-- PDF p.7, Remark 2: the explicit computed mixture of the product components.
Third factors may coincide or differ, and zero labels are ignored. -/
def productMixture (e : DiscreteDensityEnsemble (Fin 2 → ℕ)) (tau : ℕ → (Matrix ℕ ℕ ℂ)) :
    (Matrix (Fin 3 → ℕ) (Fin 3 → ℕ) ℂ) :=
  fun x y => ∑' k, (e.weight k : ℂ) * attachThird (e.component k) (tau k) x y

/-! Sanity checks of local definitions, rather than proofs of source results. -/

/-- A one-term spectral expansion exists on any unit ray, including finite spaces. -/
def pureExpansion {I : Type*} (v : QAlgorithms.RealRegisters.Hilbert I)
    (hv : ‖v‖ = 1) : SpectralExpansion I where
  weight := fun k => if k = 0 then 1 else 0
  vector := fun _ => v
  nonnegative := by intro k; split_ifs <;> norm_num
  summable := (hasSum_ite_eq (0 : ℕ) (1 : ℝ)).summable
  mass := by simp
  unit := by intro k hk; exact hv
  orthogonal := by
    intro k l hkl hk hl
    have hk0 : k = 0 := by by_contra h; simp [h] at hk
    have hl0 : l = 0 := by by_contra h; simp [h] at hl
    exact False.elim (hkl (hk0.trans hl0.symm))

example {I : Type*} (v : QAlgorithms.RealRegisters.Hilbert I) (hv : ‖v‖ = 1) :
    spectralKernel (pureExpansion v hv) = fun i j => v i * star (v j) := by
  funext i j
  change (∑' k : ℕ, ((if k = 0 then 1 else 0 : ℝ) : ℂ) *
    v i * star (v j)) = v i * star (v j)
  have hterm : (fun k : ℕ => ((if k = 0 then 1 else 0 : ℝ) : ℂ) *
      v i * star (v j)) =
      (fun k : ℕ => if k = 0 then v i * star (v j) else 0) := by
    funext k
    by_cases hk : k = 0 <;> simp [hk]
  rw [hterm]
  simp

/-- The concrete vacuum ray meets the unit-vector premise of pureExpansion. -/
example : ‖(lp.single (E := fun _ : ℕ => ℂ) 2 (0 : ℕ) (1 : ℂ))‖ = 1 := by
  simp

/-- Tensor coordinates reassemble exactly, so retained and traced factors are not swapped. -/
example {n : ℕ} (C : Cluster n) (x : Fin n → ℕ) :
    joinCoordinate C (restrictCoordinate C x) (restrictCoordinate Cᶜ x) = x := by
  funext i
  by_cases h : i ∈ C <;> simp [joinCoordinate, restrictCoordinate, h]

example : ENNReal.ofReal (-(0 : ℝ) * Real.log 0) = 0 := by simp

example : ENNReal.ofReal (-(1 : ℝ) * Real.log 1) = 0 := by simp

example {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (i : Fin n) :
    clusterInformation rho {i} = 0 := by simp [clusterInformation]

example {n : ℕ} (rho : (Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ)) (i j : Fin n) :
    (BinaryTree.split (.leaf i) (.leaf j)).information rho = clusterMutualInformation rho {i} {j} := by
  simp [BinaryTree.information, BinaryTree.parties]

example {I : Type*} (rho : (Matrix I I ℂ)) : kernelProduct rho 0 = 0 := by
  funext i j
  simp [kernelProduct]

example {m : ℕ} (X Y : Matrix (Fin m) (Fin m) ℂ)
    (h : ¬ LinearMap.ker Y.mulVecLin ≤ LinearMap.ker X.mulVecLin) :
    finiteRelativeEntropy X Y = ⊤ := by simp [QAlgorithms.Papers.Independent.Dependency50.P13.finiteRelativeEntropy, h]

end QAlgorithms.Papers.Independent.Intakethird100.P63
Source
On Mutual Information in Multipartite Quantum States and Equality in Strong Subadditivity of Entropy, arXiv:quant-ph/0311193v1, https://arxiv.org/abs/quant-ph/0311193v1
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P63.BinaryTree.Valid

For every natural number nnn, this definition assigns a validity proposition to binary trees whose leaves are labelled by {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}. Every leaf is valid. A node with children l,rl,rl,r is valid exactly when both children are valid and their sets of leaf labels are disjoint, where a leaf's label set is its singleton and a node's label set is the union of its children's label sets. Thus validity requires that no label occur in two leaves, but does not require that every possible label occur. For n=0n=0n=0, there are no trees to which the definition applies.

QAlgorithms.Papers.Independent.Intakethird100.P63.BinaryTree.information

For every natural number nnn, every complex kernel ρ\rhoρ indexed by pairs of configurations in Nn\mathbb N^nNn, and every binary tree ttt with leaves labelled by {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, this definition assigns a real number Jρ(t)J_\rho(t)Jρ​(t) recursively by Jρ(leaf⁡(i))=0J_\rho(\operatorname{leaf}(i))=0Jρ​(leaf(i))=0 and

Jρ(split⁡(l,r))=Sρ(Al)+Sρ(Ar)−Sρ(Al∪Ar)+Jρ(l)+Jρ(r),J_\rho(\operatorname{split}(l,r)) =S_\rho(A_l)+S_\rho(A_r)-S_\rho(A_l\cup A_r)+J_\rho(l)+J_\rho(r),Jρ​(split(l,r))=Sρ​(Al​)+Sρ​(Ar​)−Sρ​(Al​∪Ar​)+Jρ​(l)+Jρ​(r),

where AtA_tAt​ is the set of labels occurring in ttt. Here, for a finite coordinate set CCC, define the reduced kernel by RC(x,y)=∑z∈NCcρ(x⊕Cz,y⊕Cz)R_C(x,y)=\sum_{z\in\mathbb N^{C^c}}\rho(x\oplus_C z,y\oplus_C z)RC​(x,y)=∑z∈NCc​ρ(x⊕C​z,y⊕C​z), where x⊕Czx\oplus_C zx⊕C​z agrees with xxx on CCC and with zzz on its complement, and put Sρ(C)=H(RC).toRealS_\rho(C)=H(R_C).\mathrm{toReal}Sρ​(C)=H(RC​).toReal. For a kernel KKK on a label set III, H(K)H(K)H(K) is the infimum of ∑k=0∞max⁡{−wklog⁡wk,0}\sum_{k=0}^\infty\max\{-w_k\log w_k,0\}∑k=0∞​max{−wk​logwk​,0}, interpreted in [0,∞][0,\infty][0,∞], over all representations K(i,j)=∑k=0∞wkuk(i)uk(j)‾K(i,j)=\sum_{k=0}^\infty w_k u_k(i)\overline{u_k(j)}K(i,j)=∑k=0∞​wk​uk​(i)uk​(j)​, with nonnegative summable real weights summing to 111, vectors uk∈ℓ2(I;C)u_k\in\ell^2(I;\mathbb C)uk​∈ℓ2(I;C) of norm 111 whenever wk>0w_k>0wk​>0, and ⟨uk,ul⟩=0\langle u_k,u_l\rangle=0⟨uk​,ul​⟩=0 for distinct indices of positive weight. The infimum is +∞+\infty+∞ if there is no representation, and conversion to a real number sends +∞+\infty+∞ to 000; kernel sums use the total infinite-sum convention, which returns 000 for a nonsummable family. No density, finite-entropy, tree-validity, or coverage assumption is required, so the children's label sets may overlap.

QAlgorithms.Papers.Independent.Intakethird100.P63.BinaryTree

For each natural number nnn, this defines the type of finite trees generated by two constructors: a leaf carrying an index i∈{0,…,n−1}i\in\{0,\ldots,n-1\}i∈{0,…,n−1}, and a split carrying an ordered pair of trees of the same type, designated left and right. Every tree is built using these constructors. Leaf labels may repeat, and there is no requirement that every index occur or that the labels in the two children be disjoint. There is no constructor for an empty tree; consequently, when n=0n=0n=0, this type has no elements.

QAlgorithms.Papers.Independent.Intakethird100.P63.BinaryTree.parties

For every natural number nnn, a binary tree has leaves labelled by elements of {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} and internal nodes with two children. This definition assigns a finite set of labels to each tree: a leaf labelled iii is assigned {i}\{i\}{i}, and a node with children l,rl,rl,r is assigned the union of their assigned sets. Repeated occurrences of a label contribute only one element to this set. No validity or coverage condition is imposed. When n=0n=0n=0, there are no such trees.

QAlgorithms.Papers.Independent.Intakethird100.P63.Biorthogonal

For every natural number nnn, let LLL be the configuration set {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and let eee consist of real weights wkw_kwk​ and complex matrices ρk\rho_kρk​ indexed by LLL, for k∈Nk\in\mathbb Nk∈N, with nonnegative summable weights satisfying ∑k′wk=1\sum'_k w_k=1∑k′​wk​=1. At every positive weight, ρk\rho_kρk​ is required to have a representation ρk(x,y)=∑r′ur(x)vr(y)‾\rho_k(x,y)=\sum'_r u_r(x)\overline{v_r(y)}ρk​(x,y)=∑r′​ur​(x)vr​(y)​ with ur,vr∈ℓ2(L;C)u_r,v_r\in\ell^2(L;\mathbb C)ur​,vr​∈ℓ2(L;C) and ∑r′∥ur∥∥vr∥<∞\sum'_r\|u_r\|\|v_r\|<\infty∑r′​∥ur​∥∥vr​∥<∞, to have every finite principal quadratic form ∑x,y∈sax‾ρk(x,y)ay\sum_{x,y\in s}\overline{a_x}\rho_k(x,y)a_y∑x,y∈s​ax​​ρk​(x,y)ay​ real and nonnegative for every finite s⊆Ls\subseteq Ls⊆L and every choice of complex coefficients, and to have summable real diagonal entries with ∑x′Re⁡ρk(x,x)=1\sum'_x\operatorname{Re}\rho_k(x,x)=1∑x′​Reρk​(x,x)=1; matrices at zero weights have no such requirement. For any subsets A,B⊆{0,…,n−1}A,B\subseteq\{0,\ldots,n-1\}A,B⊆{0,…,n−1}, the defined property means that for every ordered pair k≠lk\ne lk=l with wk>0w_k>0wk​>0 and wl>0w_l>0wl​>0, and for each D∈{A,B}D\in\{A,B\}D∈{A,B}, the matrix whose (a,b)(a,b)(a,b) entry is ∑z:D→N′Rk,D(a,z)Rl,D(z,b)\sum'_{z:D\to\mathbb N}R_{k,D}(a,z)R_{l,D}(z,b)∑z:D→N′​Rk,D​(a,z)Rl,D​(z,b) is zero, where Rk,D(a,b)=∑t:Dc→N′ρk(a⊔t,b⊔t)R_{k,D}(a,b)=\sum'_{t:D^c\to\mathbb N}\rho_k(a\sqcup t,b\sqcup t)Rk,D​(a,b)=∑t:Dc→N′​ρk​(a⊔t,b⊔t) and ⊔\sqcup⊔ joins assignments on complementary subsets. These complex sums are assigned 000 when nonsummable. The subsets need not be nonempty, disjoint, or covering, and n=0n=0n=0 is allowed. Only distinct positive-weight components are tested; if there is just one positive-weight component, the property is vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P63.Cluster

For every natural number nnn, a cluster is a finite subset of the index set {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}. The empty subset is allowed, and when n=0n=0n=0 it is the only cluster.

QAlgorithms.Papers.Independent.Intakethird100.P63.Coordinates

For any natural number nnn and any finite subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, the coordinates on CCC are all functions C→NC\to\mathbb NC→N, assigning an arbitrary natural number to each member of CCC. If CCC is empty, there is exactly one such function, the empty assignment.

QAlgorithms.Papers.Independent.Intakethird100.P63.DiscreteDensityEnsemble

For any type III, a discrete density ensemble consists of real weights wkw_kwk​ and complex matrices Xk:I×I→CX_k:I\times I\to\mathbb CXk​:I×I→C, indexed by k∈Nk\in\mathbb Nk∈N, with wk≥0w_k\ge0wk​≥0, summable weights, and ∑kwk=1\sum_k w_k=1∑k​wk​=1. For every kkk with wk>0w_k>0wk​>0, XkX_kXk​ must satisfy all of the following: there exist sequences ar,br∈ℓ2(I;C)a_r,b_r\in\ell^2(I;\mathbb C)ar​,br​∈ℓ2(I;C) with ∑r∥ar∥∥br∥<∞\sum_r\|a_r\|\|b_r\|<\infty∑r​∥ar​∥∥br​∥<∞ and Xk(i,j)=∑rar(i)br(j)‾X_k(i,j)=\sum_r a_r(i)\overline{b_r(j)}Xk​(i,j)=∑r​ar​(i)br​(j)​ for every i,ji,ji,j; for every finite s⊆Is\subseteq Is⊆I and every z:I→Cz:I\to\mathbb Cz:I→C, the number ∑i∈s∑j∈sz(i)‾Xk(i,j)z(j)\sum_{i\in s}\sum_{j\in s}\overline{z(i)}X_k(i,j)z(j)∑i∈s​∑j∈s​z(i)​Xk​(i,j)z(j) has imaginary part zero and nonnegative real part; and the real diagonal family is summable with ∑i∈IRe⁡Xk(i,i)=1\sum_{i\in I}\operatorname{Re}X_k(i,i)=1∑i∈I​ReXk​(i,i)=1. Components at zero-weight indices are arbitrary matrices. No countability assumption is imposed on III; if III is empty, no ensemble exists because a positive-weight component would have to have diagonal sum 111.

QAlgorithms.Papers.Independent.Intakethird100.P63.FiniteSubsystemEntropies

For every natural number nnn and every complex matrix ρ\rhoρ indexed by configurations {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, this property means that H(RC)≠∞H(R_C)\ne\inftyH(RC​)=∞ for every subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, including the empty and full subsets. Here RC(x,y)=∑z:Cc→N′ρ(x⊔z,y⊔z)R_C(x,y)=\sum'_{z:C^c\to\mathbb N}\rho(x\sqcup z,y\sqcup z)RC​(x,y)=∑z:Cc→N′​ρ(x⊔z,y⊔z) for assignments x,y:C→Nx,y:C\to\mathbb Nx,y:C→N, where x⊔zx\sqcup zx⊔z uses xxx on CCC and zzz on its complement. For a matrix RRR indexed by a label set LLL, H(R)H(R)H(R) is the infimum in [0,∞][0,\infty][0,∞] of ∑k′max⁡{0,−wklog⁡wk}\sum'_k\max\{0,-w_k\log w_k\}∑k′​max{0,−wk​logwk​} over all representations R(a,b)=∑k′wkuk(a)uk(b)‾R(a,b)=\sum'_k w_k u_k(a)\overline{u_k(b)}R(a,b)=∑k′​wk​uk​(a)uk​(b)​, where k∈Nk\in\mathbb Nk∈N, uk∈ℓ2(L;C)u_k\in\ell^2(L;\mathbb C)uk​∈ℓ2(L;C), the real weights are nonnegative and summable with sum 111, and the vectors at positive weights are unit vectors and mutually orthogonal for distinct indices. Vectors at zero weights are unrestricted within ℓ2\ell^2ℓ2. The logarithm is natural, with log⁡0=0\log 0=0log0=0, and an empty infimum is ∞\infty∞. No density or convergence assumption on ρ\rhoρ is part of this property; nonsummable complex sums are assigned 000. When n=0n=0n=0, the property still tests the unique subset, namely ∅\varnothing∅.

QAlgorithms.Papers.Independent.Intakethird100.P63.IsPartition

For every natural number nnn and every finite set PPP of subsets of {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, this property means that every C∈PC\in PC∈P is nonempty, any two distinct members C,D∈PC,D\in PC,D∈P are disjoint, and every i∈{0,…,n−1}i\in\{0,\ldots,n-1\}i∈{0,…,n−1} belongs to some C∈PC\in PC∈P. Thus the members cover the entire index set and are pairwise disjoint. The case n=0n=0n=0 is included: then the property holds exactly when P=∅P=\varnothingP=∅.

QAlgorithms.Papers.Independent.Intakethird100.P63.Monoorthogonal

For every natural number nnn, let LLL be the configuration set {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and let eee consist of real weights wkw_kwk​ and complex matrices ρk\rho_kρk​ indexed by LLL, for k∈Nk\in\mathbb Nk∈N, with nonnegative summable weights satisfying ∑k′wk=1\sum'_k w_k=1∑k′​wk​=1. For every kkk with wk>0w_k>0wk​>0, the matrix ρk\rho_kρk​ must admit a representation ρk(x,y)=∑r′ur(x)vr(y)‾\rho_k(x,y)=\sum'_r u_r(x)\overline{v_r(y)}ρk​(x,y)=∑r′​ur​(x)vr​(y)​ with ur,vr∈ℓ2(L;C)u_r,v_r\in\ell^2(L;\mathbb C)ur​,vr​∈ℓ2(L;C) and summable products ∥ur∥∥vr∥\|u_r\|\|v_r\|∥ur​∥∥vr​∥, every finite principal quadratic form ∑x,y∈sax‾ρk(x,y)ay\sum_{x,y\in s}\overline{a_x}\rho_k(x,y)a_y∑x,y∈s​ax​​ρk​(x,y)ay​ must be real and nonnegative for every finite s⊆Ls\subseteq Ls⊆L and every choice of complex coefficients, and its real diagonal entries must be summable with sum 111; zero-weight components have no such requirement. For any subset B⊆{0,…,n−1}B\subseteq\{0,\ldots,n-1\}B⊆{0,…,n−1}, this property means that for all distinct k,lk,lk,l with wk>0w_k>0wk​>0 and wl>0w_l>0wl​>0, ∑z:B→N′Rk,B(a,z)Rl,B(z,b)=0\sum'_{z:B\to\mathbb N}R_{k,B}(a,z)R_{l,B}(z,b)=0∑z:B→N′​Rk,B​(a,z)Rl,B​(z,b)=0 for every pair of assignments a,b:B→Na,b:B\to\mathbb Na,b:B→N, where Rk,B(a,b)=∑t:Bc→N′ρk(a⊔t,b⊔t)R_{k,B}(a,b)=\sum'_{t:B^c\to\mathbb N}\rho_k(a\sqcup t,b\sqcup t)Rk,B​(a,b)=∑t:Bc→N′​ρk​(a⊔t,b⊔t) and ⊔\sqcup⊔ joins assignments on complementary subsets. Nonsummable complex sums are assigned 000. The empty subset and n=0n=0n=0 are included. No condition is tested between a component and itself or on zero-weight components; with just one positive-weight component, the property is vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P63.SpectralExpansion

For any type III, a spectral expansion consists of real weights wkw_kwk​ and vectors uk∈ℓ2(I;C)u_k\in\ell^2(I;\mathbb C)uk​∈ℓ2(I;C), indexed by k∈Nk\in\mathbb Nk∈N, such that wk≥0w_k\ge 0wk​≥0 for every kkk, the sequence of weights is summable, ∑k=0∞wk=1\sum_{k=0}^{\infty}w_k=1∑k=0∞​wk​=1, every vector with wk>0w_k>0wk​>0 has norm 111, and ⟨uk,ul⟩=0\langle u_k,u_l\rangle=0⟨uk​,ul​⟩=0 whenever k≠lk\ne lk=l and both weights are positive. Vectors at zero-weight indices are arbitrary elements of ℓ2(I;C)\ell^2(I;\mathbb C)ℓ2(I;C); neither normalization nor orthogonality is required of them. There is no countability assumption on III. If III is empty, no such expansion exists, since the mass condition requires a positive weight, whereas ℓ2(I;C)\ell^2(I;\mathbb C)ℓ2(I;C) has no unit vector.

QAlgorithms.Papers.Independent.Intakethird100.P63.amongCorrelation

For every natural number nnn, every complex matrix ρ\rhoρ indexed by configurations {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and every finite set PPP of subsets of {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, this defines ∑C∈PSρ(C)−Sρ({0,…,n−1})\sum_{C\in P}S_\rho(C)-S_\rho(\{0,\ldots,n-1\})∑C∈P​Sρ​(C)−Sρ​({0,…,n−1}). Here Sρ(D)=H(RD).toRealS_\rho(D)=H(R_D).\mathrm{toReal}Sρ​(D)=H(RD​).toReal, where RD(x,y)=∑z:Dc→N′ρ(x⊔z,y⊔z)R_D(x,y)=\sum'_{z:D^c\to\mathbb N}\rho(x\sqcup z,y\sqcup z)RD​(x,y)=∑z:Dc→N′​ρ(x⊔z,y⊔z) for assignments x,y:D→Nx,y:D\to\mathbb Nx,y:D→N, with ⊔\sqcup⊔ joining assignments on complementary subsets. For a matrix RRR indexed by LLL, H(R)H(R)H(R) is the infimum in [0,∞][0,\infty][0,∞] of ∑k′max⁡{0,−wklog⁡wk}\sum'_k\max\{0,-w_k\log w_k\}∑k′​max{0,−wk​logwk​} over representations R(a,b)=∑k′wkuk(a)uk(b)‾R(a,b)=\sum'_k w_k u_k(a)\overline{u_k(b)}R(a,b)=∑k′​wk​uk​(a)uk​(b)​, where uk∈ℓ2(L;C)u_k\in\ell^2(L;\mathbb C)uk​∈ℓ2(L;C) and the nonnegative summable real weights have sum 111, with unit vectors at positive weights and orthogonality between distinct positive-weight vectors; zero-weight vectors are unrestricted within ℓ2\ell^2ℓ2. The logarithm is natural, with log⁡0=0\log 0=0log0=0. An empty infimum is ∞\infty∞, and conversion to a real number sends ∞\infty∞ to 000. Neither a density or finite-entropy assumption on ρ\rhoρ nor a partition condition on PPP is required; nonsummable complex sums are assigned 000. In particular, P=∅P=\varnothingP=∅ gives the negative of the full-subset marginal entropy, and n=0n=0n=0 is allowed.

QAlgorithms.Papers.Independent.Intakethird100.P63.attachThird

For any complex kernel σ\sigmaσ indexed by pairs of configurations in N2\mathbb N^2N2 and any complex kernel τ\tauτ indexed by pairs of natural numbers, this definition produces a kernel on N3\mathbb N^3N3 whose entries are

attachThird⁡(σ,τ)(x,y)=σ((x0,x1),(y0,y1))τ(x2,y2).\operatorname{attachThird}(\sigma,\tau)(x,y) =\sigma\bigl((x_0,x_1),(y_0,y_1)\bigr)\tau(x_2,y_2).attachThird(σ,τ)(x,y)=σ((x0​,x1​),(y0​,y1​))τ(x2​,y2​).

There are no positivity, normalization, summability, or density assumptions on either input.

QAlgorithms.Papers.Independent.Intakethird100.P63.clusterInformation

For every natural number nnn, every complex matrix ρ\rhoρ indexed by configurations {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and every subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, this defines ∑i∈CSρ({i})−Sρ(C)\sum_{i\in C}S_\rho(\{i\})-S_\rho(C)∑i∈C​Sρ​({i})−Sρ​(C). For any subset DDD, put RD(x,y)=∑z:Dc→N′ρ(x⊔z,y⊔z)R_D(x,y)=\sum'_{z:D^c\to\mathbb N}\rho(x\sqcup z,y\sqcup z)RD​(x,y)=∑z:Dc→N′​ρ(x⊔z,y⊔z) for assignments x,y:D→Nx,y:D\to\mathbb Nx,y:D→N, with ⊔\sqcup⊔ joining assignments on complementary subsets, and let Sρ(D)=H(RD).toRealS_\rho(D)=H(R_D).\mathrm{toReal}Sρ​(D)=H(RD​).toReal. For a matrix RRR indexed by LLL, H(R)H(R)H(R) is the infimum in [0,∞][0,\infty][0,∞] of ∑k′max⁡{0,−wklog⁡wk}\sum'_k\max\{0,-w_k\log w_k\}∑k′​max{0,−wk​logwk​} over representations R(a,b)=∑k′wkuk(a)uk(b)‾R(a,b)=\sum'_k w_k u_k(a)\overline{u_k(b)}R(a,b)=∑k′​wk​uk​(a)uk​(b)​ by vectors uk∈ℓ2(L;C)u_k\in\ell^2(L;\mathbb C)uk​∈ℓ2(L;C) and nonnegative summable real weights of sum 111, with unit vectors at positive weights and orthogonality between distinct positive-weight vectors; zero-weight vectors have no such restrictions. The logarithm is natural and log⁡0=0\log 0=0log0=0. An empty infimum is ∞\infty∞, and conversion to a real number sends ∞\infty∞ to 000. No density, finiteness-of-entropy, or convergence assumption on ρ\rhoρ is imposed; nonsummable complex sums are assigned 000. In particular, C=∅C=\varnothingC=∅ gives −Sρ(∅)-S_\rho(\varnothing)−Sρ​(∅), and n=0n=0n=0 is allowed.

QAlgorithms.Papers.Independent.Intakethird100.P63.clusterMutualInformation

For every natural number nnn, every complex matrix ρ\rhoρ indexed by configurations {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and any two subsets A,B⊆{0,…,n−1}A,B\subseteq\{0,\ldots,n-1\}A,B⊆{0,…,n−1}, this defines Sρ(A)+Sρ(B)−Sρ(A∪B)S_\rho(A)+S_\rho(B)-S_\rho(A\cup B)Sρ​(A)+Sρ​(B)−Sρ​(A∪B). Here Sρ(D)=H(RD).toRealS_\rho(D)=H(R_D).\mathrm{toReal}Sρ​(D)=H(RD​).toReal, where RD(x,y)=∑z:Dc→N′ρ(x⊔z,y⊔z)R_D(x,y)=\sum'_{z:D^c\to\mathbb N}\rho(x\sqcup z,y\sqcup z)RD​(x,y)=∑z:Dc→N′​ρ(x⊔z,y⊔z) for assignments x,y:D→Nx,y:D\to\mathbb Nx,y:D→N, and ⊔\sqcup⊔ joins assignments on complementary subsets. For a matrix RRR indexed by LLL, H(R)H(R)H(R) is the infimum in [0,∞][0,\infty][0,∞] of ∑k′max⁡{0,−wklog⁡wk}\sum'_k\max\{0,-w_k\log w_k\}∑k′​max{0,−wk​logwk​} over representations R(a,b)=∑k′wkuk(a)uk(b)‾R(a,b)=\sum'_k w_k u_k(a)\overline{u_k(b)}R(a,b)=∑k′​wk​uk​(a)uk​(b)​ with uk∈ℓ2(L;C)u_k\in\ell^2(L;\mathbb C)uk​∈ℓ2(L;C) and nonnegative summable real weights of sum 111, such that every positive-weight vector is unit and vectors at distinct positive weights are orthogonal; zero-weight vectors need not be unit or orthogonal. The logarithm is natural, with log⁡0=0\log 0=0log0=0; an empty infimum is ∞\infty∞, whose conversion to a real number is 000. There is no density, finite-entropy, or convergence hypothesis on ρ\rhoρ, and nonsummable complex sums are assigned 000. The subsets may be empty, overlap, or coincide, and n=0n=0n=0 is included.

QAlgorithms.Papers.Independent.Intakethird100.P63.frameCompression

For any type III, any natural number mmm, any complex matrix ρ:I×I→C\rho:I\times I\to\mathbb Cρ:I×I→C, and any family v0,…,vm−1∈ℓ2(I;C)v_0,\ldots,v_{m-1}\in\ell^2(I;\mathbb C)v0​,…,vm−1​∈ℓ2(I;C), the defined compression is the m×mm\times mm×m matrix

Cρ,v(a,b)=∑i∈Iva(i)‾(∑j∈Iρ(i,j)vb(j))(0≤a,b<m).C_{\rho,v}(a,b) = \sum_{i\in I}\overline{v_a(i)} \left(\sum_{j\in I}\rho(i,j)v_b(j)\right) \qquad(0\le a,b<m).Cρ,v​(a,b)=i∈I∑​va​(i)​​j∈I∑​ρ(i,j)vb​(j)​(0≤a,b<m).

There is no assumption that the family is orthonormal, that ρ\rhoρ is a density matrix, or that either displayed sum is summable. These are nested total infinite sums in the displayed order; a nonsummable complex family has total sum zero. The index type III need not be countable. The value m=0m=0m=0 is allowed and gives the matrix with no entries.

QAlgorithms.Papers.Independent.Intakethird100.P63.herbutEntropy

For any type III and any complex matrix ρ:I×I→C\rho:I\times I\to\mathbb Cρ:I×I→C, its defined entropy is the extended nonnegative number

inf⁡{∑k=0∞max⁡{−wklog⁡wk,0}:ρ(i,j)=∑k=0∞wkuk(i)uk(j)‾ for every i,j∈I},\inf\left\{ \sum_{k=0}^{\infty}\max\{-w_k\log w_k,0\}: \rho(i,j)=\sum_{k=0}^{\infty}w_k u_k(i)\overline{u_k(j)} \text{ for every }i,j\in I \right\},inf{k=0∑∞​max{−wk​logwk​,0}:ρ(i,j)=k=0∑∞​wk​uk​(i)uk​(j)​ for every i,j∈I},

where the infimum ranges over all summable sequences wk≥0w_k\ge0wk​≥0 with ∑kwk=1\sum_k w_k=1∑k​wk​=1 and all sequences uk∈ℓ2(I;C)u_k\in\ell^2(I;\mathbb C)uk​∈ℓ2(I;C) satisfying ∥uk∥=1\|u_k\|=1∥uk​∥=1 when wk>0w_k>0wk​>0 and ⟨uk,ul⟩=0\langle u_k,u_l\rangle=0⟨uk​,ul​⟩=0 for distinct indices with positive weights. Each nonnegative summand is regarded as an element of [0,∞][0,\infty][0,∞], and log⁡\loglog is the natural logarithm. Zero weights contribute zero, using the total real logarithm with log⁡0=0\log 0=0log0=0, and impose no normalization or orthogonality requirements on their vectors. The sum may be infinite. No density hypothesis is imposed on ρ\rhoρ; if it has no such expansion, the infimum is over the empty set and equals +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P63.herbutMarginalEntropy

For every natural number nnn, every complex matrix ρ\rhoρ indexed by pairs of configurations in N{0,…,n−1}\mathbb N^{\{0,\ldots,n-1\}}N{0,…,n−1}, and every subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, this defines a real number as follows. Write CcC^cCc for the complement in {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1} and define the matrix RCR_CRC​, indexed by assignments x,y:C→Nx,y:C\to\mathbb Nx,y:C→N, by RC(x,y)=∑z:Cc→N′ρ(x⊔z,y⊔z)R_C(x,y)=\sum'_{z:C^c\to\mathbb N}\rho(x\sqcup z,y\sqcup z)RC​(x,y)=∑z:Cc→N′​ρ(x⊔z,y⊔z), where x⊔zx\sqcup zx⊔z agrees with xxx on CCC and with zzz on CcC^cCc. For any matrix RRR with label set LLL, define H(R)H(R)H(R) to be the infimum, in [0,∞][0,\infty][0,∞], of ∑k∈N′max⁡{0,−wklog⁡wk}\sum'_{k\in\mathbb N}\max\{0,-w_k\log w_k\}∑k∈N′​max{0,−wk​logwk​} over all representations R(a,b)=∑k′wkuk(a)uk(b)‾R(a,b)=\sum'_k w_k u_k(a)\overline{u_k(b)}R(a,b)=∑k′​wk​uk​(a)uk​(b)​ with uk∈ℓ2(L;C)u_k\in\ell^2(L;\mathbb C)uk​∈ℓ2(L;C), nonnegative summable real weights satisfying ∑k′wk=1\sum'_k w_k=1∑k′​wk​=1, unit vectors at every positive weight, and orthogonality between vectors at distinct positive weights; vectors at zero weights have no normalization or orthogonality requirement. The defined marginal entropy is H(RC).toRealH(R_C).\mathrm{toReal}H(RC​).toReal: its finite value if H(RC)<∞H(R_C)<\inftyH(RC​)<∞, and 000 if H(RC)=∞H(R_C)=\inftyH(RC​)=∞, including when there is no such representation and the infimum is over the empty set. Here log⁡\loglog is the natural logarithm, with log⁡0=0\log 0=0log0=0. No density, positivity, normalization, or convergence hypothesis on ρ\rhoρ is imposed; a nonsummable complex sum ∑′\sum'∑′ has value 000. Both C=∅C=\varnothingC=∅ and n=0n=0n=0 are included.

QAlgorithms.Papers.Independent.Intakethird100.P63.herbutMixture

For any type III and a discrete density ensemble eee with summable nonnegative weights wkw_kwk​ totaling 111, the mixture is the complex matrix

Me(i,j)=∑k=0∞wkXk(i,j)(i,j∈I).M_e(i,j)=\sum_{k=0}^{\infty}w_kX_k(i,j) \qquad(i,j\in I).Me​(i,j)=k=0∑∞​wk​Xk​(i,j)(i,j∈I).

Here, for each positive-weight index kkk, the component XkX_kXk​ has a representation Xk(i,j)=∑rar(i)br(j)‾X_k(i,j)=\sum_r a_r(i)\overline{b_r(j)}Xk​(i,j)=∑r​ar​(i)br​(j)​ with ar,br∈ℓ2(I;C)a_r,b_r\in\ell^2(I;\mathbb C)ar​,br​∈ℓ2(I;C) and summable ∥ar∥∥br∥\|a_r\|\|b_r\|∥ar​∥∥br​∥; every finite quadratic form ∑i,j∈sz(i)‾Xk(i,j)z(j)\sum_{i,j\in s}\overline{z(i)}X_k(i,j)z(j)∑i,j∈s​z(i)​Xk​(i,j)z(j), for finite s⊆Is\subseteq Is⊆I and arbitrary z:I→Cz:I\to\mathbb Cz:I→C, is real and nonnegative; and its real diagonal is summable with total 111. Zero-weight components are unrestricted and contribute zero. The displayed entrywise series uses the total complex sum, whose value is zero if its summands are not summable; no separate summability hypothesis on this series is supplied.

QAlgorithms.Papers.Independent.Intakethird100.P63.herbutReduction

For any natural number nnn, any arbitrary complex matrix ρ\rhoρ indexed in both entries by assignments {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and any finite subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, the reduction is the matrix indexed by assignments x,y:C→Nx,y:C\to\mathbb Nx,y:C→N defined by

RC(x,y)=∑z:Cc→Nρ(JC(x,z),JC(y,z)),R_C(x,y)= \sum_{z:C^{\mathrm c}\to\mathbb N} \rho\bigl(J_C(x,z),J_C(y,z)\bigr),RC​(x,y)=z:Cc→N∑​ρ(JC​(x,z),JC​(y,z)),

where CcC^{\mathrm c}Cc is the complement within {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, and JC(x,z)J_C(x,z)JC​(x,z) equals xxx on CCC and zzz on its complement. No positivity, density, trace-class, or summability hypothesis is imposed on ρ\rhoρ; the displayed sum is the total complex sum and equals zero if its summands are not summable. If CCC is empty, the reduction has one entry, given by the sum of all diagonal entries of ρ\rhoρ. If CCC is the full index set, there is one complement assignment, so the reduction returns ρ\rhoρ under the natural identification of coordinates; these conventions also cover n=0n=0n=0.

QAlgorithms.Papers.Independent.Intakethird100.P63.joinCoordinate

For any natural number nnn, any finite subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, and assignments x:C→Nx:C\to\mathbb Nx:C→N and z:Cc→Nz:C^{\mathrm c}\to\mathbb Nz:Cc→N, where the complement is taken within {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, the joined assignment JC(x,z):{0,…,n−1}→NJ_C(x,z):\{0,\ldots,n-1\}\to\mathbb NJC​(x,z):{0,…,n−1}→N is defined by

JC(x,z)(i)={x(i),i∈C,z(i),i∉C.J_C(x,z)(i)= \begin{cases} x(i),&i\in C,\\ z(i),&i\notin C. \end{cases}JC​(x,z)(i)={x(i),z(i),​i∈C,i∈/C.​

Empty clusters and empty complements are allowed; in particular, for n=0n=0n=0 the result is the unique empty assignment.

QAlgorithms.Papers.Independent.Intakethird100.P63.kernelProduct

For any type III and any two complex matrices ρ,σ\rho,\sigmaρ,σ indexed by I×II\times II×I, this defines the matrix K(i,j)=∑z∈I′ρ(i,z)σ(z,j)K(i,j)=\sum'_{z\in I}\rho(i,z)\sigma(z,j)K(i,j)=∑z∈I′​ρ(i,z)σ(z,j) for every i,j∈Ii,j\in Ii,j∈I. There is no finiteness, countability, positivity, or summability assumption on III or the matrices. The sum is the unconditional sum over III, assigned the value 000 when its summands are not summable. The empty label set is also allowed.

QAlgorithms.Papers.Independent.Intakethird100.P63.productAcross

For every natural number nnn, every complex matrix ρ\rhoρ indexed by configurations {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, and every subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, this defines a matrix on the same configurations by PC(x,y)=RC(x∣C,y∣C) RCc(x∣Cc,y∣Cc)P_C(x,y)=R_C(x|_C,y|_C)\,R_{C^c}(x|_{C^c},y|_{C^c})PC​(x,y)=RC​(x∣C​,y∣C​)RCc​(x∣Cc​,y∣Cc​). Here the complement is taken in {0,…,n−1}\{0,\ldots,n-1\}{0,…,n−1}, restrictions retain the indicated coordinates, and, for each subset DDD, RD(a,b)=∑z:Dc→N′ρ(a⊔z,b⊔z)R_D(a,b)=\sum'_{z:D^c\to\mathbb N}\rho(a\sqcup z,b\sqcup z)RD​(a,b)=∑z:Dc→N′​ρ(a⊔z,b⊔z) for assignments a,b:D→Na,b:D\to\mathbb Na,b:D→N, with ⊔\sqcup⊔ joining assignments on complementary subsets. No density, positivity, normalization, or convergence hypothesis on ρ\rhoρ is imposed; nonsummable complex sums are assigned 000. Empty and full choices of CCC are permitted, as is n=0n=0n=0; assignments on an empty subset form a singleton.

QAlgorithms.Papers.Independent.Intakethird100.P63.productMixture

Let eee consist of real weights wkw_kwk​ and complex kernels KkK_kKk​ on N2\mathbb N^2N2, indexed by k∈Nk\in\mathbb Nk∈N, with wk≥0w_k\ge0wk​≥0, the weights summable, and ∑kwk=1\sum_k w_k=1∑k​wk​=1. Whenever wk>0w_k>0wk​>0, KkK_kKk​ must have a representation Kk(i,j)=∑rar(i)br(j)‾K_k(i,j)=\sum_r a_r(i)\overline{b_r(j)}Kk​(i,j)=∑r​ar​(i)br​(j)​ with ar,br∈ℓ2(N2;C)a_r,b_r\in\ell^2(\mathbb N^2;\mathbb C)ar​,br​∈ℓ2(N2;C) and ∑r∥ar∥∥br∥<∞\sum_r\|a_r\|\|b_r\|<\infty∑r​∥ar​∥∥br​∥<∞; every finite quadratic form ∑i,j∈Fzi‾Kk(i,j)zj\sum_{i,j\in F}\overline{z_i}K_k(i,j)z_j∑i,j∈F​zi​​Kk​(i,j)zj​, for finite F⊆N2F\subseteq\mathbb N^2F⊆N2 and arbitrary z:N2→Cz:\mathbb N^2\to\mathbb Cz:N2→C, must be real and nonnegative; and its real diagonal must be summable with sum 111. For every family of arbitrary complex kernels τk\tau_kτk​ on N\mathbb NN, the definition produces the kernel

productMixture⁡(e,τ)(x,y)=∑k=0∞wkKk((x0,x1),(y0,y1))τk(x2,y2),x,y∈N3.\operatorname{productMixture}(e,\tau)(x,y) =\sum_{k=0}^\infty w_k K_k\bigl((x_0,x_1),(y_0,y_1)\bigr)\tau_k(x_2,y_2), \qquad x,y\in\mathbb N^3.productMixture(e,τ)(x,y)=k=0∑∞​wk​Kk​((x0​,x1​),(y0​,y1​))τk​(x2​,y2​),x,y∈N3.

The components KkK_kKk​ with zero weight have no density requirement, and there are no assumptions on the τk\tau_kτk​. Each displayed series is the total complex infinite sum, equal to 000 if that entry's series is not summable; the definition itself asserts no density property of the resulting kernel.

QAlgorithms.Papers.Independent.Intakethird100.P63.pureExpansion

For every label type III, every vector v∈ℓ2(I;C)v\in\ell^2(I;\mathbb C)v∈ℓ2(I;C), and an assumption ∥v∥=1\|v\|=1∥v∥=1, this definition constructs a spectral expansion with weights w0=1w_0=1w0​=1 and wk=0w_k=0wk​=0 for every k≠0k\ne0k=0, and with vector vvv at every index. The resulting weights are nonnegative and summable with sum 111; its vectors have norm 111 at every positive-weight index; and vectors at distinct positive-weight indices are required to be orthogonal. That last requirement is vacuous here because 000 is the only positive-weight index. No countability or nonemptiness assumption is imposed on III; if III is empty, the required unit vector cannot exist.

QAlgorithms.Papers.Independent.Intakethird100.P63.restrictCoordinate

For any natural number nnn, any finite subset C⊆{0,…,n−1}C\subseteq\{0,\ldots,n-1\}C⊆{0,…,n−1}, and any assignment x:{0,…,n−1}→Nx:\{0,\ldots,n-1\}\to\mathbb Nx:{0,…,n−1}→N, the restricted coordinate assignment is the function on CCC sending each i∈Ci\in Ci∈C to x(i)x(i)x(i). This includes the unique empty assignment when C=∅C=\varnothingC=∅, and includes n=0n=0n=0.

QAlgorithms.Papers.Independent.Intakethird100.P63.spectralKernel

For any type III and any spectral expansion eee, consisting of summable nonnegative weights wkw_kwk​ of total mass 111 and vectors uk∈ℓ2(I;C)u_k\in\ell^2(I;\mathbb C)uk​∈ℓ2(I;C) that are unit vectors and pairwise orthogonal at distinct positive-weight indices, the associated matrix Ke:I×I→CK_e:I\times I\to\mathbb CKe​:I×I→C is defined entrywise by

Ke(i,j)=∑k=0∞wk uk(i) uk(j)‾(i,j∈I).K_e(i,j)=\sum_{k=0}^{\infty}w_k\,u_k(i)\,\overline{u_k(j)} \qquad(i,j\in I).Ke​(i,j)=k=0∑∞​wk​uk​(i)uk​(j)​(i,j∈I).

Zero-weight vectors remain unrestricted apart from belonging to ℓ2(I;C)\ell^2(I;\mathbb C)ℓ2(I;C), and their summands are zero.

QAlgorithms.Papers.Independent.Intakethird100.P63.tensorMarginals

For every natural number nnn and every complex matrix ρ\rhoρ indexed by configurations {0,…,n−1}→N\{0,\ldots,n-1\}\to\mathbb N{0,…,n−1}→N, this defines another matrix on the same configurations by T(x,y)=∏i=0n−1R{i}(x∣{i},y∣{i})T(x,y)=\prod_{i=0}^{n-1}R_{\{i\}}(x|_{\{i\}},y|_{\{i\}})T(x,y)=∏i=0n−1​R{i}​(x∣{i}​,y∣{i}​). For any subset DDD, the matrix appearing here is RD(a,b)=∑z:Dc→N′ρ(a⊔z,b⊔z)R_D(a,b)=\sum'_{z:D^c\to\mathbb N}\rho(a\sqcup z,b\sqcup z)RD​(a,b)=∑z:Dc→N′​ρ(a⊔z,b⊔z), where a,b:D→Na,b:D\to\mathbb Na,b:D→N and the joined configuration uses the first assignment on DDD and zzz on its complement; thus the singleton assignments have values xix_ixi​ and yiy_iyi​. No positivity, normalization, density, or summability assumption on ρ\rhoρ is required, and a nonsummable complex sum is assigned 000. When n=0n=0n=0, the product is empty and equals 111, so the resulting matrix has its unique entry equal to 111.

QAlgorithms.Papers.Independent.Intakethird100.P63.traceClassDivergence

For any type III and arbitrary complex matrices ρ,σ:I×I→C\rho,\sigma:I\times I\to\mathbb Cρ,σ:I×I→C, the defined divergence is

sup⁡m∈N sup⁡v0,…,vm−1∈ℓ2(I;C)⟨va,vb⟩=δabDfin(Cρ,v,Cσ,v),\sup_{m\in\mathbb N}\ \sup_{\substack{v_0,\ldots,v_{m-1}\in\ell^2(I;\mathbb C)\\ \langle v_a,v_b\rangle=\delta_{ab}}} D_{\mathrm{fin}}(C_{\rho,v},C_{\sigma,v}),m∈Nsup​ v0​,…,vm−1​∈ℓ2(I;C)⟨va​,vb​⟩=δab​​sup​Dfin​(Cρ,v​,Cσ,v​),

with values in [0,∞][0,\infty][0,∞], where

Cρ,v(a,b)=∑i∈Iva(i)‾(∑j∈Iρ(i,j)vb(j)),C_{\rho,v}(a,b)= \sum_{i\in I}\overline{v_a(i)} \left(\sum_{j\in I}\rho(i,j)v_b(j)\right),Cρ,v​(a,b)=i∈I∑​va​(i)​​j∈I∑​ρ(i,j)vb​(j)​,

and similarly for σ\sigmaσ. For finite square matrices X,YX,YX,Y, Dfin(X,Y)=+∞D_{\mathrm{fin}}(X,Y)=+\inftyDfin​(X,Y)=+∞ if ker⁡Y⊈ker⁡X\ker Y\not\subseteq\ker XkerY⊆kerX; otherwise it is

max⁡ ⁣{Re⁡tr⁡(f(X)−Xg(Y)−X+Y), 0},f(t)=tlog⁡t,g(t)=log⁡t,\max\!\left\{ \operatorname{Re}\operatorname{tr} \bigl(f(X)-Xg(Y)-X+Y\bigr),\,0 \right\}, \qquad f(t)=t\log t,\quad g(t)=\log t,max{Retr(f(X)−Xg(Y)−X+Y),0},f(t)=tlogt,g(t)=logt,

regarded as an extended nonnegative number, with f(X)f(X)f(X) and g(Y)g(Y)g(Y) evaluated by the real continuous functional calculus used in the definition. No trace-class, positivity, normalization, or countability hypotheses are imposed on the input matrices or on III. The compression sums are total sums, taking value zero for nonsummable complex families. The supremum includes m=0m=0m=0, whose empty frame contributes zero; for any mmm admitting no orthonormal frame, the corresponding inner supremum is zero.

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