Theorem 2: Lieb concavity of trace exp(K+log A)
OpenQAlgorithms.Papers.Independent.Intakethird100.P51.liebConcavityFix a self-adjoint matrix and write for positive definite . For positive definite and ,
The matrices need not be normalized. This is non-strict concavity on the entire positive definite cone.
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P51
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)
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]
/-- Theorem 2, PDF p.5: non-strict Lieb concavity for arbitrary positive definite matrices. -/
theorem liebConcavity [Nonempty A] (L X Y : Mat A) (hL : L.IsHermitian) (hX : X.PosDef) (hY : Y.PosDef)
(t : ℝ) (ht0 : 0 ≤ t) (ht1 : t ≤ 1) :
t * (NormedSpace.exp (L + logMatrix X)).trace.re +
(1 - t) * (NormedSpace.exp (L + logMatrix Y)).trace.re ≤
(NormedSpace.exp (L + logMatrix ((t : ℂ) • X + ((1 - t : ℝ) : ℂ) • Y))).trace.re := by sorry
end Statements
end QAlgorithms.Papers.Independent.Intakethird100.P51
end
Read-back
What the Lean code literally says, in plain math · gpt-6.1-sol
For every nonempty finite index type with decidable equality, every Hermitian complex matrix indexed by , every pair of positive definite matrices indexed by , and every real satisfying , the theorem asserts
Here means the continuous-functional-calculus application of the natural logarithm to , and is the matrix exponential. No trace normalization is imposed on , no commutation hypothesis is imposed, and both endpoints and one-dimensional index types are included.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.