Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Pure-state transformations and catalysis under operations that completely preserve positivity of partial transpose

Definition
QAlgorithms_Papers_Independent_Intakethird100_P95

by Gohan · Nov 1, 2026 · Mathlib 0df444a (Lean v4.33.1)

prove2quantum

The 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.
Definition code
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
Source
Pure-state transformations and catalysis under operations that completely preserve positivity of partial transpose, arXiv:0801.4322v2, https://arxiv.org/abs/0801.4322v2
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 ddd, this abbreviates the type of ordered pairs (i,j)(i,j)(i,j) of integers with 0≤i,j<d0\le i,j<d0≤i,j<d and j<ij<ij<i. Diagonal pairs are excluded, and the type is empty when d=0d=0d=0 or d=1d=1d=1.

QAlgorithms.Papers.Independent.Intakethird100.P95.CatalyticPreparation

For any finite type III, any natural number KKK, and any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, this predicate means that there exist a natural number C≥2C\ge2C≥2 and a complex-linear map L\mathcal LL from matrices on {0,…,C−1}×{0,…,K−1}\{0,\ldots,C-1\}\times\{0,\ldots,K-1\}{0,…,C−1}×{0,…,K−1} to matrices on {0,…,C−1}×I\{0,\ldots,C-1\}\times I{0,…,C−1}×I such that L\mathcal LL is completely positive, preserves the trace of every matrix, its conjugation by partial transpose on the second subsystem is completely positive, and

L(ΩC⊗ΩK)=ΩC⊗Dλ.\mathcal L(\Omega_C\otimes\Omega_K)=\Omega_C\otimes D_\lambda.L(ΩC​⊗ΩK​)=ΩC​⊗Dλ​.

The products here group the two first subsystems together and the two second subsystems together. Explicitly, (Ωm)(i,j),(k,l)=1/m(\Omega_m)_{(i,j),(k,l)}=1/m(Ωm​)(i,j),(k,l)​=1/m if i=ji=ji=j and k=lk=lk=l, and is zero otherwise; (Dλ)(i,j),(k,l)=λiλk(D_\lambda)_{(i,j),(k,l)}=\sqrt{\lambda_i\lambda_k}(Dλ​)(i,j),(k,l)​=λi​λk​​ if i=ji=ji=j and k=lk=lk=l, and is zero otherwise. Partial transpose sends X(a,b),(c,d)X_{(a,b),(c,d)}X(a,b),(c,d)​ to X(a,d),(c,b)X_{(a,d),(c,b)}X(a,d),(c,b)​; complete positivity requires every finite ancillary block amplification, including size zero, to preserve positive semidefiniteness. No nonnegativity or normalization of λ\lambdaλ is assumed. Empty III and K=0K=0K=0 are included, and negative real square-root arguments are evaluated as zero.

QAlgorithms.Papers.Independent.Intakethird100.P95.CostRates

For a finite type III, an arbitrary function λ:I→R\lambda:I\to\mathbb Rλ:I→R, and a Boolean flag zeroError\mathrm{zeroError}zeroError, this defines a set of extended real numbers containing precisely those rrr for which there exist a finite real number RRR, sequences n,K:N→Nn,K:\mathbb N\to\mathbb Nn,K:N→N, and maps Φj\Phi_jΦj​ from matrices on Fin⁡(Kj)×Fin⁡(Kj)\operatorname{Fin}(K_j)\times\operatorname{Fin}(K_j)Fin(Kj​)×Fin(Kj​) to matrices on Inj×InjI^{n_j}\times I^{n_j}Inj​×Inj​, such that r=Rr=Rr=R, every njn_jnj​ and KjK_jKj​ is positive, nj→+∞n_j\to+\inftynj​→+∞, and log⁡(Kj)/(njlog⁡2)→R\log(K_j)/(n_j\log2)\to Rlog(Kj​)/(nj​log2)→R. Each Φj\Phi_jΦj​ must be complex-linear and trace preserving on every matrix, and both Φj\Phi_jΦj​ 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 000. Define

ΩK((a,b),(c,d))={1/K,a=b and c=d,0,otherwise,wn(x)=∏i=0n−1λxi,\Omega_K\bigl((a,b),(c,d)\bigr)= \begin{cases} 1/K,&a=b\ \text{and}\ c=d,\\ 0,&\text{otherwise}, \end{cases} \qquad w_n(x)=\prod_{i=0}^{n-1}\lambda_{x_i},ΩK​((a,b),(c,d))={1/K,0,​a=b and c=d,otherwise,​wn​(x)=i=0∏n−1​λxi​​,

and

ρλ,n((x,y),(u,v))={wn(x)wn(u),x=y and u=v,0,otherwise.\rho_{\lambda,n}\bigl((x,y),(u,v)\bigr)= \begin{cases} \sqrt{w_n(x)w_n(u)},&x=y\ \text{and}\ u=v,\\ 0,&\text{otherwise}. \end{cases}ρλ,n​((x,y),(u,v))={wn​(x)wn​(u)​,0,​x=y and u=v,otherwise.​

If the flag is true, the requirement is Φj(ΩKj)=ρλ,nj\Phi_j(\Omega_{K_j})=\rho_{\lambda,n_j}Φj​(ΩKj​​)=ρλ,nj​​ for every jjj; if it is false, the requirement is ∥Φj(ΩKj)−ρλ,nj∥1→0\|\Phi_j(\Omega_{K_j})-\rho_{\lambda,n_j}\|_1\to0∥Φj​(ΩKj​​)−ρλ,nj​​∥1​→0, where ∥M∥1=Re⁡tr⁡M†M\|M\|_1=\operatorname{Re}\operatorname{tr}\sqrt{M^\dagger M}∥M∥1​=RetrM†M​. No nonnegativity or normalization assumption is imposed on λ\lambdaλ, no monotonicity of the sequences is required, and neither infinite extended-real values nor an explicit restriction R≥0R\ge0R≥0 appears among the witnesses. If III is empty, the set is empty because positive njn_jnj​ 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 III, any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, and any Boolean zzz, this defines a set of extended real numbers consisting precisely of the finite numbers R∈RR\in\mathbb RR∈R for which there are sequences of positive integers nj,Kjn_j,K_jnj​,Kj​, with nj→∞n_j\to\inftynj​→∞ and log⁡2(Kj)/nj→R\log_2(K_j)/n_j\to Rlog2​(Kj​)/nj​→R, and maps Lj\mathcal L_jLj​ from matrices on (Inj)×(Inj)(I^{n_j})\times(I^{n_j})(Inj​)×(Inj​) to matrices on {0,…,Kj−1}2\{0,\ldots,K_j-1\}^2{0,…,Kj​−1}2, satisfying the following requirements. Each Lj\mathcal L_jLj​ 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 λn(u)=∏i=0n−1λ(ui)\lambda_n(u)=\prod_{i=0}^{n-1}\lambda(u_i)λn​(u)=∏i=0n−1​λ(ui​), and define Dλ,nD_{\lambda,n}Dλ,n​ by

(Dλ,n)(u,v),(u′,v′)={λn(u)λn(u′),u=v and u′=v′,0,otherwise.(D_{\lambda,n})_{(u,v),(u',v')}= \begin{cases} \sqrt{\lambda_n(u)\lambda_n(u')},&u=v\text{ and }u'=v',\\ 0,&\text{otherwise}. \end{cases}(Dλ,n​)(u,v),(u′,v′)​={λn​(u)λn​(u′)​,0,​u=v and u′=v′,otherwise.​

Define BKB_KBK​ to have entry 1/K1/K1/K at ((a,b),(c,d))((a,b),(c,d))((a,b),(c,d)) when a=ba=ba=b and c=dc=dc=d, and zero otherwise. If zzz is true, Lj(Dλ,nj)=BKj\mathcal L_j(D_{\lambda,n_j})=B_{K_j}Lj​(Dλ,nj​​)=BKj​​ for every jjj; if zzz is false, tr⁡Hj†Hj→0\operatorname{tr}\sqrt{H_j^\dagger H_j}\to0trHj†​Hj​​→0, where Hj=BKj−Lj(Dλ,nj)H_j=B_{K_j}-\mathcal L_j(D_{\lambda,n_j})Hj​=BKj​​−Lj​(Dλ,nj​​), with the square root interpreted through functional calculus. No nonnegativity or normalization is assumed for λ\lambdaλ, and real square roots of negative numbers are zero. The type III may be empty, although all njn_jnj​ and KjK_jKj​ must be positive; neither +∞+\infty+∞ nor −∞-\infty−∞ belongs to this set.

QAlgorithms.Papers.Independent.Intakethird100.P95.FaceStationary

For natural numbers d,Kd,Kd,K, any real vector λ=(λi)i=0d−1\lambda=(\lambda_i)_{i=0}^{d-1}λ=(λi​)i=0d−1​, and any subset C⊆{0,…,d−1}C\subseteq\{0,\ldots,d-1\}C⊆{0,…,d−1}, this predicate means that there exists a real vector xxx such that ∣xi∣≤λi|x_i|\le\sqrt{\lambda_i}∣xi​∣≤λi​​ for every iii, xi=λix_i=\sqrt{\lambda_i}xi​=λi​​ for every i∈Ci\in Ci∈C, and

∑i=0d−1xi−Kxj=0for every j∉C.\sum_{i=0}^{d-1}x_i-Kx_j=0 \qquad\text{for every }j\notin C.i=0∑d−1​xi​−Kxj​=0for every j∈/C.

There are no positivity or normalization assumptions on λ\lambdaλ and no restriction on KKK. Real square roots of negative inputs are defined to be zero. If CCC is the entire index set, the final equations impose no conditions; if d=0d=0d=0, the predicate holds using the empty vector.

QAlgorithms.Papers.Independent.Intakethird100.P95.IsUniformRank

For any finite type III, equipped with a finite enumeration, any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, and any natural number KKK, this predicate means

#{i∈I:λi≠0}=Kand∀i∈I, λi=0 or λi=K−1.\#\{i\in I:\lambda_i\ne0\}=K \quad\text{and}\quad \forall i\in I,\ \lambda_i=0\ \text{or}\ \lambda_i=K^{-1}.#{i∈I:λi​=0}=Kand∀i∈I, λi​=0 or λi​=K−1.

For K=0K=0K=0, the inverse is zero and the predicate holds exactly when λ\lambdaλ is identically zero. For K>0K>0K>0, exactly KKK entries equal 1/K1/K1/K and every other entry is zero. Empty III is allowed and then the predicate holds exactly for K=0K=0K=0.

QAlgorithms.Papers.Independent.Intakethird100.P95.OnFace

For every natural number ddd, arbitrary real vectors λ,x:{0,…,d−1}→R\lambda,x:\{0,\ldots,d-1\}\to\mathbb Rλ,x:{0,…,d−1}→R, and any subset C⊆{0,…,d−1}C\subseteq\{0,\ldots,d-1\}C⊆{0,…,d−1}, this predicate means

∣xi∣≤λifor every i,xi=λifor every i∈C.|x_i|\le\sqrt{\lambda_i}\quad\text{for every }i, \qquad x_i=\sqrt{\lambda_i}\quad\text{for every }i\in C.∣xi​∣≤λi​​for every i,xi​=λi​​for every i∈C.

No nonnegativity or normalization of λ\lambdaλ is assumed. The real square root of a negative number is zero, so λi≤0\lambda_i\le0λi​≤0 forces xi=0x_i=0xi​=0. If CCC is empty, the equality requirements are vacuous; if d=0d=0d=0, the entire predicate holds vacuously.

QAlgorithms.Papers.Independent.Intakethird100.P95.PPTChannel

For finite types A,B,C,DA,B,C,DA,B,C,D, this structure consists of a complex-linear map Φ\PhiΦ from matrices on A×BA\times BA×B to matrices on C×DC\times DC×D, together with proofs that Φ\PhiΦ is completely positive, preserves the complex trace of every matrix, and that ΓD∘Φ∘ΓB\Gamma_D\circ\Phi\circ\Gamma_BΓD​∘Φ∘ΓB​ is completely positive. Here (ΓBX)(a,b),(a′,b′)=X(a,b′),(a′,b)(\Gamma_BX)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}(ΓB​X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, and ΓD\Gamma_DΓD​ is defined analogously on the output. Complete positivity means that, for every natural number hhh, including 000, the amplification acting on each matrix block by the given map sends every positive semidefinite matrix on Fin⁡(h)×(A×B)\operatorname{Fin}(h)\times(A\times B)Fin(h)×(A×B) to a positive semidefinite matrix on Fin⁡(h)×(C×D)\operatorname{Fin}(h)\times(C\times D)Fin(h)×(C×D); the same condition applies to the transposed composite. The four finite types may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P95.PreparationFeasible

For finite types A,BA,BA,B, possibly empty, a natural number KKK, and arbitrary complex matrices X,PX,PX,P indexed by A×BA\times BA×B, this predicate requires PPP to be Hermitian and positive semidefinite and requires

XΓ+(K−1)PΓ⪰0,(K+1)PΓ−XΓ⪰0,X^\Gamma+(K-1)P^\Gamma\succeq0, \qquad (K+1)P^\Gamma-X^\Gamma\succeq0,XΓ+(K−1)PΓ⪰0,(K+1)PΓ−XΓ⪰0,

where M(a,b),(c,d)Γ=M(a,d),(c,b)M^\Gamma_{(a,b),(c,d)}=M_{(a,d),(c,b)}M(a,b),(c,d)Γ​=M(a,d),(c,b)​, and ⪰0\succeq0⪰0 means positive semidefinite. The subtraction K−1K-1K−1 is performed in the real numbers: K=0K=0K=0 gives the two conditions XΓ−PΓ⪰0X^\Gamma-P^\Gamma\succeq0XΓ−PΓ⪰0 and PΓ−XΓ⪰0P^\Gamma-X^\Gamma\succeq0PΓ−XΓ⪰0, while K=1K=1K=1 gives XΓ⪰0X^\Gamma\succeq0XΓ⪰0 and 2PΓ−XΓ⪰02P^\Gamma-X^\Gamma\succeq02PΓ−XΓ⪰0. No trace normalization or separate hypothesis on XXX is included.

QAlgorithms.Papers.Independent.Intakethird100.P95.RateLimit

For sequences n,K:N→Nn,K:\mathbb N\to\mathbb Nn,K:N→N and a real number RRR, this predicate means

lim⁡j→∞log⁡(Kj)/log⁡2nj=R.\lim_{j\to\infty}\frac{\log(K_j)/\log 2}{n_j}=R.j→∞lim​nj​log(Kj​)/log2​=R.

This predicate itself imposes no positivity, divergence, or monotonicity conditions on either sequence. The total-function conventions assign log⁡0=0\log 0=0log0=0 and division by zero the value 000, so terms with Kj=0K_j=0Kj​=0 or nj=0n_j=0nj​=0 are included and have value 000.

QAlgorithms.Papers.Independent.Intakethird100.P95.Reach

For finite types A,B,C,DA,B,C,DA,B,C,D and arbitrary complex matrices XXX on A×BA\times BA×B and YYY on C×DC\times DC×D, this predicate means that there exists a complex-linear map Φ\PhiΦ with Φ(X)=Y\Phi(X)=YΦ(X)=Y which preserves the trace of every matrix and for which both Φ\PhiΦ and ΓD∘Φ∘ΓB\Gamma_D\circ\Phi\circ\Gamma_BΓD​∘Φ∘ΓB​ are completely positive. Here ΓB\Gamma_BΓB​ exchanges the two BBB indices of a matrix entry, and ΓD\Gamma_DΓD​ exchanges the two DDD indices. Complete positivity requires positivity preservation after amplification by the identity on Fin⁡(h)\operatorname{Fin}(h)Fin(h) for every h≥0h\ge0h≥0. Neither XXX nor YYY is assumed to be positive semidefinite or normalized, and empty indexing types are allowed.

QAlgorithms.Papers.Independent.Intakethird100.P95.SchmidtDualFeasible

For every natural number ddd, and real arrays μij\mu_{ij}μij​, 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and tijt_{ij}tij​, 0≤j<i<d0\le j<i<d0≤j<i<d, this predicate means that every μij≤1\mu_{ij}\le1μij​≤1, every tij≤1t_{ij}\le1tij​≤1, and the following matrix on the standard basis ∣ij⟩|ij\rangle∣ij⟩ of Cd×d\mathbb C^{d\times d}Cd×d is positive semidefinite:

C(μ,t)=∑iμii∣ii⟩⟨ii∣+12∑j<i[(μij+tij)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij−tij)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)].C(\mu,t)= \sum_i\mu_{ii}|ii\rangle\langle ii| +\frac12\sum_{j<i}\left[ (\mu_{ij}+t_{ij})(|ij\rangle\langle ij|+|ji\rangle\langle ji|) +(\mu_{ij}-t_{ij})(|ii\rangle\langle jj|+|jj\rangle\langle ii|) \right].C(μ,t)=i∑​μii​∣ii⟩⟨ii∣+21​j<i∑​[(μij​+tij​)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij​−tij​)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)].

The coordinate bounds contain no explicit lower bounds. For d=0d=0d=0, the arrays and matrix are empty and the predicate holds.

QAlgorithms.Papers.Independent.Intakethird100.P95.SchmidtPrimalFeasible

For any natural numbers d,Kd,Kd,K, arbitrary real function λ:{0,…,d−1}→R\lambda:\{0,\ldots,d-1\}\to\mathbb Rλ:{0,…,d−1}→R, real coefficients sijs_{ij}sij​ for 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and real coefficients aija_{ij}aij​ for 0≤j<i<d0\le j<i<d0≤j<i<d, this predicate requires

sij≥λiλjK+1(j≤i),aij≥λiλjK−1(j<i),s_{ij}\ge\frac{\sqrt{\lambda_i\lambda_j}}{K+1}\quad(j\le i), \qquad a_{ij}\ge\frac{\sqrt{\lambda_i\lambda_j}}{K-1}\quad(j<i),sij​≥K+1λi​λj​​​(j≤i),aij​≥K−1λi​λj​​​(j<i),

and requires the following matrix to be positive semidefinite:

∑isiieiieii†+12∑j<i[(sij+aij)(eijeij†+ejieji†)+(sij−aij)(eiiejj†+ejjeii†)].\sum_i s_{ii}e_{ii}e_{ii}^\dagger+ \frac12\sum_{j<i} \left[ (s_{ij}+a_{ij})(e_{ij}e_{ij}^\dagger+e_{ji}e_{ji}^\dagger) +(s_{ij}-a_{ij})(e_{ii}e_{jj}^\dagger+e_{jj}e_{ii}^\dagger) \right].i∑​sii​eii​eii†​+21​j<i∑​[(sij​+aij​)(eij​eij†​+eji​eji†​)+(sij​−aij​)(eii​ejj†​+ejj​eii†​)].

Here eije_{ij}eij​ is the coordinate vector at (i,j)(i,j)(i,j), and uv†uv^\daggeruv† has entries upvq‾u_p\overline{v_q}up​vq​​. No nonnegativity or normalization assumption is imposed on λ\lambdaλ, and a negative real square-root argument gives zero. The denominators use real subtraction and total division: at K=1K=1K=1, every antisymmetric lower bound is zero; at K=0K=0K=0, it is −λiλj-\sqrt{\lambda_i\lambda_j}−λi​λj​​. When d=0d=0d=0, the coefficient conditions are vacuous and the empty matrix is positive semidefinite.

QAlgorithms.Papers.Independent.Intakethird100.P95.SymIndex

For each natural number ddd, this abbreviates the type of ordered pairs (i,j)(i,j)(i,j) of integers with 0≤i,j<d0\le i,j<d0≤i,j<d and j≤ij\le ij≤i. Diagonal pairs are included, and the type is empty when d=0d=0d=0.

QAlgorithms.Papers.Independent.Intakethird100.P95.ValidLengths

For two sequences n,K:N→Nn,K:\mathbb N\to\mathbb Nn,K:N→N, this predicate means that nj>0n_j>0nj​>0 and Kj>0K_j>0Kj​>0 for every jjj, and that nj→+∞n_j\to+\inftynj​→+∞ as j→+∞j\to+\inftyj→+∞: for every natural bound, all sufficiently late njn_jnj​ meet or exceed that bound. Neither sequence is required to be monotone, and no growth condition is imposed on KKK.

QAlgorithms.Papers.Independent.Intakethird100.P95.antisymmetricSchmidtProjector

For any natural number ddd and indices 0≤i,j<d0\le i,j<d0≤i,j<d, this defines

Aij=12(eij−eji)(eij−eji)†,A_{ij}=\frac12(e_{ij}-e_{ji})(e_{ij}-e_{ji})^\dagger,Aij​=21​(eij​−eji​)(eij​−eji​)†,

where eije_{ij}eij​ is the coordinate vector equal to one at (i,j)(i,j)(i,j) and zero elsewhere, and vv†vv^\daggervv† has entries vpvq‾v_p\overline{v_q}vp​vq​​. There is no assumption that i≠ji\ne ji=j, so Aii=0A_{ii}=0Aii​=0; when d=0d=0d=0, there are no possible indices.

QAlgorithms.Papers.Independent.Intakethird100.P95.bobTransposeLinear

For arbitrary types AAA and BBB, this defines a complex-linear endomorphism ΓB\Gamma_BΓB​ of the complex matrices indexed by A×BA\times BA×B, given by

(ΓBX)((a,b),(a′,b′))=X((a,b′),(a′,b)).(\Gamma_BX)\bigl((a,b),(a',b')\bigr) =X\bigl((a,b'),(a',b)\bigr).(ΓB​X)((a,b),(a′,b′))=X((a,b′),(a′,b)).

Thus it exchanges the two BBB indices while retaining the AAA indices. No finiteness assumptions are required.

QAlgorithms.Papers.Independent.Intakethird100.P95.coordinateMatrix

For any natural number ddd, arbitrary real coefficients sijs_{ij}sij​ indexed by 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and arbitrary real coefficients aija_{ij}aij​ indexed by 0≤j<i<d0\le j<i<d0≤j<i<d, this defines the complex matrix

C(s,a)=∑i=0d−1siieiieii†+12∑0≤j<i<d[(sij+aij)(eijeij†+ejieji†)+(sij−aij)(eiiejj†+ejjeii†)],C(s,a)= \sum_{i=0}^{d-1}s_{ii}e_{ii}e_{ii}^\dagger+ \frac12\sum_{0\le j<i<d} \left[ (s_{ij}+a_{ij})(e_{ij}e_{ij}^\dagger+e_{ji}e_{ji}^\dagger) +(s_{ij}-a_{ij})(e_{ii}e_{jj}^\dagger+e_{jj}e_{ii}^\dagger) \right],C(s,a)=i=0∑d−1​sii​eii​eii†​+21​0≤j<i<d∑​[(sij​+aij​)(eij​eij†​+eji​eji†​)+(sij​−aij​)(eii​ejj†​+ejj​eii†​)],

where eije_{ij}eij​ is the coordinate vector at (i,j)(i,j)(i,j), and uv†uv^\daggeruv† has entries upvq‾u_p\overline{v_q}up​vq​​. 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 d=0d=0d=0, both sums are empty and give the zero matrix.

QAlgorithms.Papers.Independent.Intakethird100.P95.exampleSpectrum

This defines the real vector indexed by {0,…,5}\{0,\ldots,5\}{0,…,5} whose entries, in index order, are

(120,120,120,420,420,920).\left(\frac1{20},\frac1{20},\frac1{20}, \frac4{20},\frac4{20},\frac9{20}\right).(201​,201​,201​,204​,204​,209​).

QAlgorithms.Papers.Independent.Intakethird100.P95.faceQuadratic

For arbitrary natural numbers d,Kd,Kd,K and an arbitrary real vector x=(x0,…,xd−1)x=(x_0,\ldots,x_{d-1})x=(x0​,…,xd−1​), this function is

∑0≤j<i<dxixjK−1−∑0≤j≤i<dxixjK+1.\sum_{0\le j<i<d}\frac{x_i x_j}{K-1} -\sum_{0\le j\le i<d}\frac{x_i x_j}{K+1}.0≤j<i<d∑​K−1xi​xj​​−0≤j≤i<d∑​K+1xi​xj​​.

There are no sign assumptions on xxx and no lower bound on KKK. Division by zero gives zero, so when K=1K=1K=1 the entire first sum is zero. When d=0d=0d=0, both sums are empty and the value is zero.

QAlgorithms.Papers.Independent.Intakethird100.P95.faceValue

For natural numbers d,Kd,Kd,K, any real vector λ=(λi)i=0d−1\lambda=(\lambda_i)_{i=0}^{d-1}λ=(λi​)i=0d−1​, and any subset C⊆{0,…,d−1}C\subseteq\{0,\ldots,d-1\}C⊆{0,…,d−1}, put D=K+∣C∣−dD=K+|C|-dD=K+∣C∣−d, with this subtraction performed in R\mathbb RR. The defined real number is

K(K2−1)D[(∑i∈Cλi)2−D∑i∈Cλi].\frac{K}{(K^2-1)D} \left[ \left(\sum_{i\in C}\sqrt{\lambda_i}\right)^2 -D\sum_{i\in C}\lambda_i \right].(K2−1)DK​​(i∈C∑​λi​​)2−Di∈C∑​λi​​.

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 K=1K=1K=1 or D=0D=0D=0, and also when K=0K=0K=0 or CCC is empty. The definition includes d=0d=0d=0.

QAlgorithms.Papers.Independent.Intakethird100.P95.largestProbability

For any finite type III and any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, this defines the real supremum of {λi:i∈I}\{\lambda_i:i\in I\}{λi​:i∈I}. When III is nonempty, this is its maximum; when III is empty, the real supremum convention gives 000. The entries are not assumed to be nonnegative or to sum to 111.

QAlgorithms.Papers.Independent.Intakethird100.P95.pptDistillableEntanglement

For any finite type III, arbitrary λ:I→R\lambda:I\to\mathbb Rλ:I→R, and Boolean zzz, this defines the extended real supremum of all finite real numbers RRR for which positive integer sequences nj,Kjn_j,K_jnj​,Kj​ satisfy nj→∞n_j\to\inftynj​→∞ and log⁡2(Kj)/nj→R\log_2(K_j)/n_j\to Rlog2​(Kj​)/nj​→R, and there exist complex linear, trace-preserving maps Lj\mathcal L_jLj​ from matrices on Inj×InjI^{n_j}\times I^{n_j}Inj​×Inj​ to matrices on {0,…,Kj−1}2\{0,\ldots,K_j-1\}^2{0,…,Kj​−1}2, with both Lj\mathcal L_jLj​ 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 λn(u)=∏i=0n−1λ(ui)\lambda_n(u)=\prod_{i=0}^{n-1}\lambda(u_i)λn​(u)=∏i=0n−1​λ(ui​), define Dλ,nD_{\lambda,n}Dλ,n​ to have entry λn(u)λn(u′)\sqrt{\lambda_n(u)\lambda_n(u')}λn​(u)λn​(u′)​ when its indices are ((u,u),(u′,u′))((u,u),(u',u'))((u,u),(u′,u′)) and zero elsewhere, and define BKB_KBK​ to have entry 1/K1/K1/K when both its row pair and its column pair have equal coordinates and zero elsewhere. If zzz is true, require Lj(Dλ,nj)=BKj\mathcal L_j(D_{\lambda,n_j})=B_{K_j}Lj​(Dλ,nj​​)=BKj​​ for every jjj; otherwise require Re⁡tr⁡Hj†Hj→0\operatorname{Re}\operatorname{tr}\sqrt{H_j^\dagger H_j}\to0RetrHj†​Hj​​→0, where Hj=BKj−Lj(Dλ,nj)H_j=B_{K_j}-\mathcal L_j(D_{\lambda,n_j})Hj​=BKj​​−Lj​(Dλ,nj​​) and the square root uses functional calculus. There is no nonnegativity or normalization assumption on λ\lambdaλ; negative real square-root arguments give zero, and III may be empty. An empty set of rates has supremum −∞-\infty−∞, while an unbounded set has supremum +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P95.pptEntanglementCost

For any finite type III, arbitrary λ:I→R\lambda:I\to\mathbb Rλ:I→R, and Boolean zzz, this defines an extended real number as the infimum of all finite real rates RRR with the following property: there are positive integer sequences nj,Kjn_j,K_jnj​,Kj​, with nj→∞n_j\to\inftynj​→∞ and log⁡2(Kj)/nj→R\log_2(K_j)/n_j\to Rlog2​(Kj​)/nj​→R, and complex linear maps Lj\mathcal L_jLj​ from matrices on {0,…,Kj−1}2\{0,\ldots,K_j-1\}^2{0,…,Kj​−1}2 to matrices on Inj×InjI^{n_j}\times I^{n_j}Inj​×Inj​, 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 BKB_KBK​ have entries 1/K1/K1/K when both its row pair and its column pair have equal coordinates, and zero otherwise; let λn(u)=∏i=0n−1λ(ui)\lambda_n(u)=\prod_{i=0}^{n-1}\lambda(u_i)λn​(u)=∏i=0n−1​λ(ui​), and let Dλ,nD_{\lambda,n}Dλ,n​ have entry λn(u)λn(u′)\sqrt{\lambda_n(u)\lambda_n(u')}λn​(u)λn​(u′)​ at ((u,v),(u′,v′))((u,v),(u',v'))((u,v),(u′,v′)) when u=vu=vu=v and u′=v′u'=v'u′=v′, and zero otherwise. For zzz true, Lj(BKj)=Dλ,nj\mathcal L_j(B_{K_j})=D_{\lambda,n_j}Lj​(BKj​​)=Dλ,nj​​ for every jjj; for zzz false, the trace norm of Lj(BKj)−Dλ,nj\mathcal L_j(B_{K_j})-D_{\lambda,n_j}Lj​(BKj​​)−Dλ,nj​​ tends to zero, where the trace norm of HHH is Re⁡tr⁡H†H\operatorname{Re}\operatorname{tr}\sqrt{H^\dagger H}RetrH†H​, using functional calculus. No normalization or sign condition is imposed on λ\lambdaλ, III may be empty, and negative real square-root arguments give zero. The infimum is +∞+\infty+∞ if there are no such rates.

QAlgorithms.Papers.Independent.Intakethird100.P95.preparationValue

For finite types A,BA,BA,B, possibly empty, a natural number KKK, and an arbitrary complex matrix XXX indexed by A×BA\times BA×B, this defines a real number as the infimum of

{Re⁡tr⁡P:P=P†, P⪰0, XΓ+(K−1)PΓ⪰0, (K+1)PΓ−XΓ⪰0},\left\{\operatorname{Re}\operatorname{tr}P: P=P^\dagger,\ P\succeq0,\ X^\Gamma+(K-1)P^\Gamma\succeq0,\ (K+1)P^\Gamma-X^\Gamma\succeq0\right\},{RetrP:P=P†, P⪰0, XΓ+(K−1)PΓ⪰0, (K+1)PΓ−XΓ⪰0},

where M(a,b),(c,d)Γ=M(a,d),(c,b)M^\Gamma_{(a,b),(c,d)}=M_{(a,d),(c,b)}M(a,b),(c,d)Γ​=M(a,d),(c,b)​, and K−1K-1K−1 is real subtraction. There are no normalization assumptions on XXX or PPP, and K=0K=0K=0 and K=1K=1K=1 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 A,BA,BA,B, possibly empty, a natural number KKK, and an arbitrary complex matrix XXX indexed by A×BA\times BA×B, this defines the set of real numbers vvv for which there exists a Hermitian positive semidefinite matrix PPP on the same index set satisfying

XΓ+(K−1)PΓ⪰0,(K+1)PΓ−XΓ⪰0,Re⁡tr⁡P=v.X^\Gamma+(K-1)P^\Gamma\succeq0, \qquad (K+1)P^\Gamma-X^\Gamma\succeq0, \qquad \operatorname{Re}\operatorname{tr}P=v.XΓ+(K−1)PΓ⪰0,(K+1)PΓ−XΓ⪰0,RetrP=v.

Here M(a,b),(c,d)Γ=M(a,d),(c,b)M^\Gamma_{(a,b),(c,d)}=M_{(a,d),(c,b)}M(a,b),(c,d)Γ​=M(a,d),(c,b)​, the inequalities mean positive semidefiniteness, and K−1K-1K−1 is real subtraction, including the value −1-1−1 at K=0K=0K=0. Neither XXX nor PPP is required to have trace one; the set can be empty.

QAlgorithms.Papers.Independent.Intakethird100.P95.primalValue

For arbitrary natural numbers d,Kd,Kd,K and arbitrary real numbers λi\lambda_iλi​, 0≤i<d0\le i<d0≤i<d, this is the real infimum

inf⁡{∑j≤isij+∑j<iaij:λiλjK+1≤sij, λiλjK−1≤aij, C(s,a)⪰0},\inf\left\{ \sum_{j\le i}s_{ij}+\sum_{j<i}a_{ij}: \frac{\sqrt{\lambda_i\lambda_j}}{K+1}\le s_{ij},\ \frac{\sqrt{\lambda_i\lambda_j}}{K-1}\le a_{ij},\ C(s,a)\succeq0 \right\},inf{j≤i∑​sij​+j<i∑​aij​:K+1λi​λj​​​≤sij​, K−1λi​λj​​​≤aij​, C(s,a)⪰0},

where sijs_{ij}sij​ ranges over real arrays with 0≤j≤i<d0\le j\le i<d0≤j≤i<d, aija_{ij}aij​ ranges over real arrays with 0≤j<i<d0\le j<i<d0≤j<i<d, and

C(s,a)=∑isii∣ii⟩⟨ii∣+12∑j<i[(sij+aij)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(sij−aij)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]C(s,a)= \sum_i s_{ii}|ii\rangle\langle ii| +\frac12\sum_{j<i}\left[ (s_{ij}+a_{ij})(|ij\rangle\langle ij|+|ji\rangle\langle ji|) +(s_{ij}-a_{ij})(|ii\rangle\langle jj|+|jj\rangle\langle ii|) \right]C(s,a)=i∑​sii​∣ii⟩⟨ii∣+21​j<i∑​[(sij​+aij​)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(sij​−aij​)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]

on the standard basis ∣ij⟩|ij\rangle∣ij⟩ of Cd×d\mathbb C^{d\times d}Cd×d. No restriction on λ\lambdaλ or KKK 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, d=0d=0d=0 gives value zero.

QAlgorithms.Papers.Independent.Intakethird100.P95.primalValues

For arbitrary natural numbers d,Kd,Kd,K and an arbitrary function λ:{0,…,d−1}→R\lambda:\{0,\ldots,d-1\}\to\mathbb Rλ:{0,…,d−1}→R, this is the set of real numbers vvv for which there exist real arrays sijs_{ij}sij​, 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and aija_{ij}aij​, 0≤j<i<d0\le j<i<d0≤j<i<d, satisfying

λiλjK+1≤sij,λiλjK−1≤aij,C(s,a)⪰0,v=∑j≤isij+∑j<iaij.\frac{\sqrt{\lambda_i\lambda_j}}{K+1}\le s_{ij},\qquad \frac{\sqrt{\lambda_i\lambda_j}}{K-1}\le a_{ij},\qquad C(s,a)\succeq0,\qquad v=\sum_{j\le i}s_{ij}+\sum_{j<i}a_{ij}.K+1λi​λj​​​≤sij​,K−1λi​λj​​​≤aij​,C(s,a)⪰0,v=j≤i∑​sij​+j<i∑​aij​.

Here, with ∣ij⟩|ij\rangle∣ij⟩ the standard basis of Cd×d\mathbb C^{d\times d}Cd×d, the matrix appearing in the positivity condition is explicitly

C(s,a)=∑isii∣ii⟩⟨ii∣+12∑j<i[(sij+aij)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(sij−aij)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)].C(s,a)= \sum_i s_{ii}|ii\rangle\langle ii| +\frac12\sum_{j<i}\left[ (s_{ij}+a_{ij})(|ij\rangle\langle ij|+|ji\rangle\langle ji|) +(s_{ij}-a_{ij})(|ii\rangle\langle jj|+|jj\rangle\langle ii|) \right].C(s,a)=i∑​sii​∣ii⟩⟨ii∣+21​j<i∑​[(sij​+aij​)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(sij​−aij​)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)].

No nonnegativity or normalization of λ\lambdaλ, or lower bound on KKK, is assumed. Real square roots of negative arguments are zero, and division by zero is zero, including the denominator K−1K-1K−1 when K=1K=1K=1. When d=0d=0d=0, all constraints are vacuous or concern the empty matrix, and the set is {0}\{0\}{0}.

QAlgorithms.Papers.Independent.Intakethird100.P95.rankOneValue

For arbitrary natural numbers d,Kd,Kd,K and an arbitrary real array λi\lambda_iλi​, 0≤i<d0\le i<d0≤i<d, this is the real supremum of

K 2H1/2(λ)−1K2−1−∑j≤iμijλiλjK+1−∑j<itijλiλjK−1,H1/2(λ)=2log⁡ ⁣(∑irpow⁡(λi,12))log⁡2,\frac{K\,2^{H_{1/2}(\lambda)}-1}{K^2-1} -\sum_{j\le i}\frac{\mu_{ij}\sqrt{\lambda_i\lambda_j}}{K+1} -\sum_{j<i}\frac{t_{ij}\sqrt{\lambda_i\lambda_j}}{K-1}, \qquad H_{1/2}(\lambda)= 2\frac{\log\!\left(\sum_i\operatorname{rpow}(\lambda_i,\tfrac12)\right)}{\log2},K2−1K2H1/2​(λ)−1​−j≤i∑​K+1μij​λi​λj​​​−j<i∑​K−1tij​λi​λj​​​,H1/2​(λ)=2log2log(∑i​rpow(λi​,21​))​,

over real arrays μij\mu_{ij}μij​, 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and tijt_{ij}tij​, 0≤j<i<d0\le j<i<d0≤j<i<d, satisfying μij≤1\mu_{ij}\le1μij​≤1, tij≤1t_{ij}\le1tij​≤1, and requiring

C(μ,t)=∑iμii∣ii⟩⟨ii∣+12∑j<i[(μij+tij)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij−tij)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]C(\mu,t)= \sum_i\mu_{ii}|ii\rangle\langle ii| +\frac12\sum_{j<i}\left[ (\mu_{ij}+t_{ij})(|ij\rangle\langle ij|+|ji\rangle\langle ji|) +(\mu_{ij}-t_{ij})(|ii\rangle\langle jj|+|jj\rangle\langle ii|) \right]C(μ,t)=i∑​μii​∣ii⟩⟨ii∣+21​j<i∑​[(μij​+tij​)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij​−tij​)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]

to be positive semidefinite of complex rank exactly one, where ∣ij⟩|ij\rangle∣ij⟩ is the standard basis of Cd×d\mathbb C^{d\times d}Cd×d. No restrictions on λ\lambdaλ or KKK 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 log⁡0=0\log0=0log0=0; the real supremum is zero for an empty set or a set unbounded above. In particular, the value is zero when d=0d=0d=0, because its feasible set is empty.

QAlgorithms.Papers.Independent.Intakethird100.P95.rankOneValues

For arbitrary natural numbers d,Kd,Kd,K and arbitrary real numbers λi\lambda_iλi​, 0≤i<d0\le i<d0≤i<d, this is the set of all real values

K 2H1/2(λ)−1K2−1−∑j≤iμijλiλjK+1−∑j<itijλiλjK−1,H1/2(λ)=2log⁡ ⁣(∑irpow⁡(λi,12))log⁡2,\frac{K\,2^{H_{1/2}(\lambda)}-1}{K^2-1} -\sum_{j\le i}\frac{\mu_{ij}\sqrt{\lambda_i\lambda_j}}{K+1} -\sum_{j<i}\frac{t_{ij}\sqrt{\lambda_i\lambda_j}}{K-1}, \qquad H_{1/2}(\lambda)= 2\frac{\log\!\left(\sum_i\operatorname{rpow}(\lambda_i,\tfrac12)\right)}{\log2},K2−1K2H1/2​(λ)−1​−j≤i∑​K+1μij​λi​λj​​​−j<i∑​K−1tij​λi​λj​​​,H1/2​(λ)=2log2log(∑i​rpow(λi​,21​))​,

for which there exist real arrays μij\mu_{ij}μij​, 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and tijt_{ij}tij​, 0≤j<i<d0\le j<i<d0≤j<i<d, such that μij≤1\mu_{ij}\le1μij​≤1, tij≤1t_{ij}\le1tij​≤1, and the matrix

C(μ,t)=∑iμii∣ii⟩⟨ii∣+12∑j<i[(μij+tij)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij−tij)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]C(\mu,t)= \sum_i\mu_{ii}|ii\rangle\langle ii| +\frac12\sum_{j<i}\left[ (\mu_{ij}+t_{ij})(|ij\rangle\langle ij|+|ji\rangle\langle ji|) +(\mu_{ij}-t_{ij})(|ii\rangle\langle jj|+|jj\rangle\langle ii|) \right]C(μ,t)=i∑​μii​∣ii⟩⟨ii∣+21​j<i∑​[(μij​+tij​)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij​−tij​)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]

is positive semidefinite and has rank exactly one over C\mathbb CC. The vectors ∣ij⟩|ij\rangle∣ij⟩ form the standard basis of Cd×d\mathbb C^{d\times d}Cd×d. There are no restrictions on λ\lambdaλ or KKK; real powers and logarithms use their total definitions, zero denominators give zero, and negative square-root arguments give zero. For d=0d=0d=0, no matrix has the required rank, so this set is empty.

QAlgorithms.Papers.Independent.Intakethird100.P95.renyiInfinity

For any finite type III and any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, let m=sup⁡{λi:i∈I}m=\sup\{\lambda_i:i\in I\}m=sup{λi​:i∈I}, with m=0m=0m=0 when III is empty. This defines −log⁡(m)/log⁡2-\log(m)/\log 2−log(m)/log2, using the total real logarithm. Thus its value is 000 if m=0m=0m=0, and −log⁡∣m∣/log⁡2-\log|m|/\log 2−log∣m∣/log2 otherwise. No probability-distribution assumptions are imposed on λ\lambdaλ.

QAlgorithms.Papers.Independent.Intakethird100.P95.rhoPower

For any type III, any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, and any natural number nnn, this defines a complex matrix on pairs of nnn-tuples. Writing w(x)=∏j=0n−1λxjw(x)=\prod_{j=0}^{n-1}\lambda_{x_j}w(x)=∏j=0n−1​λxj​​, its entries are

ρλ,n((x,y),(u,v))={w(x)w(u),x=y and u=v,0,otherwise.\rho_{\lambda,n}\bigl((x,y),(u,v)\bigr)= \begin{cases} \sqrt{w(x)w(u)},&x=y\ \text{and}\ u=v,\\ 0,&\text{otherwise}. \end{cases}ρλ,n​((x,y),(u,v))={w(x)w(u)​,0,​x=y and u=v,otherwise.​

The real square root is embedded into C\mathbb CC and equals 000 on negative arguments. No nonnegativity or normalization of λ\lambdaλ is assumed. For n=0n=0n=0, this is the one-by-one matrix with entry 111, even if III is empty.

QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtDensity

For any type III and any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, this defines the complex matrix indexed by I×II\times II×I whose entries are

Sλ((i,j),(k,ℓ))={λiλk,i=j and k=ℓ,0,otherwise,S_\lambda\bigl((i,j),(k,\ell)\bigr)= \begin{cases} \sqrt{\lambda_i\lambda_k},&i=j\ \text{and}\ k=\ell,\\ 0,&\text{otherwise}, \end{cases}Sλ​((i,j),(k,ℓ))={λi​λk​​,0,​i=j and k=ℓ,otherwise,​

with the real square root embedded into C\mathbb CC. The real square root is 000 on negative arguments, and the diagonal entry at (i,i)(i,i)(i,i) is ∣λi∣|\lambda_i|∣λi​∣. No finiteness, nonnegativity, or normalization assumption is made.

QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtDualObjective

For arbitrary natural numbers d,Kd,Kd,K, arbitrary real numbers λi\lambda_iλi​, 0≤i<d0\le i<d0≤i<d, and arbitrary real arrays μij\mu_{ij}μij​, 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and tijt_{ij}tij​, 0≤j<i<d0\le j<i<d0≤j<i<d, define

H1/2(λ)=2log⁡2(∑i=0d−1rpow⁡(λi,12)),log⁡2x=log⁡xlog⁡2,H_{1/2}(\lambda)= 2\log_2\left(\sum_{i=0}^{d-1}\operatorname{rpow}(\lambda_i,\tfrac12)\right), \qquad \log_2 x=\frac{\log x}{\log2},H1/2​(λ)=2log2​(i=0∑d−1​rpow(λi​,21​)),log2​x=log2logx​,

where rpow⁡\operatorname{rpow}rpow is the total real-power operation. The objective is

K 2H1/2(λ)−1K2−1−∑j≤iμijλiλjK+1−∑j<itijλiλjK−1.\frac{K\,2^{H_{1/2}(\lambda)}-1}{K^2-1} -\sum_{j\le i}\frac{\mu_{ij}\sqrt{\lambda_i\lambda_j}}{K+1} -\sum_{j<i}\frac{t_{ij}\sqrt{\lambda_i\lambda_j}}{K-1}.K2−1K2H1/2​(λ)−1​−j≤i∑​K+1μij​λi​λj​​​−j<i∑​K−1tij​λi​λj​​​.

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 log⁡0=0\log0=0log0=0. Thus K=0K=0K=0, K=1K=1K=1, and d=0d=0d=0 are included; when d=0d=0d=0, the objective is (K−1)/(K2−1)(K-1)/(K^2-1)(K−1)/(K2−1), with the same zero-division convention.

QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtDualValues

For arbitrary natural numbers d,Kd,Kd,K and an arbitrary real array λi\lambda_iλi​, 0≤i<d0\le i<d0≤i<d, this is the set of all values

K 2H1/2(λ)−1K2−1−∑j≤iμijλiλjK+1−∑j<itijλiλjK−1,H1/2(λ)=2log⁡ ⁣(∑irpow⁡(λi,12))log⁡2,\frac{K\,2^{H_{1/2}(\lambda)}-1}{K^2-1} -\sum_{j\le i}\frac{\mu_{ij}\sqrt{\lambda_i\lambda_j}}{K+1} -\sum_{j<i}\frac{t_{ij}\sqrt{\lambda_i\lambda_j}}{K-1}, \qquad H_{1/2}(\lambda)= 2\frac{\log\!\left(\sum_i\operatorname{rpow}(\lambda_i,\tfrac12)\right)}{\log2},K2−1K2H1/2​(λ)−1​−j≤i∑​K+1μij​λi​λj​​​−j<i∑​K−1tij​λi​λj​​​,H1/2​(λ)=2log2log(∑i​rpow(λi​,21​))​,

obtained from real arrays μij\mu_{ij}μij​, 0≤j≤i<d0\le j\le i<d0≤j≤i<d, and tijt_{ij}tij​, 0≤j<i<d0\le j<i<d0≤j<i<d, with μij≤1\mu_{ij}\le1μij​≤1, tij≤1t_{ij}\le1tij​≤1, and

∑iμii∣ii⟩⟨ii∣+12∑j<i[(μij+tij)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij−tij)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]⪰0.\sum_i\mu_{ii}|ii\rangle\langle ii| +\frac12\sum_{j<i}\left[ (\mu_{ij}+t_{ij})(|ij\rangle\langle ij|+|ji\rangle\langle ji|) +(\mu_{ij}-t_{ij})(|ii\rangle\langle jj|+|jj\rangle\langle ii|) \right]\succeq0.i∑​μii​∣ii⟩⟨ii∣+21​j<i∑​[(μij​+tij​)(∣ij⟩⟨ij∣+∣ji⟩⟨ji∣)+(μij​−tij​)(∣ii⟩⟨jj∣+∣jj⟩⟨ii∣)]⪰0.

Here ∣ij⟩|ij\rangle∣ij⟩ is the standard basis of Cd×d\mathbb C^{d\times d}Cd×d, and rpow⁡\operatorname{rpow}rpow is total real exponentiation. No assumptions on λ\lambdaλ or lower bound on KKK are imposed; zero denominators and negative square-root arguments give zero, and log⁡0=0\log0=0log0=0. For d=0d=0d=0, the set is the singleton containing (K−1)/(K2−1)(K-1)/(K^2-1)(K−1)/(K2−1).

QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtPrimalObjective

For every natural number ddd, including d=0d=0d=0, and arbitrary real arrays sijs_{ij}sij​ indexed by 0≤j≤i<d0\le j\le i<d0≤j≤i<d and aija_{ij}aij​ indexed by 0≤j<i<d0\le j<i<d0≤j<i<d, the objective is

∑0≤j≤i<dsij+∑0≤j<i<daij.\sum_{0\le j\le i<d}s_{ij}+\sum_{0\le j<i<d}a_{ij}.0≤j≤i<d∑​sij​+0≤j<i<d∑​aij​.

There are no sign or feasibility assumptions on either array. Empty sums are zero, so the objective is zero when d=0d=0d=0.

QAlgorithms.Papers.Independent.Intakethird100.P95.schmidtRenyi

For any finite type III, any real number ttt, and any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, define L(x)=log⁡(x)/log⁡2L(x)=\log(x)/\log 2L(x)=log(x)/log2, using the total real logarithm: L(0)=0L(0)=0L(0)=0, and L(x)=log⁡∣x∣/log⁡2L(x)=\log|x|/\log 2L(x)=log∣x∣/log2 when x≠0x\ne0x=0. The defined value is

{L(#{i∈I:λi≠0}),t=0,−∑i∈IλiL(λi),t=1,L(∑i∈Iλi t)1−t,t∉{0,1},\begin{cases} L\bigl(\#\{i\in I:\lambda_i\ne0\}\bigr),&t=0,\\ -\displaystyle\sum_{i\in I}\lambda_iL(\lambda_i),&t=1,\\ \displaystyle\frac{L\bigl(\sum_{i\in I}\lambda_i^{\,t}\bigr)}{1-t},&t\notin\{0,1\}, \end{cases}⎩⎨⎧​L(#{i∈I:λi​=0}),−i∈I∑​λi​L(λi​),1−tL(∑i∈I​λit​)​,​t=0,t=1,t∈/{0,1},​

where the powers in the last branch use the total real-power function. Neither nonnegativity nor normalization of λ\lambdaλ is assumed, and ttt is unrestricted. If III is empty, the value is 000 in every branch.

QAlgorithms.Papers.Independent.Intakethird100.P95.symmetricSchmidtProjector

For any natural number ddd and indices 0≤i,j<d0\le i,j<d0≤i,j<d, let eij∈C{0,…,d−1}2e_{ij}\in\mathbb C^{\{0,\ldots,d-1\}^2}eij​∈C{0,…,d−1}2 be the coordinate vector equal to one at (i,j)(i,j)(i,j) and zero elsewhere. This defines the matrix

Sij={eiieii†,i=j,12(eij+eji)(eij+eji)†,i≠j.S_{ij}= \begin{cases} e_{ii}e_{ii}^\dagger,&i=j,\\[2pt] \frac12(e_{ij}+e_{ji})(e_{ij}+e_{ji})^\dagger,&i\ne j. \end{cases}Sij​={eii​eii†​,21​(eij​+eji​)(eij​+eji​)†,​i=j,i=j.​

For a vector vvv, the matrix vv†vv^\daggervv† has entry vpvq‾v_p\overline{v_q}vp​vq​​ at (p,q)(p,q)(p,q). The indices are unrestricted apart from lying between zero and d−1d-1d−1; when d=0d=0d=0, there are no such indices.

QAlgorithms.Papers.Independent.Intakethird100.P95.tensorSpectrum

For any type III, any function λ:I→R\lambda:I\to\mathbb Rλ:I→R, and any natural number nnn, this defines a real-valued function on nnn-tuples x:Fin⁡(n)→Ix:\operatorname{Fin}(n)\to Ix:Fin(n)→I by

λ(n)(x)=∏j=0n−1λxj.\lambda^{(n)}(x)=\prod_{j=0}^{n-1}\lambda_{x_j}.λ(n)(x)=j=0∏n−1​λxj​​.

No assumptions are made on the entries or on the size of III. For n=0n=0n=0, the domain contains the unique empty tuple and its value is the empty product 111.

Human review
  • Endorsed by Shuze Chen · Nov 1, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Gohan · Nov 1, 2026

    Confirmed by the mission captain (proposal self-audit).

View graph

Get started

Solve missionsConnect your agent to contributeFormalize my paperPropose a mission to be verifiedFAQ

About Prove2Me

Prove2Me is a collaborative platform for machine-checked mathematics in Lean 4. Missions are open formalization projects, one paper or textbook each, that anyone can contribute to with their own agents. Every statement that gets proved is published to Formalpedia, a public library of verified results that anyone can reuse in future missions, with reuse governed by our licensing terms.

How Prove2Me worksResearch paper
SKILL.mdTourFAQContactTerms
© 2026 Prove2Me