Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Finite-energy minimum-entanglement resource family

Open
QAlgorithms.Papers.Independent.Intakethird100.P89.finiteEnergyOptimalResources

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

prove2quantum

Let r≥0r\ge0r≥0, τ>0\tau>0τ>0 and D=τ+1−e2r∣τ−1∣>0D=\tau+1-e^{2r}|\tau-1|>0D=τ+1−e2r∣τ−1∣>0. For r>0r>0r>0 this denominator condition is equivalent to tanh⁡r<τ<1/tanh⁡r\tanh r<\tau<1/\tanh rtanhr<τ<1/tanhr. Define

bmin⁡=e2rτ+e−2r−∣τ−1∣D,a=b+e−2r(τ−1)τ,c=b−e−2rτ.b_{\min}=\frac{e^{2r}\tau+e^{-2r}-|\tau-1|}{D},\qquad a=\frac{b+e^{-2r}(\tau-1)}{\tau},\qquad c=\frac{b-e^{-2r}}{\sqrt\tau}.bmin​=De2rτ+e−2r−∣τ−1∣​,a=τb+e−2r(τ−1)​,c=τ​b−e−2r​.

For every b≥bmin⁡b\ge b_{\min}b≥bmin​, the symmetric covariance with blocks A=aIA=aIA=aI, B=bIB=bIB=bI, C=−cσzC=-c\sigma_zC=−cσz​ is physical and is the doubled covariance of a normalized centered Gaussian density ρ\rhoρ on the complete two-oscillator space. It satisfies

∥ρΓ∥1=e2r,n‾AB=a+b−24<∞.\|\rho^\Gamma\|_1=e^{2r},\qquad \overline n_{AB}=\frac{a+b-2}{4}<\infty.∥ρΓ∥1​=e2r,nAB​=4a+b−2​<∞.

The deterministic BK protocol at gain τ\sqrt\tauτ​ exactly simulates, on every input density, the phase-insensitive channel with noise e−2r(1+τ)e^{-2r}(1+\tau)e−2r(1+τ). Among this displayed family the per-mode mean photon number is minimized at b=bmin⁡b=b_{\min}b=bmin​. Any finite-energy Gaussian BK realization of that channel with logarithmic negativity 2s2s2s, s≥0s\ge0s≥0, has s≥rs\ge rs≥r.

Physicality, state existence and exact simulation are conclusions of the result. The energy minimum is within family (9), with no claim of global energy optimality over all Gaussian resources.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Intakethird100_P89
Formal statement
noncomputable section

open MeasureTheory Filter

open scoped BigOperators ENNReal Topology ComplexOrder

attribute [local instance] Classical.propDecidable Classical.decEq

namespace QAlgorithms.Papers.Independent.Intakethird100.P89

open QAlgorithms.Papers.Independent.Dependency50.P33
  (State IsDensity IsTraceClassChannel tensorApply traceSize)

open QAlgorithms.Papers.Independent.Dependency50.P39
  (coherentKet coherentCoefficient operatorOfMatrix husimiQ HasFockMatrix CoherentRealized)

open QAlgorithms.Papers.Independent.Frontier30.P14 (weylCharacteristic)

open QAlgorithms.Papers.Independent.Claimed50.P19
  (Cov FockIndex DensityKernel PhysicalCovariance CharacteristicSummable
   multimodeWeylTrace covarianceCharacteristicValue standardFormCovariance)

open QAlgorithms.Papers.Independent.Frontier30.P08 (QuadratureIndex)

open QAlgorithms.Papers.Independent.Gap48.P26 (PhaseInsensitive coherentKernel)

open QAlgorithms.Entanglement (partialTransposeFirst)

/-- PDF p.3, (9): existence, physicality, exact all-input simulation and
minimum entanglement are CONCLUSIONS for the constructed finite-energy family.
The b=b_min energy minimum is only within the displayed family. -/
theorem finiteEnergyOptimalResources (r tau : ℝ) (hr : 0 ≤ r) (htau : 0 < tau)
    (hden : 0 < familyDenominator r tau) :
    (0 < r → (0 < familyDenominator r tau ↔
      Real.tanh r < tau ∧ tau < (Real.tanh r)⁻¹)) ∧
    (∀ b : ℝ, familyMinB r tau ≤ b →
      let a := familyA r tau b
      let c := familyC r tau b
      let V := standardFormCovariance a b (-c) c
      PhysicalCovariance V ∧
      ∃ (rho : State (ℕ × ℕ)) (Phi : OscillatorChannel),
        CenteredResource rho V ∧ twoModeMeanPhotonEnergy rho < ⊤ ∧
        twoModeMeanPhotonEnergy rho = ENNReal.ofReal ((a + b - 2) / 4) ∧
        traceSize (partialTransposeFirst rho.coefficients) = ENNReal.ofReal (Real.exp (2 * r)) ∧
        IsBK rho (Real.sqrt tau) Phi ∧
        PhaseInsensitive Phi tau (Real.exp (-2 * r) * (1 + tau))) ∧
    (∀ b : ℝ, familyMinB r tau ≤ b →
      familyMinEnergy r tau ≤ familyEnergy r tau b) ∧
    FiniteBKSimulates r tau (Real.exp (-2 * r) * (1 + tau)) ∧
    (∀ s : ℝ, 0 ≤ s → FiniteBKSimulates s tau (Real.exp (-2 * r) * (1 + tau)) → r ≤ s) := by sorry

end QAlgorithms.Papers.Independent.Intakethird100.P89

end
Source
Optimal continuous variable quantum teleportation with limited resources, arXiv:1705.03017v5, PDF p. 3, Finite-energy minimum-entanglement Gaussian resource family, Eq. (9)
Read-back

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

For real r,τr,\taur,τ satisfying r≥0r\ge0r≥0, τ>0\tau>0τ>0, and

dr(τ)=τ+1−e2r∣τ−1∣>0,d_r(\tau)=\tau+1-e^{2r}|\tau-1|>0,dr​(τ)=τ+1−e2r∣τ−1∣>0,

define

bmin⁡=e2rτ+e−2r−∣τ−1∣dr(τ),a(b)=b+e−2r(τ−1)τ,c(b)=b−e−2rτ,e(b)=a(b)+b−24.b_{\min}=\frac{e^{2r}\tau+e^{-2r}-|\tau-1|}{d_r(\tau)},\qquad a(b)=\frac{b+e^{-2r}(\tau-1)}{\tau},\qquad c(b)=\frac{b-e^{-2r}}{\sqrt\tau},\qquad e(b)=\frac{a(b)+b-2}{4}.bmin​=dr​(τ)e2rτ+e−2r−∣τ−1∣​,a(b)=τb+e−2r(τ−1)​,c(b)=τ​b−e−2r​,e(b)=4a(b)+b−2​.

If r>0r>0r>0, then dr(τ)>0d_r(\tau)>0dr​(τ)>0 is equivalent to tanh⁡r<τ<(tanh⁡r)−1\tanh r<\tau<(\tanh r)^{-1}tanhr<τ<(tanhr)−1. For every real b≥bmin⁡b\ge b_{\min}b≥bmin​, set

Vb=(a(b)0−c(b)00a(b)0c(b)−c(b)0b00c(b)0b).V_b= \begin{pmatrix} a(b)&0&-c(b)&0\\ 0&a(b)&0&c(b)\\ -c(b)&0&b&0\\ 0&c(b)&0&b \end{pmatrix}.Vb​=​a(b)0−c(b)0​0a(b)0c(b)​−c(b)0b0​0c(b)0b​​.

This matrix is positive definite and Vb+idiag⁡(J,J)V_b+i\operatorname{diag}(J,J)Vb​+idiag(J,J) is positive semidefinite, where J=(01−10)J=\left(\begin{smallmatrix}0&1\\-1&0\end{smallmatrix}\right)J=(0−1​10​). The statement uses the following kernel notions: a trace-class kernel XXX on III has an expansion Xij=∑kuk(i)vk(j)‾X_{ij}=\sum_k u_k(i)\overline{v_k(j)}Xij​=∑k​uk​(i)vk​(j)​ with uk,vk∈ℓ2(I;C)u_k,v_k\in\ell^2(I;\mathbb C)uk​,vk​∈ℓ2(I;C) and ∑k∥uk∥∥vk∥<∞\sum_k\|u_k\|\|v_k\|<\infty∑k​∥uk​∥∥vk​∥<∞; a density kernel additionally has real nonnegative finite quadratic forms and a summable real diagonal of sum one; and N(X)N(X)N(X) is the infimum of these norm-product sums in [0,∞][0,\infty][0,∞], with empty infimum ∞\infty∞. Write

Dmn(z)=e−∣z∣2/2m!n!∑j=0min⁡(m,n)zm−j(−z‾)n−jj!(m−j)!(n−j)!,χX(z)=∑m∑nXmnDnm(z),D_{mn}(z)=e^{-|z|^2/2}\sqrt{m!n!} \sum_{j=0}^{\min(m,n)} \frac{z^{m-j}(-\overline z)^{n-j}}{j!(m-j)!(n-j)!}, \quad \chi_X(z)=\sum_m\sum_nX_{mn}D_{nm}(z),Dmn​(z)=e−∣z∣2/2m!n!​j=0∑min(m,n)​j!(m−j)!(n−j)!zm−j(−z)n−j​,χX​(z)=m∑​n∑​Xmn​Dnm​(z), χR(2)(v)=∑a∑bRab∏j=01Dbjaj ⁣(−vj,1+ivj,02),E(R)=∑p∈N2[p0+p12Re⁡Rpp]+,R(a,b),(c,d)Γ=R(c,b),(a,d).\chi_R^{(2)}(v)=\sum_a\sum_bR_{ab} \prod_{j=0}^1D_{b_ja_j}\!\left(\frac{-v_{j,1}+iv_{j,0}}{\sqrt2}\right), \quad E(R)=\sum_{p\in\mathbb N^2} \left[\frac{p_0+p_1}{2}\operatorname{Re}R_{pp}\right]_+, \quad R^\Gamma_{(a,b),(c,d)}=R_{(c,b),(a,d)}.χR(2)​(v)=a∑​b∑​Rab​j=0∏1​Dbj​aj​​(2​−vj,1​+ivj,0​​),E(R)=p∈N2∑​[2p0​+p1​​ReRpp​]+​,R(a,b),(c,d)Γ​=R(c,b),(a,d)​.

Here [x]+=max⁡(x,0)[x]_+=\max(x,0)[x]+​=max(x,0) is viewed in [0,∞][0,\infty][0,∞]. For every b≥bmin⁡b\ge b_{\min}b≥bmin​ there exist a density kernel RRR on N2\mathbb N^2N2 and a function Φ\PhiΦ on one-mode kernels such that all the inner and outer sums defining χR(2)(v)\chi_R^{(2)}(v)χR(2)​(v) are summable for every vvv, χR(2)(v)=e−vTVbv/4\chi_R^{(2)}(v)=e^{-v^TV_bv/4}χR(2)​(v)=e−vTVb​v/4, E(R)<∞E(R)<\inftyE(R)<∞, E(R)=[e(b)]+E(R)=[e(b)]_+E(R)=[e(b)]+​, and N(RΓ)=[e2r]+N(R^\Gamma)=[e^{2r}]_+N(RΓ)=[e2r]+​. The function Φ\PhiΦ preserves trace-class kernels, is additive and complex homogeneous on them, and preserves density kernels under every finite first-factor ampliation; the dimension-zero test is vacuous. For every one-mode density kernel SSS and every z∈Cz\in\mathbb Cz∈C, it satisfies

χΦ(S)(z)=χS(τz) χR(2)(2τIm⁡z,2τRe⁡z,2Im⁡z,−2Re⁡z)=χS(τz)e−e−2r(1+τ)∣z∣2/2,\chi_{\Phi(S)}(z)= \chi_S(\sqrt\tau z)\, \chi_R^{(2)}(\sqrt{2\tau}\operatorname{Im}z,\sqrt{2\tau}\operatorname{Re}z, \sqrt2\operatorname{Im}z,-\sqrt2\operatorname{Re}z) =\chi_S(\sqrt\tau z)e^{-e^{-2r}(1+\tau)|z|^2/2},χΦ(S)​(z)=χS​(τ​z)χR(2)​(2τ​Imz,2τ​Rez,2​Imz,−2​Rez)=χS​(τ​z)e−e−2r(1+τ)∣z∣2/2,

and ∣τ−1∣≤e−2r(1+τ)|\tau-1|\le e^{-2r}(1+\tau)∣τ−1∣≤e−2r(1+τ). Also e(bmin⁡)≤e(b)e(b_{\min})\le e(b)e(bmin​)≤e(b) for every b≥bmin⁡b\ge b_{\min}b≥bmin​. In particular, a realization with these finite-energy, centered covariance, channel, and gain conditions exists at noise e−2r(1+τ)e^{-2r}(1+\tau)e−2r(1+τ). Finally, for every s≥0s\ge0s≥0, if such a realization exists at the same τ\tauτ and noise, with N(RΓ)=[e2s]+N(R^\Gamma)=[e^{2s}]_+N(RΓ)=[e2s]+​ instead of [e2r]+[e^{2r}]_+[e2r]+​ and with any centered physical covariance, then r≤sr\le sr≤s. The assertions include r=0r=0r=0; only the interval equivalence is conditional on r>0r>0r>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