Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for I-concurrence and tangle for isotropic states

Definition
QAlgorithms_Papers_Independent_Intakethird100_P21

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

prove2quantum

The objects the formalized statements of I-concurrence and tangle for isotropic states (arXiv:quant-ph/0208002) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • marginalTangle: PDF pp.3--5, (1.3),(2.3): the marginal linear-entropy cost.
  • pureTangle: PDF pp.3--4, (1.3),(2.2): pure tangle, with Alice the first factor. Only normalized vectors are used as physical states.
  • tangle: PDF p.3, (1.5): convex roof of pure tangle, not squared mixed concurrence. The attainment conclusion below pins the source's minimum convention.
  • IsPureProductState: PDF pp.2--4: a normalized pure product state, including arbitrary phase.
  • IsMaximallyEntangled: PDF p.4: maximal pure entanglement has uniform Schmidt probabilities on min(dA,dB) entries, also for unequal local dimensions.
  • IsOptimalPureEnsemble: PDF p.3, (1.4)--(1.5): an attaining normalized finite decomposition. A zero-weight member is unrestricted beyond being normalized.
  • isotropicReferenceVector: PDF p.5, (3.1): the accepted maximally entangled vector, with a total zero-dimensional fallback used nowhere in physical statements.
  • isotropicDensityAtFidelity: PDF p.5, (3.2): rho_(d,F); physical domain d>=2 and 0<=F<=1.
  • maximallyEntangledOverlap: PDF pp.5--6: squared overlap with Phi, not root fidelity.
  • SchmidtProbabilities: PDF pp.4,8: nonnegative, normalized squared Schmidt probabilities. Zero entries are explicitly permitted.
  • schmidtCost: PDF p.4, (2.2): pure tangle from squared Schmidt probabilities.
  • fidelityProfile: PDF p.6, (3.6): no conjugation is applied to the diagonal of V.
  • canonicalVector: PDF pp.4,6: canonical Schmidt vector before Alice's and Bob's unitaries.
  • profileValues: PDF pp.7--8, (3.15),(3.18)--(3.19): the constrained cost values. The nonempty physical domain of this profile is [1/d,1].
  • constrainedSchmidtTangle: PDF p.7, (3.15): raw squared-concurrence profile. Theorems always restrict its arguments to the nonempty feasible domain [1/d,1].
  • rawC: PDF p.7: raw concurrence profile c_d=sqrt(q_d).
  • minimizingLargeAmplitude: PDF p.9, (3.31): sqrt(F/d)(1+w sqrt(d-1)), written without w.
  • minimizingSmallAmplitude: PDF p.9, (3.32): sqrt(F/d)(1-w/sqrt(d-1)), written without w.
  • explicitQ: PDF p.9, (3.30): factor 2 multiplies the entire parenthesis.
  • minimizingProbabilities: PDF p.9, (3.30)--(3.32): one gamma-squared entry, d-1 delta-squared entries.
  • ScalarMixture: PDF p.7, (3.14): finite normalized scalar mixtures in a fixed interval.
  • hullValues: PDF p.7, (3.14): attainable average costs with mean fidelity F.
  • scalarConvexRoofAtFidelity: PDF p.7: scalar lower convex hull, only evaluated on [1/d,1].
  • IsConvexEnvelope: PDF p.7: the largest convex minorant, stated as the full universal property.
  • relaxedValues: PDF p.6, (3.9): relaxed pure ensembles constrain only average fidelity, not the full density. Via Schmidt decomposition these are exactly the source's (p_j,mu_j,V_j); component fidelities can be anywhere in [0,1].
  • relaxedRoof: PDF p.6, (3.9)--(3.10): fidelity-only minimum; attainment is stated below.
  • tangentPoint: PDF p.13: supporting tangent breakpoint, with d^2 in the denominator.
  • tangleFormula: PDF p.13, (3.41): the three branches of isotropic tangle (d>=3).
  • concurrenceFormula: PDF p.13, (3.42): only 2d/(d-1) is under the radical.
  • spectralTrace: PDF p.4, (2.3): real spectral trace for a Hermitian matrix. Finite spectra allow arbitrary scalar f; no operator-concavity is assumed.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_FiniteRegisters
import Definitions.Def_QAlgorithms_Defs_BipartiteSeparability
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P03
import Definitions.Def_QAlgorithms_Papers_Independent_Frontier30_P30
import Definitions.Def_QAlgorithms_Papers_Independent_Bridge50_P37

/-!
Statements from Rungta--Caves, arXiv:quant-ph/0208002v1.
All page references are PDF pages. Only pinned accepted modules are imported.
The tangle roof is distinct from the square of the mixed concurrence.
Source corrections on pp.6,7,10 are documented in MODERATION_NOTES.md;
independent moderation and mechanical compilation remain pending.
-/

noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P21

open scoped Matrix ComplexConjugate
open Filter
open QAlgorithms
open QAlgorithms.Papers.Independent.Information100.P03
  (PureEnsemble convexRoof pureConcurrence concurrence SchmidtData schmidtVector)

attribute [local instance] Classical.propDecidable Classical.decEq

/-- PDF pp.3--5, (1.3),(2.3): the marginal linear-entropy cost. -/
def marginalTangle {I : Type*} [Fintype I] (A : Matrix I I ℂ) : ℝ :=
  2 * (1 - (A * A).trace.re)

/-- PDF pp.3--4, (1.3),(2.2): pure tangle, with Alice the first factor.
Only normalized vectors are used as physical states. -/
def pureTangle {I J : Type*} [Fintype I] [Fintype J]
    (psi : EuclideanSpace ℂ (I × J)) : ℝ :=
  marginalTangle (traceRight (pureDensity psi))

/-- PDF p.3, (1.5): convex roof of pure tangle, not squared mixed concurrence.
The attainment conclusion below pins the source's minimum convention. -/
def tangle {I J : Type*} [Fintype I] [Fintype J]
    (rho : Matrix (I × J) (I × J) ℂ) : ℝ :=
  convexRoof pureTangle rho

/-- PDF pp.2--4: a normalized pure product state, including arbitrary phase. -/
def IsPureProductState {I J : Type*} [Fintype I] [Fintype J]
    (psi : EuclideanSpace ℂ (I × J)) : Prop :=
  ∃ a : EuclideanSpace ℂ I, ∃ b : EuclideanSpace ℂ J,
    IsState a ∧ IsState b ∧ psi = tensorVec a b

/-- PDF p.4: maximal pure entanglement has uniform Schmidt probabilities
on min(dA,dB) entries, also for unequal local dimensions. -/
def IsMaximallyEntangled {a b : ℕ} (psi : EuclideanSpace ℂ (Fin a × Fin b)) : Prop :=
  ∃ s : SchmidtData a b, psi = schmidtVector s ∧
    ∀ i, s.probability i = 1 / (min a b : ℝ)

/-- PDF p.3, (1.4)--(1.5): an attaining normalized finite decomposition.
A zero-weight member is unrestricted beyond being normalized. -/
def IsOptimalPureEnsemble {I : Type*} [Fintype I] (cost : EuclideanSpace ℂ I → ℝ)
    (rho : Matrix I I ℂ) (e : PureEnsemble I) : Prop :=
  ensembleDensity e.weight e.vector = rho ∧
    (∑ k, e.weight k * cost (e.vector k)) = convexRoof cost rho

/-- PDF p.5, (3.1): the accepted maximally entangled vector, with a total
zero-dimensional fallback used nowhere in physical statements. -/
def isotropicReferenceVector (d : ℕ) : EuclideanSpace ℂ (Fin d × Fin d) :=
  if hd : 0 < d then
    letI : Nonempty (Fin d) := ⟨⟨0, hd⟩⟩
    finiteMaxEntangled (Fin d)
  else 0

/-- PDF p.5, (3.2): rho_(d,F); physical domain d>=2 and 0<=F<=1. -/
def isotropicDensityAtFidelity (d : ℕ) (F : ℝ) : Matrix (Fin d × Fin d) (Fin d × Fin d) ℂ :=
  if hd : 0 < d then
    letI : NeZero d := ⟨Nat.ne_of_gt hd⟩
    Bridge50.P37.isotropicState d F
  else 0

/-- PDF pp.5--6: squared overlap with Phi, not root fidelity. -/
def maximallyEntangledOverlap {d : ℕ} (rho : Matrix (Fin d × Fin d) (Fin d × Fin d) ℂ) : ℝ :=
  effectProb (pureDensity (isotropicReferenceVector d)) rho

/-- PDF pp.4,8: nonnegative, normalized squared Schmidt probabilities.
Zero entries are explicitly permitted. -/
structure SchmidtProbabilities (d : ℕ) where
  probability : Fin d → ℝ
  nonnegative : ∀ i, 0 ≤ probability i
  total : ∑ i, probability i = 1

/-- PDF p.4, (2.2): pure tangle from squared Schmidt probabilities. -/
def schmidtCost {d : ℕ} (mu : Fin d → ℝ) : ℝ :=
  2 * (1 - ∑ i, mu i ^ 2)

/-- PDF p.6, (3.6): no conjugation is applied to the diagonal of V. -/
def fidelityProfile {d : ℕ} (mu : SchmidtProbabilities d)
    (V : Matrix (Fin d) (Fin d) ℂ) : ℝ :=
  Complex.normSq (∑ k, (Real.sqrt (mu.probability k) : ℂ) * V k k) / d

/-- PDF pp.4,6: canonical Schmidt vector before Alice's and Bob's unitaries. -/
def canonicalVector {d : ℕ} (mu : SchmidtProbabilities d) :
    EuclideanSpace ℂ (Fin d × Fin d) :=
  WithLp.toLp 2 fun p => if p.1 = p.2 then (Real.sqrt (mu.probability p.1) : ℂ) else 0

/-- PDF pp.7--8, (3.15),(3.18)--(3.19): the constrained cost values.
The nonempty physical domain of this profile is [1/d,1]. -/
def profileValues (d : ℕ) (F : ℝ) : Set ℝ :=
  {v | ∃ mu : SchmidtProbabilities d,
    (∑ k, Real.sqrt (mu.probability k)) = Real.sqrt ((d : ℝ) * F) ∧
    v = schmidtCost mu.probability}

/-- PDF p.7, (3.15): raw squared-concurrence profile. Theorems always
restrict its arguments to the nonempty feasible domain [1/d,1]. -/
def constrainedSchmidtTangle (d : ℕ) (F : ℝ) : ℝ := sInf (profileValues d F)

/-- PDF p.7: raw concurrence profile c_d=sqrt(q_d). -/
def rawC (d : ℕ) (F : ℝ) : ℝ := Real.sqrt (constrainedSchmidtTangle d F)

/-- PDF p.9, (3.31): sqrt(F/d)(1+w sqrt(d-1)), written without w. -/
def minimizingLargeAmplitude (d : ℕ) (F : ℝ) : ℝ :=
  Real.sqrt (F / d) + Real.sqrt (((d : ℝ) - 1) * (1 - F) / d)

/-- PDF p.9, (3.32): sqrt(F/d)(1-w/sqrt(d-1)), written without w. -/
def minimizingSmallAmplitude (d : ℕ) (F : ℝ) : ℝ :=
  Real.sqrt (F / d) - Real.sqrt ((1 - F) / ((d : ℝ) * ((d : ℝ) - 1)))

/-- PDF p.9, (3.30): factor 2 multiplies the entire parenthesis. -/
def explicitQ (d : ℕ) (F : ℝ) : ℝ :=
  2 * (1 - minimizingLargeAmplitude d F ^ 4 - ((d : ℝ) - 1) * minimizingSmallAmplitude d F ^ 4)

/-- PDF p.9, (3.30)--(3.32): one gamma-squared entry, d-1 delta-squared entries. -/
def minimizingProbabilities (d : ℕ) (F : ℝ) : Fin d → ℝ :=
  fun i => if i.val = 0 then minimizingLargeAmplitude d F ^ 2 else minimizingSmallAmplitude d F ^ 2

/-- PDF p.7, (3.14): finite normalized scalar mixtures in a fixed interval. -/
structure ScalarMixture (lo hi : ℝ) where
  count : ℕ
  weight : Fin count → ℝ
  point : Fin count → ℝ
  nonnegative : ∀ k, 0 ≤ weight k
  total : ∑ k, weight k = 1
  domain : ∀ k, point k ∈ Set.Icc lo hi

/-- PDF p.7, (3.14): attainable average costs with mean fidelity F. -/
def hullValues (d : ℕ) (f : ℝ → ℝ) (F : ℝ) : Set ℝ :=
  {v | ∃ m : ScalarMixture (1 / (d : ℝ)) 1,
    (∑ k, m.weight k * m.point k) = F ∧ v = ∑ k, m.weight k * f (m.point k)}

/-- PDF p.7: scalar lower convex hull, only evaluated on [1/d,1]. -/
def scalarConvexRoofAtFidelity (d : ℕ) (f : ℝ → ℝ) (F : ℝ) : ℝ := sInf (hullValues d f F)

/-- PDF p.7: the largest convex minorant, stated as the full universal property. -/
def IsConvexEnvelope (D : Set ℝ) (f g : ℝ → ℝ) : Prop :=
  ConvexOn ℝ D g ∧ (∀ x ∈ D, g x ≤ f x) ∧
    ∀ h : ℝ → ℝ, ConvexOn ℝ D h → (∀ x ∈ D, h x ≤ f x) → ∀ x ∈ D, h x ≤ g x

/-- PDF p.6, (3.9): relaxed pure ensembles constrain only average fidelity,
not the full density. Via Schmidt decomposition these are exactly the source's
(p_j,mu_j,V_j); component fidelities can be anywhere in [0,1]. -/
def relaxedValues (d : ℕ) (cost : EuclideanSpace ℂ (Fin d × Fin d) → ℝ)
    (F : ℝ) : Set ℝ :=
  {v | ∃ e : PureEnsemble (Fin d × Fin d),
    (∑ k, e.weight k * maximallyEntangledOverlap (pureDensity (e.vector k))) = F ∧
    v = ∑ k, e.weight k * cost (e.vector k)}

/-- PDF p.6, (3.9)--(3.10): fidelity-only minimum; attainment is stated below. -/
def relaxedRoof (d : ℕ) (cost : EuclideanSpace ℂ (Fin d × Fin d) → ℝ)
    (F : ℝ) : ℝ := sInf (relaxedValues d cost F)

/-- PDF p.13: supporting tangent breakpoint, with d^2 in the denominator. -/
def tangentPoint (d : ℕ) : ℝ := 4 * ((d : ℝ) - 1) / (d : ℝ)^2

/-- PDF p.13, (3.41): the three branches of isotropic tangle (d>=3). -/
def tangleFormula (d : ℕ) (F : ℝ) : ℝ :=
  if F ≤ 1 / (d : ℝ) then 0
  else if F ≤ tangentPoint d then explicitQ d F
  else 2 * d * (F - 1) / ((d : ℝ) - 1) + 2 * ((d : ℝ) - 1) / d

/-- PDF p.13, (3.42): only 2d/(d-1) is under the radical. -/
def concurrenceFormula (d : ℕ) (F : ℝ) : ℝ :=
  if F ≤ 1 / (d : ℝ) then 0
  else Real.sqrt (2 * d / ((d : ℝ) - 1)) * (F - 1 / d)

/-- PDF p.4, (2.3): real spectral trace for a Hermitian matrix.
Finite spectra allow arbitrary scalar f; no operator-concavity is assumed. -/
def spectralTrace {d : ℕ} (f : ℝ → ℝ) (rho : Matrix (Fin d) (Fin d) ℂ) : ℝ :=
  (cfc f rho).trace.re

/-! Sanity examples for the new paper-local substrate. -/

example : isotropicReferenceVector 1 (0, 0) = 1 := by
  simp [isotropicReferenceVector, finiteMaxEntangled]

example : tangentPoint 3 = (8 : ℝ) / 9 := by
  norm_num [tangentPoint]

example : concurrenceFormula 2 1 = 1 := by
  norm_num [concurrenceFormula]

example : tangleFormula 3 1 = (4 : ℝ) / 3 := by
  norm_num [tangleFormula, tangentPoint]

example : schmidtCost (fun _ : Fin 1 => (1 : ℝ)) = 0 := by
  simp [schmidtCost]

example {d : ℕ} (rho : Matrix (Fin d) (Fin d) ℂ) :
    spectralTrace (fun _ => 0) rho = 0 := by
  simp [spectralTrace]

example : pureTangle (ket (ι := Fin 1 × Fin 1) (0, 0)) = 0 := by
  simp [pureTangle, marginalTangle, traceRight, pureDensity, ket, Matrix.mul_apply, Matrix.trace]

example {I J : Type*} [Fintype I] [Fintype J]
    (a : EuclideanSpace ℂ I) (b : EuclideanSpace ℂ J) (ha : IsState a) (hb : IsState b) :
    IsPureProductState (tensorVec a b) := ⟨a, b, ha, hb, rfl⟩

example : (0 : ℝ) ∈ profileValues 1 1 := by
  refine ⟨⟨fun _ => 1, by intro k; norm_num, by simp⟩, ?_, ?_⟩
  · norm_num
  · simp [schmidtCost]

example (d : ℕ) : minimizingLargeAmplitude d 1 = minimizingSmallAmplitude d 1 := by
  simp [minimizingLargeAmplitude, minimizingSmallAmplitude]

example (d : ℕ) (f : ℝ → ℝ) (F : ℝ) (hF : F ∈ Set.Icc (1/(d : ℝ)) 1) :
    f F ∈ hullValues d f F := by
  let m : ScalarMixture (1/(d : ℝ)) 1 :=
    ⟨1, fun _ => 1, fun _ => F, by intro k; norm_num, by simp, fun _ => hF⟩
  refine ⟨m, ?_, ?_⟩ <;> simp [m]

end QAlgorithms.Papers.Independent.Intakethird100.P21
Source
I-concurrence and tangle for isotropic states, arXiv:quant-ph/0208002v1, https://arxiv.org/abs/quant-ph/0208002v1
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P21.IsConvexEnvelope

For any set D⊆RD\subseteq\mathbb RD⊆R and functions f,g:R→Rf,g:\mathbb R\to\mathbb Rf,g:R→R, this predicate requires that DDD be convex, that ggg be convex on DDD, that g(x)≤f(x)g(x)\le f(x)g(x)≤f(x) for every x∈Dx\in Dx∈D, and that every function h:R→Rh:\mathbb R\to\mathbb Rh:R→R which is convex on DDD and satisfies h(x)≤f(x)h(x)\le f(x)h(x)≤f(x) throughout DDD also satisfy h(x)≤g(x)h(x)\le g(x)h(x)≤g(x) throughout DDD. It imposes no restrictions on function values outside DDD. For D=∅D=\varnothingD=∅, the predicate holds for every fff and ggg.

QAlgorithms.Papers.Independent.Intakethird100.P21.IsMaximallyEntangled

For any natural numbers a,ba,ba,b and any vector ψ∈C{0,…,a−1}×{0,…,b−1}\psi\in\mathbb C^{\{0,\ldots,a-1\}\times\{0,\ldots,b-1\}}ψ∈C{0,…,a−1}×{0,…,b−1}, put m=min⁡(a,b)m=\min(a,b)m=min(a,b). This predicate means that there exist real numbers pi≥0p_i\ge0pi​≥0 for 0≤i<m0\le i<m0≤i<m, summing to one, and orthonormal families ui∈Cau_i\in\mathbb C^aui​∈Ca, vi∈Cbv_i\in\mathbb C^bvi​∈Cb, such that

ψ(x,y)=∑i=0m−1pi (ui)x(vi)y,pi=1m(0≤i<m).\psi_{(x,y)}=\sum_{i=0}^{m-1}\sqrt{p_i}\,(u_i)_x(v_i)_y, \qquad p_i=\frac1m\quad(0\le i<m).ψ(x,y)​=i=0∑m−1​pi​​(ui​)x​(vi​)y​,pi​=m1​(0≤i<m).

If m=0m=0m=0, the empty probability sum cannot equal one, so the predicate is false.

QAlgorithms.Papers.Independent.Intakethird100.P21.IsOptimalPureEnsemble

For any finite type III, any real-valued function ccc on CI\mathbb C^ICI, any complex matrix ρ\rhoρ indexed by III, and any supplied finite ensemble e=(K,w,ψ)e=(K,w,\psi)e=(K,w,ψ) whose weights are nonnegative and sum to one and whose vectors all have norm one, this predicate means

ρ=∑k=0K−1wkψkψk†,∑k=0K−1wkc(ψk)=inf⁡(L,p,ϕ): pℓ≥0, ∑ℓpℓ=1∥ϕℓ∥=1, ρ=∑ℓpℓϕℓϕℓ†∑ℓ=0L−1pℓc(ϕℓ),\rho=\sum_{k=0}^{K-1}w_k\psi_k\psi_k^\dagger, \qquad \sum_{k=0}^{K-1}w_k c(\psi_k) = \inf_{\substack{(L,p,\phi):\ p_\ell\ge0,\ \sum_\ell p_\ell=1\\ \|\phi_\ell\|=1,\ \rho=\sum_\ell p_\ell\phi_\ell\phi_\ell^\dagger}} \sum_{\ell=0}^{L-1}p_\ell c(\phi_\ell),ρ=k=0∑K−1​wk​ψk​ψk†​,k=0∑K−1​wk​c(ψk​)=(L,p,ϕ): pℓ​≥0, ∑ℓ​pℓ​=1∥ϕℓ​∥=1, ρ=∑ℓ​pℓ​ϕℓ​ϕℓ†​​inf​ℓ=0∑L−1​pℓ​c(ϕℓ​),

where each matrix ϕℓϕℓ†\phi_\ell\phi_\ell^\daggerϕℓ​ϕℓ†​ has entries (ϕℓ)i(ϕℓ)j‾(\phi_\ell)_i\overline{(\phi_\ell)_j}(ϕℓ​)i​(ϕℓ​)j​​. Thus the supplied ensemble must reproduce ρ\rhoρ and have cost equal to this real infimum. No boundedness assumption is imposed on ccc; if the set of competing costs is unbounded below, the real infimum convention used here gives 000. The supplied ensemble necessarily has K>0K>0K>0, and no such ensemble exists when III is empty.

QAlgorithms.Papers.Independent.Intakethird100.P21.IsPureProductState

For any finite types I,JI,JI,J and any vector ψ∈CI×J\psi\in\mathbb C^{I\times J}ψ∈CI×J, this predicate means that there exist vectors a∈CIa\in\mathbb C^Ia∈CI and b∈CJb\in\mathbb C^Jb∈CJ, each of norm one, such that

ψ(i,j)=aibjfor every i∈I, j∈J.\psi_{(i,j)}=a_i b_j \quad\text{for every }i\in I,\ j\in J.ψ(i,j)​=ai​bj​for every i∈I, j∈J.

If either type is empty, the predicate cannot hold because its corresponding vector space has no norm-one vector.

QAlgorithms.Papers.Independent.Intakethird100.P21.ScalarMixture

For any real endpoints lo\mathrm{lo}lo and hi\mathrm{hi}hi, an object of this structure consists of a natural number qqq, real weights wiw_iwi​, and real points xix_ixi​, indexed by i=0,…,q−1i=0,\ldots,q-1i=0,…,q−1, together with the requirements wi≥0w_i\ge0wi​≥0, ∑iwi=1\sum_iw_i=1∑i​wi​=1, and lo≤xi≤hi\mathrm{lo}\le x_i\le\mathrm{hi}lo≤xi​≤hi for every iii. Zero weights and repeated points are permitted, and the point restriction applies even at zero weights. The normalization excludes q=0q=0q=0; if lo>hi\mathrm{lo}>\mathrm{hi}lo>hi, no such object exists.

QAlgorithms.Papers.Independent.Intakethird100.P21.SchmidtProbabilities

For every natural number ddd, this structure consists of a function p:{0,…,d−1}→Rp:\{0,\ldots,d-1\}\to\mathbb Rp:{0,…,d−1}→R, together with the requirements

pi≥0for every i,∑i=0d−1pi=1.p_i\ge0\quad\text{for every }i, \qquad \sum_{i=0}^{d-1}p_i=1.pi​≥0for every i,i=0∑d−1​pi​=1.

There are no such structures for d=0d=0d=0, since the empty sum is zero.

QAlgorithms.Papers.Independent.Intakethird100.P21.canonicalVector

For every natural number ddd and every family μ=(μi)i=0d−1\mu=(\mu_i)_{i=0}^{d-1}μ=(μi​)i=0d−1​ of nonnegative real numbers with ∑iμi=1\sum_i\mu_i=1∑i​μi​=1, this defines the complex Euclidean vector indexed by pairs i,j∈{0,…,d−1}i,j\in\{0,\ldots,d-1\}i,j∈{0,…,d−1} whose (i,j)(i,j)(i,j)-coordinate is μi\sqrt{\mu_i}μi​​ when i=ji=ji=j, and 000 otherwise. For d=0d=0d=0, no such family μ\muμ exists, since its empty sum would have to equal 111.

QAlgorithms.Papers.Independent.Intakethird100.P21.concurrenceFormula

For every natural number ddd and every real FFF, this defines

{0,F≤1/d,2dd−1(F−1d),F>1/d.\begin{cases} 0,&F\le 1/d,\\[2pt] \displaystyle\sqrt{\frac{2d}{d-1}}\left(F-\frac1d\right),&F>1/d. \end{cases}⎩⎨⎧​0,d−12d​​(F−d1​),​F≤1/d,F>1/d.​

No positivity assumption on ddd and no interval restriction on FFF is imposed. Subtraction is in R\mathbb RR, division by zero is zero, and real square roots use their total-function convention. Consequently, for d=0d=0d=0 or d=1d=1d=1, the value is zero for every FFF.

QAlgorithms.Papers.Independent.Intakethird100.P21.constrainedSchmidtTangle

For every natural number ddd and real number FFF, this defines the real infimum of all numbers 2(1−∑i=0d−1μi2)2(1-\sum_{i=0}^{d-1}\mu_i^2)2(1−∑i=0d−1​μi2​) obtained from nonnegative real families (μi)i=0d−1(\mu_i)_{i=0}^{d-1}(μi​)i=0d−1​ satisfying ∑iμi=1\sum_i\mu_i=1∑i​μi​=1 and ∑iμi=dF\sum_i\sqrt{\mu_i}=\sqrt{dF}∑i​μi​​=dF​. No range assumption on FFF is required. If there is no such family—including when d=0d=0d=0, or when d>0d>0d>0 and F∉[1/d,1]F\notin[1/d,1]F∈/[1/d,1]—the real infimum used here is 000.

QAlgorithms.Papers.Independent.Intakethird100.P21.explicitQ

For every natural number ddd and real number FFF, this defines

2[1−(Fd+(d−1)(1−F)d)4−(d−1)(Fd−1−Fd(d−1))4].2\left[ 1-\left(\sqrt{\frac Fd}+\sqrt{\frac{(d-1)(1-F)}d}\right)^4 -(d-1)\left(\sqrt{\frac Fd}-\sqrt{\frac{1-F}{d(d-1)}}\right)^4 \right].2​1−(dF​​+d(d−1)(1−F)​​)4−(d−1)(dF​​−d(d−1)1−F​​)4​.

All occurrences of ddd in this expression are real-valued, division by zero is 000, and square roots of negative numbers are 000. No range or sign condition is imposed on the result; in particular, its value at d=0d=0d=0 is 222.

QAlgorithms.Papers.Independent.Intakethird100.P21.fidelityProfile

For every natural number ddd, every supplied probability vector pi≥0p_i\ge0pi​≥0 with ∑i=0d−1pi=1\sum_{i=0}^{d-1}p_i=1∑i=0d−1​pi​=1, and every complex d×dd\times dd×d matrix VVV, this defines

1d∣∑k=0d−1pk Vkk∣2.\frac1d\left|\sum_{k=0}^{d-1}\sqrt{p_k}\,V_{kk}\right|^2.d1​​k=0∑d−1​pk​​Vkk​​2.

The matrix VVV is not required to be unitary. Although d=0d=0d=0 is included in the natural-number parameter, no probability vector satisfying the required conditions exists in that case.

QAlgorithms.Papers.Independent.Intakethird100.P21.hullValues

For every natural number ddd, function f:R→Rf:\mathbb R\to\mathbb Rf:R→R, and real number FFF, this defines the set of real numbers vvv for which there exist a natural number qqq, nonnegative real weights wiw_iwi​ summing to 111, and points xi∈[1/d,1]x_i\in[1/d,1]xi​∈[1/d,1], indexed by i=0,…,q−1i=0,\ldots,q-1i=0,…,q−1, such that ∑iwixi=F\sum_iw_ix_i=F∑i​wi​xi​=F and v=∑iwif(xi)v=\sum_iw_if(x_i)v=∑i​wi​f(xi​). No regularity assumption is imposed on fff. The weight normalization excludes q=0q=0q=0. Since 1/0=01/0=01/0=0, the permitted interval is [0,1][0,1][0,1] when d=0d=0d=0; when d=1d=1d=1, it is the singleton {1}\{1\}{1}. The set is empty if FFF lies outside the permitted interval.

QAlgorithms.Papers.Independent.Intakethird100.P21.isotropicDensityAtFidelity

For every natural number ddd and every real number FFF, this defines a complex matrix indexed by {0,…,d−1}2\{0,\ldots,d-1\}^2{0,…,d−1}2. For d>0d>0d>0, let PdP_dPd​ have entry 1/d1/d1/d between (i,j)(i,j)(i,j) and (k,l)(k,l)(k,l) when i=ji=ji=j and k=lk=lk=l, and entry zero otherwise. The defined matrix is

FPd+1−Fd2−1(I−Pd).F P_d+\frac{1-F}{d^2-1}(I-P_d).FPd​+d2−11−F​(I−Pd​).

For d=0d=0d=0, it is the zero matrix. No restriction on FFF is imposed. Division is totalized with division by zero equal to zero, so for d=1d=1d=1 the second term vanishes and the matrix is FIF IFI.

QAlgorithms.Papers.Independent.Intakethird100.P21.isotropicReferenceVector

For every natural number ddd, this defines a vector Ωd∈C{0,…,d−1}2\Omega_d\in\mathbb C^{\{0,\ldots,d-1\}^2}Ωd​∈C{0,…,d−1}2. For d>0d>0d>0, its coordinates are

(Ωd)(i,j)={1/d,i=j,0,i≠j.(\Omega_d)_{(i,j)}= \begin{cases} 1/\sqrt d,&i=j,\\ 0,&i\ne j. \end{cases}(Ωd​)(i,j)​={1/d​,0,​i=j,i=j.​

For d=0d=0d=0, it is the zero vector in the space indexed by the empty set.

QAlgorithms.Papers.Independent.Intakethird100.P21.marginalTangle

For any finite type III and any complex matrix AAA indexed by I×II\times II×I, this defines the real number

2(1−Re⁡tr⁡(A2)).2\bigl(1-\operatorname{Re}\operatorname{tr}(A^2)\bigr).2(1−Retr(A2)).

No positivity, Hermitian, or trace normalization assumption is imposed on AAA. If III is empty, the trace is zero and the value is 222.

QAlgorithms.Papers.Independent.Intakethird100.P21.maximallyEntangledOverlap

For every natural number ddd and every complex matrix ρ\rhoρ indexed by {0,…,d−1}2\{0,\ldots,d-1\}^2{0,…,d−1}2, this defines Re⁡tr⁡(Pdρ)\operatorname{Re}\operatorname{tr}(P_d\rho)Retr(Pd​ρ), where, for d>0d>0d>0, PdP_dPd​ has entry 1/d1/d1/d between (i,i)(i,i)(i,i) and (j,j)(j,j)(j,j) and zero elsewhere. Equivalently, for d>0d>0d>0, the value is

1dRe⁡∑i,j=0d−1ρ(i,i),(j,j).\frac1d\operatorname{Re}\sum_{i,j=0}^{d-1}\rho_{(i,i),(j,j)}.d1​Rei,j=0∑d−1​ρ(i,i),(j,j)​.

For d=0d=0d=0, the value is zero. No positivity or normalization assumption is imposed on ρ\rhoρ.

QAlgorithms.Papers.Independent.Intakethird100.P21.minimizingLargeAmplitude

For every natural number ddd and real number FFF, this defines

Fd+(d−1)(1−F)d,\sqrt{\frac{F}{d}}+\sqrt{\frac{(d-1)(1-F)}{d}},dF​​+d(d−1)(1−F)​​,

where ddd is interpreted as a real number, so d−1d-1d−1 is real subtraction. There are no positivity or range assumptions on ddd or FFF; division by zero is defined to be 000, and the real square root of a negative number is 000. Thus the value is 000 when d=0d=0d=0, and F\sqrt FF​ when d=1d=1d=1.

QAlgorithms.Papers.Independent.Intakethird100.P21.minimizingProbabilities

For every natural number ddd and real number FFF, this defines a real-valued function on i∈{0,…,d−1}i\in\{0,\ldots,d-1\}i∈{0,…,d−1} by

i⟼{(Fd+(d−1)(1−F)d)2,i=0,(Fd−1−Fd(d−1))2,i≠0.i\longmapsto \begin{cases} \left(\sqrt{\frac Fd}+\sqrt{\frac{(d-1)(1-F)}d}\right)^2,&i=0,\\ \left(\sqrt{\frac Fd}-\sqrt{\frac{1-F}{d(d-1)}}\right)^2,&i\ne0. \end{cases}i⟼⎩⎨⎧​(dF​​+d(d−1)(1−F)​​)2,(dF​​−d(d−1)1−F​​)2,​i=0,i=0.​

Here subtraction involving ddd is real subtraction, division by zero is 000, and square roots of negative numbers are 000. No normalization condition is included in the definition. For d=0d=0d=0 this is the function on the empty index set; for d=1d=1d=1, its sole value is (F)2(\sqrt F)^2(F​)2.

QAlgorithms.Papers.Independent.Intakethird100.P21.minimizingSmallAmplitude

For every natural number ddd and real number FFF, this defines

Fd−1−Fd(d−1),\sqrt{\frac{F}{d}}-\sqrt{\frac{1-F}{d(d-1)}},dF​​−d(d−1)1−F​​,

where ddd is interpreted as a real number and d−1d-1d−1 is real subtraction. There are no positivity or range assumptions, division by zero is 000, and real square roots of negative numbers are 000. The value is therefore 000 when d=0d=0d=0, and F\sqrt FF​ when d=1d=1d=1; no nonnegativity of the resulting difference is required.

QAlgorithms.Papers.Independent.Intakethird100.P21.profileValues

For every natural number ddd and real number FFF, this defines the set

{ 2(1−∑i=0d−1μi2)  |  μi≥0,∑i=0d−1μi=1,∑i=0d−1μi=dF }.\left\{\,2\left(1-\sum_{i=0}^{d-1}\mu_i^2\right)\;\middle|\; \mu_i\ge0,\quad \sum_{i=0}^{d-1}\mu_i=1,\quad \sum_{i=0}^{d-1}\sqrt{\mu_i}=\sqrt{dF}\,\right\}.{2(1−i=0∑d−1​μi2​)​μi​≥0,i=0∑d−1​μi​=1,i=0∑d−1​μi​​=dF​}.

The parameter FFF has no imposed range, and the real square root is 000 on negative arguments. The set is empty when d=0d=0d=0; for d>0d>0d>0, it is empty unless 1/d≤F≤11/d\le F\le11/d≤F≤1.

QAlgorithms.Papers.Independent.Intakethird100.P21.pureTangle

For any finite types I,JI,JI,J and any vector ψ∈CI×J\psi\in\mathbb C^{I\times J}ψ∈CI×J, define the matrix RRR indexed by I×II\times II×I by

Rij=∑h∈Jψ(i,h)ψ(j,h)‾.R_{ij}=\sum_{h\in J}\psi_{(i,h)}\overline{\psi_{(j,h)}}.Rij​=h∈J∑​ψ(i,h)​ψ(j,h)​​.

This definition assigns to ψ\psiψ the real number 2(1−Re⁡tr⁡(R2))2\bigl(1-\operatorname{Re}\operatorname{tr}(R^2)\bigr)2(1−Retr(R2)). The vector ψ\psiψ need not have norm one. If either index type is empty, RRR is zero and the value is 222.

QAlgorithms.Papers.Independent.Intakethird100.P21.rawC

For every natural number ddd and real number FFF, this defines the real square root of the real infimum of

{ 2(1−∑i=0d−1μi2)  |  μi≥0,∑iμi=1,∑iμi=dF }.\left\{\,2\left(1-\sum_{i=0}^{d-1}\mu_i^2\right)\;\middle|\; \mu_i\ge0,\quad \sum_i\mu_i=1,\quad \sum_i\sqrt{\mu_i}=\sqrt{dF}\,\right\}.{2(1−i=0∑d−1​μi2​)​μi​≥0,i∑​μi​=1,i∑​μi​​=dF​}.

The real infimum of an empty set is taken to be 000, and the real square root of a negative number is also 000. In particular, this value is 000 when d=0d=0d=0 or when d>0d>0d>0 and F∉[1/d,1]F\notin[1/d,1]F∈/[1/d,1].

QAlgorithms.Papers.Independent.Intakethird100.P21.relaxedRoof

For every natural number ddd, every real-valued cost function ccc on Cd⊗Cd\mathbb C^{d}\otimes\mathbb C^{d}Cd⊗Cd, and every real FFF, this defines the real infimum of all averages ∑kwkc(ψk)\sum_kw_kc(\psi_k)∑k​wk​c(ψk​) over finite families of unit vectors ψk\psi_kψk​ with nonnegative weights summing to 111 and satisfying ∑kwkRe⁡tr⁡(∣Φd⟩⟨Φd∣ ∣ψk⟩⟨ψk∣)=F\sum_kw_k\operatorname{Re}\operatorname{tr}(|\Phi_d\rangle\langle\Phi_d|\,|\psi_k\rangle\langle\psi_k|)=F∑k​wk​Retr(∣Φd​⟩⟨Φd​∣∣ψk​⟩⟨ψk​∣)=F, where Φd=d−1/2∑i<d∣i,i⟩\Phi_d=d^{-1/2}\sum_{i<d}|i,i\rangleΦd​=d−1/2∑i<d​∣i,i⟩ for d>0d>0d>0, and Φ0=0\Phi_0=0Φ0​=0. There is no constraint on the ensemble’s average density matrix, and no regularity or boundedness assumption on ccc. This is a real-valued infimum: the convention gives 000 if the set is empty or unbounded below. In particular, it is 000 for d=0d=0d=0, for F∉[0,1]F\notin[0,1]F∈/[0,1], and for d=1d=1d=1 with F≠1F\ne1F=1.

QAlgorithms.Papers.Independent.Intakethird100.P21.relaxedValues

For every natural number ddd, every function c:Cd⊗Cd→Rc:\mathbb C^{d}\otimes\mathbb C^{d}\to\mathbb Rc:Cd⊗Cd→R on the corresponding Euclidean vector space, and every real number FFF, this defines the set of real numbers vvv for which there exist a natural number qqq, weights w0,…,wq−1≥0w_0,\ldots,w_{q-1}\ge0w0​,…,wq−1​≥0 summing to 111, and unit vectors ψ0,…,ψq−1∈Cd⊗Cd\psi_0,\ldots,\psi_{q-1}\in\mathbb C^{d}\otimes\mathbb C^{d}ψ0​,…,ψq−1​∈Cd⊗Cd, such that

∑k=0q−1wk⟨Φd,∣ψk⟩⟨ψk∣Φd⟩=F,v=∑k=0q−1wkc(ψk),\sum_{k=0}^{q-1}w_k\langle\Phi_d,|\psi_k\rangle\langle\psi_k|\Phi_d\rangle=F, \qquad v=\sum_{k=0}^{q-1}w_kc(\psi_k),k=0∑q−1​wk​⟨Φd​,∣ψk​⟩⟨ψk​∣Φd​⟩=F,v=k=0∑q−1​wk​c(ψk​),

where, for d>0d>0d>0, Φd=d−1/2∑i=0d−1∣i,i⟩\Phi_d=d^{-1/2}\sum_{i=0}^{d-1}|i,i\rangleΦd​=d−1/2∑i=0d−1​∣i,i⟩, and for d=0d=0d=0 the reference vector is zero. Here ∣ψ⟩⟨ψ∣|\psi\rangle\langle\psi|∣ψ⟩⟨ψ∣ has entries ψiψj‾\psi_i\overline{\psi_j}ψi​ψj​​, and the overlap means the real part of the trace of its product with ∣Φd⟩⟨Φd∣|\Phi_d\rangle\langle\Phi_d|∣Φd​⟩⟨Φd​∣. No particular average density matrix is prescribed. Zero weights are allowed, but q=0q=0q=0 cannot satisfy the weight sum. For d=0d=0d=0 the set is empty because there are no unit vectors; for d>0d>0d>0, it is empty when F∉[0,1]F\notin[0,1]F∈/[0,1], and for d=1d=1d=1 it is empty unless F=1F=1F=1.

QAlgorithms.Papers.Independent.Intakethird100.P21.scalarConvexRoofAtFidelity

For every natural number ddd, function f:R→Rf:\mathbb R\to\mathbb Rf:R→R, and real number FFF, this defines the real infimum of all finite weighted sums ∑iwif(xi)\sum_iw_if(x_i)∑i​wi​f(xi​) with wi≥0w_i\ge0wi​≥0, ∑iwi=1\sum_iw_i=1∑i​wi​=1, xi∈[1/d,1]x_i\in[1/d,1]xi​∈[1/d,1], and ∑iwixi=F\sum_iw_ix_i=F∑i​wi​xi​=F. No continuity, convexity, or boundedness assumption is imposed on fff. Division by zero gives 1/0=01/0=01/0=0, so the permitted interval for d=0d=0d=0 is [0,1][0,1][0,1]. If the set of weighted sums is empty or unbounded below, the real infimum used here is 000.

QAlgorithms.Papers.Independent.Intakethird100.P21.schmidtCost

For every natural number ddd and every function μ:{0,…,d−1}→R\mu:\{0,\ldots,d-1\}\to\mathbb Rμ:{0,…,d−1}→R, this defines

2(1−∑i=0d−1μi2).2\left(1-\sum_{i=0}^{d-1}\mu_i^2\right).2(1−i=0∑d−1​μi2​).

No nonnegativity or sum-to-one condition is imposed on μ\muμ. For d=0d=0d=0, the value is 222.

QAlgorithms.Papers.Independent.Intakethird100.P21.spectralTrace

For every natural number ddd, every function f:R→Rf:\mathbb R\to\mathbb Rf:R→R, and every complex d×dd\times dd×d matrix ρ\rhoρ, this defines Re⁡tr⁡(f(ρ))\operatorname{Re}\operatorname{tr}(f(\rho))Retr(f(ρ)), where f(ρ)f(\rho)f(ρ) is obtained using the real continuous functional calculus. For Hermitian ρ\rhoρ, this is ∑i=0d−1f(λi)\sum_{i=0}^{d-1}f(\lambda_i)∑i=0d−1​f(λi​), with eigenvalues counted with multiplicity. The definition assumes neither that ρ\rhoρ is Hermitian, positive semidefinite, or trace one, nor that fff is globally continuous; the functional calculus’s total-function conventions apply outside its applicable hypotheses. The value is zero when d=0d=0d=0.

QAlgorithms.Papers.Independent.Intakethird100.P21.tangentPoint

For every natural number ddd, this defines the real number

4(d−1)d2,\frac{4(d-1)}{d^2},d24(d−1)​,

with all arithmetic performed in R\mathbb RR, so d−1d-1d−1 is real subtraction. The definition includes d=0d=0d=0 and d=1d=1d=1; its value is 000 in both cases, using the convention that division by zero is zero.

QAlgorithms.Papers.Independent.Intakethird100.P21.tangle

For any finite types I,JI,JI,J and any complex matrix ρ\rhoρ indexed by I×JI\times JI×J, this defines the infimum of

∑k=0K−1wk 2(1−Re⁡tr⁡(Rk2)),(Rk)ij=∑h∈J(ψk)(i,h)(ψk)(j,h)‾,\sum_{k=0}^{K-1}w_k\,2\bigl(1-\operatorname{Re}\operatorname{tr}(R_k^2)\bigr), \qquad (R_k)_{ij}=\sum_{h\in J}(\psi_k)_{(i,h)}\overline{(\psi_k)_{(j,h)}},k=0∑K−1​wk​2(1−Retr(Rk2​)),(Rk​)ij​=h∈J∑​(ψk​)(i,h)​(ψk​)(j,h)​​,

over all natural numbers KKK, real weights wk≥0w_k\ge0wk​≥0 summing to one, and vectors ψk∈CI×J\psi_k\in\mathbb C^{I\times J}ψk​∈CI×J with ∥ψk∥=1\|\psi_k\|=1∥ψk​∥=1, subject to ρ=∑kwkψkψk†\rho=\sum_k w_k\psi_k\psi_k^\daggerρ=∑k​wk​ψk​ψk†​, where (ψkψk†)xy=(ψk)x(ψk)y‾(\psi_k\psi_k^\dagger)_{xy}=(\psi_k)_x\overline{(\psi_k)_y}(ψk​ψk†​)xy​=(ψk​)x​(ψk​)y​​. Every listed vector must be normalized, including those with zero weight; K=0K=0K=0 cannot satisfy the weight condition. No density matrix assumption is imposed on ρ\rhoρ. If there is no such ensemble, the real infimum of the empty set used here is 000; in particular, this occurs if either index type is empty.

QAlgorithms.Papers.Independent.Intakethird100.P21.tangleFormula

For every natural number ddd and every real FFF, define

u=Fd+(d−1)(1−F)d,v=Fd−1−Fd(d−1).u=\sqrt{\frac Fd}+\sqrt{\frac{(d-1)(1-F)}d}, \qquad v=\sqrt{\frac Fd}-\sqrt{\frac{1-F}{d(d-1)}}.u=dF​​+d(d−1)(1−F)​​,v=dF​​−d(d−1)1−F​​.

The declaration defines the real-valued function

{0,F≤1/d,2(1−u4−(d−1)v4),F>1/d and F≤4(d−1)/d2,2d(F−1)d−1+2(d−1)d,F>1/d and F>4(d−1)/d2.\begin{cases} 0,&F\le 1/d,\\[2pt] 2\bigl(1-u^4-(d-1)v^4\bigr), &F>1/d\ \text{and}\ F\le 4(d-1)/d^2,\\[2pt] \dfrac{2d(F-1)}{d-1}+\dfrac{2(d-1)}d, &F>1/d\ \text{and}\ F>4(d-1)/d^2. \end{cases}⎩⎨⎧​0,2(1−u4−(d−1)v4),d−12d(F−1)​+d2(d−1)​,​F≤1/d,F>1/d and F≤4(d−1)/d2,F>1/d and F>4(d−1)/d2.​

There is no restriction that F∈[0,1]F\in[0,1]F∈[0,1] or that ddd be positive. All subtractions here are real, division by zero is zero, and the real square root of a negative argument is zero. In particular, for d=0d=0d=0 or d=1d=1d=1, the resulting function is zero for every FFF.

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