
level. Geiger-mode single-photon avalanche diodes (SPADs) have found widespread use over the
past two decades [5–8].
All single-photon detectors exhibit dark counts in the absence of light: in avalanche photodiodes,
dark counts occur when avalanches are triggered by electrical carriers that are thermally generated
or emitted by trapping levels in the semiconductor [9, 10]. The rate at which these avalanches
are triggered, called the dark count rate, is heavily dependent on detector settings such as the
quantum efficiency and dead time [6, 7, 11]. Dark counts are a source of noise in any application
and increase with the bias voltage needed to enhance the quantum efficiency of the detector.
There have been several studies of dark current in avalanche photodiodes [5,12] and numerical
models have been developed to compare with experimental measurements [13].
However, a complete statistical characterization of the dark counts in 1550 nm SPADs has
not been reported. The purpose of our research is to explore the dark counts of SPADs by
using time-tagged photodetector measurements that can be used to generate histograms of
their interarrival time distributions, probability distributions, and entropy rates. A simplifying
assumption often made about the statistics of dark counts is that they follow a Poisson distribution.
We examine this assumption and show that it is not generally valid without proper processing of
the data. This processing is essential when one wants to accurately characterize the statistics of
single photon emitters or novel light sources.
It is long understood that Poisson statistics are not always an appropriate model for counts
from photomultiplier tubes [4, 14], image sensors such as CCD or CMOS sensors [15], general
purpose digital pulse processing systems [16], and SPAD photocounts [13,17]. However, it is still
often assumed that their dark counts follow Poisson statistics [9] whose probability distribution
is given by eqn. 1 [18, 19]
𝑃(𝑛)=(𝑟𝑇)𝑛𝑒(𝑟𝑇 )
𝑛!,(1)
where
𝑃(𝑛)
is the probability of
𝑛
counts being detected in the time interval
𝑇
, assuming
𝑟
is the
average rate of detected counts. We will define
𝜆=𝑟𝑇
as the average number of counts in the
interval
𝑇
. The detector quantum efficiency,
𝜂
, relates the optical power to the observed photon
count rate, 𝑟, as:
𝑟=𝜂(𝑃/ℎ𝜈),(2)
where 𝑃is the average optical power incident on the detector, and ℎ𝜈 is the energy per photon.
In many experiments the dark count rates are relatively small in comparison with the photon
count rates, therefore their effect on the distribution is negligible. However, there are other
practical cases where the dark count rate is comparable to or only slightly less than the photon
count rate. In these instances, we aim to quantitatively explore the dark count distributions and
their influence on a weak attenuated coherent light source which is known to be Poisson [1, 2].
We study the deviations of dark current statistics from Poisson statistics due to dead time and
afterpulsing [4, 20] and propose a method for correcting the histograms. First, we recommend
extending the detector dead time to eliminate the effects of afterpulsing as observed from the
interarrival distributions. Second, we suggest calculating the adjusted occurrence rate of dark
counts from the measured rate by using a standard formula. This enables us to determine the
correct slope for the interarrival time distribution of a Poisson process. Lastly, we will show that
a further step is necessary to correct the probability distributions for dead time which involves
implementing an iterative algorithm [21, 22].
2. Instrumentation
Single-photon avalanche diodes are named as such because when the reverse bias of its p-n
junction is raised above the breakdown voltage, just a single carrier can trigger an electrical
avalanche process, leading to a measurable current [9, 10]. To detect a subsequent photon, the