Generalized Quantum Photodetection

• Generalized quantum photodetection is a comprehensive framework that extends standard photodetection by mapping optical field states to outcomes via POVMs and quantum instruments.
• It incorporates realistic imperfections such as inefficiencies, dark counts, and incomplete click data, allowing controlled estimation within a defined operator basis.
• The approach is applied to full-chain detection models, quantum illumination, and advanced interference experiments, providing unified methods for observable reconstruction and validation.

Generalized quantum photodetection framework denotes a family of formalisms that extends standard electric-field, normal-order photodetection into a broader measurement theory in which detector response is represented by positive operator-valued measures (POVMs), quantum instruments, scattering networks, or fully coupled light–matter–amplification dynamics. Across this literature, photodetection is not treated as a primitive count observable alone, but as a device-dependent map from optical field states to outcomes, conditioned states, and inferred observables. The resulting frameworks accommodate informationally incomplete click detectors, realistic inefficiencies and dark counts, multimode Gaussian optics, ultrastrong light–matter coupling, and engineered detector responses involving electric and magnetic field amplitudes (Kovalenko et al., 2018, Young et al., 2018, Stefano et al., 2017, Hatifi et al., 16 May 2026).

1. Measurement geometry and generalized observables

A central operator-theoretic formulation represents a measurement by a POVM

$$\{\hat{\Pi}_n:n\in\mathcal I\},\qquad \hat{\Pi}_n\ge 0,\qquad \sum_{n\in\mathcal I}\hat{\Pi}_n=\hat 1,$$

with Born probabilities

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$

In the geometric construction developed for photocounting, this probability is interpreted as a Hilbert–Schmidt scalar product,

$$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$

so the POVM elements play the role of basis vectors for a measurement-defined subspace of observables. An observable $\hat B$ is expanded as

$\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$

where $\hat R$ is orthogonal to the POVM basis and quantifies the systematic mismatch introduced by incomplete detection. Expectation values then take the form

$$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$

This establishes generalized photodetection as a coordinate problem in operator space rather than a direct identification of raw counts with ideal photon-number statistics (Kovalenko et al., 2018).

The dual geometric object is the contravariant operator-valued measure (COVM) $\{\hat\Pi^n\}$, defined by

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$

With the POVM metric tensor

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$

covariant coordinates correspond to measured probabilities, whereas contravariant coordinates specify the generalized observable reconstructed from those probabilities. For array detectors with finitely many on–off channels, the POVM is finite-dimensional in an infinite-dimensional Fock space and is therefore informationally incomplete. In that setting, the orthogonal remainder obeys

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$0

with the bound

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$1

The framework therefore does not identify finite click data with full state information, but with controlled estimation inside the POVM span (Kovalenko et al., 2018).

A related line of work formulates imperfect photodetection explicitly as a nonorthogonal POVM rather than as an ideal projection followed by ad hoc corrections. In a cavity-QED model with a two-level atom-pointer and an imperfect ionization-chamber readout, the outcomes $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$2 are described by

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$3

with $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$4. In that description, detector inefficiency and wrong-state identification are built directly into the measurement operators and the induced field-space conditional maps (Trifanov et al., 2012).

2. End-to-end detector models and retrodictive POVMs

A major strand of the subject models the detector as a complete quantum chain comprising transmission into the device, internal conversion, amplification, and final readout. In that program, the basic outcome object is again a POVM on the external optical Hilbert space,

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$5

but the POVM is derived from a microscopic detection architecture rather than postulated. A particularly explicit result is the retrodictive wavepacket projector generated by a time-dependent two-level detector,

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$6

where

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$7

and

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$8

with $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$9. By choosing $$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$0 and $$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$1, the detector can project onto an arbitrary smooth single-photon wavepacket, and in the realistic three-stage model the external click POVM assumes the compact form

$$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$2

Arbitrary wavepacket projection remains possible provided the transmission function has no zeros on the relevant support, with matching condition

$$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$3

These results recast photodetection as a retrodictive mode-selective measurement whose temporal and spectral structure is fixed by the detector’s own resonance and filtering properties (Propp et al., 2020, Propp, 2022).

Before amplification, the same logic can be expressed as a two-port quantum scattering problem. In that setting, the first irreversible stage of single-photon detection is summarized by a complex transmission amplitude $$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$4. The corresponding long-time click POVM is

$$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$5

so $$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$6 is the maximum detection probability for a monochromatic component at frequency $$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$7. Two additional input-independent performance measures are the spectral bandwidth

$$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$8

and the group delay

$$\langle \hat\varrho,\hat\Pi_n\rangle_{\rm HS}={\rm Tr}(\hat\varrho\hat\Pi_n),$$9

For a single discrete-state detector network,

$\hat B$0

with perfect transmission at resonance when $\hat B$1. This network description is then extended to parallel, series, and hybrid topologies, all encoded by the structure of $\hat B$2 (Propp et al., 2019).

A Heisenberg-picture realization of the full chain introduces one Hamiltonian for photon generation, propagation to an absorber, absorption by a $\hat B$3-type molecule, irreversible transition, and amplification into a separate output field. The output operator can then be written schematically as signal plus noise, and the steady-state amplified output obeys a nonlinear relation of the form

$\hat B$4

with

$\hat B$5

at resonance in the prototype model. In this framework, amplification is not appended after absorption; it is triggered by an internal state change of the detector and can occur at a wavelength different from that of the absorbed photon (Biswas et al., 2020).

A more abstract open-system formulation treats the incoming field, absorption dynamics, and amplification dynamics as one coupled quantum system. The matter sector evolves through auxiliary density operators $\hat B$6 indexed by field Fock sectors, with internal Liouvillian, light–matter coupling, and amplification terms all included in one hierarchy. Ideal performance is defined by reproduction of the arrival-time statistics,

$\hat B$7

with efficiency $\hat B$8, zero latency $\hat B$9, and jitter equal to the intrinsic pulse-width contribution. This defines an explicit mathematical notion of ideal photodetector performance inside the generalized framework (Young et al., 2018).

3. Imperfections, multiplexing, and realistic count models

Generalized photodetection frameworks are especially prominent when the hardware does not deliver exact photon-number projections. A canonical example is spatial multiplexing with standard on–off detectors. In the integrated divide-and-conquer architecture, a single field is distributed approximately uniformly over $\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$0 channels and the observable is the number $\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$1 of clicks. The exact click-counting model is

$\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$2

For a coherent state, the click statistics are binomial, whereas nonclassical inputs produce deviations quantified by

$\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$3

This formalism replaces proportional photon-to-electron conversion by an exact operator model of multiplexed threshold detection (Heilmann et al., 2015).

In generalized Hong–Ou–Mandel interferometry and boson-sampling-type experiments, realistic photon-number resolving detectors are modeled by Fock-diagonal POVMs

$\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$4

so the measured multimode distribution is

$\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$5

Under the no-false-count assumption $\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$6 for $\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$7 and postselection on $\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$8, the realistic probabilities collapse to a multiplicative correction,

$\hat B=\sum_{n\in\mathcal I} B^n\hat\Pi_n+\hat R,$9

The paper derives explicit $\hat R$0 factors for pure loss, arrays of on–off detectors, and detectors with finite dead time. In monochromatic dead-time detection, for example,

$\hat R$1

This preserves access to ideal permanental interference probabilities after postselection, but only after detector response is folded into the model (Len et al., 2021).

Imperfect readout can also enter through state misidentification rather than limited number resolution. In the atom-pointer model, the empirical coefficients $\hat R$2, $\hat R$3, $\hat R$4, and $\hat R$5 describe false assignment and missed-detection channels, and the measurement backaction is encoded by outcome-conditioned quantum operations

$\hat R$6

The induced field evolution therefore depends on the actual detector error model, not merely on a scalar efficiency parameter (Trifanov et al., 2012).

4. Reconstruction, conditional dynamics, and inference

A distinctive application of generalized photodetection is the direct estimation of observables that are not themselves measured by the detector. Given a detector POVM $\hat R$7, one computes covariant coordinates

$\hat R$8

forms the dual metric $\hat R$9, and obtains contravariant coefficients

$$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$0

Expectation values are then estimated from outcome probabilities as

$$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$1

with the orthogonal remainder $$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$2 providing a systematic-error estimate. The construction is applied to photon-number moments $$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$3, normal-ordered moments $$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$4, generating functions such as $$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$5 and $$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$6, and more general functions $$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$7. The same framework introduces a generalized star prescription,

$$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$8

which replaces ordinary functional calculus by a rule tied to the actual generalized measurement outcomes (Kovalenko et al., 2018).

The geometric method also yields detector-aware state reconstruction in unbalanced homodyne detection. The Cahill–Glauber distribution is written as

$$\langle \hat B\rangle=\sum_{n\in\mathcal I} B^n\varrho_n+\langle \hat R\rangle.$$9

where

$\{\hat\Pi^n\}$0

For array detectors, the coefficients

$\{\hat\Pi^n\}$1

enter the click-basis representation used for reconstruction. The operator norm is Hilbert–Schmidt only for $\{\hat\Pi^n\}$2,

$\{\hat\Pi^n\}$3

so for $\{\hat\Pi^n\}$4 finite photon-number truncation is used to obtain a well-defined approximation and an error estimate (Kovalenko et al., 2018).

Sequential photodetection can likewise be formulated as a joint generalized measurement process rather than exclusively through the quantum regression theorem. For two measurements at $\{\hat\Pi^n\}$5,

$\{\hat\Pi^n\}$6

with $\{\hat\Pi^n\}$7. Photon counting appears as a quantum-jump POVM,

$\{\hat\Pi^n\}$8

and homodyne amplitude readout is represented by a Gaussian weak-measurement POVM. For reversed temporal ordering, the formalism introduces the past quantum state effect matrix $\{\hat\Pi^n\}$9, evolving under the adjoint master equation, and the smoothing formula

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$0

This places intensity–intensity and intensity–amplitude correlations inside a unified conditioned-dynamics framework (Xu et al., 2015).

5. Beyond electric-only Glauber response

Several generalized frameworks explicitly depart from the standard assumption that photodetection is governed by bare electric-field normal ordering. In arbitrarily strong light–matter coupling, the detector couples weakly to the interacting system through

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$1

while the system itself is diagonalized in dressed eigenstates ${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$2. The correct positive-frequency operator is then defined in the interacting basis,

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$3

and the detection rate becomes

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$4

For a broadband electric-field detector,

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$5

This dressed-basis construction resolves the ultrastrong-coupling failure of bare normal-order formulas such as ${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$6. In the quantum Rabi model, the ground state has nonzero bare-photon population, yet

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$7

so the ground-state bare excitations are not detected as real photons (Stefano et al., 2017).

A different extension keeps the detector microscopic and shows that vacuum modes act as an additional reservoir coupled to the detector’s induced polarization current. After a Markov reduction, the detector amplitude obeys

${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$8

where ${\rm Tr}(\hat\Pi^n\hat\Pi_m)=\delta^n_m.$9 is decay into the electronic reservoir and $$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$0 is decay into the radiative or vacuum reservoir. The mean photocurrent for an incident coherent state is

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$1

and at resonance

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$2

Current fluctuations satisfy

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$3

and in the Purcell regime

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$4

Vacuum modes therefore modify both quantum efficiency and shot noise through detector dynamics rather than through direct signal counts (Wadood et al., 2018).

A more recent generalization replaces the electric-only Glauber response by a coherent superposition of electric and magnetic field amplitudes,

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$5

The detection probability is

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$6

In a far-field two-source geometry, complete detector-amplitude cancellation occurs at $$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$7. In a two-mode single-photon setting, the detector measures

$$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$8

so $$g_{nm}={\rm Tr}(\hat\Pi_n\hat\Pi_m),\qquad \hat\Pi^n=\sum_m g^{nm}\hat\Pi_m,$$9 continuously rotates the effective measurement basis, with exact first-order visibility

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$00

In the lossy resonant realization, a monitored radiative channel can be dark when $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$01, while the absorption

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$02

reaches

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$03

at critical coupling $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$04. This generalization interprets engineered electric–magnetic interference as a generalized measurement operator rather than only as a scattering effect (Hatifi et al., 16 May 2026).

6. Protocol-level applications: illumination, interference, and tomography

The generalized framework has been applied directly to quantum illumination with simple Geiger-mode photodetection. In the Gaussian-state treatment, the detector is described by click/no-click POVMs,

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$05

For a Gaussian state, the perfect-detector no-click probability is computed from the overlap with vacuum in phase space, and imperfect detection is reduced to loss plus a perfect detector. Heralding an idler click from a TMSV source produces a vacuum-suppressed thermal state (VST), while idler no-click produces a photon-number suppressed thermal state (PNST). The framework yields target-detection click probabilities, single-shot posterior probabilities, and Monte Carlo repeated-shot Bayesian updates, and it shows that heralded TMSV illumination can outperform coherent-state illumination in low-energy, low-reflectivity, high-background regimes. It also identifies click-probability matching, rather than equal mean photon number, as a useful comparison criterion (Yang et al., 2020).

A detector-level microwave quantum-illumination analysis studies the experimentally relevant parametric-mixer receiver under finite efficiency, dark counts, and limited photon-number resolution. The return and idler are mixed and counted at two outputs with observables

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$06

and correlated photon counting uses

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$07

The covariance terms

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$08

and

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$09

drive the CPC enhancement. The two outputs are strongly asymmetric: $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$10 and remains near binary, so $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$11 is close to optimal for PC1, while $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$12 is background dominated and remains worse than the classical reference even as $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$13. The paper concludes that PC1 is the practical workhorse, CPC yields only modest improvement with strong calibration sensitivity, and dark counts are more damaging than moderate inefficiency (Kronowetter et al., 2023).

Realistic interference experiments with spectrally multimode Gaussian sources are likewise treated at the detector level. For a subset $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$14 of modes with reduced covariance matrix $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$15, the off probability is

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$16

threshold probabilities are obtained by inclusion–exclusion, and exact number-resolving probabilities follow from

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$17

Applied to heralded Hong–Ou–Mandel interference, this formalism separates spectral impurity from photon-number impurity and shows that increasing source brightness decreases visibility for any degree of spectral impurity, even with number-resolving detectors. Tight filtering can raise the visibility to about $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$18 at low power for the nonseparable source studied, but reduces the maximum heralding rate from about $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$19 to $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$20 and the heralding efficiency from $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$21 to about $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$22 at the rate optimum (Thomas et al., 2020).

In the generalized multimode Hong–Ou–Mandel setting with realistic PNR detectors, the ideal permanent-based probabilities

$\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$23

are convolved with detector response coefficients $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$24. Under postselection on the total photon number, the detector-smeared probabilities are proportional to the ideal ones, which preserves access to the bosonic-interference structure needed in validation and nonclassicality witnessing (Len et al., 2021).

7. Scope, limitations, and unresolved problems

The generalized framework does not define one unique formalism. Rather, the literature develops several complementary generalizations: operator-space geometry for incomplete POVMs, full-chain POVM engineering for mode-selective single-photon detection, open-system models coupling absorption and amplification, detector-response convolutions for realistic counting, and non-Glauber measurement operators for ultrastrong coupling or electric–magnetic selectivity. A plausible implication is that the common object across these approaches is not a universal count operator, but a detector-specific mapping from field states to probabilities, conditioned states, or inferred observables (Kovalenko et al., 2018, Young et al., 2018, Propp et al., 2019).

Several limits recur across the subject. Finite click-detector arrays are informationally incomplete and cannot reconstruct arbitrary quantum states exactly; the inaccessible component is represented by the orthogonal remainder $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$25. In network-based detection, if $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$26 at relevant frequencies, those spectral components are irretrievably lost and no downstream processing can recover them. In shelving-type and band detectors, reset, latency, and jitter trade off against long-lived monitored populations and count rate. In parametric-mixer quantum illumination, the advantage is highly sensitive to the weights $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$27, to dark counts, and to the requirement that the mixer and PC2 remain linear and unsaturated at large thermal background. These are device-architectural constraints rather than merely statistical nuisances (Propp et al., 2020, Young et al., 2018, Kronowetter et al., 2023, Kovalenko et al., 2018).

At the most speculative end of the literature, late-time photodetection has been proposed as the basis for a realist completion of relativistic quantum theory. In that construction, hypothetical asymptotic electromagnetic-field detections on a late hypersurface are used to define local beables at spacetime point $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$28, conditioned only on outcomes outside the future light cone of $\varrho_n={\rm Tr}(\hat\varrho\,\hat\Pi_n).$29. The scheme is explicitly idealized and the author identifies several open problems: lack of rigorous QFT foundations, the infinite soft-photon problem, detector idealization, Lorentz-invariance issues tied to preferred frames or cutoffs, dependence on detector geometry and thresholds, and the absence of a fully satisfactory theory of approximately localized measurements in QFT. This suggests that generalized quantum photodetection has become not only a detector theory, but also a testing ground for questions about inference, realism, and relativistic measurement (Kent, 2016).