Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Irreversibility of Asymptotic Entanglement Manipulation Under Quantum Operations Completely Preserving Positivity of Partial Transpose

Definition
QAlgorithms_Papers_Independent_Intakethird100_P90

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

prove2quantum

The objects the formalized statements of Irreversibility of Asymptotic Entanglement Manipulation Under Quantum Operations Completely Preserving Positivity of Partial Transpose (arXiv:1606.09421) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • pptRelativeEntropy: PDF p.2: PPT-relative entropy in bits, with extended Umegaki entropy. Support failure in relativeEntropy gives +infinity, including singular inputs.
  • regularizedPptRelativeEntropy: PDF p.2: infimum over positive tensor powers; n=0 is never divided by.
  • regularizedRains: PDF p.2: regularized Rains bound, using the independent accepted Rains cone.
  • preparationRate: PDF pp.1--2: operational dilution at rate r. The input target(2^m) is m standard Bell pairs in their Schmidt product basis. All sufficiently large positive n admit consumption at most n(r+delta) and error at most epsilon. A full trace norm is used; its scale does not change vanishing-error rates.
  • asymptoticEntanglementCost: PDF pp.1--2: infimum of operational Bell-pair consumption rates. Empty rate sets give +infinity, not an artificial real zero. On the theorem's finite-dimensional density inputs the value is finite and nonnegative.
  • loccDistillationRate: PDF pp.1--2,5: operational two-way LOCC Bell-pair yield, with all branches retained and arbitrarily small error. No bounded number of rounds is imposed.
  • loccDistillable: PDF pp.1--2,5: supremum of vanishing-error LOCC distillation rates.
  • etaFeasible: PDF p.3, Eq.(4): both support sandwich inequalities, in Loewner order.
  • supportEta: PDF p.3, Eq.(4): infimum of the partially transposed spectral norm over the support sandwich. For a density matrix the set is nonempty and its minimum is strictly positive; eta_sdp_duality below explicitly includes this fact and attainment.
  • etaEntanglement: PDF pp.3,6: the additive support bound in bits. Used only for normalized states, whose nonzero support ensures supportEta>0, so log(0) is never a physical value.
  • etaPrimalFeasible: PDF p.6, Eq.(12): primal SDP with two two-sided operator sandwiches.
  • etaPrimal: PDF p.6, Eq.(12): primal value, defined independently of the dual.
  • etaDualFeasible: PDF p.6, Eq.(13): four PSD dual variables and the printed two constraints.
  • etaDualObjective: PDF p.6, Eq.(13): real trace pairing with the support projector.
  • etaDual: PDF p.6, Eq.(13): dual supremum, not defined to equal the primal.
  • supportOverlap: PDF pp.3,7, Eq.(16): PPT-state overlap with the actual support projection.
  • supportDualFeasible: PDF p.7, Eqs.(16),(18): exactly the two one-sided constraints.
  • supportDual: PDF p.7, Eq.(16): one-sided dual infimum, independently defined.
  • wangDuanBound: PDF p.4: E_W, the improved logarithmic-negativity upper bound from Ref.[30]. X>=rho is Loewner order; no normalization of X is imposed. On density inputs traceNorm(gamma X)>=Tr X>=1, so logarithms have positive arguments.
  • antisymmetricSupported: PDF pp.1--2,7: range contained in the two-qutrit antisymmetric subspace. Each column has swap eigenvalue -1, which is exactly that range inclusion.
  • stateRank: PDF pp.1--2: rank is the complex dimension of the matrix's linear-map range.
  • antisymmetricKet: PDF p.2: the normalized antisymmetric ket on basis coordinates 0 and j.
  • rhoV: PDF p.2: rho_v is the specified equally weighted mixture, not a free state.
  • exteriorKet: PDF pp.4,7: arbitrary complex exterior-product spectral vector. Only its unit norm and orthogonality to u1 are assumed in Proposition 5.
  • alphaKet: PDF p.4: the two disjoint-support unit vectors in the alpha family.
  • alphaState: PDF pp.4--5: rho^(alpha); theorems use the exact interval [21/50,1/2].
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P36
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P32
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P06

/-!
arXiv:1606.09421v3 (Wang--Duan), PDF pp.1--7. All matrices use fixed
orthonormal product bases. Only pinned accepted-library modules are imported;
the independent-paper workflow permits their reuse. New objects stay here.
The source's PPT class is completely PPT-preserving CPTP, and its operational
currency is the normalized Bell pair. Rates use deterministic unconditional
outputs and vanishing trace-norm error (equivalently, fidelity tending to one).
No assertion is made about the separable infinite-dimensional interfaces in
the motivating supporting-paper gap list.
-/

noncomputable section

open scoped Matrix ComplexConjugate ComplexOrder

namespace QAlgorithms.Papers.Independent.Intakethird100.P90

attribute [local instance] Classical.propDecidable Classical.decEq

open QAlgorithms
open QAlgorithms.Papers.Independent.Dependency50.P36
  (relativeEntropy rains target canonicalMatrix allowed pptDistillable)
open QAlgorithms.Papers.Independent.Information100.P04 (gamma copies IsPPT)
open QAlgorithms.Papers.Independent.Next46.P42 (regroupProduct)
open QAlgorithms.Papers.Independent.Information100.P06 (supportProjection)
open QAlgorithms.Papers.Independent.Information100.P21 (FiniteLOCC)
open QAlgorithms.Papers.Independent.Gap48.P32 (antisymmetricQutrit)

/-- PDF p.2: PPT-relative entropy in bits, with extended Umegaki entropy.
Support failure in relativeEntropy gives +infinity, including singular inputs. -/
def pptRelativeEntropy {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : EReal :=
  sInf {e | ∃ sigma : Matrix (A × B) (A × B) ℂ,
    IsPPT sigma ∧ e = relativeEntropy rho sigma}

/-- PDF p.2: infimum over positive tensor powers; n=0 is never divided by. -/
def regularizedPptRelativeEntropy {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : EReal :=
  ⨅ n : ℕ, pptRelativeEntropy (copies rho (n+1)) / ((n+1 : ℕ) : EReal)

/-- PDF p.2: regularized Rains bound, using the independent accepted Rains cone. -/
def regularizedRains {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : EReal :=
  ⨅ n : ℕ, rains (copies rho (n+1)) / ((n+1 : ℕ) : EReal)

/-- PDF pp.1--2: operational dilution at rate r. The input target(2^m) is m
standard Bell pairs in their Schmidt product basis. All sufficiently large
positive n admit consumption at most n(r+delta) and error at most epsilon.
A full trace norm is used; its scale does not change vanishing-error rates. -/
def preparationRate {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (locc : Bool)
    (rho : Matrix (A × B) (A × B) ℂ) (r : ℝ) : Prop :=
  0 ≤ r ∧ ∀ epsilon > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
    ∀ n : ℕ, N ≤ n → 0 < n → ∃ m : ℕ,
      (m : ℝ) ≤ (n : ℝ) * (r + delta) ∧
      ∃ Phi : Superop (Fin (2^m) × Fin (2^m))
        (Fin (Fintype.card (Fin n → A)) × Fin (Fintype.card (Fin n → B))),
        allowed (if locc then .locc else .ppt) Phi ∧
        traceNorm (Phi (target (2^m)) - canonicalMatrix (copies rho n)) ≤ epsilon

/-- PDF pp.1--2: infimum of operational Bell-pair consumption rates.
Empty rate sets give +infinity, not an artificial real zero. On the theorem's
finite-dimensional density inputs the value is finite and nonnegative. -/
def asymptoticEntanglementCost {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (locc : Bool)
    (rho : Matrix (A × B) (A × B) ℂ) : EReal :=
  sInf {e | ∃ r : ℝ, preparationRate locc rho r ∧ e = (r : EReal)}

/-- PDF pp.1--2,5: operational two-way LOCC Bell-pair yield, with all branches
retained and arbitrarily small error. No bounded number of rounds is imposed. -/
def loccDistillationRate {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) (r : ℝ) : Prop :=
  0 ≤ r ∧ ∀ epsilon > (0 : ℝ), ∀ delta > (0 : ℝ), ∃ N : ℕ,
    ∀ n : ℕ, N ≤ n → 0 < n → ∃ m : ℕ,
      (n : ℝ) * (r - delta) ≤ (m : ℝ) ∧
      ∃ Phi : Superop
        (Fin (Fintype.card (Fin n → A)) × Fin (Fintype.card (Fin n → B)))
        (Fin (2^m) × Fin (2^m)),
        FiniteLOCC Phi ∧
        traceNorm (Phi (canonicalMatrix (copies rho n)) - target (2^m)) ≤ epsilon

/-- PDF pp.1--2,5: supremum of vanishing-error LOCC distillation rates. -/
def loccDistillable {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : EReal :=
  sSup {e | ∃ r : ℝ, loccDistillationRate rho r ∧ e = (r : EReal)}

/-- PDF p.3, Eq.(4): both support sandwich inequalities, in Loewner order. -/
def etaFeasible {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho Y : Matrix (A × B) (A × B) ℂ) : Prop :=
  Y.IsHermitian ∧ (Y - gamma (supportProjection rho)).PosSemidef ∧
    (Y + gamma (supportProjection rho)).PosSemidef

/-- PDF p.3, Eq.(4): infimum of the partially transposed spectral norm over the
support sandwich. For a density matrix the set is nonempty and its minimum is
strictly positive; eta_sdp_duality below explicitly includes this fact and attainment. -/
def supportEta {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  sInf {t | ∃ Y, etaFeasible rho Y ∧ t = specNorm (gamma Y)}

/-- PDF pp.3,6: the additive support bound in bits. Used only for normalized
states, whose nonzero support ensures supportEta>0, so log(0) is never a physical value. -/
def etaEntanglement {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  -Real.logb 2 (supportEta rho)

/-- PDF p.6, Eq.(12): primal SDP with two two-sided operator sandwiches. -/
def etaPrimalFeasible {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ)
    (t : ℝ) (Y : Matrix (A × B) (A × B) ℂ) : Prop :=
  etaFeasible rho Y ∧
    ((t : ℂ) • (1 : Matrix (A × B) (A × B) ℂ) - gamma Y).PosSemidef ∧
    (gamma Y + (t : ℂ) • (1 : Matrix (A × B) (A × B) ℂ)).PosSemidef

/-- PDF p.6, Eq.(12): primal value, defined independently of the dual. -/
def etaPrimal {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  sInf {t | ∃ Y, etaPrimalFeasible rho t Y}

/-- PDF p.6, Eq.(13): four PSD dual variables and the printed two constraints. -/
def etaDualFeasible {A B : Type} [Fintype A] [Fintype B]
    (V F W X : Matrix (A × B) (A × B) ℂ) : Prop :=
  V.PosSemidef ∧ F.PosSemidef ∧ W.PosSemidef ∧ X.PosSemidef ∧
    (gamma (W-X) - (V+F)).PosSemidef ∧ (W+X).trace.re ≤ 1

/-- PDF p.6, Eq.(13): real trace pairing with the support projector. -/
def etaDualObjective {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B]
    (rho V F : Matrix (A × B) (A × B) ℂ) : ℝ :=
  (supportProjection rho * gamma (V-F)).trace.re

/-- PDF p.6, Eq.(13): dual supremum, not defined to equal the primal. -/
def etaDual {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  sSup {v | ∃ V F W X : Matrix (A × B) (A × B) ℂ,
    etaDualFeasible V F W X ∧ v = etaDualObjective rho V F}

/-- PDF pp.3,7, Eq.(16): PPT-state overlap with the actual support projection. -/
def supportOverlap {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  sSup {v | ∃ sigma : Matrix (A × B) (A × B) ℂ,
    IsPPT sigma ∧ v = (supportProjection rho * sigma).trace.re}

/-- PDF p.7, Eqs.(16),(18): exactly the two one-sided constraints. -/
def supportDualFeasible {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ)
    (t : ℝ) (Y : Matrix (A × B) (A × B) ℂ) : Prop :=
  Y.IsHermitian ∧ (Y - gamma (supportProjection rho)).PosSemidef ∧
    ((t : ℂ) • (1 : Matrix (A × B) (A × B) ℂ) - gamma Y).PosSemidef

/-- PDF p.7, Eq.(16): one-sided dual infimum, independently defined. -/
def supportDual {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  sInf {t | ∃ Y, supportDualFeasible rho t Y}

/-- PDF p.4: E_W, the improved logarithmic-negativity upper bound from Ref.[30].
X>=rho is Loewner order; no normalization of X is imposed. On density inputs
traceNorm(gamma X)>=Tr X>=1, so logarithms have positive arguments. -/
def wangDuanBound {A B : Type} [Fintype A] [Fintype B]
    [DecidableEq A] [DecidableEq B] (rho : Matrix (A × B) (A × B) ℂ) : ℝ :=
  sInf {e | ∃ X : Matrix (A × B) (A × B) ℂ,
    X.IsHermitian ∧ (X-rho).PosSemidef ∧ e = Real.logb 2 (traceNorm (gamma X))}

/-- PDF pp.1--2,7: range contained in the two-qutrit antisymmetric subspace.
Each column has swap eigenvalue -1, which is exactly that range inclusion. -/
def antisymmetricSupported (rho : Matrix (Fin 3 × Fin 3) (Fin 3 × Fin 3) ℂ) : Prop :=
  ∀ i j, rho (i.2,i.1) j = -rho i j

/-- PDF pp.1--2: rank is the complex dimension of the matrix's linear-map range. -/
def stateRank {I : Type} [Fintype I] [DecidableEq I] (rho : Matrix I I ℂ) : ℕ :=
  Module.finrank ℂ (LinearMap.range rho.mulVecLin)

/-- PDF p.2: the normalized antisymmetric ket on basis coordinates 0 and j. -/
def antisymmetricKet (j : Fin 3) : EuclideanSpace ℂ (Fin 3 × Fin 3) :=
  ((Real.sqrt 2 : ℂ)⁻¹) • (ket (0,j) - ket (j,0))

/-- PDF p.2: rho_v is the specified equally weighted mixture, not a free state. -/
def rhoV : Matrix (Fin 3 × Fin 3) (Fin 3 × Fin 3) ℂ :=
  (1/2 : ℂ) • (pureDensity (antisymmetricKet 1) + pureDensity (antisymmetricKet 2))

/-- PDF pp.4,7: arbitrary complex exterior-product spectral vector. Only its
unit norm and orthogonality to u1 are assumed in Proposition 5. -/
def exteriorKet (a b : EuclideanSpace ℂ (Fin 3)) : EuclideanSpace ℂ (Fin 3 × Fin 3) :=
  WithLp.toLp 2 fun i => (a i.1 * b i.2 - b i.1 * a i.2) / (Real.sqrt 2 : ℂ)

/-- PDF p.4: the two disjoint-support unit vectors in the alpha family. -/
def alphaKet (alpha : ℝ) (j : Fin 3) : EuclideanSpace ℂ (Fin 3 × Fin 3) :=
  (Real.sqrt alpha : ℂ) • ket (0,j) - (Real.sqrt (1-alpha) : ℂ) • ket (j,0)

/-- PDF pp.4--5: rho^(alpha); theorems use the exact interval [21/50,1/2]. -/
def alphaState (alpha : ℝ) : Matrix (Fin 3 × Fin 3) (Fin 3 × Fin 3) ℂ :=
  (1/2 : ℂ) • (pureDensity (alphaKet alpha 1) + pureDensity (alphaKet alpha 2))

/-- Product-basis check: Bob's transpose sends |01><12| to |02><11|. -/
example (X : Matrix (Fin 3 × Fin 3) (Fin 3 × Fin 3) ℂ) :
    gamma X (0,1) (1,2) = X (0,2) (1,1) := rfl

/-- The Bell currency has the correct 00:11 coherence and zero cross-sector weight. -/
example : target 2 (0,0) (1,1) = (1/2 : ℂ) ∧ target 2 (0,1) (0,1) = 0 := by
  norm_num [target]

/-- Zero consumed Bell pairs mean a scalar trace-one state, not a zero matrix. -/
example : target (2^0) (0,0) (0,0) = (1 : ℂ) := by norm_num [target]

/-- The tensor power at block length zero is the scalar identity. -/
example {A B : Type} (rho : Matrix (A × B) (A × B) ℂ)
    (i j : (Fin 0 → A) × (Fin 0 → B)) : copies rho 0 i j = 1 := by simp [copies]

/-- The finite zero state has zero support projection. Physical theorems exclude it. -/
example {I : Type} [Fintype I] [DecidableEq I] :
    supportProjection (0 : Matrix I I ℂ) = 0 := by
  simp [supportProjection, cfc_zero]

/-- The support sandwich orientation is the paper's -Y<=P^TB<=Y. -/
example {A B : Type} [Fintype A] [Fintype B] [DecidableEq A] [DecidableEq B]
    (rho Y : Matrix (A × B) (A × B) ℂ) (h : etaFeasible rho Y) :
    (Y - gamma (supportProjection rho)).PosSemidef ∧
      (Y + gamma (supportProjection rho)).PosSemidef := ⟨h.2.1,h.2.2⟩

/-- Regrouping keeps both input coherences. -/
example {A B C D : Type} (rho : Matrix (A × B) (A × B) ℂ)
    (sigma : Matrix (C × D) (C × D) ℂ) (a a1 : A) (b b1 : B) (c c1 : C) (d d1 : D) :
    regroupProduct rho sigma ((a,c),(b,d)) ((a1,c1),(b1,d1)) =
      rho (a,b) (a1,b1) * sigma (c,d) (c1,d1) := rfl

/-- Antisymmetric means that equal local factors give the zero exterior vector. -/
example (a : EuclideanSpace ℂ (Fin 3)) : exteriorKet a a = 0 := by
  ext i
  simp [exteriorKet]

/-- The accepted rank-three example is (I-swap)/6, with the correct sign. -/
example : antisymmetricQutrit (0,1) (1,0) = -(1/6 : ℂ) := by
  norm_num [antisymmetricQutrit]

end QAlgorithms.Papers.Independent.Intakethird100.P90
Source
Irreversibility of Asymptotic Entanglement Manipulation Under Quantum Operations Completely Preserving Positivity of Partial Transpose, arXiv:1606.09421v3, https://arxiv.org/abs/1606.09421v3
Read-back

What the Lean code literally says, in plain math · gpt-6.1-sol

QAlgorithms.Papers.Independent.Intakethird100.P90.alphaKet

For every real number α\alphaα and every j∈{0,1,2}j\in\{0,1,2\}j∈{0,1,2}, this defines a vector in C{0,1,2}×{0,1,2}\mathbb C^{\{0,1,2\}\times\{0,1,2\}}C{0,1,2}×{0,1,2} by

vα,j=α ∣0,j⟩−1−α ∣j,0⟩,v_{\alpha,j}=\sqrt{\alpha}\,|0,j\rangle-\sqrt{1-\alpha}\,|j,0\rangle,vα,j​=α​∣0,j⟩−1−α​∣j,0⟩,

where ∣a,b⟩|a,b\rangle∣a,b⟩ has coordinate 111 at (a,b)(a,b)(a,b) and coordinate 000 everywhere else. There is no restriction that α∈[0,1]\alpha\in[0,1]α∈[0,1]; the real square root used here equals 000 on negative inputs. The index j=0j=0j=0 is included, in which case the two basis vectors coincide and vα,0=(α−1−α)∣0,0⟩v_{\alpha,0}=(\sqrt{\alpha}-\sqrt{1-\alpha})|0,0\ranglevα,0​=(α​−1−α​)∣0,0⟩. No unit-norm condition is part of this definition.

QAlgorithms.Papers.Independent.Intakethird100.P90.alphaState

For every real number α\alphaα, this defines the complex matrix

ρα=12(vα,1vα,1†+vα,2vα,2†),vα,j=α ∣0,j⟩−1−α ∣j,0⟩,\rho_\alpha=\frac12\left(v_{\alpha,1}v_{\alpha,1}^{\dagger}+v_{\alpha,2}v_{\alpha,2}^{\dagger}\right), \qquad v_{\alpha,j}=\sqrt{\alpha}\,|0,j\rangle-\sqrt{1-\alpha}\,|j,0\rangle,ρα​=21​(vα,1​vα,1†​+vα,2​vα,2†​),vα,j​=α​∣0,j⟩−1−α​∣j,0⟩,

indexed by pairs from {0,1,2}\{0,1,2\}{0,1,2}. Here ∣a,b⟩|a,b\rangle∣a,b⟩ is the coordinate basis vector, and vv†vv^\daggervv† means the matrix with entry vuvv‾v_u\overline{v_v}vu​vv​​ at row uuu, column vvv. Equivalently,

ρα=12∑j=12[(α)2∣0,j⟩⟨0,j∣+(1−α)2∣j,0⟩⟨j,0∣−α1−α(∣0,j⟩⟨j,0∣+∣j,0⟩⟨0,j∣)].\rho_\alpha=\frac12\sum_{j=1}^{2}\left[ (\sqrt{\alpha})^2|0,j\rangle\langle0,j| +(\sqrt{1-\alpha})^2|j,0\rangle\langle j,0| -\sqrt{\alpha}\sqrt{1-\alpha}\bigl(|0,j\rangle\langle j,0|+|j,0\rangle\langle0,j|\bigr) \right].ρα​=21​j=1∑2​[(α​)2∣0,j⟩⟨0,j∣+(1−α​)2∣j,0⟩⟨j,0∣−α​1−α​(∣0,j⟩⟨j,0∣+∣j,0⟩⟨0,j∣)].

The parameter α\alphaα is unrestricted, including values outside [0,1][0,1][0,1], and real square roots of negative numbers are 000. The definition itself imposes no trace-one hypothesis.

QAlgorithms.Papers.Independent.Intakethird100.P90.antisymmetricKet

For every j∈{0,1,2}j\in\{0,1,2\}j∈{0,1,2}, this definition gives the vector in C{0,1,2}2\mathbb C^{\{0,1,2\}^2}C{0,1,2}2

vj=e(0,j)−e(j,0)2,v_j=\frac{e_{(0,j)}-e_{(j,0)}}{\sqrt2},vj​=2​e(0,j)​−e(j,0)​​,

where e(a,b)e_{(a,b)}e(a,b)​ has coordinate one at (a,b)(a,b)(a,b) and zero elsewhere. The allowed value j=0j=0j=0 gives the zero vector; the definition imposes no condition excluding it.

QAlgorithms.Papers.Independent.Intakethird100.P90.antisymmetricSupported

For an arbitrary complex matrix ρ\rhoρ indexed by {0,1,2}×{0,1,2}\{0,1,2\}\times\{0,1,2\}{0,1,2}×{0,1,2}, this predicate requires

ρ(b,a),j=−ρ(a,b),j\rho_{(b,a),j}=-\rho_{(a,b),j}ρ(b,a),j​=−ρ(a,b),j​

for every a,b∈{0,1,2}a,b\in\{0,1,2\}a,b∈{0,1,2} and every column index j∈{0,1,2}2j\in\{0,1,2\}^2j∈{0,1,2}2. The condition swaps only the two components of the row index. In particular, every row indexed by (a,a)(a,a)(a,a) must be zero. No positivity, Hermiticity, trace normalization, or separate column-swap condition is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P90.asymptoticEntanglementCost

For finite types A,BA,BA,B with decidable equality, a Boolean selector bbb, and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the extended real infimum of all real r≥0r\ge0r≥0 with the following property: for every ϵ,δ>0\epsilon,\delta>0ϵ,δ>0, all sufficiently large positive integers nnn admit m∈Nm\in\mathbb Nm∈N and a map Φ\PhiΦ such that

m≤n(r+δ),∥Φ(T2m)−ρ[n]^∥1≤ϵ.m\le n(r+\delta), \qquad \left\|\Phi(T_{2^m})-\widehat{\rho^{[n]}}\right\|_1\le\epsilon.m≤n(r+δ),​Φ(T2m​)−ρ[n]​​1​≤ϵ.

The input matrix TkT_kTk​ has entry 1/k1/k1/k when both its row and column pairs have equal components, and zero otherwise; ρ[n]\rho^{[n]}ρ[n] has entries ∏t=0n−1ρ(at,bt),(at′,bt′)\prod_{t=0}^{n-1}\rho_{(a_t,b_t),(a'_t,b'_t)}∏t=0n−1​ρ(at​,bt​),(at′​,bt′​)​, and the hat reindexes An,BnA^n,B^nAn,Bn by their chosen finite ordinal enumerations. The norm is Re⁡Tr⁡X†X\operatorname{Re}\operatorname{Tr}\sqrt{X^\dagger X}ReTrX†X​. If bbb is true, Φ\PhiΦ must be built recursively from identities and complete finite local Kraus instruments on either party, with an outcome-dependent recursive continuation and summed branches. If bbb is false, Φ\PhiΦ and its conjugation by the second-register partial transpose Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​ must both have finite Kraus representations X↦∑jKjXKj†X\mapsto\sum_jK_jXK_j^\daggerX↦∑j​Kj​XKj†​ with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. The quantifiers allow m=0m=0m=0, but require n>0n>0n>0; no density-matrix condition is imposed on ρ\rhoρ, and A,BA,BA,B may be empty. If there are no such rates, the infimum is +∞+\infty+∞; attainment is not asserted.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaDual

For finite types A,BA,BA,B equipped with decidable equality and an arbitrary complex matrix ρ\rhoρ indexed by A×BA\times BA×B, this definition is the real supremum

sup⁡{Re⁡Tr⁡ ⁣(Pρ Γ(V−F)):V,F,W,X⪰0, Γ(W−X)−(V+F)⪰0, Re⁡Tr⁡(W+X)≤1},\sup\left\{ \operatorname{Re}\operatorname{Tr}\!\left(P_\rho\,\Gamma(V-F)\right): V,F,W,X\succeq0,\ \Gamma(W-X)-(V+F)\succeq0,\ \operatorname{Re}\operatorname{Tr}(W+X)\le1 \right\},sup{ReTr(Pρ​Γ(V−F)):V,F,W,X⪰0, Γ(W−X)−(V+F)⪰0, ReTr(W+X)≤1},

where all four optimization variables are complex matrices on A×BA\times BA×B, (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​, and Pρ=cfc⁡(f,ρ)P_\rho=\operatorname{cfc}(f,\rho)Pρ​=cfc(f,ρ) for f(0)=0f(0)=0f(0)=0 and f(t)=1f(t)=1f(t)=1 when t≠0t\ne0t=0, using real continuous functional calculus. No state or Hermiticity assumption is imposed on ρ\rhoρ. The all-zero choice is feasible and gives objective zero. Empty index types are included and give value zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaDualFeasible

For finite types A,BA,BA,B, possibly empty, and complex matrices V,F,W,XV,F,W,XV,F,W,X indexed by A×BA\times BA×B, this predicate means that all four matrices are positive semidefinite and

Γ(W−X)−(V+F)⪰0,Re⁡Tr⁡(W+X)≤1,\Gamma(W-X)-(V+F)\succeq0, \qquad \operatorname{Re}\operatorname{Tr}(W+X)\le1,Γ(W−X)−(V+F)⪰0,ReTr(W+X)≤1,

where the operation Γ\GammaΓ is defined entrywise by (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​. Here ⪰0\succeq0⪰0 denotes positive semidefiniteness. There are no trace-one requirements on these matrices; in particular, taking all four to be zero satisfies the predicate.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaDualObjective

For finite types A,BA,BA,B equipped with decidable equality, and arbitrary complex matrices ρ,V,F\rho,V,Fρ,V,F indexed by A×BA\times BA×B, this definition returns

Re⁡Tr⁡ ⁣(Pρ Γ(V−F)),\operatorname{Re}\operatorname{Tr}\!\left(P_\rho\,\Gamma(V-F)\right),ReTr(Pρ​Γ(V−F)),

where Pρ=cfc⁡(f,ρ)P_\rho=\operatorname{cfc}(f,\rho)Pρ​=cfc(f,ρ), with f(0)=0f(0)=0f(0)=0 and f(t)=1f(t)=1f(t)=1 for every nonzero real ttt, using the real continuous functional calculus, and (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​. No density, positivity, Hermiticity, or feasibility hypothesis is imposed on ρ,V,F\rho,V,Fρ,V,F. Empty index types are permitted, in which case the trace is zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaEntanglement

For finite types A,BA,BA,B with decidable equality and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the real number

−log⁡2 ⁣(inf⁡Y=Y†Y−Γ(Pρ)⪰0Y+Γ(Pρ)⪰0∥Γ(Y)∥op).-\log_2\!\left( \inf_{\substack{Y=Y^\dagger\\ Y-\Gamma(P_\rho)\succeq0\\ Y+\Gamma(P_\rho)\succeq0}} \|\Gamma(Y)\|_{\mathrm{op}} \right).−log2​​Y=Y†Y−Γ(Pρ​)⪰0Y+Γ(Pρ​)⪰0​inf​∥Γ(Y)∥op​​.

Here Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, the norm is the Euclidean operator norm, and Pρ=f(ρ)P_\rho=f(\rho)Pρ​=f(ρ) uses real continuous functional calculus with f(0)=0f(0)=0f(0)=0 and f(x)=1f(x)=1f(x)=1 for x≠0x\ne0x=0; it selects nonzero eigenspaces for Hermitian ρ\rhoρ and gives zero for non-Hermitian ρ\rhoρ. No positivity or normalization is assumed for ρ\rhoρ. The logarithm is the totalized real base-two logarithm, so its value at zero is zero: in particular, ρ=0\rho=0ρ=0 and empty underlying types give value zero, rather than an infinite value.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaFeasible

For finite types A,BA,BA,B with decidable equality and any complex matrices ρ,Y\rho,Yρ,Y on A×BA\times BA×B, this predicate means

Y=Y†,Y−Γ(Pρ)⪰0,Y+Γ(Pρ)⪰0,Y=Y^\dagger,\qquad Y-\Gamma(P_\rho)\succeq0,\qquad Y+\Gamma(P_\rho)\succeq0,Y=Y†,Y−Γ(Pρ​)⪰0,Y+Γ(Pρ​)⪰0,

where Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​ and Pρ=f(ρ)P_\rho=f(\rho)Pρ​=f(ρ) is defined by real continuous functional calculus with f(0)=0f(0)=0f(0)=0 and f(x)=1f(x)=1f(x)=1 for x≠0x\ne0x=0. For Hermitian ρ\rhoρ, PρP_\rhoPρ​ projects onto its nonzero eigenspaces; for non-Hermitian ρ\rhoρ, this real functional calculus gives zero. No positivity or normalization hypothesis is imposed on ρ\rhoρ. Empty AAA or BBB are permitted, in which case all three conditions hold for the unique matrices.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaPrimal

For finite types A,BA,BA,B with decidable equality and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the real infimum

inf⁡{t∈R:∃Y, Y=Y†, Y±Γ(Pρ)⪰0, tI±Γ(Y)⪰0},\inf\left\{t\in\mathbb R: \exists Y,\ Y=Y^\dagger,\ Y\pm\Gamma(P_\rho)\succeq0,\ tI\pm\Gamma(Y)\succeq0 \right\},inf{t∈R:∃Y, Y=Y†, Y±Γ(Pρ​)⪰0, tI±Γ(Y)⪰0},

where each ±\pm± requires both signs, III is the identity, Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, and Pρ=f(ρ)P_\rho=f(\rho)Pρ​=f(ρ) uses real continuous functional calculus with f(0)=0f(0)=0f(0)=0 and f(x)=1f(x)=1f(x)=1 for x≠0x\ne0x=0. For Hermitian ρ\rhoρ, PρP_\rhoPρ​ projects onto its nonzero eigenspaces; for non-Hermitian ρ\rhoρ, it is zero. No density-matrix assumption is imposed on ρ\rhoρ, and the defining set does not separately restrict ttt to be nonnegative. If ρ=0\rho=0ρ=0, the value is zero. If either underlying type is empty, every real ttt is feasible, and the real infimum operation assigns zero to this unbounded-below set.

QAlgorithms.Papers.Independent.Intakethird100.P90.etaPrimalFeasible

For finite types A,BA,BA,B with decidable equality, arbitrary complex matrices ρ,Y\rho,Yρ,Y on A×BA\times BA×B, and any real ttt, this predicate means

Y=Y†,Y−Γ(Pρ)⪰0,Y+Γ(Pρ)⪰0,tI−Γ(Y)⪰0,Γ(Y)+tI⪰0.Y=Y^\dagger,\qquad Y-\Gamma(P_\rho)\succeq0,\qquad Y+\Gamma(P_\rho)\succeq0,\qquad tI-\Gamma(Y)\succeq0,\qquad \Gamma(Y)+tI\succeq0.Y=Y†,Y−Γ(Pρ​)⪰0,Y+Γ(Pρ​)⪰0,tI−Γ(Y)⪰0,Γ(Y)+tI⪰0.

Here III is the identity matrix, Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, and Pρ=f(ρ)P_\rho=f(\rho)Pρ​=f(ρ) is real continuous functional calculus with f(0)=0f(0)=0f(0)=0 and f(x)=1f(x)=1f(x)=1 for x≠0x\ne0x=0. It projects onto nonzero eigenspaces when ρ\rhoρ is Hermitian and is zero otherwise. There is no explicit assumption t≥0t\ge0t≥0 and no positivity or normalization assumption on ρ\rhoρ. If either underlying type is empty, the conditions hold for every real ttt.

QAlgorithms.Papers.Independent.Intakethird100.P90.exteriorKet

For arbitrary vectors a,b∈C{0,1,2}a,b\in\mathbb C^{\{0,1,2\}}a,b∈C{0,1,2}, this definition gives the vector v∈C{0,1,2}2v\in\mathbb C^{\{0,1,2\}^2}v∈C{0,1,2}2 with coordinates

v(i,j)=aibj−biaj2(i,j∈{0,1,2}).v_{(i,j)}=\frac{a_i b_j-b_i a_j}{\sqrt2} \qquad(i,j\in\{0,1,2\}).v(i,j)​=2​ai​bj​−bi​aj​​(i,j∈{0,1,2}).

There are no normalization, orthogonality, or nonzero assumptions on a,ba,ba,b; in particular, linearly dependent inputs give the zero vector.

QAlgorithms.Papers.Independent.Intakethird100.P90.loccDistillable

For finite types A,BA,BA,B with decidable equality and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the extended real supremum of all real r≥0r\ge0r≥0 such that, for every ϵ,δ>0\epsilon,\delta>0ϵ,δ>0, every sufficiently large positive integer nnn admits m∈Nm\in\mathbb Nm∈N and a finite recursive local-operation map Φ\PhiΦ with

n(r−δ)≤m,∥Φ(ρ[n]^)−T2m∥1≤ϵ.n(r-\delta)\le m, \qquad \left\|\Phi(\widehat{\rho^{[n]}})-T_{2^m}\right\|_1\le\epsilon.n(r−δ)≤m,​Φ(ρ[n]​)−T2m​​1​≤ϵ.

Here ρ(a,b),(a′,b′)[n]=∏t=0n−1ρ(at,bt),(at′,bt′)\rho^{[n]}_{(a,b),(a',b')}=\prod_{t=0}^{n-1}\rho_{(a_t,b_t),(a'_t,b'_t)}ρ(a,b),(a′,b′)[n]​=∏t=0n−1​ρ(at​,bt​),(at′​,bt′​)​; the hat reindexes An,BnA^n,B^nAn,Bn through their chosen finite ordinal enumerations; TkT_kTk​ has entry 1/k1/k1/k when both index pairs have equal components and zero otherwise; and ∥X∥1=Re⁡Tr⁡X†X\|X\|_1=\operatorname{Re}\operatorname{Tr}\sqrt{X^\dagger X}∥X∥1​=ReTrX†X​. A permitted map is an identity or is obtained by applying a finite complete local Kraus family on either party, following each outcome by another permitted map, and summing branches; completeness means ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. No density-matrix assumption is imposed on ρ\rhoρ, empty underlying types are included, and m=0m=0m=0 is allowed. An empty set of achievable rates has supremum −∞-\infty−∞, while an unbounded set has supremum +∞+\infty+∞; attainment is not asserted.

QAlgorithms.Papers.Independent.Intakethird100.P90.loccDistillationRate

For finite types A,BA,BA,B with decidable equality, any complex matrix ρ\rhoρ on A×BA\times BA×B, and r∈Rr\in\mathbb Rr∈R, this predicate means that r≥0r\ge0r≥0 and, for every ϵ,δ>0\epsilon,\delta>0ϵ,δ>0, there exists N∈NN\in\mathbb NN∈N such that every integer n≥Nn\ge Nn≥N with n>0n>0n>0 admits m∈Nm\in\mathbb Nm∈N and a map Φ\PhiΦ satisfying

n(r−δ)≤m,∥Φ(ρ[n]^)−T2m∥1≤ϵ.n(r-\delta)\le m, \qquad \left\|\Phi(\widehat{\rho^{[n]}})-T_{2^m}\right\|_1\le\epsilon.n(r−δ)≤m,​Φ(ρ[n]​)−T2m​​1​≤ϵ.

The input is indexed by finite ordinal enumerations of AnA^nAn and BnB^nBn, and the output by two registers of size 2m2^m2m. Here ρ(a,b),(a′,b′)[n]=∏t=0n−1ρ(at,bt),(at′,bt′)\rho^{[n]}_{(a,b),(a',b')}=\prod_{t=0}^{n-1}\rho_{(a_t,b_t),(a'_t,b'_t)}ρ(a,b),(a′,b′)[n]​=∏t=0n−1​ρ(at​,bt​),(at′​,bt′​)​, the hat denotes reindexing by the chosen bijections, and (Tk)(i,j),(i′,j′)=1/k(T_k)_{(i,j),(i',j')}=1/k(Tk​)(i,j),(i′,j′)​=1/k when i=ji=ji=j and i′=j′i'=j'i′=j′, and zero otherwise. The norm is ∥X∥1=Re⁡Tr⁡X†X\|X\|_1=\operatorname{Re}\operatorname{Tr}\sqrt{X^\dagger X}∥X∥1​=ReTrX†X​. The map Φ\PhiΦ must be built by the following finite recursion: identities are allowed; either party may apply a finite local Kraus family with ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I, follow each outcome by another recursively allowed map, and sum the resulting branches. No positivity or normalization assumption is made on ρ\rhoρ, and empty A,BA,BA,B are permitted. The case m=0m=0m=0 is allowed; when r−δ<0r-\delta<0r−δ<0, the displayed rate inequality alone imposes no positive lower bound on mmm.

QAlgorithms.Papers.Independent.Intakethird100.P90.pptRelativeEntropy

For finite types A,BA,BA,B with decidable equality, and any complex matrix ρ\rhoρ indexed by A×BA\times BA×B, this defines the extended real number

inf⁡σ⪰0, Tr⁡σ=1Γ(σ)⪰0D(ρ∥σ),Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b).\inf_{\substack{\sigma\succeq0,\ \operatorname{Tr}\sigma=1\\ \Gamma(\sigma)\succeq0}}D(\rho\Vert\sigma), \qquad \Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}.σ⪰0, Trσ=1Γ(σ)⪰0​inf​D(ρ∥σ),Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​.

Here

D(ρ∥σ)={Re⁡Tr⁡ ⁣(f(ρ)−ρ ℓ(σ))log⁡2,ker⁡σ⊆ker⁡ρ,+∞,otherwise,D(\rho\Vert\sigma)= \begin{cases} \dfrac{\operatorname{Re}\operatorname{Tr}\!\left(f(\rho)-\rho\,\ell(\sigma)\right)}{\log 2}, &\ker\sigma\subseteq\ker\rho,\\[2mm] +\infty,&\text{otherwise}, \end{cases}D(ρ∥σ)=⎩⎨⎧​log2ReTr(f(ρ)−ρℓ(σ))​,+∞,​kerσ⊆kerρ,otherwise,​

where f(x)=xlog⁡xf(x)=x\log xf(x)=xlogx and ℓ(x)=log⁡x\ell(x)=\log xℓ(x)=logx are applied through real continuous functional calculus. The real logarithm is totalized, with log⁡0=0\log0=0log0=0 and log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣ for nonzero real xxx; this real functional calculus gives zero on non-Hermitian matrices. No positivity or normalization is required of ρ\rhoρ, and finite values of DDD are not truncated to be nonnegative. Empty AAA or BBB are permitted; then no trace-one σ\sigmaσ exists, so the infimum is +∞+\infty+∞.

QAlgorithms.Papers.Independent.Intakethird100.P90.preparationRate

For finite types A,BA,BA,B with decidable equality, a Boolean selector bbb, any complex matrix ρ\rhoρ on A×BA\times BA×B, and r∈Rr\in\mathbb Rr∈R, this predicate means that r≥0r\ge0r≥0 and, for every ϵ>0\epsilon>0ϵ>0 and δ>0\delta>0δ>0, there is N∈NN\in\mathbb NN∈N such that every integer n≥Nn\ge Nn≥N with n>0n>0n>0 admits an integer m≥0m\ge0m≥0 and a map Φ\PhiΦ from matrices on {0,…,2m−1}2\{0,\ldots,2^m-1\}^2{0,…,2m−1}2 to matrices on {0,…,∣An∣−1}×{0,…,∣Bn∣−1}\{0,\ldots,|A^n|-1\}\times\{0,\ldots,|B^n|-1\}{0,…,∣An∣−1}×{0,…,∣Bn∣−1} satisfying

m≤n(r+δ),∥Φ(T2m)−ρ[n]^∥1≤ϵ.m\le n(r+\delta), \qquad \left\|\Phi(T_{2^m})-\widehat{\rho^{[n]}}\right\|_1\le\epsilon.m≤n(r+δ),​Φ(T2m​)−ρ[n]​​1​≤ϵ.

Here (Tk)(i,j),(i′,j′)=1/k(T_k)_{(i,j),(i',j')}=1/k(Tk​)(i,j),(i′,j′)​=1/k when i=ji=ji=j and i′=j′i'=j'i′=j′, and is zero otherwise; ρ(a,b),(a′,b′)[n]=∏t=0n−1ρ(at,bt),(at′,bt′)\rho^{[n]}_{(a,b),(a',b')}=\prod_{t=0}^{n-1}\rho_{(a_t,b_t),(a'_t,b'_t)}ρ(a,b),(a′,b′)[n]​=∏t=0n−1​ρ(at​,bt​),(at′​,bt′​)​; and the hat reindexes the two tuple registers using the chosen bijections with their finite ordinal enumerations. The norm is ∥X∥1=Re⁡Tr⁡X†X\|X\|_1=\operatorname{Re}\operatorname{Tr}\sqrt{X^\dagger X}∥X∥1​=ReTrX†X​. When bbb is true, Φ\PhiΦ must be obtainable by finitely many recursive local operations: either an identity map, or a complete finite local Kraus family on either party followed, for each outcome, by another such map, with the branches summed. When bbb is false, both Φ\PhiΦ and X↦Γ(Φ(Γ(X)))X\mapsto\Gamma(\Phi(\Gamma(X)))X↦Γ(Φ(Γ(X))) must be channels, where Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, and a channel means a map X↦∑jKjXKj†X\mapsto\sum_jK_jXK_j^\daggerX↦∑j​Kj​XKj†​ with finitely many KjK_jKj​ satisfying ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. No positivity or normalization hypothesis is imposed on ρ\rhoρ; empty underlying types are included. The case m=0m=0m=0, whose input is the one-dimensional matrix T1=(1)T_1=(1)T1​=(1), is allowed.

QAlgorithms.Papers.Independent.Intakethird100.P90.regularizedPptRelativeEntropy

For finite types A,BA,BA,B with decidable equality, and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines

inf⁡k≥11kinf⁡σ⪰0, Tr⁡σ=1Γ(σ)⪰0D(ρ[k]∥σ)\inf_{k\ge1}\frac1k \inf_{\substack{\sigma\succeq0,\ \operatorname{Tr}\sigma=1\\\Gamma(\sigma)\succeq0}} D(\rho^{[k]}\Vert\sigma)k≥1inf​k1​σ⪰0, Trσ=1Γ(σ)⪰0​inf​D(ρ[k]∥σ)

in the extended reals. The matrix ρ[k]\rho^{[k]}ρ[k], indexed by Ak×BkA^k\times B^kAk×Bk, has entries

ρ(a,b),(a′,b′)[k]=∏i=0k−1ρ(ai,bi),(ai′,bi′),\rho^{[k]}_{(a,b),(a',b')}=\prod_{i=0}^{k-1}\rho_{(a_i,b_i),(a'_i,b'_i)},ρ(a,b),(a′,b′)[k]​=i=0∏k−1​ρ(ai​,bi​),(ai′​,bi′​)​,

and Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, with the same formula for tuple indices. The quantity D(X∥Y)D(X\Vert Y)D(X∥Y) is +∞+\infty+∞ unless ker⁡Y⊆ker⁡X\ker Y\subseteq\ker XkerY⊆kerX; otherwise it is Re⁡Tr⁡(f(X)−Xℓ(Y))/log⁡2\operatorname{Re}\operatorname{Tr}(f(X)-X\ell(Y))/\log2ReTr(f(X)−Xℓ(Y))/log2, using real continuous functional calculus with f(x)=xlog⁡xf(x)=x\log xf(x)=xlogx, ℓ(x)=log⁡x\ell(x)=\log xℓ(x)=logx, the totalized real logarithm, and zero functional-calculus values on non-Hermitian matrices. No density-matrix assumption is made on ρ\rhoρ. Only positive copy numbers occur, so there is no division by zero and no zero-copy term. If either underlying type is empty, the value is +∞+\infty+∞. The definition takes an infimum and does not assert convergence of the normalized quantities.

QAlgorithms.Papers.Independent.Intakethird100.P90.regularizedRains

For finite types A,BA,BA,B with decidable equality, and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the extended real number

inf⁡k≥11kinf⁡σ⪰0∥Γ(σ)∥1≤1D(ρ[k]∥σ).\inf_{k\ge1}\frac1k \inf_{\substack{\sigma\succeq0\\\|\Gamma(\sigma)\|_1\le1}} D(\rho^{[k]}\Vert\sigma).k≥1inf​k1​σ⪰0∥Γ(σ)∥1​≤1​inf​D(ρ[k]∥σ).

Here ρ(a,b),(a′,b′)[k]=∏i=0k−1ρ(ai,bi),(ai′,bi′)\rho^{[k]}_{(a,b),(a',b')}=\prod_{i=0}^{k-1}\rho_{(a_i,b_i),(a'_i,b'_i)}ρ(a,b),(a′,b′)[k]​=∏i=0k−1​ρ(ai​,bi​),(ai′​,bi′​)​ on Ak×BkA^k\times B^kAk×Bk, Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, and ∥X∥1=Re⁡Tr⁡X†X\|X\|_1=\operatorname{Re}\operatorname{Tr}\sqrt{X^\dagger X}∥X∥1​=ReTrX†X​. The quantity D(X∥Y)D(X\Vert Y)D(X∥Y) equals +∞+\infty+∞ when ker⁡Y⊈ker⁡X\ker Y\not\subseteq\ker XkerY⊆kerX, and otherwise equals Re⁡Tr⁡(f(X)−Xℓ(Y))/log⁡2\operatorname{Re}\operatorname{Tr}(f(X)-X\ell(Y))/\log2ReTr(f(X)−Xℓ(Y))/log2, where f(x)=xlog⁡xf(x)=x\log xf(x)=xlogx and ℓ(x)=log⁡x\ell(x)=\log xℓ(x)=logx use real continuous functional calculus, the totalized real logarithm, and zero functional-calculus values on non-Hermitian matrices. Neither ρ\rhoρ nor the feasible σ\sigmaσ is required to have trace one; ρ\rhoρ is not required to be positive. Only positive copy numbers occur, and no limit is asserted. Empty AAA or BBB are allowed and give value zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.rhoV

This definition gives the complex matrix on {0,1,2}2\{0,1,2\}^2{0,1,2}2

ρV=14[(e(0,1)−e(1,0))(e(0,1)−e(1,0))†+(e(0,2)−e(2,0))(e(0,2)−e(2,0))†],\rho_V= \frac14\left[ (e_{(0,1)}-e_{(1,0)})(e_{(0,1)}-e_{(1,0)})^\dagger+ (e_{(0,2)}-e_{(2,0)})(e_{(0,2)}-e_{(2,0)})^\dagger \right],ρV​=41​[(e(0,1)​−e(1,0)​)(e(0,1)​−e(1,0)​)†+(e(0,2)​−e(2,0)​)(e(0,2)​−e(2,0)​)†],

where e(a,b)e_{(a,b)}e(a,b)​ is the coordinate vector at (a,b)(a,b)(a,b), and a vector outer product vv†vv^\daggervv† has entry vivj‾v_i\overline{v_j}vi​vj​​. Equivalently, it is one half the sum of the outer products of the two vectors (e(0,j)−e(j,0))/2(e_{(0,j)}-e_{(j,0)})/\sqrt2(e(0,j)​−e(j,0)​)/2​ for j=1,2j=1,2j=1,2.

QAlgorithms.Papers.Independent.Intakethird100.P90.stateRank

For a finite type III equipped with decidable equality and an arbitrary complex matrix ρ\rhoρ indexed by III, this definition returns the complex dimension of the range of the linear map

v⟼ρv,(ρv)i=∑j∈Iρijvj.v\longmapsto \rho v,\qquad (\rho v)_i=\sum_{j\in I}\rho_{ij}v_j.v⟼ρv,(ρv)i​=j∈I∑​ρij​vj​.

It is therefore the matrix rank as a natural number, without any assumption that ρ\rhoρ is a density matrix. The empty type III is included and gives rank zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.supportDual

For finite types A,BA,BA,B equipped with decidable equality and an arbitrary complex matrix ρ\rhoρ indexed by A×BA\times BA×B, this definition is the real infimum

inf⁡{t∈R:∃ Y, Y=Y†, Y−ΓPρ⪰0, tI−ΓY⪰0},\inf\left\{ t\in\mathbb R: \exists\,Y,\ Y=Y^\dagger,\ Y-\Gamma P_\rho\succeq0,\ tI-\Gamma Y\succeq0 \right\},inf{t∈R:∃Y, Y=Y†, Y−ΓPρ​⪰0, tI−ΓY⪰0},

where YYY is a complex matrix on A×BA\times BA×B, (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​, and Pρ=cfc⁡(f,ρ)P_\rho=\operatorname{cfc}(f,\rho)Pρ​=cfc(f,ρ), with f(0)=0f(0)=0f(0)=0 and f(t)=1f(t)=1f(t)=1 for nonzero real ttt, using real continuous functional calculus. No state assumption on ρ\rhoρ or explicit restriction t≥0t\ge0t≥0 is imposed. If either index type is empty, every real ttt is feasible; the real infimum of this unbounded-below set is defined to be zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.supportDualFeasible

For finite types A,BA,BA,B equipped with decidable equality, arbitrary complex matrices ρ,Y\rho,Yρ,Y indexed by A×BA\times BA×B, and a real number ttt, this predicate means

Y=Y†,Y−ΓPρ⪰0,tI−ΓY⪰0,Y=Y^\dagger,\qquad Y-\Gamma P_\rho\succeq0,\qquad tI-\Gamma Y\succeq0,Y=Y†,Y−ΓPρ​⪰0,tI−ΓY⪰0,

where III is the identity matrix, (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​, and Pρ=cfc⁡(f,ρ)P_\rho=\operatorname{cfc}(f,\rho)Pρ​=cfc(f,ρ) for f(0)=0f(0)=0f(0)=0 and f(t)=1f(t)=1f(t)=1 for nonzero real ttt, using real continuous functional calculus. There is no explicit nonnegativity condition on ttt, nor any state assumption on ρ\rhoρ. Empty index types are allowed; then the predicate holds for every real ttt and the unique matrix YYY.

QAlgorithms.Papers.Independent.Intakethird100.P90.supportEta

For finite types A,BA,BA,B with decidable equality and any complex matrix ρ\rhoρ on A×BA\times BA×B, this defines the real infimum

inf⁡Y=Y†Y−Γ(Pρ)⪰0Y+Γ(Pρ)⪰0∥Γ(Y)∥op,\inf_{\substack{Y=Y^\dagger\\ Y-\Gamma(P_\rho)\succeq0\\ Y+\Gamma(P_\rho)\succeq0}} \|\Gamma(Y)\|_{\mathrm{op}},Y=Y†Y−Γ(Pρ​)⪰0Y+Γ(Pρ​)⪰0​inf​∥Γ(Y)∥op​,

where ∥⋅∥op\|\cdot\|_{\mathrm{op}}∥⋅∥op​ is the operator norm for Euclidean norms, Γ(X)(a,b),(a′,b′)=X(a,b′),(a′,b)\Gamma(X)_{(a,b),(a',b')}=X_{(a,b'),(a',b)}Γ(X)(a,b),(a′,b′)​=X(a,b′),(a′,b)​, and Pρ=f(ρ)P_\rho=f(\rho)Pρ​=f(ρ) uses real continuous functional calculus with f(0)=0f(0)=0f(0)=0 and f(x)=1f(x)=1f(x)=1 for x≠0x\ne0x=0. For Hermitian ρ\rhoρ, this selects its nonzero eigenspaces; for non-Hermitian ρ\rhoρ, it gives zero. No density-matrix assumption is made on ρ\rhoρ, and no minimizing YYY is specified by the definition. If ρ=0\rho=0ρ=0, or either underlying type is empty, the value is zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.supportOverlap

For finite types A,BA,BA,B equipped with decidable equality and an arbitrary complex matrix ρ\rhoρ indexed by A×BA\times BA×B, this definition is

sup⁡σ⪰0, Tr⁡σ=1Γσ⪰0Re⁡Tr⁡(Pρσ),\sup_{\substack{\sigma\succeq0,\ \operatorname{Tr}\sigma=1\\ \Gamma\sigma\succeq0}} \operatorname{Re}\operatorname{Tr}(P_\rho\sigma),σ⪰0, Trσ=1Γσ⪰0​sup​ReTr(Pρ​σ),

where σ\sigmaσ ranges over complex matrices on A×BA\times BA×B, (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​, and Pρ=cfc⁡(f,ρ)P_\rho=\operatorname{cfc}(f,\rho)Pρ​=cfc(f,ρ), with f(0)=0f(0)=0f(0)=0 and f(t)=1f(t)=1f(t)=1 for nonzero real ttt, using real continuous functional calculus. No density or Hermiticity assumption is made on ρ\rhoρ. The supremum is taken in R\mathbb RR; if either index type is empty, no trace-one matrix σ\sigmaσ exists, and the real supremum of the empty set is zero.

QAlgorithms.Papers.Independent.Intakethird100.P90.wangDuanBound

For finite types A,BA,BA,B equipped with decidable equality and an arbitrary complex matrix ρ\rhoρ indexed by A×BA\times BA×B, this definition is

inf⁡X=X†X−ρ⪰0log⁡2 ⁣(Re⁡Tr⁡[cfc⁡ ⁣(⋅,(ΓX)†(ΓX))]),\inf_{\substack{X=X^\dagger\\X-\rho\succeq0}} \log_2\!\left( \operatorname{Re}\operatorname{Tr} \left[ \operatorname{cfc}\!\left(\sqrt{\cdot},(\Gamma X)^\dagger(\Gamma X)\right) \right]\right),X=X†X−ρ⪰0​inf​log2​(ReTr[cfc(⋅​,(ΓX)†(ΓX))]),

where XXX ranges over complex matrices on A×BA\times BA×B, and (ΓZ)(a,b),(a′,b′)=Z(a,b′),(a′,b)(\Gamma Z)_{(a,b),(a',b')}=Z_{(a,b'),(a',b)}(ΓZ)(a,b),(a′,b′)​=Z(a,b′),(a′,b)​. Thus the logarithm is applied to the explicitly defined trace norm of ΓX\Gamma XΓX. No density or Hermiticity hypothesis is imposed on ρ\rhoρ, and XXX itself is required to be Hermitian but is not separately required to be positive semidefinite. The logarithm is total, with log⁡20=0\log_2 0=0log2​0=0; the infimum is real-valued, with value zero for an empty or unbounded-below set of objective values. Empty index types are permitted and give value zero.

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