Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Tight bounds for private communication over bosonic Gaussian channels based on teleportation simulation with optimal finite resources

Definition
QAlgorithms_Papers_Independent_Intakethird100_P17

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

prove2quantum

The objects the formalized statements of Tight bounds for private communication over bosonic Gaussian channels based on teleportation simulation with optimal finite resources (arXiv:1808.00608) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • GaussianInConvention: Appendix A, PDF p.7, (A1): rescale the pinned [q,p]=2i coordinates. Covariance here is the centered half anticommutator, and the vacuum is I/(2κ). Finite photon energy is explicit; no moment identity is a field.
  • quadratureKernel: Appendix A, PDF p.7, before (A1): actual quadrature coefficients on the finite-particle core. The momentum has +i on creation and -i on annihilation.
  • polynomialKernelProduct: Appendix A, PDF p.7, (A9): composition of finite-band polynomial kernels. The theorem asserts convergence of the internal sums, not a default tsum value.
  • PolynomialProductSummable: Domain guard for the unbounded polynomial product on its common core.
  • quadraticKernel: Appendix A, PDF p.7, (A7): coefficients of the genuine q^T A q.
  • densityCoefficients: Actual Fock matrix coefficients of a bounded oscillator density.
  • PolynomialPairingSummable: Absolute convergence of the density-weighted polynomial trace pairing.
  • polynomialExpectation: Appendix A, PDF p.7, (A8)--(A12): actual polynomial expectation. Only the summable domain is used, and the source theorem concludes that domain.
  • quadraticVariance: Appendix A, PDF p.7, (A12): variance formed from actual first/second moments.
  • complexCovariance: Entrywise real-to-complex conversion of a finite covariance matrix.
  • StrictGibbsCovariance: Appendix A, PDF p.7, finite Gibbs domain of (A2)--(A4). Strict uncertainty is equivalent to all symplectic eigenvalues >1/(2κ). No pure-mode log or inverse-coth endpoint is assigned a finite value.
  • gibbsSeriesTerm: Appendix A, PDF p.7, (A4): analytic arcoth series in inverse powers. The strict Gibbs domain has spectral radius of the inverse <1; convergence is concluded in each theorem using this definition. No normality of sigma*iOmega is assumed, and no covariance commutativity is required.
  • kappaGibbsMatrix: Appendix A, PDF p.7, (A4): G=2iOmega arcoth(2κ sigma iOmega). The result is real on the stated physical domain.
  • gaussianSigma: Appendix A, PDF p.7, (A19): the covariance expression, independently of operator entropy. Strict uncertainty makes the determinant positive real.
  • oscillatorLogCutoff: Appendix A, PDF pp.7--8: bounded support-log regularization in BITS. At faithful but not bounded-below densities the unbounded logarithm is obtained only through the explicitly convergent density-weighted limits below.
  • crossEntropyCutoff: Appendix A, PDF p.7, (A18): independently computed regularized cross entropy.
  • oscillatorCrossEntropy: The actual cross entropy on its convergent finite-log-moment domain. Each identity below also states convergence, excluding off-domain limUnder defaults.
  • oscillatorRelativeEntropy: Independently defined Umegaki relative entropy, in bits, retaining +infinity. The accepted finite-projection definition is used, not the Gaussian formula.
  • relativeVarianceCutoff: Appendix A, PDF p.7, (A14)--(A15): square of the ACTUAL centered operator log difference at a bounded cutoff, without assuming the density operators commute.
  • oscillatorRelativeVariance: Relative entropy variance on the domain where the density-weighted squared support-log difference converges. Convergence is part of the theorem.
  • gaussianVarianceExpression: Appendix A, PDF p.8, final unnumbered formula. The Omega term is positive as printed: replacing the last Omega by its transpose would reverse its sign.
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P30
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P13
import Definitions.Def_QAlgorithms_Papers_Independent_Next46_P17

/-!
arXiv:1808.00608v4, Appendix A, PDF pp.7--8.
Actual oscillator densities use the pinned complete Fock Weyl system. All
paper-local definitions have bodies; no entropy or moment identity is assumed.
Operational nodes and source ambiguities are recorded individually in HARD.md.
Compilation and independent moderation are pending the separate handoff.
-/

noncomputable section

open Filter
open scoped BigOperators ENNReal ComplexOrder Topology Matrix

attribute [local instance] Classical.propDecidable Classical.decEq

namespace QAlgorithms.Papers.Independent.Intakethird100.P17

open QAlgorithms.Papers.Independent.Claimed50.P30
  (ModeOperator Quadrature PhaseVector numberKet trace TraceSummable
   FockWeylSystem GaussianMoments canonicalOmega)
open QAlgorithms.Papers.Independent.Dependency50.P13 (IsTraceClass relativeEntropy)
open QAlgorithms.Papers.Independent.Next46.P17 (raisingCoefficient)

/-- Appendix A, PDF p.7, (A1): rescale the pinned [q,p]=2i coordinates.
Covariance here is the centered half anticommutator, and the vacuum is I/(2κ).
Finite photon energy is explicit; no moment identity is a field. -/
def GaussianInConvention {n : ℕ} (W : FockWeylSystem n) (kappa : ℝ)
    (rho : ModeOperator n) (u : PhaseVector n)
    (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) : Prop :=
  GaussianMoments W rho (Real.sqrt (2 * kappa) • u) ((2 * kappa) • sigma) ∧
    Summable (fun j : Fin n → ℕ =>
      (∑ m, (j m : ℝ)) * (inner ℂ (numberKet j) (rho (numberKet j))).re)

/-- Appendix A, PDF p.7, before (A1): actual quadrature coefficients on the
finite-particle core. The momentum has +i on creation and -i on annihilation. -/
def quadratureKernel {n : ℕ} (kappa : ℝ) (a : Quadrature n)
    (i j : Fin n → ℕ) : ℂ :=
  if a.2 = 0 then
    (raisingCoefficient a.1 i j + raisingCoefficient a.1 j i) /
      (Real.sqrt (2 * kappa) : ℂ)
  else Complex.I *
    (raisingCoefficient a.1 i j - raisingCoefficient a.1 j i) /
      (Real.sqrt (2 * kappa) : ℂ)

/-- Appendix A, PDF p.7, (A9): composition of finite-band polynomial kernels.
The theorem asserts convergence of the internal sums, not a default tsum value. -/
def polynomialKernelProduct {n : ℕ}
    (A B : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) :
    Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ :=
  fun i j => ∑' k, A i k * B k j

/-- Domain guard for the unbounded polynomial product on its common core. -/
def PolynomialProductSummable {n : ℕ}
    (A B : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) : Prop :=
  ∀ i j, Summable (fun k => A i k * B k j)

/-- Appendix A, PDF p.7, (A7): coefficients of the genuine q^T A q. -/
def quadraticKernel {n : ℕ} (kappa : ℝ)
    (A : Matrix (Quadrature n) (Quadrature n) ℝ) :
    Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ :=
  fun i j => ∑ a, ∑ b, (A a b : ℂ) *
    polynomialKernelProduct (quadratureKernel kappa a) (quadratureKernel kappa b) i j

/-- Actual Fock matrix coefficients of a bounded oscillator density. -/
def densityCoefficients {n : ℕ} (rho : ModeOperator n) :
    Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ :=
  fun i j => rho (numberKet j) i

/-- Absolute convergence of the density-weighted polynomial trace pairing. -/
def PolynomialPairingSummable {n : ℕ} (rho : ModeOperator n)
    (A : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) : Prop :=
  Summable (fun z : (Fin n → ℕ) × (Fin n → ℕ) =>
    densityCoefficients rho z.1 z.2 * A z.2 z.1)

/-- Appendix A, PDF p.7, (A8)--(A12): actual polynomial expectation.
Only the summable domain is used, and the source theorem concludes that domain. -/
def polynomialExpectation {n : ℕ} (rho : ModeOperator n)
    (A : Matrix (Fin n → ℕ) (Fin n → ℕ) ℂ) : ℝ :=
  (∑' z : (Fin n → ℕ) × (Fin n → ℕ),
    densityCoefficients rho z.1 z.2 * A z.2 z.1).re

/-- Appendix A, PDF p.7, (A12): variance formed from actual first/second moments. -/
def quadraticVariance {n : ℕ} (rho : ModeOperator n) (kappa : ℝ)
    (A : Matrix (Quadrature n) (Quadrature n) ℝ) : ℝ :=
  let Q := quadraticKernel kappa A
  polynomialExpectation rho (polynomialKernelProduct Q Q) - polynomialExpectation rho Q ^ 2

/-- Entrywise real-to-complex conversion of a finite covariance matrix. -/
def complexCovariance {n : ℕ} (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) :
    Matrix (Quadrature n) (Quadrature n) ℂ := fun i j => (sigma i j : ℂ)

/-- Appendix A, PDF p.7, finite Gibbs domain of (A2)--(A4).
Strict uncertainty is equivalent to all symplectic eigenvalues >1/(2κ).
No pure-mode log or inverse-coth endpoint is assigned a finite value. -/
def StrictGibbsCovariance {n : ℕ} (kappa : ℝ)
    (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) : Prop :=
  sigma.transpose = sigma ∧ sigma.PosDef ∧
    ((kappa : ℂ) • complexCovariance sigma +
      (Complex.I / 2) • complexCovariance (canonicalOmega n)).PosDef

/-- Appendix A, PDF p.7, (A4): analytic arcoth series in inverse powers.
The strict Gibbs domain has spectral radius of the inverse <1; convergence is
concluded in each theorem using this definition. No normality of sigma*iOmega
is assumed, and no covariance commutativity is required. -/
def gibbsSeriesTerm {n : ℕ} (kappa : ℝ)
    (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) (k : ℕ) :
    Matrix (Quadrature n) (Quadrature n) ℂ :=
  let B := ((2 * kappa : ℝ) : ℂ) •
    (complexCovariance sigma * (Complex.I • complexCovariance (canonicalOmega n)))
  (((2 * k + 1 : ℕ) : ℂ)⁻¹) • (B⁻¹) ^ (2 * k + 1)

/-- Appendix A, PDF p.7, (A4): G=2iOmega arcoth(2κ sigma iOmega).
The result is real on the stated physical domain. -/
def kappaGibbsMatrix {n : ℕ} (kappa : ℝ)
    (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) :
    Matrix (Quadrature n) (Quadrature n) ℝ :=
  let G := (2 * Complex.I) •
    (complexCovariance (canonicalOmega n) * ∑' k, gibbsSeriesTerm kappa sigma k)
  fun i j => (G i j).re

/-- Appendix A, PDF p.7, (A19): the covariance expression, independently of
operator entropy. Strict uncertainty makes the determinant positive real. -/
def gaussianSigma {n : ℕ} (kappa : ℝ)
    (sigma1 sigma2 : Matrix (Quadrature n) (Quadrature n) ℝ)
    (delta : PhaseVector n) : ℝ :=
  let G := kappaGibbsMatrix kappa sigma2
  let Zsquare := ((kappa : ℂ) • complexCovariance sigma2 +
    (Complex.I / 2) • complexCovariance (canonicalOmega n)).det.re
  (Real.log Zsquare + kappa * (sigma1 * G).trace +
    kappa * ∑ i, ∑ j, delta i * G i j * delta j) / (2 * Real.log 2)

/-- Appendix A, PDF pp.7--8: bounded support-log regularization in BITS.
At faithful but not bounded-below densities the unbounded logarithm is obtained
only through the explicitly convergent density-weighted limits below. -/
def oscillatorLogCutoff {n : ℕ} (rho : ModeOperator n) (k : ℕ) : ModeOperator n := by
  letI : ContinuousFunctionalCalculus ℝ (ModeOperator n) IsSelfAdjoint :=
    IsSelfAdjoint.instContinuousFunctionalCalculus (A := ModeOperator n)
  exact cfc (fun x : ℝ => Real.log (max x (1 / ((k : ℝ) + 1))) / Real.log 2) rho

/-- Appendix A, PDF p.7, (A18): independently computed regularized cross entropy. -/
def crossEntropyCutoff {n : ℕ} (rho sigma : ModeOperator n) (k : ℕ) : ℝ :=
  -(trace (rho.comp (oscillatorLogCutoff sigma k))).re

/-- The actual cross entropy on its convergent finite-log-moment domain.
Each identity below also states convergence, excluding off-domain limUnder defaults. -/
def oscillatorCrossEntropy {n : ℕ} (rho sigma : ModeOperator n) : ℝ :=
  Filter.limUnder atTop (crossEntropyCutoff rho sigma)

/-- Independently defined Umegaki relative entropy, in bits, retaining +infinity.
The accepted finite-projection definition is used, not the Gaussian formula. -/
def oscillatorRelativeEntropy {n : ℕ} (rho sigma : ModeOperator n)
    (hrho : IsTraceClass rho) (hsigma : IsTraceClass sigma) : ℝ≥0∞ :=
  relativeEntropy ⟨rho, hrho⟩ ⟨sigma, hsigma⟩ / ENNReal.ofReal (Real.log 2)

/-- Appendix A, PDF p.7, (A14)--(A15): square of the ACTUAL centered operator
log difference at a bounded cutoff, without assuming the density operators commute. -/
def relativeVarianceCutoff {n : ℕ} (rho sigma : ModeOperator n)
    (mean : ℝ) (k : ℕ) : ℝ :=
  let L := oscillatorLogCutoff rho k - oscillatorLogCutoff sigma k -
    (mean : ℂ) • (1 : ModeOperator n)
  (trace (rho.comp (L.comp L))).re

/-- Relative entropy variance on the domain where the density-weighted squared
support-log difference converges. Convergence is part of the theorem. -/
def oscillatorRelativeVariance {n : ℕ} (rho sigma : ModeOperator n) (mean : ℝ) : ℝ :=
  Filter.limUnder atTop (relativeVarianceCutoff rho sigma mean)

/-- Appendix A, PDF p.8, final unnumbered formula. The Omega term is positive
as printed: replacing the last Omega by its transpose would reverse its sign. -/
def gaussianVarianceExpression {n : ℕ} (kappa : ℝ)
    (sigma1 sigma2 : Matrix (Quadrature n) (Quadrature n) ℝ)
    (delta : PhaseVector n) : ℝ :=
  let G1 := kappaGibbsMatrix kappa sigma1
  let G2 := kappaGibbsMatrix kappa sigma2
  let D := G1 - G2
  let Omega := canonicalOmega n
  let B := (8 * kappa ^ 2) • (G2 * sigma1 * G2)
  (4 * kappa ^ 2 * (sigma1 * D * sigma1 * D).trace +
    (D * Omega * D * Omega).trace + ∑ i, ∑ j, delta i * B i j * delta j) /
    (2 * (2 * Real.log 2) ^ 2)

example {n : ℕ} (m : Fin n) (j : Fin n → ℕ) :
    raisingCoefficient m (Function.update j m (j m + 1)) j =
      (Real.sqrt ((j m : ℝ) + 1) : ℂ) := by
  simp [raisingCoefficient]

example {n : ℕ} (kappa : ℝ) (m : Fin n) (i j : Fin n → ℕ) :
    quadratureKernel kappa (m, 1) i j = Complex.I *
      (raisingCoefficient m i j - raisingCoefficient m j i) /
      (Real.sqrt (2 * kappa) : ℂ) := by
  simp [quadratureKernel]

example {n : ℕ} (rho : ModeOperator n) (k : ℕ) :
    relativeVarianceCutoff rho rho 0 k = 0 := by
  simp [relativeVarianceCutoff, trace,
    QAlgorithms.Papers.Independent.Dependency50.P13.operatorTrace]

example {n : ℕ} (kappa : ℝ) (sigma : Matrix (Quadrature n) (Quadrature n) ℝ) :
    gaussianVarianceExpression kappa sigma sigma 0 = 0 := by
  simp [gaussianVarianceExpression]

end QAlgorithms.Papers.Independent.Intakethird100.P17
Source
Tight bounds for private communication over bosonic Gaussian channels based on teleportation simulation with optimal finite resources, arXiv:1808.00608v4, https://arxiv.org/abs/1808.00608v4
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P17.GaussianInConvention

For any natural number nnn, including 000, put J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}, H=ℓ2(J;C)H=\ell^2(J;\mathbb C)H=ℓ2(J;C), and Q={0,…,n−1}×{0,1}Q=\{0,\ldots,n-1\}\times\{0,1\}Q={0,…,n−1}×{0,1}, and let ej∈He_j\in Hej​∈H have coordinate 111 at jjj and 000 elsewhere. The argument WWW is a family of bounded unitary operators Wξ:H→HW_\xi:H\to HWξ​:H→H, indexed by ξ∈RQ\xi\in\mathbb R^Qξ∈RQ, satisfying

(Wξej)i=∏m=0n−1D ⁣(−ξ(m,1)+iξ(m,0);im,jm),(W_\xi e_j)_i=\prod_{m=0}^{n-1}D\!\left(-\xi_{(m,1)}+\mathrm i\xi_{(m,0)};i_m,j_m\right),(Wξ​ej​)i​=m=0∏n−1​D(−ξ(m,1)​+iξ(m,0)​;im​,jm​),

where, for z∈Cz\in\mathbb Cz∈C and r,s∈Nr,s\in\mathbb Nr,s∈N,

D(z;r,s)=e−∣z∣2/2r! s!∑h=0min⁡(r,s)zr−h(−z‾)s−hh! (r−h)! (s−h)!.D(z;r,s)=e^{-|z|^2/2}\sqrt{r!\,s!} \sum_{h=0}^{\min(r,s)} \frac{z^{r-h}(-\overline z)^{s-h}}{h!\,(r-h)!\,(s-h)!}.D(z;r,s)=e−∣z∣2/2r!s!​h=0∑min(r,s)​h!(r−h)!(s−h)!zr−h(−z)s−h​.

For any real κ\kappaκ, bounded complex-linear operator ρ:H→H\rho:H\to Hρ:H→H, vector u∈RQu\in\mathbb R^Qu∈RQ, and real matrix σ\sigmaσ indexed by QQQ, this predicate means all of the following: there exist sequences vr,wr∈Hv_r,w_r\in Hvr​,wr​∈H with ∑r∥vr∥∥wr∥<∞\sum_r\|v_r\|\|w_r\|<\infty∑r​∥vr​∥∥wr​∥<∞ such that ρx=∑r⟨wr,x⟩vr\rho x=\sum_r\langle w_r,x\rangle v_rρx=∑r​⟨wr​,x⟩vr​ for every x∈Hx\in Hx∈H; ρ∗=ρ\rho^*=\rhoρ∗=ρ and Re⁡⟨x,ρx⟩≥0\operatorname{Re}\langle x,\rho x\rangle\ge0Re⟨x,ρx⟩≥0 for every xxx; the family (⟨ej,ρej⟩)j∈J(\langle e_j,\rho e_j\rangle)_{j\in J}(⟨ej​,ρej​⟩)j∈J​ is summable and has sum 111; the real matrix V=2κσV=2\kappa\sigmaV=2κσ is symmetric and positive definite, and V+iΩV+\mathrm i\OmegaV+iΩ is positive semidefinite, where

Ω(m,b),(ℓ,c)={1m=ℓ, b=0, c=1,−1m=ℓ, b=1, c=0,0otherwise;\Omega_{(m,b),(\ell,c)} =\begin{cases} 1&m=\ell,\ b=0,\ c=1,\\ -1&m=\ell,\ b=1,\ c=0,\\ 0&\text{otherwise}; \end{cases}Ω(m,b),(ℓ,c)​=⎩⎨⎧​1−10​m=ℓ, b=0, c=1,m=ℓ, b=1, c=0,otherwise;​

for every ξ∈RQ\xi\in\mathbb R^Qξ∈RQ, the family (⟨ej,ρWξej⟩)j∈J(\langle e_j,\rho W_\xi e_j\rangle)_{j\in J}(⟨ej​,ρWξ​ej​⟩)j∈J​ is summable and

∑j∈J⟨ej,ρWξej⟩=exp⁡ ⁣(−12∑a,b∈Qξa(2κσab)ξb+i∑a∈Q2κ uaξa);\sum_{j\in J}\langle e_j,\rho W_\xi e_j\rangle = \exp\!\left( -\frac12\sum_{a,b\in Q}\xi_a(2\kappa\sigma_{ab})\xi_b +\mathrm i\sum_{a\in Q}\sqrt{2\kappa}\,u_a\xi_a \right);j∈J∑​⟨ej​,ρWξ​ej​⟩=exp​−21​a,b∈Q∑​ξa​(2κσab​)ξb​+ia∈Q∑​2κ​ua​ξa​​;

and the real family

((∑m=0n−1jm)Re⁡⟨ej,ρej⟩)j∈J\left(\left(\sum_{m=0}^{n-1}j_m\right) \operatorname{Re}\langle e_j,\rho e_j\rangle\right)_{j\in J}((m=0∑n−1​jm​)Re⟨ej​,ρej​⟩)j∈J​

is summable. There is no explicit sign restriction on κ\kappaκ: the real square root is 000 when 2κ≤02\kappa\le02κ≤0; when n>0n>0n>0 and κ=0\kappa=0κ=0, the required positive definiteness of VVV cannot hold. When n=0n=0n=0, JJJ is a singleton, QQQ is empty, and the last summability condition concerns a single zero term.

QAlgorithms.Papers.Independent.Intakethird100.P17.PolynomialPairingSummable

For any natural number nnn, let J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}, H=ℓ2(J;C)H=\ell^2(J;\mathbb C)H=ℓ2(J;C), and eje_jej​ be the coordinate unit vector at jjj. For any bounded complex-linear operator ρ:H→H\rho:H\to Hρ:H→H and arbitrary complex matrix A:J×J→CA:J\times J\to\mathbb CA:J×J→C, this predicate means that the single family indexed by pairs

((ρej)iAji)(i,j)∈J×J\bigl((\rho e_j)_i A_{ji}\bigr)_{(i,j)\in J\times J}((ρej​)i​Aji​)(i,j)∈J×J​

is summable in C\mathbb CC, equivalently

∑(i,j)∈J×J∣(ρej)iAji∣<∞.\sum_{(i,j)\in J\times J}|(\rho e_j)_i A_{ji}|<\infty.(i,j)∈J×J∑​∣(ρej​)i​Aji​∣<∞.

No positivity or normalization of ρ\rhoρ, or symmetry of AAA, is required. When n=0n=0n=0, the indexing set is a singleton and this predicate holds automatically.

QAlgorithms.Papers.Independent.Intakethird100.P17.PolynomialProductSummable

For any natural number nnn, with J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}, and arbitrary complex matrices A,B:J×J→CA,B:J\times J\to\mathbb CA,B:J×J→C, this predicate means that for every i,j∈Ji,j\in Ji,j∈J, the family

(AikBkj)k∈J\bigl(A_{ik}B_{kj}\bigr)_{k\in J}(Aik​Bkj​)k∈J​

is summable in C\mathbb CC, equivalently ∑k∈J∣AikBkj∣<∞\sum_{k\in J}|A_{ik}B_{kj}|<\infty∑k∈J​∣Aik​Bkj​∣<∞. This is an entrywise requirement and imposes no uniform bound over i,ji,ji,j. When n=0n=0n=0, JJJ is a singleton and the predicate holds for every A,BA,BA,B.

QAlgorithms.Papers.Independent.Intakethird100.P17.StrictGibbsCovariance

For any natural number nnn, real κ\kappaκ, and real matrix σ\sigmaσ indexed by Q={0,…,n−1}×{0,1}Q=\{0,\ldots,n-1\}\times\{0,1\}Q={0,…,n−1}×{0,1}, this predicate means that σT=σ\sigma^{\mathsf T}=\sigmaσT=σ, that σ\sigmaσ is positive definite over R\mathbb RR, and that the complex matrix

κσC+i2ΩC\kappa\sigma_{\mathbb C}+\frac{\mathrm i}{2}\Omega_{\mathbb C}κσC​+2i​ΩC​

is positive definite over C\mathbb CC. Here the subscript C\mathbb CC means entrywise embedding of a real matrix into C\mathbb CC, and

Ω(m,b),(ℓ,c)={1m=ℓ, b=0, c=1,−1m=ℓ, b=1, c=0,0otherwise.\Omega_{(m,b),(\ell,c)} =\begin{cases} 1&m=\ell,\ b=0,\ c=1,\\ -1&m=\ell,\ b=1,\ c=0,\\ 0&\text{otherwise}. \end{cases}Ω(m,b),(ℓ,c)​=⎩⎨⎧​1−10​m=ℓ, b=0, c=1,m=ℓ, b=1, c=0,otherwise.​

Positive definiteness requires the respective Hermitian quadratic forms to be strictly positive on every nonzero vector. There is no separate sign hypothesis on κ\kappaκ. When n=0n=0n=0, both index sets are empty and these conditions hold vacuously for the unique empty matrix, for every κ\kappaκ.

QAlgorithms.Papers.Independent.Intakethird100.P17.complexCovariance

For any natural number nnn, including 000, and arbitrary real matrix σ\sigmaσ indexed by Q={0,…,n−1}×{0,1}Q=\{0,\ldots,n-1\}\times\{0,1\}Q={0,…,n−1}×{0,1}, this definition gives the complex matrix with entries

(σC)ab=σab+i 0(a,b∈Q).(\sigma_{\mathbb C})_{ab}=\sigma_{ab}+\mathrm i\,0 \qquad(a,b\in Q).(σC​)ab​=σab​+i0(a,b∈Q).

It imposes no symmetry, positivity, or covariance condition. When n=0n=0n=0, it is the unique matrix with empty row and column index sets.

QAlgorithms.Papers.Independent.Intakethird100.P17.crossEntropyCutoff

For every natural number nnn, bounded complex linear operators ρ,σ\rho,\sigmaρ,σ on Hn=ℓ2({0,…,n−1}→N;C)\mathcal H_n=\ell^2(\{0,\ldots,n-1\}\to\mathbb N;\mathbb C)Hn​=ℓ2({0,…,n−1}→N;C), and natural number kkk, this defines

−Re⁡∑j:{0,…,n−1}→N⟨ej,ρ(fk(σ)ej)⟩,fk(x)=log⁡(max⁡{x,1/(k+1)})log⁡2,-\operatorname{Re}\sum_{j:\{0,\ldots,n-1\}\to\mathbb N} \left\langle e_j,\rho\bigl(f_k(\sigma)e_j\bigr)\right\rangle, \qquad f_k(x)=\frac{\log(\max\{x,1/(k+1)\})}{\log 2},−Rej:{0,…,n−1}→N∑​⟨ej​,ρ(fk​(σ)ej​)⟩,fk​(x)=log2log(max{x,1/(k+1)})​,

where eje_jej​ is the standard unit vector at jjj, and fk(σ)f_k(\sigma)fk​(σ) is the real continuous functional calculus for self-adjoint operators, with value zero if σ\sigmaσ is not self-adjoint. No density, trace-class, or summability hypotheses are required. The displayed infinite sum is total and has value zero if its complex-valued summands are not summable. The cases k=0k=0k=0 and n=0n=0n=0 are included.

QAlgorithms.Papers.Independent.Intakethird100.P17.densityCoefficients

For any natural number nnn, let J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}, let H=ℓ2(J;C)H=\ell^2(J;\mathbb C)H=ℓ2(J;C), and let eje_jej​ denote the vector with coordinate 111 at jjj and 000 elsewhere. For any bounded complex-linear operator ρ:H→H\rho:H\to Hρ:H→H, this definition gives the complex matrix

ρij=(ρej)i(i,j∈J).\rho_{ij}=(\rho e_j)_i \qquad(i,j\in J).ρij​=(ρej​)i​(i,j∈J).

No positivity, trace-class, or normalization assumption on ρ\rhoρ is imposed. When n=0n=0n=0, JJJ is a singleton, so this matrix has one entry.

QAlgorithms.Papers.Independent.Intakethird100.P17.gaussianSigma

For every natural number nnn, real number κ\kappaκ, real matrices σ1,σ2\sigma_1,\sigma_2σ1​,σ2​ indexed by J={0,…,n−1}×{0,1}J=\{0,\ldots,n-1\}\times\{0,1\}J={0,…,n−1}×{0,1}, and real vector δ\deltaδ indexed by JJJ, this defines

log⁡Z2+κ tr⁡(σ1G)+κ∑a,b∈JδaGabδb2log⁡2,Z2=Re⁡det⁡ ⁣(κσ2+i2Ω),\frac{ \log Z^2+\kappa\,\operatorname{tr}(\sigma_1G) +\kappa\sum_{a,b\in J}\delta_aG_{ab}\delta_b }{2\log 2}, \qquad Z^2=\operatorname{Re}\det\!\left(\kappa\sigma_2+\frac{i}{2}\Omega\right),2log2logZ2+κtr(σ1​G)+κ∑a,b∈J​δa​Gab​δb​​,Z2=Redet(κσ2​+2i​Ω),

where Ω(r,0),(r,1)=1\Omega_{(r,0),(r,1)}=1Ω(r,0),(r,1)​=1, Ω(r,1),(r,0)=−1\Omega_{(r,1),(r,0)}=-1Ω(r,1),(r,0)​=−1, all its other entries are zero, and real matrices inside complex expressions are embedded entrywise. Here

B=2κ σ2(iΩ),Gab=Re⁡[2iΩ∑k=0∞(B−1)2k+12k+1]ab.B=2\kappa\,\sigma_2(i\Omega),\qquad G_{ab}=\operatorname{Re}\left[ 2i\Omega\sum_{k=0}^{\infty}\frac{(B^{-1})^{2k+1}}{2k+1} \right]_{ab}.B=2κσ2​(iΩ),Gab​=Re[2iΩk=0∑∞​2k+1(B−1)2k+1​]ab​.

No covariance, positivity, invertibility, or summability conditions are imposed. Singular matrix inverses and nonsummable matrix series have value zero. The real logarithm is total, with log⁡0=0\log 0=0log0=0 and log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣ for negative xxx, so Z2Z^2Z2 need not be positive. The case n=0n=0n=0 is included and gives zero.

QAlgorithms.Papers.Independent.Intakethird100.P17.gaussianVarianceExpression

For every natural number nnn, real number κ\kappaκ, real matrices σ1,σ2\sigma_1,\sigma_2σ1​,σ2​ indexed by J={0,…,n−1}×{0,1}J=\{0,\ldots,n-1\}\times\{0,1\}J={0,…,n−1}×{0,1}, and real vector δ\deltaδ indexed by JJJ, let Ω(r,0),(r,1)=1\Omega_{(r,0),(r,1)}=1Ω(r,0),(r,1)​=1, Ω(r,1),(r,0)=−1\Omega_{(r,1),(r,0)}=-1Ω(r,1),(r,0)​=−1, with all other entries zero. Regarding real matrices as complex matrices entrywise, set

Br=2κ σr(iΩ),(Gr)ab=Re⁡[2iΩ∑k=0∞(Br−1)2k+12k+1]ab(r=1,2),D=G1−G2.B_r=2\kappa\,\sigma_r(i\Omega),\qquad (G_r)_{ab}=\operatorname{Re}\left[ 2i\Omega\sum_{k=0}^{\infty}\frac{(B_r^{-1})^{2k+1}}{2k+1} \right]_{ab}\quad(r=1,2), \qquad D=G_1-G_2.Br​=2κσr​(iΩ),(Gr​)ab​=Re[2iΩk=0∑∞​2k+1(Br−1​)2k+1​]ab​(r=1,2),D=G1​−G2​.

The definition returns the real number

4κ2tr⁡(σ1Dσ1D)+tr⁡(DΩDΩ)+∑a,b∈Jδa(8κ2G2σ1G2)abδb2(2log⁡2)2.\frac{ 4\kappa^2\operatorname{tr}(\sigma_1D\sigma_1D) +\operatorname{tr}(D\Omega D\Omega) +\displaystyle\sum_{a,b\in J}\delta_a \bigl(8\kappa^2G_2\sigma_1G_2\bigr)_{ab}\delta_b }{ 2(2\log 2)^2 }.2(2log2)24κ2tr(σ1​Dσ1​D)+tr(DΩDΩ)+a,b∈J∑​δa​(8κ2G2​σ1​G2​)ab​δb​​.

No positivity, symmetry, invertibility, convergence, or sign conditions on the result are required. Singular matrix inverses and nonsummable matrix series have value zero. Both κ=0\kappa=0κ=0 and n=0n=0n=0 are allowed; for n=0n=0n=0, the result is zero.

QAlgorithms.Papers.Independent.Intakethird100.P17.gibbsSeriesTerm

For any natural number nnn, real κ\kappaκ, arbitrary real matrix σ\sigmaσ indexed by Q={0,…,n−1}×{0,1}Q=\{0,\ldots,n-1\}\times\{0,1\}Q={0,…,n−1}×{0,1}, and natural number kkk, including 000, embed σ\sigmaσ entrywise into the complex matrix σC\sigma_{\mathbb C}σC​, and define Ω\OmegaΩ by

Ω(m,b),(ℓ,c)={1m=ℓ, b=0, c=1,−1m=ℓ, b=1, c=0,0otherwise.\Omega_{(m,b),(\ell,c)} =\begin{cases} 1&m=\ell,\ b=0,\ c=1,\\ -1&m=\ell,\ b=1,\ c=0,\\ 0&\text{otherwise}. \end{cases}Ω(m,b),(ℓ,c)​=⎩⎨⎧​1−10​m=ℓ, b=0, c=1,m=ℓ, b=1, c=0,otherwise.​

With B=2κ σC(iΩC)B=2\kappa\,\sigma_{\mathbb C}(\mathrm i\Omega_{\mathbb C})B=2κσC​(iΩC​), this definition gives the complex matrix

Tk=12k+1(B−1)2k+1.T_k=\frac{1}{2k+1}(B^{-1})^{2k+1}.Tk​=2k+11​(B−1)2k+1.

The exponent 2k+12k+12k+1 is a positive natural number, and its reciprocal is taken in C\mathbb CC. No symmetry, positivity, strict Gibbs condition, or invertibility hypothesis is imposed. The matrix inverse is total: for a singular BBB, B−1=0B^{-1}=0B−1=0, hence Tk=0T_k=0Tk​=0. When n=0n=0n=0, the result is the unique empty matrix.

QAlgorithms.Papers.Independent.Intakethird100.P17.kappaGibbsMatrix

For every natural number nnn, every real number κ\kappaκ, and every real matrix σ\sigmaσ indexed by J={0,…,n−1}×{0,1}J=\{0,\ldots,n-1\}\times\{0,1\}J={0,…,n−1}×{0,1}, this defines a real matrix Gκ(σ)G_\kappa(\sigma)Gκ​(σ). Let Ω(r,0),(r,1)=1\Omega_{(r,0),(r,1)}=1Ω(r,0),(r,1)​=1, Ω(r,1),(r,0)=−1\Omega_{(r,1),(r,0)}=-1Ω(r,1),(r,0)​=−1, and let all other entries of Ω\OmegaΩ be zero; regard real matrices as complex matrices entrywise. Then

B=2κ σ(iΩ),Gκ(σ)ab=Re⁡[2iΩ∑k=0∞(B−1)2k+12k+1]ab.B=2\kappa\,\sigma(i\Omega),\qquad G_\kappa(\sigma)_{ab} =\operatorname{Re}\left[ 2i\Omega\sum_{k=0}^{\infty} \frac{(B^{-1})^{2k+1}}{2k+1} \right]_{ab}.B=2κσ(iΩ),Gκ​(σ)ab​=Re[2iΩk=0∑∞​2k+1(B−1)2k+1​]ab​.

There are no positivity, symmetry, invertibility, or convergence assumptions. Matrix inversion is total and gives the zero matrix for a singular matrix; the infinite matrix sum is the total mathematical sum, which is zero when the series is not summable. In particular, κ=0\kappa=0κ=0 is allowed. When n=0n=0n=0, the result is the unique matrix on the empty index set.

QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorCrossEntropy

For every natural number nnn and bounded complex linear operators ρ,σ\rho,\sigmaρ,σ on Hn=ℓ2({0,…,n−1}→N;C)\mathcal H_n=\ell^2(\{0,\ldots,n-1\}\to\mathbb N;\mathbb C)Hn​=ℓ2({0,…,n−1}→N;C), this defines the real value selected by the total filter-limit operation, as k→∞k\to\inftyk→∞ through natural numbers, for the sequence

−Re⁡∑j⟨ej,ρ(fk(σ)ej)⟩,fk(x)=log⁡(max⁡{x,1/(k+1)})log⁡2.-\operatorname{Re}\sum_j \left\langle e_j,\rho\bigl(f_k(\sigma)e_j\bigr)\right\rangle, \qquad f_k(x)=\frac{\log(\max\{x,1/(k+1)\})}{\log 2}.−Rej∑​⟨ej​,ρ(fk​(σ)ej​)⟩,fk​(x)=log2log(max{x,1/(k+1)})​.

Here eje_jej​ is the standard number basis, indexed by all functions j:{0,…,n−1}→Nj:\{0,\ldots,n-1\}\to\mathbb Nj:{0,…,n−1}→N, and fk(σ)f_k(\sigma)fk​(σ) uses the real continuous functional calculus for self-adjoint operators, giving zero outside that condition. Nonsummable complex sums have value zero. No convergence, self-adjointness, density, or trace-class hypothesis is imposed: if the real sequence converges, the selected value is its limit, but the definition itself does not assert convergence. The case n=0n=0n=0 is included.

QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorLogCutoff

For every natural number nnn, every bounded complex linear operator ρ\rhoρ on Hn=ℓ2({0,…,n−1}→N;C)\mathcal H_n=\ell^2(\{0,\ldots,n-1\}\to\mathbb N;\mathbb C)Hn​=ℓ2({0,…,n−1}→N;C), and every natural number kkk, this defines the operator obtained by applying the real continuous functional calculus for self-adjoint operators to ρ\rhoρ with the scalar function

fk(x)=log⁡ ⁣(max⁡{x,1k+1})log⁡2.f_k(x)=\frac{\log\!\left(\max\left\{x,\frac1{k+1}\right\}\right)}{\log 2}.fk​(x)=log2log(max{x,k+11​})​.

No self-adjointness, positivity, normalization, or trace-class assumption on ρ\rhoρ is supplied; the total functional-calculus operation gives zero when its self-adjointness condition fails. The scalar logarithm always has a positive argument here. Both k=0k=0k=0, where the lower cutoff is 111, and n=0n=0n=0, where H0\mathcal H_0H0​ is one-dimensional, are included.

QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorRelativeEntropy

For every natural number nnn and bounded complex linear operators ρ,σ\rho,\sigmaρ,σ on Hn=ℓ2({0,…,n−1}→N;C)\mathcal H_n=\ell^2(\{0,\ldots,n-1\}\to\mathbb N;\mathbb C)Hn​=ℓ2({0,…,n−1}→N;C), supplied with proofs that each admits a representation

Tx=∑r=0∞⟨vr,x⟩urfor every x,∑r=0∞∥ur∥ ∥vr∥<∞Tx=\sum_{r=0}^{\infty}\langle v_r,x\rangle u_r \quad\text{for every }x,\qquad \sum_{r=0}^{\infty}\|u_r\|\,\|v_r\|<\inftyTx=r=0∑∞​⟨vr​,x⟩ur​for every x,r=0∑∞​∥ur​∥∥vr​∥<∞

for some sequences ur,vr∈Hnu_r,v_r\in\mathcal H_nur​,vr​∈Hn​, this defines an element of [0,∞][0,\infty][0,∞]. For every natural number mmm and every orthonormal family v0,…,vm−1v_0,\ldots,v_{m-1}v0​,…,vm−1​ in Hn\mathcal H_nHn​, form

Xab=⟨va,ρvb⟩,Yab=⟨va,σvb⟩.X_{ab}=\langle v_a,\rho v_b\rangle,\qquad Y_{ab}=\langle v_a,\sigma v_b\rangle.Xab​=⟨va​,ρvb​⟩,Yab​=⟨va​,σvb​⟩.

Assign this pair the value

R(X,Y)={max⁡ ⁣{0,Re⁡tr⁡(h(X)−Xlog⁡(Y)−X+Y)},ker⁡Y⊆ker⁡X,+∞,otherwise,h(x)=xlog⁡x,R(X,Y)= \begin{cases} \max\!\left\{0,\operatorname{Re}\operatorname{tr} \bigl(h(X)-X\log(Y)-X+Y\bigr)\right\}, &\ker Y\subseteq\ker X,\\ +\infty,&\text{otherwise}, \end{cases} \qquad h(x)=x\log x,R(X,Y)={max{0,Retr(h(X)−Xlog(Y)−X+Y)},+∞,​kerY⊆kerX,otherwise,​h(x)=xlogx,

where h(X)h(X)h(X) and log⁡(Y)\log(Y)log(Y) use the real continuous functional calculus for self-adjoint matrices, including its total default outside the required conditions. The defined value is

sup⁡m,vR(X,Y)log⁡2,\frac{\displaystyle\sup_{m,v}R(X,Y)}{\log 2},log2m,vsup​R(X,Y)​,

with the denominator interpreted in [0,∞][0,\infty][0,∞]. No positivity or trace-one assumption is imposed on either operator. Empty orthonormal families m=0m=0m=0, and the case n=0n=0n=0, are included.

QAlgorithms.Papers.Independent.Intakethird100.P17.oscillatorRelativeVariance

For every natural number nnn, bounded complex linear operators ρ,σ\rho,\sigmaρ,σ on Hn=ℓ2({0,…,n−1}→N;C)\mathcal H_n=\ell^2(\{0,\ldots,n-1\}\to\mathbb N;\mathbb C)Hn​=ℓ2({0,…,n−1}→N;C), and real number mmm, this defines the real value selected by the total filter-limit operation, as k→∞k\to\inftyk→∞, for

Re⁡∑j⟨ej,ρ(Lk(Lkej))⟩,Lk=fk(ρ)−fk(σ)−mI,fk(x)=log⁡(max⁡{x,1/(k+1)})log⁡2.\operatorname{Re}\sum_j\langle e_j,\rho(L_k(L_ke_j))\rangle, \qquad L_k=f_k(\rho)-f_k(\sigma)-mI,\qquad f_k(x)=\frac{\log(\max\{x,1/(k+1)\})}{\log 2}.Rej∑​⟨ej​,ρ(Lk​(Lk​ej​))⟩,Lk​=fk​(ρ)−fk​(σ)−mI,fk​(x)=log2log(max{x,1/(k+1)})​.

Here eje_jej​ is the standard number basis indexed by j:{0,…,n−1}→Nj:\{0,\ldots,n-1\}\to\mathbb Nj:{0,…,n−1}→N; the real functional calculus gives zero outside self-adjoint operators, and a nonsummable complex trace series has value zero. No convergence, density, trace-class, or positivity hypotheses are imposed, and mmm is arbitrary. If the displayed real sequence converges, the definition gives its limit; no such convergence is asserted by the definition. The case n=0n=0n=0 is included.

QAlgorithms.Papers.Independent.Intakethird100.P17.polynomialExpectation

For any natural number nnn, let J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}, H=ℓ2(J;C)H=\ell^2(J;\mathbb C)H=ℓ2(J;C), and eje_jej​ be the coordinate unit vector at jjj. For any bounded complex-linear operator ρ:H→H\rho:H\to Hρ:H→H and arbitrary complex matrix A:J×J→CA:J\times J\to\mathbb CA:J×J→C, this definition assigns the real number

Re⁡(∑(i,j)∈J×J(ρej)iAji).\operatorname{Re}\left( \sum_{(i,j)\in J\times J}(\rho e_j)_i A_{ji} \right).Re​(i,j)∈J×J∑​(ρej​)i​Aji​​.

The sum is over pairs, rather than an iterated sum. There is no summability, positivity, normalization, or symmetry hypothesis; a nonsummable complex family has total sum 000, so the defined value is then 000. When n=0n=0n=0, the expression is the real part of the product of the sole operator coefficient and the sole entry of AAA.

QAlgorithms.Papers.Independent.Intakethird100.P17.polynomialKernelProduct

For any natural number nnn, let J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}. For arbitrary complex matrices A,B:J×J→CA,B:J\times J\to\mathbb CA,B:J×J→C, their defined kernel product is the matrix

(A⋆B)ij=∑k∈JAikBkj(i,j∈J).(A\mathbin{\star}B)_{ij}=\sum_{k\in J}A_{ik}B_{kj} \qquad(i,j\in J).(A⋆B)ij​=k∈J∑​Aik​Bkj​(i,j∈J).

No summability hypothesis is required: each sum is the total infinite-sum operation, which returns 000 if its complex-valued family is not summable. When n=0n=0n=0, JJJ is a singleton and this is ordinary multiplication of the two sole entries.

QAlgorithms.Papers.Independent.Intakethird100.P17.quadraticKernel

For any natural number nnn, real κ\kappaκ, and arbitrary real matrix AAA indexed by Q={0,…,n−1}×{0,1}Q=\{0,\ldots,n-1\}\times\{0,1\}Q={0,…,n−1}×{0,1}, this definition gives a complex matrix on J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1} by

Qκ,A(i,j)=∑a∈Q∑b∈QAab(∑k∈JKκ,a(i,k)Kκ,b(k,j)).Q_{\kappa,A}(i,j) =\sum_{a\in Q}\sum_{b\in Q}A_{ab} \left(\sum_{k\in J}K_{\kappa,a}(i,k)K_{\kappa,b}(k,j)\right).Qκ,A​(i,j)=a∈Q∑​b∈Q∑​Aab​(k∈J∑​Kκ,a​(i,k)Kκ,b​(k,j)).

Here, writing a=(m,c)a=(m,c)a=(m,c), let Rm(i,j)R_m(i,j)Rm​(i,j) be jm+1\sqrt{j_m+1}jm​+1​ if iii is obtained from jjj by increasing coordinate mmm by 111, and 000 otherwise; then

Kκ,(m,0)(i,j)=Rm(i,j)+Rm(j,i)2κ,Kκ,(m,1)(i,j)=i(Rm(i,j)−Rm(j,i))2κ.K_{\kappa,(m,0)}(i,j) =\frac{R_m(i,j)+R_m(j,i)}{\sqrt{2\kappa}}, \qquad K_{\kappa,(m,1)}(i,j) =\frac{\mathrm i(R_m(i,j)-R_m(j,i))}{\sqrt{2\kappa}}.Kκ,(m,0)​(i,j)=2κ​Rm​(i,j)+Rm​(j,i)​,Kκ,(m,1)​(i,j)=2κ​i(Rm​(i,j)−Rm​(j,i))​.

All real factors are embedded in C\mathbb CC. No symmetry or positivity of AAA, or positivity of κ\kappaκ, is assumed. Infinite sums use the total sum operation, returning 000 for a nonsummable family. For κ≤0\kappa\le0κ≤0, division by the zero real square root makes every Kκ,aK_{\kappa,a}Kκ,a​, and hence Qκ,AQ_{\kappa,A}Qκ,A​, zero. For n=0n=0n=0, the outer sums are empty and the resulting matrix is zero on the singleton JJJ.

QAlgorithms.Papers.Independent.Intakethird100.P17.quadraticVariance

For any natural number nnn, real κ\kappaκ, bounded complex-linear operator ρ\rhoρ on H=ℓ2(J;C)H=\ell^2(J;\mathbb C)H=ℓ2(J;C), where J=N{0,…,n−1}J=\mathbb N^{\{0,\ldots,n-1\}}J=N{0,…,n−1}, and arbitrary real matrix AAA indexed by Q={0,…,n−1}×{0,1}Q=\{0,\ldots,n-1\}\times\{0,1\}Q={0,…,n−1}×{0,1}, let eje_jej​ be the coordinate unit vector at jjj. Define Rm(i,j)=jm+1R_m(i,j)=\sqrt{j_m+1}Rm​(i,j)=jm​+1​ when iii is obtained by increasing coordinate mmm of jjj by 111, and Rm(i,j)=0R_m(i,j)=0Rm​(i,j)=0 otherwise, and set

Kκ,(m,0)(i,j)=Rm(i,j)+Rm(j,i)2κ,Kκ,(m,1)(i,j)=i(Rm(i,j)−Rm(j,i))2κ,K_{\kappa,(m,0)}(i,j)=\frac{R_m(i,j)+R_m(j,i)}{\sqrt{2\kappa}}, \qquad K_{\kappa,(m,1)}(i,j)= \frac{\mathrm i(R_m(i,j)-R_m(j,i))}{\sqrt{2\kappa}},Kκ,(m,0)​(i,j)=2κ​Rm​(i,j)+Rm​(j,i)​,Kκ,(m,1)​(i,j)=2κ​i(Rm​(i,j)−Rm​(j,i))​, Qij=∑a,b∈QAab∑k∈JKκ,a(i,k)Kκ,b(k,j).Q_{ij}=\sum_{a,b\in Q}A_{ab} \sum_{k\in J}K_{\kappa,a}(i,k)K_{\kappa,b}(k,j).Qij​=a,b∈Q∑​Aab​k∈J∑​Kκ,a​(i,k)Kκ,b​(k,j).

The defined real number is

Re⁡(∑(i,j)∈J×J(ρej)i(∑k∈JQjkQki))−[Re⁡(∑(i,j)∈J×J(ρej)iQji)]2.\operatorname{Re}\left( \sum_{(i,j)\in J\times J}(\rho e_j)_i \left(\sum_{k\in J}Q_{jk}Q_{ki}\right) \right) - \left[ \operatorname{Re}\left( \sum_{(i,j)\in J\times J}(\rho e_j)_iQ_{ji} \right) \right]^2.Re​(i,j)∈J×J∑​(ρej​)i​(k∈J∑​Qjk​Qki​)​−​Re​(i,j)∈J×J∑​(ρej​)i​Qji​​​2.

No summability, normalization, positivity, or symmetry assumptions are imposed, and no nonnegativity of this value is asserted. Each infinite sum returns 000 if its family is not summable. Real square roots of nonpositive numbers and division by zero have value 000; consequently the expression is 000 when κ≤0\kappa\le0κ≤0. It is also 000 when n=0n=0n=0, because then QQQ is empty and the matrix QijQ_{ij}Qij​ is zero.

QAlgorithms.Papers.Independent.Intakethird100.P17.quadratureKernel

For any natural number nnn, real κ\kappaκ, quadrature index a=(m,b)∈{0,…,n−1}×{0,1}a=(m,b)\in\{0,\ldots,n-1\}\times\{0,1\}a=(m,b)∈{0,…,n−1}×{0,1}, and occupation indices i,j∈N{0,…,n−1}i,j\in\mathbb N^{\{0,\ldots,n-1\}}i,j∈N{0,…,n−1}, define

Rm(i,j)={jm+1i equals j with coordinate m replaced by jm+1,0otherwise.R_m(i,j)= \begin{cases} \sqrt{j_m+1}&i\text{ equals }j\text{ with coordinate }m\text{ replaced by }j_m+1,\\ 0&\text{otherwise}. \end{cases}Rm​(i,j)={jm​+1​0​i equals j with coordinate m replaced by jm​+1,otherwise.​

The complex-valued kernel is

Kκ,(m,b)(i,j)={Rm(i,j)+Rm(j,i)2κb=0,i(Rm(i,j)−Rm(j,i))2κb=1,K_{\kappa,(m,b)}(i,j)= \begin{cases} \dfrac{R_m(i,j)+R_m(j,i)}{\sqrt{2\kappa}}&b=0,\\[6pt] \dfrac{\mathrm i\bigl(R_m(i,j)-R_m(j,i)\bigr)}{\sqrt{2\kappa}}&b=1, \end{cases}Kκ,(m,b)​(i,j)=⎩⎨⎧​2κ​Rm​(i,j)+Rm​(j,i)​2κ​i(Rm​(i,j)−Rm​(j,i))​​b=0,b=1,​

with real quantities embedded in C\mathbb CC. No positivity assumption on κ\kappaκ is imposed. If κ≤0\kappa\le0κ≤0, its real square-root denominator is 000, and complex division by 000 makes the kernel 000. For n=0n=0n=0, there is no quadrature index aaa, so there are no such evaluations.

QAlgorithms.Papers.Independent.Intakethird100.P17.relativeVarianceCutoff

For every natural number nnn, bounded complex linear operators ρ,σ\rho,\sigmaρ,σ on Hn=ℓ2({0,…,n−1}→N;C)\mathcal H_n=\ell^2(\{0,\ldots,n-1\}\to\mathbb N;\mathbb C)Hn​=ℓ2({0,…,n−1}→N;C), real number mmm, and natural number kkk, this defines

Re⁡∑j⟨ej,ρ(Lk(Lkej))⟩,Lk=fk(ρ)−fk(σ)−mI,fk(x)=log⁡(max⁡{x,1/(k+1)})log⁡2.\operatorname{Re}\sum_j\langle e_j,\rho(L_k(L_ke_j))\rangle, \qquad L_k=f_k(\rho)-f_k(\sigma)-mI,\qquad f_k(x)=\frac{\log(\max\{x,1/(k+1)\})}{\log 2}.Rej∑​⟨ej​,ρ(Lk​(Lk​ej​))⟩,Lk​=fk​(ρ)−fk​(σ)−mI,fk​(x)=log2log(max{x,1/(k+1)})​.

The vectors eje_jej​ are the standard unit vectors indexed by all j:{0,…,n−1}→Nj:\{0,\ldots,n-1\}\to\mathbb Nj:{0,…,n−1}→N. Each functional calculus is the real continuous functional calculus for self-adjoint operators and gives zero when its self-adjointness condition fails; a nonsummable complex trace series has value zero. There are no density, trace-class, positivity, summability, or restrictions on mmm, and mmm is not required to be an expectation. The product is Lk∘LkL_k\circ L_kLk​∘Lk​. The cases k=0k=0k=0 and n=0n=0n=0 are included.

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