Definitions for Inequalities for Quantum Entropy: A Review with Conditions for Equality
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P51The 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.
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
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 and any complex matrix indexed by , this defines a faithful density matrix to mean that is positive definite and its complex trace is exactly . No nonemptiness assumption is imposed on ; when is empty, its trace is , so this predicate cannot hold.
QAlgorithms.Papers.Independent.Intakethird100.P51.JC
This proposition says that for every , every family of real weights with and , and every pair of families of positive semidefinite complex matrices satisfying for every , the matrices and satisfy
For a supported pair, , where and the logarithm are applied by continuous functional calculus, with their total real values at zero. The relative entropy definition assigns 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; and are included.
QAlgorithms.Papers.Independent.Intakethird100.P51.MONO
This defines the following proposition: for every pair of natural numbers , every channel from complex matrices to complex matrices, and every positive semidefinite pair of input matrices satisfying , one has
Here a channel means a map admitting a finite representation with , and
with matrix functions given by real continuous functional calculus. This is the real conversion of an extended-real relative entropy whose unsupported value is , converted to zero. Neither nor must have trace one or be nonzero. The cases and 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 , including zero, and all positive semidefinite complex matrices indexed by satisfying , their partial traces over the first factor satisfy
Here , and, for a supported pair ,
with both matrix functions defined by continuous functional calculus using the total real logarithm, for which . The underlying relative entropy is extended-real-valued and equals 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 nor 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 , including zero, and every positive semidefinite complex matrix of trace one on registers , , and ,
The marginals are defined entrywise by
Here , where is applied through continuous functional calculus and . Thus the entropy uses natural logarithms. Singular density matrices and one-dimensional registers are included.
QAlgorithms.Papers.Independent.Intakethird100.P51.SSAlogCondition
For finite types , each with decidable equality, and any complex matrix on , let , , and have entries
Writing for the real natural logarithm applied through continuous functional calculus—using , for negative real , and the zero matrix for non-Hermitian —this predicate means precisely that, for every , , and ,
Here is when holds and 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 with decidable equality and any complex matrices indexed by , this defines “ is supported by ” to mean : every vector satisfying also satisfies . Neither matrix is assumed positive or Hermitian. If is empty, the condition holds automatically.
QAlgorithms.Papers.Independent.Intakethird100.P51.UEP
For finite nonempty types , both with decidable equality, and a supplied channel from complex matrices on to complex matrices on , this defines the real supremum of
over every positive integer , every positive semidefinite with , and every family of positive semidefinite matrices on satisfying . In this expression,
where is the matrix unit and , , using real continuous functional calculus. Zero reciprocals are zero, and terms with are omitted. The supplied channel hypothesis means that some finite family of rectangular matrices satisfies and for every . The definition takes a supremum and does not itself assert attainment.
QAlgorithms.Papers.Independent.Intakethird100.P51.UpperLogContour
For any finite index type with decidable equality, any complex matrix indexed by , any , and any , this predicate means precisely
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 , and 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 , each with decidable equality, and any complex matrix indexed by , this defines the proposition
where
Here , with applied by real continuous functional calculus. No positivity or trace normalization of 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 , any complex square matrix indexed by , and any , this defines the matrix
the outer product of column of with itself. No finiteness, unitarity, orthogonality, or normalization assumption on is imposed; a zero column gives the zero matrix.
QAlgorithms.Papers.Independent.Intakethird100.P51.coordinatePVM
For types with decidable equality and a chosen , this defines the complex matrix indexed by , explicitly
Neither type is required to be finite; may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.dualEffects
For any type , a finite type with decidable equality, any complex matrices and indexed by , and any real weights , this defines the family
where is real continuous functional calculus, with scalar value zero for and with value zero on non-Hermitian matrices. No probability, density-matrix, or measurement conditions are required of . The type need not be finite, and or may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.dualStates
For any type , a finite type with decidable equality, any complex matrix indexed by , any family of such matrices, and , this defines
The real square root is applied by continuous functional calculus, which returns zero on non-Hermitian matrices. The reciprocal is total, so gives ; negative is also permitted. No density-matrix or measurement assumptions are imposed, need not be finite, and may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.dualUnnormalized
For any type , a finite type with decidable equality, any complex matrix indexed by , any family of such matrices, and a chosen , this defines
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 need not be finite, and may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.epsteinTrace
For any finite index type with decidable equality, arbitrary complex matrices indexed by , and any , this definition returns the complex number
Here, if a matrix has its complex spectrum inside some open disk with , then is , using one such disk selected by classical choice; otherwise . The exponential is the matrix exponential. Consequently the returned value is whenever no qualifying disk exists for . No Hermitian, positivity, or nonemptiness hypothesis is imposed; for empty , the trace is zero.
QAlgorithms.Papers.Independent.Intakethird100.P51.hallBound
For finite types , with decidable equality on , any complex matrix indexed by , and any family of such matrices, this defines the real number
Here , , 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 contribute zero; if is empty, the value is . The type may also be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.hallRelative
For finite types , with decidable equality on , any complex matrix indexed by , and any family of such matrices, this defines
For arbitrary matrices , the real number is
All matrix functions use real continuous functional calculus. The second branch arises because the underlying extended-real relative entropy is when the kernel inclusion fails, and its conversion to a real number is zero. No density-matrix or measurement assumptions are imposed; terms with are zero, empty gives zero, and may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.imaginaryExp
For any finite type with decidable equality, any real number , and any complex matrix indexed by , this defines
where is the matrix exponential. The matrix need not be Hermitian, so unitarity is not part of the definition. The case gives the identity matrix, and the empty index type is also allowed.
QAlgorithms.Papers.Independent.Intakethird100.P51.inverseRoot
For a finite type with decidable equality and any complex square matrix indexed by , this defines
using real continuous functional calculus. The scalar square root and inverse are total: for , including . 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 and any family of complex matrices with rows indexed by and columns indexed by , this defines the complex matrix , with rows indexed by and columns indexed by , by
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 , with decidable equality on , and any complex matrix indexed by , this defines a matrix on by
No finiteness, positivity, or nonemptiness is required of these types or matrices.
QAlgorithms.Papers.Independent.Intakethird100.P51.liftB
For arbitrary types , with decidable equality on and , and any complex matrix indexed by , this defines a matrix on by
No finiteness, positivity, or nonemptiness is assumed.
QAlgorithms.Papers.Independent.Intakethird100.P51.liftBC
For arbitrary types , with decidable equality on , and any complex matrix indexed by , this defines a matrix on by
No finiteness, positivity, or nonemptiness is assumed.
QAlgorithms.Papers.Independent.Intakethird100.P51.logMatrix
For any finite type with decidable equality and any complex matrix indexed by , this defines as the continuous functional calculus of the real natural logarithm applied to . For Hermitian , this applies to each nonzero eigenvalue and assigns to a zero eigenvalue; for non-Hermitian , the functional calculus returns the zero matrix. No positivity, invertibility, or normalization is assumed, and may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P51.marginalABC
For arbitrary types , with finite, and any complex matrix on , this defines a matrix on by
No density or positivity condition is assumed, and need not be finite. If is empty, the result is the zero matrix.
QAlgorithms.Papers.Independent.Intakethird100.P51.marginalABD
For arbitrary types , with finite, and any complex matrix on , this defines a matrix on by
No density or positivity condition is assumed, and need not be finite. If is empty, the result is the zero matrix.
QAlgorithms.Papers.Independent.Intakethird100.P51.principalLog
For any finite index type with decidable equality and any complex matrix indexed by , this definition first asks whether there exist and such that and . If such pairs exist, it chooses one by classical choice and returns
where the circle has the standard positive orientation and 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 and imposes no positivity or Hermitian condition on .
QAlgorithms.Papers.Independent.Intakethird100.P51.splitProduct
For arbitrary types , any complex matrix on , and any complex matrix on , this defines a matrix on whose entries are
No finiteness, decidable equality, positivity, or nonemptiness assumptions are imposed.
QAlgorithms.Papers.Independent.Intakethird100.P51.ssaGap
For finite types , each with decidable equality, and any complex matrix on , define matrices
This declaration defines the real number
where uses continuous functional calculus and is the real natural logarithm, with . No density, Hermiticity, or positivity condition is imposed, and no nonnegativity is asserted. The types may be empty, with empty sums equal to zero.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.