Section 3.3: purification and equality of complementary entropies
OpenQAlgorithms.Papers.Independent.Intakethird100.P51.purificationEvery finite normalized positive semidefinite state has a finite nonzero auxiliary system and a unit vector such that
Moreover, for every bipartite unit vector, its two reduced states have equal entropy. Singular states are permitted. The existence of the purification is part of the conclusion and supplies the complementary-entropy method.
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]
/-- Section 3.3, PDF p.9: purification exists; complementary entropies agree for
all unit bipartite pure states, including singular density matrices. -/
theorem purification [Nonempty A] [Nonempty B] (rho : Mat A) (h : IsDensityMatrix rho) :
(∃ m : ℕ, 0 < m ∧ ∃ psi : EuclideanSpace ℂ (A × Fin m),
‖psi‖ = 1 ∧ traceRight (pureDensity psi) = rho ∧
entropy (traceLeft (pureDensity psi)) = entropy rho) ∧
(∀ psi : EuclideanSpace ℂ (A × B), ‖psi‖ = 1 →
entropy (traceRight (pureDensity psi)) = entropy (traceLeft (pureDensity psi))) := 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 finite nonempty index types with decidable equality and every positive semidefinite complex matrix on with trace , there exist a natural number and a vector of Euclidean norm such that and . In addition, every vector of Euclidean norm satisfies . Here has entries , partial traces sum over the discarded index, and , with applied by continuous functional calculus using natural logarithms and . Singular density matrices and registers of cardinality are included.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.