Definitions for Quantum Markov chains, sufficiency of quantum channels, and Renyi information measures
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P52The 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.
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
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 and , with decidable equality on , and an arbitrary function , this predicate means that has a finite Kraus representation and sends every positive definite matrix to a positive definite matrix. Explicitly, there exist and matrices , indexed by , such that and for every complex matrix on ; additionally, implies for every such . 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 with decidable equality, arbitrary complex matrices on , and an arbitrary function , this predicate means that there exists a function satisfying and , where has a finite Kraus representation: for some and matrices , and for every . No positivity or normalization is required of , and no channel or linearity assumption is imposed on .
QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalMaxInformation
For finite, possibly empty types , each with decidable equality, and any complex matrix on , set , where , and put , , with partial trace summing equal -indices. Let and be continuous functional calculus for and , respectively, with inverse zero at zero, and define
The declaration defines
as an extended-real infimum over the embedded real numbers. An empty feasible set gives , and a feasible set unbounded below gives . No positivity or normalization assumption on is required.
QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalMinInformation
For finite, possibly empty types , each with decidable equality, and any complex matrix on , put , where , and let , , with partial trace summing equal -indices. Write for continuous functional calculus for , and for that for , whose value at zero is zero, and define
The declaration defines the extended-real quantity
where for and for . In particular, a zero squared expression gives . There is no positivity or normalization assumption on .
QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalPetzReference
For finite, possibly empty types , each with decidable equality, and any complex matrix on , define , where , and let , , where . If denotes continuous functional calculus for , and denotes continuous functional calculus for , the declaration defines the matrix
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 , each with decidable equality and possibly empty, any real , and any complex matrix indexed by , define , where
Write and , with . For any square matrix , let denote continuous functional calculus applied to the real function that is at and otherwise, and set
The declaration defines the extended-real value
There are no positivity, trace, or parameter restrictions: are included, real division by zero is defined as zero, extended-real arithmetic is used, and assigns zero to a zero spectral value even when .
QAlgorithms.Papers.Independent.Intakethird100.P52.conditionalRenyiTilde
For finite, possibly empty types , each with decidable equality, any real , and any complex matrix on , set , where , and set , , where partial trace sums over equal -indices. Let be continuous functional calculus for the function equal to at zero and otherwise. Define
The declaration defines
as an extended real, where for and otherwise. No positivity or normalization of , or exclusion of , 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 finite, possibly empty type with decidable equality, and any complex matrix indexed by , this declaration defines a matrix on the same index set by
Thus each diagonal block in the register is the matrix obtained by summing over equal indices of . Neither nor is required to be finite, and no positivity, Hermiticity, or normalization is assumed. If is empty, the matrix has an empty index set.
QAlgorithms.Papers.Independent.Intakethird100.P52.deltaMax
For finite, possibly empty types with decidable equality, arbitrary complex matrices on , and any function , write for continuous functional calculus and define using matrix units. Set , , and . This declaration defines
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 , whereas a feasible set unbounded below gives ; in particular, an empty gives .
QAlgorithms.Papers.Independent.Intakethird100.P52.deltaMin
For finite, possibly empty types with decidable equality, arbitrary complex matrices on , and any function , write for continuous functional calculus and define using matrix units. Set , , and . The declaration defines , where
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 whenever , including the empty- case.
QAlgorithms.Papers.Independent.Intakethird100.P52.deltaTilde
For finite, possibly empty types with decidable equality, any , arbitrary complex matrices on , and any function , let be continuous functional calculus for the function equal to at zero and otherwise, and define using matrix units. Set
The declaration defines in the extended reals, where for and 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 is included, and at the resulting extended-real value is zero.
QAlgorithms.Papers.Independent.Intakethird100.P52.jointTrace
For any finite, possibly empty type with decidable equality, any real , and any complex square matrices indexed by , let denote continuous functional calculus for the real function equal to zero at zero and otherwise. The declaration defines the real number
No positivity or normalization hypotheses are imposed. The parameter is included through real division by zero being zero, and the function defining has value zero at zero.
QAlgorithms.Papers.Independent.Intakethird100.P52.minRelative
For a finite, possibly empty type with decidable equality and arbitrary complex matrices on , let denote continuous functional calculus with the real function , and put
This declaration defines in the extended reals, where when and when . Thus the value is if and if . No positivity or trace normalization is assumed; in particular, the empty-index case gives .
QAlgorithms.Papers.Independent.Intakethird100.P52.outputReference
For finite, possibly empty types , each with decidable equality, any real , any complex square matrices indexed by , and any function from such matrices to complex square matrices indexed by , let denote continuous functional calculus for the function equal to at zero and otherwise. The declaration defines the matrix on
There are no assumptions that is linear or a channel, or that either matrix is positive or normalized. The definition includes , using ; still assigns zero at a zero spectral value.
QAlgorithms.Papers.Independent.Intakethird100.P52.renyiDifference
For finite, possibly empty types with decidable equality, any , arbitrary complex matrices on , and any function , let be continuous functional calculus for the function equal to at zero and otherwise, and define using matrix units. With for and for , this declaration defines the extended-real value
The real prefactor is embedded into the extended reals before multiplication. There are no positivity, normalization, channel, linearity, or parameter restrictions; at the prefactor is zero, so the defined value is zero, including when the logarithm is .
QAlgorithms.Papers.Independent.Intakethird100.P52.renyiKernel
For finite, possibly empty types with decidable equality, any , arbitrary complex matrices on , and any function , define by continuous functional calculus with the real function that is at and otherwise. Define entrywise by , where has entry at and zero elsewhere. The declaration defines the matrix
There are no positivity, normalization, channel, or linearity assumptions. Every real , including , is included; uses the function at zero and elsewhere, rather than the constant function .
QAlgorithms.Papers.Independent.Intakethird100.P52.renyiReference
For finite, possibly empty types with decidable equality, any , arbitrary complex matrices on , and any function , let denote continuous functional calculus for the function equal to at zero and otherwise, and define using matrix units . This declaration defines
There are no positivity, normalization, channel, or linearity assumptions. In particular is allowed: the outer exponent is then zero under the convention , and uses the function at zero and elsewhere.
QAlgorithms.Papers.Independent.Intakethird100.P52.sandwichedKernel
For finite, possibly empty types with decidable equality, any , arbitrary complex matrices on , and any function , let be continuous functional calculus for the real function equal to at zero and elsewhere, and put , with the matrix unit. The declaration defines
No positivity, normalization, channel, or linearity assumptions are imposed. The definition includes , where each displayed exponent is zero because real division by zero is defined as zero; uses the function at zero and 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 with decidable equality, any , arbitrary complex matrices on , and any function , let denote continuous functional calculus for the function equal to at zero and otherwise, and set using matrix units . The defined matrix is
No positivity, normalization, channel, or linearity assumptions are imposed. Both and are included, with every division by zero evaluated as zero; an exponent zero invokes the function at zero and elsewhere.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.