Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Quantum Markov chains, sufficiency of quantum channels, and Renyi information measures

Definition
QAlgorithms_Papers_Independent_Intakethird100_P52

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

prove2quantum

The objects the formalized statements of Quantum Markov chains, sufficiency of quantum channels, and Renyi information measures (arXiv:1501.05636) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • IsStrictChannel: Section 3, PDF p.7: strictness on every positive definite input.
  • IsSufficient: Definition 1, PDF p.3, including its extension to unnormalized operators: one common CPTP recovery must recover both members of the pair.
  • renyiKernel: Equation (1.16), PDF p.4: the operator inside the trace after rho^alpha.
  • sandwichedKernel: Equation (1.17), PDF p.4: the sandwiched reference kernel, before rho factors.
  • renyiReference: Equations (3.23),(3.57), PDF pp.10,14: the computed reference R_alpha.
  • sandwichedReference: Equations (3.24),(3.59), PDF pp.10,14: the computed reference S_alpha.
  • renyiDifference: Equation (1.16), PDF p.4. The extended logarithm is -infinity at raw trace zero, so below order one the difference is +infinity, not zero.
  • deltaTilde: Equations (1.17),(A.17), PDF pp.4,21: for the positive operator X inside the alpha-norm, alpha log ||X||_alpha = log Tr X^alpha. This retains the primary normalization, rather than the erroneous auxiliary prefactor in (A.19), p.22.
  • minRelative: Equations (3.2),(3.8), PDF pp.7–8: -log of squared root fidelity. The square-root sandwich formula equals the squared trace norm in (3.2).
  • deltaMin: Definition 16, equation (5.1), PDF p.17: actual Petz-recovered reference.
  • deltaMax: Definition 16, equation (5.2), PDF p.17: Dmax as in primary (3.9), in bits.
  • conditionalSigma: Equation (1.21), PDF p.5: sigma=rho_AC tensor I_B, in coordinates A x (C x B). The AC marginal is computed from this same rho; the physical middle register is C.
  • conditionalRenyiInformation: Equations (1.14),(1.22), PDF pp.4–5, in the exact channel specialization (1.21).
  • conditionalRenyiTilde: Equations (1.15),(1.23), PDF pp.4–5, with the same computed marginal reference.
  • conditionalPetzReference: Equations (5.4),(5.7), PDF pp.17–18: the reference equals rho_AC^(1/2) rho_C^(-1/2) rho_BC rho_C^(-1/2) rho_AC^(1/2), with identities on complementary registers. It is constructed by (1.11),(1.21).
  • conditionalMinInformation: Equation (5.7), PDF p.18: primary min-information definition.
  • conditionalMaxInformation: Equation (5.4), PDF p.17: primary max-information definition, using 2^lambda from (3.9), rather than the inconsistent exp(lambda) in the alternative (5.5).
  • outputReference: Appendix A, equation (A.1), PDF p.19: output-space reference P_alpha. Both outer N(sigma) exponents are (alpha-1)/2. No trace normalization is applied.
  • jointTrace: Appendix A, equations (A.10)–(A.12), PDF p.20: power stays inside the trace.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P21
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P35
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P71
import Definitions.Def_QAlgorithms_Papers_Independent_Renew46_P35

/-!
Datta–Wilde, arXiv:1501.05636v4. Finite nonzero physical registers throughout.
Information is in bits; operator logs in the exponential inequality are natural.
All real spectral powers, including exponent zero, vanish on the kernel.
Compilation is pending the requested batch handoff.
-/

noncomputable section

open scoped Matrix ComplexOrder

namespace QAlgorithms.Papers.Independent.Intakethird100.P52

attribute [local instance] Classical.propDecidable Classical.decEq

open QAlgorithms.Papers.Independent.Information100.P63 (Mat)
open QAlgorithms.Papers.Independent.Dependency50.P21 (channelAdjoint petz marginalAB)
open QAlgorithms.Papers.Independent.Dependency50.P35 (tr elog maxRelative)
open QAlgorithms.Papers.Independent.Information100.P71 (supportPower)
open QAlgorithms.Papers.Independent.Renew46.P35 (HasMarkovBlocks)

/-- Section 3, PDF p.7: strictness on every positive definite input. -/
def IsStrictChannel {H K : Type} [Fintype H] [Fintype K] [DecidableEq H]
    (N : Superop H K) : Prop :=
  IsChannel N ∧ ∀ X : Mat H, X.PosDef → (N X).PosDef

/-- Definition 1, PDF p.3, including its extension to unnormalized operators:
one common CPTP recovery must recover both members of the pair. -/
def IsSufficient {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (rho sigma : Mat H) (N : Superop H K) : Prop :=
  ∃ R : Superop K H, IsChannel R ∧ R (N rho) = rho ∧ R (N sigma) = sigma

/-- Equation (1.16), PDF p.4: the operator inside the trace after rho^alpha. -/
def renyiKernel {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : Mat H :=
  supportPower sigma ((1 - a) / 2) *
    channelAdjoint N (supportPower (N sigma) ((a - 1) / 2) *
      supportPower (N rho) (1 - a) * supportPower (N sigma) ((a - 1) / 2)) *
    supportPower sigma ((1 - a) / 2)

/-- Equation (1.17), PDF p.4: the sandwiched reference kernel, before rho factors. -/
def sandwichedKernel {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : Mat H :=
  supportPower sigma ((1 - a) / (2 * a)) *
    channelAdjoint N (supportPower (N sigma) ((a - 1) / (2 * a)) *
      supportPower (N rho) ((1 - a) / a) * supportPower (N sigma) ((a - 1) / (2 * a))) *
    supportPower sigma ((1 - a) / (2 * a))

/-- Equations (3.23),(3.57), PDF pp.10,14: the computed reference R_alpha. -/
def renyiReference {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : Mat H :=
  supportPower (renyiKernel a rho sigma N) (1 / (1 - a))

/-- Equations (3.24),(3.59), PDF pp.10,14: the computed reference S_alpha. -/
def sandwichedReference {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : Mat H :=
  supportPower (sandwichedKernel a rho sigma N) (a / (1 - a))

/-- Equation (1.16), PDF p.4. The extended logarithm is -infinity at raw trace
zero, so below order one the difference is +infinity, not zero. -/
def renyiDifference {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : EReal :=
  ((1 / (a - 1) : ℝ) : EReal) *
    elog (tr (supportPower rho a * renyiKernel a rho sigma N))

/-- Equations (1.17),(A.17), PDF pp.4,21: for the positive operator X inside
the alpha-norm, alpha log ||X||_alpha = log Tr X^alpha. This retains the primary
normalization, rather than the erroneous auxiliary prefactor in (A.19), p.22. -/
def deltaTilde {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : EReal :=
  ((1 / (a - 1) : ℝ) : EReal) * elog (tr (supportPower
    (supportPower rho (1 / 2) * sandwichedKernel a rho sigma N * supportPower rho (1 / 2)) a))

/-- Equations (3.2),(3.8), PDF pp.7–8: -log of squared root fidelity.
The square-root sandwich formula equals the squared trace norm in (3.2). -/
def minRelative {H : Type} [Fintype H] [DecidableEq H]
    (rho tau : Mat H) : EReal :=
  -elog (rootFidelity rho tau ^ 2)

/-- Definition 16, equation (5.1), PDF p.17: actual Petz-recovered reference. -/
def deltaMin {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (rho sigma : Mat H) (N : Superop H K) : EReal :=
  minRelative rho (petz sigma N (N rho))

/-- Definition 16, equation (5.2), PDF p.17: Dmax as in primary (3.9), in bits. -/
def deltaMax {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (rho sigma : Mat H) (N : Superop H K) : EReal :=
  maxRelative rho (petz sigma N (N rho))

/-- Equation (1.21), PDF p.5: sigma=rho_AC tensor I_B, in coordinates A x (C x B).
The AC marginal is computed from this same rho; the physical middle register is C. -/
def conditionalSigma {A C B : Type} [Fintype B] [DecidableEq B]
    (rho : Mat (A × (C × B))) : Mat (A × (C × B)) :=
  fun i j => if i.2.2 = j.2.2 then marginalAB rho (i.1, i.2.1) (j.1, j.2.1) else 0

/-- Equations (1.14),(1.22), PDF pp.4–5, in the exact channel specialization (1.21). -/
def conditionalRenyiInformation {A C B : Type} [Fintype A] [Fintype C] [Fintype B]
    [DecidableEq A] [DecidableEq C] [DecidableEq B]
    (a : ℝ) (rho : Mat (A × (C × B))) : EReal :=
  renyiDifference a rho (conditionalSigma rho) traceLeft

/-- Equations (1.15),(1.23), PDF pp.4–5, with the same computed marginal reference. -/
def conditionalRenyiTilde {A C B : Type} [Fintype A] [Fintype C] [Fintype B]
    [DecidableEq A] [DecidableEq C] [DecidableEq B]
    (a : ℝ) (rho : Mat (A × (C × B))) : EReal :=
  deltaTilde a rho (conditionalSigma rho) traceLeft

/-- Equations (5.4),(5.7), PDF pp.17–18: the reference equals
rho_AC^(1/2) rho_C^(-1/2) rho_BC rho_C^(-1/2) rho_AC^(1/2),
with identities on complementary registers. It is constructed by (1.11),(1.21). -/
def conditionalPetzReference {A C B : Type} [Fintype A] [Fintype C] [Fintype B]
    [DecidableEq A] [DecidableEq C] [DecidableEq B]
    (rho : Mat (A × (C × B))) : Mat (A × (C × B)) :=
  petz (conditionalSigma rho) traceLeft (traceLeft rho)

/-- Equation (5.7), PDF p.18: primary min-information definition. -/
def conditionalMinInformation {A C B : Type} [Fintype A] [Fintype C] [Fintype B]
    [DecidableEq A] [DecidableEq C] [DecidableEq B]
    (rho : Mat (A × (C × B))) : EReal :=
  minRelative rho (conditionalPetzReference rho)

/-- Equation (5.4), PDF p.17: primary max-information definition, using 2^lambda
from (3.9), rather than the inconsistent exp(lambda) in the alternative (5.5). -/
def conditionalMaxInformation {A C B : Type} [Fintype A] [Fintype C] [Fintype B]
    [DecidableEq A] [DecidableEq C] [DecidableEq B]
    (rho : Mat (A × (C × B))) : EReal :=
  maxRelative rho (conditionalPetzReference rho)

/-- Appendix A, equation (A.1), PDF p.19: output-space reference P_alpha.
Both outer N(sigma) exponents are (alpha-1)/2. No trace normalization is applied. -/
def outputReference {H K : Type} [Fintype H] [Fintype K]
    [DecidableEq H] [DecidableEq K] (a : ℝ) (rho sigma : Mat H)
    (N : Superop H K) : Mat K :=
  supportPower (supportPower (N sigma) ((a - 1) / 2) *
    N (supportPower sigma ((1 - a) / 2) * supportPower rho a * supportPower sigma ((1 - a) / 2)) *
    supportPower (N sigma) ((a - 1) / 2)) (1 / a)

/-- Appendix A, equations (A.10)–(A.12), PDF p.20: power stays inside the trace. -/
def jointTrace {H : Type} [Fintype H] [DecidableEq H]
    (a : ℝ) (P Q : Mat H) : ℝ :=
  tr (supportPower (supportPower Q ((a - 1) / 2) * P * supportPower Q ((a - 1) / 2)) (1 / a))

/-! Sanity checks for the kernel convention, exponential convention, zero references,
and the conditioning-register order. These prove only definition sanity, not the paper results. -/

example {H : Type} [Fintype H] [DecidableEq H] (t : ℝ) :
    supportPower (0 : Mat H) t = 0 := by
  simp [supportPower]

example {H : Type} [Fintype H] [DecidableEq H] :
    cfc Real.exp (0 : Mat H) = 1 := by
  simp

example {H K : Type} [Fintype H] [Fintype K] [DecidableEq H] [DecidableEq K]
    (a : ℝ) (rho : Mat H) (N : Superop H K) :
    renyiReference a rho 0 N = 0 ∧ sandwichedReference a rho 0 N = 0 := by
  simp [renyiReference, sandwichedReference, renyiKernel, sandwichedKernel,
    supportPower]

example {H : Type} [Fintype H] [DecidableEq H] (rho : Mat H) :
    minRelative rho 0 = ⊤ := by
  simp [minRelative, rootFidelity, elog]

example (rho : Mat (Unit × (Unit × Unit))) :
    conditionalSigma rho ((), ((), ())) ((), ((), ())) =
      rho ((), ((), ())) ((), ((), ())) := by
  simp [conditionalSigma, marginalAB]

example {H : Type} [Fintype H] [DecidableEq H] (a : ℝ) (Q : Mat H) :
    jointTrace a 0 Q = 0 := by
  simp [jointTrace, supportPower, tr]

example {H : Type} [Fintype H] [DecidableEq H]
    (hid : IsChannel (id : Superop H H)) : IsStrictChannel (id : Superop H H) := by
  exact ⟨hid, fun _ h => h⟩

example {H : Type} [Fintype H] [DecidableEq H]
    (rho sigma : Mat H) (hid : IsChannel (id : Superop H H)) :
    IsSufficient rho sigma (id : Superop H H) := by
  exact ⟨id, hid, rfl, rfl⟩

example {H : Type} [Fintype H] [DecidableEq H] (a : ℝ) (P : Mat H) :
    jointTrace a P 0 = 0 := by
  simp [jointTrace, supportPower, tr]

end QAlgorithms.Papers.Independent.Intakethird100.P52
Source
Quantum Markov chains, sufficiency of quantum channels, and Renyi information measures, arXiv:1501.05636v4, https://arxiv.org/abs/1501.05636v4
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P52.IsStrictChannel

For finite, possibly empty types HHH and KKK, with decidable equality on HHH, and an arbitrary function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), this predicate means that NNN has a finite Kraus representation and sends every positive definite matrix to a positive definite matrix. Explicitly, there exist m∈Nm\in\mathbb Nm∈N and matrices Ej∈CK×HE_j\in\mathbb C^{K\times H}Ej​∈CK×H, indexed by 0≤j<m0\le j<m0≤j<m, such that ∑jEj†Ej=IH\sum_jE_j^\dagger E_j=I_H∑j​Ej†​Ej​=IH​ and N(X)=∑jEjXEj†N(X)=\sum_jE_jXE_j^\daggerN(X)=∑j​Ej​XEj†​ for every complex matrix XXX on HHH; additionally, X≻0X\succ0X≻0 implies N(X)≻0N(X)\succ0N(X)≻0 for every such XXX. The quantified matrices need not have trace one, and the Kraus family is allowed to be empty.

QAlgorithms.Papers.Independent.Intakethird100.P52.IsSufficient

For finite, possibly empty types H,KH,KH,K with decidable equality, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and an arbitrary function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), this predicate means that there exists a function R:Mat⁡K(C)→Mat⁡H(C)R:\operatorname{Mat}_K(\mathbb C)\to\operatorname{Mat}_H(\mathbb C)R:MatK​(C)→MatH​(C) satisfying R(N(ρ))=ρR(N(\rho))=\rhoR(N(ρ))=ρ and R(N(σ))=σR(N(\sigma))=\sigmaR(N(σ))=σ, where RRR has a finite Kraus representation: for some m∈Nm\in\mathbb Nm∈N and matrices Ej∈CH×KE_j\in\mathbb C^{H\times K}Ej​∈CH×K, ∑jEj†Ej=IK\sum_jE_j^\dagger E_j=I_K∑j​Ej†​Ej​=IK​ and R(Y)=∑jEjYEj†R(Y)=\sum_jE_jYE_j^\daggerR(Y)=∑j​Ej​YEj†​ for every YYY. No positivity or normalization is required of ρ,σ\rho,\sigmaρ,σ, and no channel or linearity assumption is imposed on NNN.

QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalMaxInformation

For finite, possibly empty types A,C,BA,C,BA,C,B, each with decidable equality, and any complex matrix ρ\rhoρ on A×(C×B)A\times(C\times B)A×(C×B), set Σ=ρAC⊗IB\Sigma=\rho_{AC}\otimes I_BΣ=ρAC​⊗IB​, where (ρAC)(i,c),(j,d)=∑bρ(i,c,b),(j,d,b)(\rho_{AC})_{(i,c),(j,d)}=\sum_b\rho_{(i,c,b),(j,d,b)}(ρAC​)(i,c),(j,d)​=∑b​ρ(i,c,b),(j,d,b)​, and put D=Tr⁡AΣD=\operatorname{Tr}_A\SigmaD=TrA​Σ, Y=Tr⁡AρY=\operatorname{Tr}_A\rhoY=TrA​ρ, with partial trace summing equal AAA-indices. Let Q(X)Q(X)Q(X) and J(X)J(X)J(X) be continuous functional calculus for x\sqrt{x}x​ and (x)−1(\sqrt{x})^{-1}(x​)−1, respectively, with inverse zero at zero, and define

Rρ=Q(Σ)[IA⊗J(D)YJ(D)]Q(Σ).R_\rho=Q(\Sigma)\left[I_A\otimes J(D)YJ(D)\right]Q(\Sigma).Rρ​=Q(Σ)[IA​⊗J(D)YJ(D)]Q(Σ).

The declaration defines

inf⁡{λ∈R:2λRρ−ρ is positive semidefinite}\inf\left\{\lambda\in\mathbb R: 2^\lambda R_\rho-\rho\ \text{is positive semidefinite}\right\}inf{λ∈R:2λRρ​−ρ is positive semidefinite}

as an extended-real infimum over the embedded real numbers. An empty feasible set gives +∞+\infty+∞, and a feasible set unbounded below gives −∞-\infty−∞. No positivity or normalization assumption on ρ\rhoρ is required.

QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalMinInformation

For finite, possibly empty types A,C,BA,C,BA,C,B, each with decidable equality, and any complex matrix ρ\rhoρ on A×(C×B)A\times(C\times B)A×(C×B), put Σ=ρAC⊗IB\Sigma=\rho_{AC}\otimes I_BΣ=ρAC​⊗IB​, where (ρAC)(i,c),(j,d)=∑bρ(i,c,b),(j,d,b)(\rho_{AC})_{(i,c),(j,d)}=\sum_b\rho_{(i,c,b),(j,d,b)}(ρAC​)(i,c),(j,d)​=∑b​ρ(i,c,b),(j,d,b)​, and let D=Tr⁡AΣD=\operatorname{Tr}_A\SigmaD=TrA​Σ, Y=Tr⁡AρY=\operatorname{Tr}_A\rhoY=TrA​ρ, with partial trace summing equal AAA-indices. Write Q(X)Q(X)Q(X) for continuous functional calculus for x↦xx\mapsto\sqrt{x}x↦x​, and J(X)J(X)J(X) for that for x↦(x)−1x\mapsto(\sqrt{x})^{-1}x↦(x​)−1, whose value at zero is zero, and define

Rρ=Q(Σ)[IA⊗J(D)YJ(D)]Q(Σ).R_\rho=Q(\Sigma)\left[I_A\otimes J(D)YJ(D)\right]Q(\Sigma).Rρ​=Q(Σ)[IA​⊗J(D)YJ(D)]Q(Σ).

The declaration defines the extended-real quantity

−ℓ2 ⁣([Re⁡Tr⁡Q ⁣(Q(ρ)RρQ(ρ))]2),-\ell_2\!\left( \left[\operatorname{Re}\operatorname{Tr} Q\!\left(Q(\rho)R_\rho Q(\rho)\right)\right]^2 \right),−ℓ2​([ReTrQ(Q(ρ)Rρ​Q(ρ))]2),

where ℓ2(x)=log⁡2x\ell_2(x)=\log_2xℓ2​(x)=log2​x for x>0x>0x>0 and −∞-\infty−∞ for x≤0x\le0x≤0. In particular, a zero squared expression gives +∞+\infty+∞. There is no positivity or normalization assumption on ρ\rhoρ.

QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalPetzReference

For finite, possibly empty types A,C,BA,C,BA,C,B, each with decidable equality, and any complex matrix ρ\rhoρ on A×(C×B)A\times(C\times B)A×(C×B), define Σ=ρAC⊗IB\Sigma=\rho_{AC}\otimes I_BΣ=ρAC​⊗IB​, where (ρAC)(i,c),(j,d)=∑bρ(i,c,b),(j,d,b)(\rho_{AC})_{(i,c),(j,d)}=\sum_b\rho_{(i,c,b),(j,d,b)}(ρAC​)(i,c),(j,d)​=∑b​ρ(i,c,b),(j,d,b)​, and let D=Tr⁡AΣD=\operatorname{Tr}_A\SigmaD=TrA​Σ, Y=Tr⁡AρY=\operatorname{Tr}_A\rhoY=TrA​ρ, where (Tr⁡AX)(c,b),(d,e)=∑iX(i,c,b),(i,d,e)(\operatorname{Tr}_A X)_{(c,b),(d,e)}=\sum_iX_{(i,c,b),(i,d,e)}(TrA​X)(c,b),(d,e)​=∑i​X(i,c,b),(i,d,e)​. If Q(X)Q(X)Q(X) denotes continuous functional calculus for x↦xx\mapsto\sqrt{x}x↦x​, and J(X)J(X)J(X) denotes continuous functional calculus for x↦(x)−1x\mapsto(\sqrt{x})^{-1}x↦(x​)−1, the declaration defines the matrix

Rρ=Q(Σ)[IA⊗(J(D)YJ(D))]Q(Σ).R_\rho=Q(\Sigma)\left[I_A\otimes\bigl(J(D)YJ(D)\bigr)\right]Q(\Sigma).Rρ​=Q(Σ)[IA​⊗(J(D)YJ(D))]Q(Σ).

No density-matrix or positivity hypothesis is imposed. The scalar square root and inverse are their total real operations, so the inverse-square-root function takes value zero at zero.

QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalRenyiInformation

For finite types A,C,BA,C,BA,C,B, each with decidable equality and possibly empty, any real aaa, and any complex matrix ρ\rhoρ indexed by A×(C×B)A\times(C\times B)A×(C×B), define Σ=ρAC⊗IB\Sigma=\rho_{AC}\otimes I_BΣ=ρAC​⊗IB​, where

(ρAC)(i,c),(j,d)=∑bρ(i,c,b),(j,d,b).(\rho_{AC})_{(i,c),(j,d)}=\sum_b\rho_{(i,c,b),(j,d,b)}.(ρAC​)(i,c),(j,d)​=b∑​ρ(i,c,b),(j,d,b)​.

Write ρCB=Tr⁡Aρ\rho_{CB}=\operatorname{Tr}_A\rhoρCB​=TrA​ρ and ΣCB=Tr⁡AΣ\Sigma_{CB}=\operatorname{Tr}_A\SigmaΣCB​=TrA​Σ, with (Tr⁡AX)(c,b),(d,e)=∑iX(i,c,b),(i,d,e)(\operatorname{Tr}_A X)_{(c,b),(d,e)}=\sum_i X_{(i,c,b),(i,d,e)}(TrA​X)(c,b),(d,e)​=∑i​X(i,c,b),(i,d,e)​. For any square matrix XXX, let St(X)S_t(X)St​(X) denote continuous functional calculus applied to the real function that is 000 at x=0x=0x=0 and xtx^txt otherwise, and set

Ka= ⁣S(1−a)/2(Σ)[IA⊗(S(a−1)/2(ΣCB)S1−a(ρCB)S(a−1)/2(ΣCB))]S(1−a)/2(Σ).K_a=\! S_{(1-a)/2}(\Sigma) \left[I_A\otimes \left( S_{(a-1)/2}(\Sigma_{CB})S_{1-a}(\rho_{CB})S_{(a-1)/2}(\Sigma_{CB}) \right)\right] S_{(1-a)/2}(\Sigma).Ka​=S(1−a)/2​(Σ)[IA​⊗(S(a−1)/2​(ΣCB​)S1−a​(ρCB​)S(a−1)/2​(ΣCB​))]S(1−a)/2​(Σ).

The declaration defines the extended-real value

1a−1 ℓ2 ⁣(Re⁡Tr⁡(Sa(ρ)Ka)),ℓ2(x)={log⁡2x,x>0,−∞,x≤0.\frac{1}{a-1}\, \ell_2\!\left(\operatorname{Re}\operatorname{Tr}(S_a(\rho)K_a)\right), \qquad \ell_2(x)= \begin{cases}\log_2x,&x>0,\\-\infty,&x\le0.\end{cases}a−11​ℓ2​(ReTr(Sa​(ρ)Ka​)),ℓ2​(x)={log2​x,−∞,​x>0,x≤0.​

There are no positivity, trace, or parameter restrictions: a=0,1a=0,1a=0,1 are included, real division by zero is defined as zero, extended-real arithmetic is used, and StS_tSt​ assigns zero to a zero spectral value even when t=0t=0t=0.

QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalRenyiTilde

For finite, possibly empty types A,C,BA,C,BA,C,B, each with decidable equality, any real aaa, and any complex matrix ρ\rhoρ on A×(C×B)A\times(C\times B)A×(C×B), set Σ=ρAC⊗IB\Sigma=\rho_{AC}\otimes I_BΣ=ρAC​⊗IB​, where (ρAC)(i,c),(j,d)=∑bρ(i,c,b),(j,d,b)(\rho_{AC})_{(i,c),(j,d)}=\sum_b\rho_{(i,c,b),(j,d,b)}(ρAC​)(i,c),(j,d)​=∑b​ρ(i,c,b),(j,d,b)​, and set ρCB=Tr⁡Aρ\rho_{CB}=\operatorname{Tr}_A\rhoρCB​=TrA​ρ, ΣCB=Tr⁡AΣ\Sigma_{CB}=\operatorname{Tr}_A\SigmaΣCB​=TrA​Σ, where partial trace sums over equal AAA-indices. Let St(X)S_t(X)St​(X) be continuous functional calculus for the function equal to 000 at zero and xtx^txt otherwise. Define

K~a=S(1−a)/(2a)(Σ)[IA⊗(S(a−1)/(2a)(ΣCB)S(1−a)/a(ρCB)S(a−1)/(2a)(ΣCB))]S(1−a)/(2a)(Σ).\widetilde K_a= S_{(1-a)/(2a)}(\Sigma) \left[I_A\otimes \left( S_{(a-1)/(2a)}(\Sigma_{CB}) S_{(1-a)/a}(\rho_{CB}) S_{(a-1)/(2a)}(\Sigma_{CB}) \right)\right] S_{(1-a)/(2a)}(\Sigma).Ka​=S(1−a)/(2a)​(Σ)[IA​⊗(S(a−1)/(2a)​(ΣCB​)S(1−a)/a​(ρCB​)S(a−1)/(2a)​(ΣCB​))]S(1−a)/(2a)​(Σ).

The declaration defines

1a−1 ℓ2 ⁣(Re⁡Tr⁡Sa ⁣(S1/2(ρ)K~aS1/2(ρ)))\frac{1}{a-1}\, \ell_2\!\left( \operatorname{Re}\operatorname{Tr} S_a\!\left(S_{1/2}(\rho)\widetilde K_aS_{1/2}(\rho)\right) \right)a−11​ℓ2​(ReTrSa​(S1/2​(ρ)Ka​S1/2​(ρ)))

as an extended real, where ℓ2(x)=log⁡2x\ell_2(x)=\log_2xℓ2​(x)=log2​x for x>0x>0x>0 and ℓ2(x)=−∞\ell_2(x)=-\inftyℓ2​(x)=−∞ otherwise. No positivity or normalization of ρ\rhoρ, or exclusion of a=0,1a=0,1a=0,1, is assumed; real division by zero is zero, and the functional calculus function is zero at a zero spectral value for every exponent.

QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalSigma

For arbitrary types A,CA,CA,C, a finite, possibly empty type BBB with decidable equality, and any complex matrix ρ\rhoρ indexed by A×(C×B)A\times(C\times B)A×(C×B), this declaration defines a matrix on the same index set by

Σ(a,c,b),(a′,c′,b′)={∑u∈Bρ(a,c,u),(a′,c′,u),b=b′,0,b≠b′.\Sigma_{(a,c,b),(a',c',b')}= \begin{cases} \displaystyle\sum_{u\in B} \rho_{(a,c,u),(a',c',u)},&b=b',\\ 0,&b\ne b'. \end{cases}Σ(a,c,b),(a′,c′,b′)​=⎩⎨⎧​u∈B∑​ρ(a,c,u),(a′,c′,u)​,0,​b=b′,b=b′.​

Thus each diagonal block in the BBB register is the matrix obtained by summing over equal BBB indices of ρ\rhoρ. Neither AAA nor CCC is required to be finite, and no positivity, Hermiticity, or normalization is assumed. If BBB is empty, the matrix has an empty index set.

QAlgorithms.Papers.Independent.Intakethird100.P52.deltaMax

For finite, possibly empty types H,KH,KH,K with decidable equality, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), write Ff(X)F_f(X)Ff​(X) for continuous functional calculus and define (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​ using matrix units. Set S=Fx(σ)S=F_{\sqrt{\phantom{x}}}(\sigma)S=Fx​​(σ), T=Fx↦(x)−1(N(σ))T=F_{x\mapsto(\sqrt{x})^{-1}}(N(\sigma))T=Fx↦(x​)−1​(N(σ)), and τ=SDN(TN(ρ)T)S\tau=S D_N(TN(\rho)T)Sτ=SDN​(TN(ρ)T)S. This declaration defines

inf⁡{ℓ∈R:  2ℓτ−ρ⪰0},\inf\left\{\ell\in\mathbb R:\;2^\ell\tau-\rho\succeq0\right\},inf{ℓ∈R:2ℓτ−ρ⪰0},

where the real numbers in this set are embedded into the extended reals before taking the infimum. No positivity, normalization, channel, or linearity assumptions are imposed; the inverse square-root function takes value zero at zero. An empty feasible set gives +∞+\infty+∞, whereas a feasible set unbounded below gives −∞-\infty−∞; in particular, an empty HHH gives −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P52.deltaMin

For finite, possibly empty types H,KH,KH,K with decidable equality, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), write Ff(X)F_f(X)Ff​(X) for continuous functional calculus and define (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​ using matrix units. Set S=Fx(σ)S=F_{\sqrt{\phantom{x}}}(\sigma)S=Fx​​(σ), T=Fx↦(x)−1(N(σ))T=F_{x\mapsto(\sqrt{x})^{-1}}(N(\sigma))T=Fx↦(x​)−1​(N(σ)), and τ=SDN(TN(ρ)T)S\tau=S D_N(TN(\rho)T)Sτ=SDN​(TN(ρ)T)S. The declaration defines −L(r2)-L(r^2)−L(r2), where

r=Re⁡Tr⁡Fx ⁣(Fx(ρ) τ Fx(ρ)),L(x)={log⁡2x,x>0,−∞,x≤0.r=\operatorname{Re}\operatorname{Tr} F_{\sqrt{\phantom{x}}}\!\left( F_{\sqrt{\phantom{x}}}(\rho)\,\tau\,F_{\sqrt{\phantom{x}}}(\rho) \right), \qquad L(x)= \begin{cases} \log_2x,&x>0,\\ -\infty,&x\le0. \end{cases}r=ReTrFx​​(Fx​​(ρ)τFx​​(ρ)),L(x)={log2​x,−∞,​x>0,x≤0.​

No positivity, normalization, channel, or linearity assumptions are imposed. The inverse square-root function takes value zero at zero under real inversion, and the defined extended-real value is +∞+\infty+∞ whenever r=0r=0r=0, including the empty-HHH case.

QAlgorithms.Papers.Independent.Intakethird100.P52.deltaTilde

For finite, possibly empty types H,KH,KH,K with decidable equality, any a∈Ra\in\mathbb Ra∈R, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), let St(X)S_t(X)St​(X) be continuous functional calculus for the function equal to 000 at zero and rpow⁡(x,t)\operatorname{rpow}(x,t)rpow(x,t) otherwise, and define (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​ using matrix units. Set

B=S(1−a)/(2a)(σ) DN ⁣(S(a−1)/(2a)(N(σ)) S(1−a)/a(N(ρ)) S(a−1)/(2a)(N(σ))) S(1−a)/(2a)(σ).B= S_{(1-a)/(2a)}(\sigma)\, D_N\!\left( S_{(a-1)/(2a)}(N(\sigma))\,S_{(1-a)/a}(N(\rho))\, S_{(a-1)/(2a)}(N(\sigma)) \right)\, S_{(1-a)/(2a)}(\sigma).B=S(1−a)/(2a)​(σ)DN​(S(a−1)/(2a)​(N(σ))S(1−a)/a​(N(ρ))S(a−1)/(2a)​(N(σ)))S(1−a)/(2a)​(σ).

The declaration defines 1a−1L ⁣(Re⁡Tr⁡Sa(S1/2(ρ)BS1/2(ρ)))\frac1{a-1}L\!\left(\operatorname{Re}\operatorname{Tr}S_a(S_{1/2}(\rho)BS_{1/2}(\rho))\right)a−11​L(ReTrSa​(S1/2​(ρ)BS1/2​(ρ))) in the extended reals, where L(x)=log⁡2xL(x)=\log_2xL(x)=log2​x for x>0x>0x>0 and L(x)=−∞L(x)=-\inftyL(x)=−∞ otherwise, and the real prefactor is embedded before multiplication. No positivity, normalization, channel, linearity, or parameter restrictions are imposed. Division by zero is zero, so a=0a=0a=0 is included, and at a=1a=1a=1 the resulting extended-real value is zero.

QAlgorithms.Papers.Independent.Intakethird100.P52.jointTrace

For any finite, possibly empty type HHH with decidable equality, any real aaa, and any complex square matrices P,QP,QP,Q indexed by HHH, let St(X)S_t(X)St​(X) denote continuous functional calculus for the real function equal to zero at zero and xtx^txt otherwise. The declaration defines the real number

Re⁡Tr⁡S1/a ⁣(S(a−1)/2(Q) P S(a−1)/2(Q)).\operatorname{Re}\operatorname{Tr} S_{1/a}\!\left(S_{(a-1)/2}(Q)\,P\,S_{(a-1)/2}(Q)\right).ReTrS1/a​(S(a−1)/2​(Q)PS(a−1)/2​(Q)).

No positivity or normalization hypotheses are imposed. The parameter a=0a=0a=0 is included through real division by zero being zero, and the function defining S0S_0S0​ has value zero at zero.

QAlgorithms.Papers.Independent.Intakethird100.P52.minRelative

For a finite, possibly empty type HHH with decidable equality and arbitrary complex matrices ρ,τ\rho,\tauρ,τ on HHH, let Ff(X)F_f(X)Ff​(X) denote continuous functional calculus with the real function fff, and put

r=Re⁡Tr⁡Fx ⁣(Fx(ρ) τ Fx(ρ)).r=\operatorname{Re}\operatorname{Tr} F_{\sqrt{\phantom{x}}}\!\left( F_{\sqrt{\phantom{x}}}(\rho)\,\tau\,F_{\sqrt{\phantom{x}}}(\rho) \right).r=ReTrFx​​(Fx​​(ρ)τFx​​(ρ)).

This declaration defines −L(r2)-L(r^2)−L(r2) in the extended reals, where L(x)=log⁡2xL(x)=\log_2xL(x)=log2​x when x>0x>0x>0 and L(x)=−∞L(x)=-\inftyL(x)=−∞ when x≤0x\le0x≤0. Thus the value is −log⁡2(r2)-\log_2(r^2)−log2​(r2) if r≠0r\ne0r=0 and +∞+\infty+∞ if r=0r=0r=0. No positivity or trace normalization is assumed; in particular, the empty-index case gives +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P52.outputReference

For finite, possibly empty types H,KH,KH,K, each with decidable equality, any real aaa, any complex square matrices ρ,σ\rho,\sigmaρ,σ indexed by HHH, and any function NNN from such matrices to complex square matrices indexed by KKK, let St(X)S_t(X)St​(X) denote continuous functional calculus for the function equal to 000 at zero and xtx^txt otherwise. The declaration defines the matrix on KKK

S1/a ⁣(S(a−1)/2(N(σ)) N ⁣(S(1−a)/2(σ)Sa(ρ)S(1−a)/2(σ)) S(a−1)/2(N(σ))).S_{1/a}\!\left( S_{(a-1)/2}(N(\sigma))\, N\!\left(S_{(1-a)/2}(\sigma)S_a(\rho)S_{(1-a)/2}(\sigma)\right)\, S_{(a-1)/2}(N(\sigma)) \right).S1/a​(S(a−1)/2​(N(σ))N(S(1−a)/2​(σ)Sa​(ρ)S(1−a)/2​(σ))S(a−1)/2​(N(σ))).

There are no assumptions that NNN is linear or a channel, or that either matrix is positive or normalized. The definition includes a=0a=0a=0, using 1/0=01/0=01/0=0; S0S_0S0​ still assigns zero at a zero spectral value.

QAlgorithms.Papers.Independent.Intakethird100.P52.renyiDifference

For finite, possibly empty types H,KH,KH,K with decidable equality, any a∈Ra\in\mathbb Ra∈R, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), let St(X)S_t(X)St​(X) be continuous functional calculus for the function equal to 000 at zero and rpow⁡(x,t)\operatorname{rpow}(x,t)rpow(x,t) otherwise, and define (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​ using matrix units. With L(x)=log⁡2xL(x)=\log_2xL(x)=log2​x for x>0x>0x>0 and L(x)=−∞L(x)=-\inftyL(x)=−∞ for x≤0x\le0x≤0, this declaration defines the extended-real value

1a−1 L ⁣(Re⁡Tr⁡ ⁣[Sa(ρ) S(1−a)/2(σ) DN ⁣(S(a−1)/2(N(σ)) S1−a(N(ρ)) S(a−1)/2(N(σ))) S(1−a)/2(σ)]).\frac1{a-1}\, L\!\left( \operatorname{Re}\operatorname{Tr}\!\left[ S_a(\rho)\,S_{(1-a)/2}(\sigma)\, D_N\!\left( S_{(a-1)/2}(N(\sigma))\,S_{1-a}(N(\rho))\,S_{(a-1)/2}(N(\sigma)) \right)\, S_{(1-a)/2}(\sigma) \right]\right).a−11​L(ReTr[Sa​(ρ)S(1−a)/2​(σ)DN​(S(a−1)/2​(N(σ))S1−a​(N(ρ))S(a−1)/2​(N(σ)))S(1−a)/2​(σ)]).

The real prefactor is embedded into the extended reals before multiplication. There are no positivity, normalization, channel, linearity, or parameter restrictions; at a=1a=1a=1 the prefactor is zero, so the defined value is zero, including when the logarithm is −∞-\infty−∞.

QAlgorithms.Papers.Independent.Intakethird100.P52.renyiKernel

For finite, possibly empty types H,KH,KH,K with decidable equality, any a∈Ra\in\mathbb Ra∈R, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), define St(X)S_t(X)St​(X) by continuous functional calculus with the real function that is 000 at x=0x=0x=0 and rpow⁡(x,t)\operatorname{rpow}(x,t)rpow(x,t) otherwise. Define DN(Y)D_N(Y)DN​(Y) entrywise by (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​, where EijE_{ij}Eij​ has entry 111 at (i,j)(i,j)(i,j) and zero elsewhere. The declaration defines the matrix

S(1−a)/2(σ) DN ⁣(S(a−1)/2(N(σ)) S1−a(N(ρ)) S(a−1)/2(N(σ))) S(1−a)/2(σ).S_{(1-a)/2}(\sigma)\, D_N\!\left( S_{(a-1)/2}(N(\sigma))\,S_{1-a}(N(\rho))\,S_{(a-1)/2}(N(\sigma)) \right)\, S_{(1-a)/2}(\sigma).S(1−a)/2​(σ)DN​(S(a−1)/2​(N(σ))S1−a​(N(ρ))S(a−1)/2​(N(σ)))S(1−a)/2​(σ).

There are no positivity, normalization, channel, or linearity assumptions. Every real aaa, including a=1a=1a=1, is included; S0S_0S0​ uses the function 000 at zero and 111 elsewhere, rather than the constant function 111.

QAlgorithms.Papers.Independent.Intakethird100.P52.renyiReference

For finite, possibly empty types H,KH,KH,K with decidable equality, any a∈Ra\in\mathbb Ra∈R, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), let St(X)S_t(X)St​(X) denote continuous functional calculus for the function equal to 000 at zero and rpow⁡(x,t)\operatorname{rpow}(x,t)rpow(x,t) otherwise, and define (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​ using matrix units EijE_{ij}Eij​. This declaration defines

S1/(1−a) ⁣(S(1−a)/2(σ) DN ⁣(S(a−1)/2(N(σ)) S1−a(N(ρ)) S(a−1)/2(N(σ))) S(1−a)/2(σ)).S_{1/(1-a)}\!\left( S_{(1-a)/2}(\sigma)\, D_N\!\left( S_{(a-1)/2}(N(\sigma))\,S_{1-a}(N(\rho))\,S_{(a-1)/2}(N(\sigma)) \right)\, S_{(1-a)/2}(\sigma) \right).S1/(1−a)​(S(1−a)/2​(σ)DN​(S(a−1)/2​(N(σ))S1−a​(N(ρ))S(a−1)/2​(N(σ)))S(1−a)/2​(σ)).

There are no positivity, normalization, channel, or linearity assumptions. In particular a=1a=1a=1 is allowed: the outer exponent is then zero under the convention 1/0=01/0=01/0=0, and S0S_0S0​ uses the function 000 at zero and 111 elsewhere.

QAlgorithms.Papers.Independent.Intakethird100.P52.sandwichedKernel

For finite, possibly empty types H,KH,KH,K with decidable equality, any a∈Ra\in\mathbb Ra∈R, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), let St(X)S_t(X)St​(X) be continuous functional calculus for the real function equal to 000 at zero and rpow⁡(x,t)\operatorname{rpow}(x,t)rpow(x,t) elsewhere, and put (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​, with EijE_{ij}Eij​ the matrix unit. The declaration defines

S(1−a)/(2a)(σ) DN ⁣(S(a−1)/(2a)(N(σ)) S(1−a)/a(N(ρ)) S(a−1)/(2a)(N(σ))) S(1−a)/(2a)(σ).S_{(1-a)/(2a)}(\sigma)\, D_N\!\left( S_{(a-1)/(2a)}(N(\sigma))\, S_{(1-a)/a}(N(\rho))\, S_{(a-1)/(2a)}(N(\sigma)) \right)\, S_{(1-a)/(2a)}(\sigma).S(1−a)/(2a)​(σ)DN​(S(a−1)/(2a)​(N(σ))S(1−a)/a​(N(ρ))S(a−1)/(2a)​(N(σ)))S(1−a)/(2a)​(σ).

No positivity, normalization, channel, or linearity assumptions are imposed. The definition includes a=0a=0a=0, where each displayed exponent is zero because real division by zero is defined as zero; S0S_0S0​ uses the function 000 at zero and 111 elsewhere. Negative exponents likewise retain the explicitly assigned value zero at a zero spectral argument.

QAlgorithms.Papers.Independent.Intakethird100.P52.sandwichedReference

For finite, possibly empty types H,KH,KH,K with decidable equality, any a∈Ra\in\mathbb Ra∈R, arbitrary complex matrices ρ,σ\rho,\sigmaρ,σ on HHH, and any function N:Mat⁡H(C)→Mat⁡K(C)N:\operatorname{Mat}_H(\mathbb C)\to\operatorname{Mat}_K(\mathbb C)N:MatH​(C)→MatK​(C), let St(X)S_t(X)St​(X) denote continuous functional calculus for the function equal to 000 at zero and rpow⁡(x,t)\operatorname{rpow}(x,t)rpow(x,t) otherwise, and set (DN(Y))ij=∑k,l∈KN(Eij)kl‾Ykl(D_N(Y))_{ij}=\sum_{k,l\in K}\overline{N(E_{ij})_{kl}}Y_{kl}(DN​(Y))ij​=∑k,l∈K​N(Eij​)kl​​Ykl​ using matrix units EijE_{ij}Eij​. The defined matrix is

Sa/(1−a) ⁣(S(1−a)/(2a)(σ) DN ⁣(S(a−1)/(2a)(N(σ)) S(1−a)/a(N(ρ)) S(a−1)/(2a)(N(σ))) S(1−a)/(2a)(σ)).S_{a/(1-a)}\!\left( S_{(1-a)/(2a)}(\sigma)\, D_N\!\left( S_{(a-1)/(2a)}(N(\sigma))\,S_{(1-a)/a}(N(\rho))\, S_{(a-1)/(2a)}(N(\sigma)) \right)\, S_{(1-a)/(2a)}(\sigma) \right).Sa/(1−a)​(S(1−a)/(2a)​(σ)DN​(S(a−1)/(2a)​(N(σ))S(1−a)/a​(N(ρ))S(a−1)/(2a)​(N(σ)))S(1−a)/(2a)​(σ)).

No positivity, normalization, channel, or linearity assumptions are imposed. Both a=0a=0a=0 and a=1a=1a=1 are included, with every division by zero evaluated as zero; an exponent zero invokes the function 000 at zero and 111 elsewhere.

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