Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Quantum and private capacities of low-noise channels

Definition
QAlgorithms_Papers_Independent_Intakethird100_P50

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

prove2quantum

The objects the formalized statements of Quantum and private capacities of low-noise channels (arXiv:1705.04335) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • dilationQuantumCapacity: §2.3, PDF p.6: the accepted operational quantum capacity applied to the receiver of this dilation; this adapter does not define capacity by an entropy formula.
  • channelDifference: §2.2, PDF p.4: pointwise difference of channel maps on every input matrix. Distinct from the accepted entropy-difference scalar function.
  • maxEntry: §2.2, PDF p.4: entrywise maximum. Used only on nonempty finite index types.
  • dilationChoiRank: §2.1, PDF p.4: actual Choi rank; redundant environment coordinates do not count.
  • degradability: Definition 5, PDF p.7: minimize full diamond error over all CPTP degraders. On physical nonempty finite systems the admissible set is nonempty in [0,2].
  • quantumCorrection: §3.1, PDF p.8, (9): quantum correction, used for Choi rank e≥2 near eta=0.
  • privateCorrection: §3.1, PDF p.10: full private correction f2, with coefficients 1,4,2,4. Only used on e≥2 and physical eta near zero.
  • quantumPolynomialLeading: Lemma 9, PDF p.9: retained terms, evaluated only for c>0, p>0, e≥2.
  • privatePolynomialLeading: Lemma 10, PDF p.10: retained private-capacity terms.
  • quantumLowNoiseLeading: Theorem 12, PDF p.12: quantum expression; the exponent in the first term is 1.5.
  • privateLowNoiseLeading: Theorem 12, PDF p.12: private expression, with all rank-dependent coefficients.
  • LowNoiseCapacityBounds: Theorem 12, PDF p.12: rank-one and zero-noise branches are exact; log(0) is never used as a physical logarithm. Both eventual remainders are retained.
  • Probabilities: §4, PDF p.12: all probabilities are nonnegative and sum to one.
  • pauliKrausOperator: §4, PDF pp.12–13: canonical Kraus labels I,X,Y,Z.
  • pauliDilation: Appendix A, PDF p.20: four-dimensional canonical Pauli dilation. Zero-probability labels remain present. Physicality is the probability simplex.
  • depolarizingDilation: §4.1, PDF p.13, (15): p is TOTAL nonidentity Pauli error.
  • xzDilation: §4.2, PDF p.17, (25): independent equal X and Z errors, each with probability p.
  • depolarizingHashing: §4.1, PDF p.15: hashing expression, used only eventually near p=0.
  • polynomialWeights: §4.2, PDF p.16: THREE noise polynomials; p0 is their complement.
  • PolynomialNoise: §4.2, PDF p.16: zero constant terms and physical probabilities near zero. Eventual validity suffices for asymptotic conclusions; no constant-rank assumption.
  • linearCoefficient: Theorem 15, PDF p.16: coefficient of the linear term of each polynomial.
  • pairCoefficient: Corollary 16, PDF p.17: coefficient C_p; includes the zero case.
  • modifiedWeights: §4.2, PDF p.16, (22) and footnote 2: EVALUATE each polynomial at p+a_i p².
  • positiveCyclicShift: §4, PDF p.12: positive cyclic shift X|i>=|i+1>.
  • quditPauliDilation: §4, PDF p.12: canonical generalized Pauli dilation in X^k Z^l ordering.
  • composedDilation: §5, PDF pp.18–19: compatible sequential dilation; environment order E1,E2.
  • recoveryDegrader: §5, PDF p.19: the explicitly constructed degrader for the ORIGINAL channel.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_Channel
import Definitions.Def_QAlgorithms_Defs_Ensemble
import Definitions.Def_QAlgorithms_Defs_MatrixDistances
import Definitions.Def_QAlgorithms_Defs_FiniteRegisters
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P67
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P80
import Definitions.Def_QAlgorithms_Papers_Independent_Information100_P69
import Definitions.Def_QAlgorithms_Papers_Independent_Support50_P07
import Definitions.Def_QAlgorithms_Papers_Independent_Frontier18_P17

noncomputable section

open scoped BigOperators Matrix ComplexConjugate ComplexOrder
open Filter

namespace QAlgorithms.Papers.Independent.Intakethird100.P50

variable {A B E F : Type}
variable [Fintype A] [DecidableEq A] [Fintype B] [DecidableEq B]
variable [Fintype E] [DecidableEq E] [Fintype F] [DecidableEq F]

/-! Exact accepted definitions are imported from the pinned library. The independent
paper workflow permits these accepted modules; no unreviewed module is imported.
The source role here is arXiv:1705.04335v3, §2.1–2.3, PDF pp.3–7.
Only genuine paper-specific constructions remain local. -/

open QAlgorithms.Papers.Independent.Information100.P67
  (IsIsometry receiverChannel complementaryChannel maximumCoherentInformation
    dilationPower idealJoint)
open QAlgorithms.Papers.Independent.Information100.P80 (log2)
open QAlgorithms.Papers.Independent.Information100.P69 (diamond choi)
open QAlgorithms.Papers.Independent.Support50.P07 (fullTracePrivateCapacity)
open QAlgorithms.Papers.Independent.Frontier18.P17 (clock)

/-- §2.3, PDF p.6: the accepted operational quantum capacity applied to the
receiver of this dilation; this adapter does not define capacity by an entropy formula. -/
def dilationQuantumCapacity (V : Matrix (B × E) A ℂ) : ℝ :=
  QAlgorithms.Papers.Independent.Information100.P80.quantumCapacity (receiverChannel V)

/-- §2.2, PDF p.4: pointwise difference of channel maps on every input matrix.
Distinct from the accepted entropy-difference scalar function. -/
def channelDifference (Phi Psi : Superop A B) : Superop A B := fun X => Phi X - Psi X

/-- §2.2, PDF p.4: entrywise maximum. Used only on nonempty finite index types. -/
def maxEntry (X : Matrix A B ℂ) : ℝ := sSup {v | ∃ i j, v = ‖X i j‖}

/-- §2.1, PDF p.4: actual Choi rank; redundant environment coordinates do not count. -/
def dilationChoiRank (V : Matrix (B × E) A ℂ) : ℕ := (choi (receiverChannel V)).rank

/-- Definition 5, PDF p.7: minimize full diamond error over all CPTP degraders.
On physical nonempty finite systems the admissible set is nonempty in [0,2]. -/
def degradability (V : Matrix (B × E) A ℂ) : ℝ :=
  sInf {eta : ℝ | 0 ≤ eta ∧ eta ≤ 2 ∧ ∃ D : Superop B E, IsChannel D ∧
    diamond (channelDifference (complementaryChannel V) (D ∘ receiverChannel V)) ≤ eta}

/-- §3.1, PDF p.8, (9): quantum correction, used for Choi rank e≥2 near eta=0. -/
def quantumCorrection (e : ℕ) (eta : ℝ) : ℝ :=
  eta / 2 * log2 ((e : ℝ) - 1) + eta * log2 e + binEntropy2 (eta / 2) +
    (1 + eta / 2) * binEntropy2 (eta / (2 + eta))

/-- §3.1, PDF p.10: full private correction f2, with coefficients 1,4,2,4.
Only used on e≥2 and physical eta near zero. -/
def privateCorrection (e : ℕ) (eta : ℝ) : ℝ :=
  eta * log2 ((e : ℝ) - 1) + 4 * eta * log2 e + 2 * binEntropy2 (eta / 2) +
    4 * (1 + eta / 2) * binEntropy2 (eta / (2 + eta))

/-- Lemma 9, PDF p.9: retained terms, evaluated only for c>0, p>0, e≥2. -/
def quantumPolynomialLeading (e : ℕ) (c r p : ℝ) : ℝ :=
  c * r * p ^ (r - 1) * (-p * log2 p) + c * p ^ r *
    (-log2 c + 1 + (Real.log 2)⁻¹ + log2 e + log2 ((e : ℝ) - 1) / 2)

/-- Lemma 10, PDF p.10: retained private-capacity terms. -/
def privatePolynomialLeading (e : ℕ) (c r p : ℝ) : ℝ :=
  3 * c * r * p ^ (r - 1) * (-p * log2 p) + c * p ^ r *
    (-3 * log2 c + 3 + 3 * (Real.log 2)⁻¹ + log2 ((e : ℝ) - 1) + 4 * log2 e)

/-- Theorem 12, PDF p.12: quantum expression; the exponent in the first term is 1.5. -/
def quantumLowNoiseLeading (e : ℕ) (eps : ℝ) : ℝ :=
  -3 * eps ^ (3 / 2 : ℝ) * log2 eps +
    ((Real.log 2)⁻¹ + log2 e + log2 ((e : ℝ) - 1) / 2) * 2 * eps ^ (3 / 2 : ℝ)

/-- Theorem 12, PDF p.12: private expression, with all rank-dependent coefficients. -/
def privateLowNoiseLeading (e : ℕ) (eps : ℝ) : ℝ :=
  -9 * eps ^ (3 / 2 : ℝ) * log2 eps +
    (3 * (Real.log 2)⁻¹ + 4 * log2 e + log2 ((e : ℝ) - 1)) * 2 * eps ^ (3 / 2 : ℝ)

/-- Theorem 12, PDF p.12: rank-one and zero-noise branches are exact; log(0) is
never used as a physical logarithm. Both eventual remainders are retained. -/
def LowNoiseCapacityBounds (V : Matrix (B × E) A ℂ) (eps Cq Cp : ℝ) : Prop :=
  if dilationChoiRank V = 1 ∨ eps = 0 then
    dilationQuantumCapacity V = maximumCoherentInformation V ∧ fullTracePrivateCapacity V = maximumCoherentInformation V
  else
    |dilationQuantumCapacity V - maximumCoherentInformation V| ≤
      quantumLowNoiseLeading (dilationChoiRank V) eps + Cq * eps ^ (9 / 2 : ℝ) ∧
    |fullTracePrivateCapacity V - maximumCoherentInformation V| ≤
      privateLowNoiseLeading (dilationChoiRank V) eps + Cp * eps ^ (3 : ℕ)

/-- §4, PDF p.12: all probabilities are nonnegative and sum to one. -/
def Probabilities {I : Type} [Fintype I] (w : I → ℝ) : Prop :=
  (∀ i, 0 ≤ w i) ∧ ∑ i, w i = 1

/-- §4, PDF pp.12–13: canonical Kraus labels I,X,Y,Z. -/
def pauliKrausOperator : Fin 4 → Matrix (Fin 2) (Fin 2) ℂ :=
  ![1, !![0, 1; 1, 0], !![0, -Complex.I; Complex.I, 0], !![1, 0; 0, -1]]

/-- Appendix A, PDF p.20: four-dimensional canonical Pauli dilation.
Zero-probability labels remain present. Physicality is the probability simplex. -/
def pauliDilation (w : Fin 4 → ℝ) : Matrix (Fin 2 × Fin 4) (Fin 2) ℂ :=
  fun be a => (Real.sqrt (w be.2) : ℂ) * pauliKrausOperator be.2 be.1 a

/-- §4.1, PDF p.13, (15): p is TOTAL nonidentity Pauli error. -/
def depolarizingDilation (p : ℝ) : Matrix (Fin 2 × Fin 4) (Fin 2) ℂ :=
  pauliDilation ![1 - p, p / 3, p / 3, p / 3]

/-- §4.2, PDF p.17, (25): independent equal X and Z errors, each with probability p. -/
def xzDilation (p : ℝ) : Matrix (Fin 2 × Fin 4) (Fin 2) ℂ :=
  pauliDilation ![(1 - p) ^ 2, p - p ^ 2, p ^ 2, p - p ^ 2]

/-- §4.1, PDF p.15: hashing expression, used only eventually near p=0. -/
def depolarizingHashing (p : ℝ) : ℝ := 1 - binEntropy2 p - p * log2 3

/-- §4.2, PDF p.16: THREE noise polynomials; p0 is their complement. -/
def polynomialWeights (P : Fin 3 → Polynomial ℝ) (p : ℝ) : Fin 4 → ℝ :=
  ![1 - ∑ i, (P i).eval p, (P 0).eval p, (P 1).eval p, (P 2).eval p]

/-- §4.2, PDF p.16: zero constant terms and physical probabilities near zero.
Eventual validity suffices for asymptotic conclusions; no constant-rank assumption. -/
def PolynomialNoise (P : Fin 3 → Polynomial ℝ) : Prop :=
  (∀ i, (P i).eval 0 = 0) ∧ ∃ delta : ℝ, 0 < delta ∧
    ∀ p ∈ Set.Icc 0 delta, Probabilities (polynomialWeights P p)

/-- Theorem 15, PDF p.16: coefficient of the linear term of each polynomial. -/
def linearCoefficient (P : Fin 3 → Polynomial ℝ) (i : Fin 3) : ℝ := (P i).coeff 1

/-- Corollary 16, PDF p.17: coefficient C_p; includes the zero case. -/
def pairCoefficient (P : Fin 3 → Polynomial ℝ) : ℝ :=
  |linearCoefficient P 0 * linearCoefficient P 1 +
    linearCoefficient P 0 * linearCoefficient P 2 +
    linearCoefficient P 1 * linearCoefficient P 2|

/-- §4.2, PDF p.16, (22) and footnote 2: EVALUATE each polynomial at p+a_i p². -/
def modifiedWeights (P : Fin 3 → Polynomial ℝ) (a : Fin 3 → ℝ) (p : ℝ) : Fin 4 → ℝ :=
  let q := fun i : Fin 3 => (P i).eval (p + a i * p ^ 2)
  ![1 - ∑ i, q i, q 0, q 1, q 2]

/-- §4, PDF p.12: positive cyclic shift X|i>=|i+1>. -/
def positiveCyclicShift (d : ℕ) : Matrix (Fin d) (Fin d) ℂ :=
  fun i j => if i.val = (j.val + 1) % d then 1 else 0

/-- §4, PDF p.12: canonical generalized Pauli dilation in X^k Z^l ordering. -/
def quditPauliDilation (d : ℕ) (w : (Fin d × Fin d) → ℝ) :
    Matrix (Fin d × (Fin d × Fin d)) (Fin d) ℂ :=
  fun be a => (Real.sqrt (w be.2) : ℂ) *
    (positiveCyclicShift d ^ be.2.1.val * clock d ^ be.2.2.val) be.1 a

/-! Sanity checks of local conventions, not proofs of source results. -/

example : log2 1 = 0 := by simp [log2]

example (V : Matrix (B × E) A ℂ) (rho : Matrix A A ℂ) :
    complementaryChannel V rho = traceLeft (V * rho * Vᴴ) := rfl

example (Phi : Superop A B) (a c : A) (b d : B) :
    choi Phi (a,b) (c,d) = Phi (Matrix.single a c 1) b d := rfl

example (V : Matrix (B × E) A ℂ)
    (be : (Fin 0 → B) × (Fin 0 → E)) (a : Fin 0 → A) :
    dilationPower V 0 be a = 1 := by simp [dilationPower]

example (m : ℕ) (sigma : Matrix E E ℂ) (i : Fin m) (e f : E) :
    idealJoint m sigma (i,(i,e)) (i,(i,f)) = (m : ℂ)⁻¹ * sigma e f := by
  simp [idealJoint]

example : pauliKrausOperator 1 0 1 = 1 ∧ pauliKrausOperator 2 0 1 = -Complex.I := by constructor <;> rfl

example (P : Fin 3 → Polynomial ℝ) (a : Fin 3 → ℝ) (p : ℝ) :
    modifiedWeights P a p 1 = (P 0).eval (p + a 0 * p ^ 2) := rfl

example : positiveCyclicShift 2 1 0 = 1 := by norm_num [positiveCyclicShift]

/-- The three-dimensional case distinguishes the forward shift from its inverse. -/
example : positiveCyclicShift 3 1 0 = 1 ∧ positiveCyclicShift 3 2 0 = 0 := by
  norm_num [positiveCyclicShift]

example (p : ℝ) : xzDilation p =
    pauliDilation ![(1 - p) ^ 2, p - p ^ 2, p ^ 2, p - p ^ 2] := rfl

example (e : ℕ) : quantumCorrection e 0 = 0 := by
  simp [quantumCorrection, binEntropy2]

example (e : ℕ) : privateCorrection e 0 = 0 := by
  simp [privateCorrection, binEntropy2]

example : depolarizingDilation 0 (0,0) 0 = 1 := by
  norm_num [depolarizingDilation, pauliDilation, pauliKrausOperator]

/-- §5, PDF pp.18–19: compatible sequential dilation; environment order E1,E2. -/
def composedDilation (V : Matrix (B × E) A ℂ) (W : Matrix (A × F) B ℂ) :
    Matrix (A × (E × F)) A ℂ :=
  fun ae a => ∑ b, W (ae.1, ae.2.2) b * V (b, ae.2.1) a

/-- §5, PDF p.19: the explicitly constructed degrader for the ORIGINAL channel. -/
def recoveryDegrader (V : Matrix (B × E) A ℂ) (W : Matrix (A × F) B ℂ) : Superop B E :=
  fun rho => traceRight (complementaryChannel (composedDilation V W) (receiverChannel W rho))

example (V : Matrix (B × E) A ℂ) (W : Matrix (A × F) B ℂ)
    (a c : A) (e : E) (f : F) :
    composedDilation V W (a,(e,f)) c = ∑ b, W (a,f) b * V (b,e) c := rfl

end QAlgorithms.Papers.Independent.Intakethird100.P50
Source
Quantum and private capacities of low-noise channels, arXiv:1705.04335v3, https://arxiv.org/abs/1705.04335v3
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P50.LowNoiseCapacityBounds

For finite types A,B,EA,B,EA,B,E with decidable equality, any complex matrix V:(A→B×E)V:(A\to B\times E)V:(A→B×E), and arbitrary real numbers ε,Cq,Cp\varepsilon,C_q,C_pε,Cq​,Cp​, define Φ(X)=Tr⁡E(VXV†)\Phi(X)=\operatorname{Tr}_E(VXV^\dagger)Φ(X)=TrE​(VXV†), Φc(X)=Tr⁡B(VXV†)\Phi^c(X)=\operatorname{Tr}_B(VXV^\dagger)Φc(X)=TrB​(VXV†), and e=rank⁡Je=\operatorname{rank}Je=rankJ, where J(a,b),(a′,b′)=[Φ(∣a⟩⟨a′∣)]b,b′J_{(a,b),(a',b')}=[\Phi(|a\rangle\langle a'|)]_{b,b'}J(a,b),(a′,b′)​=[Φ(∣a⟩⟨a′∣)]b,b′​. Put S(X)=−Re⁡tr⁡[f(X)]S(X)=-\operatorname{Re}\operatorname{tr}[f(X)]S(X)=−Retr[f(X)], where f(x)=xlog⁡2xf(x)=x\log_2xf(x)=xlog2​x is applied by functional calculus, and put

Ic=sup⁡ρ≥0, tr⁡ρ=1(S(Φ(ρ))−S(Φc(ρ))).I_c=\sup_{\rho\ge0,\ \operatorname{tr}\rho=1} \bigl(S(\Phi(\rho))-S(\Phi^c(\rho))\bigr).Ic​=ρ≥0, trρ=1sup​(S(Φ(ρ))−S(Φc(ρ))).

Let QQQ be the real supremum of nonnegative rates rrr admitting a sequence, indexed by every natural nnn, of positive message dimensions MnM_nMn​ and encoding and decoding channels En,Dn\mathcal E_n,\mathcal D_nEn​,Dn​, with

Re⁡tr⁡ ⁣[ΩMn(id⊗(Dn∘Φ⊗n∘En))(ΩMn)]→1,∀δ>0,r−δ≤log⁡2Mnn eventually,\operatorname{Re}\operatorname{tr}\!\left[ \Omega_{M_n}(\mathrm{id}\otimes(\mathcal D_n\circ\Phi^{\otimes n}\circ\mathcal E_n))(\Omega_{M_n}) \right]\to1, \qquad \forall\delta>0,\quad r-\delta\le\frac{\log_2M_n}{n}\ \text{eventually},Retr[ΩMn​​(id⊗(Dn​∘Φ⊗n∘En​))(ΩMn​​)]→1,∀δ>0,r−δ≤nlog2​Mn​​ eventually,

where ΩM=M−1∑i,j∣i,i⟩⟨j,j∣\Omega_M=M^{-1}\sum_{i,j}|i,i\rangle\langle j,j|ΩM​=M−1∑i,j​∣i,i⟩⟨j,j∣. Let PPP be the real supremum of nonnegative rates rrr such that, for every ϵ0>0\epsilon_0>0ϵ0​>0 and slack s>0s>0s>0, there is a positive n0n_0n0​ such that every n≥n0n\ge n_0n≥n0​ admits M>0M>0M>0, density matrices ρm\rho_mρm​ on AnA^nAn, and positive semidefinite measurement matrices DjD_jDj​ on BnB^nBn satisfying ∑jDj=I\sum_jD_j=I∑j​Dj​=I, for which r−s≤log⁡2M/nr-s\le\log_2M/nr−s≤log2​M/n and ∥ω−KM,σ∥1<ϵ0\|\omega-\mathcal K_{M,\sigma}\|_1<\epsilon_0∥ω−KM,σ​∥1​<ϵ0​. In this condition, Vn(b,e;a)=∏i=0n−1V(bi,ei;ai)V_n(b,e;a)=\prod_{i=0}^{n-1}V(b_i,e_i;a_i)Vn​(b,e;a)=∏i=0n−1​V(bi​,ei​;ai​),

ω(m,j,e),(m′,j′,e′)={1M∑b,cDj(c,b)(VnρmVn†)(b,e),(c,e′),m=m′, j=j′,0,otherwise,σ=Tr⁡m,jω,KM,σ=1M∑m∣m,m⟩⟨m,m∣⊗σ.\omega_{(m,j,e),(m',j',e')} = \begin{cases} \displaystyle\frac1M\sum_{b,c}D_j(c,b) (V_n\rho_mV_n^\dagger)_{(b,e),(c,e')}, &m=m',\ j=j',\\ 0,&\text{otherwise}, \end{cases} \quad \sigma=\operatorname{Tr}_{m,j}\omega, \quad \mathcal K_{M,\sigma}=\frac1M\sum_m|m,m\rangle\langle m,m|\otimes\sigma.ω(m,j,e),(m′,j′,e′)​=⎩⎨⎧​M1​b,c∑​Dj​(c,b)(Vn​ρm​Vn†​)(b,e),(c,e′)​,0,​m=m′, j=j′,otherwise,​σ=Trm,j​ω,KM,σ​=M1​m∑​∣m,m⟩⟨m,m∣⊗σ.

A density matrix is positive semidefinite with trace 111; a channel has a finite Kraus representation 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; partial traces sum diagonal indices of the discarded registers; ∥X∥1=Re⁡tr⁡X†X\|X\|_1=\operatorname{Re}\operatorname{tr}\sqrt{X^\dagger X}∥X∥1​=RetrX†X​; and tensor powers act on matrix units by Φ⊗n(∣x⟩⟨y∣)=⨂iΦ(∣xi⟩⟨yi∣)\Phi^{\otimes n}(|x\rangle\langle y|)=\bigotimes_i\Phi(|x_i\rangle\langle y_i|)Φ⊗n(∣x⟩⟨y∣)=⨂i​Φ(∣xi​⟩⟨yi​∣). With these quantities, the proposition defined here is Q=IcQ=I_cQ=Ic​ and P=IcP=I_cP=Ic​ when e=1e=1e=1 or ε=0\varepsilon=0ε=0; otherwise it is the conjunction

∣Q−Ic∣≤−3ε3/2log⁡2ε+2(1ln⁡2+log⁡2e+12log⁡2(e−1))ε3/2+Cqε9/2,∣P−Ic∣≤−9ε3/2log⁡2ε+2(3ln⁡2+4log⁡2e+log⁡2(e−1))ε3/2+Cpε3.\begin{aligned} |Q-I_c| &\le -3\varepsilon^{3/2}\log_2\varepsilon +2\left(\frac1{\ln2}+\log_2e+\frac12\log_2(e-1)\right)\varepsilon^{3/2} +C_q\varepsilon^{9/2},\\ |P-I_c| &\le -9\varepsilon^{3/2}\log_2\varepsilon +2\left(\frac3{\ln2}+4\log_2e+\log_2(e-1)\right)\varepsilon^{3/2} +C_p\varepsilon^3. \end{aligned}∣Q−Ic​∣∣P−Ic​∣​≤−3ε3/2log2​ε+2(ln21​+log2​e+21​log2​(e−1))ε3/2+Cq​ε9/2,≤−9ε3/2log2​ε+2(ln23​+4log2​e+log2​(e−1))ε3/2+Cp​ε3.​

No isometry condition on VVV, nonemptiness condition on the registers, or sign condition on ε,Cq,Cp\varepsilon,C_q,C_pε,Cq​,Cp​ is imposed. Real suprema use their total convention, including 000 for an empty defining set; log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2 uses ln⁡0=0\ln0=0ln0=0 and the absolute-value convention for negative arguments; fractional powers are total real powers. Zero-length tensor products and division by zero also use their total definitions.

QAlgorithms.Papers.Independent.Intakethird100.P50.PolynomialNoise

For three real polynomials P0,P1,P2P_0,P_1,P_2P0​,P1​,P2​, this predicate means that Pi(0)=0P_i(0)=0Pi​(0)=0 for every i∈{0,1,2}i\in\{0,1,2\}i∈{0,1,2}, and that there exists a real δ>0\delta>0δ>0 such that, for every real p∈[0,δ]p\in[0,\delta]p∈[0,δ], all four entries of

w(p)=(1−∑i=02Pi(p), P0(p), P1(p), P2(p))w(p)=\left(1-\sum_{i=0}^{2}P_i(p),\ P_0(p),\ P_1(p),\ P_2(p)\right)w(p)=(1−i=0∑2​Pi​(p), P0​(p), P1​(p), P2​(p))

are nonnegative and ∑k=03wk(p)=1\sum_{k=0}^{3}w_k(p)=1∑k=03​wk​(p)=1. Both endpoints of the interval are included; no condition is required outside that interval.

QAlgorithms.Papers.Independent.Intakethird100.P50.Probabilities

For a finite type III and a function w:I→Rw:I\to\mathbb Rw:I→R, this predicate means that wi≥0w_i\ge 0wi​≥0 for every i∈Ii\in Ii∈I and ∑i∈Iwi=1\sum_{i\in I}w_i=1∑i∈I​wi​=1. If III is empty, the sum is 000, so the predicate cannot hold.

QAlgorithms.Papers.Independent.Intakethird100.P50.channelDifference

For arbitrary types A,BA,BA,B and arbitrary maps Φ,Ψ\Phi,\PsiΦ,Ψ from complex A×AA\times AA×A matrices to complex B×BB\times BB×B matrices, their defined difference is the map

X⟼Φ(X)−Ψ(X).X\longmapsto\Phi(X)-\Psi(X).X⟼Φ(X)−Ψ(X).

There are no assumptions of linearity, positivity, normalization, finiteness, or nonemptiness.

QAlgorithms.Papers.Independent.Intakethird100.P50.composedDilation

For types A,B,E,FA,B,E,FA,B,E,F, with BBB finite so that the following sum is defined, and arbitrary complex matrices VVV with rows B×EB\times EB×E and columns AAA, and WWW with rows A×FA\times FA×F and columns BBB, this defines the complex matrix ZZZ with rows A×(E×F)A\times(E\times F)A×(E×F) and columns AAA by

Z(aout,(e,f)),ain=∑b∈BW(aout,f),b V(b,e),ain.Z_{(a_{\mathrm{out}},(e,f)),a_{\mathrm{in}}} =\sum_{b\in B}W_{(a_{\mathrm{out}},f),b}\,V_{(b,e),a_{\mathrm{in}}}.Z(aout​,(e,f)),ain​​=b∈B∑​W(aout​,f),b​V(b,e),ain​​.

Neither input matrix is assumed to be an isometry. Empty index types are allowed; in particular, if BBB is empty, every entry of ZZZ is zero.

QAlgorithms.Papers.Independent.Intakethird100.P50.degradability

For finite types A,B,EA,B,EA,B,E with decidable equality and any complex matrix V:(A→B×E)V:(A\to B\times E)V:(A→B×E), put ΦV(X)=Tr⁡E(VXV†)\Phi_V(X)=\operatorname{Tr}_E(VXV^\dagger)ΦV​(X)=TrE​(VXV†) and ΦVc(X)=Tr⁡B(VXV†)\Phi_V^c(X)=\operatorname{Tr}_B(VXV^\dagger)ΦVc​(X)=TrB​(VXV†). This defines

inf⁡{η∈[0,2]:∃D,D(X)=∑jKjXKj†,∑jKj†Kj=IB,∥ΦVc−D∘ΦV∥⋄≤η},\inf\left\{\eta\in[0,2]: \exists\mathcal D,\quad \mathcal D(X)=\sum_jK_jXK_j^\dagger,\quad \sum_jK_j^\dagger K_j=I_B,\quad \|\Phi_V^c-\mathcal D\circ\Phi_V\|_\diamond\le\eta \right\},inf{η∈[0,2]:∃D,D(X)=j∑​Kj​XKj†​,j∑​Kj†​Kj​=IB​,∥ΦVc​−D∘ΦV​∥⋄​≤η},

where the Kraus family is finite, each KjK_jKj​ has rows indexed by EEE and columns indexed by BBB, and the particular norm used is

∥Λ∥⋄=sup⁡X on A×A∥X∥1≤1∥(idA⊗Λ)(X)∥1,∥X∥1=Re⁡tr⁡X†X.\|\Lambda\|_\diamond = \sup_{\substack{X\text{ on }A\times A\\\|X\|_1\le1}} \|(\mathrm{id}_A\otimes\Lambda)(X)\|_1, \qquad \|X\|_1=\operatorname{Re}\operatorname{tr}\sqrt{X^\dagger X}.∥Λ∥⋄​=X on A×A∥X∥1​≤1​sup​∥(idA​⊗Λ)(X)∥1​,∥X∥1​=RetrX†X​.

The square root is defined by functional calculus, and idA⊗Λ\mathrm{id}_A\otimes\LambdaidA​⊗Λ applies Λ\LambdaΛ to each reference-indexed matrix block. The two partial traces sum the corresponding diagonal register indices. No isometry condition on VVV is assumed, and the set of admissible η\etaη need not be nonempty; an empty defining set gives value 000 under the real infimum convention.

QAlgorithms.Papers.Independent.Intakethird100.P50.depolarizingDilation

For every real ppp, this defines a complex matrix with entries

V(b,k),a=wk (Pk)b,a,(w0,w1,w2,w3)=(1−p,p3,p3,p3),V_{(b,k),a}=\sqrt{w_k}\,(P_k)_{b,a}, \qquad (w_0,w_1,w_2,w_3)=\left(1-p,\frac p3,\frac p3,\frac p3\right),V(b,k),a​=wk​​(Pk​)b,a​,(w0​,w1​,w2​,w3​)=(1−p,3p​,3p​,3p​),

where a,b∈{0,1}a,b\in\{0,1\}a,b∈{0,1}, k∈{0,1,2,3}k\in\{0,1,2,3\}k∈{0,1,2,3}, and

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\\mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

There is no restriction 0≤p≤10\le p\le10≤p≤1; any negative weight contributes a zero square-root factor.

QAlgorithms.Papers.Independent.Intakethird100.P50.depolarizingHashing

For every real ppp, this defines the real number

1−−plog⁡p−(1−p)log⁡(1−p)log⁡2−p log⁡3log⁡2.1-\frac{-p\log p-(1-p)\log(1-p)}{\log2} -p\,\frac{\log3}{\log2}.1−log2−plogp−(1−p)log(1−p)​−plog2log3​.

Here log⁡\loglog is the total real logarithm used by the code: log⁡0=0\log0=0log0=0, and for negative arguments it equals the logarithm of the absolute value. Thus the definition includes p=0p=0p=0, p=1p=1p=1, and values outside [0,1][0,1][0,1].

QAlgorithms.Papers.Independent.Intakethird100.P50.dilationChoiRank

For finite types A,B,EA,B,EA,B,E, with decidable equality on AAA, and any complex matrix V:(A→B×E)V:(A\to B\times E)V:(A→B×E), this defines the matrix rank of JVJ_VJV​, indexed by A×BA\times BA×B, whose entries are

(JV)(a,b),(a′,b′)=[Tr⁡E ⁣(V∣a⟩⟨a′∣V†)]b,b′.(J_V)_{(a,b),(a',b')} = \left[\operatorname{Tr}_E\!\left(V|a\rangle\langle a'|V^\dagger\right)\right]_{b,b'}.(JV​)(a,b),(a′,b′)​=[TrE​(V∣a⟩⟨a′∣V†)]b,b′​.

Here ∣a⟩⟨a′∣|a\rangle\langle a'|∣a⟩⟨a′∣ is the matrix with its only nonzero entry, equal to 111, at (a,a′)(a,a')(a,a′), and (Tr⁡EY)b,b′=∑eY(b,e),(b′,e)(\operatorname{Tr}_EY)_{b,b'}=\sum_eY_{(b,e),(b',e)}(TrE​Y)b,b′​=∑e​Y(b,e),(b′,e)​. There is no normalization factor in JVJ_VJV​, and no isometry or nonemptiness hypothesis is imposed on VVV or the registers; empty indexing cases therefore include rank 000.

QAlgorithms.Papers.Independent.Intakethird100.P50.dilationQuantumCapacity

For finite types A,B,EA,B,EA,B,E, with decidable equality on AAA and BBB, and any complex matrix VVV with rows indexed by B×EB\times EB×E and columns indexed by AAA, define ΦV(X)=Tr⁡E(VXV†)\Phi_V(X)=\operatorname{Tr}_E(VXV^\dagger)ΦV​(X)=TrE​(VXV†), where (Tr⁡EY)b,b′=∑eY(b,e),(b′,e)(\operatorname{Tr}_E Y)_{b,b'}=\sum_eY_{(b,e),(b',e)}(TrE​Y)b,b′​=∑e​Y(b,e),(b′,e)​. The value defined here is the real supremum of all r≥0r\ge0r≥0 for which there is, for every natural number nnn, a positive integer MnM_nMn​, an encoding channel En\mathcal E_nEn​ from MnM_nMn​-dimensional matrices to matrices on AnA^nAn, and a decoding channel Dn\mathcal D_nDn​ from matrices on BnB^nBn to MnM_nMn​-dimensional matrices, such that

Re⁡tr⁡ ⁣[ΩMn (id⊗(Dn∘ΦV⊗n∘En))(ΩMn)]⟶1,∀δ>0,r−δ≤log⁡2Mnn eventually.\operatorname{Re}\operatorname{tr}\!\left[ \Omega_{M_n}\, (\mathrm{id}\otimes(\mathcal D_n\circ\Phi_V^{\otimes n}\circ\mathcal E_n))(\Omega_{M_n}) \right]\longrightarrow1, \qquad \forall\delta>0,\quad r-\delta\le\frac{\log_2M_n}{n}\ \text{eventually}.Retr[ΩMn​​(id⊗(Dn​∘ΦV⊗n​∘En​))(ΩMn​​)]⟶1,∀δ>0,r−δ≤nlog2​Mn​​ eventually.

Here AnA^nAn means the functions {0,…,n−1}→A\{0,\ldots,n-1\}\to A{0,…,n−1}→A, ΩM=M−1∑i,j=0M−1∣i,i⟩⟨j,j∣\Omega_M=M^{-1}\sum_{i,j=0}^{M-1}|i,i\rangle\langle j,j|ΩM​=M−1∑i,j=0M−1​∣i,i⟩⟨j,j∣, and a channel means a map X↦∑jKjXKj†X\mapsto\sum_jK_jXK_j^\daggerX↦∑j​Kj​XKj†​ with finitely many Kraus matrices satisfying ∑jKj†Kj=I\sum_jK_j^\dagger K_j=I∑j​Kj†​Kj​=I. The tensor power is determined by ΦV⊗n(∣x⟩⟨x′∣)=⨂iΦV(∣xi⟩⟨xi′∣)\Phi_V^{\otimes n}(|x\rangle\langle x'|)=\bigotimes_i\Phi_V(|x_i\rangle\langle x'_i|)ΦV⊗n​(∣x⟩⟨x′∣)=⨂i​ΦV​(∣xi​⟩⟨xi′​∣), and id⊗Φ\mathrm{id}\otimes\Phiid⊗Φ applies Φ\PhiΦ to each matrix block indexed by the reference register. The logarithm is log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2. The sequence includes n=0n=0n=0, with the empty tensor product and total division x/0=0x/0=0x/0=0, although the conditions concern n→∞n\to\inftyn→∞. No isometry condition on VVV is assumed; the supremum is taken in R\mathbb RR, with its total supremum convention, including value 000 for the empty set.

QAlgorithms.Papers.Independent.Intakethird100.P50.linearCoefficient

For any three real polynomials P0,P1,P2P_0,P_1,P_2P0​,P1​,P2​ and any index i∈{0,1,2}i\in\{0,1,2\}i∈{0,1,2}, this defines the coefficient of the degree-one monomial in PiP_iPi​. It is 000 when that monomial is absent; no noise or probability condition is assumed.

QAlgorithms.Papers.Independent.Intakethird100.P50.maxEntry

For types A,BA,BA,B and any complex matrix XXX indexed by A×BA\times BA×B, this defines

sup⁡{ ∣Xi,j∣:i∈A, j∈B }\sup\{\,|X_{i,j}|:i\in A,\ j\in B\,\}sup{∣Xi,j​∣:i∈A, j∈B}

as a real supremum, where ∣⋅∣|\cdot|∣⋅∣ is the complex modulus. For finite nonempty A,BA,BA,B, this is the largest entry modulus; if either indexing type is empty, the defining set is empty and the value is 000. The definition itself imposes no finiteness or boundedness hypothesis and uses the total real supremum operation.

QAlgorithms.Papers.Independent.Intakethird100.P50.modifiedWeights

For any three real polynomials P0,P1,P2P_0,P_1,P_2P0​,P1​,P2​, any three real numbers a0,a1,a2a_0,a_1,a_2a0​,a1​,a2​, and any real ppp, set qi=Pi(p+aip2)q_i=P_i(p+a_ip^2)qi​=Pi​(p+ai​p2). This defines the four real weights

(1−q0−q1−q2, q0, q1, q2),(1-q_0-q_1-q_2,\ q_0,\ q_1,\ q_2),(1−q0​−q1​−q2​, q0​, q1​, q2​),

indexed in that order by 0,1,2,30,1,2,30,1,2,3. No sign, normalization, or noise hypotheses are imposed on the inputs; individual weights may be negative.

QAlgorithms.Papers.Independent.Intakethird100.P50.pairCoefficient

For any three real polynomials P0,P1,P2P_0,P_1,P_2P0​,P1​,P2​, let cic_ici​ be the coefficient of the degree-one monomial in PiP_iPi​, with absent coefficients equal to 000. This defines

∣c0c1+c0c2+c1c2∣.\left|c_0c_1+c_0c_2+c_1c_2\right|.∣c0​c1​+c0​c2​+c1​c2​∣.

The absolute value applies to the entire sum, and no assumptions are imposed on the polynomials.

QAlgorithms.Papers.Independent.Intakethird100.P50.pauliDilation

For any four real numbers w0,w1,w2,w3w_0,w_1,w_2,w_3w0​,w1​,w2​,w3​, this defines a complex matrix VVV whose rows are indexed by pairs (b,k)∈{0,1}×{0,1,2,3}(b,k)\in\{0,1\}\times\{0,1,2,3\}(b,k)∈{0,1}×{0,1,2,3} and whose columns are indexed by a∈{0,1}a\in\{0,1\}a∈{0,1}, by

V(b,k),a=wk (Pk)b,a,(P0,P1,P2,P3)=((1001),(0110),(0−ii0),(100−1)).V_{(b,k),a}=\sqrt{w_k}\,(P_k)_{b,a}, \qquad (P_0,P_1,P_2,P_3)= \left( \begin{pmatrix}1&0\\0&1\end{pmatrix}, \begin{pmatrix}0&1\\1&0\end{pmatrix}, \begin{pmatrix}0&-\mathrm i\\\mathrm i&0\end{pmatrix}, \begin{pmatrix}1&0\\0&-1\end{pmatrix} \right).V(b,k),a​=wk​​(Pk​)b,a​,(P0​,P1​,P2​,P3​)=((10​01​),(01​10​),(0i​−i0​),(10​0−1​)).

The real square root is regarded as a complex number. No nonnegativity or normalization assumption is imposed on www; the real square root used here is 000 for negative inputs.

QAlgorithms.Papers.Independent.Intakethird100.P50.pauliKrausOperator

This defines four complex 2×22\times22×2 matrices, indexed by 0,1,2,30,1,2,30,1,2,3, with rows and columns indexed by 0,10,10,1:

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1),P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\qquad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\qquad P_2=\begin{pmatrix}0&-\mathrm i\\ \mathrm i&0\end{pmatrix},\qquad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix},P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​),

where i2=−1\mathrm i^2=-1i2=−1.

QAlgorithms.Papers.Independent.Intakethird100.P50.polynomialWeights

For any three real polynomials P0,P1,P2P_0,P_1,P_2P0​,P1​,P2​ and any real ppp, this defines the four real weights

(1−P0(p)−P1(p)−P2(p), P0(p), P1(p), P2(p)),\left(1-P_0(p)-P_1(p)-P_2(p),\ P_0(p),\ P_1(p),\ P_2(p)\right),(1−P0​(p)−P1​(p)−P2​(p), P0​(p), P1​(p), P2​(p)),

indexed in that order by 0,1,2,30,1,2,30,1,2,3. No assumptions on the polynomials or on ppp are imposed, and the weights need not be nonnegative.

QAlgorithms.Papers.Independent.Intakethird100.P50.positiveCyclicShift

For every natural number ddd, this defines a complex d×dd\times dd×d matrix SSS, with row and column indices i,j∈{0,…,d−1}i,j\in\{0,\ldots,d-1\}i,j∈{0,…,d−1}, by

Si,j={1,i=(j+1) mod d,0,otherwise.S_{i,j}= \begin{cases} 1,&i=(j+1)\bmod d,\\ 0,&\text{otherwise}. \end{cases}Si,j​={1,0,​i=(j+1)modd,otherwise.​

For d=0d=0d=0, the matrix has no rows or columns, so no remainder expression is evaluated at an entry; for d=1d=1d=1, it is the 1×11\times11×1 matrix with entry 111.

QAlgorithms.Papers.Independent.Intakethird100.P50.privateCorrection

For every natural number eee and every real number η\etaη, this defines

ηlog⁡2(e−1)+4ηlog⁡2e+2H2 ⁣(η2)+4(1+η2)H2 ⁣(η2+η),\eta\log_2(e-1)+4\eta\log_2e +2H_2\!\left(\frac{\eta}{2}\right) +4\left(1+\frac{\eta}{2}\right) H_2\!\left(\frac{\eta}{2+\eta}\right),ηlog2​(e−1)+4ηlog2​e+2H2​(2η​)+4(1+2η​)H2​(2+ηη​),

where log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2 and H2(x)=[−xln⁡x−(1−x)ln⁡(1−x)]/ln⁡2H_2(x)=[-x\ln x-(1-x)\ln(1-x)]/\ln2H2​(x)=[−xlnx−(1−x)ln(1−x)]/ln2, using the total real logarithm. There are no restrictions on η\etaη and no lower bound on eee. In particular, e=0,1e=0,1e=0,1 and η=−2\eta=-2η=−2 are included: e−1e-1e−1 is real subtraction, ln⁡0=0\ln0=0ln0=0, negative logarithm arguments are interpreted through absolute value, and division by zero is 000.

QAlgorithms.Papers.Independent.Intakethird100.P50.privateLowNoiseLeading

For every natural number eee and every real number ε\varepsilonε, this defines

−9ε3/2log⁡2ε+2(3ln⁡2+4log⁡2e+log⁡2(e−1))ε3/2.-9\varepsilon^{3/2}\log_2\varepsilon +2\left( \frac3{\ln2}+4\log_2e+\log_2(e-1) \right)\varepsilon^{3/2}.−9ε3/2log2​ε+2(ln23​+4log2​e+log2​(e−1))ε3/2.

The exponent 3/23/23/2 uses the total real-power operation, and log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2 uses the total real logarithm, including ln⁡0=0\ln0=0ln0=0 and ln⁡x=ln⁡∣x∣\ln x=\ln|x|lnx=ln∣x∣ for negative xxx. All real ε\varepsilonε and all natural eee, including e=0,1e=0,1e=0,1, are allowed; e−1e-1e−1 denotes real subtraction.

QAlgorithms.Papers.Independent.Intakethird100.P50.privatePolynomialLeading

For every natural number eee and all real numbers c,r,pc,r,pc,r,p, this defines

3crp r−1(−plog⁡2p)+cpr(−3log⁡2c+3+3ln⁡2+log⁡2(e−1)+4log⁡2e).3crp^{\,r-1}(-p\log_2p) +cp^r\left( -3\log_2c+3+\frac3{\ln2} +\log_2(e-1)+4\log_2e \right).3crpr−1(−plog2​p)+cpr(−3log2​c+3+ln23​+log2​(e−1)+4log2​e).

The logarithm means log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2, with ln⁡0=0\ln0=0ln0=0 and ln⁡x=ln⁡∣x∣\ln x=\ln|x|lnx=ln∣x∣ for negative xxx; the powers with exponents r−1r-1r−1 and rrr are total real powers. The quantities r−1r-1r−1 and e−1e-1e−1 are real differences. The definition includes arbitrary real c,r,pc,r,pc,r,p, including nonpositive values, and every natural eee, including 000 and 111.

QAlgorithms.Papers.Independent.Intakethird100.P50.quantumCorrection

For every natural number eee and every real number η\etaη, this defines

η2log⁡2(e−1)+ηlog⁡2e+H2 ⁣(η2)+(1+η2)H2 ⁣(η2+η),\frac{\eta}{2}\log_2(e-1)+\eta\log_2e +H_2\!\left(\frac{\eta}{2}\right) +\left(1+\frac{\eta}{2}\right) H_2\!\left(\frac{\eta}{2+\eta}\right),2η​log2​(e−1)+ηlog2​e+H2​(2η​)+(1+2η​)H2​(2+ηη​),

where log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2 and H2(x)=[−xln⁡x−(1−x)ln⁡(1−x)]/ln⁡2H_2(x)=[-x\ln x-(1-x)\ln(1-x)]/\ln2H2​(x)=[−xlnx−(1−x)ln(1−x)]/ln2, using the total real logarithm. No conditions such as e≥2e\ge2e≥2 or 0≤η≤20\le\eta\le20≤η≤2 are imposed. The subtraction e−1e-1e−1 occurs in R\mathbb RR; logarithms at zero have value 000, logarithms at negative arguments use ln⁡x=ln⁡∣x∣\ln x=\ln|x|lnx=ln∣x∣, and division by zero is 000, including when η=−2\eta=-2η=−2.

QAlgorithms.Papers.Independent.Intakethird100.P50.quantumLowNoiseLeading

For every natural number eee and every real number ε\varepsilonε, this defines

−3ε3/2log⁡2ε+2(1ln⁡2+log⁡2e+12log⁡2(e−1))ε3/2.-3\varepsilon^{3/2}\log_2\varepsilon +2\left( \frac1{\ln2}+\log_2e+\frac12\log_2(e-1) \right)\varepsilon^{3/2}.−3ε3/2log2​ε+2(ln21​+log2​e+21​log2​(e−1))ε3/2.

Here 3/23/23/2 is a real exponent and the power is the total real-power operation; log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2 uses the total real logarithm, with ln⁡0=0\ln0=0ln0=0 and ln⁡x=ln⁡∣x∣\ln x=\ln|x|lnx=ln∣x∣ for negative xxx. There is no nonnegativity or smallness condition on ε\varepsilonε, and no restriction excluding e=0,1e=0,1e=0,1; e−1e-1e−1 is computed in R\mathbb RR.

QAlgorithms.Papers.Independent.Intakethird100.P50.quantumPolynomialLeading

For every natural number eee and all real numbers c,r,pc,r,pc,r,p, this defines

crp r−1(−plog⁡2p)+cpr(−log⁡2c+1+1ln⁡2+log⁡2e+12log⁡2(e−1)).crp^{\,r-1}(-p\log_2p) +cp^r\left( -\log_2c+1+\frac1{\ln2} +\log_2e+\frac12\log_2(e-1) \right).crpr−1(−plog2​p)+cpr(−log2​c+1+ln21​+log2​e+21​log2​(e−1)).

Here log⁡2x=ln⁡x/ln⁡2\log_2x=\ln x/\ln2log2​x=lnx/ln2, with the total real logarithm, so ln⁡0=0\ln0=0ln0=0 and ln⁡x=ln⁡∣x∣\ln x=\ln|x|lnx=ln∣x∣ for negative xxx. The exponents r−1r-1r−1 and rrr use the total real-power operation, and both subtractions r−1r-1r−1 and e−1e-1e−1 occur in R\mathbb RR. No positivity, integrality, or smallness assumptions are imposed on c,r,pc,r,pc,r,p, and e=0,1e=0,1e=0,1 are included.

QAlgorithms.Papers.Independent.Intakethird100.P50.quditPauliDilation

For every natural number ddd and every real function www on {0,…,d−1}2\{0,\ldots,d-1\}^2{0,…,d−1}2, this defines a complex matrix VwV_wVw​ with columns indexed by a∈{0,…,d−1}a\in\{0,\ldots,d-1\}a∈{0,…,d−1} and rows indexed by (b,(j,k))∈{0,…,d−1}3(b,(j,k))\in\{0,\ldots,d-1\}^3(b,(j,k))∈{0,…,d−1}3, by

(Vw)(b,(j,k)),a=w(j,k) (SdjCdk)b,a,(Sd)u,v={1u=(v+1) mod d,0otherwise,(Cd)u,v={e2πiu/du=v,0u≠v.(V_w)_{(b,(j,k)),a}=\sqrt{w(j,k)}\,(S_d^jC_d^k)_{b,a}, \qquad (S_d)_{u,v}= \begin{cases}1&u=(v+1)\bmod d,\\0&\text{otherwise},\end{cases} \qquad (C_d)_{u,v}= \begin{cases}e^{2\pi i u/d}&u=v,\\0&u\ne v.\end{cases}(Vw​)(b,(j,k)),a​=w(j,k)​(Sdj​Cdk​)b,a​,(Sd​)u,v​={10​u=(v+1)modd,otherwise,​(Cd​)u,v​={e2πiu/d0​u=v,u=v.​

The square root is the real square root, subsequently regarded as complex, and the exponents j,kj,kj,k are the natural representatives of the indices. There is no assumption that the weights are nonnegative or sum to one, and no isometry assertion; negative weights contribute zero because the real square root of a negative number is zero. The definition includes d=0d=0d=0, when both matrix index sets are empty, and d=1d=1d=1.

QAlgorithms.Papers.Independent.Intakethird100.P50.recoveryDegrader

For finite types A,B,E,FA,B,E,FA,B,E,F and arbitrary complex matrices V:B×E←AV:B\times E\leftarrow AV:B×E←A and W:A×F←BW:A\times F\leftarrow BW:A×F←B, this defines a map from complex matrices on BBB to complex matrices on EEE. Set

Z(a,(e,f)),a′=∑b∈BW(a,f),bV(b,e),a′,σa,a′=∑f∈F(WρW†)(a,f),(a′,f).Z_{(a,(e,f)),a'}=\sum_{b\in B}W_{(a,f),b}V_{(b,e),a'}, \qquad \sigma_{a,a'}=\sum_{f\in F}(W\rho W^\dagger)_{(a,f),(a',f)}.Z(a,(e,f)),a′​=b∈B∑​W(a,f),b​V(b,e),a′​,σa,a′​=f∈F∑​(WρW†)(a,f),(a′,f)​.

Its value at any complex matrix ρ\rhoρ on BBB is the matrix

DV,W(ρ)e,e′=∑f∈F∑a∈A(ZσZ†)(a,(e,f)),(a,(e′,f)).\mathcal D_{V,W}(\rho)_{e,e'} =\sum_{f\in F}\sum_{a\in A} (Z\sigma Z^\dagger)_{(a,(e,f)),(a,(e',f))}.DV,W​(ρ)e,e′​=f∈F∑​a∈A∑​(ZσZ†)(a,(e,f)),(a,(e′,f))​.

Thus the definition first traces out FFF from WρW†W\rho W^\daggerWρW†, applies conjugation by the displayed ZZZ, and traces out its AAA and FFF coordinates. No isometry, positivity, normalization, or channel hypothesis is imposed on V,W,ρV,W,\rhoV,W,ρ, and the definition itself makes no channel assertion. Empty finite types are included, with empty sums equal to zero.

QAlgorithms.Papers.Independent.Intakethird100.P50.xzDilation

For every real ppp, this defines a complex matrix with entries

V(b,k),a=wk (Pk)b,a,(w0,w1,w2,w3)=((1−p)2, p−p2, p2, p−p2),V_{(b,k),a}=\sqrt{w_k}\,(P_k)_{b,a}, \qquad (w_0,w_1,w_2,w_3)=\bigl((1-p)^2,\ p-p^2,\ p^2,\ p-p^2\bigr),V(b,k),a​=wk​​(Pk​)b,a​,(w0​,w1​,w2​,w3​)=((1−p)2, p−p2, p2, p−p2),

where a,b∈{0,1}a,b\in\{0,1\}a,b∈{0,1}, k∈{0,1,2,3}k\in\{0,1,2,3\}k∈{0,1,2,3}, and

P0=(1001),P1=(0110),P2=(0−ii0),P3=(100−1).P_0=\begin{pmatrix}1&0\\0&1\end{pmatrix},\quad P_1=\begin{pmatrix}0&1\\1&0\end{pmatrix},\quad P_2=\begin{pmatrix}0&-\mathrm i\\\mathrm i&0\end{pmatrix},\quad P_3=\begin{pmatrix}1&0\\0&-1\end{pmatrix}.P0​=(10​01​),P1​=(01​10​),P2​=(0i​−i0​),P3​=(10​0−1​).

No restriction on ppp is imposed; the real square root of a negative weight is taken to be 000.

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