Definitions for Pure-state transformations and catalysis under operations that completely preserve positivity of partial transpose
DefinitionQAlgorithms_Papers_Independent_Intakethird100_P95The objects the formalized statements of Pure-state transformations and catalysis under operations that completely preserve positivity of partial transpose (arXiv:0801.4322) are stated in. Each was checked against the paper by a blind read-back and a moderator.
schmidtRenyi: PDF p.2, equation (5): finite orders, with separate support and Shannon branches.largestProbability: PDF p.2, equation (5): maximum Schmidt probability. Its range is finite; normalization ensures it is nonempty and has a strictly positive maximum.renyiInfinity: PDF p.2, equation (5): the infinity order is a separate branch.schmidtDensity: PDF p.2, equation (1): the normalized Schmidt-form density when lam is a probability vector.bobTransposeLinear: PDF p.2, equation (4): complex-linear packaging of the accepted Bob partial transpose. The entry function is reused, rather than copied.PPTChannel: PDF p.2: CPTP and complete positivity of Gamma_out L Gamma_in. The accepted IsCompletelyPositive quantifies over every finite ancilla.Reach: PDF p.2: deterministic exact operational reachability, independently of all SDPs.tensorSpectrum: PDF p.3: tensor Schmidt probabilities. The empty tensor is the rank-one product state.rhoPower: PDF p.3: n copies, with all of Alice's copies grouped opposite all of Bob's.ValidLengths: PDF p.3, Definitions 1–2: positive block lengths tending to infinity.RateLimit: PDF p.3, Definitions 1–2: normalized log resource rank, independent of the protocols.CostRates: PDF p.3, Definition 2: achievable formation rates. Exact means equality at EVERY block; otherwise the trace norm error tends to zero.DistillationRates: PDF p.3, Definition 1: achievable distillation rates with zero error at every block in the exact case, and vanishing trace norm error otherwise.pptEntanglementCost: PDF p.3, Definition 2: operational infimum, retaining infinite values.pptDistillableEntanglement: PDF p.3, Definition 1: operational supremum, retaining infinite values.PreparationFeasible: PDF p.5, equation (14): one common Hermitian positive variable P in BOTH inequalities.preparationValues: PDF p.5, equation (14): objective set of the general mixed-state preparation SDP.preparationValue: PDF p.5, equation (14): infimum defined from the SDP, independently of reachability. Lemma 8 additionally asserts a minimum, excluding empty/unbounded default infima.SymIndex: PDF p.5: symmetric coordinates, one for each j ≤ i.AntiIndex: PDF p.5: antisymmetric coordinates, one for each j < i.symmetricSchmidtProjector: PDF p.5: sigma_ii is |ii><ii|; off diagonal symmetric projectors include 1/2.antisymmetricSchmidtProjector: PDF p.5: antisymmetric projectors include 1/2.coordinateMatrix: PDF pp.5–6, equations (20),(29): partially transposed coordinate matrix.SchmidtPrimalFeasible: PDF p.5, equation (20): the primal constraints.schmidtPrimalObjective: PDF p.5, equation (20): each real coordinate is counted once.primalValues: PDF p.5, equation (20): the set of feasible objective values.primalValue: PDF p.5, equation (20): T is defined independently as an SDP infimum; Proposition 9 records its attainment.SchmidtDualFeasible: PDF p.6, equation (29): upper bounds only, with no nonnegativity constraints on coordinates.schmidtDualObjective: PDF p.6, equation (28): dual objective with both weighted sums subtracted.schmidtDualValues: PDF p.6, equation (28): the dual's objective set.rankOneValues: PDF p.7, equation (36): rank EXACTLY one, without the spurious zero witness in (A3).rankOneValue: PDF p.7, Remark 14: supremum version of the restricted value; no assertion of attainment. Its nonemptiness and upper boundedness are conclusions of Remark 14 below.IsUniformRank: PDF pp.2,7,9: equality to Phi_K up to local unitaries and zero padding.CatalyticPreparation: PDF p.9, Theorem 21: exact, uncorrelated return of a finite maximally entangled catalyst.faceQuadratic: PDF p.11, (A4): the actual quadratic objective, with all diagonal terms in the second sum.OnFace: PDF pp.11–12: a point on the CLOSED face F_C; negative free coordinates are allowed.FaceStationary: PDF p.12, (A6): stationary for the restriction to the affine hull H_C and lying on F_C. The free-coordinate derivative of faceQuadratic is 2 / (K^2-1) times (sum_i x_i - K*x_j); this is not a constrained local extremum.faceValue: PDF p.12, (A7): stationary value; cast BEFORE subtraction. Under Lemma 23's hypotheses the denominator cannot vanish. Negative denominators are retained.exampleSpectrum: PDF p.4, Example 7: the displayed multiset, independent of its erroneous down-arrow label.
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Information
import Definitions.Def_QAlgorithms_Defs_MatrixMapPositivity
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P28
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P27
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P63
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P19
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P04
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P18
import Definitions.Def_QAlgorithms_Papers_Independent_Support50_P11
/-!
Matthews and Winter, arXiv:0801.4322v2. All page numbers are PDF pages.
Paper-specific finite-dimensional definitions are local to this paper; exact
accepted primitives are imported without local copies. They do not implement the
arbitrary-separable-space interfaces requested by the separate 2007.14270 paper.
Theorems 13 and 20 and Conjectures 19 and 22 are deferred in HARD.md.
-/
noncomputable section
open scoped BigOperators Matrix ComplexOrder
open Filter
namespace QAlgorithms.Papers.Independent.Intakethird100.P95
attribute [local instance] Classical.propDecidable Classical.decEq
-- Exact primitives from the pinned accepted library; no local aliases.
open QAlgorithms.Papers.Independent.Claimed50.P28 (maximallyEntangledDensity)
open QAlgorithms.Papers.Independent.Information100.P27 (IsProbability)
open QAlgorithms.Papers.Independent.Information100.P63 (logTwo)
open QAlgorithms.Papers.Independent.Information100.P19 (outer)
open QAlgorithms.Papers.Independent.Dependency50.P18 (bipartiteTensor)
open QAlgorithms.Papers.Independent.Support50.P11 (Loewner)
/-- PDF p.2, equation (5): finite orders, with separate support and Shannon branches. -/
def schmidtRenyi {I : Type*} [Fintype I] (t : ℝ) (lam : I → ℝ) : ℝ :=
if t = 0 then logTwo ((Finset.univ.filter fun i => lam i ≠ 0).card : ℝ)
else if t = 1 then - ∑ i, lam i * logTwo (lam i)
else logTwo (∑ i, Real.rpow (lam i) t) / (1 - t)
/-- PDF p.2, equation (5): maximum Schmidt probability. Its range is finite;
normalization ensures it is nonempty and has a strictly positive maximum. -/
def largestProbability {I : Type*} [Fintype I] (lam : I → ℝ) : ℝ :=
sSup (Set.range lam)
/-- PDF p.2, equation (5): the infinity order is a separate branch. -/
def renyiInfinity {I : Type*} [Fintype I] (lam : I → ℝ) : ℝ :=
- logTwo (largestProbability lam)
/-- PDF p.2, equation (1): the normalized Schmidt-form density when lam is a probability vector. -/
def schmidtDensity {I : Type*} (lam : I → ℝ) : Matrix (I × I) (I × I) ℂ :=
fun p q => if p.1 = p.2 ∧ q.1 = q.2 then (Real.sqrt (lam p.1 * lam q.1) : ℂ) else 0
/-- PDF p.2, equation (4): complex-linear packaging of the accepted Bob partial
transpose. The entry function is reused, rather than copied. -/
def bobTransposeLinear {A B : Type*} :
Matrix (A × B) (A × B) ℂ →ₗ[ℂ] Matrix (A × B) (A × B) ℂ where
toFun := QAlgorithms.Papers.Independent.Information100.P04.gamma
map_add' _ _ := rfl
map_smul' _ _ := rfl
/-- PDF p.2: CPTP and complete positivity of Gamma_out L Gamma_in.
The accepted IsCompletelyPositive quantifies over every finite ancilla. -/
structure PPTChannel (A B C D : Type*) [Fintype A] [Fintype B]
[Fintype C] [Fintype D] where
map : Matrix (A × B) (A × B) ℂ →ₗ[ℂ] Matrix (C × D) (C × D) ℂ
completelyPositive : QAlgorithms.Entanglement.IsCompletelyPositive map
tracePreserving : ∀ X, (map X).trace = X.trace
transposeCompletelyPositive :
QAlgorithms.Entanglement.IsCompletelyPositive (bobTransposeLinear.comp (map.comp bobTransposeLinear))
/-- PDF p.2: deterministic exact operational reachability, independently of all SDPs. -/
def Reach {A B C D : Type*} [Fintype A] [Fintype B] [Fintype C] [Fintype D]
(X : Matrix (A × B) (A × B) ℂ) (Y : Matrix (C × D) (C × D) ℂ) : Prop :=
∃ L : PPTChannel A B C D, L.map X = Y
/-- PDF p.3: tensor Schmidt probabilities. The empty tensor is the rank-one product state. -/
def tensorSpectrum {I : Type*} (lam : I → ℝ) (n : ℕ) : (Fin n → I) → ℝ :=
fun x => ∏ j, lam (x j)
/-- PDF p.3: n copies, with all of Alice's copies grouped opposite all of Bob's. -/
def rhoPower {I : Type*} (lam : I → ℝ) (n : ℕ) := schmidtDensity (tensorSpectrum lam n)
/-- PDF p.3, Definitions 1–2: positive block lengths tending to infinity. -/
def ValidLengths (n K : ℕ → ℕ) : Prop :=
(∀ j, 0 < n j) ∧ Tendsto n atTop atTop ∧ (∀ j, 0 < K j)
/-- PDF p.3, Definitions 1–2: normalized log resource rank, independent of the protocols. -/
def RateLimit (n K : ℕ → ℕ) (R : ℝ) : Prop :=
Tendsto (fun j => logTwo (K j : ℝ) / (n j : ℝ)) atTop (nhds R)
/-- PDF p.3, Definition 2: achievable formation rates. Exact means equality at
EVERY block; otherwise the trace norm error tends to zero. -/
def CostRates {I : Type*} [Fintype I] (lam : I → ℝ) (zeroError : Bool) : Set EReal :=
{r | ∃ R : ℝ, r = (R : EReal) ∧ ∃ n K : ℕ → ℕ,
ValidLengths n K ∧ RateLimit n K R ∧
∃ L : (j : ℕ) → PPTChannel (Fin (K j)) (Fin (K j)) (Fin (n j) → I) (Fin (n j) → I),
if zeroError then ∀ j, (L j).map (maximallyEntangledDensity (K j)) = rhoPower lam (n j)
else Tendsto (fun j => QAlgorithms.traceNorm
((L j).map (maximallyEntangledDensity (K j)) - rhoPower lam (n j))) atTop (nhds 0)}
/-- PDF p.3, Definition 1: achievable distillation rates with zero error at every
block in the exact case, and vanishing trace norm error otherwise. -/
def DistillationRates {I : Type*} [Fintype I] (lam : I → ℝ) (zeroError : Bool) : Set EReal :=
{r | ∃ R : ℝ, r = (R : EReal) ∧ ∃ n K : ℕ → ℕ,
ValidLengths n K ∧ RateLimit n K R ∧
∃ L : (j : ℕ) → PPTChannel (Fin (n j) → I) (Fin (n j) → I) (Fin (K j)) (Fin (K j)),
if zeroError then ∀ j, (L j).map (rhoPower lam (n j)) = maximallyEntangledDensity (K j)
else Tendsto (fun j => QAlgorithms.traceNorm
(maximallyEntangledDensity (K j) - (L j).map (rhoPower lam (n j)))) atTop (nhds 0)}
/-- PDF p.3, Definition 2: operational infimum, retaining infinite values. -/
def pptEntanglementCost {I : Type*} [Fintype I] (lam : I → ℝ) (zeroError : Bool) : EReal :=
sInf (CostRates lam zeroError)
/-- PDF p.3, Definition 1: operational supremum, retaining infinite values. -/
def pptDistillableEntanglement {I : Type*} [Fintype I] (lam : I → ℝ) (zeroError : Bool) : EReal :=
sSup (DistillationRates lam zeroError)
/-- PDF p.5, equation (14): one common Hermitian positive variable P in BOTH inequalities. -/
def PreparationFeasible {A B : Type*} [Fintype A] [Fintype B] (K : ℕ)
(X P : Matrix (A × B) (A × B) ℂ) : Prop :=
P.IsHermitian ∧ P.PosSemidef ∧
Loewner ((-((K : ℝ) - 1) : ℂ) • bobTransposeLinear P) (bobTransposeLinear X) ∧
Loewner (bobTransposeLinear X) (((K : ℝ) + 1 : ℂ) • bobTransposeLinear P)
/-- PDF p.5, equation (14): objective set of the general mixed-state preparation SDP. -/
def preparationValues {A B : Type*} [Fintype A] [Fintype B]
(K : ℕ) (X : Matrix (A × B) (A × B) ℂ) : Set ℝ :=
{v | ∃ P, PreparationFeasible K X P ∧ P.trace.re = v}
/-- PDF p.5, equation (14): infimum defined from the SDP, independently of reachability.
Lemma 8 additionally asserts a minimum, excluding empty/unbounded default infima. -/
def preparationValue {A B : Type*} [Fintype A] [Fintype B]
(K : ℕ) (X : Matrix (A × B) (A × B) ℂ) : ℝ := sInf (preparationValues K X)
/-- PDF p.5: symmetric coordinates, one for each j ≤ i. -/
abbrev SymIndex (d : ℕ) := {p : Fin d × Fin d // p.2 ≤ p.1}
/-- PDF p.5: antisymmetric coordinates, one for each j < i. -/
abbrev AntiIndex (d : ℕ) := {p : Fin d × Fin d // p.2 < p.1}
/-- PDF p.5: sigma_ii is |ii><ii|; off diagonal symmetric projectors include 1/2. -/
def symmetricSchmidtProjector {d : ℕ} (i j : Fin d) : Matrix (Fin d × Fin d) (Fin d × Fin d) ℂ :=
if i = j then outer (WithLp.ofLp (QAlgorithms.ket (i, i)))
else (1 / 2 : ℂ) • outer (WithLp.ofLp (QAlgorithms.ket (i, j)) + WithLp.ofLp (QAlgorithms.ket (j, i)))
/-- PDF p.5: antisymmetric projectors include 1/2. -/
def antisymmetricSchmidtProjector {d : ℕ} (i j : Fin d) : Matrix (Fin d × Fin d) (Fin d × Fin d) ℂ :=
(1 / 2 : ℂ) • outer (WithLp.ofLp (QAlgorithms.ket (i, j)) - WithLp.ofLp (QAlgorithms.ket (j, i)))
/-- PDF pp.5–6, equations (20),(29): partially transposed coordinate matrix. -/
def coordinateMatrix {d : ℕ} (s : SymIndex d → ℝ) (a : AntiIndex d → ℝ) :
Matrix (Fin d × Fin d) (Fin d × Fin d) ℂ :=
(∑ p, (s p : ℂ) • bobTransposeLinear (symmetricSchmidtProjector p.val.1 p.val.2)) +
∑ p, (a p : ℂ) • bobTransposeLinear (antisymmetricSchmidtProjector p.val.1 p.val.2)
/-- PDF p.5, equation (20): the primal constraints. -/
def SchmidtPrimalFeasible {d : ℕ} (K : ℕ) (lam : Fin d → ℝ)
(s : SymIndex d → ℝ) (a : AntiIndex d → ℝ) : Prop :=
(∀ p, Real.sqrt (lam p.val.1 * lam p.val.2) / ((K : ℝ) + 1) ≤ s p) ∧
(∀ p, Real.sqrt (lam p.val.1 * lam p.val.2) / ((K : ℝ) - 1) ≤ a p) ∧
(coordinateMatrix s a).PosSemidef
/-- PDF p.5, equation (20): each real coordinate is counted once. -/
def schmidtPrimalObjective {d : ℕ} (s : SymIndex d → ℝ) (a : AntiIndex d → ℝ) : ℝ :=
(∑ p, s p) + ∑ p, a p
/-- PDF p.5, equation (20): the set of feasible objective values. -/
def primalValues {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) : Set ℝ :=
{v | ∃ s a, SchmidtPrimalFeasible K lam s a ∧ schmidtPrimalObjective s a = v}
/-- PDF p.5, equation (20): T is defined independently as an SDP infimum;
Proposition 9 records its attainment. -/
def primalValue {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) : ℝ := sInf (primalValues K lam)
/-- PDF p.6, equation (29): upper bounds only, with no nonnegativity constraints on coordinates. -/
def SchmidtDualFeasible {d : ℕ} (mu : SymIndex d → ℝ) (t : AntiIndex d → ℝ) : Prop :=
(∀ p, mu p ≤ 1) ∧ (∀ p, t p ≤ 1) ∧ (coordinateMatrix mu t).PosSemidef
/-- PDF p.6, equation (28): dual objective with both weighted sums subtracted. -/
def schmidtDualObjective {d : ℕ} (K : ℕ) (lam : Fin d → ℝ)
(mu : SymIndex d → ℝ) (t : AntiIndex d → ℝ) : ℝ :=
((K : ℝ) * Real.rpow 2 (schmidtRenyi (1 / 2) lam) - 1) / ((K : ℝ)^2 - 1) -
((∑ p, mu p * Real.sqrt (lam p.val.1 * lam p.val.2) / ((K : ℝ) + 1)) +
∑ p, t p * Real.sqrt (lam p.val.1 * lam p.val.2) / ((K : ℝ) - 1))
/-- PDF p.6, equation (28): the dual's objective set. -/
def schmidtDualValues {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) : Set ℝ :=
{v | ∃ mu t, SchmidtDualFeasible mu t ∧ schmidtDualObjective K lam mu t = v}
/-- PDF p.7, equation (36): rank EXACTLY one, without the spurious zero witness in (A3). -/
def rankOneValues {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) : Set ℝ :=
{v | ∃ mu t, SchmidtDualFeasible mu t ∧ (coordinateMatrix mu t).rank = 1 ∧
schmidtDualObjective K lam mu t = v}
/-- PDF p.7, Remark 14: supremum version of the restricted value; no assertion
of attainment. Its nonemptiness and upper boundedness are conclusions of Remark 14 below. -/
def rankOneValue {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) : ℝ := sSup (rankOneValues K lam)
/-- PDF pp.2,7,9: equality to Phi_K up to local unitaries and zero padding. -/
def IsUniformRank {I : Type*} [Fintype I] (lam : I → ℝ) (K : ℕ) : Prop :=
(Finset.univ.filter fun i => lam i ≠ 0).card = K ∧
∀ i, lam i = 0 ∨ lam i = (K : ℝ)⁻¹
/-- PDF p.9, Theorem 21: exact, uncorrelated return of a finite maximally entangled catalyst. -/
def CatalyticPreparation {I : Type} [Fintype I] (K : ℕ) (lam : I → ℝ) : Prop :=
∃ C : ℕ, 2 ≤ C ∧ Reach (bipartiteTensor (maximallyEntangledDensity C) (maximallyEntangledDensity K))
(bipartiteTensor (maximallyEntangledDensity C) (schmidtDensity lam))
/-- PDF p.11, (A4): the actual quadratic objective, with all diagonal terms in the second sum. -/
def faceQuadratic {d : ℕ} (K : ℕ) (x : Fin d → ℝ) : ℝ :=
(∑ p : AntiIndex d, x p.val.1 * x p.val.2 / ((K : ℝ) - 1)) -
∑ p : SymIndex d, x p.val.1 * x p.val.2 / ((K : ℝ) + 1)
/-- PDF pp.11–12: a point on the CLOSED face F_C; negative free coordinates are allowed. -/
def OnFace {d : ℕ} (lam : Fin d → ℝ) (C : Finset (Fin d)) (x : Fin d → ℝ) : Prop :=
(∀ i, |x i| ≤ Real.sqrt (lam i)) ∧ (∀ i ∈ C, x i = Real.sqrt (lam i))
/-- PDF p.12, (A6): stationary for the restriction to the affine hull H_C and
lying on F_C. The free-coordinate derivative of faceQuadratic is
2 / (K^2-1) times (sum_i x_i - K*x_j); this is not a constrained local extremum. -/
def FaceStationary {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) (C : Finset (Fin d)) : Prop :=
∃ x : Fin d → ℝ, OnFace lam C x ∧
∀ j, j ∉ C → (∑ i, x i) - (K : ℝ) * x j = 0
/-- PDF p.12, (A7): stationary value; cast BEFORE subtraction. Under Lemma 23's
hypotheses the denominator cannot vanish. Negative denominators are retained. -/
def faceValue {d : ℕ} (K : ℕ) (lam : Fin d → ℝ) (C : Finset (Fin d)) : ℝ :=
let D : ℝ := (K : ℝ) + (C.card : ℝ) - (d : ℝ)
(K : ℝ) / (((K : ℝ)^2 - 1) * D) *
((∑ i ∈ C, Real.sqrt (lam i))^2 - D * ∑ i ∈ C, lam i)
/-- PDF p.4, Example 7: the displayed multiset, independent of its erroneous down-arrow label. -/
def exampleSpectrum : Fin 6 → ℝ := ![1/20, 1/20, 1/20, 4/20, 4/20, 9/20]
/-! Sanity examples for the complete local definitions (not theorem proofs). -/
example {A B : Type*} (X : Matrix (A × B) (A × B) ℂ) : bobTransposeLinear (bobTransposeLinear X) = X := rfl
example : bobTransposeLinear (fun p q : Fin 2 × Fin 2 =>
if p = (0, 1) ∧ q = (1, 0) then (1 : ℂ) else 0) (0, 0) (1, 1) = 1 := by
change (if ((0, 1) : Fin 2 × Fin 2) = (0, 1) ∧
((1, 0) : Fin 2 × Fin 2) = (1, 0) then (1 : ℂ) else 0) = 1
simp
example : schmidtDensity (fun _ : Fin 1 => (1 : ℝ)) (0, 0) (0, 0) = 1 := by
norm_num [schmidtDensity]
example {I : Type*} (lam : I → ℝ) (x : Fin 0 → I) : tensorSpectrum lam 0 x = 1 := by
simp [tensorSpectrum]
example : maximallyEntangledDensity 2 (0, 0) (1, 1) = (1 / 2 : ℂ) := by
norm_num [maximallyEntangledDensity]
example {A B C D : Type} (X : Matrix (A × B) (A × B) ℂ)
(Y : Matrix (C × D) (C × D) ℂ) (a a' : A) (b b' : B) (c c' : C) (d d' : D) :
bipartiteTensor X Y ((a,c),(b,d)) ((a',c'),(b',d')) =
X (a,b) (a',b') * Y (c,d) (c',d') := rfl
example : symmetricSchmidtProjector (0 : Fin 2) 1 (0, 1) (1, 0) = (1 / 2 : ℂ) := by
norm_num [symmetricSchmidtProjector, outer, QAlgorithms.ket, PiLp.ofLp_single, Pi.single_apply]
example : antisymmetricSchmidtProjector (0 : Fin 2) 1 (0, 1) (1, 0) = (-1 / 2 : ℂ) := by
norm_num [antisymmetricSchmidtProjector, outer, QAlgorithms.ket, PiLp.ofLp_single, Pi.single_apply]
example : schmidtRenyi 1 (fun _ : Fin 1 => (1 : ℝ)) = 0 := by
simp [schmidtRenyi, logTwo]
example : (∑ i, exampleSpectrum i) = 1 := by
norm_num [exampleSpectrum, Fin.sum_univ_succ]
example {A B : Type*} [Fintype A] [Fintype B] (X P : Matrix (A × B) (A × B) ℂ) :
PreparationFeasible 2 X P ↔ P.IsHermitian ∧ P.PosSemidef ∧
Loewner ((-1 : ℂ) • bobTransposeLinear P) (bobTransposeLinear X) ∧
Loewner (bobTransposeLinear X) ((3 : ℂ) • bobTransposeLinear P) := by
norm_num [PreparationFeasible]
end QAlgorithms.Papers.Independent.Intakethird100.P95
Read-back
What the Lean code literally says, in plain math · gpt-6.1-sol
QAlgorithms.Papers.Independent.Intakethird100.P95.AntiIndex
For each natural number , this abbreviates the type of ordered pairs of integers with and . Diagonal pairs are excluded, and the type is empty when or .
QAlgorithms.Papers.Independent.Intakethird100.P95.CatalyticPreparation
For any finite type , any natural number , and any function , this predicate means that there exist a natural number and a complex-linear map from matrices on to matrices on such that is completely positive, preserves the trace of every matrix, its conjugation by partial transpose on the second subsystem is completely positive, and
The products here group the two first subsystems together and the two second subsystems together. Explicitly, if and , and is zero otherwise; if and , and is zero otherwise. Partial transpose sends to ; complete positivity requires every finite ancillary block amplification, including size zero, to preserve positive semidefiniteness. No nonnegativity or normalization of is assumed. Empty and are included, and negative real square-root arguments are evaluated as zero.
QAlgorithms.Papers.Independent.Intakethird100.P95.CostRates
For a finite type , an arbitrary function , and a Boolean flag , this defines a set of extended real numbers containing precisely those for which there exist a finite real number , sequences , and maps from matrices on to matrices on , such that , every and is positive, , and . Each must be complex-linear and trace preserving on every matrix, and both and its composite with partial transposition on the second input and second output factors must be completely positive, meaning positivity preserving under every finite identity amplification, including dimension . Define
and
If the flag is true, the requirement is for every ; if it is false, the requirement is , where . No nonnegativity or normalization assumption is imposed on , no monotonicity of the sequences is required, and neither infinite extended-real values nor an explicit restriction appears among the witnesses. If is empty, the set is empty because positive would make the output matrix space have identically zero trace, preventing trace preservation from the nonempty input space.
QAlgorithms.Papers.Independent.Intakethird100.P95.DistillationRates
For any finite type , any function , and any Boolean , this defines a set of extended real numbers consisting precisely of the finite numbers for which there are sequences of positive integers , with and , and maps from matrices on to matrices on , satisfying the following requirements. Each is complex linear, trace preserving on every matrix, and completely positive, and its conjugation by partial transpose on the second subsystem is also completely positive; complete positivity means preservation of positive semidefiniteness after adjoining an identity map of every finite dimension, including dimension zero. Write , and define by
Define to have entry at when and , and zero otherwise. If is true, for every ; if is false, , where , with the square root interpreted through functional calculus. No nonnegativity or normalization is assumed for , and real square roots of negative numbers are zero. The type may be empty, although all and must be positive; neither nor belongs to this set.
QAlgorithms.Papers.Independent.Intakethird100.P95.FaceStationary
For natural numbers , any real vector , and any subset , this predicate means that there exists a real vector such that for every , for every , and
There are no positivity or normalization assumptions on and no restriction on . Real square roots of negative inputs are defined to be zero. If is the entire index set, the final equations impose no conditions; if , the predicate holds using the empty vector.
QAlgorithms.Papers.Independent.Intakethird100.P95.IsUniformRank
For any finite type , equipped with a finite enumeration, any function , and any natural number , this predicate means
For , the inverse is zero and the predicate holds exactly when is identically zero. For , exactly entries equal and every other entry is zero. Empty is allowed and then the predicate holds exactly for .
QAlgorithms.Papers.Independent.Intakethird100.P95.OnFace
For every natural number , arbitrary real vectors , and any subset , this predicate means
No nonnegativity or normalization of is assumed. The real square root of a negative number is zero, so forces . If is empty, the equality requirements are vacuous; if , the entire predicate holds vacuously.
QAlgorithms.Papers.Independent.Intakethird100.P95.PPTChannel
For finite types , this structure consists of a complex-linear map from matrices on to matrices on , together with proofs that is completely positive, preserves the complex trace of every matrix, and that is completely positive. Here , and is defined analogously on the output. Complete positivity means that, for every natural number , including , the amplification acting on each matrix block by the given map sends every positive semidefinite matrix on to a positive semidefinite matrix on ; the same condition applies to the transposed composite. The four finite types may be empty.
QAlgorithms.Papers.Independent.Intakethird100.P95.PreparationFeasible
For finite types , possibly empty, a natural number , and arbitrary complex matrices indexed by , this predicate requires to be Hermitian and positive semidefinite and requires
where , and means positive semidefinite. The subtraction is performed in the real numbers: gives the two conditions and , while gives and . No trace normalization or separate hypothesis on is included.
QAlgorithms.Papers.Independent.Intakethird100.P95.RateLimit
For sequences and a real number , this predicate means
This predicate itself imposes no positivity, divergence, or monotonicity conditions on either sequence. The total-function conventions assign and division by zero the value , so terms with or are included and have value .
QAlgorithms.Papers.Independent.Intakethird100.P95.Reach
For finite types and arbitrary complex matrices on and on , this predicate means that there exists a complex-linear map with which preserves the trace of every matrix and for which both and are completely positive. Here exchanges the two indices of a matrix entry, and exchanges the two indices. Complete positivity requires positivity preservation after amplification by the identity on for every . Neither nor is assumed to be positive semidefinite or normalized, and empty indexing types are allowed.
QAlgorithms.Papers.Independent.Intakethird100.P95.SchmidtDualFeasible
For every natural number , and real arrays , , and , , this predicate means that every , every , and the following matrix on the standard basis of is positive semidefinite:
The coordinate bounds contain no explicit lower bounds. For , the arrays and matrix are empty and the predicate holds.
QAlgorithms.Papers.Independent.Intakethird100.P95.SchmidtPrimalFeasible
For any natural numbers , arbitrary real function , real coefficients for , and real coefficients for , this predicate requires
and requires the following matrix to be positive semidefinite:
Here is the coordinate vector at , and has entries . No nonnegativity or normalization assumption is imposed on , and a negative real square-root argument gives zero. The denominators use real subtraction and total division: at , every antisymmetric lower bound is zero; at , it is . When , the coefficient conditions are vacuous and the empty matrix is positive semidefinite.
QAlgorithms.Papers.Independent.Intakethird100.P95.SymIndex
For each natural number , this abbreviates the type of ordered pairs of integers with and . Diagonal pairs are included, and the type is empty when .
QAlgorithms.Papers.Independent.Intakethird100.P95.ValidLengths
For two sequences , this predicate means that and for every , and that as : for every natural bound, all sufficiently late meet or exceed that bound. Neither sequence is required to be monotone, and no growth condition is imposed on .
QAlgorithms.Papers.Independent.Intakethird100.P95.antisymmetricSchmidtProjector
For any natural number and indices , this defines
where is the coordinate vector equal to one at and zero elsewhere, and has entries . There is no assumption that , so ; when , there are no possible indices.
QAlgorithms.Papers.Independent.Intakethird100.P95.bobTransposeLinear
For arbitrary types and , this defines a complex-linear endomorphism of the complex matrices indexed by , given by
Thus it exchanges the two indices while retaining the indices. No finiteness assumptions are required.
QAlgorithms.Papers.Independent.Intakethird100.P95.coordinateMatrix
For any natural number , arbitrary real coefficients indexed by , and arbitrary real coefficients indexed by , this defines the complex matrix
where is the coordinate vector at , and has entries . Equivalently, the definition takes partial transpose on the second subsystem of the weighted symmetric and antisymmetric vector outer products. No sign conditions are imposed on the coefficients; for , both sums are empty and give the zero matrix.
QAlgorithms.Papers.Independent.Intakethird100.P95.exampleSpectrum
This defines the real vector indexed by whose entries, in index order, are
QAlgorithms.Papers.Independent.Intakethird100.P95.faceQuadratic
For arbitrary natural numbers and an arbitrary real vector , this function is
There are no sign assumptions on and no lower bound on . Division by zero gives zero, so when the entire first sum is zero. When , both sums are empty and the value is zero.
QAlgorithms.Papers.Independent.Intakethird100.P95.faceValue
For natural numbers , any real vector , and any subset , put , with this subtraction performed in . The defined real number is
No positivity, normalization, or stationarity condition is required. Real square roots of negative numbers are zero, and division by zero is defined to be zero; consequently the value is zero when or , and also when or is empty. The definition includes .
QAlgorithms.Papers.Independent.Intakethird100.P95.largestProbability
For any finite type and any function , this defines the real supremum of . When is nonempty, this is its maximum; when is empty, the real supremum convention gives . The entries are not assumed to be nonnegative or to sum to .
QAlgorithms.Papers.Independent.Intakethird100.P95.pptDistillableEntanglement
For any finite type , arbitrary , and Boolean , this defines the extended real supremum of all finite real numbers for which positive integer sequences satisfy and , and there exist complex linear, trace-preserving maps from matrices on to matrices on , with both and its conjugation by partial transpose on the second subsystem completely positive. Complete positivity requires preservation of positive semidefiniteness after adjoining an identity map of every finite dimension. Put , define to have entry when its indices are and zero elsewhere, and define to have entry when both its row pair and its column pair have equal coordinates and zero elsewhere. If is true, require for every ; otherwise require , where and the square root uses functional calculus. There is no nonnegativity or normalization assumption on ; negative real square-root arguments give zero, and may be empty. An empty set of rates has supremum , while an unbounded set has supremum .
QAlgorithms.Papers.Independent.Intakethird100.P95.pptEntanglementCost
For any finite type , arbitrary , and Boolean , this defines an extended real number as the infimum of all finite real rates with the following property: there are positive integer sequences , with and , and complex linear maps from matrices on to matrices on , such that each map preserves trace on every matrix and both the map and its conjugation by partial transpose on the second subsystem are completely positive. Here complete positivity requires positive semidefiniteness to be preserved after adjoining an identity map of every finite dimension. Let have entries when both its row pair and its column pair have equal coordinates, and zero otherwise; let , and let have entry at when and , and zero otherwise. For true, for every ; for false, the trace norm of tends to zero, where the trace norm of is , using functional calculus. No normalization or sign condition is imposed on , may be empty, and negative real square-root arguments give zero. The infimum is if there are no such rates.
QAlgorithms.Papers.Independent.Intakethird100.P95.preparationValue
For finite types , possibly empty, a natural number , and an arbitrary complex matrix indexed by , this defines a real number as the infimum of
where , and is real subtraction. There are no normalization assumptions on or , and and are included. This is the real-valued infimum operation, which returns zero when the displayed set is empty; the definition does not assert attainment.
QAlgorithms.Papers.Independent.Intakethird100.P95.preparationValues
For finite types , possibly empty, a natural number , and an arbitrary complex matrix indexed by , this defines the set of real numbers for which there exists a Hermitian positive semidefinite matrix on the same index set satisfying
Here , the inequalities mean positive semidefiniteness, and is real subtraction, including the value at . Neither nor is required to have trace one; the set can be empty.
QAlgorithms.Papers.Independent.Intakethird100.P95.primalValue
For arbitrary natural numbers and arbitrary real numbers , , this is the real infimum
where ranges over real arrays with , ranges over real arrays with , and
on the standard basis of . No restriction on or is imposed; negative square-root arguments and zero denominators are evaluated as zero. This definition does not assert that a minimizing array exists. The real infimum operation returns zero for an empty set or a set unbounded below; in particular, gives value zero.
QAlgorithms.Papers.Independent.Intakethird100.P95.primalValues
For arbitrary natural numbers and an arbitrary function , this is the set of real numbers for which there exist real arrays , , and , , satisfying
Here, with the standard basis of , the matrix appearing in the positivity condition is explicitly
No nonnegativity or normalization of , or lower bound on , is assumed. Real square roots of negative arguments are zero, and division by zero is zero, including the denominator when . When , all constraints are vacuous or concern the empty matrix, and the set is .
QAlgorithms.Papers.Independent.Intakethird100.P95.rankOneValue
For arbitrary natural numbers and an arbitrary real array , , this is the real supremum of
over real arrays , , and , , satisfying , , and requiring
to be positive semidefinite of complex rank exactly one, where is the standard basis of . No restrictions on or are imposed, and no maximizing arrays are asserted to exist. The total conventions give zero for zero-denominator divisions and negative square-root arguments, with ; the real supremum is zero for an empty set or a set unbounded above. In particular, the value is zero when , because its feasible set is empty.
QAlgorithms.Papers.Independent.Intakethird100.P95.rankOneValues
For arbitrary natural numbers and arbitrary real numbers , , this is the set of all real values
for which there exist real arrays , , and , , such that , , and the matrix
is positive semidefinite and has rank exactly one over . The vectors form the standard basis of . There are no restrictions on or ; real powers and logarithms use their total definitions, zero denominators give zero, and negative square-root arguments give zero. For , no matrix has the required rank, so this set is empty.
QAlgorithms.Papers.Independent.Intakethird100.P95.renyiInfinity
For any finite type and any function , let , with when is empty. This defines , using the total real logarithm. Thus its value is if , and otherwise. No probability-distribution assumptions are imposed on .
QAlgorithms.Papers.Independent.Intakethird100.P95.rhoPower
For any type , any function , and any natural number , this defines a complex matrix on pairs of -tuples. Writing , its entries are
The real square root is embedded into and equals on negative arguments. No nonnegativity or normalization of is assumed. For , this is the one-by-one matrix with entry , even if is empty.
QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtDensity
For any type and any function , this defines the complex matrix indexed by whose entries are
with the real square root embedded into . The real square root is on negative arguments, and the diagonal entry at is . No finiteness, nonnegativity, or normalization assumption is made.
QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtDualObjective
For arbitrary natural numbers , arbitrary real numbers , , and arbitrary real arrays , , and , , define
where is the total real-power operation. The objective is
No feasibility, sign, or normalization conditions are assumed. All operations are total: division by zero and square roots of negative real numbers give zero, and . Thus , , and are included; when , the objective is , with the same zero-division convention.
QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtDualValues
For arbitrary natural numbers and an arbitrary real array , , this is the set of all values
obtained from real arrays , , and , , with , , and
Here is the standard basis of , and is total real exponentiation. No assumptions on or lower bound on are imposed; zero denominators and negative square-root arguments give zero, and . For , the set is the singleton containing .
QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtPrimalObjective
For every natural number , including , and arbitrary real arrays indexed by and indexed by , the objective is
There are no sign or feasibility assumptions on either array. Empty sums are zero, so the objective is zero when .
QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtRenyi
For any finite type , any real number , and any function , define , using the total real logarithm: , and when . The defined value is
where the powers in the last branch use the total real-power function. Neither nonnegativity nor normalization of is assumed, and is unrestricted. If is empty, the value is in every branch.
QAlgorithms.Papers.Independent.Intakethird100.P95.symmetricSchmidtProjector
For any natural number and indices , let be the coordinate vector equal to one at and zero elsewhere. This defines the matrix
For a vector , the matrix has entry at . The indices are unrestricted apart from lying between zero and ; when , there are no such indices.
QAlgorithms.Papers.Independent.Intakethird100.P95.tensorSpectrum
For any type , any function , and any natural number , this defines a real-valued function on -tuples by
No assumptions are made on the entries or on the size of . For , the domain contains the unique empty tuple and its value is the empty product .
Confirmed by the mission captain (proposal self-audit).
Confirmed by the moderator at approval.