Comparative study of afterpulsing behavior and models in single photon counting avalanche photo diode detectors

TL;DR

This study compares afterpulsing behaviors in avalanche photodiodes from different manufacturers, revealing significant differences and the presence of high-order afterpulses, which impact quantum optics applications and require individual detector calibration.

Contribution

It provides a comparative analysis of afterpulsing in APDs from multiple manufacturers, highlighting the need for detector-specific calibration and introducing a procedure for it.

Findings

Significant differences in afterpulsing behavior between detectors.

Detection of high-order afterpulses not covered by standard models.

Calibration procedure for individual detectors demonstrated.

Abstract

Single-photon detectors, like Avalanche Photo Diodes (APDs), have a great importance in many fields like quantum key distribution, laser ranging, florescence microscopy, etc. Afterpulsing is a typical non ideal behavior of APDs, operated in Geiger mode, which adversely affects any application which measures the number or timing of detection events. Several conflicting studies have tried to model afterpulsing behavior and link it to fundamental semiconductor physics. Here we describe experiments we performed on at least 2 different detectors from 3 manufacturers, all commonly used in quantum optics setups. We report on the timing structure of the counts. By fitting different models to these data, we found statistically significant differences between the detectors used. Furthermore, we report on the presence of high order afterpulses which are not accounted for in any of the standard models of afterpulse behavior. This might have significant implication e.g. in quantum key distribution systems, as the modeling of the devices in use will have to be part of the security analysis for some systems. Significant differences between manufacturers and even between seemingly identical detectors of the same make and batch suggest that every individual detector needs to be calibrated. We describe a reliable procedure to perform such a calibration and present our findings on three different commercial products.

Tables

Table 1. Table 1 : Table comparing various makes of commercially available single photon detectors we tested. The total after pulse probility shows a piece to piece variability as expected. More importantly, the piece to piece variation extends to the best fitting theoretical model. This should not happen for theories based on fundamental properties of semiconductors. Clearly the applicability of all 3 mathematical models discussed here is very limited and ill defined.

Make of detector	Serial number	Factory date	P(AP) (%)	Best model
Excelitas SPCM-NIR	29860	27 Aug. 15	1.29043 $\pm$0.00071	Multiple exponential
	29864	27 Aug. 15	0.01130 $\pm$0.00025
	32403	26 Oct. 16	0.87774 $\pm$0.00028	Power
	32404	26 Oct. 16	2.01987 $\pm$0.00039
	32405	26 Oct. 16	1.38750 $\pm$0.00033
	32406	26 Oct. 16	1.56231 $\pm$0.00035
	32424	02 Nov. 16	2.91185 $\pm$0.00048	Hyperbolic sinc
	32425	02 Nov. 16	3.96425 $\pm$0.00057
	32426	02 Nov. 16	4.30335 $\pm$0.00059
	32427	02 Nov. 16	4.69090 $\pm$0.00060
	32428	02 Nov. 16	3.31100 $\pm$0.00051
	32429	02 Nov. 16	2.98385 $\pm$0.00049
PerkinElmer SPCM-AQ4C	195 (Ch 1)	25 Jan. 07	0.28542 $\pm$0.00079	Multiple exponential
	195 (Ch 2)	25 Jan. 07	0.28542 $\pm$0.00079
Laser Components $\tau$-SPAD-fast	1019917	01 Nov. 13	5.10601 $\pm$0.00310	Multiple exponential
	1019920	01 Nov. 13	8.52410 $\pm$0.00213

Equations

$$P_{exp}(t)=\sum _k{A}_k{e}^{-\frac{t}{{\tau }_k}}+d,$$

$$P_{pow}(t)=A\cdot t^{-\lambda }+d,$$

$$P_{sinc}(t)=2\cdot A\cdot \frac{\text{sinh}(\Delta \cdot t)}{t}\cdot e^{-\gamma \cdot t}+d,$$

Full text

Comparative study of afterpulsing behavior and models in single photon counting avalanche photo diode detectors

Abdul Waris Ziarkash,1, † Siddarth Koduru Joshi,1,,† Mario Stipčević,2 and Rupert Ursin3,*

*1**Institute for Quantum Optics and Quantum Information, Austrian Academy of Sciences Vienna A-1090

*2**Ruder Boskovic Institute, Center of Excellence for Advanced Materials and Sensors and Division of Experimental Physics, Zagreb 10000, Croatia

*3**Institute for Quantum Optics and Quantum Information, Austrian Academy of Sciences, Vienna A-1090

Vienna Center for Quantum Science and Technology (VCQ).*

† These authors contributed equally and are joint first authors.

***Siddarth.Koduru.Joshi@oeaw.ac.at.

****Rupert.Ursin@oeaw.ac.at.

OCIS codes: (270.0270) Quantum optics; (270.5568) Quantum cryptography; (060.5565) Quantum communications; (250.1345) Avalanche photodiodes; (040.6070) Solid state detectors; (120.4800) Optical standards and testing; (120.4640) Optical instruments; (120.6200) Spectrometers and spectroscopic; instrumentation; (180.2520) Fluorescence microscopy; (170.0170) Medical optics and biotechnology .

References and links

• [1]

M.D. Eisaman, J. Fan, A. Migdall, and S.V. Polyakov.

Invited review article: Single-photon sources and detectors.

Review of scientific instruments, 82(7):071101, 2011.

• [2]

Robert H. Hadfield.

Single-photon detectors for optical quantum information applications.

Nature photonics, 3(12):696–705, 2009.

• [3]

Winfried Denk, James H. Strickler, Watt W. Webb, et al.

Two-photon laser scanning fluorescence microscopy.

Science, 248(4951):73–76, 1990.

• [4]

David E. Kuhl, Jorge R. Barrio, Sung-Cheng Huang, Carl Selin, Robert F Ackermann, James L. Lear, J.L. Wu, TH Lin, and Michael E. Phelps.

Quantifying local cerebral blood flow by N-isopropyl-p-[123I]iodoamphetamine (IMP) tomography.

Journal of nuclear medicine: official publication, Society of Nuclear Medicine, 23(3):196–203, 1982.

• [5]

Anja Wagner, Heiko Mahrholdt, Thomas A. Holly, Michael D. Elliott, Matthias Regenfus, Michele Parker, Francis J. Klocke, Robert O. Bonow, Raymond J. Kim, and Robert M. Judd.

Contrast-enhanced MRI and routine single photon emission computed tomography (SPECT) perfusion imaging for detection of subendocardial myocardial infarcts: An imaging study.

The Lancet (London, England), 361(9355):374–379, 2003.

• [6]

Haifeng Wang, Terry B. Huff, Daniel A. Zweifel, Wei He, Philip S. Low, Alexander Wei, and Ji-Xin Cheng.

In vitro and in vivo two-photon luminescence imaging of single gold nanorods.

Proceedings of the National Academy of Sciences of the United States of America, 102(44):15752–15756, 2005.

• [7]

Jamie Holder, R.W. Atkins, H.M. Badran, G. Blaylock, S.M. Bradbury, J.H. Buckley, K.L. Byrum, D.A. Carter-Lewis, O. Celik, Y.C.K. Chow, et al.

The first veritas telescope.

Astroparticle Physics, 25(6):391–401, 2006.

• [8]

R. G. W. Brown, K. D. Ridley, and J. G. Rarity.

Characterization of silicon avalanche photodiodes for photon correlation measurements 1: Passive quenching.

Applied Optics, 25(22):4122, 1986.

• [9]

Thiago Ferreira Da Silva, Guilherme B. Xavier, and Jean Pierre Von Der Weid.

Real-time characterization of gated-mode single-photon detectors.

IEEE Journal of Quantum Electronics, 47(9):1251–1256, 2011.

• [10]

Sergio Cova, A. Lacaita, and G. Ripamonti.

Trapping Phenomena in Avalanche Photodiodes on Nanosecond Scale.

IEEE Electron Device Letters, 12(12):685–687, 1991.

• [11]

M. Stipčević.

A novel active quenching circuit for single photon detection with geiger mode avalanche photodiodes.

Applied Optics, 48:1705–1714, 2009.

• [12]

S. V. Polyakov and A. L. Migdall.

High accuracy verification of a correlated-photon-based method for determining photoncounting detection efficiency.

Optics Express, 15:1390–1407, 2007.

• [13]

Mario Stipčević and Daniel J. Gauthier.

Precise Monte Carlo Simulation of Single-Photon Detectors.

5:1–5, 2014.

• [14]

Michael Ware, Alan Migdall, Joshua C. Bienfang, and Sergey V. Polyakov.

Calibrating photon-counting detectors to high accuracy: background and deadtime issues.

Journal of Modern Optics, 54(2-3):361–372, 2007.

• [15]

Nitin Jain, Birgit Stiller, Imran Khan, Vadim Makarov, Christoph Marquardt, and Gerd Leuchs.

Risk analysis of Trojan - horse attacks on practical quantum key distribution systems.

IEEE Journal on Selected Topic in Quatum Electronics, 21(3):1077–260X, 2014.

• [16]

H. T. Yen, S. D. Lin, and C. M. Tsai.

A simple method to characterize the afterpulsing effect in single photon avalanche photodiode.

Journal of Applied Physics, 104(5):1–5, 2008.

• [17]

Boris Korzh, Tommaso Lunghi, Kateryna Kuzmenko, Gianluca Boso, and Hugo Zbinden.

Afterpulsing studies of low-noise InGaAs/InP single-photon negative-feedback avalanche diodes.

Journal of Modern Optics, (July):1–7, 2015.

• [18]

A. C. Giudice, M. Ghioni, S. Cova, and F. Zappa.

A process and deep level evaluation tool: Afterpulsing in avalanche junctions.

In European Solid-State Device Research Conference, pages 347–350, 2003.

• [19]

Christoph Schaeff Sven Ramelow Mario Stipčević Gerhard Humer, Momtchil Peev and Rupert Ursin.

A simple and robust method for estimating afterpulsing in single photon detectors.

Journal of Lightwave Technology, 33(14):3098–3107, 2015.

• [20]

Mark A. Itzler, Xudong Jiang, and Mark Entwistle.

Power law temporal dependence of InGaAs/InP SPAD afterpulsing.

Journal of Modern Optics, 59(17):1–9, 2012.

• [21]

Mark A. Itzler, Mark Entwistle, and Xudong Jiang.

High-rate photon counting with Geiger-mode APDs.

IEEE Photonic Society 24th Annual Meeting, PHO 2011, (April 2016):348–349, 2011.

• [22]

Alessandro Restelli, Joshua C. Bienfang, and Alan L. Migdall.

Time-domain measurements of afterpulsing in InGaAs/InP SPAD gated with sub-nanosecond pulses.

Journal of Modern Optics, 59(August):1–7, 2012.

• [23]

We could not illuminate the SPCM-AQ4C fiber coupled detector in the same manner as the other free space coupled detectors, so we have left out these results for consistency.

• [24]

Fang-Xiang Wang, Wei Chen, Ya-Ping Li, De-Yong He, Chao Wang, Yun-Guang Han, Shuang Wang, Zhen-Qiang Yin, and Zheng-Fu Han.

Non-Markov property of afterpulsing effect in single-photon avalanche detector.

pages 18–20, 2016.

• [25]

M. D. Eisaman, J. Fan, A. Migdall and S. V. Polyakov.

Invited Review Article: Single-photon sources and detectors.

Acta Medica Okayama, 67(4):259–263, 2013.

• [26]

D. B. Horoshko, V. N. Chizhevsky, and S. Ya. Kilin.

Afterpulsing model based on the quasi-continuous distribution of deep levels in single-photon avalanche diodes.

Journal of Modern Optics, (September).1–5, 2016.

• [27]

Mario Stipčević, Daqing Wang, and Rupert Ursin.

Characterization of a Commercially Available Large Area, High Detection Efficiency Single-Photon Avalanche Diode.

Journal of Lightwave Technology, 31(23):3591–3596, 2013.

• [28]

K. E. Jensen, P. I. Hopman, E. K. Duerr, E. A. Dauler, J. P. Donnelly, S. H. Groves, L. J. Mahoney, K. a. McIntosh, K. M. Molvar, A. Napoleone, D. C. Oakley, S. Verghese, C. J. Vineis, and R. D. Younger.

Afterpulsing in Geiger-mode avalanche photodiodes for 1.06 $\mu$m wavelength.

Applied Physics Letters, 88(13):10–13, 2006.

• [29]

Daniel Gottesman, H-K Lo, Norbert Lutkenhaus, and John Preskill.

Security of quantum key distribution with imperfect devices.

In Information Theory, 2004. ISIT 2004. Proceedings. International Symposium on, page 136. IEEE, 2004.

• [30]

Daniel Gottesman, Hoi-Kwong Lo, Norbert Lütkenhaus, and John Preskill.

Security of quantum key distribution with imperfect devices.

arXiv preprint quant-ph/0212066, 2002.

• [31]

M. Stipčević, B.G. Christensen, P.G.Kwiat, PG and D.J. Gauthier.

Advanced active quenching circuits for single-photon avalanche photodiodes.

International Society for Optics and Photonics. Proceedings. page 98580R–98580R, 2016.

• [32]

Daniel J. Gauthier, Christoph F. Wildfeuer, Hannah Guilbert, Mario Stipcevic, Bradley G. Christensen, Daniel Kumor, Paul Kwiat, Kevin T. McCusker, Thomas Brougham, and Stephen Barnett.

Quantum key distribution using hyperentangled time-bin states.

Invited lecture. Conference on Coherence and Quantum Optics. Optical Society of America, 2013.

• [33]

We chose a value of 900 ns such that we include a significant amount of the afterpulse tail while avoiding any potential effects due to the earliest possible subsequent pulse 1 $\mu$s later. In silicon based SPADs the afterpulse tail dies off some hundred nanoseconds, whereas in InGaAs detectors it can last for several $\mu$s.

1 Introduction

Applications of single-photon detectors, which focus on the timing of a very weak optical signal, mostly use single-photon avalanche diodes (SPADs) operated above their breakdown voltage in Geiger mode. Such as in photonics, quantum processing tasks, laser ranging, fluorescence microscopy, neural imaging with blood flow tomography, contrast-enhanced MRI, two-photon luminescence imaging, astronomical telescopes, etc. [1, 2, 3, 4, 5, 6, 7]. They are widely manufactured and sold by many different companies and as such exhibit different properties/behavior. To efficiently and accurately perform these experiments we must account for all non-ideal behavior of the detectors used [8, 9]. This is especially true for sensitive applications like quantum communication, because the security of any real world implementation (i.e. with a high transmission loss) depends on the devices used. Hence the precise modeling/characterization of the non-ideal behavior of single photon detectors (and all other components of the quantum communication device) are critical for practical security proofs.

An ideal single-photon detector generates one and only one electric pulse for every incident photon. However, in a real detector, it is possible that a single incident photon results in more than one electrical pulse per incident photon. This is known as afterpulsing. [8] In this work we, approach this problem from the perspective of an end user. Consequently, we use an ipso-facto definition of an “afterpulse” as any pulse in addition to and following an isolated detection event and its subsequent dead time, regardless of its etiology. This behavior has been extensively studied due to its importance in semiconductor physics in general. It has been suggested that afterpulsing can be linked to charges trapped in the deep levels of the semiconductor’s band structure and released at a later time [10] as well as to an other causes (such as those described in Refs. [11, 12, 13, 14]

Afterpulsing has different implications, depending on the application of the detector used. It can result in an overestimation of the total count rate by up to 10 % as well as a reduction of the duty cycle and detection efficiency due to the increased dead times. In fluorescence microscopy it could lead to an overestimation of the concentration of fluorophores. In quantum communication the overestimation of coincidence events leads to a larger Quantum Bit Error Rate (QBER). Afterpulsing can also adversely affect security [15]. Afterpulsing also poses a significant problem to the measurement of photon arrival times both in quantum communication protocols as well as for laser ranging. Hence a proper characterization of the afterpulse behaviour could improve the accuracy of several measurements in metrology in general.

Several previous works have attempted to characterize the afterpulse behavior and fit the results to various models [18, 10, 20, 19, 26, 28, 27]. Worryingly, these papers do not agree on the most suitable model describing the statistics of the arrival times. Most of these studies have focused on one make/manufacture of detector at a time. In this work we take a more comprehensive approach and perform a comparative study on 3 different makes of detectors (for each make we compare up to 12 individual detectors with almost consecutive serial numbers). We test the universality of the various theoretical models and attempt to simply resolve the conflicts raised by previous works. To do so we eschew the common practice of studying the timing auto-correlation of the detector signals in favor of the more comprehensive timing cross-correlation between the emission of a photon and all the following detector signals. We describe this procedure in Section 3. We then fit to the various standard models described in Section 2 and compare these fits in Section  4.1. Due to our use of cross-correlations we were also able to detect higher orders of afterpulses which we present in Section 4.2 and finally in Section 4.3 we discuss the corrections that can be applied.

2 Canonical models of afterpulsing

Much of the previous work on afterpulsing in SPADs focused on modeling their behavior. In this section we present the three most common models. We go on to show the inadequacy of all of these models in Section 4.1. The characteristic decay of the afterpulse probability was sometimes thought to depend on the deep level in which the charge is trapped. Initial models [18, 10, 16, 19, 17] considered the decay from several distinct deep levels and proposed the “multiple exponential model” where the afterpulsing probability ($P_{exp}(t)$) at a time t is given by:

[TABLE]

where $\tau_{k}$ and $A_{k}$ is the de-trapping lifetime and amplitude factor for the $k^{\mathrm{th}}$ deep level, and d is the offset due to noise counts.

Furthermore [20, 21] considered a continuum of deep levels, in InGaAs/InP detectors, that could trap a charge and found that the power law model was a good empirical fit. Here the afterpulse probability ($P_{pow}(t)$) is given by:

[TABLE]

where $\lambda$ is a effective decay constant and $A$ is the initial afterpulse probability.

In an another attempt to create a more physically meaningful model Ref. [26] derived the “hyperbolic Sinc model” from the Arrhenius law (once more assuming a continuum of levels), where the afterpulse probability ($P_{sinc}(t)$) is given by:

[TABLE]

where $\Delta$ and $\gamma$ are both functions of the minimum and maximum energies of the deep levels in which charges may be trapped.

In this paper, we compare the above three standard mathematical models (see Equations 1, 2 and 3) with the measured behavior of several different single photon counting detectors (SPCM-AQ4C from PerkinElmer (see Figure 2a), SPCM-NIR from Excelitas (see Figure 2c) and $\tau$-SPAD-fast from Laser Components (see Figure 2e). We show that none of these canonical models are universal and vary between individual detectors as well as makes of detectors. We also report on the presence of several higher order afterpulses. An investigation of the region of reduced detection efficiency called the “dead-time” found between the signal pulse and the afterpulse shows yet another characteristic variation between makes of detectors, which is not part of the models used so far.

3 Experimental Setup

The general scheme of the setup used for the experiments is shown in Figure 1. A function-generator triggers a 798 nm laser with a repetition rate of 0.25 to 1 MHz. The emitted light pulse of $\approx$ 1 ns length gets attenuated in a neutral attenuator to $<$0.1  photons per pulse on average and arrives at the active area of the respective detector in use. The arrival time at the detector and the trigger from the function generator are individually sent to a time-tagging-module (TTM). We then calculated the temporal cross-correlation histogram ($g^{(2)}$) between those two signals using a bin width of 1 ns. We normalize to the number of detected laser pulses to obtain the pulse probability per bin. To ensure that our results are not skewed by an individual “faulty” detector we tested at least 2 detectors of each type (same manufacturer’s model number with nearly consecutive serial numbers) all of similar ages.

We consider the first event arriving at the TTM, after a given laser pulse to be the ”detection pulse” which stems from the real photon impinging on the active area of the detector. Any detection event after that, but before the next laser pulse, should be due to either afterpulsing or noise. From an end user perspective it is impossible to distinguish a true afterpulse from another count that just happened to occur at a similar time delay. Consequently, we do not differentiate between an afterpulse due to processes in the semiconductor diode from those due to secondary avalanches caused by any other means (such as faults in the quenching circuit). This is typical of most applications and previous studies of these detectors ranging from (quantum) optics experiments to fluorescence microscopy.

We fit each of the standard models (see Equations 1, 2 and Equation 3) to the tail of the afterpulse as seen in the $g^{(2)}$ obtained for each different SPAD detectors. Figure 2 shows the resultant fits as well as the residuals.

All current models of afterpulsing discussed, concern the timing structure of the afterpulse tail. An ideal model of a detector should include its entire behaviour including the duration of the dead time, the jitter and the rising edge of the afterpulse peak. As in most previous previous studies, we then isolate the peak (due to afterpulses) and its tail in time. The tail of the afterpulse of all detectors are shown in Figure 2 after normalization and accidentals correction have been applied. The accidental rate is estimated using a region well after the end of the afterpulse tail and before the next incident photon (this method of correction has been empirically verified in Section 4.3). We define the total afterpulse probability to be the sum total of the afterpulse probability in each bin of the corrected $g^{(2)}$ histogram for a duration of 900 ns. [33]. Regardless of the model, in order to obtain a good fit we had to ignore the first two points of the afterpulse tail.

4 Results

4.1 Comparison of afterpulsing models

We observed a wide variation between different tested detectors. We tested at least 2 and up to 12 identical detectors (in terms of make, model number and year of purchase) from each of three manufacturers. We saw a drastic change in the total afterpulse probability between otherwise identical detectors(see Table 1).

On the whole, the power law model had the most consistent good fit and the hyperbolic Sinc model the worst. Unlike the hyperbolic Sinc model, the power law model approximates the experimental data more accurately independent of the detector used AQ4C (see Figure 2a), SPCM-NIR (see Figure 2c), $\tau$-SPAD-fast (see Figure 2e). Ref. [26] showed that the hyperbolic sinc model was better than the power law model using id100-MMF50 SPAD module from idQuantique, we were able to duplicate those results only with 6 of the 16 detectors we tested.

Figure 2 shows sample fits of the three models for one detector from each manufacturer. In this case the power law model is fitting the entire range of the tail better than the Hyperbolic sinc model. The fits can be quantitatively compared using the residuals seen in Figure 2 which are well within the $\pm$ 2 standard deviation range (shown in red) for the power law model but not for the hyperbolic sinc model.

Power law model as a simple analytic model function leads to some good results and proved to be the best model for 4 of the 16 detectors tested. The power law model was primarily shown to be effective for InGaAS/InP detectors [21, 20, 22] and we now show that it is a reasonably good model for some Si detectors. However, it’s estimated parameters gives a rather loose description of the physical origin of afterpulses. A truly good model should be able to predict the behavior of the entire tail region. The power law model only provides a suitable fit if the beginning two points of the afterpulse tail are ignored.

The multiple exponential model provides a good fit when $k\geqslant 5$ and is slightly better than the power law model for the SPCM-AQ4C detector. It proved to be the best fit for the behaviour of 6 of the 16 tested detectors. It is nearly commensurate with the power law model for the $\tau$-SPAD. For the SPCM-NIR, the exponential model gives a good fit only when we ignore the higher order afterpulse peaks (which are discussed in Section 4.2). We observed that it was possible to get different fits of almost the same quality with different characteristic times ($\tau_{k}$ also called the de-trapping times) depending on the number of exponentials used. Thus, upon examining the fitting procedure we concur with [20] who state that “It is evident from this whole fitting procedure that the extracted values for the de-trapping times depend entirely on number of exponentials in the model function and the range of hold-off times used in the data set.”

In case of the hyperbolic sinc model, we consistently obtain unsuitable fits to all but 6 detectors. In figure 2 we can see that the residuals for several fits are much larger than those for other models. However, for 6 detectors, the hyperbolic sinc model was found to be the most suitable.

The total probability of obtaining an afterpulse ($P(AP)$) for each detector is shown in Table 1. Clearly, afterpulse behavior and probabilities vary drastically between brands and between individual specimens of the same brand.Various detectors of the same make from the same company (even if manufactured 1 week appart) provide evidence in support of different and contradictory mathematical models. As do detectors from different brands. Clearly a universal model to describe afterpulsing behaviour does not exist and every detector needs to be characterized individually.

4.2 Higher order afterpulses

Any detection event may cause an afterpulse which, being a detection event, may induce secondary and further afterpulses called “higher order afterpulses”. Usually afterpulsing is a small effect and the afterpulsing probability distribution function is smeared such that the higher order afterpulses are improbable. However, in the presence of strong enough twilighting [12], photon detections and afterpulses accumulate in a narrow peak that appear just after the dead time. We note this behavior in SPCM-NIR and to a lesser extent in SPCM-AQ4C, as can be seen in Figure 3. The higher order afterpulses are clearly visible as a series of peaks after the main peak with a period exactly equal to the dead time of the detector. The presence of such higher order afterpulsing was first speculated in Ref. [20] but, to the best of our knowledge, never reported.

For the SPCM-NIR, the signal pulse and the first order of afterpulse are 22 ns apart, which exactly corresponds to the duration of the dead time of this particular detector module. The measured time intervals between all following higher orders of afterpulses (as seen in the inset of the Figure 3) have the same time delay of 22 ns. We have obtained a similar plot for three detectors with the same model number: in each case, peaks appear separated by the dead time of the particular detector. Such a behavior is clearly undesirable, notably in time-resolved spectroscopy where higher order peaks could be mistaken for, or mask the true signal.

We verified that the higher order afterpulses seen are not the result of stray light, the shape of the laser pulse, electronic noise, impedance mismatch or optical reflections. We also processed the data with different bin widths and laser pulse frequencies to ensure that the observed higher order afterpulse peaks are not due to digitization noise.

Further, we note that it is possible to explain the area of the $n^{\mathrm{th}}$ higher order peak based on the probability of the first afterpulse ($P(AP)$) as $P(AP)^{(n+1)}+$ the probability of an afterpulse in the bin just before the $n^{\mathrm{th}}$ higher order peak $\times$ the number of bins in the peak. This geometric progression agrees to within 4 to 6 % for the $2^{nd}$ to $5^{th}$ order afterpulses.

4.3 Background and accidentals corrections during the dead time

In this subsection we discuss the corrections that we can apply to the cross correlated histograms we used in the previous subsections. In typical quantum optics experiments, there is a probability that a coincidence is detected between two different detectors erroneously, we call these coincidences “accidentals”. Typically, they can be estimated from Poissonian statistics as: $r_{acc}=r_{1}r_{2}t_{c}$, where $r_{acc}$ is the rate of the accidental counts, $r_{1}$ and $r_{2}$ are the count rates of the individual detectors and $t_{c}$ is the coincidence time window used. By introducing a controlled amount of continuous background illumination we experimentally verify this standard practice. In the duration well after a detection event, this provides a very good estimate (with a maximum variation of $\ll$ 3 %) of the behavior (see Figures 4(a) and 4(b)).

An interesting consequence of using the cross-correlation technique described above is the non-zero probability of a count during the “dead-time”. We attribute these detection events to the probability that the detector did not click because of the laser pulse but did click within the dead-time region.

When used for some tasks (like quantum communication), the end user may not be able to correct for the accidentals obtained while computing the cross correlation histogram. In these cases it is important, for a complete model of the detector’s behavior, that these accidentals during the dead-time be correctly accounted for.

We use a continuous wave battery powered LED as a steady and controllable source of background illumination in addition to the attenuated laser pulses (see Figure 1). We then measured the $g^{(2)}$ histograms for different background count rates. Figure 4(a) shows these histograms for the SPCM-NIR while Figure 4(b) shows the same for the $\tau$-SPAD-fast modules [23]. We show that the increase in click probability during the dead-time scales linearly with the amount of additional background illumination (see Figures 4(a) and 4(b)).

In the regions before the detection event, well after the detection event and even during the afterpulse, the measured accidental corrections agreed with the value calculated assuming Poissonian statistics.

5 Conclusion

We have clearly demonstrated that different standard models are required to appropriately describe the distribution of electrical signals generated by different detectors. This explains the conflicting nature of several previous studies; for example Refs. [20, 21, 22] show strong evidence for the power law model while Refs. [18, 10] show equally compelling evidence for the multiple exponential model and Ref. [26] provides evidence for the hyperbolic sinc model. Recently [24] showed that the afterpulse probability is dependent on past events — a property not considered in the exponential, power, or hyperbolic sinc models. By comparing previously reported results to our own, we realize that there is a large variation between the different commonly used manufactures/makes of detectors and between individual detectors. The suppression of afterpulses by different quenching methods [25] leads us to believe that the afterpulsing behavior is more dependent on the electronic quenching circuit used rather than the properties of the semiconductor (such as the presence and distribution of discrete/continuous/quasi-continuous deep levels). Our data clearly proves that none of the current theoretical models are universal which makes it hard to draw conclusions about the underlying mechanism based on fundamental semiconductor physics.

As seen in Table 1 supposedly identical detectors manufactured within a week of each other, provide evidence to support one model over the other. Several previous studies have shown that one or the other mathematical model fits the afterpulsing behaviour of a few individual detectors. Further many of these studies use such evidence to confirm or disprove hypotheses/assumptions about semiconductor physics. Clearly, with such extreme individual variation no conclusion about the physics behind these mathematical models can be drawn from any similar test. Especially tests with a small number of individual detectors or with detectors of only one kind.

We also report on the presence of higher order afterpulses in one of the tested detector models. To make this possible, unlike in several previous measurements [26, 19, 27] who used an auto-correlation signal, we use a cross-correlation histogram between the detector and the trigger. This allows us to look for both higher order afterpulses as well as the behavior during the dead time. These higher order afterpulses can cause large errors in measurements of the arrival times of photons and must be carefully accounted for. We were able to exclude frequency dependence in our detectors since we repeated all the experiments described above for several different repetition rates (all far from detector saturation) of the laser pulses ranging from 10 kHz to 1.2 MHz and found no significant variation.

Further, the total afterpulsing probability for several individual detectors all of the same manufacturer and part number varied drastically (in some cases by a factor of 460!) implies that every individual detector needs to be calibrated for all applications that need to accurately measure count rates or arrival times of photons.

For many years, afterpulsing has been extensively studied from a semiconductor physics based perspective, where it is important to understanding how trapped charges/energy levels decay. However, to correctly study this, one must separate the effects due to the behaviour of the diode or seimiconductor junction from the effects due to the electronics and quenching circuits. For example, based on modeling the semiconductor junction, a longer dead time is thought to lead to a lower afterpulsing probability [19]. However, as seen in Table 1 and Figure 3 both the SPCM-AQ4C and the $\tau$-SPAD-fast have a significantly longer dead time than the SPCM-NIR but the former has a lower afterpulse probability while the latter has a larger one. In this case it is impossible to draw a conclusion about the relationship between dead time and afterpulse probability without carefully considering the quenching circuits used. This type of analysis requires proprietary and confidential information about the diode and circuits used by the manufacturer; which is inaccessible to a typical end-user. Instead of drawing potentially erroneous conclusions about semiconductor behavior based on which model had a better fit, we focus on an application oriented perspective. For most practical purposes, it is sufficient to understand the statistical nature of the afterpulses rather than their causal mechanisms.

We also would like to mention, that we did not consider a possible aging effect of the detectors. All measurements were performed in the time span of only a few months. Commercial products, may have to be characterized repeatedly during their long operational lifetimes. The study of any aging effects on the statistical behaviour of SPAD’s in general is an interesting avenue for further exploration.

Most applications of SPADs are hindered by afterpulsing, in many cases these ill effects can be corrected for if each individual detector is properly characterized. The individual characterization is necessary due to the large variation in both the total afterpulse probability and the mathematical form of the probability distribution between detectors of the same model number, age and manufacturer under nearly identical laboratory test conditions. It is possible, although inadvisable, to ignore all clicks for several hundred ns after any detection event. This effectively increases the dead time and avoids the bulk of afterpulses but this severely limits the maximum count rates and detection efficiency due to saturation-like effects. To correctly account for the ill effects of afterpulses it is sufficient to characterise each individual detector prior to/during use, using the method described in this work.

From a quantum communication perspective, our method of characterization of every individual detector can be included into the overall device dependent security analysis [29, 30]. These characterizations can also be used to improve the accuracy of results in quantum meterology.

Acknowledgements

The Authors thank FFG-ALR (contract Nr. 844360), ESA (contract Nr. 4000112591/14/NL/US) and MoSES 533-19-14-0002, the Austrian Academy of Sciences as well as Austrian Science Fund (FWF): SFB F40 (FOQUS) for their financial support. We would like to thank Dmitri Horoshko, Vyacheslav Chizhevsky and Sergei Kilin for discussing their methods and sharing the results of their paper [26].