Research portal

Mathematical note

Operator-qualified optical inference contract

math/operator-qualified-optical-inference.md

Edition
Site v0.3.0 · continuous main snapshot
Source revision
ec2865b0eac15148675c629981a545632b3571c5
Extent
2,146 words
Public route
https://www.cordana.dev/math/operator-qualified-optical-inference/
Mapped records2 mapped records

Direct repository links only; no document-level evidence status is implied.

This note defines the quantitative boundary for Fixture F-007. It operationalizes the measurement-operator-aware residual from the optics, photonics, and inverse-sensing audit. The contract binds every decoded output to a versioned physical operator, separates measurement information from prior selection, and compares optical, digital, and hybrid routes at equal task information and lifecycle budget.

Episode and operator identity

For episode ee and acquisition tt, seal

Ie,t=(Xe,Ae,t,νe,t,Ce,t,Re,t,Qe,Be),I_{e,t}=(X_e,A_{e,t},\nu_{e,t},C_{e,t},R_{e,t},Q_e,B_e),

where XeX_e identifies the physical scene or latent-state generator, Ae,tA_{e,t} identifies the acquisition action, νe,t\nu_{e,t} is an immutable operator-version identifier, Ce,tC_{e,t} is the calibration record, Re,tR_{e,t} is the hidden regime record, QeQ_e is the registered downstream-query set, and BeB_e is the resource-budget record. Identifiers and hashes are byte strings; timestamps within the records are seconds [s] from a declared clock origin.

Let X\mathcal X be the admissible latent-state space and let xeXx_e\in\mathcal X be the latent physical state. Each component of xex_e retains its native unit, such as radiance [W sr1^{-1} m2^{-2}], range [m], or concentration [mol m3^{-3}]. The observation is

ye,t=ge,t ⁣(Hνe,t(ae,t,ce,t)xe)+ne,t,y_{e,t}=g_{e,t}\!\left(\mathcal H_{\nu_{e,t}} (a_{e,t},c_{e,t})x_e\right)+n_{e,t},

where ae,ta_{e,t} is the acquisition action, ce,tc_{e,t} is the calibrated parameter vector, Hνe,t\mathcal H_{\nu_{e,t}} is the versioned physical forward operator, ge,tg_{e,t} is the detector response including conversion and clipping, ne,tn_{e,t} is additive or conditionally specified noise, and ye,ty_{e,t} is the raw observation in detector counts [count] or another declared sensor unit. The operator carries the conversion units required to map components of xex_e to the input unit of ge,tg_{e,t}. Any non-additive, coherent, correlated, or signal-dependent noise is part of the likelihood rather than being hidden in ne,tn_{e,t}.

The submitted observation contract is

Oe,t=(ye,t,ae,t,νe,t,ce,t,Σe,tc,Me,tsat,τe,t,Ve,t),O_{e,t}=(y_{e,t},a_{e,t},\nu_{e,t},c_{e,t}, \Sigma^{c}_{e,t},M^{\mathrm{sat}}_{e,t}, \tau_{e,t},\mathcal V_{e,t}),

where Σe,tc\Sigma^{c}_{e,t} is calibration-parameter covariance in the squared native parameter units, Me,tsatM^{\mathrm{sat}}_{e,t} is a binary saturation or dead-time mask, τe,t\tau_{e,t} is acquisition time [s], and Ve,t\mathcal V_{e,t} is the declared validity envelope. The contract is invalid outside Ve,t\mathcal V_{e,t} until recalibrated or explicitly downgraded.

Aperture, diffraction, and recoverable modes

For wavelength λ\lambda [m] and numerical aperture NA\mathrm{NA} [dimensionless], the conventional Rayleigh lateral scale for two incoherent point sources is

dR=0.61λNA[m].d_{\mathrm R}=0.61\frac{\lambda}{\mathrm{NA}} \quad [\mathrm m].

dRd_{\mathrm R} is a criterion under stated imaging assumptions, not a universal task-resolution value. Aperture diameter DapD_{\mathrm{ap}} [m], focal length ff [m], field of view Ω\Omega [sr], coherence, sampling pitch psampp_{\mathrm{samp}} [m], exposure, and noise must be reported separately. A reconstructed pixel pitch below dRd_{\mathrm R} does not by itself establish additional measured information.

For a linearized forward operator He,tH_{e,t} with singular-value decomposition

He,t=Ue,tΣe,tVe,t,H_{e,t}=U_{e,t}\Sigma_{e,t}V_{e,t}^{*},

Ue,tU_{e,t} and Ve,tV_{e,t} are unitary bases, Ve,tV_{e,t}^{*} is the conjugate transpose, and diagonal entry σe,t,j\sigma_{e,t,j} of Σe,t\Sigma_{e,t} has the units of He,tH_{e,t}. A state perturbation ve,t,jv_{e,t,j} in the corresponding right singular direction is unidentifiable from that acquisition when σe,t,j=0\sigma_{e,t,j}=0. For a preregistered tolerance ϵH\epsilon_H with the same units as a singular value, define the effective measured rank

rϵH(He,t)=j1[σe,t,j>ϵH][mode],r_{\epsilon_H}(H_{e,t})= \sum_j\mathbb 1[\sigma_{e,t,j}>\epsilon_H] \quad [\mathrm{mode}],

where 1[]\mathbb 1[\cdot] is the indicator function and jj indexes singular modes. Report the full singular spectrum or a validated task-relevant summary; rϵHr_{\epsilon_H} is threshold-qualified.

For two task-distinct states x1x_1 and x2x_2, likelihood separation is

D12=DKL ⁣(p(yx1,a,c,ν)p(yx2,a,c,ν))[nat],D_{12}=D_{\mathrm{KL}}\!\left( p(y\mid x_1,a,c,\nu)\,\|\,p(y\mid x_2,a,c,\nu) \right)\quad [\mathrm{nat}],

where DKLD_{\mathrm{KL}} is Kullback--Leibler divergence in nats. The null-space honesty track treats x1x_1 and x2x_2 as measurement-indistinguishable when D12D_{12} lies below a preregistered discrimination threshold supported by a power calculation. A method must then return calibrated alternatives, a bound, or abstention unless it acquires additional evidence.

The Fixture F-007 analytical likelihood plot visualizes one such indistinguishable base operator and a separating added measurement. It is not an empirical performance result.

Photons, detector response, and dynamic range

For detector element ii, use the photon-counting model when its assumptions hold:

kiPoisson(μi),μi=ηiΦiτi+bi,k_i\sim\operatorname{Poisson}(\mu_i), \qquad \mu_i=\eta_i\Phi_i\tau_i+b_i,

where kik_i is detected count [count], μi\mu_i is expected count [count], ηi\eta_i is quantum or detection efficiency [dimensionless], Φi\Phi_i is incident photon rate [photon/s], τi\tau_i is exposure [s], and bib_i is expected background plus dark count [count]. For the ideal background-free Poisson case,

SNRshot=Nγ,\operatorname{SNR}_{\mathrm{shot}}=\sqrt{N_\gamma},

where NγN_\gamma is expected detected photon count [count] and the signal-to- noise ratio is dimensionless. Read noise [electron rms], fixed-pattern error, coherent receiver noise, afterpulsing, pile-up, and dead time receive explicit terms whenever present.

For full-well or count-rate limit KimaxK_i^{\max} [count], a simplified clipped detector output is

yi=min(ki,Kimax),si=1[kiKimax],y_i=\min(k_i,K_i^{\max}), \qquad s_i=\mathbb 1[k_i\ge K_i^{\max}],

where yiy_i is recorded count [count] and sis_i is a dimensionless saturation indicator. The saturation fraction is

Fsat=1Ndeti=1Ndetsi,F_{\mathrm{sat}}= \frac{1}{N_{\mathrm{det}}} \sum_{i=1}^{N_{\mathrm{det}}}s_i,

where NdetN_{\mathrm{det}} is detector-element count [element] and FsatF_{\mathrm{sat}} is dimensionless. Report full well, count-rate ceiling, read noise, dark signal, analog-to-digital converter range, and dead-time or pile-up policy independently; nominal bit depth is not dynamic range.

Phase ambiguity and prior-qualified reconstruction

For coherent intensity measurement,

y=Ax2+n,y=|Ax|^2+n,

where AA is a declared complex-valued propagation and sampling operator, xx is a complex field amplitude in a declared native unit, yy is intensity or detector count in its native unit, 2|\cdot|^2 is elementwise squared magnitude, and nn is measurement noise in the same unit as yy. The ambiguity class is

E(y;A)={xX:Ax2=Ax2},\mathcal E(y;A)= \{x'\in\mathcal X:|Ax'|^2=|Ax|^2\},

where E\mathcal E is a set of physically admissible fields. Global phase, translation, conjugate inversion, and geometry-specific ambiguities are scored as separate equivalence relations when applicable.

For reconstruction method mm with prior πm(x)\pi_m(x) and likelihood pm(yx,O)p_m(y\mid x,O), the posterior is

pm(xy,O)=pm(yx,O)πm(x)Xpm(yx,O)πm(x)dx,p_m(x\mid y,O)= \frac{p_m(y\mid x,O)\pi_m(x)} {\int_{\mathcal X}p_m(y\mid x',O)\pi_m(x')\,\mathrm dx'},

where OO is the observation contract, xx' is an integration variable with the same native units as xx, and the posterior density carries the reciprocal units implied by the measure dx\mathrm dx'. Method mm must label information origin as measurement, prior, calibration, or active intervention.

For hidden truth xex_e and a nominal (1α)(1-\alpha) credible set Cm,e,1α\mathcal C_{m,e,1-\alpha}, empirical coverage over NN independent episodes is

Cov^m,1α=1Ne=1N1[xeCm,e,1α],\widehat{\operatorname{Cov}}_{m,1-\alpha}= \frac{1}{N}\sum_{e=1}^{N} \mathbb 1[x_e\in\mathcal C_{m,e,1-\alpha}],

where α\alpha and coverage are dimensionless and NN is episode count [episode]. Report coverage after source-family, texture, sparsity, positivity, motion, and noise-model shifts. Perceptual quality and truth fidelity remain separate outcomes.

Active sensing and illumination safety

Let btb_t be the belief state before action aa, θ\theta the uncertain task state, dd a downstream decision, U(d,θ)U(d,\theta) task utility in a declared native unit, and yy the prospective observation. Expected value of information is

EVI(abt)=Eyp(ya,bt)[maxdE[U(d,θ)bt,a,y]]maxdE[U(d,θ)bt].\operatorname{EVI}(a\mid b_t)= \mathbb E_{y\sim p(y\mid a,b_t)} \left[\max_d\mathbb E[U(d,\theta)\mid b_t,a,y]\right] -\max_d\mathbb E[U(d,\theta)\mid b_t].

EVI\operatorname{EVI} has the same unit as UU. Action admissibility is the componentwise condition

c(a)=(Nγ,Ea,La,Da,Wa,Ra)(Nγmax,Eamax,Lamax,Damax,Wamax,Ramax),\mathbf c(a)= (N_\gamma,E_a,L_a,D_a,W_a,R_a) \preceq (N_\gamma^{\max},E_a^{\max},L_a^{\max},D_a^{\max},W_a^{\max},R_a^{\max}),

where NγN_\gamma is incident or detected photon count [photon] as explicitly labeled, EaE_a is energy [J], LaL_a is latency [s], DaD_a is dose in the task-native safety unit, WaW_a is actuator wear [cycle], and RaR_a is risk on a declared scale. Superscript max\max denotes the preregistered ceiling in the same unit, and \preceq means every component is within its ceiling. A scalarized cost may guide a policy only after the component ceilings are enforced and its weights are published.

Multiplex, fusion, and calibration crossover

For route rr in regime ρ\rho, define the protected outcome vector

Yr,ρ=(Ltask,Ucal,Nγ,Elife,Twall,Fsat,Ccross),\mathbf Y_{r,\rho}= (L_{\mathrm{task}},U_{\mathrm{cal}},N_\gamma,E_{\mathrm{life}}, T_{\mathrm{wall}},F_{\mathrm{sat}},C_{\mathrm{cross}}),

where LtaskL_{\mathrm{task}} is task loss in its native unit, UcalU_{\mathrm{cal}} is dimensionless uncertainty-calibration error, NγN_\gamma is photon count [photon], ElifeE_{\mathrm{life}} is lifecycle energy [J], TwallT_{\mathrm{wall}} is wall latency [s], FsatF_{\mathrm{sat}} is saturation fraction [dimensionless], and CcrossC_{\mathrm{cross}} is dimensionless crosstalk. Direct and multiplexed routes are compared componentwise across photon flux, background, detector noise, occupancy, crosstalk, and saturation regimes; no single favorable point establishes an advantage.

For two sensor estimates x^1\hat x_1 and x^2\hat x_2 with errors e1=x^1xe_1=\hat x_1-x and e2=x^2xe_2=\hat x_2-x, retain

P12=E[e1e2],P_{12}=\mathbb E[e_1e_2^{\top}],

where P12P_{12} is cross-covariance in squared native state units and e2e_2^{\top} is transpose. Setting P12=0P_{12}=0 is a tested assumption, not a default. Fused estimates must report marginal covariance, cross-covariance or a justified bound, alignment error in native spatial and temporal units, and failure under common-mode perturbations.

For reference measurement rtrefr_t^{\mathrm{ref}} and calibrated prediction r^ref(ct)\widehat r^{\mathrm{ref}}(c_t) in the same observation unit, define

zt2=(rtrefr^ref(ct))St1(rtrefr^ref(ct)),z_t^2=(r_t^{\mathrm{ref}}-\widehat r^{\mathrm{ref}}(c_t))^{\top} S_t^{-1}(r_t^{\mathrm{ref}}-\widehat r^{\mathrm{ref}}(c_t)),

where StS_t is residual covariance in squared observation units and zt2z_t^2 is dimensionless. The monitor declares a threshold zmax2z_{\max}^2, a window length ww [sample], a false-alarm target [dimensionless probability], and a fallback. Detection delay is seconds [s], recovery time is seconds [s], calibration cost is samples [sample] and joules [J], and pre-detection task loss is reported in the task's native unit. Scene shift, source drift, alignment drift, detector gain, and thermal drift are separate hidden causes.

Optical, digital, and hybrid route accounting

For route rr, end-to-end latency is

Tr=Tintegrate+Tencode+Tpropagate+Tdetect+Tconvert+Ttransfer+Tdigital+Tcontrol,T_r=T_{\mathrm{integrate}}+T_{\mathrm{encode}}+T_{\mathrm{propagate}}+ T_{\mathrm{detect}}+T_{\mathrm{convert}}+T_{\mathrm{transfer}}+ T_{\mathrm{digital}}+T_{\mathrm{control}},

where every term is seconds [s]. Optical propagation time cannot replace TrT_r. The route identity fixes input origin, operator shape, batch, sparsity, effective precision, output dimension, programming frequency, operator reuse, utilization, and duty cycle.

Let the intended transform be u=Wvu=Wv, where vv is an input vector in declared native units, WW is a linear operator with corresponding conversion units, and uu is the desired output vector. Device dd at temperature ϑ\vartheta [K] and age \ell [s] realizes

u^d,ϑ,=(W+ΔWd,ϑ,)v+ϵd,ϑ,analog+ϵd,ϑ,read,\widehat u_{d,\vartheta,\ell}= (W+\Delta W_{d,\vartheta,\ell})v+ \epsilon^{\mathrm{analog}}_{d,\vartheta,\ell}+ \epsilon^{\mathrm{read}}_{d,\vartheta,\ell},

where ΔWd,ϑ,\Delta W_{d,\vartheta,\ell} has the units of WW, ϵanalog\epsilon^{\mathrm{analog}} is analog transform error in output units, and ϵread\epsilon^{\mathrm{read}} is detector, conversion, and readout error in output units. Report bias, covariance, tails, effective precision [bit], and task loss separately across fan-in, depth, device, temperature, age, and workload.

For device population D\mathcal D with NDN_D fabricated devices, fabrication yield is

Yfab=1NDdD1[d meets the preregistered envelope],Y_{\mathrm{fab}}= \frac{1}{N_D}\sum_{d\in\mathcal D} \mathbb 1[d\text{ meets the preregistered envelope}],

where YfabY_{\mathrm{fab}} is dimensionless and NDN_D is device count [device]. The envelope includes transfer-function error, task quality, trimming time [s], tuning energy [J], steady thermal power [W], thermal crosstalk [dimensionless or a declared transfer unit], and stability over the declared interval. Failed dies remain in the denominator.

Query-registered physical compaction

Let encoder hh transform raw observation record OO into retained artifact z=h(O)z=h(O) with size SzS_z [byte]. After the encoder is frozen, query qq from registered set QQ produces answer fq(O)f_q(O) from raw data and reconstructed answer f^q(z)\widehat f_q(z) from the retained artifact. Query recovery is

RQ(z)=1QqQ1 ⁣[dq ⁣(f^q(z),fq(O))ϵq],R_Q(z)=\frac{1}{|Q|}\sum_{q\in Q} \mathbb 1\!\left[ d_q\!\left(\widehat f_q(z),f_q(O)\right)\le\epsilon_q \right],

where Q|Q| is query count [query], dqd_q is error in the native unit of query qq, ϵq\epsilon_q is a tolerance in the same unit, and RQR_Q is dimensionless. The confirmatory evaluator adds sealed future queries and an incident- investigation query after route and retention policies are frozen. Failure is reported as lost query classes, not only as an average score.

Lifecycle, labor, and matched budgets

For one accepted output, complete lifecycle energy is

Elife=Esource+Emodulate+Epropagate+Edetect+EADC+EDAC+Econtrol+Edigital+Ethermal+Efacility+Ecalibrate+Emaintain+Eembodied,\begin{aligned} E_{\mathrm{life}}={}&E_{\mathrm{source}}+E_{\mathrm{modulate}}+ E_{\mathrm{propagate}}+E_{\mathrm{detect}}+E_{\mathrm{ADC}}+E_{\mathrm{DAC}}\\ &+E_{\mathrm{control}}+E_{\mathrm{digital}}+E_{\mathrm{thermal}}+ E_{\mathrm{facility}}+E_{\mathrm{calibrate}}+E_{\mathrm{maintain}}+ E_{\mathrm{embodied}}, \end{aligned}

where every term is joules [J] and analog-to-digital and digital-to-analog conversion are denoted ADC and DAC. Powered propagation elements are charged to EpropagateE_{\mathrm{propagate}}; passive loss appears through increased source or amplifier demand. Amortized embodied energy is

Eembodied=Efabricate+Epackage+EreplaceErecoverNaccepted,life,E_{\mathrm{embodied}}= \frac{E_{\mathrm{fabricate}}+E_{\mathrm{package}}+E_{\mathrm{replace}}- E_{\mathrm{recover}}}{N_{\mathrm{accepted,life}}},

where numerator terms are joules [J] and Naccepted,lifeN_{\mathrm{accepted,life}} is the accepted-output count [output] over measured or conservatively modeled service life. Negative recovery credit must be independently substantiated.

The equal-budget vector for method mm is

Bm=(Ntrain,Nscene,Nγ,Ddose,Nact,Ncal,Nfab,Nquery,Sstate,Twall,Thuman,Elife),\mathbf B_m=(N_{\mathrm{train}},N_{\mathrm{scene}},N_\gamma,D_{\mathrm{dose}}, N_{\mathrm{act}},N_{\mathrm{cal}},N_{\mathrm{fab}},N_{\mathrm{query}}, S_{\mathrm{state}},T_{\mathrm{wall}},T_{\mathrm{human}},E_{\mathrm{life}}),

where NtrainN_{\mathrm{train}} is training-example count [example], NsceneN_{\mathrm{scene}} is physical-scene count [scene], NγN_\gamma is photon count [photon], DdoseD_{\mathrm{dose}} is dose in a declared task-native unit, NactN_{\mathrm{act}} is actuator-cycle count [cycle], NcalN_{\mathrm{cal}} is calibration-sample count [sample], NfabN_{\mathrm{fab}} is fabricated-device count [device], NqueryN_{\mathrm{query}} is evaluator-query count [query], SstateS_{\mathrm{state}} is retained and working storage [byte], TwallT_{\mathrm{wall}} is wall time [s], ThumanT_{\mathrm{human}} is labor [person-hour], and ElifeE_{\mathrm{life}} is lifecycle energy [J]. Human design, alignment, labeling, tuning, calibration, inspection, safety review, maintenance, and incident response remain role-stratified entries.

An arm is budget-matched only when every binding component is within its preregistered tolerance or when the comparison is explicitly a Pareto frontier. Freed budget from an ablation remains unused.

Protected outcome vector and retirement estimand

Keep the following outcome families separate:

Zm=(Zaperture,Zphoton,Zphase,Zprior,Zdrift,Zrange,Zfusion,Ztransform,Zconversion,Zanalog,Zfabrication,Zthermal,Zsafety,Zlifecycle),\mathbf Z_m=(Z_{\mathrm{aperture}},Z_{\mathrm{photon}},Z_{\mathrm{phase}}, Z_{\mathrm{prior}},Z_{\mathrm{drift}},Z_{\mathrm{range}},Z_{\mathrm{fusion}}, Z_{\mathrm{transform}},Z_{\mathrm{conversion}},Z_{\mathrm{analog}}, Z_{\mathrm{fabrication}},Z_{\mathrm{thermal}},Z_{\mathrm{safety}}, Z_{\mathrm{lifecycle}}),

where the components respectively contain aperture/diffraction results, photon/shot-noise results, phase-ambiguity results, prior-mismatch results, calibration/drift results, saturation/dynamic-range results, fusion-covariance results, optical-transform results, conversion/readout results, analog-error results, fabrication results, thermal-work results, active-illumination safety results, and full lifecycle energy and labor results. Each ZZ is a structured record with the native units defined above; the vector is not scalarized for acceptance.

For protected component jj, the paired effect of mechanism kk is

Δk,j=Zj(mfull)Zj(mk),\Delta_{k,j}=Z_j(m_{\mathrm{full}})-Z_j(m_{-k}),

where mkm_{-k} removes only mechanism kk and receives no replacement resource. Report Δk,j\Delta_{k,j} in the native unit of outcome jj with a 95% uncertainty interval over independent scene, operator, regime, device, site, and seed strata.

Retire the residual when an equal-budget mature-null composition matches the protected vector; when a gain depends on hidden prior, operator, calibration, or query leakage; when it fails held-out devices or regimes; or when no preregistered ablation isolates value beyond the registered inverse, uncertainty, sensing, control, acceleration, optical, hybrid-design, and calibration baselines. The editable source for the system diagram is operator-qualified-physical-inference.mmd.