Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Attained one-copy deterministic SDP characterizations

Open
QAlgorithms.Papers.Independent.Intakethird100.P83.deterministicSDPs

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

prove2quantum

Let ρ\rhoρ be a normalized state on arbitrary nonzero finite-dimensional local spaces, let PρP_\rhoPρ​ be its orthogonal support projector, and let FΓF_\GammaFΓ​ be the positive-real-target SDP (3). Independently define E1(ρ)=sup⁡{log⁡2k:k>0, FΓ(ρ,k)=1}E_1(\rho)=\sup\{\log_2 k:k>0,\ F_\Gamma(\rho,k)=1\}E1​(ρ)=sup{log2​k:k>0, FΓ​(ρ,k)=1}. Define

W0(ρ)=max⁡R{Re⁡Tr⁡(ρR):0≤R≤(Re⁡Tr⁡(ρR))I, −I≤TB(R)≤I},m(ρ)=min⁡Pρ≤R≤I∥TB(R)∥∞.\begin{aligned}W_0(\rho)&=\max_R\{\operatorname{Re}\operatorname{Tr}(\rho R):0\le R\le (\operatorname{Re}\operatorname{Tr}(\rho R))I,\ -I\le T_B(R)\le I\},\\m(\rho)&=\min_{P_\rho\le R\le I}\|T_B(R)\|_\infty.\end{aligned}W0​(ρ)m(ρ)​=Rmax​{ReTr(ρR):0≤R≤(ReTr(ρR))I, −I≤TB​(R)≤I},=Pρ​≤R≤Imin​∥TB​(R)∥∞​.​

Both extrema are attained, W0(ρ)≥1W_0(\rho)\ge1W0​(ρ)≥1, m(ρ)>0m(\rho)>0m(ρ)>0, and

E1(ρ)=log⁡2W0(ρ)=−log⁡2m(ρ).E_1(\rho)=\log_2W_0(\rho)=-\log_2m(\rho).E1​(ρ)=log2​W0​(ρ)=−log2​m(ρ).

The defining supremum for E1E_1E1​ is also attained at some positive real kkk. This is equivalent to (15), maximizing the negative logarithm over the support interval; the norm minimization itself is the SDP.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P83
Formal statement
noncomputable section

open scoped Matrix ComplexConjugate ComplexOrder

open Filter

namespace QAlgorithms.Papers.Independent.Intakethird100.P83

open QAlgorithms.Papers.Independent.Support50.P11 (Loewner)

open QAlgorithms.Papers.Independent.Dependency50.P32 (partialTranspose)

open QAlgorithms.Papers.Independent.Information100.P34 (supportProjector)

variable {A B C D : Type*}

variable [Fintype A] [Fintype B] [Fintype C] [Fintype D]

variable [DecidableEq A] [DecidableEq B] [DecidableEq C] [DecidableEq D]

/-- Eqs. (14)–(15), PDF p.5, using the independent F_Gamma definition of
oneCopyExact on p.4. Both SDP optima and the positive-real target size are attained;
the norm is minimized before applying the negative logarithm. -/
theorem deterministicSDPs [Nonempty A] [Nonempty B]
    (rho : Matrix (A × B) (A × B) ℂ) (hrho : IsDensityMatrix rho) :
    oneCopyExact rho = Real.logb 2 (wZeroValue rho) ∧
    oneCopyExact rho = -Real.logb 2 (supportNormMinimum rho) ∧
    1 ≤ wZeroValue rho ∧ 0 < supportNormMinimum rho ∧
    (∃ R, WZeroFeasible rho R ∧ effectProb rho R = wZeroValue rho) ∧
    (∃ R, Loewner (supportProjector rho) R ∧ Loewner R 1 ∧
      complexSpectralNorm (partialTranspose R) = supportNormMinimum rho) ∧
    (∃ k : ℝ, 0 < k ∧ pptDistillationFidelity rho k = 1 ∧ oneCopyExact rho = Real.logb 2 k) := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P83

end
Source
Improved Semidefinite Programming Upper Bound on Distillable Entanglement, arXiv:1601.07940v4, PDF p. 5, Unnumbered main result: one-copy deterministic SDPs (14)–(15)
Read-back

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

For any finite nonempty types A,BA,BA,B, each equipped with decidable equality, and any complex matrix ρ\rhoρ indexed by A×BA\times BA×B that is positive semidefinite with trace 111, write X(a,b),(a′,b′)Γ=X(a,b′),(a′,b)X^\Gamma_{(a,b),(a',b')}=X_{(a,b'),(a',b)}X(a,b),(a′,b′)Γ​=X(a,b′),(a′,b)​, let PρP_\rhoPρ​ be the orthogonal projector onto the range of ρ\rhoρ, and let ∥⋅∥op\|\cdot\|_{\mathrm{op}}∥⋅∥op​ denote the Euclidean operator norm. Write X⪯YX\preceq YX⪯Y when Y−XY-XY−X is positive semidefinite. Define

Fρ(k)=sup⁡{Re⁡tr⁡(ρQ):k>0, 0⪯Q⪯I, −k−1I⪯QΓ⪯k−1I},D(ρ)=sup⁡{log⁡2k:k>0, Fρ(k)=1},F_\rho(k)=\sup\{\operatorname{Re}\operatorname{tr}(\rho Q):k>0,\ 0\preceq Q\preceq I,\ -k^{-1}I\preceq Q^\Gamma\preceq k^{-1}I\},\qquad D(\rho)=\sup\{\log_2 k:k>0,\ F_\rho(k)=1\},Fρ​(k)=sup{Retr(ρQ):k>0, 0⪯Q⪯I, −k−1I⪯QΓ⪯k−1I},D(ρ)=sup{log2​k:k>0, Fρ​(k)=1},

where kkk is real, not required to be an integer. Define w0(ρ)w_0(\rho)w0​(ρ) as the supremum of Re⁡tr⁡(ρR)\operatorname{Re}\operatorname{tr}(\rho R)Retr(ρR) over matrices RRR satisfying R⪰0R\succeq0R⪰0, R⪯Re⁡tr⁡(ρR) IR\preceq\operatorname{Re}\operatorname{tr}(\rho R)\,IR⪯Retr(ρR)I, and −I⪯RΓ⪯I-I\preceq R^\Gamma\preceq I−I⪯RΓ⪯I; define s(ρ)=inf⁡{∥RΓ∥op:Pρ⪯R⪯I}s(\rho)=\inf\{\|R^\Gamma\|_{\mathrm{op}}:P_\rho\preceq R\preceq I\}s(ρ)=inf{∥RΓ∥op​:Pρ​⪯R⪯I}. All these extrema are real-valued suprema or infima. Then

D(ρ)=log⁡2w0(ρ)=−log⁡2s(ρ),1≤w0(ρ),0<s(ρ).D(\rho)=\log_2 w_0(\rho)=-\log_2 s(\rho),\qquad 1\le w_0(\rho),\qquad 0<s(\rho).D(ρ)=log2​w0​(ρ)=−log2​s(ρ),1≤w0​(ρ),0<s(ρ).

There exists an RRR satisfying the displayed constraints defining w0(ρ)w_0(\rho)w0​(ρ) and attaining Re⁡tr⁡(ρR)=w0(ρ)\operatorname{Re}\operatorname{tr}(\rho R)=w_0(\rho)Retr(ρR)=w0​(ρ); there exists, possibly separately, an RRR with Pρ⪯R⪯IP_\rho\preceq R\preceq IPρ​⪯R⪯I attaining ∥RΓ∥op=s(ρ)\|R^\Gamma\|_{\mathrm{op}}=s(\rho)∥RΓ∥op​=s(ρ); and there exists a real k>0k>0k>0 such that Fρ(k)=1F_\rho(k)=1Fρ​(k)=1 and D(ρ)=log⁡2kD(\rho)=\log_2 kD(ρ)=log2​k.

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