Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

IV.C.1–2: calibrated statistical identities, Eqs.(44)–(48)

Proved
QAlgorithms.Papers.Independent.Gap48.P25.calibratedStatistics

by Goku · Oct 21, 2026 · Mathlib 0df444a (Lean v4.33.1)

prove2quantum

Let pAp_ApA​ be a normalized photon law, ν~S>0\widetilde\nu_S>0νS​>0, calibrated apparatus transmission b=tBη∈[0,1]b=t_B\eta\in[0,1]b=tB​η∈[0,1], 0≤pd<1/20\le p_d<1/20≤pd​<1/2, and photon forwarding and line-error probabilities fn,εn∈[0,1]f_n,\varepsilon_n\in[0,1]fn​,εn​∈[0,1]. Assume f0=0f_0=0f0​=0, excluding photon creation on vacuum inputs as in footnote 57. In the restricted one-forwarded-photon model define Rn,p=ν~SpA(n)fnbR_{n,p}=\widetilde\nu_Sp_A(n)f_nbRn,p​=νS​pA​(n)fn​b, Rn,d=ν~SpA(n)(1−fnb)2pdR_{n,d}=\widetilde\nu_Sp_A(n)(1-f_nb)2p_dRn,d​=νS​pA​(n)(1−fn​b)2pd​ and R=∑n(Rn,p+Rn,d)>0R=\sum_n(R_{n,p}+R_{n,d})>0R=∑n​(Rn,p​+Rn,d​)>0. For Y=∑nRn,p/RY=\sum_nR_{n,p}/RY=∑n​Rn,p​/R and δ=(1−Y)/2\delta=(1-Y)/2δ=(1−Y)/2,

0≤Y≤1,Y=1−2pdν~S/R1−2pd,∑nRn,dR=2δ.0\le Y\le1,\qquad Y={1-2p_d\widetilde\nu_S/R\over1-2p_d},\qquad {\sum_nR_{n,d}\over R}=2\delta.0≤Y≤1,Y=1−2pd​1−2pd​νS​/R​,R∑n​Rn,d​​=2δ.

The observed error is Q=(∑nRn,pεn)/R+δQ=(\sum_nR_{n,p}\varepsilon_n)/R+\deltaQ=(∑n​Rn,p​εn​)/R+δ, or Q=Yε+δQ=Y\varepsilon+\deltaQ=Yε+δ if Y>0Y>0Y>0 and ε\varepsilonε is the normalized photon-error average. These are calibrated statistical identities, without an unconditional security conclusion.

Preamble
import Definitions.Def_QAlgorithms_Papers_Independent_Gap48_P25
Formal statement
noncomputable section

namespace QAlgorithms.Papers.Independent.Gap48.P25

open scoped BigOperators ComplexConjugate

open Filter

open QAlgorithms.Papers.Independent.Claimed50.P08 (classicalEntropy)

/-- Eqs.(44)–(48), PDF p.31: calibrated statistical ratio, with IV.C's distinct
Y_n convention and independent dark-count law. Footnote 57 excludes photon
creation on vacuum inputs, so f 0 = 0; no operational security conclusion. -/
theorem calibratedStatistics (p : PhotonLaw) (nu b pd : ℝ) (f e : ℕ → ℝ)
    (hnu : 0 < nu) (hb : b ∈ Set.Icc (0 : ℝ) 1) (hpd : 0 ≤ pd ∧ pd < 1/2)
    (hf : ∀ n, f n ∈ Set.Icc (0 : ℝ) 1) (hf0 : f 0 = 0)
    (he : ∀ n, e n ∈ Set.Icc (0 : ℝ) 1)
    (hR : 0 < calibratedTotal p nu b pd f) :
    let R := calibratedTotal p nu b pd f
    let Y := calibratedFraction p nu b pd f
    let delta := (1-Y)/2
    let photonWrong := ∑' n, calibratedPhoton p nu b f n*e n
    let dark := ∑' n, calibratedDark p nu b pd f n
    Y ∈ Set.Icc (0 : ℝ) 1 ∧ Y = (1-2*pd*nu/R)/(1-2*pd) ∧
    dark/R = 2*delta ∧ (photonWrong+dark/2)/R = (photonWrong/R)+delta ∧
    (0 < Y → (photonWrong+dark/2)/R = Y*(photonWrong/(R*Y))+delta) := by sorry

end QAlgorithms.Papers.Independent.Gap48.P25

end
Source
The Security of Practical Quantum Key Distribution, arXiv:0802.4155v3, PDF p. 31, IV.C.1–2: calibrated statistical identities, Eqs.(44)–(48)
Read-back

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

Let ppp be a sequence of real masses pn≥0p_n\ge0pn​≥0, indexed by all natural numbers including zero, whose series sums to 111. Let ν,b,pd\nu,b,p_dν,b,pd​ be real numbers with ν>0\nu>0ν>0, 0≤b≤10\le b\le10≤b≤1, and 0≤pd<1/20\le p_d<1/20≤pd​<1/2, and let f,e:N→Rf,e:\mathbb N\to\mathbb Rf,e:N→R satisfy 0≤fn,en≤10\le f_n,e_n\le10≤fn​,en​≤1 for every nnn and f0=0f_0=0f0​=0. Define the photon and dark contributions, their total, the photon fraction, and the photon error contribution by

Pn=νpnfnb,Dn=νpn(1−fnb)(2pd),R=∑n=0∞(Pn+Dn),Y=∑n=0∞PnR,δ=1−Y2,W=∑n=0∞Pnen,D=∑n=0∞Dn.P_n=\nu p_n f_n b,\qquad D_n=\nu p_n(1-f_n b)(2p_d),\qquad R=\sum_{n=0}^{\infty}(P_n+D_n),\qquad Y=\frac{\sum_{n=0}^{\infty}P_n}{R},\qquad \delta=\frac{1-Y}{2},\qquad W=\sum_{n=0}^{\infty}P_ne_n,\qquad D=\sum_{n=0}^{\infty}D_n.Pn​=νpn​fn​b,Dn​=νpn​(1−fn​b)(2pd​),R=n=0∑∞​(Pn​+Dn​),Y=R∑n=0∞​Pn​​,δ=21−Y​,W=n=0∑∞​Pn​en​,D=n=0∑∞​Dn​.

Assume R>0R>0R>0. Then

0≤Y≤1,Y=1−2pdν/R1−2pd,DR=2δ,W+D/2R=WR+δ.0\le Y\le1,\qquad Y=\frac{1-2p_d\nu/R}{1-2p_d},\qquad \frac{D}{R}=2\delta,\qquad \frac{W+D/2}{R}=\frac WR+\delta.0≤Y≤1,Y=1−2pd​1−2pd​ν/R​,RD​=2δ,RW+D/2​=RW​+δ.

If additionally Y>0Y>0Y>0, then (W+D/2)/R=Y(W/(RY))+δ(W+D/2)/R=Y\bigl(W/(RY)\bigr)+\delta(W+D/2)/R=Y(W/(RY))+δ. The unconditional conclusions allow Y=0Y=0Y=0; the final identity is asserted only under the additional strict positivity condition.

Human review
  • Endorsed by Shuze Chen · Oct 21, 2026

    Confirmed by the moderator at approval.

  • Endorsed by Goku · Oct 21, 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