Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Inequalities for Quantum Entropy: A Review with Conditions for Equality

Definition
QAlgorithms_Papers_Independent_Intakethird100_P51

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

prove2quantum

The objects the formalized statements of Inequalities for Quantum Entropy: A Review with Conditions for Equality (arXiv:quant-ph/0205064) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • logMatrix: Natural Hermitian logarithm (PDF pp.4–5). At a zero eigenvalue the value is zero; this is used only inside x log x, or as zero extension of a common-support logarithm difference. It does not assert that a singular matrix has a finite logarithm.
  • SupportedBy: Support inclusion supp X ⊆ supp Y, PDF p.4.
  • FaithfulDensity: Default strictly positive normalized density convention, PDF pp.3,5.
  • liftAB: AB operator lifted by the unnormalized identity on C (PDF p.4).
  • liftBC: BC operator lifted by the identity on A (PDF p.4).
  • liftB: B operator lifted by identities on A and C (PDF p.4).
  • ssaGap: Right minus left side of SSA (3), PDF p.3, with actual marginals.
  • SSAlogCondition: The literal lifted logarithm equality (4), PDF p.4.
  • imaginaryExp: exp(itL), with the natural-log imaginary-power convention, PDF p.7.
  • splitProduct: Regrouping of rho_(1,2a) tensor rho_(2b,3), PDF p.7.
  • marginalABC: ABC marginal of the four-register purification, PDF p.9.
  • marginalABD: ABD marginal of the same four-register purification, PDF p.9.
  • WeakMonotonicity: Inequality (22) on the actual 124 marginal (PDF pp.9,15).
  • krausV: The computed rectangular Kraus isometry in (26), PDF pp.10–11.
  • coordinatePVM: Coordinate projections I tensor |b><b| in Section 3.5, PDF p.11.
  • basisPVM: Common-basis rank-one spectral measurement, PDF p.17.
  • inverseRoot: Inverse square root for the faithful ensemble–POVM duality, PDF p.21.
  • dualEffects: Dual POVM, PDF p.21. Its uses require a faithful reference and positive weights.
  • dualUnnormalized: Unnormalized dual ensemble component sqrt(rho) E_b sqrt(rho), PDF p.21.
  • dualStates: Dual normalized state at an active outcome tau_b>0, PDF p.21. Its zero-mass value is zero and is never treated as a density matrix or placed inside a log.
  • hallBound: Hall's right side (59), PDF p.21, explicitly omitting zero-mass outcomes.
  • hallRelative: Equation (60), PDF p.21, retaining the factor tau_b.
  • UEP: The optimizing expression (61), PDF p.21. Only the definition is retained; no strict comparison between global optima is inferred from one maximizing pair.
  • MONO: MONO as a global finite PSD/support principle (PDF pp.14–15). Unnormalized positive matrices are included, so embedded references used in the reductions are covered.
  • MPT: MPT as a global principle, with genuine partial traces, PDF pp.13–15.
  • SSA: SSA as a global normalized PSD inequality, including block and pure states, PDF p.15.
  • JC: JC as a global finite PSD/support principle, PDF pp.14–15.
  • UpperLogContour: Admissible positively oriented circle for the principal analytic matrix logarithm, equation (63), PDF p.24. Its closed disk is in the scalar upper half-plane and its open disk contains the entire spectrum, ensuring winding number one and no resolvent pole.
  • principalLog: Principal holomorphic matrix logarithm on upper-half-plane spectrum, using exactly the contour-resolvent calculus (63), PDF p.24. Different admissible circles give the same value by holomorphic functional calculus. Outside this domain the zero value is a totalization only, never used by the Pick theorem.
  • epsteinTrace: Epstein's nonnormal matrix pencil trace function (62), PDF p.24.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_Ensemble
import Definitions.Def_QAlgorithms_Defs_Measurement
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P27
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P56
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P63
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P21
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P35
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P19
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P33

noncomputable section

open scoped Matrix ComplexOrder ComplexConjugate Matrix.Norms.L2Operator
open MeasureTheory

namespace QAlgorithms.Papers.Independent.Intakethird100.P51

attribute [local instance] Classical.propDecidable Classical.decEq

-- Exact interfaces from the frozen accepted library; no local redeclarations.
open QAlgorithms.Papers.Independent.Information100.P63 (Mat)
open QAlgorithms.Papers.Independent.Information100.P56 (entropy holevo IsLaw)
open QAlgorithms.Papers.Independent.Information100.P27 (average)
open QAlgorithms.Papers.Independent.Dependency50.P21 (relativeEntropy marginalAB channelAdjoint)
open QAlgorithms.Papers.Independent.Dependency50.P35 (marginalAC)
open QAlgorithms.Papers.Independent.Gap48.P19 (shannonNats)
open QAlgorithms.Papers.Independent.Gap48.P33 (measurementInformation)

/- The accepted relativeEntropy is EReal-valued in nats. Every use below is finite
by positive definiteness, explicit support inclusion, common channel-output support,
or the active ensemble construction, so .toReal is its finite real value.
The accepted measurementInformation uses bits; multiplying by Real.log 2 gives
exactly the natural-logarithm entropy difference (47). -/

/-- Natural Hermitian logarithm (PDF pp.4–5). At a zero eigenvalue the value is
zero; this is used only inside x log x, or as zero extension of a common-support
logarithm difference. It does not assert that a singular matrix has a finite logarithm. -/
def logMatrix {A : Type} [Fintype A] [DecidableEq A] (X : Mat A) : Mat A :=
  cfc Real.log X

/-- Support inclusion supp X ⊆ supp Y, PDF p.4. -/
def SupportedBy {A : Type} [Fintype A] [DecidableEq A] (X Y : Mat A) : Prop :=
  LinearMap.ker Y.mulVecLin ≤ LinearMap.ker X.mulVecLin

/-- Default strictly positive normalized density convention, PDF pp.3,5. -/
def FaithfulDensity {A : Type} [Fintype A] (X : Mat A) : Prop :=
  X.PosDef ∧ X.trace = 1

/-- AB operator lifted by the unnormalized identity on C (PDF p.4). -/
def liftAB {A B C : Type} [DecidableEq C] (X : Mat (A × B)) : Mat (A × (B × C)) :=
  fun p q => if p.2.2 = q.2.2 then X (p.1, p.2.1) (q.1, q.2.1) else 0

/-- BC operator lifted by the identity on A (PDF p.4). -/
def liftBC {A B C : Type} [DecidableEq A] (X : Mat (B × C)) : Mat (A × (B × C)) :=
  fun p q => if p.1 = q.1 then X p.2 q.2 else 0

/-- B operator lifted by identities on A and C (PDF p.4). -/
def liftB {A B C : Type} [DecidableEq A] [DecidableEq C] (X : Mat B) : Mat (A × (B × C)) :=
  fun p q => if p.1 = q.1 ∧ p.2.2 = q.2.2 then X p.2.1 q.2.1 else 0

/-- Right minus left side of SSA (3), PDF p.3, with actual marginals. -/
def ssaGap {A B C : Type} [Fintype A] [Fintype B] [Fintype C]
    [DecidableEq A] [DecidableEq B] [DecidableEq C] (X : Mat (A × (B × C))) : ℝ :=
  entropy (marginalAB X) + entropy (traceLeft X) - entropy X -
    entropy (traceRight (traceLeft X))

/-- The literal lifted logarithm equality (4), PDF p.4. -/
def SSAlogCondition {A B C : Type} [Fintype A] [Fintype B] [Fintype C]
    [DecidableEq A] [DecidableEq B] [DecidableEq C] (X : Mat (A × (B × C))) : Prop :=
  logMatrix X - liftAB (logMatrix (marginalAB X)) =
    liftBC (logMatrix (traceLeft X)) - liftB (logMatrix (traceRight (traceLeft X)))

/-- exp(itL), with the natural-log imaginary-power convention, PDF p.7. -/
def imaginaryExp {A : Type} [Fintype A] [DecidableEq A] (t : ℝ) (L : Mat A) : Mat A :=
  NormedSpace.exp (((t : ℂ) * Complex.I) • L)

/-- Regrouping of rho_(1,2a) tensor rho_(2b,3), PDF p.7. -/
def splitProduct {A B D C : Type} (X : Mat (A × B)) (Y : Mat (D × C)) :
    Mat (A × ((B × D) × C)) :=
  fun p q => X (p.1, p.2.1.1) (q.1, q.2.1.1) * Y (p.2.1.2, p.2.2) (q.2.1.2, q.2.2)

/-- ABC marginal of the four-register purification, PDF p.9. -/
def marginalABC {A B C D : Type} [Fintype D] (X : Mat (A × (B × (C × D)))) :
    Mat (A × (B × C)) :=
  fun p q => ∑ d, X (p.1, (p.2.1, (p.2.2, d))) (q.1, (q.2.1, (q.2.2, d)))

/-- ABD marginal of the same four-register purification, PDF p.9. -/
def marginalABD {A B C D : Type} [Fintype C] (X : Mat (A × (B × (C × D)))) :
    Mat (A × (B × D)) :=
  fun p q => ∑ c, X (p.1, (p.2.1, (c, p.2.2))) (q.1, (q.2.1, (c, q.2.2)))

/-- Inequality (22) on the actual 124 marginal (PDF pp.9,15). -/
def WeakMonotonicity {A B D : Type} [Fintype A] [Fintype B] [Fintype D]
    [DecidableEq A] [DecidableEq B] [DecidableEq D] (X : Mat (A × (B × D))) : Prop :=
  entropy (traceLeft (traceLeft X)) + entropy (traceRight (traceLeft X)) ≤
    entropy (marginalAB X) + entropy (marginalAC X)

/-- The computed rectangular Kraus isometry in (26), PDF pp.10–11. -/
def krausV {A B K : Type} (F : K → Matrix B A ℂ) : Matrix (B × K) A ℂ :=
  fun bk a => F bk.2 bk.1 a

/-- Coordinate projections I tensor |b><b| in Section 3.5, PDF p.11. -/
def coordinatePVM {A O : Type} [DecidableEq A] [DecidableEq O] (b : O) : Mat (A × O) :=
  Matrix.kronecker (1 : Mat A) (Matrix.single b b 1)

/-- Common-basis rank-one spectral measurement, PDF p.17. -/
def basisPVM {A : Type} (U : Mat A) (a : A) : Mat A :=
  fun i j => U i a * conj (U j a)

/-- Inverse square root for the faithful ensemble–POVM duality, PDF p.21. -/
def inverseRoot {A : Type} [Fintype A] [DecidableEq A] (rho : Mat A) : Mat A :=
  cfc (fun x : ℝ => (Real.sqrt x)⁻¹) rho

/-- Dual POVM, PDF p.21. Its uses require a faithful reference and positive weights. -/
def dualEffects {J A : Type} [Fintype A] [DecidableEq A]
    (rho : Mat A) (w : J → ℝ) (sigma : J → Mat A) : J → Mat A :=
  fun j => (w j : ℂ) • (inverseRoot rho * sigma j * inverseRoot rho)

/-- Unnormalized dual ensemble component sqrt(rho) E_b sqrt(rho), PDF p.21. -/
def dualUnnormalized {O A : Type} [Fintype A] [DecidableEq A]
    (rho : Mat A) (E : O → Mat A) (b : O) : Mat A :=
  cfc Real.sqrt rho * E b * cfc Real.sqrt rho

/-- Dual normalized state at an active outcome tau_b>0, PDF p.21. Its zero-mass
value is zero and is never treated as a density matrix or placed inside a log. -/
def dualStates {O A : Type} [Fintype A] [DecidableEq A]
    (rho : Mat A) (E : O → Mat A) (b : O) : Mat A :=
  ((effectProb (E b) rho : ℂ)⁻¹) • dualUnnormalized rho E b

/-- Hall's right side (59), PDF p.21, explicitly omitting zero-mass outcomes. -/
def hallBound {O A : Type} [Fintype O] [Fintype A] [DecidableEq A]
    (rho : Mat A) (E : O → Mat A) : ℝ :=
  entropy rho - ∑ b, if 0 < effectProb (E b) rho then
    effectProb (E b) rho * entropy (dualStates rho E b) else 0

/-- Equation (60), PDF p.21, retaining the factor tau_b. -/
def hallRelative {O A : Type} [Fintype O] [Fintype A] [DecidableEq A]
    (rho : Mat A) (E : O → Mat A) : ℝ :=
  ∑ b, if 0 < effectProb (E b) rho then
    effectProb (E b) rho * (relativeEntropy (dualStates rho E b) rho).toReal else 0

/-- The optimizing expression (61), PDF p.21. Only the definition is retained;
no strict comparison between global optima is inferred from one maximizing pair. -/
def UEP {A B : Type} [Fintype A] [Fintype B] [DecidableEq A] [DecidableEq B]
    (Phi : Superop A B) (_hPhi : IsChannel Phi) [Nonempty A] [Nonempty B] : ℝ :=
  sSup {v : ℝ | ∃ rho : Mat A, IsDensityMatrix rho ∧
    ∃ m : ℕ, 0 < m ∧ ∃ E : Fin m → Mat B, IsPOVM E ∧
      v = hallBound rho (fun b => channelAdjoint Phi (E b))}

/-- MONO as a global finite PSD/support principle (PDF pp.14–15). Unnormalized
positive matrices are included, so embedded references used in the reductions are covered. -/
def MONO : Prop :=
  ∀ n m : ℕ, ∀ Phi : Superop (Fin (n + 1)) (Fin (m + 1)), IsChannel Phi →
    ∀ X Y : Mat (Fin (n + 1)), X.PosSemidef → Y.PosSemidef → SupportedBy X Y →
      SupportedBy (Phi X) (Phi Y) ∧ (relativeEntropy (Phi X) (Phi Y)).toReal ≤ (relativeEntropy X Y).toReal

/-- MPT as a global principle, with genuine partial traces, PDF pp.13–15. -/
def MPT : Prop :=
  ∀ n m : ℕ, ∀ X Y : Mat (Fin (n + 1) × Fin (m + 1)),
    X.PosSemidef → Y.PosSemidef → SupportedBy X Y →
      SupportedBy (traceLeft X) (traceLeft Y) ∧
      (relativeEntropy (traceLeft X) (traceLeft Y)).toReal ≤ (relativeEntropy X Y).toReal

/-- SSA as a global normalized PSD inequality, including block and pure states, PDF p.15. -/
def SSA : Prop :=
  ∀ n m k : ℕ, ∀ X : Mat (Fin (n + 1) × (Fin (m + 1) × Fin (k + 1))),
    IsDensityMatrix X → 0 ≤ ssaGap X

/-- JC as a global finite PSD/support principle, PDF pp.14–15. -/
def JC : Prop :=
  ∀ n k : ℕ, ∀ w : Fin (k + 1) → ℝ, IsLaw w →
    ∀ X Y : Fin (k + 1) → Mat (Fin (n + 1)),
      (∀ j, (X j).PosSemidef ∧ (Y j).PosSemidef ∧ SupportedBy (X j) (Y j)) →
      SupportedBy (average w X) (average w Y) ∧
      (relativeEntropy (average w X) (average w Y)).toReal ≤ ∑ j, w j * (relativeEntropy (X j) (Y j)).toReal

/-- Admissible positively oriented circle for the principal analytic matrix logarithm,
equation (63), PDF p.24. Its closed disk is in the scalar upper half-plane and
its open disk contains the entire spectrum, ensuring winding number one and no resolvent pole. -/
def UpperLogContour {A : Type} [Fintype A] [DecidableEq A]
    (M : Mat A) (c : ℂ) (r : ℝ) : Prop :=
  0 < r ∧ r < c.im ∧ spectrum ℂ M ⊆ Metric.ball c r

/-- Principal holomorphic matrix logarithm on upper-half-plane spectrum, using
exactly the contour-resolvent calculus (63), PDF p.24. Different admissible circles
give the same value by holomorphic functional calculus. Outside this domain the
zero value is a totalization only, never used by the Pick theorem. -/
def principalLog {A : Type} [Fintype A] [DecidableEq A] (M : Mat A) : Mat A :=
  if h : ∃ cr : ℂ × ℝ, UpperLogContour M cr.1 cr.2 then
    let cr := Classical.choose h
    ((2 * (Real.pi : ℂ) * Complex.I)⁻¹) •
      circleIntegral (fun z => Complex.log z • (z • (1 : Mat A) - M)⁻¹) cr.1 cr.2
  else 0

/-- Epstein's nonnormal matrix pencil trace function (62), PDF p.24. -/
def epsteinTrace {A : Type} [Fintype A] [DecidableEq A]
    (K R B : Mat A) (z : ℂ) : ℂ :=
  (NormedSpace.exp (K + principalLog (z • R + B))).trace

/-! Sanity checks pinning the local definitions, without proving the paper's theorems. -/

example {A : Type} [Fintype A] [DecidableEq A] : entropy (1 : Mat A) = 0 := by
  simp [entropy, cfc_one]

example {A : Type} [Fintype A] [DecidableEq A] : logMatrix (1 : Mat A) = 0 := by
  simp [logMatrix, cfc_one]

example : shannonNats (fun _ : Unit => (1 : ℝ)) = 0 := by
  simp [shannonNats]

example {A B K : Type} (F : K → Matrix B A ℂ) (b : B) (k : K) (a : A) :
    krausV F (b, k) a = F k b a := rfl

example {A B C : Type} [Fintype C] (X : Mat (A × (B × C))) (a a2 : A) (b b2 : B) :
    marginalAB X (a, b) (a2, b2) = ∑ c, X (a, (b, c)) (a2, (b2, c)) := rfl

example {A : Type} [Fintype A] [DecidableEq A] (X : Mat A)
    (h : ¬ ∃ cr : ℂ × ℝ, UpperLogContour X cr.1 cr.2) : principalLog X = 0 := by
  simp [principalLog, h]

example {A : Type} [Fintype A] [DecidableEq A] : (relativeEntropy (1 : Mat A) 1).toReal = 0 := by
  simp [relativeEntropy, QAlgorithms.Papers.Independent.Dependency50.P21.IsSupported, cfc_one]

example {A : Type} [Fintype A] [DecidableEq A] (X : Mat A) : SupportedBy X X := le_refl _

example : FaithfulDensity (1 : Mat (Fin 1)) := by
  exact ⟨Matrix.PosDef.one, by simp [Matrix.trace]⟩

example : IsLaw (fun _ : Unit => (1 : ℝ)) := by
  exact ⟨fun _ => by norm_num, by simp⟩

example {A : Type} (X : Mat A) : average (fun _ : Unit => (1 : ℝ)) (fun _ => X) = X := by
  simp [QAlgorithms.Papers.Independent.Information100.P27.average]

example {A B C : Type} [DecidableEq C] (X : Mat (A × B)) (a a2 : A) (b b2 : B) (c : C) :
    liftAB X (a, (b, c)) (a2, (b2, c)) = X (a, b) (a2, b2) := by
  simp [liftAB]

example {A B C : Type} [DecidableEq A] (X : Mat (B × C)) (a : A) (b b2 : B) (c c2 : C) :
    liftBC X (a, (b, c)) (a, (b2, c2)) = X (b, c) (b2, c2) := by
  simp [liftBC]

example {A B C : Type} [DecidableEq A] [DecidableEq C] (X : Mat B) (a : A) (b b2 : B) (c : C) :
    liftB X (a, (b, c)) (a, (b2, c)) = X b b2 := by
  simp [liftB]

example {A : Type} [Fintype A] [DecidableEq A] (L : Mat A) : imaginaryExp 0 L = 1 := by
  simp [imaginaryExp]

example {A B D C : Type} (X : Mat (A × B)) (Y : Mat (D × C))
    (a a2 : A) (b b2 : B) (d d2 : D) (c c2 : C) :
    splitProduct X Y (a, ((b, d), c)) (a2, ((b2, d2), c2)) =
      X (a, b) (a2, b2) * Y (d, c) (d2, c2) := rfl

example {A : Type} [Fintype A] [DecidableEq A] (X : Mat A) :
    holevo (fun _ : Unit => (1 : ℝ)) (fun _ => X) = 0 := by
  simp [holevo, QAlgorithms.Papers.Independent.Information100.P56.average]

example {A : Type} [Fintype A] [DecidableEq A] : inverseRoot (1 : Mat A) = 1 := by
  simp [inverseRoot, cfc_one]

example {O A : Type} [Fintype A] [DecidableEq A] (E : O → Mat A) (b : O) :
    dualUnnormalized (1 : Mat A) E b = E b := by
  simp [dualUnnormalized, cfc_one]

example {A : Type} [DecidableEq A] (a : A) : basisPVM (1 : Mat A) a a a = 1 := by
  simp [basisPVM]

section Statements

variable {A B C D J K O : Type}
variable [Fintype A] [Fintype B] [Fintype C] [Fintype D] [Fintype J] [Fintype K] [Fintype O]
variable [DecidableEq A] [DecidableEq B] [DecidableEq C] [DecidableEq D]
variable [DecidableEq J] [DecidableEq K] [DecidableEq O]

end Statements

end QAlgorithms.Papers.Independent.Intakethird100.P51
Source
Inequalities for Quantum Entropy: A Review with Conditions for Equality, arXiv:quant-ph/0205064v2, https://arxiv.org/abs/quant-ph/0205064v2
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P51.FaithfulDensity

For any finite type AAA and any complex matrix XXX indexed by A×AA\times AA×A, this defines a faithful density matrix to mean that XXX is positive definite and its complex trace is exactly 111. No nonemptiness assumption is imposed on AAA; when AAA is empty, its trace is 000, so this predicate cannot hold.

QAlgorithms.Papers.Independent.Intakethird100.P51.JC

This proposition says that for every n,k∈Nn,k\in\mathbb Nn,k∈N, every family of real weights w0,…,wkw_0,\ldots,w_kw0​,…,wk​ with wj≥0w_j\ge0wj​≥0 and ∑jwj=1\sum_jw_j=1∑j​wj​=1, and every pair of families of positive semidefinite complex (n+1)×(n+1)(n+1)\times(n+1)(n+1)×(n+1) matrices Xj,YjX_j,Y_jXj​,Yj​ satisfying ker⁡Yj⊆ker⁡Xj\ker Y_j\subseteq\ker X_jkerYj​⊆kerXj​ for every jjj, the matrices X‾=∑jwjXj\overline X=\sum_jw_jX_jX=∑j​wj​Xj​ and Y‾=∑jwjYj\overline Y=\sum_jw_jY_jY=∑j​wj​Yj​ satisfy

ker⁡Y‾⊆ker⁡X‾,D(X‾∥Y‾)≤∑jwjD(Xj∥Yj).\ker\overline Y\subseteq\ker\overline X, \qquad D(\overline X\|\overline Y)\le\sum_jw_jD(X_j\|Y_j).kerY⊆kerX,D(X∥Y)≤j∑​wj​D(Xj​∥Yj​).

For a supported pair, D(P∥Q)=Re⁡tr⁡(f(P)−Plog⁡Q)D(P\|Q)=\operatorname{Re}\operatorname{tr}(f(P)-P\log Q)D(P∥Q)=Retr(f(P)−PlogQ), where f(x)=xln⁡xf(x)=x\ln xf(x)=xlnx and the logarithm are applied by continuous functional calculus, with their total real values at zero. The relative entropy definition assigns +∞+\infty+∞ outside the support condition, and the displayed inequality uses its real conversion, which sends infinities to zero; every pair used here is required or asserted to be supported. The matrices need not be normalized or nonzero. Zero weights are allowed, but their corresponding matrices must still satisfy all hypotheses; n=0n=0n=0 and k=0k=0k=0 are included.

QAlgorithms.Papers.Independent.Intakethird100.P51.MONO

This defines the following proposition: for every pair of natural numbers n,mn,mn,m, every channel Φ\PhiΦ from complex (n+1)×(n+1)(n+1)\times(n+1)(n+1)×(n+1) matrices to complex (m+1)×(m+1)(m+1)\times(m+1)(m+1)×(m+1) matrices, and every positive semidefinite pair X,YX,YX,Y of input matrices satisfying ker⁡Y⊆ker⁡X\ker Y\subseteq\ker XkerY⊆kerX, one has

ker⁡Φ(Y)⊆ker⁡Φ(X)andR(Φ(X),Φ(Y))≤R(X,Y).\ker\Phi(Y)\subseteq\ker\Phi(X) \quad\text{and}\quad \mathcal R(\Phi(X),\Phi(Y))\le\mathcal R(X,Y).kerΦ(Y)⊆kerΦ(X)andR(Φ(X),Φ(Y))≤R(X,Y).

Here a channel means a map admitting a finite representation Φ(Z)=∑jKjZKj†\Phi(Z)=\sum_jK_jZK_j^\daggerΦ(Z)=∑j​Kj​ZKj†​ with ∑jKj†Kj=In+1\sum_jK_j^\dagger K_j=I_{n+1}∑j​Kj†​Kj​=In+1​, and

R(P,Q)={Re⁡tr⁡ ⁣(f(P)−Pℓ(Q)),ker⁡Q⊆ker⁡P,0,otherwise,f(x)=xln⁡x,ℓ(x)=ln⁡x,\mathcal R(P,Q)= \begin{cases} \operatorname{Re}\operatorname{tr}\!\bigl(f(P)-P\ell(Q)\bigr),&\ker Q\subseteq\ker P,\\ 0,&\text{otherwise}, \end{cases} \qquad f(x)=x\ln x,\quad\ell(x)=\ln x,R(P,Q)={Retr(f(P)−Pℓ(Q)),0,​kerQ⊆kerP,otherwise,​f(x)=xlnx,ℓ(x)=lnx,

with matrix functions given by real continuous functional calculus. This is the real conversion of an extended-real relative entropy whose unsupported value is +∞+\infty+∞, converted to zero. Neither XXX nor YYY must have trace one or be nonzero. The cases n=0n=0n=0 and m=0m=0m=0 give one-dimensional registers, and zero matrices are included.

QAlgorithms.Papers.Independent.Intakethird100.P51.MPT

This proposition says that for every pair of natural numbers n,mn,mn,m, including zero, and all positive semidefinite complex matrices X,YX,YX,Y indexed by {0,…,n}×{0,…,m}\{0,\ldots,n\}\times\{0,\ldots,m\}{0,…,n}×{0,…,m} satisfying ker⁡Y⊆ker⁡X\ker Y\subseteq\ker XkerY⊆kerX, their partial traces over the first factor satisfy

ker⁡(Tr⁡1Y)⊆ker⁡(Tr⁡1X),D(Tr⁡1X∥Tr⁡1Y)≤D(X∥Y).\ker(\operatorname{Tr}_1Y)\subseteq\ker(\operatorname{Tr}_1X), \qquad D(\operatorname{Tr}_1X\|\operatorname{Tr}_1Y)\le D(X\|Y).ker(Tr1​Y)⊆ker(Tr1​X),D(Tr1​X∥Tr1​Y)≤D(X∥Y).

Here (Tr⁡1X)bb′=∑aX(a,b),(a,b′)(\operatorname{Tr}_1X)_{bb'}=\sum_aX_{(a,b),(a,b')}(Tr1​X)bb′​=∑a​X(a,b),(a,b′)​, and, for a supported pair P,QP,QP,Q,

D(P∥Q)=Re⁡tr⁡ ⁣(f(P)−Plog⁡Q),f(x)=xln⁡x,D(P\|Q)=\operatorname{Re}\operatorname{tr}\!\left(f(P)-P\log Q\right), \qquad f(x)=x\ln x,D(P∥Q)=Retr(f(P)−PlogQ),f(x)=xlnx,

with both matrix functions defined by continuous functional calculus using the total real logarithm, for which ln⁡0=0\ln0=0ln0=0. The underlying relative entropy is extended-real-valued and equals +∞+\infty+∞ when the kernel inclusion fails; the statement compares its real conversions, which send either infinity to zero, although the stated support conditions make the quantities here finite. Neither XXX nor YYY is required to have trace one or to be nonzero; the two register dimensions are always positive.

QAlgorithms.Papers.Independent.Intakethird100.P51.SSA

This proposition says that for every n,m,k∈Nn,m,k\in\mathbb Nn,m,k∈N, including zero, and every positive semidefinite complex matrix XXX of trace one on registers A={0,…,n}A=\{0,\ldots,n\}A={0,…,n}, B={0,…,m}B=\{0,\ldots,m\}B={0,…,m}, and C={0,…,k}C=\{0,\ldots,k\}C={0,…,k},

0≤S(XAB)+S(XBC)−S(X)−S(XB).0\le S(X_{AB})+S(X_{BC})-S(X)-S(X_B).0≤S(XAB​)+S(XBC​)−S(X)−S(XB​).

The marginals are defined entrywise by

(XAB)(a,b),(a′,b′)=∑cX(a,b,c),(a′,b′,c),(XBC)(b,c),(b′,c′)=∑aX(a,b,c),(a,b′,c′),(XB)b,b′=∑a,cX(a,b,c),(a,b′,c).(X_{AB})_{(a,b),(a',b')}=\sum_cX_{(a,b,c),(a',b',c)},\quad (X_{BC})_{(b,c),(b',c')}=\sum_aX_{(a,b,c),(a,b',c')},\quad (X_B)_{b,b'}=\sum_{a,c}X_{(a,b,c),(a,b',c)}.(XAB​)(a,b),(a′,b′)​=c∑​X(a,b,c),(a′,b′,c)​,(XBC​)(b,c),(b′,c′)​=a∑​X(a,b,c),(a,b′,c′)​,(XB​)b,b′​=a,c∑​X(a,b,c),(a,b′,c)​.

Here S(Z)=−Re⁡tr⁡f(Z)S(Z)=-\operatorname{Re}\operatorname{tr}f(Z)S(Z)=−Retrf(Z), where f(x)=xln⁡xf(x)=x\ln xf(x)=xlnx is applied through continuous functional calculus and f(0)=0f(0)=0f(0)=0. Thus the entropy uses natural logarithms. Singular density matrices and one-dimensional registers are included.

QAlgorithms.Papers.Independent.Intakethird100.P51.SSAlogCondition

For finite types A,B,CA,B,CA,B,C, each with decidable equality, and any complex matrix XXX on A×(B×C)A\times(B\times C)A×(B×C), let XABX_{AB}XAB​, XBCX_{BC}XBC​, and XBX_BXB​ have entries

(XAB)(a,b),(a′,b′)=∑cX(a,(b,c)),(a′,(b′,c)),(XBC)(b,c),(b′,c′)=∑aX(a,(b,c)),(a,(b′,c′)),(XB)b,b′=∑c(XBC)(b,c),(b′,c).(X_{AB})_{(a,b),(a',b')}=\sum_cX_{(a,(b,c)),(a',(b',c))},\qquad (X_{BC})_{(b,c),(b',c')}=\sum_aX_{(a,(b,c)),(a,(b',c'))},\qquad (X_B)_{b,b'}=\sum_c(X_{BC})_{(b,c),(b',c)}.(XAB​)(a,b),(a′,b′)​=c∑​X(a,(b,c)),(a′,(b′,c))​,(XBC​)(b,c),(b′,c′)​=a∑​X(a,(b,c)),(a,(b′,c′))​,(XB​)b,b′​=c∑​(XBC​)(b,c),(b′,c)​.

Writing L(M)L(M)L(M) for the real natural logarithm applied through continuous functional calculus—using log⁡0=0\log0=0log0=0, log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣ for negative real xxx, and the zero matrix for non-Hermitian MMM—this predicate means precisely that, for every a,a′∈Aa,a'\in Aa,a′∈A, b,b′∈Bb,b'\in Bb,b′∈B, and c,c′∈Cc,c'\in Cc,c′∈C,

L(X)(a,(b,c)),(a′,(b′,c′))−1c=c′L(XAB)(a,b),(a′,b′)=1a=a′L(XBC)(b,c),(b′,c′)−1a=a′1c=c′L(XB)b,b′.L(X)_{(a,(b,c)),(a',(b',c'))} -\mathbf1_{c=c'}L(X_{AB})_{(a,b),(a',b')} = \mathbf1_{a=a'}L(X_{BC})_{(b,c),(b',c')} -\mathbf1_{a=a'}\mathbf1_{c=c'}L(X_B)_{b,b'}.L(X)(a,(b,c)),(a′,(b′,c′))​−1c=c′​L(XAB​)(a,b),(a′,b′)​=1a=a′​L(XBC​)(b,c),(b′,c′)​−1a=a′​1c=c′​L(XB​)b,b′​.

Here 1P\mathbf1_P1P​ is 111 when PPP holds and 000 otherwise. No positivity, normalization, invertibility, or nonemptiness is assumed; if the full index type is empty, this matrix equality holds vacuously.

QAlgorithms.Papers.Independent.Intakethird100.P51.SupportedBy

For any finite type AAA with decidable equality and any complex matrices X,YX,YX,Y indexed by A×AA\times AA×A, this defines “XXX is supported by YYY” to mean ker⁡Y⊆ker⁡X\ker Y\subseteq\ker XkerY⊆kerX: every vector v:A→Cv:A\to\mathbb Cv:A→C satisfying Yv=0Yv=0Yv=0 also satisfies Xv=0Xv=0Xv=0. Neither matrix is assumed positive or Hermitian. If AAA is empty, the condition holds automatically.

QAlgorithms.Papers.Independent.Intakethird100.P51.UEP

For finite nonempty types A,BA,BA,B, both with decidable equality, and a supplied channel Φ\PhiΦ from complex matrices on AAA to complex matrices on BBB, this defines the real supremum of

H(ρ)−∑b=0,…,m−1pb>0pbH(Zb)H(\rho)-\sum_{\substack{b=0,\ldots,m-1\\p_b>0}}p_bH(Z_b)H(ρ)−b=0,…,m−1pb​>0​∑​pb​H(Zb​)

over every positive integer mmm, every positive semidefinite ρ\rhoρ with tr⁡ρ=1\operatorname{tr}\rho=1trρ=1, and every family of positive semidefinite matrices EbE_bEb​ on BBB satisfying ∑bEb=IB\sum_bE_b=I_B∑b​Eb​=IB​. In this expression,

Qb(i,j)=∑u,v∈BΦ(EijA)u,v‾ (Eb)u,v,pb=Re⁡tr⁡(Qbρ),Zb=pb−1SρQbSρ,Sρ=(⋅)(ρ),\begin{aligned} Q_b(i,j)&=\sum_{u,v\in B}\overline{\Phi(E^{A}_{ij})_{u,v}}\,(E_b)_{u,v},& p_b&=\operatorname{Re}\operatorname{tr}(Q_b\rho),\\ Z_b&=p_b^{-1}S_\rho Q_bS_\rho,& S_\rho&=(\sqrt{\cdot})(\rho), \end{aligned}Qb​(i,j)Zb​​=u,v∈B∑​Φ(EijA​)u,v​​(Eb​)u,v​,=pb−1​Sρ​Qb​Sρ​,​pb​Sρ​​=Retr(Qb​ρ),=(⋅​)(ρ),​

where EijAE^{A}_{ij}EijA​ is the matrix unit and H(M)=−Re⁡tr⁡(f(M))H(M)=-\operatorname{Re}\operatorname{tr}(f(M))H(M)=−Retr(f(M)), f(x)=xln⁡xf(x)=x\ln xf(x)=xlnx, using real continuous functional calculus. Zero reciprocals are zero, and terms with pb≤0p_b\le0pb​≤0 are omitted. The supplied channel hypothesis means that some finite family of rectangular matrices KjK_jKj​ satisfies ∑jKj†Kj=IA\sum_jK_j^\dagger K_j=I_A∑j​Kj†​Kj​=IA​ and Φ(X)=∑jKjXKj†\Phi(X)=\sum_jK_jXK_j^\daggerΦ(X)=∑j​Kj​XKj†​ for every XXX. The definition takes a supremum and does not itself assert attainment.

QAlgorithms.Papers.Independent.Intakethird100.P51.UpperLogContour

For any finite index type AAA with decidable equality, any complex matrix MMM indexed by AAA, any c∈Cc\in\mathbb Cc∈C, and any r∈Rr\in\mathbb Rr∈R, this predicate means precisely

0<r<Im⁡candσC(M)⊆{z∈C:∣z−c∣<r}.0<r<\operatorname{Im}c \quad\text{and}\quad \sigma_{\mathbb C}(M)\subseteq\{z\in\mathbb C:|z-c|<r\}.0<r<ImcandσC​(M)⊆{z∈C:∣z−c∣<r}.

Thus the entire complex spectrum must lie strictly inside an open disk whose closed boundary also lies in the open upper half-plane. No Hermitian or positivity assumption is imposed on MMM, and AAA is allowed to be empty, in which case its spectrum is empty and the spectral inclusion is vacuous.

QAlgorithms.Papers.Independent.Intakethird100.P51.WeakMonotonicity

For finite types A,B,DA,B,DA,B,D, each with decidable equality, and any complex matrix XXX indexed by A×(B×D)A\times(B\times D)A×(B×D), this defines the proposition

H(XD)+H(XB)≤H(XAB)+H(XAD),H(X_D)+H(X_B)\le H(X_{AB})+H(X_{AD}),H(XD​)+H(XB​)≤H(XAB​)+H(XAD​),

where

XD(d,d′)=∑a,bX(a,b,d),(a,b,d′),XB(b,b′)=∑a,dX(a,b,d),(a,b′,d),XAB((a,b),(a′,b′))=∑dX(a,b,d),(a′,b′,d),XAD((a,d),(a′,d′))=∑bX(a,b,d),(a′,b,d′).\begin{aligned} X_D(d,d')&=\sum_{a,b}X_{(a,b,d),(a,b,d')},& X_B(b,b')&=\sum_{a,d}X_{(a,b,d),(a,b',d)},\\ X_{AB}((a,b),(a',b'))&=\sum_d X_{(a,b,d),(a',b',d)},& X_{AD}((a,d),(a',d'))&=\sum_b X_{(a,b,d),(a',b,d')}. \end{aligned}XD​(d,d′)XAB​((a,b),(a′,b′))​=a,b∑​X(a,b,d),(a,b,d′)​,=d∑​X(a,b,d),(a′,b′,d)​,​XB​(b,b′)XAD​((a,d),(a′,d′))​=a,d∑​X(a,b,d),(a,b′,d)​,=b∑​X(a,b,d),(a′,b,d′)​.​

Here H(M)=−Re⁡tr⁡(f(M))H(M)=-\operatorname{Re}\operatorname{tr}(f(M))H(M)=−Retr(f(M)), with f(x)=xln⁡xf(x)=x\ln xf(x)=xlnx applied by real continuous functional calculus. No positivity or trace normalization of XXX is assumed. The types may be empty; if any is empty, all four entropy expressions are zero.

QAlgorithms.Papers.Independent.Intakethird100.P51.basisPVM

For any type AAA, any complex square matrix UUU indexed by AAA, and any a∈Aa\in Aa∈A, this defines the matrix

(Pa)i,j=Ui,a Uj,a‾(i,j∈A),(P_a)_{i,j}=U_{i,a}\,\overline{U_{j,a}}\qquad(i,j\in A),(Pa​)i,j​=Ui,a​Uj,a​​(i,j∈A),

the outer product of column aaa of UUU with itself. No finiteness, unitarity, orthogonality, or normalization assumption on UUU is imposed; a zero column gives the zero matrix.

QAlgorithms.Papers.Independent.Intakethird100.P51.coordinatePVM

For types A,OA,OA,O with decidable equality and a chosen b∈Ob\in Ob∈O, this defines the complex matrix Pb=IA⊗EbbP_b=I_A\otimes E_{bb}Pb​=IA​⊗Ebb​ indexed by A×OA\times OA×O, explicitly

(Pb)(a,o),(a′,o′)={1,a=a′, o=b, o′=b,0,otherwise.(P_b)_{(a,o),(a',o')}= \begin{cases} 1,&a=a',\ o=b,\ o'=b,\\ 0,&\text{otherwise}. \end{cases}(Pb​)(a,o),(a′,o′)​={1,0,​a=a′, o=b, o′=b,otherwise.​

Neither type is required to be finite; AAA may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.dualEffects

For any type JJJ, a finite type AAA with decidable equality, any complex matrices ρ\rhoρ and σj\sigma_jσj​ indexed by AAA, and any real weights wjw_jwj​, this defines the family

Dj=wjRρσjRρ,Rρ=f(ρ),f(x)=(x)−1,D_j=w_jR_\rho\sigma_jR_\rho,\qquad R_\rho=f(\rho),\quad f(x)=(\sqrt{x})^{-1},Dj​=wj​Rρ​σj​Rρ​,Rρ​=f(ρ),f(x)=(x​)−1,

where f(ρ)f(\rho)f(ρ) is real continuous functional calculus, with scalar value zero for x≤0x\le0x≤0 and with value zero on non-Hermitian matrices. No probability, density-matrix, or measurement conditions are required of w,ρ,σw,\rho,\sigmaw,ρ,σ. The type JJJ need not be finite, and JJJ or AAA may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.dualStates

For any type OOO, a finite type AAA with decidable equality, any complex matrix ρ\rhoρ indexed by AAA, any family EbE_bEb​ of such matrices, and b∈Ob\in Ob∈O, this defines

Zb=pb−1SρEbSρ,pb=Re⁡tr⁡(Ebρ),Sρ=(⋅)(ρ).Z_b=p_b^{-1}S_\rho E_bS_\rho,\qquad p_b=\operatorname{Re}\operatorname{tr}(E_b\rho),\qquad S_\rho=(\sqrt{\cdot})(\rho).Zb​=pb−1​Sρ​Eb​Sρ​,pb​=Retr(Eb​ρ),Sρ​=(⋅​)(ρ).

The real square root is applied by continuous functional calculus, which returns zero on non-Hermitian matrices. The reciprocal is total, so pb=0p_b=0pb​=0 gives Zb=0Z_b=0Zb​=0; negative pbp_bpb​ is also permitted. No density-matrix or measurement assumptions are imposed, OOO need not be finite, and AAA may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.dualUnnormalized

For any type OOO, a finite type AAA with decidable equality, any complex matrix ρ\rhoρ indexed by AAA, any family EbE_bEb​ of such matrices, and a chosen b∈Ob\in Ob∈O, this defines

Tb=SρEbSρ,Sρ=(⋅)(ρ),T_b=S_\rho E_bS_\rho,\qquad S_\rho=(\sqrt{\cdot})(\rho),Tb​=Sρ​Eb​Sρ​,Sρ​=(⋅​)(ρ),

where the square root is applied through real continuous functional calculus. Scalar negative inputs have square root zero, and this functional calculus returns zero on non-Hermitian matrices. No positivity, normalization, or measurement assumptions are imposed. The type OOO need not be finite, and AAA may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.epsteinTrace

For any finite index type AAA with decidable equality, arbitrary complex matrices K,R,BK,R,BK,R,B indexed by AAA, and any z∈Cz\in\mathbb Cz∈C, this definition returns the complex number

tr⁡exp⁡ ⁣(K+L(zR+B)).\operatorname{tr}\exp\!\bigl(K+\mathcal L(zR+B)\bigr).trexp(K+L(zR+B)).

Here, if a matrix MMM has its complex spectrum inside some open disk ∣ζ−c∣<r|\zeta-c|<r∣ζ−c∣<r with 0<r<Im⁡c0<r<\operatorname{Im}c0<r<Imc, then L(M)\mathcal L(M)L(M) is (2πi)−1∮∣ζ−c∣=rLog⁡(ζ)(ζI−M)−1 dζ(2\pi i)^{-1}\oint_{|\zeta-c|=r}\operatorname{Log}(\zeta)(\zeta I-M)^{-1}\,d\zeta(2πi)−1∮∣ζ−c∣=r​Log(ζ)(ζI−M)−1dζ, using one such disk selected by classical choice; otherwise L(M)=0\mathcal L(M)=0L(M)=0. The exponential is the matrix exponential. Consequently the returned value is tr⁡eK\operatorname{tr}e^KtreK whenever no qualifying disk exists for zR+BzR+BzR+B. No Hermitian, positivity, or nonemptiness hypothesis is imposed; for empty AAA, the trace is zero.

QAlgorithms.Papers.Independent.Intakethird100.P51.hallBound

For finite types O,AO,AO,A, with decidable equality on AAA, any complex matrix ρ\rhoρ indexed by AAA, and any family EbE_bEb​ of such matrices, this defines the real number

H(ρ)−∑b∈Opb>0pbH(Zb),pb=Re⁡tr⁡(Ebρ),Zb=pb−1SρEbSρ,Sρ=(⋅)(ρ).H(\rho)-\sum_{\substack{b\in O\\p_b>0}}p_bH(Z_b), \qquad p_b=\operatorname{Re}\operatorname{tr}(E_b\rho),\quad Z_b=p_b^{-1}S_\rho E_bS_\rho,\quad S_\rho=(\sqrt{\cdot})(\rho).H(ρ)−b∈Opb​>0​∑​pb​H(Zb​),pb​=Retr(Eb​ρ),Zb​=pb−1​Sρ​Eb​Sρ​,Sρ​=(⋅​)(ρ).

Here H(M)=−Re⁡tr⁡(f(M))H(M)=-\operatorname{Re}\operatorname{tr}(f(M))H(M)=−Retr(f(M)), f(x)=xln⁡xf(x)=x\ln xf(x)=xlnx, and both matrix functions use real continuous functional calculus, returning zero on non-Hermitian inputs. No positivity, normalization, or measurement assumptions are imposed. Terms with pb≤0p_b\le0pb​≤0 contribute zero; if OOO is empty, the value is H(ρ)H(\rho)H(ρ). The type AAA may also be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.hallRelative

For finite types O,AO,AO,A, with decidable equality on AAA, any complex matrix ρ\rhoρ indexed by AAA, and any family EbE_bEb​ of such matrices, this defines

∑b∈Opb>0pb R(Zb,ρ),pb=Re⁡tr⁡(Ebρ),Zb=pb−1SρEbSρ,Sρ=(⋅)(ρ).\sum_{\substack{b\in O\\p_b>0}}p_b\,\mathcal R(Z_b,\rho), \qquad p_b=\operatorname{Re}\operatorname{tr}(E_b\rho),\quad Z_b=p_b^{-1}S_\rho E_bS_\rho,\quad S_\rho=(\sqrt{\cdot})(\rho).b∈Opb​>0​∑​pb​R(Zb​,ρ),pb​=Retr(Eb​ρ),Zb​=pb−1​Sρ​Eb​Sρ​,Sρ​=(⋅​)(ρ).

For arbitrary matrices X,YX,YX,Y, the real number R(X,Y)\mathcal R(X,Y)R(X,Y) is

R(X,Y)={Re⁡tr⁡ ⁣(f(X)−Xℓ(Y)),ker⁡Y⊆ker⁡X,0,otherwise,f(x)=xln⁡x,ℓ(x)=ln⁡x.\mathcal R(X,Y)= \begin{cases} \operatorname{Re}\operatorname{tr}\!\bigl(f(X)-X\ell(Y)\bigr), &\ker Y\subseteq\ker X,\\ 0,&\text{otherwise}, \end{cases} \qquad f(x)=x\ln x,\quad \ell(x)=\ln x.R(X,Y)={Retr(f(X)−Xℓ(Y)),0,​kerY⊆kerX,otherwise,​f(x)=xlnx,ℓ(x)=lnx.

All matrix functions use real continuous functional calculus. The second branch arises because the underlying extended-real relative entropy is +∞+\infty+∞ when the kernel inclusion fails, and its conversion to a real number is zero. No density-matrix or measurement assumptions are imposed; terms with pb≤0p_b\le0pb​≤0 are zero, empty OOO gives zero, and AAA may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.imaginaryExp

For any finite type AAA with decidable equality, any real number ttt, and any complex matrix LLL indexed by A×AA\times AA×A, this defines

imaginaryExp⁡(t,L)=exp⁡(itL),\operatorname{imaginaryExp}(t,L)=\exp(itL),imaginaryExp(t,L)=exp(itL),

where exp⁡\expexp is the matrix exponential. The matrix LLL need not be Hermitian, so unitarity is not part of the definition. The case t=0t=0t=0 gives the identity matrix, and the empty index type is also allowed.

QAlgorithms.Papers.Independent.Intakethird100.P51.inverseRoot

For a finite type AAA with decidable equality and any complex square matrix ρ\rhoρ indexed by AAA, this defines

Rρ=f(ρ),f(x)=(x)−1,R_\rho=f(\rho),\qquad f(x)=(\sqrt{x})^{-1},Rρ​=f(ρ),f(x)=(x​)−1,

using real continuous functional calculus. The scalar square root and inverse are total: f(x)=0f(x)=0f(x)=0 for x≤0x\le0x≤0, including f(0)=0f(0)=0f(0)=0. Thus zero eigenvalues are assigned zero, rather than an infinite inverse. No positivity, invertibility, or trace normalization is assumed; the real functional calculus returns zero for a non-Hermitian matrix. The index type may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.krausV

For arbitrary types A,B,KA,B,KA,B,K and any family FkF_kFk​ of complex matrices with rows indexed by BBB and columns indexed by AAA, this defines the complex matrix VVV, with rows indexed by B×KB\times KB×K and columns indexed by AAA, by

V(b,k),a=(Fk)b,a.V_{(b,k),a}=(F_k)_{b,a}.V(b,k),a​=(Fk​)b,a​.

None of the types is required to be finite or nonempty, and no normalization, positivity, or completeness condition on the family is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P51.liftAB

For arbitrary types A,B,CA,B,CA,B,C, with decidable equality on CCC, and any complex matrix XXX indexed by (A×B)×(A×B)(A\times B)\times(A\times B)(A×B)×(A×B), this defines a matrix on A×(B×C)A\times(B\times C)A×(B×C) by

lift⁡AB(X)(a,(b,c)),(a′,(b′,c′))={X(a,b),(a′,b′),c=c′,0,c≠c′.\operatorname{lift}_{AB}(X)_{(a,(b,c)),(a',(b',c'))} = \begin{cases} X_{(a,b),(a',b')},&c=c',\\ 0,&c\ne c'. \end{cases}liftAB​(X)(a,(b,c)),(a′,(b′,c′))​={X(a,b),(a′,b′)​,0,​c=c′,c=c′.​

No finiteness, positivity, or nonemptiness is required of these types or matrices.

QAlgorithms.Papers.Independent.Intakethird100.P51.liftB

For arbitrary types A,B,CA,B,CA,B,C, with decidable equality on AAA and CCC, and any complex matrix XXX indexed by B×BB\times BB×B, this defines a matrix on A×(B×C)A\times(B\times C)A×(B×C) by

lift⁡B(X)(a,(b,c)),(a′,(b′,c′))={Xb,b′,a=a′ and c=c′,0,otherwise.\operatorname{lift}_{B}(X)_{(a,(b,c)),(a',(b',c'))} = \begin{cases} X_{b,b'},&a=a'\ \text{and}\ c=c',\\ 0,&\text{otherwise}. \end{cases}liftB​(X)(a,(b,c)),(a′,(b′,c′))​={Xb,b′​,0,​a=a′ and c=c′,otherwise.​

No finiteness, positivity, or nonemptiness is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P51.liftBC

For arbitrary types A,B,CA,B,CA,B,C, with decidable equality on AAA, and any complex matrix XXX indexed by (B×C)×(B×C)(B\times C)\times(B\times C)(B×C)×(B×C), this defines a matrix on A×(B×C)A\times(B\times C)A×(B×C) by

lift⁡BC(X)(a,(b,c)),(a′,(b′,c′))={X(b,c),(b′,c′),a=a′,0,a≠a′.\operatorname{lift}_{BC}(X)_{(a,(b,c)),(a',(b',c'))} = \begin{cases} X_{(b,c),(b',c')},&a=a',\\ 0,&a\ne a'. \end{cases}liftBC​(X)(a,(b,c)),(a′,(b′,c′))​={X(b,c),(b′,c′)​,0,​a=a′,a=a′.​

No finiteness, positivity, or nonemptiness is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P51.logMatrix

For any finite type AAA with decidable equality and any complex matrix XXX indexed by A×AA\times AA×A, this defines \logMatrix(X)\logMatrix(X)\logMatrix(X) as the continuous functional calculus of the real natural logarithm applied to XXX. For Hermitian XXX, this applies log⁡∣λ∣\log|\lambda|log∣λ∣ to each nonzero eigenvalue λ\lambdaλ and assigns 000 to a zero eigenvalue; for non-Hermitian XXX, the functional calculus returns the zero matrix. No positivity, invertibility, or normalization is assumed, and AAA may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P51.marginalABC

For arbitrary types A,B,C,DA,B,C,DA,B,C,D, with DDD finite, and any complex matrix XXX on A×(B×(C×D))A\times(B\times(C\times D))A×(B×(C×D)), this defines a matrix on A×(B×C)A\times(B\times C)A×(B×C) by

marginal⁡ABC(X)(a,(b,c)),(a′,(b′,c′))=∑d∈DX(a,(b,(c,d))),(a′,(b′,(c′,d))).\operatorname{marginal}_{ABC}(X)_{(a,(b,c)),(a',(b',c'))} = \sum_{d\in D}X_{(a,(b,(c,d))),(a',(b',(c',d)))}.marginalABC​(X)(a,(b,c)),(a′,(b′,c′))​=d∈D∑​X(a,(b,(c,d))),(a′,(b′,(c′,d)))​.

No density or positivity condition is assumed, and A,B,CA,B,CA,B,C need not be finite. If DDD is empty, the result is the zero matrix.

QAlgorithms.Papers.Independent.Intakethird100.P51.marginalABD

For arbitrary types A,B,C,DA,B,C,DA,B,C,D, with CCC finite, and any complex matrix XXX on A×(B×(C×D))A\times(B\times(C\times D))A×(B×(C×D)), this defines a matrix on A×(B×D)A\times(B\times D)A×(B×D) by

marginal⁡ABD(X)(a,(b,d)),(a′,(b′,d′))=∑c∈CX(a,(b,(c,d))),(a′,(b′,(c,d′))).\operatorname{marginal}_{ABD}(X)_{(a,(b,d)),(a',(b',d'))} = \sum_{c\in C}X_{(a,(b,(c,d))),(a',(b',(c,d')))}.marginalABD​(X)(a,(b,d)),(a′,(b′,d′))​=c∈C∑​X(a,(b,(c,d))),(a′,(b′,(c,d′)))​.

No density or positivity condition is assumed, and A,B,DA,B,DA,B,D need not be finite. If CCC is empty, the result is the zero matrix.

QAlgorithms.Papers.Independent.Intakethird100.P51.principalLog

For any finite index type AAA with decidable equality and any complex matrix MMM indexed by AAA, this definition first asks whether there exist c∈Cc\in\mathbb Cc∈C and r∈Rr\in\mathbb Rr∈R such that 0<r<Im⁡c0<r<\operatorname{Im}c0<r<Imc and σC(M)⊆{z:∣z−c∣<r}\sigma_{\mathbb C}(M)\subseteq\{z:|z-c|<r\}σC​(M)⊆{z:∣z−c∣<r}. If such pairs exist, it chooses one by classical choice and returns

12πi∮∣z−c∣=rLog⁡(z) (zI−M)−1 dz,\frac{1}{2\pi i}\oint_{|z-c|=r}\operatorname{Log}(z)\,(zI-M)^{-1}\,dz,2πi1​∮∣z−c∣=r​Log(z)(zI−M)−1dz,

where the circle has the standard positive orientation and Log⁡\operatorname{Log}Log is the principal scalar complex logarithm. If no such pair exists, it returns the zero matrix. The inverse in the integrand is the total matrix inverse, agreeing with the usual inverse on invertible matrices and returning zero on singular matrices; on the chosen circle, the spectral inclusion ensures invertibility. The definition allows empty AAA and imposes no positivity or Hermitian condition on MMM.

QAlgorithms.Papers.Independent.Intakethird100.P51.splitProduct

For arbitrary types A,B,D,CA,B,D,CA,B,D,C, any complex matrix XXX on A×BA\times BA×B, and any complex matrix YYY on D×CD\times CD×C, this defines a matrix on A×((B×D)×C)A\times((B\times D)\times C)A×((B×D)×C) whose entries are

splitProduct⁡(X,Y)(a,((b,d),c)),(a′,((b′,d′),c′))=X(a,b),(a′,b′) Y(d,c),(d′,c′).\operatorname{splitProduct}(X,Y)_{(a,((b,d),c)),(a',((b',d'),c'))} = X_{(a,b),(a',b')}\,Y_{(d,c),(d',c')}.splitProduct(X,Y)(a,((b,d),c)),(a′,((b′,d′),c′))​=X(a,b),(a′,b′)​Y(d,c),(d′,c′)​.

No finiteness, decidable equality, positivity, or nonemptiness assumptions are imposed.

QAlgorithms.Papers.Independent.Intakethird100.P51.ssaGap

For finite types A,B,CA,B,CA,B,C, each with decidable equality, and any complex matrix XXX on A×(B×C)A\times(B\times C)A×(B×C), define matrices

(XAB)(a,b),(a′,b′)=∑c∈CX(a,(b,c)),(a′,(b′,c)),(XBC)(b,c),(b′,c′)=∑a∈AX(a,(b,c)),(a,(b′,c′)),(XB)b,b′=∑bC∈C(XBC)(b,bC),(b′,bC).(X_{AB})_{(a,b),(a',b')}=\sum_{c\in C}X_{(a,(b,c)),(a',(b',c))},\qquad (X_{BC})_{(b,c),(b',c')}=\sum_{a\in A}X_{(a,(b,c)),(a,(b',c'))},\qquad (X_B)_{b,b'}=\sum_{b_C\in C}(X_{BC})_{(b,b_C),(b',b_C)}.(XAB​)(a,b),(a′,b′)​=c∈C∑​X(a,(b,c)),(a′,(b′,c))​,(XBC​)(b,c),(b′,c′)​=a∈A∑​X(a,(b,c)),(a,(b′,c′))​,(XB​)b,b′​=bC​∈C∑​(XBC​)(b,bC​),(b′,bC​)​.

This declaration defines the real number

ssaGap⁡(X)=S(XAB)+S(XBC)−S(X)−S(XB),S(M)=−Re⁡tr⁡ ⁣[f(M)],f(x)=xlog⁡x,\operatorname{ssaGap}(X)=S(X_{AB})+S(X_{BC})-S(X)-S(X_B), \qquad S(M)=-\operatorname{Re}\operatorname{tr}\!\bigl[f(M)\bigr],\quad f(x)=x\log x,ssaGap(X)=S(XAB​)+S(XBC​)−S(X)−S(XB​),S(M)=−Retr[f(M)],f(x)=xlogx,

where f(M)f(M)f(M) uses continuous functional calculus and log⁡\loglog is the real natural logarithm, with log⁡0=0\log0=0log0=0. No density, Hermiticity, or positivity condition is imposed, and no nonnegativity is asserted. The types may be empty, with empty sums equal to 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