Resolving Photon Numbers Using Ultra-High-Resolution Timing of a Single Low-Jitter Superconducting Nanowire Detector
Abstract
Photon-number-resolving (PNR) detectors are a key enabling technology in photonic quantum information processing. Here, we demonstrate the PNR capacity of conventional superconducting nanowire single-photon detectors by performing ultra-high-resolution time-tagging of the detector-generated electrical pulses. This method provides a viable approach for PNR with high detection efficiency and a high operational repetition rate. We present the implementation of such a PNR detector in the telecom C-band and its characterization by measuring the photon-number statistics of coherent light with tunable intensity. Additionally, we demonstrate the capabilities of the detection method by measuring photon-number correlations of non-classical states.
I Introduction
Efficient detectors play a crucial role in photonic quantum technology, and many recent milestones in the field can be traced back to the advancement of detector technology over the past decades [1, 2, 3]. Two types of detectors can be distinguished: threshold detectors, which provide a ”click” in presence of one or more photons, and detectors with photon number resolution (PNR).
PNR is an essential functionality for protocols like quantum teleportation and entanglement swapping [4], photonic quantum computing [5], and quantum-enhanced sensing [6]. Transition-edge sensors (TES) can be used, which intrinsically provide energy- and photon-number resolution encoded in the pulse form of the detector response [7, 8, 9, 10]. This comes at the cost of long recovery times. Presently, a commonly employed alternative approach is to achieve PNR functionality by multiplexing the photon flux onto several threshold detectors (e.g. superconducting nanowire single-photon detectors (SNSPDs)), where spatial-, temporal-, and multi-pixel approaches have been demonstrated [11, 12, 13, 14, 15]. This approaches are already established, but continuous further development in terms of efficiency and repetition rate are a key challenge.
Recent studies [16, 17], however, have provided findings that challenged the conventional perception of SNSPDs as threshold detectors. These studies demonstrated that the waveform of SNSPD detector pulses can provide information about the number of photons absorbed within the thermal relaxation time of the nanowire [18]. In a comprehensive detector tomography of an SNSPD, Endo et al. [17] showed that the photon number can be inferred by matching waveform traces recorded in a high-resolution oscilloscope. Moreover, the turn-on dynamics of SNSPDs, as modeled by Nicolich et al. [19] indicate that the rising slope of that waveform increases with a higher photon number of a detection event. Correspondingly, an earlier trigger at a particular threshold value should also allow for a distinction between photon numbers. This suggests that PNR detection should be possible using time-tagging devices with exceptional timing resolution.
In this work, we empirically study the timing statistics of SNSPD detector waveforms for photon wavepackets with different photon number distributions using ultra-fast time-tagging devices. Our experimental findings reveal that for our detection system and amplification circuit, simple threshold detection of the rising edge offers basic photon number resolution (PNR) functionality up to approximately N=2. Moreover, we find that substantial improvement of the PNR functionality can be achieved by not only considering the rising edge but also classifying detection events by the timing of the falling edge of the SNSPD-generated electrical pulse. This is also in good agreement with recent findings from an extensive principle component analysis of SNSPD electrical pulse shapes [20].
We experimentally implement this PNR detection scheme in the telecom C-band and characterize it by measuring the photon-number statistics of a ps-pulsed coherent light source with a tunable average photon number. We demonstrate that high-resolution time-tagging of the rising and falling edge of pulses enables PNR measurements with a single conventional low-jitter SNSPD, resolving up to N=5 photons in a single-shot measurement. To further demonstrate the applicability to quantum applications, we employ the approach for measurements of the joint photon-number distributions (JPND) of N00N states (N=2) generated from type-II spontaneous parametric down-conversion (SPDC) in a Hong-Ou-Mandel interferometer.
The manuscript is organized as follows. In Section II.1, we explain the detection approach and then evaluate the performance of the detection scheme using coherent states of light with different mean photon numbers. Further, we demonstrate the effectiveness of the method by measuring joint photon-number distributions of non-classical states in Section II.2. The main findings are summarized in Section III, together with an outlook on future steps towards broader impact.

II Results
II.1 PNR Characterization using coherent states
Experimental Setup
The experimental setup for PNR detection is shown in Fig. 2. We produce coherent light pulses with low mean photon numbers by attenuating laser pulses at nm with a duration of ps. To avoid saturation of the SNSPD we reduce the repetition rate of the laser system from MHz to kHz via pulse picking. For count rates up to approx. Mcps, the detection efficiency of the SNSPD (Single Quantum, EOS Series) is specified as .
For each measurement run, the detection efficiency is optimized by alignment of the polarization of the light with a fiber polarization controller placed before the SNSPD.

BS: beam splitter, PD: photodiode, VOA: variable optical attenuator, PC: polarization controller, SNSPD: superconducting nanowire single-photon detector, TTX: ultra-high-resolution time-tagging unit.
The arrival time of the wave packet on the SNSPD is precisely determined by splitting off a part of the laser with a beam splitter and detecting it with a high-bandwidth photodiode. The electric pulses from the photodiode and the SNSPD are then converted to time stamps using an ultra-high-resolution time-tagging unit (TTX, Swabian Instruments, TimeTagger X).
Results and Discussion
Fig. 1 A shows the two-dimensional histogram of the difference between the arrival time of the wave packet on the detector and the timing of the rising and falling edge of the electrical pulses of the SNSPD. The histogram exhibits distinct modes that correspond to the photon number in the coherent state, with a detected mean photon number of .
In order to determine the photon number of individual detection events, a simple analysis of the time difference between the trigger and the rising edge (which corresponds to projection along vertical lines in Fig. 1 A) allows to distinguish between N=1 and N>1 with high accuracy. Essentially, this can be understood as a consequence of the turn-on dynamics and the steeper rising edge of the detection waveform [19]. It is, however, evident that a clearer distinction of the photon-number modes can be found by considering not only the rising edge, but also the falling edge to optimally cluster detection events. A simple ad-hoc approach is to separate the clusters with parallel lines at an angle, i.e., by comparing different projections of the 2D histogram, along the falling edge timing and the found distinction lines. The optimal position and angle for the separation can be found by numerical optimization, minimizing the cross-talk between photon-number modes. A comparison in Fig. 1 B shows the significantly improved distinction between photon-number modes enabled by consideration of the falling edge timing. Zero-photon detection events are inferred from photodiode detections events without corresponding detection events on the SNSPD.
Note that this clustering of detector modes is only possible with low timing jitter of the measurement. The main contributors to the measured jitter is the SNSPD with a low jitter of ps root mean squared (RMS), enabled by cryogenic amplification, and the TTX with a jitter of ps RMS per channel.
To confirm the correspondence between clusters and photon numbers, we characterize the photon number statistics of coherent states of light. We vary the mean photon number of the coherent light by means of attenuation and compare the measured photon-number statistics with the expected Poissonian photon-number distributions.

The measured photon-number statistics shown in Fig. 3 are in good agreement with the fitted Poissonian photon-number distributions (solid line, truncated at N=5). Values for were obtained by fitting the truncated photon-number distribution (summing up occurrences of photon numbers N4). Graphical inspection shows the clear distinguishability of photon-number modes up to N=4.
To better quantify the benefit of optimal clustering, we model the marginal timing distribution (along the optimal projection axis) as a sum of the individual cluster distributions. We observe excellent agreement with our experimental data when modeling the individual clusters using a Voigt profile fit, as shown in Fig. 4.

By calculating the overlap between the fitted Voigt profiles and the identified photon-number bins, we can determine the cross talk between the photon-number channels. Fig. 5 displays the cross talk calculated for photon-number measurements conducted using only the rising edge timing and optimal clustering methods. Comparing these cross talk matrices reveals a significant reduction in cross talk when utilizing the optimal clustering method. This reduction enables a more precise identification of photon numbers in measurements.

II.2 Joint PNR Detection using non-classical light

Experimental Setup
To generate and characterize the joint photon-number distribution of N00N states [21], the experimental setup was modified as shown in Fig. 6. The N00N states with N=2 are generated by interfering photon-pairs from a type-II SPDC process in a Hong-Ou-Mandel interferometer. The photon pairs are generated in a periodically poled potassium titanyl phosphate (ppKTP) crystal, pumped by the second harmonic of the pump laser at nm. The two-step non-linear process to generate the photon pairs is pumped with a pulsed laser with a repetition rate of MHz and a transform-limited pulse duration of approx. ps. In this experiment, no pulse picking was required due to the low average count rates produced by the SPDC process. Every detection event resulting in a photon number measurement is associated with a time tag. Thus, we can identify correlated detection events across several detectors, enabling the characterization of joint photon-number distributions.
The transformation of the polarization state in the common-path Hong-Ou-Mandel interferometer allows switching between creating N00N states, and deterministically splitting photon pairs created in the SPDC process. The generation of N00N states is achieved by minimizing coincidence counts between the two SNSPD channels. Respectively, a maximization of those coincidence counts results in the deterministic splitting of photon pairs generated in the SPDC process.
The second detector channel employed for these measurements has an estimated detection probability of and a timing jitter of ps RMS.
Results and Discussion
The obtained JPNDs in Fig. 7 show the expected photon-number correlations between the two PNR detectors depending on the polarization transformation in the Hong-Ou-Mandel interferometer. The difference between the two generated states becomes especially evident when one examines the occurrences of two-photon events depicted in Fig. 8. This confirms the accurate identification of correlations of individual PNR measurements.

Furthermore, the observed JPNDs provide further information on the system’s performance parameters. For example, by comparing the rates of detecting single photons with the rates of detecting photon pairs, it is possible to estimate the total system efficiency of approx. (including fiber coupling efficiency, system transmission, and detection efficiency). The extension to higher-order photon number contributions could provide further characterization metrics, such as multi-pair generation rates.

III Conclusions
In this work, we have explored the implementation of a PNR detector using a single conventional low-jitter superconducting nanowire detector and threshold analysis of the electric waveform with ultra-high-resolution timing.
The method requires precise knowledge of the arrival time of the wave packet on the detector, either by splitting off a part of the wave packet or by analysis of a pump pulse involved in a non-linear optical process needed for the generation of a certain state. Nevertheless, this approach can dramatically reduce resource overhead compared to approaches based on multiplexing and promises substantially shorter recovery rates than PNR detectors based on TES.
We have shown a successful implementation and characterization of the detection method with clear distinguishability of photon modes up to photons. This is confirmed by showing good agreement of measured photon-number distributions of a coherent light source with Poissonian photon-number statistics. It also indicates a potential extension to a simple PNR detection tomography method using coherent states of light. Coherent states offer several advantages for the characterization of the PNR detection system. They are easily and reliably generated and offer non-vanishing probabilities for higher photon numbers, even at low intensities. Additionally, they are resilient to the inevitable losses in the light propagation and detection and thus offer a good reference state. However, this property makes it impossible to assess the loss characteristics of the detection scheme. The method thus needs to be employed with accurately characterized attenuators. Alternatively a well-characterized threshold detector could provide a necessary reference. A detailed discussion of absolute detector tomography, however, is beyond the scope of this initial work.
Furthermore, we demonstrate the usefulness of the shot-by-shot PNR detection approach by characterizing the joint photon-number distributions of N00N states with N=2. This would be a daunting task for methods that rely on the analysis of the complete waveform, in particular when scaling to higher photon numbers.
In conclusion, we have demonstrated that it is possible to resolve high photon numbers with minimal overhead and infrastructure already available to researchers in many labs. We note that the approaches outlined herein are heavily dependent on the specifics of the detection system and hope that our results will inspire further studies into the prediction of the turn-on and turn-off behaviour of superconducting nanowires. Although a more comprehensive understanding of the observed phenomena needs to be explored in further work, we are confident that our empirical results already provide a simple means for leveraging PNR capability of existing detectors that are widely available in the research community.
Moreover, we believe that further enhancement of the approach could be accomplished either by suitable pre-conditioning of electronic detector signals or by evaluating more than two threshold levels.
Acknowledgments
We acknowledge support from the European Union’s Horizon 2020 Research and Innovation Action under Grant Agreement No. 899824 (SURQUID). GS acknowledges funding of this project from the Carl-Zeiss-Stiftung. GS is a member of the Max Planck School of Photonics, supported by the German Federal Ministry of Education and Research, the Max Planck Society, and the Fraunhofer Society. The authors gratefully acknowledge the support of René Sondenheimer, Oskar Kohout, Cristina Amaya, Shreya Gouravaram, and Helmut Fedder for their invaluable contributions to discussions and extend their thanks to Markus Leipe and Cristina Amaya, for their assistance in operating the SPDC source. We also acknowledge Single Quantum B.V. for providing a preliminary experimental setup and their detector systems for first tests.
References
- Christ et al. [2013] A. Christ, A. Fedrizzi, H. Hübel, T. Jennewein, and C. Silberhorn, in Single-Photon Generation and Detection, Experimental Methods in the Physical Sciences, Vol. 45, edited by A. Migdall, S. V. Polyakov, J. Fan, and J. C. Bienfang (Academic Press, 2013) pp. 351–410.
- Marsili et al. [2013] F. Marsili, V. B. Verma, J. A. Stern, S. Harrington, A. E. Lita, T. Gerrits, I. Vayshenker, B. Baek, M. D. Shaw, R. P. Mirin, and S. W. Nam, Nature Photonics 7, 210 (2013).
- Chang et al. [2021] J. Chang, J. W. N. Los, J. O. Tenorio-Pearl, N. Noordzij, R. Gourgues, A. Guardiani, J. R. Zichi, S. F. Pereira, H. P. Urbach, V. Zwiller, S. N. Dorenbos, and I. Esmaeil Zadeh, APL Photonics 6, 036114 (2021).
- Żukowski et al. [1993] M. Żukowski, A. Zeilinger, M. A. Horne, and A. K. Ekert, Physical Review Letters 71, 4287 (1993).
- Kok et al. [2007] P. Kok, W. J. Munro, K. Nemoto, T. C. Ralph, J. P. Dowling, and G. J. Milburn, Reviews of Modern Physics 79, 135 (2007).
- Edamatsu et al. [2002] K. Edamatsu, R. Shimizu, and T. Itoh, Physical Review Letters 89, 213601 (2002).
- Rosenberg et al. [2005] D. Rosenberg, A. E. Lita, A. J. Miller, and S. W. Nam, Phys. Rev. A 71, 061803 (2005).
- Lita et al. [2008] A. E. Lita, A. J. Miller, and S. W. Nam, Optics Express 16, 3032 (2008).
- Lita et al. [2010] A. E. Lita, B. Calkins, L. A. Pellouchoud, A. J. Miller, and S. Nam, in Advanced Photon Counting Techniques IV, Vol. 7681, edited by M. A. Itzler and J. C. Campbell, International Society for Optics and Photonics (SPIE, 2010) p. 76810D.
- Schmidt et al. [2018] M. Schmidt, M. von Helversen, M. López, F. Gericke, E. Schlottmann, T. Heindel, S. Kück, S. Reitzenstein, and J. Beyer, Journal of Low Temperature Physics 193, 1243 (2018).
- Jönsson and Björk [2020] M. Jönsson and G. Björk, Physical Review A 101, 013815 (2020).
- Achilles et al. [2004] D. Achilles, C. Silberhorn, C. Sliwa, K. Banaszek, and I. Walmsley, Optics letters 28, 2387 (2004).
- Divochiy et al. [2008] A. Divochiy, F. Marsili, D. Bitauld, A. Gaggero, R. Leoni, F. Mattioli, A. Korneev, V. Seleznev, N. Kaurova, O. Minaeva, G. Gol’tsman, K. G. Lagoudakis, M. Benkhaoul, F. Lévy, and A. Fiore, Nature Photonics 2, 302 (2008).
- Fitch et al. [2003] M. J. Fitch, B. C. Jacobs, T. B. Pittman, and J. D. Franson, Physical Review A 68, 043814 (2003).
- Stasi et al. [2022] L. Stasi, G. Gras, R. Berrazouane, M. Perrenoud, H. Zbinden, and F. Bussières, High-efficiency and fast photon-number resolving parallel superconducting nanowire single-photon detector (2022), arXiv:2207.14538 [quant-ph] .
- Cahall et al. [2017] C. Cahall, K. L. Nicolich, N. T. Islam, G. P. Lafyatis, A. J. Miller, D. J. Gauthier, and J. Kim, Optica 4, 1534 (2017).
- Endo et al. [2021] M. Endo, T. Sonoyama, M. Matsuyama, F. Okamoto, S. Miki, M. Yabuno, F. China, H. Terai, and A. Furusawa, Opt. Express 29, 11728 (2021).
- Natarajan et al. [2012] C. M. Natarajan, M. G. Tanner, and R. H. Hadfield, Superconductor Science and Technology 25, 063001 (2012).
- Nicolich et al. [2019] K. L. Nicolich, C. Cahall, N. T. Islam, G. P. Lafyatis, J. Kim, A. J. Miller, and D. J. Gauthier, Physical Review Applied 12, 034020 (2019).
- Bartley [2023] T. Bartley, private communication (2023).
- Dowling [2008] J. P. Dowling, Contemporary physics 49, 125 (2008).