Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Explicit assisted lower curves and convexified purification upper bound

Open
QAlgorithms.Papers.Independent.Intakethird100.P27.isotropic_unassisted_bounds

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

prove2quantum

For the memoryless maximally mixed qubit source with Bell purification Φ\PhiΦ and distortion observable I−ΦI-\PhiI−Φ, let Rq(D)R^q(D)Rq(D) be the operational unassisted quantum rate in qubits per copy. For 0≤D≤10\le D\le10≤D≤1, define the two explicit curves from Theorems 11 and 12 by

Lea(D)={1−12H(1−D,D/3,D/3,D/3),D≤3/4,0,D>3/4,L→(D)={h2(1/2+D(1−D)),D<1/2,0,D≥1/2.L_{\rm ea}(D)=\begin{cases}1-\tfrac12H(1-D,D/3,D/3,D/3),&D\le3/4,\\0,&D>3/4,\end{cases}\qquad L_{\to}(D)=\begin{cases}h_2(1/2+\sqrt{D(1-D)}),&D<1/2,\\0,&D\ge1/2.\end{cases}Lea​(D)={1−21​H(1−D,D/3,D/3,D/3),0,​D≤3/4,D>3/4,​L→​(D)={h2​(1/2+D(1−D)​),0,​D<1/2,D≥1/2.​

Here HHH is Shannon entropy and h2(p)=−plog⁡2p−(1−p)log⁡2(1−p)h_2(p)=-p\log_2p-(1-p)\log_2(1-p)h2​(p)=−plog2​p−(1−p)log2​(1−p), with zero logarithmic terms contributing zero. Put τd=(1−d)Φ+(d/3)(I−Φ)\tau_d=(1-d)\Phi+(d/3)(I-\Phi)τd​=(1−d)Φ+(d/3)(I−Φ) and c(d)=Ep(τd)c(d)=E_p(\tau_d)c(d)=Ep​(τd​) for 0≤d≤3/40\le d\le3/40≤d≤3/4, and c(d)=0c(d)=0c(d)=0 for 3/4<d≤13/4<d\le13/4<d≤1. Define

U(D)=inf⁡{∑i=1mpic(di):m≥1, pi≥0, ∑ipi=1, 0≤di≤1, ∑ipidi≤D}.U(D)=\inf\left\{\sum_{i=1}^m p_i c(d_i):m\ge1,\ p_i\ge0,\ \sum_i p_i=1,\ 0\le d_i\le1,\ \sum_i p_id_i\le D\right\}.U(D)=inf{i=1∑m​pi​c(di​):m≥1, pi​≥0, i∑​pi​=1, 0≤di​≤1, i∑​pi​di​≤D}.

Then

max⁡{Lea(D),L→(D)}≤Rq(D)≤U(D).\max\{L_{\rm ea}(D),L_{\to}(D)\}\le R^q(D)\le U(D).max{Lea​(D),L→​(D)}≤Rq(D)≤U(D).

EpE_pEp​ is optimized over all finite normalized purifications with exactly the prescribed marginal. The upper bound uses finite time-sharing mixtures, not an assumed numerical evaluation; no exact formula for the unassisted rate is claimed.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P27
Formal statement
noncomputable section

namespace QAlgorithms.Papers.Independent.Intakethird100.P27

open Filter

open scoped Matrix ComplexConjugate ComplexOrder ENNReal

attribute [local instance] Classical.propDecidable Classical.decEq

open QAlgorithms.Papers.Independent.Dependency50.P19
  (first marginalRB mutualInformation bipartiteBlock epr)

open QAlgorithms.Papers.Independent.Information100.P08 (formation)

open QAlgorithms.Papers.Independent.Gap48.P20 (redistributionMarginal)

open QAlgorithms.Papers.Independent.Bridge50.P37 (isotropicState normalizedEntangledProjection)

/-- PDF pp.9--10, unnumbered bounds and footnote 2: the explicit curves from
Theorems 11 and 12 lower-bound the unassisted operational rate, which is at most
the convexified exact E_p strategy cost. The Devetak--Berger curve has no factor
1/2 before binary entropy and vanishes for D >= 1/2. This is not an exact formula
for the unassisted rate. -/
theorem isotropic_unassisted_bounds
    (S : DistortionSource (Fin 2) (Fin 2))
    (hsource : pureDensity S.psi = epr 2) (hdelta : S.delta = 1 - epr 2)
    (D : ℝ) (hD0 : 0 ≤ D) (hD1 : D ≤ 1) :
    max (ENNReal.ofReal (if D ≤ 3 / 4 then 1 - isotropicEntropy D / 2 else 0))
      (ENNReal.ofReal (if D < 1 / 2 then
        binEntropy2 (1 / 2 + Real.sqrt (D * (1 - D))) else 0)) ≤
      rateDistortion S .unassisted D ∧
    rateDistortion S .unassisted D ≤ isotropicConvexUpper D := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P27

end
Source
Quantum rate distortion coding with auxiliary resources, arXiv:1212.5316v3, PDF p. 9, §V-A, unnumbered; p.10 and footnote 2
Read-back

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

Let SSS consist of a normalized vector ψ∈C2×2\psi\in\mathbb C^{2\times2}ψ∈C2×2 and a positive semidefinite distortion matrix Δ\DeltaΔ, and suppose ∣ψ⟩⟨ψ∣=Ω|\psi\rangle\langle\psi|=\Omega∣ψ⟩⟨ψ∣=Ω and Δ=I−Ω\Delta=I-\OmegaΔ=I−Ω, where Ω(i,j),(k,l)=12\Omega_{(i,j),(k,l)}=\tfrac12Ω(i,j),(k,l)​=21​ if i=ji=ji=j and k=lk=lk=l, and is zero otherwise. For 0≤D≤10\le D\le10≤D≤1, define R(D)∈[0,∞]R(D)\in[0,\infty]R(D)∈[0,∞] as the infimum of all real r≥0r\ge0r≥0 such that, for every ϵ>0\epsilon>0ϵ>0 and s>0s>0s>0, every sufficiently large integer n≥1n\ge1n≥1 admits channels encoding the nnn source registers into a positive dimensional quantum register of dimension qqq and decoding it into nnn output qubits, with log⁡2q≤n(r+s)\log_2q\le n(r+s)log2​q≤n(r+s), and with

1n∑i=1nRe⁡tr⁡[(I−Ω) ωAiBi]≤D+ϵ.\frac1n\sum_{i=1}^n\operatorname{Re}\operatorname{tr}\bigl[(I-\Omega)\,\omega_{A_iB_i}\bigr]\le D+\epsilon.n1​i=1∑n​Retr[(I−Ω)ωAi​Bi​​]≤D+ϵ.

Here ω\omegaω is obtained by applying the encoder and decoder to the second halves of Ω⊗n\Omega^{\otimes n}Ω⊗n, leaving the first halves untouched; ωAiBi\omega_{A_iB_i}ωAi​Bi​​ is its iii-th reference–output marginal. Channels 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 code has one dimensional classical and shared registers. Let h2(x)=binEntropy⁡(x)/ln⁡2h_2(x)=\operatorname{binEntropy}(x)/\ln2h2​(x)=binEntropy(x)/ln2, so on [0,1][0,1][0,1] it is binary entropy with zero endpoint terms, and set

Hiso(D)=−(1−D)log⁡2(1−D)−3D3log⁡2D3.H_{\mathrm{iso}}(D)=-(1-D)\log_2(1-D)-3\frac D3\log_2\frac D3.Hiso​(D)=−(1−D)log2​(1−D)−33D​log2​3D​.

For a bipartite matrix ρ\rhoρ, let EP(ρ)E_P(\rho)EP​(ρ) be the infimum in [0,∞][0,\infty][0,∞] of max⁡{H(Tr⁡BV∣v⟩⟨v∣),0}\max\{H(\operatorname{Tr}_{BV}|v\rangle\langle v|),0\}max{H(TrBV​∣v⟩⟨v∣),0} over all normalized vectors vvv on (A×U)×(B×V)(A\times U)\times(B\times V)(A×U)×(B×V), with finite positive dimensional U,VU,VU,V, whose marginal on A×BA\times BA×B equals ρ\rhoρ; here (∣v⟩⟨v∣)ij=vivj‾(|v\rangle\langle v|)_{ij}=v_i\overline{v_j}(∣v⟩⟨v∣)ij​=vi​vj​​ and H(X)=−Re⁡tr⁡[f(X)]H(X)=-\operatorname{Re}\operatorname{tr}[f(X)]H(X)=−Retr[f(X)], where f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x is applied by continuous functional calculus. Put

C(d)={EP ⁣((1−d)Ω+d3(I−Ω)),d≤34,0,d>34,C(d)= \begin{cases} E_P\!\left((1-d)\Omega+\dfrac d3(I-\Omega)\right),&d\le\tfrac34,\\ 0,&d>\tfrac34, \end{cases}C(d)=⎩⎨⎧​EP​((1−d)Ω+3d​(I−Ω)),0,​d≤43​,d>43​,​

and let U(D)U(D)U(D) be the infimum of ∑i=1mpiC(di)\sum_{i=1}^m p_iC(d_i)∑i=1m​pi​C(di​) over every positive integer mmm, real weights pi≥0p_i\ge0pi​≥0 with ∑ipi=1\sum_i p_i=1∑i​pi​=1, and di∈[0,1]d_i\in[0,1]di​∈[0,1] satisfying ∑ipidi≤D\sum_i p_id_i\le D∑i​pi​di​≤D. The assertion is

max⁡ ⁣{max⁡ ⁣(0,{1−12Hiso(D),D≤34,0,D>34,),max⁡ ⁣(0,{h2 ⁣(12+D(1−D)),D<12,0,D≥12)}≤R(D)≤U(D).\max\!\left\{ \max\!\left(0, \begin{cases}1-\tfrac12H_{\mathrm{iso}}(D),&D\le\tfrac34,\\0,&D>\tfrac34,\end{cases}\right), \max\!\left(0, \begin{cases}h_2\!\left(\tfrac12+\sqrt{D(1-D)}\right),&D<\tfrac12,\\0,&D\ge\tfrac12\end{cases}\right) \right\} \le R(D)\le U(D).max{max(0,{1−21​Hiso​(D),0,​D≤43​,D>43​,​),max(0,{h2​(21​+D(1−D)​),0,​D<21​,D≥21​​)}≤R(D)≤U(D).

The endpoints D=0,1D=0,1D=0,1 are included; the first lower-bound branch includes D=34D=\tfrac34D=43​, while the second excludes D=12D=\tfrac12D=21​. All costs, infima, and weighted sums use extended nonnegative arithmetic, including 0⋅(+∞)=00\cdot(+\infty)=00⋅(+∞)=0.

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