---
title: Photon-Number Resolution in Quantum Detectors
url: https://www.emergentmind.com/topics/photon-number-resolution-pnr
type: topic
---

# Photon-Number Resolution in Quantum Detectors

Photon-number resolution (PNR) is the ability of a detector to discriminate how many photons were present in an optical pulse or mode, rather than only reporting whether at least one photon was detected. In ideal form, it corresponds to measurement outcomes \(n=0,1,2,\dots\) with POVM elements \(|n\rangle\langle n|\); in practice, it is realized by intrinsically energy-resolving detectors, by spatial or temporal multiplexing of ON-OFF detectors, or by intrinsic multi-photon discrimination in superconducting nanowire devices. Because many protocols rely on post-selection conditioned on specific photon-number outcomes or on precise characterization of photon-number statistics, PNR is central to entanglement swapping, teleportation, linear-optical quantum computing, Gaussian boson sampling, heralded state preparation, quantum-enhanced metrology, source characterization, and quantum communications [2310.12472] [2411.15360] [2507.05554].

## 1. Operational meaning and measurement outcomes

In operational terms, PNR means mapping a detection event to an integer \(n\) or to an estimate \(\hat n\) that corresponds to the number of absorbed photons. Realistic devices do not perfectly resolve large \(n\); they have a maximal resolvable number \(N_{\max}\) with finite confidence, and different platforms realize this mapping through different observables [1812.05422].

One broad class is continuous-output PNR. Transition-edge sensors are described as phase-insensitive PNR detectors characterized by POVM elements diagonal in the Fock basis,
\[
\Pi_n=\sum_{m=0}^{\infty}\theta_{nm}\,|m\rangle\langle m|,
\]
where \(\theta_{nm}\) is the conditional probability of reporting \(n\) given \(m\) incident photons [2411.15360]. A related continuous-outcome picture appears in intrinsic SNSPD PNR, where the measured variable can be arrival time, amplitude, slope, or a multidimensional feature vector derived from the pulse trace. The detector then implements a classification rule on regions of that continuous outcome space rather than directly returning a digitized energy [2505.09271].

A second broad class is multiplexed PNR. Here, the detector output is inferred by aggregating clicks from multiple ON-OFF detectors within a coincidence window or across many spatial, temporal, or frequency channels. In that setting, the output is not an intrinsic energy measurement; it is a statistical inference based on click patterns, coincidence windows, or binned time tags [2110.04245] [2507.05554].

A third class comprises architectures that return only partially resolved outcomes. In arrival-time-based SNSPD PNR, events earlier than a final boundary can be grouped as “4+” when individual classes are no longer resolvable, so the outcome alphabet need not coincide with exact photon number for all \(n\) [2510.00714]. This suggests that encyclopedic use of “PNR” includes exact low-\(n\) resolution, truncated resolution, and calibrated coarse-grainings, provided the detector output remains correlated with photon number in a single shot.

## 2. Statistical models and performance criteria

The starting point for most PNR analyses is the coherent-state distribution
\[
P(n)=e^{-\mu}\frac{\mu^n}{n!},
\]
with mean photon number \(\mu\). With overall detection efficiency \(\eta\), detected counts follow a thinned distribution,
\[
P_{\mathrm{det}}(k)=\sum_{n=k}^{\infty}\binom{n}{k}\eta^k(1-\eta)^{n-k}P(n),
\]
which reduces to a Poisson law with mean \(\eta\mu\) for coherent states [2310.12472]. In detector-characterization experiments this relation is used to compare measured photon-number histograms against calibrated Poisson predictions, and in tomography it provides the known input ensemble against which a detector response matrix is reconstructed [2411.15360] [2504.02202].

A general forward model writes the measured distribution as
\[
Q(n)=\sum_m P(n\mid m)\,S(m),
\]
where \(S(m)\) is the incident photon-number distribution and \(P(n\mid m)\) is the detector response matrix or confusion matrix. This form appears in superconducting parallel nanowire detectors, in multiplexed click-detector arrays, and in amplitude-coded parallel SNSPDs [0902.4824] [1812.05422] [2207.14538]. In timestamp-based intrinsic PNR, an analogous likelihood is written over continuous kernels \(p(t\mid n)\),
\[
\mathcal{L}=\sum_i \ln\!\left[\sum_{n=0}^{N_{\max}}P(n)\,p(t_i\mid n)\right],
\]
and practical implementations often replace full maximum-likelihood estimation by decision regions and overlap integrals of fitted distributions [2310.12472].

Several non-equivalent figures of merit coexist. For multiplexed ON-OFF architectures, Jönsson and Björk define the PNR quality \(Q_n\) as the worst-case probability that the detector correctly predicts the input photon number, and show that \(Q_n\) is monotone non-increasing in \(n\) [1812.05422]. For interleaved SNSPD arrays, the photon-number-resolution SNR is defined as
\[
\mathrm{SNR}_{\mathrm{PNR}(n)}=\frac{H(n)-H(n-1)}{V_n},
\]
where \(H(n)\) is the \(n\)-peak position and \(V_n\) is its FWHM [2312.17589]. For intrinsic SNSPD timing methods, a resolvability criterion is proposed:
\[
(\mu_n-\mu_{n+1})\geq \mathrm{FWHM}_n,
\]
with a Gaussian-limit form
\[
(\mu_n-\mu_{n+1})\geq 2\sqrt{2\ln 2}\,\sigma_{\mathrm{tot},n}.
\]
This criterion is explicitly tied to the continuous distributions typical of latency- or amplitude-based SNSPD PNR [2505.09271] [2510.00714].

The shape of the fitted response function is itself part of the performance model. In one timing-based SNSPD study, Voigt profiles were used to model projected timing histograms, separating a Gaussian instrument response from Lorentzian broadening [2310.12472]. In another, exponentially-modified Gaussian (EMG) fits were used because Gaussian-only models underestimate the long tail and hence underestimate misidentification probabilities; for \(N=1\), for example, \(p_{\mathrm{misidentified}}=0.14\%\) with EMG versus \(0.0032\%\) with Gaussian fits [2510.00714]. This suggests that “PNR performance” is inseparable from the statistical model used to define overlap, threshold placement, and reconstructed POVMs.

## 3. Detector architectures and implementation classes

Conventional routes to PNR include transition-edge sensors, multiplexing with threshold detectors, and parallel or segmented SNSPD architectures [2310.12472]. These routes differ mainly in the physical observable that carries photon-number information: absorbed energy, click occupancy, amplitude step height, or latency shift.

Transition-edge sensors remain the canonical intrinsically energy-resolving architecture. They operate as ultralow-noise microcalorimeters with near-unity efficiency and high energy resolution, but their thermal recovery tail persists for microseconds, so the sensor should ideally start from the same baseline for each event. A recent machine-learning study showed that overlapping TES traces can still be classified at 500–800 kHz, extending accurate TES PNR operation to 800 kHz while maintaining accurate photon-number assignment up to at least five photons [2411.15360].

Multiplexed ON-OFF detectors realize approximate PNR by distributing photons over many bins. Spatial arrays, temporal arrays, loop architectures, and adaptive storage loops all belong to this class [1812.05422] [2009.08334] [2311.13515]. Intrinsic SNSPD PNR instead exploits measurable changes in a single detector’s response time, rising-edge slope, or amplitude when more than one photon is absorbed within a short time window [2505.09271].

| Architecture | Physical basis | Representative result |
|---|---|---|
| TES with machine-learning processing | Pulse amplitude, area, and full-trace classification | 800 kHz operation, at least a four-fold improvement, accurate photon-number assignment up to at least five photons [2411.15360] |
| Temporal multiplexing with two SNSPDs | 16 time bins measured with binary detectors | Accurate prediction of photon numbers between \(10^{-3}\) to \(10^{2}\); effective quantum efficiency \(49\%\) [2009.08334] |
| Parallel/segmented SNSPD architectures | Summed amplitude from multiple pixels | 4-pixel P-SNSPD with SDE \(92.5\%\) and \(K=0\)–4; segmented twin-layer SNSPD with system detection efficiency of \(\sim 98\%\) and photon-number resolution of 32 [2207.14538] [2504.02202] |
| Up-conversion plus SiPM | Telecom-to-visible conversion followed by multi-pixel APD readout | Overall detection efficiency \(\eta_{\mathrm{tot}}\approx 4\%\); Poissonian statistics maintained up to approximately 20 simultaneous detections [1003.5579] |
| Adaptive storage loop with a single click detector | Tunable outcoupling and Bayesian/adaptive inference | Dynamic range extended by up to an order of magnitude relative to a purely passive setup [2311.13515] |

Parallel-nanowire implementations are historically important because they provided intrinsic amplitude-coded PNR with a single electrical output proportional to photon number. The Parallel Nanowire Detector used subwavelength-scale spatial multiplexing, achieved pulse widths as low as 660 ps FWHM, showed counting performance at 80 MHz repetition rate, and exhibited no observable multiplication-noise buildup [0807.0526] [0902.4824]. A later parallel SNSPD architecture in MoSi demonstrated a system detection efficiency of \(92.5\%\), \(42\) ps FWHM jitter at the single-photon level, recovery to more than \(50\%\) efficiency after \(10\) ns, and full recovery in about \(40\) ns [2207.14538].

Array architectures can also improve the electronic separability of photon-number levels. In a 16-element interleaved SNSPD array, a level comparator circuit converted noisy analog pulse heights into digital steps, improved the PNR SNR by more than a factor of four, and reduced system timing jitter from about 90 ps to about 72 ps [2312.17589]. At the opposite end of the complexity spectrum, a room-temperature telecom detector based on up-conversion and a SiPM achieved only moderate efficiency, but preserved coherent-state Poisson statistics up to approximately 20 simultaneous detections and operated without cryogenic nanowire hardware [1003.5579].

## 4. Intrinsic SNSPD PNR: timing, amplitude, and electrothermal dynamics

Intrinsic SNSPD PNR relies on the fact that multi-photon absorption modifies the electrothermal turn-on dynamics of the nanowire. Absorption of one or more photons creates one or more resistive hotspots, triggering redistribution of the bias current and formation of a normal domain. Multi-photon absorption accelerates hotspot growth and enlarges the initial normal domain; as modeled in prior work on turn-on dynamics, the rising slope increases with photon number, so multi-photon events cross a fixed electrical threshold earlier. The same body of work emphasizes that the rising-edge separation between adjacent photon-number levels decreases roughly as \(1/\sqrt{n}\) [2310.12472] [2505.09271].

Several readout strategies exploit that physics. A differentiated slew-rate readout converts the \(k\)-dependent rising-edge slope into amplitude. Impedance-matching tapers make the analog pulse height itself strongly sensitive to the number of simultaneous hotspots. Waveform pattern matching, linear fits to the rising-edge slope, two-threshold timing, and PCA on full traces all extract related information from the same electrothermal event [1911.09485] [2505.09271]. These methods differ in hardware overhead, but they share a common dependence on low jitter, short optical pulses, and a stable mapping between hotspot dynamics and the measured trace.

A key demonstration used a single conventional SNSPD in the telecom C-band and ultra-high-resolution time-tagging of both the rising and falling edges of the detector pulse. At \(1554\) nm, with a Single Quantum EOS Series detector specified at \(86\%\) system detection efficiency and \(8.1\) ps RMS jitter, plus a time-tagger contribution of \(1.3\) ps RMS per channel, the combined time-difference jitter was about \(8.3\) ps RMS. Distinct clusters in a two-dimensional histogram of rising versus falling edge time differences were visible up to at least \(N=5\); classification by optimal projection and Voigt-profile fitting showed that using both edges significantly reduced cross-talk compared to using only the rising edge. The same detector was used to measure photon-number statistics of coherent light and joint photon-number distributions of non-classical states in a type-II SPDC experiment [2310.12472].

A different intrinsic route uses amplitude-sensitive nanowire engineering. The superconducting tapered nanowire detector resolved up to five absorbed photons at \(1550\) nm, had \(16.1\) ps timing jitter, \(<2\) c.p.s. device dark count rate, \(\sim 86\) ns reset time, and \(5.6\%\) system detection efficiency without cavity enhancement [1911.09485]. A segmented twin-layer NbN/SiO\(_2\)/NbN SNSPD on a dielectric mirror pushed this architecture further: it exhibited a single-photon system detection efficiency of about \(98\%\), photon-number resolution capability up to 32, count rate about \(41.7\) MHz at \(50\%\) of the maximum SDE, and timing jitter as low as \(40.4\) ps for \(n=32\) [2504.02202].

Intrinsic timing PNR is strongly affected by optical pulse shape and duration. A practical study using a commercial SNSPD at \(1550\) nm found that Gaussian spectral filtering yields much cleaner arrival-time histograms than a bandpass filter of equal bandwidth, that \(T_p\approx 2.9\) ps gives clean separation, that separation begins to degrade at \(T_p\approx 25\) ps, and that at \(T_p\approx 60\) ps the \(N=1\) and \(N=2\) peaks overlap strongly. In that system, EMG fits described the data over four orders of magnitude in count rate, reconstructed POVMs showed sharp photon-number-diagonal features, and the resolvability criterion yielded \(N_{\max}=3\) [2510.00714]. By contrast, NbTiN SNSPDs optimized for \(850\)–\(950\) nm demonstrated over \(94\%\) system detection efficiency at \(940\) nm, sub-11 ps timing jitter for one photon, sub-7 ps for two photons, and resolution up to 7 photons using conventional cryogenic readout circuitry [2401.07265].

Kinetic inductance introduces a direct speed–resolution trade-off. In waveguide-integrated SNSPDs, increasing the nanowire kinetic inductance from \(165\) nH to \(872\) nH improved PNR quality by \(12\%\), \(31\%\), and \(23\%\) over the first three photon numbers, because larger \(L_k\) increased the latency separation between photon-number peaks. The same change reduced the detector’s count rate from \(165\) Mcps to \(19\) Mcps because the recovery time increased linearly with \(L_k\) [2410.23162]. This suggests that intrinsic SNSPD PNR is best understood as a circuit-engineering problem coupled to hotspot physics rather than as a purely material property.

## 5. Multiplexed PNR, scaling laws, and adaptive schemes

Multiplexed PNR begins from a binary detector model and distributes the input across many modes so that the total number of clicks approximates photon number. For a balanced array of \(M\) ideal ON-OFF detectors, the probability of obtaining exactly \(k\) clicks given \(n\) photons is
\[
P(K=k\mid n)=\frac{(M)_k}{M^n}S(n,k),
\]
where \(S(n,k)\) are Stirling numbers of the second kind and \((M)_k=M(M-1)\cdots(M-k+1)\) [2507.05554]. In a temporal array measured with two SNSPDs, coherent light with mean \(\mu\) per pulse gives a single-bin click probability
\[
p(\mu)=1-(1-p_d)e^{-\eta\mu/M},
\]
and the maximum-likelihood estimate of \(\mu\) follows by inverting the binomial mean across \(M\) bins [2009.08334].

These models make the scaling limits explicit. Jönsson and Björk showed that the required quantum efficiency is very high in order to achieve even moderate photon resolution with multiplexed click detectors, and that the number of click detectors required grows quadratically with the maximal number of photons resolvable. For \(\eta=1\), their simulations gave \(M\ge 20\) elements to resolve 5 photons with \(Q_n\ge 0.5\), and \(M\ge 72\) elements to resolve 10 photons with \(Q_n\ge 0.5\) [1812.05422]. A more recent theoretical framework proves that, for multiplexed PNR detectors, the estimation error in terms of photon-number moments decreases inverse proportionally to the number of detectors, i.e. as \(\mathcal{O}(1/M)\), under the paper’s assumptions [2507.05554].

Experimentally, multiplexing is valuable even when single-shot fidelity is limited. A 16-element temporal array built from fiber couplers and two SNSPDs accurately predicted photon numbers between \(10^{-3}\) and \(10^2\) for the same number of input pulses used to estimate the click statistics, but its effective quantum efficiency of \(49\%\) prevented high-precision single-shot PNR [2009.08334]. Adaptive storage-loop PNR addresses a related limitation by changing the outcoupling fraction round by round. In that architecture, the click probability in bin \(i\) is
\[
P_{\mathrm{click},i}(N_i)=1-(1-d_i)(1-\eta_i T_i\gamma_i)^{N_i},
\]
and the control law chooses \(T_i\) to keep each bin informative. Simulations showed that adaptive feedback can extend dynamic range by up to an order of magnitude relative to a passive setup and can produce sub-shot-noise photon-number variance under favorable \(\eta\) and \(\gamma\) [2311.13515].

Multiplexing also extends beyond beamsplitters and delay lines. A proposed PNR detector based on a cascade of waveguide-coupled \(\Lambda\)-type emitters uses deterministic single-photon subtraction followed by ON-OFF detection. For a two-photon subtraction case, 7 dB of squeezing, and an array of 20 detectors of efficiency \(95\%\), the calculation predicts fidelity \(\sim 0.88\) with success probability \(\sim 3.8\%\), and the paper argues that megahertz-rate cat-state generation is achievable using an on-chip array of tens of ON-OFF detectors [2507.09034] [2507.05554].

## 6. Applications, interpretive issues, and outlook

Across architectures, PNR is used both as a detector primitive and as a characterization tool. In coherent-state characterization it measures truncated Poisson distributions and detector POVMs; in non-classical-state experiments it gives joint photon-number distributions \(P(n_s,n_i)\), conditional moments, covariance, and bunching or anti-bunching signatures [2310.12472] [2411.15360]. Timing-based SNSPD PNR has already been used to observe photon-pair bunching at a single output port of a Hong–Ou–Mandel interferometer, and segmented or parallel architectures are explicitly motivated by linear-optical quantum computing, Gaussian boson sampling, quantum communications, heralded source characterization, \(g^{(2)}\) measurements, and precision metrology [1911.09485] [2504.02202].

Calibration remains an open technical issue. Coherent states are a robust reference, but coherent-state-based characterization cannot by itself reveal absolute loss; accurate attenuator calibration or a reference threshold detector is required for absolute detector tomography [2310.12472]. Timing-based intrinsic PNR additionally requires a precise optical timing reference, re-optimization of clustering boundaries to accommodate slow drifts, and a statistical model that captures non-Gaussian broadening [2310.12472] [2510.00714]. In TES systems, overlap-aware supervised learning and unsupervised clustering reduce the penalty imposed by slow thermal recovery, but retraining or online adaptation may be needed if pulse shapes drift [2411.15360].

The interpretation of multiplexed PNR has also been debated. One analysis argues that, for realistic detector parameters and commonly used coincidence windows, accidental coincidences dominate, so low-SNR multiplexed PNR outputs can be described by a classical wave model based on amplitude-threshold detection and coincidence aggregation. That work proposes that reliable PNR consistent with corpuscular light requires \(\Delta t\lesssim \Delta \tau\), \(\mathrm{SNR}_{\mathrm{dB}}\ge 3~\mathrm{dB}\), and operation away from saturation [2110.04245]. Other works, by contrast, focus on calibrated response matrices, POVMs, confusion matrices, and task-specific fidelity rather than on using PNR as stand-alone evidence for the discrete nature of light [1812.05422] [2411.15360]. A plausible implication is that “PNR” now denotes a family of calibrated measurement schemes whose validity depends on the detector model, the regime of operation, and the figure of merit being reported.

Current development directions are consistent across platforms: further jitter reduction, optimized amplification and thresholding, more than two threshold levels, advanced statistical classifiers such as Bayesian or EM inference, on-chip timing references, integrated electronics, and combinations of intrinsic PNR with modest multiplexing to extend dynamic range [2310.12472] [2505.09271]. In SNSPDs, Gaussian pulse shaping, EMG-based histogram modeling, larger or engineered kinetic inductance, and richer use of full-trace information all improve separability, but often at the cost of count rate [2510.00714] [2410.23162]. In multiplexed systems, larger \(M\), higher efficiency, and lower dark counts directly tighten the \(\mathcal{O}(1/M)\) approximation error and expand the domain in which click statistics track photon-number moments faithfully [2507.05554]. The field therefore continues to evolve along two parallel lines: better intrinsic detectors and better inference on imperfect detectors.

Source: https://www.emergentmind.com/topics/photon-number-resolution-pnr