Conditions EQ-MEAS, EQ-LIND, EQ-HS, EQ-KL on a black-box process
Definitionpme_measurement_conditionsThe four characterizations of the spectral equilibration process associated with a Hermitian observable on , as properties of a map on matrices (only its values on density matrices matter).
- EQ-MEAS (nonselective projective measurement): for every density matrix .
- EQ-LIND (spectral equilibration): for every density matrix , the Lindblad equation with has a solution on , and for every such solution.
- EQ-HS (closest commuting ensemble): for every density matrix , is a density matrix with attaining .
- EQ-KL (closest commuting ensemble in information): the same with the relative entropy in place of .
Here is encoded as " is a minimizer"; uniqueness of the minimizer is part of what Theorem 1 proves.
import Mathlib
import Definitions.Def_pme_quantum_basics
/-!
The four characterizations EQ-MEAS, EQ-LIND, EQ-HS, EQ-KL of the black-box spectral
equilibration process `Φ_X` (Carcassi, AoP brief 017, Section 2), as predicates on a map
`Φ : Matrix n n ℂ → Matrix n n ℂ` (only its values on density matrices matter).
-/
namespace ProjectiveMeasurementEquilibration
open Matrix Filter Topology
variable {n : Type} [Fintype n] [DecidableEq n]
/-- Condition EQ-MEAS: `Φ(ρ) = ∑ₓ Pₓ ρ Pₓ` for every density matrix `ρ`. -/
def EqMeas {X : Matrix n n ℂ} (hX : X.IsHermitian) (Φ : Matrix n n ℂ → Matrix n n ℂ) : Prop :=
∀ ρ, IsDensityMatrix ρ → Φ ρ = pinching hX ρ
/-- Condition EQ-LIND: for every density matrix `ρ`, the Lindblad equation
`dρ/dt = XρX − ½{X², ρ}` started at `ρ` has a solution on `t ≥ 0`, and every such solution
converges to `Φ(ρ)` as `t → ∞`. -/
def EqLind (X : Matrix n n ℂ) (Φ : Matrix n n ℂ → Matrix n n ℂ) : Prop :=
∀ ρ, IsDensityMatrix ρ →
(∃ ρt : ℝ → Matrix n n ℂ, IsLindbladSolution X ρ ρt) ∧
∀ ρt : ℝ → Matrix n n ℂ, IsLindbladSolution X ρ ρt → Tendsto ρt atTop (𝓝 (Φ ρ))
/-- Condition EQ-HS: `Φ(ρ)` is a density matrix commuting with `X` that minimizes
`‖ρ − σ‖²_HS` over all density matrices `σ` commuting with `X`. -/
def EqHS (X : Matrix n n ℂ) (Φ : Matrix n n ℂ → Matrix n n ℂ) : Prop :=
∀ ρ, IsDensityMatrix ρ →
IsDensityMatrix (Φ ρ) ∧ X * Φ ρ = Φ ρ * X ∧
∀ σ, IsDensityMatrix σ → X * σ = σ * X → hsDistSq ρ (Φ ρ) ≤ hsDistSq ρ σ
/-- Condition EQ-KL: `Φ(ρ)` is a density matrix commuting with `X` that minimizes
`D(ρ‖σ)` over all density matrices `σ` commuting with `X`. -/
def EqKL (X : Matrix n n ℂ) (Φ : Matrix n n ℂ → Matrix n n ℂ) : Prop :=
∀ ρ, IsDensityMatrix ρ →
IsDensityMatrix (Φ ρ) ∧ X * Φ ρ = Φ ρ * X ∧
∀ σ, IsDensityMatrix σ → X * σ = σ * X → relEntropy ρ (Φ ρ) ≤ relEntropy ρ σ
end ProjectiveMeasurementEquilibration
Read-back
What the Lean code literally says, in plain math · Aristotle (Harmonic) - same agent as the drafter; non-blind
Non-blind read-back. This read-back was written by the same agent (Aristotle, by Harmonic) that drafted the Lean statements, at the explicit request of the account owner. It is not blind and not independent testimony: the author knew the intended meaning while writing it. Reviewers should not treat it as an independent audit.
For a finite index type , a matrix and an arbitrary function , four propositions are defined (all using the notions of the companion basics file). In each, and range over all density matrices (positive semidefinite, trace ).
-
(takes a proof that is Hermitian): for every density matrix , , the sum over the distinct eigenvalues of , the eigenprojections.
-
: for every density matrix , (a) there exists a function with and right-sided/within- derivative for all ; and (b) every such function satisfies entrywise as .
-
: for every density matrix : is a density matrix, , and for every density matrix with , .
-
: the same as 3 but with the inequality in the extended reals, where when and otherwise ( by functional calculus with ).
Only refers to the Hermitian proof; the others take alone. None constrains outside density matrices. If is empty all four hold vacuously.
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.