Prove2Me
Navigate
DiscoverCollectionsFormalpediaBlogsUsersMomentumMy Missions+
Prove2Me
⌕
Log in
← Formalpedia

Definitions for Discriminating quantum-optical beam-splitter channels with number-diagonal signal states: Applications to quantum reading and target detection

Definition
QAlgorithms_Papers_Independent_Intakethird100_P54

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

prove2quantum

The objects the formalized statements of Discriminating quantum-optical beam-splitter channels with number-diagonal signal states: Applications to quantum reading and target detection (arXiv:1105.4063) are stated in. Each was checked against the paper by a blind read-back and a moderator.

  • opticalKernel: Fock coefficients of an operator; validity is supplied by its construction.
  • fockVector: A full Fock vector from coefficients. Off l² only, the total extension is zero. The stated construction clauses rule this branch out for every physical input.
  • Splitter: A single-mode vacuum-environment splitter, PDF p.1, (1).
  • Splitter.r: Nonnegative field reflection magnitude.
  • Splitter.t: Nonnegative field transmission magnitude.
  • Splitter.field: Complex returned field coefficient, PDF p.1.
  • Channels: Two candidate splitters; the prior is built into the error definition.
  • ValidSplitters: The paper's ordering and distinctness assumptions, PDF pp.1–2.
  • Reading: Section IV's common phase removed by the same known unitary.
  • splitterOverlap: The common real overlap coefficient of (17),(49).
  • photons: Total signal photon number, with no idler energy charged, PDF p.4, (18).
  • lossAmplitude: Vacuum-environment binomial amplitude (12). The phase multiplies the entire input occupation, not just the returned photon count. The guard precedes subtraction.
  • sectorCoefficients: Loss-sector coefficients (11),(25); idlers are retained without change.
  • lossSector: The unnormalized retained state in a loss sector, PDF p.4, (25).
  • returned: Actual returned density after tracing inaccessible independent vacuum modes, PDF pp.1,3–5, (2),(11),(32). Every series is over full Fock labels.
  • signalLaw: Reduced signal number law, PDF p.4, before (17).
  • totalLaw: Total signal photon law (18), distinct from the multimode occupation law.
  • HasSignalEnergy: Finite mean SIGNAL energy, PDF pp.6–7. Summability rules out totalized tsum.
  • IsNDS: Source definition of NDS: the signal marginal is number diagonal, PDF p.4. For pure probes this is equivalent to orthogonal occupied idler Schmidt vectors.
  • purificationMarginal: Reduced density of a normalized countable-ancilla purification.
  • PurifiesOpticalKernel: Purification as actual norm-one amplitudes, not an assumed fidelity value.
  • squaredOpticalFidelity: Squared Uhlmann fidelity (4),(19), as the supremum over all purifications. For physical densities here this set is nonempty and bounded in [0,1]; the countable ancilla suffices because the optical Hilbert spaces are separable.
  • Resolution: A spectral resolution of a normalized positive trace-class density. Only positive-weight vectors are required to be orthonormal. The kernel equality is part of a spectral representation, not an assumption on any source theorem.
  • IsOpticalEffect: A binary measurement effect is a bounded complex-linear positive contraction. E is the decision-zero effect; 1-E is the decision-one effect, PDF pp.3,5.
  • opticalBayesRisk: Binary Bayesian error for actual spectral states and an effect; arbitrary joint signal-idler effects are admitted, including non-projective measurements.
  • bayesError: Operational minimum Bayesian error. For normalized physical densities and priors this is a nonempty infimum of numbers in [0,1]; spectral representations do not constrain the allowed receiver. Equal priors are used except in (42).
  • opticalError: Equal-prior minimum average error, PDF p.3, (5).
  • supportWeightPower: Support spectral power: the zero eigenspace contributes zero even at t=0. This implements the support-projection convention on PDF p.6.
  • chernoffAt: Tr(rho^s sigma^(1-s)), defined independently from NDS formulas (45). For densities and 0<=s<=1 every spectral representation gives the same absolutely convergent nonnegative double series, in [0,1]. Hence this supremum is nonempty and bounded; it does not optimize the probe or the receiver.
  • opticalChernoff: Chernoff coefficient WITHOUT prior factor, including both support endpoints (9).
  • outputFidelity: Pure-input output fidelity, for the actual channels.
  • outputError: Pure-input operational output error.
  • outputChernoffAt: Pure-input Chernoff-type coefficient.
  • outputChernoff: Pure-input optimized Chernoff coefficient.
  • purificationOverlapTerm: Complex overlap series in the TOTAL photon law, PDF p.4, (17).
  • purificationOverlap: Squared overlap (17); the absolute value encloses the WHOLE complex sum.
  • lowerError: Fidelity-based LOWER error envelope (6),(20),(53).
  • lossProbability: Loss-sector probability (27) directly from its actual vector.
  • lossOverlap: Loss-sector overlap (30); the first argument is conjugated.
  • lossOverlapTerm: Explicit complex binomial summand (31), with componentwise loss guard.
  • lossBlock: The at-most-two-dimensional block H_k of (33)–(35).
  • lossBlockDifference: Weighted block difference; equal priors use weights 1,1 for rho0-rho1.
  • BlockReceiver: A complete block spectral receiver, (40)–(42). E is the zero decision. It is projective, preserves every loss block, and decides the positive/negative spectra correctly. Zero eigenvectors may be assigned either decision.
  • sectorRisk: Risk of a concrete effect, directly on unnormalized output loss vectors.
  • ndsChernoffTerm: Term of (45), explicitly omitting either zero-probability sector.
  • opticalPureKernel: Rank-one full-Fock density.
  • fock: A Fock signal with no retained idlers; the empty idler occupation has one value.
  • productOpticalProbe: Independent product probes with arbitrary finite idler count J_m per signal, PDF p.6, (46)–(47). The idler register is grouped, without tracing any factors.
  • opticalCopies: Independent copies of a joint multimode probe, with all signals grouped.
  • coherentOpticalProbe: Product coherent probe; arbitrary classical idler amplitudes are retained.
  • coherentEnergy: Total signal energy of a product coherent state.
  • tmsvCoefficients: Multimode EPR coefficients (73), with N the per-pair SIGNAL energy. Natural powers preserve the vacuum convention 0^0=1.
  • tmsvProbe: Product of M identical EPR pairs, with no photon cutoff.
  • countLaw: Product law of observed RETURNED counts, PDF p.8, (62) with j=N-k.
  • totalCountLaw: Binomial law of the total returned photon count.
  • countingError: Equal-prior overlap of the total count distributions, not an assumed quantum optimum.
  • CountsAs: Photon counting followed by d(j), where true means decision one. E is the zero-decision effect and acts trivially on every retained idler.
  • countThreshold: Interior logarithmic threshold (65); endpoint rules are separate theorems.
  • fockCoefficient: One-photon classical Chernoff coefficient, PDF p.8, (66).
  • fockOptimizer: Explicit optimizer (67); theorem hypotheses restrict it to strict interior channels.
  • numberExponent: Per-photon number-state exponent, with natural logarithm, (66).
  • coherentExponent: Per-photon coherent exponent (68).
  • exponentGain: Exponent gain (69); no claim that it is large in every regime.
  • eprAlpha: Parameters (78)–(81); beta is used only when R1<1.
  • eprBeta: Parameter beta in (79).
  • eprC: Parameter C in (80).
  • eprD: Parameter D in (81).
  • eprDerivativeSign: Equivalent logarithmic case test (84)–(85): no exponent/product ambiguity.
  • eprOptimizer: Optimizer (85)–(86). N=0 and ideal memory are handled outside this definition.
  • classicalReturned: Classical positive-P probes include arbitrary probability measures and Dirac masses. Signal and idler amplitudes may be classically correlated, (3), pp.2,8.
  • vacuumEfficiency: An additional uniform independent vacuum attenuator, PDF p.12, IV.E.2.
  • efficiencyScaled: A uniform receiver efficiency multiplies energy reflectance. Eta=0 is allowed and makes the two effective returned channels identical.
  • attenuate: Vacuum attenuator on a full joint operator kernel, retaining every idler.
  • weightedTraceNorm: Dual trace norm of a weighted Hermitian difference, used in (42). The supremum ranges over ALL bounded self-adjoint contractions. Spectral resolutions merely evaluate traces of the actual two states.
  • kernelProduct: Composition of full Fock operator kernels on the physical trace-class arguments.
  • targetLower: Target-detection universal envelope (59), with fixed real signal energy N.
  • targetCoherent: Coherent target-detection error envelope (60).
Definition code
import Mathlib
import Definitions.Def_QAlgorithms_Defs_RealLabelRegisters
import Definitions.Def_QAlgorithms_Papers_Independent_Claimed50_P47
import Definitions.Def_QAlgorithms_Papers_Independent_Dependency50_P39

/-!
Statements from R. Nair, arXiv:1105.4063v3. PDF pages are cited throughout.
The optical space is full l² of occupation strings, with retained idlers first.
The accepted Hilbert, occupation labels and coherent coefficients are imported.
The remaining optical definitions below have complete bodies; their domains differ
from accepted declarations sharing their former short names. Only source results use sorry.
Fidelity is squared; Q excludes the equal-prior factor 1/2. Spectral zero powers
are support projections. Compilation is delegated to the batch handoff.
-/

noncomputable section

open scoped BigOperators ComplexConjugate
open MeasureTheory Filter

namespace QAlgorithms.Papers.Independent.Intakethird100.P54

open QAlgorithms.RealRegisters (Hilbert)
open QAlgorithms.Papers.Independent.Claimed50.P47 (Occupation)
open QAlgorithms.Papers.Independent.Dependency50.P39 (coherentCoefficient)

attribute [local instance] Classical.propDecidable
attribute [local instance] Classical.decEq

/-- Fock coefficients of an operator; validity is supplied by its construction. -/
abbrev opticalKernel (X : Type) := X → X → ℂ

/-- A full Fock vector from coefficients. Off l² only, the total extension is zero.
The stated construction clauses rule this branch out for every physical input. -/
def fockVector {X : Type} (a : X → ℂ) : Hilbert X :=
  if h : Memℓp a 2 then ⟨a, h⟩ else 0

/-- A single-mode vacuum-environment splitter, PDF p.1, (1). -/
structure Splitter where
  reflectance : ℝ
  phase : ℝ
  nonnegative : 0 ≤ reflectance
  atMostOne : reflectance ≤ 1

/-- Nonnegative field reflection magnitude. -/
def Splitter.r (a : Splitter) : ℝ := Real.sqrt a.reflectance

/-- Nonnegative field transmission magnitude. -/
def Splitter.t (a : Splitter) : ℝ := Real.sqrt (1 - a.reflectance)

/-- Complex returned field coefficient, PDF p.1. -/
def Splitter.field (a : Splitter) : ℂ :=
  (a.r : ℂ) * Complex.exp (Complex.I * (a.phase : ℂ))

/-- Two candidate splitters; the prior is built into the error definition. -/
abbrev Channels := Bool → Splitter

/-- The paper's ordering and distinctness assumptions, PDF pp.1–2. -/
def ValidSplitters (c : Channels) : Prop :=
  (c false).reflectance ≤ (c true).reflectance ∧
    (c false).field ≠ (c true).field

/-- Section IV's common phase removed by the same known unitary. -/
def Reading (c : Channels) : Prop :=
  (c false).phase = 0 ∧ (c true).phase = 0

/-- The common real overlap coefficient of (17),(49). -/
def splitterOverlap (c : Channels) : ℝ :=
  (c false).r * (c true).r + (c false).t * (c true).t

/-- Total signal photon number, with no idler energy charged, PDF p.4, (18). -/
def photons {M : ℕ} (n : Occupation (Fin M)) : ℕ := ∑ m, n m

/-- Vacuum-environment binomial amplitude (12). The phase multiplies the entire
input occupation, not just the returned photon count. The guard precedes subtraction. -/
def lossAmplitude {M : ℕ} (a : Splitter) (n k : Occupation (Fin M)) : ℂ :=
  if ∀ m, k m ≤ n m then
    ∏ m, (Real.sqrt (Nat.choose (n m) (k m) : ℝ) : ℂ) *
      Complex.exp (Complex.I * (n m : ℂ) * (a.phase : ℂ)) *
      (a.r : ℂ) ^ (n m - k m) * (a.t : ℂ) ^ k m
  else 0

/-- Loss-sector coefficients (11),(25); idlers are retained without change. -/
def sectorCoefficients {I : Type} {M : ℕ}
    (a : Splitter) (psi : Hilbert (I × Occupation (Fin M)))
    (k : Occupation (Fin M)) (x : I × Occupation (Fin M)) : ℂ :=
  psi (x.1, x.2 + k) * lossAmplitude a (x.2 + k) k

/-- The unnormalized retained state in a loss sector, PDF p.4, (25). -/
def lossSector {I : Type} {M : ℕ} (a : Splitter)
    (psi : Hilbert (I × Occupation (Fin M))) (k : Occupation (Fin M)) :
    Hilbert (I × Occupation (Fin M)) := fockVector (sectorCoefficients a psi k)

/-- Actual returned density after tracing inaccessible independent vacuum modes,
PDF pp.1,3–5, (2),(11),(32). Every series is over full Fock labels. -/
def returned {I : Type} {M : ℕ} (a : Splitter)
    (psi : Hilbert (I × Occupation (Fin M))) : opticalKernel (I × Occupation (Fin M)) :=
  fun x y => ∑' k : Occupation (Fin M),
    sectorCoefficients a psi k x * star (sectorCoefficients a psi k y)

/-- Reduced signal number law, PDF p.4, before (17). -/
def signalLaw {I : Type} {M : ℕ} (psi : Hilbert (I × Occupation (Fin M)))
    (n : Occupation (Fin M)) : ℝ := ∑' i : I, Complex.normSq (psi (i, n))

/-- Total signal photon law (18), distinct from the multimode occupation law. -/
def totalLaw {I : Type} {M : ℕ} (psi : Hilbert (I × Occupation (Fin M)))
    (n : ℕ) : ℝ :=
  ∑' q : Occupation (Fin M), if photons q = n then signalLaw psi q else 0

/-- Finite mean SIGNAL energy, PDF pp.6–7. Summability rules out totalized tsum. -/
def HasSignalEnergy {I : Type} {M : ℕ} (psi : Hilbert (I × Occupation (Fin M))) (N : ℝ) : Prop :=
  0 ≤ N ∧ Summable (fun n : ℕ => (n : ℝ) * totalLaw psi n) ∧
    (∑' n : ℕ, (n : ℝ) * totalLaw psi n) = N

/-- Source definition of NDS: the signal marginal is number diagonal, PDF p.4.
For pure probes this is equivalent to orthogonal occupied idler Schmidt vectors. -/
def IsNDS {I : Type} {M : ℕ} (psi : Hilbert (I × Occupation (Fin M))) : Prop :=
  ∀ n q : Occupation (Fin M), n ≠ q →
    (∑' i : I, psi (i, n) * star (psi (i, q))) = 0

/-- Reduced density of a normalized countable-ancilla purification. -/
def purificationMarginal {X : Type} (z : Hilbert (X × ℕ)) : opticalKernel X :=
  fun x y => ∑' e : ℕ, z (x, e) * star (z (y, e))

/-- Purification as actual norm-one amplitudes, not an assumed fidelity value. -/
def PurifiesOpticalKernel {X : Type} (rho : opticalKernel X) (z : Hilbert (X × ℕ)) : Prop :=
  ‖z‖ = 1 ∧ purificationMarginal z = rho

/-- Squared Uhlmann fidelity (4),(19), as the supremum over all purifications.
For physical densities here this set is nonempty and bounded in [0,1]; the
countable ancilla suffices because the optical Hilbert spaces are separable. -/
def squaredOpticalFidelity {X : Type} (rho sigma : opticalKernel X) : ℝ :=
  sSup {v : ℝ | ∃ z w : Hilbert (X × ℕ),
    PurifiesOpticalKernel rho z ∧ PurifiesOpticalKernel sigma w ∧
      v = Complex.normSq (inner ℂ z w)}

/-- A spectral resolution of a normalized positive trace-class density. Only
positive-weight vectors are required to be orthonormal. The kernel equality is
part of a spectral representation, not an assumption on any source theorem. -/
structure Resolution {X : Type} (rho : opticalKernel X) where
  weight : ℕ → ℝ
  ket : ℕ → Hilbert X
  nonnegative : ∀ n, 0 ≤ weight n
  summable : Summable weight
  mass : ∑' n, weight n = 1
  normalized : ∀ n, 0 < weight n → ‖ket n‖ = 1
  orthogonal : ∀ n m, n ≠ m → 0 < weight n → 0 < weight m →
    inner ℂ (ket n) (ket m) = 0
  represents : ∀ x y, rho x y =
    ∑' n : ℕ, (weight n : ℂ) * ket n x * star (ket n y)

/-- A binary measurement effect is a bounded complex-linear positive contraction.
E is the decision-zero effect; 1-E is the decision-one effect, PDF pp.3,5. -/
def IsOpticalEffect {X : Type} (E : Hilbert X →L[ℂ] Hilbert X) : Prop :=
  (∀ u v, inner ℂ u (E v) = star (inner ℂ v (E u))) ∧
    ∀ v, 0 ≤ (inner ℂ v (E v)).re ∧ (inner ℂ v (E v)).re ≤ ‖v‖ ^ 2

/-- Binary Bayesian error for actual spectral states and an effect; arbitrary
joint signal-idler effects are admitted, including non-projective measurements. -/
def opticalBayesRisk {X : Type} {rho sigma : opticalKernel X}
    (A : Resolution rho) (B : Resolution sigma)
    (E : Hilbert X →L[ℂ] Hilbert X) (pi0 pi1 : ℝ) : ℝ :=
  pi0 * (∑' n : ℕ, A.weight n * (1 - (inner ℂ (A.ket n) (E (A.ket n))).re)) +
  pi1 * (∑' n : ℕ, B.weight n * (inner ℂ (B.ket n) (E (B.ket n))).re)

/-- Operational minimum Bayesian error. For normalized physical densities and
priors this is a nonempty infimum of numbers in [0,1]; spectral representations
do not constrain the allowed receiver. Equal priors are used except in (42). -/
def bayesError {X : Type} (rho sigma : opticalKernel X) (pi0 pi1 : ℝ) : ℝ :=
  sInf {v : ℝ | ∃ (A : Resolution rho) (B : Resolution sigma)
    (E : Hilbert X →L[ℂ] Hilbert X), IsOpticalEffect E ∧ v = opticalBayesRisk A B E pi0 pi1}

/-- Equal-prior minimum average error, PDF p.3, (5). -/
def opticalError {X : Type} (rho sigma : opticalKernel X) : ℝ := bayesError rho sigma (1/2) (1/2)

/-- Support spectral power: the zero eigenspace contributes zero even at t=0.
This implements the support-projection convention on PDF p.6. -/
def supportWeightPower (p t : ℝ) : ℝ := if p = 0 then 0 else Real.rpow p t

/-- Tr(rho^s sigma^(1-s)), defined independently from NDS formulas (45).
For densities and 0<=s<=1 every spectral representation gives the same absolutely
convergent nonnegative double series, in [0,1]. Hence this supremum is nonempty
and bounded; it does not optimize the probe or the receiver. -/
def chernoffAt {X : Type} (rho sigma : opticalKernel X) (s : ℝ) : ℝ :=
  sSup {v : ℝ | ∃ (A : Resolution rho) (B : Resolution sigma),
    v = ∑' nm : ℕ × ℕ,
      supportWeightPower (A.weight nm.1) s * supportWeightPower (B.weight nm.2) (1-s) *
        Complex.normSq (inner ℂ (A.ket nm.1) (B.ket nm.2))}

/-- Chernoff coefficient WITHOUT prior factor, including both support endpoints (9). -/
def opticalChernoff {X : Type} (rho sigma : opticalKernel X) : ℝ :=
  sInf {v : ℝ | ∃ s ∈ Set.Icc (0 : ℝ) 1, v = chernoffAt rho sigma s}

/-- Pure-input output fidelity, for the actual channels. -/
def outputFidelity {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) : ℝ :=
  squaredOpticalFidelity (returned (c false) psi) (returned (c true) psi)

/-- Pure-input operational output error. -/
def outputError {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) : ℝ :=
  opticalError (returned (c false) psi) (returned (c true) psi)

/-- Pure-input Chernoff-type coefficient. -/
def outputChernoffAt {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (s : ℝ) : ℝ :=
  chernoffAt (returned (c false) psi) (returned (c true) psi) s

/-- Pure-input optimized Chernoff coefficient. -/
def outputChernoff {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) : ℝ :=
  opticalChernoff (returned (c false) psi) (returned (c true) psi)

/-- Complex overlap series in the TOTAL photon law, PDF p.4, (17). -/
def purificationOverlapTerm {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (n : ℕ) : ℂ :=
  (totalLaw psi n : ℂ) *
    Complex.exp (Complex.I * (n : ℂ) * ((c true).phase - (c false).phase : ℝ)) *
    (splitterOverlap c : ℂ) ^ n

/-- Squared overlap (17); the absolute value encloses the WHOLE complex sum. -/
def purificationOverlap {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) : ℝ :=
  Complex.normSq (∑' n : ℕ, purificationOverlapTerm c psi n)

/-- Fidelity-based LOWER error envelope (6),(20),(53). -/
def lowerError (F : ℝ) : ℝ := (1 - Real.sqrt (1-F)) / 2

/-- Loss-sector probability (27) directly from its actual vector. -/
def lossProbability {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (b : Bool) (k : Occupation (Fin M)) : ℝ :=
  ‖lossSector (c b) psi k‖ ^ 2

/-- Loss-sector overlap (30); the first argument is conjugated. -/
def lossOverlap {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (k : Occupation (Fin M)) : ℂ :=
  inner ℂ (lossSector (c false) psi k) (lossSector (c true) psi k)

/-- Explicit complex binomial summand (31), with componentwise loss guard. -/
def lossOverlapTerm {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (k n : Occupation (Fin M)) : ℂ :=
  (signalLaw psi n : ℂ) * star (lossAmplitude (c false) n k) * lossAmplitude (c true) n k

/-- The at-most-two-dimensional block H_k of (33)–(35). -/
def lossBlock {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (k : Occupation (Fin M)) :
    Submodule ℂ (Hilbert (I × Occupation (Fin M))) :=
  Submodule.span ℂ {lossSector (c false) psi k, lossSector (c true) psi k}

/-- Weighted block difference; equal priors use weights 1,1 for rho0-rho1. -/
def lossBlockDifference {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (k : Occupation (Fin M)) (w0 w1 : ℝ)
    (v : Hilbert (I × Occupation (Fin M))) : Hilbert (I × Occupation (Fin M)) :=
  ((w0 : ℂ) * inner ℂ (lossSector (c false) psi k) v) • lossSector (c false) psi k -
  ((w1 : ℂ) * inner ℂ (lossSector (c true) psi k) v) • lossSector (c true) psi k

/-- A complete block spectral receiver, (40)–(42). E is the zero decision.
It is projective, preserves every loss block, and decides the positive/negative
spectra correctly. Zero eigenvectors may be assigned either decision. -/
def BlockReceiver {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (w0 w1 : ℝ)
    (E : Hilbert (I × Occupation (Fin M)) →L[ℂ] Hilbert (I × Occupation (Fin M))) : Prop :=
  IsOpticalEffect E ∧ E.comp E = E ∧
    (∀ k v, v ∈ lossBlock c psi k → E v ∈ lossBlock c psi k) ∧
    (∀ k v (a : ℝ), v ∈ lossBlock c psi k →
      lossBlockDifference c psi k w0 w1 v = (a : ℂ) • v →
        (0 < a → E v = v) ∧ (a < 0 → E v = 0))

/-- Risk of a concrete effect, directly on unnormalized output loss vectors. -/
def sectorRisk {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M)))
    (E : Hilbert (I × Occupation (Fin M)) →L[ℂ] Hilbert (I × Occupation (Fin M)))
    (pi0 pi1 : ℝ) : ℝ :=
  pi0 * (∑' k : Occupation (Fin M), (lossProbability c psi false k -
    (inner ℂ (lossSector (c false) psi k) (E (lossSector (c false) psi k))).re)) +
  pi1 * (∑' k, (inner ℂ (lossSector (c true) psi k) (E (lossSector (c true) psi k))).re)

/-- Term of (45), explicitly omitting either zero-probability sector. -/
def ndsChernoffTerm {I : Type} {M : ℕ} (c : Channels)
    (psi : Hilbert (I × Occupation (Fin M))) (s : ℝ) (k : Occupation (Fin M)) : ℝ :=
  if lossProbability c psi false k = 0 ∨ lossProbability c psi true k = 0 then 0
  else Real.rpow (lossProbability c psi false k) (s-1) *
    Real.rpow (lossProbability c psi true k) (-s) * Complex.normSq (lossOverlap c psi k)

/-- Rank-one full-Fock density. -/
def opticalPureKernel {X : Type} (psi : Hilbert X) : opticalKernel X :=
  fun x y => psi x * star (psi y)

/-- A Fock signal with no retained idlers; the empty idler occupation has one value. -/
def fock {M : ℕ} (N : Occupation (Fin M)) : Hilbert (Occupation (Fin 0) × Occupation (Fin M)) :=
  lp.single 2 (0, N) 1

/-- Independent product probes with arbitrary finite idler count J_m per signal,
PDF p.6, (46)–(47). The idler register is grouped, without tracing any factors. -/
def productOpticalProbe {M : ℕ} (J : Fin M → ℕ)
    (psi : ∀ m, Hilbert (Occupation (Fin (J m)) × Occupation (Fin 1))) :
    Hilbert ((∀ m, Occupation (Fin (J m))) × Occupation (Fin M)) :=
  fockVector (fun x => ∏ m, psi m (x.1 m, fun _ => x.2 m))

/-- Independent copies of a joint multimode probe, with all signals grouped. -/
def opticalCopies {J M : ℕ} (psi : Hilbert (Occupation (Fin J) × Occupation (Fin M))) (L : ℕ) :
    Hilbert ((Fin L → Occupation (Fin J)) × Occupation (Fin (L*M))) :=
  fockVector (fun x => ∏ l : Fin L,
    psi (x.1 l, fun m : Fin M => x.2 ⟨l.val*M+m.val, by
      calc
        l.val*M+m.val < l.val*M+M := Nat.add_lt_add_left m.isLt _
        _ = (l.val+1)*M := (Nat.succ_mul _ _).symm
        _ ≤ L*M := Nat.mul_le_mul_right M (Nat.succ_le_of_lt l.isLt)⟩))

/-- Product coherent probe; arbitrary classical idler amplitudes are retained. -/
def coherentOpticalProbe {J M : ℕ} (beta : Fin J → ℂ) (alpha : Fin M → ℂ) :
    Hilbert (Occupation (Fin J) × Occupation (Fin M)) :=
  fockVector (fun x => (∏ j, coherentCoefficient (beta j) (x.1 j)) *
    ∏ m, coherentCoefficient (alpha m) (x.2 m))

/-- Total signal energy of a product coherent state. -/
def coherentEnergy {M : ℕ} (alpha : Fin M → ℂ) : ℝ :=
  ∑ m, Complex.normSq (alpha m)

/-- Multimode EPR coefficients (73), with N the per-pair SIGNAL energy.
Natural powers preserve the vacuum convention 0^0=1. -/
def tmsvCoefficients (M : ℕ) (N : ℝ) (x : Occupation (Fin M) × Occupation (Fin M)) : ℂ :=
  if x.1 = x.2 then
    ∏ m, ((Real.sqrt (1/(N+1)) : ℝ) : ℂ) *
      ((Real.sqrt (N/(N+1)) : ℝ) : ℂ) ^ x.2 m
  else 0

/-- Product of M identical EPR pairs, with no photon cutoff. -/
def tmsvProbe (M : ℕ) (N : ℝ) : Hilbert (Occupation (Fin M) × Occupation (Fin M)) :=
  fockVector (tmsvCoefficients M N)

/-- Product law of observed RETURNED counts, PDF p.8, (62) with j=N-k. -/
def countLaw {M : ℕ} (a : Splitter) (N j : Occupation (Fin M)) : ℝ :=
  if ∀ m, j m ≤ N m then
    ∏ m, (Nat.choose (N m) (j m) : ℝ) *
      a.reflectance ^ j m * (1-a.reflectance) ^ (N m-j m)
  else 0

/-- Binomial law of the total returned photon count. -/
def totalCountLaw (a : Splitter) (K j : ℕ) : ℝ :=
  if j ≤ K then (Nat.choose K j : ℝ) * a.reflectance ^ j *
    (1-a.reflectance) ^ (K-j) else 0

/-- Equal-prior overlap of the total count distributions, not an assumed quantum optimum. -/
def countingError (c : Channels) (K : ℕ) : ℝ :=
  (∑ j ∈ Finset.range (K+1),
    min (totalCountLaw (c false) K j) (totalCountLaw (c true) K j)) / 2

/-- Photon counting followed by d(j), where true means decision one.
E is the zero-decision effect and acts trivially on every retained idler. -/
def CountsAs {I : Type} {M : ℕ}
    (E : Hilbert (I × Occupation (Fin M)) →L[ℂ] Hilbert (I × Occupation (Fin M)))
    (d : Occupation (Fin M) → Bool) : Prop :=
  ∀ v x, (E v) x = if d x.2 then 0 else v x

/-- Interior logarithmic threshold (65); endpoint rules are separate theorems. -/
def countThreshold (c : Channels) (K : ℕ) : ℝ :=
  (K : ℝ) * Real.log ((1-(c false).reflectance)/(1-(c true).reflectance)) /
    Real.log ((c true).reflectance*(1-(c false).reflectance) /
      ((c false).reflectance*(1-(c true).reflectance)))

/-- One-photon classical Chernoff coefficient, PDF p.8, (66). -/
def fockCoefficient (c : Channels) (s : ℝ) : ℝ :=
  Real.rpow (c false).reflectance s * Real.rpow (c true).reflectance (1-s) +
  Real.rpow (1-(c false).reflectance) s * Real.rpow (1-(c true).reflectance) (1-s)

/-- Explicit optimizer (67); theorem hypotheses restrict it to strict interior channels. -/
def fockOptimizer (c : Channels) : ℝ :=
  Real.log (-(c true).reflectance * Real.log ((c false).reflectance/(c true).reflectance) /
    ((1-(c true).reflectance) * Real.log ((1-(c false).reflectance)/(1-(c true).reflectance)))) /
  Real.log ((1-(c false).reflectance)*(c true).reflectance /
    ((1-(c true).reflectance)*(c false).reflectance))

/-- Per-photon number-state exponent, with natural logarithm, (66). -/
def numberExponent (c : Channels) : ℝ := -Real.log (fockCoefficient c (fockOptimizer c))

/-- Per-photon coherent exponent (68). -/
def coherentExponent (c : Channels) : ℝ := ((c true).r-(c false).r)^2

/-- Exponent gain (69); no claim that it is large in every regime. -/
def exponentGain (c : Channels) : ℝ := numberExponent c / coherentExponent c

/-- Parameters (78)–(81); beta is used only when R1<1. -/
def eprAlpha (c : Channels) (N : ℝ) : ℝ :=
  ((1-(c false).reflectance)*N+1) / ((1-(c true).reflectance)*N+1)

/-- Parameter beta in (79). -/
def eprBeta (c : Channels) : ℝ := (1-(c false).reflectance)/(1-(c true).reflectance)

/-- Parameter C in (80). -/
def eprC (c : Channels) (N : ℝ) : ℝ :=
  ((1-(c false).r*(c true).r)*N+1)^2 / ((1-(c false).reflectance)*N+1)

/-- Parameter D in (81). -/
def eprD (c : Channels) (N : ℝ) : ℝ := (1-(c true).reflectance)*N

/-- Equivalent logarithmic case test (84)–(85): no exponent/product ambiguity. -/
def eprDerivativeSign (c : Channels) (N : ℝ) : ℝ :=
  eprD c N * Real.log (eprBeta c) - eprC c N * Real.log (eprAlpha c N)

/-- Optimizer (85)–(86). N=0 and ideal memory are handled outside this definition. -/
def eprOptimizer (c : Channels) (N : ℝ) : ℝ :=
  if 0 ≤ eprDerivativeSign c N then 0 else
    min 1 (Real.log (eprC c N * Real.log (eprAlpha c N) /
      (eprD c N * Real.log (eprBeta c))) /
        Real.log (eprBeta c / eprAlpha c N))

/-- Classical positive-P probes include arbitrary probability measures and Dirac
masses. Signal and idler amplitudes may be classically correlated, (3), pp.2,8. -/
def classicalReturned {J M : ℕ} (a : Splitter)
    (nu : Measure ((Fin J → ℂ) × (Fin M → ℂ))) :
    opticalKernel (Occupation (Fin J) × Occupation (Fin M)) :=
  fun x y => ∫ alpha,
    opticalPureKernel (coherentOpticalProbe alpha.1 (fun m => a.field * alpha.2 m)) x y ∂nu

/-- An additional uniform independent vacuum attenuator, PDF p.12, IV.E.2. -/
def vacuumEfficiency (eta : ℝ) (h0 : 0 ≤ eta) (h1 : eta ≤ 1) : Splitter :=
  ⟨eta, 0, h0, h1⟩

/-- A uniform receiver efficiency multiplies energy reflectance. Eta=0 is allowed
and makes the two effective returned channels identical. -/
def efficiencyScaled (a : Splitter) (eta : ℝ) (h0 : 0 ≤ eta) (h1 : eta ≤ 1) : Splitter :=
  ⟨eta*a.reflectance, a.phase, mul_nonneg h0 a.nonnegative,
    by
      calc
        eta*a.reflectance ≤ 1*a.reflectance := mul_le_mul_of_nonneg_right h1 a.nonnegative
        _ = a.reflectance := one_mul _
        _ ≤ 1 := a.atMostOne⟩

/-- Vacuum attenuator on a full joint operator kernel, retaining every idler. -/
def attenuate {I : Type} {M : ℕ} (a : Splitter)
    (rho : opticalKernel (I × Occupation (Fin M))) : opticalKernel (I × Occupation (Fin M)) :=
  fun x y => ∑' k : Occupation (Fin M),
    lossAmplitude a (x.2+k) k * star (lossAmplitude a (y.2+k) k) *
      rho (x.1, x.2+k) (y.1, y.2+k)

/-! Small sanity tests pin definitions without proving source results. -/

example (p : ℝ) : supportWeightPower 0 p = 0 := by simp [supportWeightPower]

example {M : ℕ} (N : Occupation (Fin M)) : fock N (0, N) = 1 := by
  simp [fock]

example {M : ℕ} (N : Occupation (Fin M)) : ‖fock N‖ = 1 := by
  simp [fock, lp.norm_single (by norm_num : (0 : ENNReal) < 2)]

example (a : Splitter) : totalCountLaw a 0 0 = 1 := by
  simp [totalCountLaw]

example (s : ℝ) : lowerError 1 = 1/2 := by norm_num [lowerError]

example (eta : ℝ) (h0 : 0 ≤ eta) (h1 : eta ≤ 1) (a : Splitter) :
    (efficiencyScaled a eta h0 h1).reflectance = eta*a.reflectance := rfl

example {X : Type} (rho : opticalKernel X) : ¬ PurifiesOpticalKernel rho (0 : Hilbert (X × ℕ)) := by
  simp [PurifiesOpticalKernel]

example {X : Type} : IsOpticalEffect (0 : Hilbert X →L[ℂ] Hilbert X) := by
  simp [IsOpticalEffect, sq_nonneg]

example (c : Channels) : countingError c 0 = 1/2 := by
  simp [countingError, totalCountLaw]

example (M : ℕ) : tmsvCoefficients M 0 (0, 0) = 1 := by
  simp [tmsvCoefficients]

example {M : ℕ} (a : Splitter) : lossAmplitude a (0 : Occupation (Fin M)) 0 = 1 := by
  simp [lossAmplitude]

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

/-- Dual trace norm of a weighted Hermitian difference, used in (42).
The supremum ranges over ALL bounded self-adjoint contractions. Spectral
resolutions merely evaluate traces of the actual two states. -/
def weightedTraceNorm {X : Type} (rho sigma : opticalKernel X) (w0 w1 : ℝ) : ℝ :=
  sSup {v : ℝ | ∃ (A : Resolution rho) (B : Resolution sigma)
    (T : Hilbert X →L[ℂ] Hilbert X),
    (∀ u v, inner ℂ u (T v) = star (inner ℂ v (T u))) ∧
    (∀ u, |(inner ℂ u (T u)).re| ≤ ‖u‖^2) ∧
    v = w0 * (∑' n, A.weight n * (inner ℂ (A.ket n) (T (A.ket n))).re) -
      w1 * (∑' n, B.weight n * (inner ℂ (B.ket n) (T (B.ket n))).re)}

/-- Composition of full Fock operator kernels on the physical trace-class arguments. -/
def kernelProduct {X : Type} (rho sigma : opticalKernel X) : opticalKernel X :=
  fun x y => ∑' z : X, rho x z * sigma z y

/-- Target-detection universal envelope (59), with fixed real signal energy N. -/
def targetLower (R N : ℝ) : ℝ := lowerError (Real.rpow (1-R) N)

/-- Coherent target-detection error envelope (60). -/
def targetCoherent (R N : ℝ) : ℝ := lowerError (Real.exp (-R*N))

end QAlgorithms.Papers.Independent.Intakethird100.P54
Source
Discriminating quantum-optical beam-splitter channels with number-diagonal signal states: Applications to quantum reading and target detection, arXiv:1105.4063v3, https://arxiv.org/abs/1105.4063v3
Read-back

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

QAlgorithms.Papers.Independent.Intakethird100.P54.BlockReceiver

For any type III, natural number M≥0M\ge0M≥0, splitter pair ccc, vector ψ∈H=ℓ2(I×NM;C)\psi\in H=\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈H=ℓ2(I×NM​;C), real weights w0,w1w_0,w_1w0​,w1​, and bounded complex-linear operator E:H→HE:H\to HE:H→H, where NM={0,…,M−1}→N\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb NNM​={0,…,M−1}→N, define Lb,kL_{b,k}Lb,k​, for b=0,1b=0,1b=0,1 and k∈NMk\in\mathcal N_Mk∈NM​, by the coordinates

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m},ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,

where Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R are splitter c(b)c(b)c(b)’s parameters, and use the zero vector if these coordinates are not square-summable. Put Vk=span⁡C{L0,k,L1,k}V_k=\operatorname{span}_{\mathbb C}\{L_{0,k},L_{1,k}\}Vk​=spanC​{L0,k​,L1,k​} and Dkv=w0⟨L0,k,v⟩L0,k−w1⟨L1,k,v⟩L1,kD_kv=w_0\langle L_{0,k},v\rangle L_{0,k}-w_1\langle L_{1,k},v\rangle L_{1,k}Dk​v=w0​⟨L0,k​,v⟩L0,k​−w1​⟨L1,k​,v⟩L1,k​. The defined proposition requires all of the following: for every u,v∈Hu,v\in Hu,v∈H, ⟨u,Ev⟩=⟨v,Eu⟩‾\langle u,Ev\rangle=\overline{\langle v,Eu\rangle}⟨u,Ev⟩=⟨v,Eu⟩​; for every v∈Hv\in Hv∈H, 0≤Re⁡⟨v,Ev⟩≤∥v∥20\le\operatorname{Re}\langle v,Ev\rangle\le\|v\|^20≤Re⟨v,Ev⟩≤∥v∥2; E2=EE^2=EE2=E; every VkV_kVk​ is invariant under EEE; and for every kkk, every v∈Vkv\in V_kv∈Vk​, and every real aaa satisfying Dkv=avD_kv=avDk​v=av, one has Ev=vEv=vEv=v if a>0a>0a>0 and Ev=0Ev=0Ev=0 if a<0a<0a<0. The last condition imposes no further requirement when a=0a=0a=0. The weights are unrestricted and ψ\psiψ need not be normalized. The case M=0M=0M=0 is included, with empty product 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.Channels

A channel pair is a function from the two Boolean labels to splitters. Writing its values as c0=c(false)c_0=c(\mathrm{false})c0​=c(false) and c1=c(true)c_1=c(\mathrm{true})c1​=c(true), each cbc_bcb​ has a reflectance Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and a real phase θb\theta_bθb​. The two splitters may coincide.

QAlgorithms.Papers.Independent.Intakethird100.P54.CountsAs

For any type III, any natural number MMM, any bounded complex-linear operator EEE on ℓ2(I×NM;C)\ell^2(I\times\mathbb N^M;\mathbb C)ℓ2(I×NM;C), and any Boolean-valued function d:NM→{false,true}d:\mathbb N^M\to\{\mathrm{false},\mathrm{true}\}d:NM→{false,true}, this is the proposition that, for every vector vvv in that Hilbert space and every coordinate (i,n)∈I×NM(i,n)\in I\times\mathbb N^M(i,n)∈I×NM,

(Ev)(i,n)={0,d(n)=true,v(i,n),d(n)=false.(Ev)(i,n)= \begin{cases} 0,&d(n)=\mathrm{true},\\ v(i,n),&d(n)=\mathrm{false}. \end{cases}(Ev)(i,n)={0,v(i,n),​d(n)=true,d(n)=false.​

There is no finiteness, countability, or nonemptiness assumption on III. If III is empty, the coordinate condition is vacuous. If M=0M=0M=0, the occupation index set is a singleton, so ddd has a single Boolean value.

QAlgorithms.Papers.Independent.Intakethird100.P54.HasSignalEnergy

For any type III, natural number MMM, square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), and real number NNN, define

Pψ(n)=∑q∈N{0,…,M−1}{∑i∈I∣ψ(i,q)∣2,∑m=0M−1qm=n,0,otherwise.P_\psi(n)= \sum_{q\in\mathbb N^{\{0,\ldots,M-1\}}} \begin{cases} \displaystyle\sum_{i\in I}|\psi(i,q)|^2, &\displaystyle\sum_{m=0}^{M-1}q_m=n,\\ 0,&\text{otherwise}. \end{cases}Pψ​(n)=q∈N{0,…,M−1}∑​⎩⎨⎧​i∈I∑​∣ψ(i,q)∣2,0,​m=0∑M−1​qm​=n,otherwise.​

The stated property means that N≥0N\ge0N≥0, the real series ∑n=0∞nPψ(n)\sum_{n=0}^{\infty}nP_\psi(n)∑n=0∞​nPψ​(n) is summable, and its sum equals NNN. No normalization of ψ\psiψ is required. The n=0n=0n=0 term contributes 000; if M=0M=0M=0, if III is empty, or if ψ=0\psi=0ψ=0, the property holds exactly when N=0N=0N=0.

QAlgorithms.Papers.Independent.Intakethird100.P54.IsNDS

For any type III, natural number MMM, and square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), this property requires that every pair of distinct occupation vectors n,q∈N{0,…,M−1}n,q\in\mathbb N^{\{0,\ldots,M-1\}}n,q∈N{0,…,M−1} satisfy

∑i∈Iψ(i,n) ψ(i,q)‾=0.\sum_{i\in I}\psi(i,n)\,\overline{\psi(i,q)}=0.i∈I∑​ψ(i,n)ψ(i,q)​=0.

There is no requirement for pairs with n=qn=qn=q, and no normalization or nonzero condition on ψ\psiψ. If M=0M=0M=0, there are no distinct occupation vectors, so the property holds vacuously. If III is empty or ψ=0\psi=0ψ=0, every displayed sum is 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.IsOpticalEffect

For any type XXX and continuous complex-linear operator E:ℓ2(X)→ℓ2(X)E:\ell^2(X)\to\ell^2(X)E:ℓ2(X)→ℓ2(X), this property means that

⟨u,Ev⟩=⟨v,Eu⟩‾for every u,v∈ℓ2(X),\langle u,Ev\rangle=\overline{\langle v,Eu\rangle} \qquad\text{for every }u,v\in\ell^2(X),⟨u,Ev⟩=⟨v,Eu⟩​for every u,v∈ℓ2(X),

and that

0≤Re⁡⟨v,Ev⟩≤∥v∥2for every v∈ℓ2(X).0\le\operatorname{Re}\langle v,Ev\rangle\le\|v\|^2 \qquad\text{for every }v\in\ell^2(X).0≤Re⟨v,Ev⟩≤∥v∥2for every v∈ℓ2(X).

The vectors range over the entire Hilbert space, including 000, without normalization. No idempotence requirement is included. If XXX is empty, the unique operator on the zero Hilbert space satisfies the property.

QAlgorithms.Papers.Independent.Intakethird100.P54.PurifiesOpticalKernel

For any type XXX, arbitrary complex-valued kernel ρ:X×X→C\rho:X\times X\to\mathbb Cρ:X×X→C, and square-summable complex family z∈ℓ2(X×N)z\in\ell^2(X\times\mathbb N)z∈ℓ2(X×N), this property means both that ∥z∥=1\|z\|=1∥z∥=1 and that

ρ(x,y)=∑e=0∞z(x,e) z(y,e)‾for every x,y∈X.\rho(x,y)=\sum_{e=0}^{\infty}z(x,e)\,\overline{z(y,e)} \qquad\text{for every }x,y\in X.ρ(x,y)=e=0∑∞​z(x,e)z(y,e)​for every x,y∈X.

There are no separate assumptions on ρ\rhoρ; the displayed requirements are the entire property. If XXX is empty, its Hilbert space is zero-dimensional and no such unit vector exists, so the property cannot hold.

QAlgorithms.Papers.Independent.Intakethird100.P54.Reading

For any channel pair ccc, consisting of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and real phases θ0,θ1\theta_0,\theta_1θ0​,θ1​, this property means exactly that θ0=0\theta_0=0θ0​=0 and θ1=0\theta_1=0θ1​=0, as equalities of real numbers. It imposes no ordering or distinctness condition on the reflectances.

QAlgorithms.Papers.Independent.Intakethird100.P54.Resolution

For any type XXX and arbitrary complex-valued kernel ρ:X×X→C\rho:X\times X\to\mathbb Cρ:X×X→C, a resolution consists of a real sequence (pn)n∈N(p_n)_{n\in\mathbb N}(pn​)n∈N​ and vectors un∈ℓ2(X)u_n\in\ell^2(X)un​∈ℓ2(X) such that pn≥0p_n\ge0pn​≥0 for every nnn, the sequence (pn)(p_n)(pn​) is summable with ∑n=0∞pn=1\sum_{n=0}^{\infty}p_n=1∑n=0∞​pn​=1, every unu_nun​ with pn>0p_n>0pn​>0 has norm 111, every distinct pair n,mn,mn,m with pn>0p_n>0pn​>0 and pm>0p_m>0pm​>0 satisfies ⟨un,um⟩=0\langle u_n,u_m\rangle=0⟨un​,um​⟩=0, and

ρ(x,y)=∑n=0∞pn un(x) un(y)‾for every x,y∈X.\rho(x,y)=\sum_{n=0}^{\infty}p_n\,u_n(x)\,\overline{u_n(y)} \qquad\text{for every }x,y\in X.ρ(x,y)=n=0∑∞​pn​un​(x)un​(y)​for every x,y∈X.

Vectors at indices of zero weight have no normalization or orthogonality requirement. The sequence may have finitely or infinitely many positive weights. Existence of a resolution is not asserted for an arbitrary kernel; if XXX is empty, no resolution exists because total weight 111 forces a positive weight but the zero Hilbert space has no unit vector.

QAlgorithms.Papers.Independent.Intakethird100.P54.Splitter.field

For any splitter aaa, with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and real phase θ\thetaθ, this defines the complex number R eiθ\sqrt R\,e^{i\theta}R​eiθ. When R=0R=0R=0, this number is zero for every phase.

QAlgorithms.Papers.Independent.Intakethird100.P54.Splitter

A splitter consists of two real numbers, a reflectance RRR and a phase θ\thetaθ, together with proofs that 0≤R0\le R0≤R and R≤1R\le1R≤1. Both endpoints R=0R=0R=0 and R=1R=1R=1 are allowed, and θ\thetaθ is unrestricted.

QAlgorithms.Papers.Independent.Intakethird100.P54.Splitter.r

For any splitter aaa, with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and real phase θ\thetaθ, this defines r(a)=Rr(a)=\sqrt Rr(a)=R​, using the nonnegative real square root; in particular, r(a)=0r(a)=0r(a)=0 when R=0R=0R=0.

QAlgorithms.Papers.Independent.Intakethird100.P54.Splitter.t

For any splitter aaa, with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and real phase θ\thetaθ, this defines t(a)=1−Rt(a)=\sqrt{1-R}t(a)=1−R​, using the nonnegative real square root; in particular, t(a)=0t(a)=0t(a)=0 when R=1R=1R=1.

QAlgorithms.Papers.Independent.Intakethird100.P54.ValidSplitters

For any channel pair ccc, whose splitter at label b∈{0,1}b\in\{0,1\}b∈{0,1} has reflectance Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and real phase θb\theta_bθb​, this property means exactly that R0≤R1R_0\le R_1R0​≤R1​ and R0 eiθ0≠R1 eiθ1\sqrt{R_0}\,e^{i\theta_0}\ne\sqrt{R_1}\,e^{i\theta_1}R0​​eiθ0​=R1​​eiθ1​. Equality of the reflectances is allowed provided the complex fields differ; if both reflectances are zero, the property fails.

QAlgorithms.Papers.Independent.Intakethird100.P54.attenuate

For any type III, any natural number MMM, a splitter aaa with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and phase θ∈R\theta\in\mathbb Rθ∈R, and any complex kernel ρ\rhoρ on I×N{0,…,M−1}I\times\mathbb N^{\{0,\ldots,M-1\}}I×N{0,…,M−1}, this defines another kernel as follows. Put r=Rr=\sqrt Rr=R​, t=1−Rt=\sqrt{1-R}t=1−R​, and, for occupation vectors p,kp,kp,k, put Aa(p,k)=∏m=0M−1(pm+kmkm) e i(pm+km)θrpmtkmA_a(p,k)=\prod_{m=0}^{M-1}\sqrt{\binom{p_m+k_m}{k_m}}\,e^{\,i(p_m+k_m)\theta}r^{p_m}t^{k_m}Aa​(p,k)=∏m=0M−1​(km​pm​+km​​)​ei(pm​+km​)θrpm​tkm​. Then

attenuate⁡(a,ρ)((i,p),(j,q))=∑k∈NMAa(p,k) Aa(q,k)‾ ρ((i,p+k),(j,q+k)).\operatorname{attenuate}(a,\rho)\bigl((i,p),(j,q)\bigr) =\sum_{k\in\mathbb N^M}A_a(p,k)\,\overline{A_a(q,k)}\, \rho\bigl((i,p+k),(j,q+k)\bigr).attenuate(a,ρ)((i,p),(j,q))=k∈NM∑​Aa​(p,k)Aa​(q,k)​ρ((i,p+k),(j,q+k)).

Addition of occupation vectors is componentwise. No positivity, normalization, or summability hypothesis is imposed on ρ\rhoρ; the infinite sum is the total complex summation operation, which equals zero if the summand is not summable. The definition includes M=0M=0M=0, when there is one empty occupation vector, all products are empty products equal to 111, and the resulting kernel equals ρ\rhoρ. The endpoint reflectances 000 and 111 are included, with natural powers satisfying 00=10^0=100=1.

QAlgorithms.Papers.Independent.Intakethird100.P54.bayesError

For any type XXX, arbitrary kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, and arbitrary real numbers π0,π1\pi_0,\pi_1π0​,π1​, put H=ℓ2(X;C)H=\ell^2(X;\mathbb C)H=ℓ2(X;C). A resolution (wn,un)(w_n,u_n)(wn​,un​) of a kernel KKK consists of summable nonnegative real weights with ∑nwn=1\sum_nw_n=1∑n​wn​=1, vectors un∈Hu_n\in Hun​∈H normalized whenever wn>0w_n>0wn​>0, orthogonality ⟨un,um⟩=0\langle u_n,u_m\rangle=0⟨un​,um​⟩=0 whenever n≠mn\ne mn=m and both weights are positive, and the identities K(x,y)=∑nwnun(x)un(y)‾K(x,y)=\sum_nw_nu_n(x)\overline{u_n(y)}K(x,y)=∑n​wn​un​(x)un​(y)​. The value bayesError⁡(ρ,σ,π0,π1)\operatorname{bayesError}(\rho,\sigma,\pi_0,\pi_1)bayesError(ρ,σ,π0​,π1​) is the real infimum, over all resolutions (an,un)(a_n,u_n)(an​,un​) of ρ\rhoρ, resolutions (bn,vn)(b_n,v_n)(bn​,vn​) of σ\sigmaσ, and continuous complex-linear operators E:H→HE:H\to HE:H→H satisfying

⟨u,Ev⟩=⟨v,Eu⟩‾(u,v∈H),0≤Re⁡⟨v,Ev⟩≤∥v∥2(v∈H),\langle u,Ev\rangle=\overline{\langle v,Eu\rangle}\quad(u,v\in H), \qquad 0\le\operatorname{Re}\langle v,Ev\rangle\le\|v\|^2\quad(v\in H),⟨u,Ev⟩=⟨v,Eu⟩​(u,v∈H),0≤Re⟨v,Ev⟩≤∥v∥2(v∈H),

of the numbers

π0∑n=0∞an(1−Re⁡⟨un,Eun⟩)+π1∑n=0∞bnRe⁡⟨vn,Evn⟩.\pi_0\sum_{n=0}^{\infty}a_n\bigl(1-\operatorname{Re}\langle u_n,Eu_n\rangle\bigr) +\pi_1\sum_{n=0}^{\infty}b_n\operatorname{Re}\langle v_n,Ev_n\rangle.π0​n=0∑∞​an​(1−Re⟨un​,Eun​⟩)+π1​n=0∑∞​bn​Re⟨vn​,Evn​⟩.

The priors have no sign or normalization assumptions, and the kernels are not assumed to admit resolutions. If there are no candidates, the real infimum of the empty set is 000; in particular, this applies when XXX is empty. Infinite sums use the total convention of value 000 for nonsummable families.

QAlgorithms.Papers.Independent.Intakethird100.P54.chernoffAt

For any type XXX, arbitrary kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, and any real sss, let H=ℓ2(X;C)H=\ell^2(X;\mathbb C)H=ℓ2(X;C) and define ht(p)=0h_t(p)=0ht​(p)=0 if p=0p=0p=0 and ht(p)=pth_t(p)=p^tht​(p)=pt otherwise, using real powers. A resolution (wn,un)(w_n,u_n)(wn​,un​) of a kernel KKK means nonnegative summable weights with sum 111, vectors in HHH normalized at positive weights and mutually orthogonal between distinct positive-weight indices, and K(x,y)=∑nwnun(x)un(y)‾K(x,y)=\sum_nw_nu_n(x)\overline{u_n(y)}K(x,y)=∑n​wn​un​(x)un​(y)​ for every x,yx,yx,y. Then

chernoffAt⁡(ρ,σ,s)=sup⁡(an,un) resolves ρ(bm,vm) resolves σ{∑(n,m)∈N2hs(an) h1−s(bm) ∣⟨un,vm⟩∣2},\operatorname{chernoffAt}(\rho,\sigma,s) = \sup_{\substack{(a_n,u_n)\text{ resolves }\rho\\(b_m,v_m)\text{ resolves }\sigma}} \left\{ \sum_{(n,m)\in\mathbb N^2} h_s(a_n)\,h_{1-s}(b_m)\,|\langle u_n,v_m\rangle|^2 \right\},chernoffAt(ρ,σ,s)=(an​,un​) resolves ρ(bm​,vm​) resolves σ​sup​⎩⎨⎧​(n,m)∈N2∑​hs​(an​)h1−s​(bm​)∣⟨un​,vm​⟩∣2⎭⎬⎫​,

where the supremum is taken in R\mathbb RR. There is no restriction 0≤s≤10\le s\le10≤s≤1 in this definition. A zero weight contributes 000 even when its exponent is zero or negative. The sum is indexed by pairs, and has value 000 if its family is nonsummable. Resolutions need not exist; the real supremum is 000 for an empty candidate set, including when XXX is empty, and also uses the total value 000 for a set unbounded above.

QAlgorithms.Papers.Independent.Intakethird100.P54.classicalReturned

For arbitrary natural numbers J,MJ,MJ,M, including zero, a splitter aaa with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and real phase θ\thetaθ, and any measure ν\nuν on CJ×CM\mathbb C^J\times\mathbb C^MCJ×CM, this defines a complex kernel on pairs of occupations x=(p,q)x=(p,q)x=(p,q), where p:{0,…,J−1}→Np:\{0,\ldots,J-1\}\to\mathbb Np:{0,…,J−1}→N and q:{0,…,M−1}→Nq:\{0,\ldots,M-1\}\to\mathbb Nq:{0,…,M−1}→N. Put f=R eiθf=\sqrt R\,e^{\mathrm i\theta}f=R​eiθ and

h(z,n)=e−∣z∣2/2znn!,gβ,α(p,q)=(∏j=0J−1h(βj,pj))(∏m=0M−1h(fαm,qm)).h(z,n)=e^{-|z|^2/2}\frac{z^n}{\sqrt{n!}},\qquad g_{\beta,\alpha}(p,q)= \left(\prod_{j=0}^{J-1}h(\beta_j,p_j)\right) \left(\prod_{m=0}^{M-1}h(f\alpha_m,q_m)\right).h(z,n)=e−∣z∣2/2n!​zn​,gβ,α​(p,q)=(j=0∏J−1​h(βj​,pj​))(m=0∏M−1​h(fαm​,qm​)).

Let ψβ,α\psi_{\beta,\alpha}ψβ,α​ be this coefficient function regarded as a vector in the complex space ℓ2\ell^2ℓ2 of occupation pairs when it belongs to that space, and the zero vector otherwise. The defined kernel is

K(x,y)=∫CJ×CMψβ,α(x) ψβ,α(y)‾ dν(β,α).K(x,y)=\int_{\mathbb C^J\times\mathbb C^M} \psi_{\beta,\alpha}(x)\, \overline{\psi_{\beta,\alpha}(y)}\,d\nu(\beta,\alpha).K(x,y)=∫CJ×CM​ψβ,α​(x)ψβ,α​(y)​dν(β,α).

Empty products are 111. The measure is not required to be a probability measure, and no integrability hypothesis is supplied; each integral uses the total Bochner integral, which is zero when its integrand is not integrable.

QAlgorithms.Papers.Independent.Intakethird100.P54.coherentEnergy

For any natural number MMM, including zero, and any complex amplitudes α=(αm)m=0M−1\alpha=(\alpha_m)_{m=0}^{M-1}α=(αm​)m=0M−1​, this defines the real number

∑m=0M−1∣αm∣2=∑m=0M−1((Re⁡αm)2+(Im⁡αm)2).\sum_{m=0}^{M-1}|\alpha_m|^2 = \sum_{m=0}^{M-1}\left((\operatorname{Re}\alpha_m)^2+(\operatorname{Im}\alpha_m)^2\right).m=0∑M−1​∣αm​∣2=m=0∑M−1​((Reαm​)2+(Imαm​)2).

For M=0M=0M=0, this is the empty sum and equals 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.coherentExponent

For any pair ccc of splitters, each having reflectance Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and an arbitrary real phase, the definition assigns

(Rtrue−Rfalse)2,\left(\sqrt{R_{\mathrm{true}}}-\sqrt{R_{\mathrm{false}}}\right)^2,(Rtrue​​−Rfalse​​)2,

where the square roots are nonnegative. No ordering or distinctness assumption is imposed; equal reflectances give zero, and the phases do not affect the value.

QAlgorithms.Papers.Independent.Intakethird100.P54.coherentOpticalProbe

For arbitrary natural numbers J,MJ,MJ,M, complex amplitudes β=(βj)j=0J−1\beta=(\beta_j)_{j=0}^{J-1}β=(βj​)j=0J−1​ and α=(αm)m=0M−1\alpha=(\alpha_m)_{m=0}^{M-1}α=(αm​)m=0M−1​, this defines a vector in ℓ2(NJ×NM;C)\ell^2(\mathbb N^J\times\mathbb N^M;\mathbb C)ℓ2(NJ×NM;C) from the coefficients

a(u,v)=(∏j=0J−1e−∣βj∣2/2βjujuj!)(∏m=0M−1e−∣αm∣2/2αmvmvm!),(u,v)∈NJ×NM.a(u,v)= \left(\prod_{j=0}^{J-1}e^{-|\beta_j|^2/2}\frac{\beta_j^{u_j}}{\sqrt{u_j!}}\right) \left(\prod_{m=0}^{M-1}e^{-|\alpha_m|^2/2}\frac{\alpha_m^{v_m}}{\sqrt{v_m!}}\right), \qquad (u,v)\in\mathbb N^J\times\mathbb N^M.a(u,v)=(j=0∏J−1​e−∣βj​∣2/2uj​!​βjuj​​​)(m=0∏M−1​e−∣αm​∣2/2vm​!​αmvm​​​),(u,v)∈NJ×NM.

If this coefficient function is square-summable, the defined vector has these coefficients; otherwise it is the zero vector. Either number of modes may be zero, and each empty product equals 111. Natural powers use 00=10^0=100=1, and every factorial denominator is positive.

QAlgorithms.Papers.Independent.Intakethird100.P54.countLaw

For any natural number MMM, any splitter aaa consisting of a real reflectance R∈[0,1]R\in[0,1]R∈[0,1] and an arbitrary real phase, and any occupation vectors N,j∈NMN,j\in\mathbb N^MN,j∈NM, this defines the real number

{∏m=0M−1(Nmjm)Rjm(1−R)Nm−jm,jm≤Nm for every m,0,otherwise.\begin{cases} \displaystyle\prod_{m=0}^{M-1}\binom{N_m}{j_m}R^{j_m}(1-R)^{N_m-j_m}, &j_m\le N_m\text{ for every }m,\\[6pt] 0,&\text{otherwise}. \end{cases}⎩⎨⎧​m=0∏M−1​(jm​Nm​​)Rjm​(1−R)Nm​−jm​,0,​jm​≤Nm​ for every m,otherwise.​

The phase does not enter this expression. The exponents use natural-number subtraction, which truncates at zero, although the branch condition ensures Nm−jmN_m-j_mNm​−jm​ is the ordinary nonnegative difference wherever the product is used. The boundary reflectances R=0R=0R=0 and R=1R=1R=1 are included, with 00=10^0=100=1. For M=0M=0M=0, the branch condition is vacuous and the empty product equals 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.countThreshold

For any pair of splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false) and c1=c(true)c_1=c(\mathrm{true})c1​=c(true), with respective reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and arbitrary real phases, and any natural number KKK, this defines the real number

Klog⁡ ⁣(1−R01−R1)log⁡ ⁣(R1(1−R0)R0(1−R1)).\frac{ K\log\!\left(\dfrac{1-R_0}{1-R_1}\right) }{ \log\!\left(\dfrac{R_1(1-R_0)}{R_0(1-R_1)}\right) }.log(R0​(1−R1​)R1​(1−R0​)​)Klog(1−R1​1−R0​​)​.

The phases do not enter the expression, and no ordering, distinctness, or strict interior condition is imposed on the reflectances. All divisions and logarithms are total: division by zero equals 000, and log⁡0=0\log 0=0log0=0. Thus the expression is defined even at boundary reflectances and when its final denominator is zero; in the latter case its value is 000. For K=0K=0K=0, its value is also 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.countingError

For any pair of splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false) and c1=c(true)c_1=c(\mathrm{true})c1​=c(true), each having a real reflectance Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and an arbitrary real phase, and any natural number KKK, this defines

12∑j=0Kmin⁡ ⁣{(Kj)R0j(1−R0)K−j,(Kj)R1j(1−R1)K−j}.\frac12\sum_{j=0}^{K} \min\!\left\{ \binom Kj R_0^j(1-R_0)^{K-j}, \binom Kj R_1^j(1-R_1)^{K-j} \right\}.21​j=0∑K​min{(jK​)R0j​(1−R0​)K−j,(jK​)R1j​(1−R1​)K−j}.

No ordering or distinctness of the splitters, and no restriction on their phases, is assumed; only their reflectances enter the expression. The endpoints Rb=0,1R_b=0,1Rb​=0,1 are included, with 00=10^0=100=1. For K=0K=0K=0, the sum has one term and the defined value is 12\tfrac1221​.

QAlgorithms.Papers.Independent.Intakethird100.P54.efficiencyScaled

For every splitter aaa, whose reflectance RRR satisfies 0≤R≤10\le R\le10≤R≤1 and whose phase θ\thetaθ is an arbitrary real number, and every real η\etaη with supplied proofs of 0≤η≤10\le\eta\le10≤η≤1, the definition constructs a splitter with reflectance ηR\eta RηR and phase θ\thetaθ, together with proofs that 0≤ηR≤10\le\eta R\le10≤ηR≤1. The cases η=0\eta=0η=0, η=1\eta=1η=1, R=0R=0R=0, and R=1R=1R=1 are all included; at η=0\eta=0η=0, the resulting reflectance is zero while its phase remains θ\thetaθ.

QAlgorithms.Papers.Independent.Intakethird100.P54.eprAlpha

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, and for every real NNN, the definition assigns

(1−R0)N+1(1−R1)N+1.\frac{(1-R_0)N+1}{(1-R_1)N+1}.(1−R1​)N+1(1−R0​)N+1​.

There is no nonnegativity assumption on NNN, nor any ordering or distinctness assumption on the splitters. If the denominator is zero, the quotient is defined to be zero; the phases do not enter the formula.

QAlgorithms.Papers.Independent.Intakethird100.P54.eprBeta

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, the definition assigns

1−R01−R1.\frac{1-R_0}{1-R_1}.1−R1​1−R0​​.

No ordering or distinctness is assumed, and the phases do not enter the formula. The endpoint R1=1R_1=1R1​=1 is included and gives zero under the convention that division by zero is zero.

QAlgorithms.Papers.Independent.Intakethird100.P54.eprC

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, and for every real NNN, the definition assigns

((1−R0R1)N+1)2(1−R0)N+1,\frac{\left((1-\sqrt{R_0}\sqrt{R_1})N+1\right)^2} {(1-R_0)N+1},(1−R0​)N+1((1−R0​​R1​​)N+1)2​,

using nonnegative square roots. No sign restriction on NNN, ordering, or distinctness is assumed. A zero denominator gives zero, and the phases do not affect the value.

QAlgorithms.Papers.Independent.Intakethird100.P54.eprD

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, and for every real NNN, the definition assigns

(1−R1)N.(1-R_1)N.(1−R1​)N.

There is no nonnegativity assumption on NNN, nor any ordering or distinctness requirement. The value is zero when N=0N=0N=0 or R1=1R_1=1R1​=1, and neither phase enters the formula.

QAlgorithms.Papers.Independent.Intakethird100.P54.eprDerivativeSign

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, and for every real NNN, the definition assigns the real number

(1−R1)Nlog⁡ ⁣(1−R01−R1)−((1−R0R1)N+1)2(1−R0)N+1log⁡ ⁣((1−R0)N+1(1−R1)N+1).(1-R_1)N\log\!\left(\frac{1-R_0}{1-R_1}\right) - \frac{\left((1-\sqrt{R_0}\sqrt{R_1})N+1\right)^2} {(1-R_0)N+1} \log\!\left(\frac{(1-R_0)N+1}{(1-R_1)N+1}\right).(1−R1​)Nlog(1−R1​1−R0​​)−(1−R0​)N+1((1−R0​​R1​​)N+1)2​log((1−R1​)N+1(1−R0​)N+1​).

No sign restriction on NNN, ordering, or distinctness is assumed. Division by zero is zero, log⁡0=0\log0=0log0=0, and log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣ for negative xxx. The phases do not enter this expression, and the definition itself asserts no derivative identity or sign property.

QAlgorithms.Papers.Independent.Intakethird100.P54.eprOptimizer

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, and for every real NNN, define the auxiliary real numbers

A=(1−R0)N+1(1−R1)N+1,B=1−R01−R1,C=((1−R0R1)N+1)2(1−R0)N+1,D=(1−R1)N.A=\frac{(1-R_0)N+1}{(1-R_1)N+1},\qquad B=\frac{1-R_0}{1-R_1},\qquad C=\frac{\left((1-\sqrt{R_0}\sqrt{R_1})N+1\right)^2}{(1-R_0)N+1}, \qquad D=(1-R_1)N.A=(1−R1​)N+1(1−R0​)N+1​,B=1−R1​1−R0​​,C=(1−R0​)N+1((1−R0​​R1​​)N+1)2​,D=(1−R1​)N.

The definition assigns

{0,0≤Dlog⁡B−Clog⁡A,min⁡ ⁣{1, log⁡ ⁣(Clog⁡ADlog⁡B)log⁡(B/A)},Dlog⁡B−Clog⁡A<0.\begin{cases} 0,&0\le D\log B-C\log A,\\[2mm] \displaystyle\min\!\left\{1,\, \frac{\log\!\left(\frac{C\log A}{D\log B}\right)} {\log(B/A)}\right\},&D\log B-C\log A<0. \end{cases}⎩⎨⎧​0,min⎩⎨⎧​1,log(B/A)log(DlogBClogA​)​⎭⎬⎫​,​0≤DlogB−ClogA,DlogB−ClogA<0.​

Equality in the branch condition selects zero. No sign restriction on NNN, ordering, or distinctness is assumed. All divisions and logarithms are total, with division by zero equal to zero, log⁡0=0\log0=0log0=0, and log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣ for negative xxx. The second branch caps the expression above by 111 without imposing a lower bound of zero, and no minimizing property is asserted.

QAlgorithms.Papers.Independent.Intakethird100.P54.exponentGain

For any pair ccc of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1], indexed respectively by false\mathrm{false}false and true\mathrm{true}true, and arbitrary real phases, set

s∗=log⁡ ⁣(−R1log⁡(R0/R1)(1−R1)log⁡((1−R0)/(1−R1)))log⁡ ⁣((1−R0)R1(1−R1)R0).s_*= \frac{\log\!\left( \frac{-R_1\log(R_0/R_1)} {(1-R_1)\log((1-R_0)/(1-R_1))} \right)} {\log\!\left(\frac{(1-R_0)R_1}{(1-R_1)R_0}\right)}.s∗​=log((1−R1​)R0​(1−R0​)R1​​)log((1−R1​)log((1−R0​)/(1−R1​))−R1​log(R0​/R1​)​)​.

The definition assigns

−log⁡ ⁣(R0s∗R11−s∗+(1−R0)s∗(1−R1)1−s∗)(R1−R0)2.\frac{-\log\!\left( R_0^{s_*}R_1^{1-s_*}+ (1-R_0)^{s_*}(1-R_1)^{1-s_*} \right)} {(\sqrt{R_1}-\sqrt{R_0})^2}.(R1​​−R0​​)2−log(R0s∗​​R11−s∗​​+(1−R0​)s∗​(1−R1​)1−s∗​)​.

No ordering, distinctness, or phase restriction is assumed. Division by zero is zero, log⁡0=0\log0=0log0=0, and log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣ for negative xxx; powers are real powers with 00=10^0=100=1 and 0t=00^t=00t=0 for t≠0t\ne0t=0. In particular, equal reflectances make the denominator zero and the defined quotient zero.

QAlgorithms.Papers.Independent.Intakethird100.P54.fock

For any natural number M≥0M\ge0M≥0 and any occupation tuple N:{0,…,M−1}→NN:\{0,\ldots,M-1\}\to\mathbb NN:{0,…,M−1}→N, this defines a vector in ℓ2(N0×NM;C)\ell^2(\mathcal N_0\times\mathcal N_M;\mathbb C)ℓ2(N0​×NM​;C), where Nj={0,…,j−1}→N\mathcal N_j=\{0,\ldots,j-1\}\to\mathbb NNj​={0,…,j−1}→N, by

fN(a,n)={1,(a,n)=(0,N),0,otherwise.f_N(a,n)= \begin{cases} 1,&(a,n)=(0,N),\\ 0,&\text{otherwise}. \end{cases}fN​(a,n)={1,0,​(a,n)=(0,N),otherwise.​

Here N0\mathcal N_0N0​ contains exactly one function, denoted 000, so this is the unit coordinate vector at the pair consisting of that empty function and NNN. When M=0M=0M=0, the entire coordinate index set is a singleton and the vector’s sole coordinate is 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.fockCoefficient

For any pair of splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false) and c1=c(true)c_1=c(\mathrm{true})c1​=c(true), with respective real reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and arbitrary real phases, and any real number sss, this defines

R0 sR1 1−s+(1−R0) s(1−R1) 1−s,R_0^{\,s}R_1^{\,1-s}+(1-R_0)^{\,s}(1-R_1)^{\,1-s},R0s​R11−s​+(1−R0​)s(1−R1​)1−s,

where the powers are real powers. There is no restriction that sss lie in [0,1][0,1][0,1], nor any ordering or distinctness assumption on the reflectances; the phases do not enter the formula. For a positive base xxx, xt=exp⁡(tlog⁡x)x^t=\exp(t\log x)xt=exp(tlogx); at a zero base, the convention is 00=10^0=100=1 and 0t=00^t=00t=0 for every nonzero real ttt, including negative ttt.

QAlgorithms.Papers.Independent.Intakethird100.P54.fockOptimizer

For any pair of splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false) and c1=c(true)c_1=c(\mathrm{true})c1​=c(true), with respective reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and arbitrary real phases, this defines the real number

log⁡ ⁣(−R1log⁡(R0/R1)(1−R1)log⁡((1−R0)/(1−R1)))log⁡ ⁣((1−R0)R1(1−R1)R0).\frac{ \log\!\left( \frac{-R_1\log(R_0/R_1)} {(1-R_1)\log((1-R_0)/(1-R_1))} \right) }{ \log\!\left( \frac{(1-R_0)R_1}{(1-R_1)R_0} \right) }.log((1−R1​)R0​(1−R0​)R1​​)log((1−R1​)log((1−R0​)/(1−R1​))−R1​log(R0​/R1​)​)​.

Only the reflectances enter this expression. There is no ordering, distinctness, strict interior, or logarithm-argument positivity assumption. Division by zero equals 000; the real logarithm satisfies log⁡0=0\log 0=0log0=0 and, for negative arguments, log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣. Consequently the formula assigns a real value even at boundary reflectances or vanishing denominators, and a zero final denominator makes the defined value 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.fockVector

For any type XXX and function a:X→Ca:X\to\mathbb Ca:X→C, this defines a vector in ℓ2(X;C)\ell^2(X;\mathbb C)ℓ2(X;C): if aaa is square-summable, the vector has coordinates a(x)a(x)a(x); otherwise, it is the zero vector. No normalization is required. If XXX is empty, the result is the unique zero vector.

QAlgorithms.Papers.Independent.Intakethird100.P54.kernelProduct

For any type XXX and any two complex kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, this defines their kernel product by

(ρ⋆σ)(x,y)=∑z∈Xρ(x,z)σ(z,y)(x,y∈X).(\rho\mathbin{\star}\sigma)(x,y)=\sum_{z\in X}\rho(x,z)\sigma(z,y) \qquad(x,y\in X).(ρ⋆σ)(x,y)=z∈X∑​ρ(x,z)σ(z,y)(x,y∈X).

There is no finiteness, countability, boundedness, positivity, or summability assumption. The sum is the total complex summation operation, equal to zero when its summand is not summable. Empty index types are included.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossAmplitude

For every natural number MMM, every splitter aaa with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and real phase θ\thetaθ, and every pair of occupation functions n,k:{0,…,M−1}→Nn,k:\{0,\ldots,M-1\}\to\mathbb Nn,k:{0,…,M−1}→N, this defines the complex number

{∏m=0M−1((nmkm) einmθ (R) nm−km (1−R) km),if km≤nm for every m,0,otherwise.\begin{cases} \displaystyle\prod_{m=0}^{M-1} \left(\sqrt{\binom{n_m}{k_m}}\, e^{i n_m\theta}\, (\sqrt R)^{\,n_m-k_m}\, (\sqrt{1-R})^{\,k_m}\right), & \text{if }k_m\le n_m\text{ for every }m,\\[4pt] 0,&\text{otherwise}. \end{cases}⎩⎨⎧​m=0∏M−1​((km​nm​​)​einm​θ(R​)nm​−km​(1−R​)km​),0,​if km​≤nm​ for every m,otherwise.​

The subtraction in each exponent is natural-number subtraction, and the condition in the first branch ensures it equals the ordinary nonnegative difference. Zero exponents give the value 111, including for a zero base, so the formula includes R=0R=0R=0 and R=1R=1R=1. For M=0M=0M=0, the condition is vacuous and the empty product equals 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossBlock

For any type III, natural number M≥0M\ge0M≥0, splitter pair ccc, vector ψ∈ℓ2(I×NM;C)\psi\in\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈ℓ2(I×NM​;C), and k∈NM={0,…,M−1}→Nk\in\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb Nk∈NM​={0,…,M−1}→N, write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for the reflectance and phase of c(b)c(b)c(b), for b=0,1b=0,1b=0,1. Let Lb,kL_{b,k}Lb,k​ be the vector whose coordinate at (i,q)(i,q)(i,q) is

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m},ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,

with the entire vector defined as zero if this coordinate function is not square-summable. The definition returns the complex linear subspace

span⁡C{L0,k,L1,k}={αL0,k+βL1,k:α,β∈C}.\operatorname{span}_{\mathbb C}\{L_{0,k},L_{1,k}\} =\{\alpha L_{0,k}+\beta L_{1,k}:\alpha,\beta\in\mathbb C\}.spanC​{L0,k​,L1,k​}={αL0,k​+βL1,k​:α,β∈C}.

The two vectors need not be nonzero or independent, so the subspace may have dimension zero, one, or two. No normalization of ψ\psiψ is required. The case M=0M=0M=0 uses the unique empty occupation tuple and product 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossBlockDifference

Let III be any type, M≥0M\ge0M≥0 a natural number, NM={0,…,M−1}→N\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb NNM​={0,…,M−1}→N, ccc a splitter pair, ψ,v∈ℓ2(I×NM;C)\psi,v\in\ell^2(I\times\mathcal N_M;\mathbb C)ψ,v∈ℓ2(I×NM​;C), k∈NMk\in\mathcal N_Mk∈NM​, and w0,w1∈Rw_0,w_1\in\mathbb Rw0​,w1​∈R. Write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for the reflectance and phase of c(b)c(b)c(b). Define Lb,kL_{b,k}Lb,k​, for b=0,1b=0,1b=0,1, from the coordinates

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m},ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,

using the zero vector if this function is not square-summable. This definition returns the vector

w0⟨L0,k,v⟩L0,k−w1⟨L1,k,v⟩L1,k,w_0\langle L_{0,k},v\rangle L_{0,k} -w_1\langle L_{1,k},v\rangle L_{1,k},w0​⟨L0,k​,v⟩L0,k​−w1​⟨L1,k​,v⟩L1,k​,

where the inner product is conjugate-linear in its first argument. The real weights are unrestricted, neither vector must be normalized, and vvv need not lie in the span of the two sector vectors. For M=0M=0M=0, occupation tuples are empty and the product is 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossOverlap

For any type III, natural number M≥0M\ge0M≥0, vector ψ∈ℓ2(I×NM;C)\psi\in\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈ℓ2(I×NM​;C), splitter pair ccc, and k∈NMk\in\mathcal N_Mk∈NM​, where NM={0,…,M−1}→N\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb NNM​={0,…,M−1}→N, write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for the reflectance and phase of splitter c(b)c(b)c(b), for b=0,1b=0,1b=0,1. For each bbb, let Lb,kL_{b,k}Lb,k​ have coordinates

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m}ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​

when this function is square-summable, and let Lb,k=0L_{b,k}=0Lb,k​=0 otherwise. This defines the complex number ⟨L0,k,L1,k⟩\langle L_{0,k},L_{1,k}\rangle⟨L0,k​,L1,k​⟩, using the Hilbert inner product conjugate-linear in its first argument. The vector ψ\psiψ need not be normalized. The case M=0M=0M=0 is included, with one empty occupation tuple and empty product 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossOverlapTerm

For any type III, natural number M≥0M\ge0M≥0, splitter pair ccc, vector ψ∈ℓ2(I×NM;C)\psi\in\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈ℓ2(I×NM​;C), and occupation tuples k,n∈NM={0,…,M−1}→Nk,n\in\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb Nk,n∈NM​={0,…,M−1}→N, put

pψ(n)=∑i∈I∣ψ(i,n)∣2.p_\psi(n)=\sum_{i\in I}|\psi(i,n)|^2.pψ​(n)=i∈I∑​∣ψ(i,n)∣2.

For b=0,1b=0,1b=0,1, write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for splitter c(b)c(b)c(b)’s reflectance and phase, and define

Ab(n,k)={∏m=0M−1(nmkm) e inmθb(Rb)nm−km(1−Rb)km,km≤nm for every m,0,otherwise.A_b(n,k)= \begin{cases} \displaystyle\prod_{m=0}^{M-1} \sqrt{\binom{n_m}{k_m}}\, e^{\,\mathrm i n_m\theta_b} (\sqrt{R_b})^{n_m-k_m}(\sqrt{1-R_b})^{k_m}, & k_m\le n_m\text{ for every }m,\\ 0,&\text{otherwise}. \end{cases}Ab​(n,k)=⎩⎨⎧​m=0∏M−1​(km​nm​​)​einm​θb​(Rb​​)nm​−km​(1−Rb​​)km​,0,​km​≤nm​ for every m,otherwise.​

The definition returns pψ(n) A0(n,k)‾A1(n,k)p_\psi(n)\,\overline{A_0(n,k)}A_1(n,k)pψ​(n)A0​(n,k)​A1​(n,k), as a complex number. In particular, it returns zero if any km>nmk_m>n_mkm​>nm​. No normalization or orthogonality condition on ψ\psiψ is imposed. For M=0M=0M=0, both tuples are the unique empty tuple, the condition is vacuous, and both amplitude products equal 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossProbability

Let III be any type, let M≥0M\ge0M≥0 be a natural number, and write NM={0,…,M−1}→N\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb NNM​={0,…,M−1}→N. Let ψ∈ℓ2(I×NM;C)\psi\in\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈ℓ2(I×NM​;C), let ccc assign a splitter to each Boolean bbb, and fix bbb and k∈NMk\in\mathcal N_Mk∈NM​. Write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for that splitter’s reflectance and phase. Define Lb,kL_{b,k}Lb,k​ to be the square-summable vector with coordinates

Lb,k(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,L_{b,k}(i,q)=\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m},Lb,k​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,

if this coordinate function is square-summable, and the zero vector otherwise. The definition returns ∥Lb,k∥2\|L_{b,k}\|^2∥Lb,k​∥2. No normalization of ψ\psiψ, ordering of the two reflectances, or relation between their phases is assumed. When M=0M=0M=0, NM\mathcal N_MNM​ contains just the empty tuple and the product is 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.lossSector

For any type III, natural number MMM, splitter with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and phase θ∈R\theta\in\mathbb Rθ∈R, square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), and occupation vector k∈N{0,…,M−1}k\in\mathbb N^{\{0,\ldots,M-1\}}k∈N{0,…,M−1}, form the complex family

sk(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θ(R)qm(1−R)km.s_k(i,q)=\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta} (\sqrt R)^{q_m}(\sqrt{1-R})^{k_m}.sk​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θ(R​)qm​(1−R​)km​.

The defined vector is this family regarded as an element of ℓ2(I×N{0,…,M−1})\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ℓ2(I×N{0,…,M−1}) if it is square-summable, and is the zero vector otherwise. No unit-norm condition is imposed on ψ\psiψ or on the resulting vector. Endpoint reflectances and zero exponents are included, with 00=10^0=100=1; for M=0M=0M=0, the product is empty and the resulting vector is ψ\psiψ.

QAlgorithms.Papers.Independent.Intakethird100.P54.lowerError

For every real number FFF, this defines

lowerError⁡(F)=1−1−F2.\operatorname{lowerError}(F)=\frac{1-\sqrt{1-F}}{2}.lowerError(F)=21−1−F​​.

There is no restriction 0≤F≤10\le F\le10≤F≤1. The real square root is defined to be zero on negative inputs, so the value is 1/21/21/2 whenever F≥1F\ge1F≥1.

QAlgorithms.Papers.Independent.Intakethird100.P54.ndsChernoffTerm

For any type III, natural number M≥0M\ge0M≥0, splitter pair ccc, vector ψ∈ℓ2(I×NM;C)\psi\in\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈ℓ2(I×NM​;C), real sss, and k∈NM={0,…,M−1}→Nk\in\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb Nk∈NM​={0,…,M−1}→N, write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for splitter c(b)c(b)c(b)’s parameters. Define Lb,kL_{b,k}Lb,k​, for b=0,1b=0,1b=0,1, from the coordinates

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m},ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,

using zero if this function is not square-summable, and put pb=∥Lb,k∥2p_b=\|L_{b,k}\|^2pb​=∥Lb,k​∥2. This defines

{0,p0=0 or p1=0,p0 s−1p1−s ∣⟨L0,k,L1,k⟩∣2,otherwise,\begin{cases} 0,&p_0=0\text{ or }p_1=0,\\ p_0^{\,s-1}p_1^{-s}\,|\langle L_{0,k},L_{1,k}\rangle|^2, &\text{otherwise}, \end{cases}{0,p0s−1​p1−s​∣⟨L0,k​,L1,k​⟩∣2,​p0​=0 or p1​=0,otherwise,​

where powers are real powers. In the second case both bases are positive. The parameter sss is not restricted to [0,1][0,1][0,1], and no normalization or occupation-sector orthogonality of ψ\psiψ is assumed. The case M=0M=0M=0 uses the unique empty tuple and product 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.numberExponent

For any pair ccc of splitters indexed by b∈{false,true}b\in\{\mathrm{false},\mathrm{true}\}b∈{false,true}, each with reflectance Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and an arbitrary real phase, write R0=RfalseR_0=R_{\mathrm{false}}R0​=Rfalse​ and R1=RtrueR_1=R_{\mathrm{true}}R1​=Rtrue​, and set

s∗=log⁡ ⁣(−R1log⁡(R0/R1)(1−R1)log⁡((1−R0)/(1−R1)))log⁡ ⁣((1−R0)R1(1−R1)R0).s_*= \frac{\log\!\left( \frac{-R_1\log(R_0/R_1)} {(1-R_1)\log((1-R_0)/(1-R_1))} \right)} {\log\!\left(\frac{(1-R_0)R_1}{(1-R_1)R_0}\right)}.s∗​=log((1−R1​)R0​(1−R0​)R1​​)log((1−R1​)log((1−R0​)/(1−R1​))−R1​log(R0​/R1​)​)​.

The definition assigns

−log⁡ ⁣(R0s∗R11−s∗+(1−R0)s∗(1−R1)1−s∗).-\log\!\left(R_0^{s_*}R_1^{1-s_*} +(1-R_0)^{s_*}(1-R_1)^{1-s_*}\right).−log(R0s∗​​R11−s∗​​+(1−R0​)s∗​(1−R1​)1−s∗​).

The phases do not enter this expression. No ordering or distinctness of the reflectances is required, and endpoint reflectances are included. All expressions use total real operations: division by zero is zero, log⁡0=0\log 0=0log0=0, and for negative xxx, log⁡x=log⁡∣x∣\log x=\log|x|logx=log∣x∣; the powers are real powers, with 00=10^0=100=1 and 0t=00^t=00t=0 for t≠0t\ne0t=0. The definition makes no assertion that s∗s_*s∗​ minimizes any function or belongs to [0,1][0,1][0,1].

QAlgorithms.Papers.Independent.Intakethird100.P54.opticalBayesRisk

For any type XXX, let H=ℓ2(X;C)H=\ell^2(X;\mathbb C)H=ℓ2(X;C) and let ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C be arbitrary kernels. A resolution of a kernel KKK consists of weights wn≥0w_n\ge0wn​≥0 and vectors un∈Hu_n\in Hun​∈H, indexed by n∈Nn\in\mathbb Nn∈N, such that the weights are summable with sum 111, every vector with positive weight has norm 111, distinct vectors with positive weights are orthogonal, and K(x,y)=∑n=0∞wnun(x)un(y)‾K(x,y)=\sum_{n=0}^{\infty}w_nu_n(x)\overline{u_n(y)}K(x,y)=∑n=0∞​wn​un​(x)un​(y)​ for every x,yx,yx,y. Given resolutions A=(an,un)A=(a_n,u_n)A=(an​,un​) of ρ\rhoρ and B=(bn,vn)B=(b_n,v_n)B=(bn​,vn​) of σ\sigmaσ, any continuous complex-linear operator E:H→HE:H\to HE:H→H, and any real numbers π0,π1\pi_0,\pi_1π0​,π1​, the definition is

opticalBayesRisk⁡(A,B,E,π0,π1)=π0∑n=0∞an(1−Re⁡⟨un,Eun⟩)+π1∑n=0∞bnRe⁡⟨vn,Evn⟩.\operatorname{opticalBayesRisk}(A,B,E,\pi_0,\pi_1) =\pi_0\sum_{n=0}^{\infty}a_n\bigl(1-\operatorname{Re}\langle u_n,Eu_n\rangle\bigr) +\pi_1\sum_{n=0}^{\infty}b_n\operatorname{Re}\langle v_n,Ev_n\rangle.opticalBayesRisk(A,B,E,π0​,π1​)=π0​n=0∑∞​an​(1−Re⟨un​,Eun​⟩)+π1​n=0∑∞​bn​Re⟨vn​,Evn​⟩.

There is no assumption here that EEE is an effect or that the real numbers π0,π1\pi_0,\pi_1π0​,π1​ are nonnegative or sum to 111. Zero-weight vectors need not be normalized or orthogonal. The infinite sums use the total convention that a nonsummable family has sum 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.opticalChernoff

For any type XXX and arbitrary kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, put H=ℓ2(X;C)H=\ell^2(X;\mathbb C)H=ℓ2(X;C) and ht(p)=0h_t(p)=0ht​(p)=0 for p=0p=0p=0, with ht(p)=pth_t(p)=p^tht​(p)=pt otherwise. For each s∈[0,1]s\in[0,1]s∈[0,1], form the real supremum

Cs=sup⁡(an,un) resolves ρ(bm,vm) resolves σ{∑(n,m)∈N2hs(an) h1−s(bm) ∣⟨un,vm⟩∣2}.C_s= \sup_{\substack{(a_n,u_n)\text{ resolves }\rho\\(b_m,v_m)\text{ resolves }\sigma}} \left\{ \sum_{(n,m)\in\mathbb N^2} h_s(a_n)\,h_{1-s}(b_m)\,|\langle u_n,v_m\rangle|^2 \right\}.Cs​=(an​,un​) resolves ρ(bm​,vm​) resolves σ​sup​⎩⎨⎧​(n,m)∈N2∑​hs​(an​)h1−s​(bm​)∣⟨un​,vm​⟩∣2⎭⎬⎫​.

Here a resolution (wn,un)(w_n,u_n)(wn​,un​) of KKK has nonnegative summable weights summing to 111, vectors in HHH normalized whenever their weights are positive, orthogonality between distinct positive-weight vectors, and K(x,y)=∑nwnun(x)un(y)‾K(x,y)=\sum_nw_nu_n(x)\overline{u_n(y)}K(x,y)=∑n​wn​un​(x)un​(y)​ for all x,yx,yx,y. The definition is opticalChernoff⁡(ρ,σ)=inf⁡{Cs:0≤s≤1}\operatorname{opticalChernoff}(\rho,\sigma)=\inf\{C_s:0\le s\le1\}opticalChernoff(ρ,σ)=inf{Cs​:0≤s≤1}, a real infimum; both endpoints are included, and zero weights retain value 000 at exponent zero. No existence of resolutions is assumed. Empty or unbounded-above candidate sets for the inner real supremum have value 000, and a nonsummable infinite family has sum 000; in particular, empty XXX gives value 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.opticalCopies

For arbitrary natural numbers J,M,LJ,M,LJ,M,L, including zero, and any vector ψ∈ℓ2(NJ×NM;C)\psi\in\ell^2(\mathbb N^J\times\mathbb N^M;\mathbb C)ψ∈ℓ2(NJ×NM;C), this defines a vector in ℓ2((NJ)L×NLM;C)\ell^2((\mathbb N^J)^L\times\mathbb N^{LM};\mathbb C)ℓ2((NJ)L×NLM;C) from the coefficient function

a(i,n)=∏ℓ=0L−1ψ ⁣(iℓ,(nℓM+m)m=0M−1),a(i,n)=\prod_{\ell=0}^{L-1}\psi\!\left(i_\ell,\bigl(n_{\ell M+m}\bigr)_{m=0}^{M-1}\right),a(i,n)=ℓ=0∏L−1​ψ(iℓ​,(nℓM+m​)m=0M−1​),

where i=(iℓ)ℓ=0L−1i=(i_\ell)_{\ell=0}^{L-1}i=(iℓ​)ℓ=0L−1​ is a list of LLL occupation vectors in NJ\mathbb N^JNJ, and n∈NLMn\in\mathbb N^{LM}n∈NLM. If aaa is square-summable, the defined vector has exactly these coefficients; otherwise it is the zero vector. There is no normalization assumption on ψ\psiψ. When L=0L=0L=0, the product is empty and equals 111, and the output index set is a singleton; J=0J=0J=0 and M=0M=0M=0 are also allowed, with their occupation vectors interpreted as functions on the empty set.

QAlgorithms.Papers.Independent.Intakethird100.P54.opticalError

For any type XXX and arbitrary kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, let H=ℓ2(X;C)H=\ell^2(X;\mathbb C)H=ℓ2(X;C). The value opticalError⁡(ρ,σ)\operatorname{opticalError}(\rho,\sigma)opticalError(ρ,σ) is the real infimum of

12∑n=0∞an(1−Re⁡⟨un,Eun⟩)+12∑n=0∞bnRe⁡⟨vn,Evn⟩\frac12\sum_{n=0}^{\infty}a_n\bigl(1-\operatorname{Re}\langle u_n,Eu_n\rangle\bigr) +\frac12\sum_{n=0}^{\infty}b_n\operatorname{Re}\langle v_n,Ev_n\rangle21​n=0∑∞​an​(1−Re⟨un​,Eun​⟩)+21​n=0∑∞​bn​Re⟨vn​,Evn​⟩

over the following data: nonnegative summable weights an,bna_n,b_nan​,bn​ each summing to 111; vectors un,vn∈Hu_n,v_n\in Hun​,vn​∈H normalized at positive weights and orthogonal within each family between distinct positive-weight indices; representations ρ(x,y)=∑nanun(x)un(y)‾\rho(x,y)=\sum_na_nu_n(x)\overline{u_n(y)}ρ(x,y)=∑n​an​un​(x)un​(y)​ and σ(x,y)=∑nbnvn(x)vn(y)‾\sigma(x,y)=\sum_nb_nv_n(x)\overline{v_n(y)}σ(x,y)=∑n​bn​vn​(x)vn​(y)​ for every x,yx,yx,y; and a continuous complex-linear operator E:H→HE:H\to HE:H→H such that ⟨u,Ev⟩=⟨v,Eu⟩‾\langle u,Ev\rangle=\overline{\langle v,Eu\rangle}⟨u,Ev⟩=⟨v,Eu⟩​ for all u,vu,vu,v and 0≤Re⁡⟨v,Ev⟩≤∥v∥20\le\operatorname{Re}\langle v,Ev\rangle\le\|v\|^20≤Re⟨v,Ev⟩≤∥v∥2 for all vvv. Thus the two weights multiplying the risks are exactly 1/21/21/2. Existence of these data is not assumed; the real infimum is 000 if the candidate set is empty, including when XXX is empty. Infinite sums have value 000 for nonsummable families.

QAlgorithms.Papers.Independent.Intakethird100.P54.opticalKernel

For any type XXX, an optical kernel on XXX is an arbitrary function X×X→CX\times X\to\mathbb CX×X→C. No finiteness, symmetry, positivity, or normalization condition is imposed; XXX may be empty.

QAlgorithms.Papers.Independent.Intakethird100.P54.opticalPureKernel

For every type XXX and every square-summable complex function ψ∈ℓ2(X;C)\psi\in\ell^2(X;\mathbb C)ψ∈ℓ2(X;C), this defines the complex-valued kernel on X×XX\times XX×X by

Kψ(x,y)=ψ(x)ψ(y)‾(x,y∈X).K_\psi(x,y)=\psi(x)\overline{\psi(y)} \qquad(x,y\in X).Kψ​(x,y)=ψ(x)ψ(y)​(x,y∈X).

No unit-norm assumption is imposed: the zero vector gives the zero kernel. Empty XXX is also allowed, in which case the kernel has no entries.

QAlgorithms.Papers.Independent.Intakethird100.P54.outputChernoff

For any type III, M∈NM\in\mathbb NM∈N, splitters cbc_bcb​ indexed by b∈{0,1}b\in\{0,1\}b∈{0,1} with reflectances Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and phases θb∈R\theta_b\in\mathbb Rθb​∈R, and any ψ∈ℓ2(X;C)\psi\in\ell^2(X;\mathbb C)ψ∈ℓ2(X;C) where Q=N{0,…,M−1}Q=\mathbb N^{\{0,\ldots,M-1\}}Q=N{0,…,M−1} and X=I×QX=I\times QX=I×Q, define

fb,k(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,Kb(x,y)=∑k∈Qfb,k(x)fb,k(y)‾.f_{b,k}(i,q)=\psi(i,q+k) \prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m}, \qquad K_b(x,y)=\sum_{k\in Q}f_{b,k}(x)\overline{f_{b,k}(y)}.fb,k​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,Kb​(x,y)=k∈Q∑​fb,k​(x)fb,k​(y)​.

For each s∈[0,1]s\in[0,1]s∈[0,1], let CsC_sCs​ be the real supremum of ∑(n,m)∈N2hs(an)h1−s(bm)∣⟨un,vm⟩∣2\sum_{(n,m)\in\mathbb N^2}h_s(a_n)h_{1-s}(b_m)|\langle u_n,v_m\rangle|^2∑(n,m)∈N2​hs​(an​)h1−s​(bm​)∣⟨un​,vm​⟩∣2 over resolutions (an,un)(a_n,u_n)(an​,un​) of K0K_0K0​ and (bn,vn)(b_n,v_n)(bn​,vn​) of K1K_1K1​, where ht(0)=0h_t(0)=0ht​(0)=0 and ht(p)=pth_t(p)=p^tht​(p)=pt for p≠0p\ne0p=0. A resolution (wn,zn)(w_n,z_n)(wn​,zn​) consists of nonnegative summable weights summing to 111, vectors in ℓ2(X;C)\ell^2(X;\mathbb C)ℓ2(X;C) normalized at positive weights and orthogonal between distinct positive-weight indices, and K(x,y)=∑nwnzn(x)zn(y)‾K(x,y)=\sum_nw_nz_n(x)\overline{z_n(y)}K(x,y)=∑n​wn​zn​(x)zn​(y)​ for every x,yx,yx,y. The definition is outputChernoff⁡(c,ψ)=inf⁡{Cs:0≤s≤1}\operatorname{outputChernoff}(c,\psi)=\inf\{C_s:0\le s\le1\}outputChernoff(c,ψ)=inf{Cs​:0≤s≤1} in R\mathbb RR, including both endpoints and retaining h0(0)=0h_0(0)=0h0​(0)=0. No normalization of ψ\psiψ, ordering or distinctness of splitters, or restriction on phases is assumed. The case M=0M=0M=0 uses a singleton QQQ and empty products equal to 111. Infinite sums are 000 for nonsummable families, and empty or unbounded-above candidate sets for the inner real supremum have value 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.outputChernoffAt

For any type III, M∈NM\in\mathbb NM∈N, splitters cbc_bcb​ indexed by b∈{0,1}b\in\{0,1\}b∈{0,1} with reflectances Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and phases θb∈R\theta_b\in\mathbb Rθb​∈R, any ψ∈ℓ2(X;C)\psi\in\ell^2(X;\mathbb C)ψ∈ℓ2(X;C) where Q=N{0,…,M−1}Q=\mathbb N^{\{0,\ldots,M-1\}}Q=N{0,…,M−1} and X=I×QX=I\times QX=I×Q, and any real sss, define

fb,k(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,Kb(x,y)=∑k∈Qfb,k(x)fb,k(y)‾.f_{b,k}(i,q)=\psi(i,q+k) \prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m}, \qquad K_b(x,y)=\sum_{k\in Q}f_{b,k}(x)\overline{f_{b,k}(y)}.fb,k​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,Kb​(x,y)=k∈Q∑​fb,k​(x)fb,k​(y)​.

With ht(p)=0h_t(p)=0ht​(p)=0 if p=0p=0p=0 and ht(p)=pth_t(p)=p^tht​(p)=pt otherwise, outputChernoffAt⁡(c,ψ,s)\operatorname{outputChernoffAt}(c,\psi,s)outputChernoffAt(c,ψ,s) is the real supremum of

∑(n,m)∈N2hs(an)h1−s(bm)∣⟨un,vm⟩∣2\sum_{(n,m)\in\mathbb N^2}h_s(a_n)h_{1-s}(b_m)|\langle u_n,v_m\rangle|^2(n,m)∈N2∑​hs​(an​)h1−s​(bm​)∣⟨un​,vm​⟩∣2

over resolutions (an,un)(a_n,u_n)(an​,un​) of K0K_0K0​ and (bn,vn)(b_n,v_n)(bn​,vn​) of K1K_1K1​. A resolution (wn,zn)(w_n,z_n)(wn​,zn​) has nonnegative summable weights summing to 111, vectors in ℓ2(X;C)\ell^2(X;\mathbb C)ℓ2(X;C) normalized at positive weights and orthogonal between distinct positive-weight indices, and the representation K(x,y)=∑nwnzn(x)zn(y)‾K(x,y)=\sum_nw_nz_n(x)\overline{z_n(y)}K(x,y)=∑n​wn​zn​(x)zn​(y)​ for every x,yx,yx,y. There is no restriction on sss, no normalization assumption on ψ\psiψ, and no ordering, distinctness, or phase restriction on the splitters. Zero weights give 000 even at zero or negative exponents. The case M=0M=0M=0 is allowed, with empty products equal to 111. Nonsummable families have sum 000, and an empty or unbounded-above candidate set has real supremum 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.outputError

For any type III, M∈NM\in\mathbb NM∈N, splitters cbc_bcb​ indexed by b∈{0,1}b\in\{0,1\}b∈{0,1} with reflectances Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and phases θb∈R\theta_b\in\mathbb Rθb​∈R, and any ψ∈H=ℓ2(X;C)\psi\in H=\ell^2(X;\mathbb C)ψ∈H=ℓ2(X;C) where Q=N{0,…,M−1}Q=\mathbb N^{\{0,\ldots,M-1\}}Q=N{0,…,M−1} and X=I×QX=I\times QX=I×Q, set

fb,k(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,Kb(x,y)=∑k∈Qfb,k(x)fb,k(y)‾.f_{b,k}(i,q)=\psi(i,q+k) \prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m}, \qquad K_b(x,y)=\sum_{k\in Q}f_{b,k}(x)\overline{f_{b,k}(y)}.fb,k​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,Kb​(x,y)=k∈Q∑​fb,k​(x)fb,k​(y)​.

The value outputError⁡(c,ψ)\operatorname{outputError}(c,\psi)outputError(c,ψ) is the real infimum of

12∑n=0∞an(1−Re⁡⟨un,Eun⟩)+12∑n=0∞bnRe⁡⟨vn,Evn⟩\frac12\sum_{n=0}^{\infty}a_n\bigl(1-\operatorname{Re}\langle u_n,Eu_n\rangle\bigr) +\frac12\sum_{n=0}^{\infty}b_n\operatorname{Re}\langle v_n,Ev_n\rangle21​n=0∑∞​an​(1−Re⟨un​,Eun​⟩)+21​n=0∑∞​bn​Re⟨vn​,Evn​⟩

over nonnegative summable weights an,bna_n,b_nan​,bn​ each summing to 111, vectors un,vn∈Hu_n,v_n\in Hun​,vn​∈H normalized at positive weights and orthogonal within each family between distinct positive-weight indices, representations K0(x,y)=∑nanun(x)un(y)‾K_0(x,y)=\sum_na_nu_n(x)\overline{u_n(y)}K0​(x,y)=∑n​an​un​(x)un​(y)​ and K1(x,y)=∑nbnvn(x)vn(y)‾K_1(x,y)=\sum_nb_nv_n(x)\overline{v_n(y)}K1​(x,y)=∑n​bn​vn​(x)vn​(y)​, and continuous complex-linear operators E:H→HE:H\to HE:H→H satisfying ⟨u,Ev⟩=⟨v,Eu⟩‾\langle u,Ev\rangle=\overline{\langle v,Eu\rangle}⟨u,Ev⟩=⟨v,Eu⟩​ for all u,vu,vu,v and 0≤Re⁡⟨v,Ev⟩≤∥v∥20\le\operatorname{Re}\langle v,Ev\rangle\le\|v\|^20≤Re⟨v,Ev⟩≤∥v∥2 for all vvv. No normalization of ψ\psiψ, ordering or distinctness of the splitters, or restriction on phases is required. The case M=0M=0M=0 uses a singleton QQQ and empty products equal to 111. If no admissible resolutions and operator exist, the real infimum is 000; nonsummable infinite families have sum 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.outputFidelity

For any type III, any natural number MMM, two splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false) and c1=c(true)c_1=c(\mathrm{true})c1​=c(true) with reflectances Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and phases θb∈R\theta_b\in\mathbb Rθb​∈R, and any ψ∈ℓ2(I×N{0,…,M−1};C)\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}};\mathbb C)ψ∈ℓ2(I×N{0,…,M−1};C), put Q=N{0,…,M−1}Q=\mathbb N^{\{0,\ldots,M-1\}}Q=N{0,…,M−1}, X=I×QX=I\times QX=I×Q, rb=Rbr_b=\sqrt{R_b}rb​=Rb​​, and tb=1−Rbt_b=\sqrt{1-R_b}tb​=1−Rb​​. Define, for k,q∈Qk,q\in Qk,q∈Q,

fb,k(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb rbqmtbkm,Kb(x,y)=∑k∈Qfb,k(x)fb,k(y)‾.f_{b,k}(i,q)=\psi(i,q+k) \prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b}\, r_b^{q_m}t_b^{k_m}, \qquad K_b(x,y)=\sum_{k\in Q}f_{b,k}(x)\overline{f_{b,k}(y)}.fb,k​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​rbqm​​tbkm​​,Kb​(x,y)=k∈Q∑​fb,k​(x)fb,k​(y)​.

Then outputFidelity⁡(c,ψ)\operatorname{outputFidelity}(c,\psi)outputFidelity(c,ψ) is the real supremum of ∣⟨z,w⟩∣2|\langle z,w\rangle|^2∣⟨z,w⟩∣2 over all z,w∈ℓ2(X×N;C)z,w\in\ell^2(X\times\mathbb N;\mathbb C)z,w∈ℓ2(X×N;C) with ∥z∥=∥w∥=1\|z\|=\|w\|=1∥z∥=∥w∥=1 and, for every x,y∈Xx,y\in Xx,y∈X,

∑e=0∞z(x,e)z(y,e)‾=K0(x,y),∑e=0∞w(x,e)w(y,e)‾=K1(x,y).\sum_{e=0}^{\infty}z(x,e)\overline{z(y,e)}=K_0(x,y), \qquad \sum_{e=0}^{\infty}w(x,e)\overline{w(y,e)}=K_1(x,y).e=0∑∞​z(x,e)z(y,e)​=K0​(x,y),e=0∑∞​w(x,e)w(y,e)​=K1​(x,y).

No normalization of ψ\psiψ, ordering or distinctness of the splitters, or restriction on their phases is assumed. The case M=0M=0M=0 is included, with QQQ a singleton and empty products equal to 111. The real supremum is 000 if no such purifications exist; infinite sums have value 000 for nonsummable families.

QAlgorithms.Papers.Independent.Intakethird100.P54.photons

For every natural number MMM, including M=0M=0M=0, and every occupation function n:{0,…,M−1}→Nn:\{0,\ldots,M-1\}\to\mathbb Nn:{0,…,M−1}→N, this defines the natural number ∑m=0M−1nm\sum_{m=0}^{M-1}n_m∑m=0M−1​nm​. For M=0M=0M=0, the sum is empty and equals zero.

QAlgorithms.Papers.Independent.Intakethird100.P54.productOpticalProbe

For any natural number M≥0M\ge0M≥0, any function J:{0,…,M−1}→NJ:\{0,\ldots,M-1\}\to\mathbb NJ:{0,…,M−1}→N, and any family

ψm∈ℓ2(NJ(m)×N1;C),Nj={0,…,j−1}→N,\psi_m\in\ell^2(\mathcal N_{J(m)}\times\mathcal N_1;\mathbb C), \qquad \mathcal N_j=\{0,\ldots,j-1\}\to\mathbb N,ψm​∈ℓ2(NJ(m)​×N1​;C),Nj​={0,…,j−1}→N,

this defines a vector in ℓ2((∏m=0M−1NJ(m))×NM;C)\ell^2((\prod_{m=0}^{M-1}\mathcal N_{J(m)})\times\mathcal N_M;\mathbb C)ℓ2((∏m=0M−1​NJ(m)​)×NM​;C) from the coordinate function

f(a,n)=∏m=0M−1ψm(am,nm‾),f(a,n)=\prod_{m=0}^{M-1}\psi_m(a_m,\underline{n_m}),f(a,n)=m=0∏M−1​ψm​(am​,nm​​),

where am∈NJ(m)a_m\in\mathcal N_{J(m)}am​∈NJ(m)​, n∈NMn\in\mathcal N_Mn∈NM​, and nm‾\underline{n_m}nm​​ is the function on the singleton {0}\{0\}{0} with value nmn_mnm​. The result has these coordinates if fff is square-summable, and is the zero vector otherwise. No normalization of the input vectors is required, and any J(m)J(m)J(m) may be zero. For M=0M=0M=0, both coordinate factors are singletons and the empty product gives the vector with sole coordinate 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.purificationMarginal

For any type XXX and square-summable complex family z∈ℓ2(X×N)z\in\ell^2(X\times\mathbb N)z∈ℓ2(X×N), this defines the complex-valued kernel on XXX by

Kz(x,y)=∑e=0∞z(x,e) z(y,e)‾(x,y∈X).K_z(x,y)=\sum_{e=0}^{\infty}z(x,e)\,\overline{z(y,e)} \qquad(x,y\in X).Kz​(x,y)=e=0∑∞​z(x,e)z(y,e)​(x,y∈X).

No unit-norm condition is imposed on zzz. In particular, z=0z=0z=0 gives the zero kernel, and if XXX is empty the resulting kernel has no entries.

QAlgorithms.Papers.Independent.Intakethird100.P54.purificationOverlap

For any type III, M∈NM\in\mathbb NM∈N, two splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false), c1=c(true)c_1=c(\mathrm{true})c1​=c(true) with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and arbitrary real phases θ0,θ1\theta_0,\theta_1θ0​,θ1​, and any ψ∈ℓ2(I×Q;C)\psi\in\ell^2(I\times Q;\mathbb C)ψ∈ℓ2(I×Q;C) with Q=N{0,…,M−1}Q=\mathbb N^{\{0,\ldots,M-1\}}Q=N{0,…,M−1}, put

Pn=∑q∈Q{∑i∈I∣ψ(i,q)∣2,∑m=0M−1qm=n,0,otherwise,κ=R0R1+1−R01−R1.P_n=\sum_{q\in Q} \begin{cases} \displaystyle\sum_{i\in I}|\psi(i,q)|^2,&\displaystyle\sum_{m=0}^{M-1}q_m=n,\\ 0,&\text{otherwise}, \end{cases} \qquad \kappa=\sqrt{R_0}\sqrt{R_1}+\sqrt{1-R_0}\sqrt{1-R_1}.Pn​=q∈Q∑​⎩⎨⎧​i∈I∑​∣ψ(i,q)∣2,0,​m=0∑M−1​qm​=n,otherwise,​κ=R0​​R1​​+1−R0​​1−R1​​.

The definition is the real number

purificationOverlap⁡(c,ψ)=∣∑n=0∞Pne in(θ1−θ0)κn∣2.\operatorname{purificationOverlap}(c,\psi) =\left|\sum_{n=0}^{\infty}P_n e^{\,\mathrm i n(\theta_1-\theta_0)}\kappa^n\right|^2.purificationOverlap(c,ψ)=​n=0∑∞​Pn​ein(θ1​−θ0​)κn​2.

No normalization or energy assumption on ψ\psiψ, or ordering or distinctness assumption on the splitters, is required. The sum includes n=0n=0n=0, whose exponential and power factors are 111, including when κ=0\kappa=0κ=0. The case M=0M=0M=0 is included: the sole occupation has photon count 000, so the expression becomes (∑i∈I∣ψ(i,0)∣2)2\bigl(\sum_{i\in I}|\psi(i,0)|^2\bigr)^2(∑i∈I​∣ψ(i,0)∣2)2, where 000 denotes the empty occupation. Infinite sums have value 000 for nonsummable families.

QAlgorithms.Papers.Independent.Intakethird100.P54.purificationOverlapTerm

For any type III, M∈NM\in\mathbb NM∈N, two splitters c0=c(false)c_0=c(\mathrm{false})c0​=c(false), c1=c(true)c_1=c(\mathrm{true})c1​=c(true) with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and arbitrary real phases θ0,θ1\theta_0,\theta_1θ0​,θ1​, any ψ∈ℓ2(I×Q;C)\psi\in\ell^2(I\times Q;\mathbb C)ψ∈ℓ2(I×Q;C) with Q=N{0,…,M−1}Q=\mathbb N^{\{0,\ldots,M-1\}}Q=N{0,…,M−1}, and any n∈Nn\in\mathbb Nn∈N, define

Pn=∑q∈Q{∑i∈I∣ψ(i,q)∣2,∑m=0M−1qm=n,0,otherwise,κ=R0R1+1−R01−R1.P_n=\sum_{q\in Q} \begin{cases} \displaystyle\sum_{i\in I}|\psi(i,q)|^2,&\displaystyle\sum_{m=0}^{M-1}q_m=n,\\ 0,&\text{otherwise}, \end{cases} \qquad \kappa=\sqrt{R_0}\sqrt{R_1}+\sqrt{1-R_0}\sqrt{1-R_1}.Pn​=q∈Q∑​⎩⎨⎧​i∈I∑​∣ψ(i,q)∣2,0,​m=0∑M−1​qm​=n,otherwise,​κ=R0​​R1​​+1−R0​​1−R1​​.

Then purificationOverlapTerm⁡(c,ψ,n)\operatorname{purificationOverlapTerm}(c,\psi,n)purificationOverlapTerm(c,ψ,n) is the complex number Pne in(θ1−θ0)κnP_ne^{\,\mathrm i n(\theta_1-\theta_0)}\kappa^nPn​ein(θ1​−θ0​)κn. No normalization or energy assumption on ψ\psiψ, or ordering or distinctness assumption on the splitters, is imposed. The index n=0n=0n=0 is included, with κ0=1\kappa^0=1κ0=1 even when κ=0\kappa=0κ=0. If M=0M=0M=0, QQQ contains just the empty occupation and its photon count is 000. All infinite sums use the total convention of value 000 for nonsummable families.

QAlgorithms.Papers.Independent.Intakethird100.P54.returned

For any type III, natural number MMM, splitter with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and phase θ∈R\theta\in\mathbb Rθ∈R, and square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), define

sk(i,q)=ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θ(R)qm(1−R)km.s_k(i,q)=\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta} (\sqrt R)^{q_m}(\sqrt{1-R})^{k_m}.sk​(i,q)=ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θ(R​)qm​(1−R​)km​.

The returned kernel is the complex-valued function of two indices x,y∈I×N{0,…,M−1}x,y\in I\times\mathbb N^{\{0,\ldots,M-1\}}x,y∈I×N{0,…,M−1} given by

K(x,y)=∑k∈N{0,…,M−1}sk(x) sk(y)‾.K(x,y)=\sum_{k\in\mathbb N^{\{0,\ldots,M-1\}}} s_k(x)\,\overline{s_k(y)}.K(x,y)=k∈N{0,…,M−1}∑​sk​(x)sk​(y)​.

These are the coefficient families themselves, without applying the square-summability test used to construct a loss-sector vector. No normalization of ψ\psiψ is assumed. The reflectances 000 and 111 are included, with 00=10^0=100=1. When M=0M=0M=0, the sum has one occupation-vector index and K(x,y)=ψ(x)ψ(y)‾K(x,y)=\psi(x)\overline{\psi(y)}K(x,y)=ψ(x)ψ(y)​.

QAlgorithms.Papers.Independent.Intakethird100.P54.sectorCoefficients

For any type III, natural number MMM, splitter aaa with reflectance R∈[0,1]R\in[0,1]R∈[0,1] and phase θ∈R\theta\in\mathbb Rθ∈R, square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), occupation vector k∈N{0,…,M−1}k\in\mathbb N^{\{0,\ldots,M-1\}}k∈N{0,…,M−1}, and x=(i,q)x=(i,q)x=(i,q) in the same index set as ψ\psiψ, this defines the complex number

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θ(R)qm(1−R)km.\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta} (\sqrt R)^{q_m}(\sqrt{1-R})^{k_m}.ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θ(R​)qm​(1−R​)km​.

Addition of occupation vectors is componentwise. The conditional in the underlying amplitude takes its product branch because km≤qm+kmk_m\le q_m+k_mkm​≤qm​+km​ for every mmm. No normalization of ψ\psiψ is required. The definition includes R=0R=0R=0, R=1R=1R=1, and zero exponents, with 00=10^0=100=1. When M=0M=0M=0, the occupation vectors are the unique empty vectors and the empty product is 111, so the coefficient is ψ(i,q)\psi(i,q)ψ(i,q).

QAlgorithms.Papers.Independent.Intakethird100.P54.sectorRisk

For any type III, natural number M≥0M\ge0M≥0, splitter pair ccc, vector ψ∈H=ℓ2(I×NM;C)\psi\in H=\ell^2(I\times\mathcal N_M;\mathbb C)ψ∈H=ℓ2(I×NM​;C), bounded complex-linear operator E:H→HE:H\to HE:H→H, and real numbers π0,π1\pi_0,\pi_1π0​,π1​, where NM={0,…,M−1}→N\mathcal N_M=\{0,\ldots,M-1\}\to\mathbb NNM​={0,…,M−1}→N, write Rb∈[0,1]R_b\in[0,1]Rb​∈[0,1] and θb∈R\theta_b\in\mathbb Rθb​∈R for splitter c(b)c(b)c(b)’s reflectance and phase. Let Lb,kL_{b,k}Lb,k​ have coordinates

ψ(i,q+k)∏m=0M−1(qm+kmkm) e i(qm+km)θb(Rb)qm(1−Rb)km,\psi(i,q+k)\prod_{m=0}^{M-1} \sqrt{\binom{q_m+k_m}{k_m}}\, e^{\,\mathrm i(q_m+k_m)\theta_b} (\sqrt{R_b})^{q_m}(\sqrt{1-R_b})^{k_m},ψ(i,q+k)m=0∏M−1​(km​qm​+km​​)​ei(qm​+km​)θb​(Rb​​)qm​(1−Rb​​)km​,

or be zero if this coordinate function is not square-summable. The defined real number is

π0∑k∈NM(∥L0,k∥2−Re⁡⟨L0,k,EL0,k⟩)+π1∑k∈NMRe⁡⟨L1,k,EL1,k⟩.\pi_0\sum_{k\in\mathcal N_M} \left(\|L_{0,k}\|^2-\operatorname{Re}\langle L_{0,k},EL_{0,k}\rangle\right) +\pi_1\sum_{k\in\mathcal N_M} \operatorname{Re}\langle L_{1,k},EL_{1,k}\rangle.π0​k∈NM​∑​(∥L0,k​∥2−Re⟨L0,k​,EL0,k​⟩)+π1​k∈NM​∑​Re⟨L1,k​,EL1,k​⟩.

There is no assumption that EEE is an effect, that ψ\psiψ is normalized, or that the weights are nonnegative or sum to one. Each infinite sum uses the total summation convention, which assigns zero to a nonsummable family. For M=0M=0M=0, the index set has one empty tuple and the coordinate product is 111.

QAlgorithms.Papers.Independent.Intakethird100.P54.signalLaw

For any type III, natural number MMM, square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), and occupation vector n∈N{0,…,M−1}n\in\mathbb N^{\{0,\ldots,M-1\}}n∈N{0,…,M−1}, this defines the real number

pψ(n)=∑i∈I∣ψ(i,n)∣2.p_\psi(n)=\sum_{i\in I}|\psi(i,n)|^2.pψ​(n)=i∈I∑​∣ψ(i,n)∣2.

The type III need not be countable or nonempty, and ψ\psiψ need not have norm 111; in particular, these values are not required by the definition to have total mass 111. If III is empty, the value is 000. If M=0M=0M=0, there is one occupation vector and its value is ∥ψ∥2\|\psi\|^2∥ψ∥2.

QAlgorithms.Papers.Independent.Intakethird100.P54.splitterOverlap

For any channel pair ccc, consisting of splitters with reflectances R0,R1∈[0,1]R_0,R_1\in[0,1]R0​,R1​∈[0,1] and arbitrary real phases, this defines the real number R0R1+1−R01−R1\sqrt{R_0}\sqrt{R_1}+\sqrt{1-R_0}\sqrt{1-R_1}R0​​R1​​+1−R0​​1−R1​​. The phases do not enter this expression, and no validity or zero-phase hypothesis is required.

QAlgorithms.Papers.Independent.Intakethird100.P54.squaredOpticalFidelity

For any type XXX and arbitrary complex-valued kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, this defines the real number

sup⁡{∣⟨z,w⟩∣2:z,w∈ℓ2(X×N),∥z∥=∥w∥=1,ρ(x,y)=∑e=0∞z(x,e)z(y,e)‾for all x,y∈X,σ(x,y)=∑e=0∞w(x,e)w(y,e)‾for all x,y∈X}.\sup\left\{ |\langle z,w\rangle|^2: \begin{array}{l} z,w\in\ell^2(X\times\mathbb N),\quad \|z\|=\|w\|=1,\\ \displaystyle\rho(x,y)=\sum_{e=0}^{\infty}z(x,e)\overline{z(y,e)} \quad\text{for all }x,y\in X,\\ \displaystyle\sigma(x,y)=\sum_{e=0}^{\infty}w(x,e)\overline{w(y,e)} \quad\text{for all }x,y\in X \end{array} \right\}.sup⎩⎨⎧​∣⟨z,w⟩∣2:z,w∈ℓ2(X×N),∥z∥=∥w∥=1,ρ(x,y)=e=0∑∞​z(x,e)z(y,e)​for all x,y∈X,σ(x,y)=e=0∑∞​w(x,e)w(y,e)​for all x,y∈X​⎭⎬⎫​.

Here ⟨⋅,⋅⟩\langle\cdot,\cdot\rangle⟨⋅,⋅⟩ is the complex Hilbert-space inner product. Existence of the two unit vectors is not assumed. If either kernel has no such representation, the set is empty and the total real supremum used by the definition has value 000; this includes the case X=∅X=\varnothingX=∅. Whenever the set is nonempty, its elements lie in [0,1][0,1][0,1].

QAlgorithms.Papers.Independent.Intakethird100.P54.supportWeightPower

For arbitrary real numbers p,tp,tp,t, the definition is

supportWeightPower⁡(p,t)={0,p=0,pt,p≠0,\operatorname{supportWeightPower}(p,t)= \begin{cases} 0,&p=0,\\ p^t,&p\ne0, \end{cases}supportWeightPower(p,t)={0,pt,​p=0,p=0,​

where the second branch uses the total real-power function. In particular, the value at p=0p=0p=0 is 000 for every exponent, including t=0t=0t=0 and negative exponents. No nonnegativity assumption on ppp or restriction on ttt is imposed.

QAlgorithms.Papers.Independent.Intakethird100.P54.targetCoherent

For arbitrary real numbers R,NR,NR,N, this defines

targetCoherent⁡(R,N)=1− 1−e−RN 2.\operatorname{targetCoherent}(R,N) =\frac{1-\sqrt{\,1-e^{-RN}\,}}{2}.targetCoherent(R,N)=21−1−e−RN​​.

No sign or range restriction is imposed on either argument. The real square root equals zero on negative arguments, so the value is 1/21/21/2 when RN<0RN<0RN<0, and it is 000 when RN=0RN=0RN=0.

QAlgorithms.Papers.Independent.Intakethird100.P54.targetLower

For arbitrary real numbers R,NR,NR,N, this defines

targetLower⁡(R,N)=1− 1−(1−R)N 2,\operatorname{targetLower}(R,N) =\frac{1-\sqrt{\,1-(1-R)^N\,}}{2},targetLower(R,N)=21−1−(1−R)N​​,

where the power is the total real-power operation. There is no assumption that R∈[0,1]R\in[0,1]R∈[0,1] or that N≥0N\ge0N≥0. The real square root is also total and equals zero on negative arguments; thus the defined value is 1/21/21/2 whenever (1−R)N>1(1-R)^N>1(1−R)N>1.

QAlgorithms.Papers.Independent.Intakethird100.P54.tmsvCoefficients

For any natural number MMM, any real number NNN, and any pair of occupation vectors u,v∈NMu,v\in\mathbb N^Mu,v∈NM, this defines the complex coefficient

aM,N(u,v)={∏m=0M−1(1N+1)(NN+1)vm,u=v,0,u≠v,a_{M,N}(u,v)= \begin{cases} \displaystyle\prod_{m=0}^{M-1} \left(\sqrt{\frac1{N+1}}\right) \left(\sqrt{\frac{N}{N+1}}\right)^{v_m},&u=v,\\[6pt] 0,&u\ne v, \end{cases}aM,N​(u,v)=⎩⎨⎧​m=0∏M−1​(N+11​​)(N+1N​​)vm​,0,​u=v,u=v,​

with the real factors regarded as complex numbers. There is no restriction that NNN be nonnegative: division by zero is assigned the value 000, the real square root of a negative number is 000, and natural powers satisfy 00=10^0=100=1. For M=0M=0M=0, both occupation vectors are the unique empty function, so they are equal and the empty product gives coefficient 111 for every real NNN.

QAlgorithms.Papers.Independent.Intakethird100.P54.tmsvProbe

For any natural number MMM and any real number NNN, this defines a vector in ℓ2(NM×NM;C)\ell^2(\mathbb N^M\times\mathbb N^M;\mathbb C)ℓ2(NM×NM;C) using the coefficient function

a(u,v)={∏m=0M−11N+1(NN+1)vm,u=v,0,u≠v.a(u,v)= \begin{cases} \displaystyle\prod_{m=0}^{M-1} \sqrt{\frac1{N+1}} \left(\sqrt{\frac{N}{N+1}}\right)^{v_m},&u=v,\\[6pt] 0,&u\ne v. \end{cases}a(u,v)=⎩⎨⎧​m=0∏M−1​N+11​​(N+1N​​)vm​,0,​u=v,u=v.​

If aaa is square-summable, the vector has exactly these coefficients; otherwise it is the zero vector. The real factors are embedded into C\mathbb CC. No nonnegativity assumption is imposed on NNN; division by zero equals 000, real square roots of negative numbers equal 000, and 00=10^0=100=1. When M=0M=0M=0, the coefficient function is 111 on its singleton index set, and hence defines the vector with that single coefficient.

QAlgorithms.Papers.Independent.Intakethird100.P54.totalCountLaw

For any splitter aaa with real reflectance R∈[0,1]R\in[0,1]R∈[0,1] and arbitrary real phase, and any natural numbers K,jK,jK,j, including zero, this defines

{(Kj)Rj(1−R)K−j,j≤K,0,j>K.\begin{cases} \displaystyle\binom Kj R^j(1-R)^{K-j},&j\le K,\\[4pt] 0,&j>K. \end{cases}⎩⎨⎧​(jK​)Rj(1−R)K−j,0,​j≤K,j>K.​

The phase does not enter the formula. Natural-number subtraction truncates at zero, but in the first branch K−jK-jK−j is the ordinary nonnegative difference. The reflectance endpoints are permitted and natural powers use 00=10^0=100=1; in particular, for K=0K=0K=0 the value is 111 at j=0j=0j=0 and 000 at every j>0j>0j>0.

QAlgorithms.Papers.Independent.Intakethird100.P54.totalLaw

For any type III, natural number MMM, square-summable complex family ψ∈ℓ2(I×N{0,…,M−1})\psi\in\ell^2(I\times\mathbb N^{\{0,\ldots,M-1\}})ψ∈ℓ2(I×N{0,…,M−1}), and natural number nnn, this defines

Pψ(n)=∑q∈N{0,…,M−1}{∑i∈I∣ψ(i,q)∣2,∑m=0M−1qm=n,0,otherwise.P_\psi(n)= \sum_{q\in\mathbb N^{\{0,\ldots,M-1\}}} \begin{cases} \displaystyle\sum_{i\in I}|\psi(i,q)|^2, &\displaystyle\sum_{m=0}^{M-1}q_m=n,\\ 0,&\text{otherwise}. \end{cases}Pψ​(n)=q∈N{0,…,M−1}∑​⎩⎨⎧​i∈I∑​∣ψ(i,q)∣2,0,​m=0∑M−1​qm​=n,otherwise.​

Thus it sums the squared coordinate magnitudes over occupation vectors whose total occupation is exactly nnn. There is no unit-norm assumption on ψ\psiψ, and n=0n=0n=0 is included. When M=0M=0M=0, the unique occupation vector has total occupation 000, so Pψ(0)=∥ψ∥2P_\psi(0)=\|\psi\|^2Pψ​(0)=∥ψ∥2 and Pψ(n)=0P_\psi(n)=0Pψ​(n)=0 for every n>0n>0n>0. If III is empty, all these values are 000.

QAlgorithms.Papers.Independent.Intakethird100.P54.vacuumEfficiency

For every real η\etaη, given proofs of 0≤η0\le\eta0≤η and η≤1\eta\le1η≤1, the definition constructs a splitter whose reflectance is η\etaη and whose phase is 000, retaining those proofs as the required reflectance bounds. Both endpoints η=0\eta=0η=0 and η=1\eta=1η=1 are included.

QAlgorithms.Papers.Independent.Intakethird100.P54.weightedTraceNorm

For any type XXX, any complex kernels ρ,σ:X×X→C\rho,\sigma:X\times X\to\mathbb Cρ,σ:X×X→C, and arbitrary real weights w0,w1w_0,w_1w0​,w1​, let H=ℓ2(X;C)H=\ell^2(X;\mathbb C)H=ℓ2(X;C). A resolution of a kernel KKK consists of weights pn≥0p_n\ge0pn​≥0, indexed by n∈Nn\in\mathbb Nn∈N, with ∑npn=1\sum_n p_n=1∑n​pn​=1 and summable weight sequence, and vectors an∈Ha_n\in Han​∈H such that ∥an∥=1\|a_n\|=1∥an​∥=1 whenever pn>0p_n>0pn​>0, distinct vectors with positive weights are orthogonal, and K(x,y)=∑npnan(x)an(y)‾K(x,y)=\sum_n p_n a_n(x)\overline{a_n(y)}K(x,y)=∑n​pn​an​(x)an​(y)​ for every x,yx,yx,y. Vectors with zero weights need not be normalized or satisfy the orthogonality condition. The definition is

sup⁡(p,a),(q,b),T[w0∑npnRe⁡⟨an,Tan⟩−w1∑nqnRe⁡⟨bn,Tbn⟩],\sup_{(p,a),(q,b),T} \left[ w_0\sum_n p_n\operatorname{Re}\langle a_n,Ta_n\rangle -w_1\sum_n q_n\operatorname{Re}\langle b_n,Tb_n\rangle \right],(p,a),(q,b),Tsup​[w0​n∑​pn​Re⟨an​,Tan​⟩−w1​n∑​qn​Re⟨bn​,Tbn​⟩],

where (p,a)(p,a)(p,a) resolves ρ\rhoρ, (q,b)(q,b)(q,b) resolves σ\sigmaσ, and T:H→HT:H\to HT:H→H is bounded and complex linear, satisfies ⟨u,Tv⟩=⟨v,Tu⟩‾\langle u,Tv\rangle=\overline{\langle v,Tu\rangle}⟨u,Tv⟩=⟨v,Tu⟩​ for every u,vu,vu,v, and satisfies ∣Re⁡⟨u,Tu⟩∣≤∥u∥2\left|\operatorname{Re}\langle u,Tu\rangle\right|\le\|u\|^2∣Re⟨u,Tu⟩∣≤∥u∥2 for every uuu. The weights w0,w1w_0,w_1w0​,w1​ are not required to be nonnegative or sum to 111, and the kernels are not assumed to possess resolutions. If either has no resolution, the supremum is over the empty set and the defined real value is 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