Optical parametric down-conversion (PDC) is a central tool in quantum optics experiments. The number of collected down-converted modes greatly affects the quality of the produced photon state. We use Silicon Photomultiplier (SiPM) number-resolving detectors in order to observe the photon-number distribution of a PDC source, and show its dependence on the number of collected modes. Additionally, we show how the stimulated emission of photons and the partition of photons into several modes determine the overall photon number. We present a novel analytical model for the optical crosstalk effect in SiPM detectors, and use it to analyze the results.
© 2012 Optical Society of America
For more than two decades, the optical process of parametric down-conversion (PDC)  has been the primary photon source in quantum optics experiments . This nonlinear process converts single photons into pairs of signal and idler photons, which possess quantum correlations between their various degrees of freedom. The statistical nature of this process is a subject of great interest, as the photon-number statistics of the process provide a direct indication of its non-classility [3, 4, 5, 6]. Additionally, the photon statistics have practical implications on state preparation schemes, where the quality of the produced state depends on the number of photons generated simultaneously within the PDC process .
Recent developments in photon-number resolving detectors [8, 9, 10, 11, 12] enable the implementation of a wide range of applications in quantum optics, including state preparation schemes [13,14,15], quantum state filtering  and quantum gates [17,18]. Number discriminating capabilities have also been shown to aid in probing eavesdropping attacks on quantum cryptography protocols . Furthermore, one of the most prominent advantages of number-resolving detectors is the ability to directly measure the photon-number statistics of an optical state.
The photon-number distribution of a PDC source can range between a thermal distribution, when a single down-converted mode is collected, and a Poisson distribution when an infinitely large number of spatial and spectral modes (multimode) is collected . The photon statistics of PDC sources has been directly measured in several works [21, 3, 4, 22], all of which report a Poisson distribution. However, the majority of quantum optics schemes take great care to ascertain that only a single mode is collected. Using an indirect homodyne detection method, a thermal distribution from a parametric amplifier was demonstrated . Transition between thermal to Poisson distributions was indirectly observed through the g(2) parameter  and through high-order PDC events . In this work, we use a photon-number resolving detector in order to directly measure the photon statistics of the signal and idler photons collected from a controlled number of modes, allowing us to observe the range of distributions between the single-mode and the multimode extremes.
The photon-number resolving detector used in this work is the Silicon Photomultiplier (SiPM) . The SiPM detector is composed of multiple avalanche-photo-diodes (APD) operating in Geiger-mode, combined on a single substrate. When a photon impinges on an APD element, an electric discharge can be generated. The photons impinging on the detector are distributed across its many elements. The output signals generated by all APDs are combined to form a single output pulse whose intensity is proportional to the number of triggered elements. In comparison to other technologies, the SiPM detector offers a better photon-number resolution (up to 20 photons, compared to 10 or less), easy integration in optical setups and operation at room temperature. On the other hand, its photon detection efficiency is 2–6 times lower, depending on the working wavelength.
The photon statistics measured by SiPM detectors may deviate from that of the original state. The three main reasons for these deviations are: 1) the non-perfect overall photon detection efficiency (loss), which is determined by the detection and coupling efficiencies, 2) false detections caused by thermal excitations (dark counts), and 3) optical crosstalk, in which secondary photons, created by carrier relaxation in one APD element, are detected by a neighboring element, falsely increasing the photon count . Another deviation is caused by events when two or more photons impinge on the same element and are detected as one. However, we have shown using a numerical model that this effect becomes significant only when the number of detected photons is higher than 20% of the total number of elements . The descriptions of loss and dark counts are quite straightforward. However, several currently available models for the crosstalk effect are either computationally difficult  or only applicable under limiting conditions [29, 30, 31]. We introduce here an analytical crosstalk model which can be applied for the entire range of experimental settings.
This paper is organized as follows: In Sec. 2 we present a model for the distortion effects in the SiPM detector, that includes a new approach to the crosstalk effect. This model is used for the interpretation of experimental data and the reconstruction of the original photon-number statistics. Sec. 3 describes the experimental setup for the generation of PDC states and their detection using an SiPM photon-number resolving detector. In Sec. 4 we present measurements of photon-number distributions of parametric down-conversion and their dependence on the number of collected modes. Conclusions are presented in Sec. 5.
2. The SiPM detection model
In order to interpret the photon-number statistics measured using SiPM detectors, the effects of loss, dark counts, and optical crosstalk must be properly modeled. If the measured photon-number distribution is represented by a vector of probabilities p⃗m, the original probability distribution p⃗o, is related to the measured distribution through the relationEq. 1, the original distribution can be reconstructed, in principle, by multiplying the measured distribution with the inverse of the distortion matrix. This method, however, is limited to relatively large probabilities and may produce non-physical results otherwise . Another approach, which we use here, is to apply a fitting algorithm which finds the original distribution based on a goodness of fit test.
The matrix M can be written as a product of three matrices, M = Mct · Mdk · Mloss, which represent the loss, dark counts (dk) and crosstalk (ct) effects. The individual matrices are constructed as follows:
The probability for loss is described by the relation 
The probability of dark counts is assumed to follow a Poisson distribution  of the form
Optical crosstalk occurs when a secondary photon generated during an electronic discharge in one APD element is detected in a neighboring APD element, thereby creating a secondary avalanche. Photons produced in the secondary avalanches can then continue to trigger further avalanches in additional elements (crosstalk-generated-crosstalk). The crosstalk probabilities in our model are calculated by dividing the crosstalk events into stages and recursively counting the number of possible crosstalk events in each stage. The recursion formula for the crosstalk matrix element is:
This novel analytical crosstalk model does not limit the number of crosstalk events and also accounts for crosstalk-generated-crosstalk. The model was developed since other available analytical models [29, 30, 31] were inconsistent with our experimental results due to limitations on the crosstalk probability values and the number of crosstalk events. Our model does not account for the geometrical arrangement of the APD elements in the detector, nor its finite size. However, we show below that for values within the experimental range, this analytical model agrees with a numerical model that has no such limitations . Additionally, this model is less computationally demanding than the numerical model. Thus, it is a useful tool in the interpretation of experimental data from photon-number resolving detectors, where crosstalk between the detection elements plays a crucial role.
3. The experimental setup
We have measured the photon-number distributions of the polarization modes of a stimulated type-II collinear PDC process. The experimental setup is shown in Fig. 1(a). Pulses of 150 fs at a repetition rate of 250 kHz are amplified, frequency doubled to 390 nm and used to pump a 2 mm long type-II collinear β – BaB2O4 (BBO) nonlinear crystal. The degenerate wavelength signal and idler photons at 780 nm, which are created with orthogonal polarizations (horizontal (H) and vertical (V)), are split at a polarizing beam splitter (PBS) and detected separately using two photon-number resolving detectors. This work focuses on measuring the photon-number distribution of a single polarization mode. Thus, the detection of both polarization modes is only used to evaluate the heralded detection efficiency of the SiPM detector. The results presented for a single polarization mode also apply to the photon pair statistics, as a collinear process produces the same statistics for pairs and for individual photons. Before the photons are coupled to the detectors, they are spectrally and spatially filtered by bandpass filters and by optical fibers, respectively. The number of spatial modes which are collected can be changed by using optical fibers with different core diameters and numerical apertures , and the number of collected spectral modes can be changed by using filters of different bandwidths.
The detection configuration is shown in Fig. 1(b). The optical signal from the photon source is coupled to an SiPM detector with 100 elements (Hamamatsu Photonics, S10362-11-100U) using optical fibers. The output electrical signal from the detector is amplified using low-noise amplifiers, digitized and then analyzed by FPGA electronics, where the number of detected photons is extracted in real time. This data is continuously transmitted to a computer, which presents the statistics. In order to minimize the effects of dark counts and after-pulsing , the 1 ns long sampling time of the analog electrical signal is synchronized with the arrival time of the photons. The dark counts are further reduced by moderately cooling the detectors using a thermoelectric cooler to ≈ −10°C, and the bias voltage level applied to the SiPM is modified accordingly.
4. Results and discussion
Figure 2 shows an example of a histogram of pulse intensities produced by a coherent state. The probability of detecting n photons is proportional to the area of the nth Gaussian. We obtain high photon-number discrimination for as high as 20 photons. The photon-number resolution error, which results from the overlapping between neighboring Gaussians, is smaller than 1% below 12 photons and approaches 12% for 20 photons. The number of resolvable photons is limited due to the decrease of the signal to noise ratio as the number of photons is increased.
We demonstrate the use of our model for reconstructing the original photon statistics with the photon-number distribution measurements of a thermal state. The thermal state was generated by collecting a single polarization mode of PDC using the setup shown in Fig. 1. The original thermal distribution can be written asEq. 5 due to optical crosstalk and dark counts. We measured the same thermal state with three different sets of values for the amount of loss, dark counts and crosstalk, obtained by operating the SiPM detector with different bias voltage values . As a result, we obtained three different measured distributions, although the original photon-number distribution was the same. The measurements are shown in Fig. 3. The main distortions in the distributions are the low values of n̄ due to losses and the deviation from a straight line in the semi-log plot when the number of measured photons exceeds 1. The latter is the effect of the crosstalk process, which falsely increases the probability of measuring high numbers of photons. The crosstalk probabilities of these measurements range from the highest probability value (maximum bias voltage), which creates the highest deviation, to the minimal probability (minimal bias voltage) which creates the lowest deviation.
The original distribution is reconstructed by fitting the measured data p⃗m to the function p⃗m = MctMdkMlossp⃗th, where p⃗th is the thermal distribution of Eq. 5 and n̄ is the only free parameter. The overall detection efficiency η, the average number of dark counts λdk, and the crosstalk probability ɛ, that are used for the construction of the distortion matrices, are measured separately. The value of η is defined as the heralded efficiency between the two down-converted polarization modes and is determined by the ratio between two-photon coincidence and single-photon counts in the limit n̄ ≪ 1. The dark counts and crosstalk parameters are obtained through a separate measurement of the Poissonian dark count statistics which are fit to the function p⃗m = MctMdk1⃗, where 1⃗ is a vector of zeros with 1 at the first position, and λdk and ɛ are the free fit parameters.
All three measurements in Fig. 3 were reconstructed to the same original state, within the margins of error, despite the different measured distributions. This result indicates the consistency of our evaluation of the dark counts, loss and crosstalk parameters. We compare these results to those obtained using the numerical crosstalk model introduced in Ref. . The numerical model is based on simulations of the detection process, and takes into account the geometrical configuration of the detector, the finite number of detection elements, and possible attempts to trigger neighboring cells which have already been triggered. The model is characterized through a parameter ɛnn, defined as the probability of generating crosstalk in one particular neighbor among the four nearest neighbors, rather than the overall probability of generating crosstalk among the nearest neighbors, which we define here as ɛ. The values of ɛnn obtained from the fit correspond to the values of ɛ through the relation ɛ = 1 −(1 − ɛnn)4 . The agreement between the two models shows that although the analytical model does not account for the geometry of the detector, it provides a good description of the crosstalk process for all crosstalk values in the experimental range.
We now consider the photon-number distributions of down-converted photons as a function of the number of collected modes. If n photons are distributed among s modes, such that n1 + n2 + ... ns = n, then the overall probability of measuring n photons is given by :Eq. 5 into Eq. 6, we obtain the following probability distribution
The dependence of the photon-number distribution on the number of modes is demonstrated by measuring the photon statistics collected from a controlled number of modes. Figure 4 shows four representative measurements of the photon-number distributions which range between a single mode (thermal) distribution, and a multimode (Poisson) distribution. The solid lines are fits to the function p⃗m = MctMdkMlossp⃗s, where p⃗s is defined in Eq. 7 with n̄ and s as the free parameters. We detected up to 14 photons and obtained good fits to the data. Notice that the difference between the distributions becomes more marked as the number of photons increases (see Figs. 4(b)–(e)). The difference is almost undetectable for low photon numbers. In fact, measurements of at least 10 photons are required for our experimental parameter range, in order to properly discriminate between distributions with different mode numbers.
The number of modes obtained for all combinations of fibers and filters is shown in table 1. The non-integer mode numbers should be interpreted as a weighted number of occupied modes and result from a non-uniform distribution of photons between the modes. We can discriminate between distributions containing up to about 10 modes, as distributions containing over 10 modes cannot be distinguished from Poisson distributions. Table 1 shows the proportional increase of the number of spectral modes with the increase in the bandwidth of the bandpass filters. The width of the phase-matched spectrum of the down-converted photons is estimated to be ≈10 nm. Due to spectral distinguishability, this width includes about 3 temporal modes [34, 25, 24]. The number of spatial modes measured with our system is also consistent with the choice of fibers. When the PDC photons, which passed through a 3 nm bandpass filter, were coupled into a single mode fiber for 780 nm, a thermal distribution was produced. Whereas, when the photons were coupled into a multimode fiber, a Poisson distribution was produced, even though the same filter was used. The relatively small number of modes obtained with the single-mode fiber for 1550 nm can be attributed to non-uniform mode coupling, which does not equally stimulate all of the modes in the fiber.
The stimulated nature of PDC can also be used to obtain the number of collected modes . The average number of down-converted photons in a single mode is determined by the stimulation parameter τi, which has a linear dependence on the pump field, the crystal length and the nonlinear coefficient of the crystal . Assuming an equal distribution of the photons between s modes, the mean number of photons n̄ is given by
In order to show the dependence of n̄ on the number of modes, we measured photon-number distributions as a function of the pump intensity for the different numbers of collected modes. The values of n̄ were extracted by fitting the different distributions to Eq. 7. The fit values of n̄ as a function of the total pump power are shown in Fig. 5. All measurements show a nonlinear dependence, indicative of a stimulated process. The number of modes s obtained from these fits corresponds well to the number of modes obtained through fits to the photon-number distributions (see Fig. 4). Furthermore, all fits resulted in similar values of the α parameter, in agreement with the model of Eq. 8.
In conclusion, we have built a detection setup, which incorporates a Silicon Photomultiplier (SiPM) as a photon-number resolving detector. Our setup enables good photon-number discrimination for as many as 20 photons. The output signal of the SiPM is analyzed in real time using FPGA electronics. The measured photon-number distribution probabilities are interpreted by modeling the effects of loss, dark counts and optical crosstalk, which distort the original photon-number probabilities. We present a useful analytical crosstalk model that is applicable for all crosstalk probability values in the experimental range and includes crosstalk-generated-crosstalk.
We have performed measurements of photon-number distributions of a single polarization mode of a collinear type-II PDC process. The high photon-number resolution of our system enables the differentiation between distributions collected from a different number of spectral and spatial modes. We present distributions which range from single-mode (thermal) to multimode (Poissonian) statistics. The results show how the photon-number distribution is determined by the number of collected modes. Furthermore, we show the dependence of the average number of photons on the number of collected modes, and the stimulated nature of parametric down-conversion.
References and links
1. L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge University Press, 1995).
2. P. G. Kwiat, K. Mattle, H. Weinfurter, A. Zeilinger, A. V. Sergienko, and Y. Shih, “New high-intensity source of polarization-entangled photon pairs,” Phys. Rev. Lett. 75, 4337–4341 (1995). [CrossRef] [PubMed]
3. E. Waks, E. Diamanti, B. C. Sanders, S. D. Bartlett, and Y. Yamamoto, “Direct observation of nonclassical photon statistics in parametric down-conversion,” Phys. Rev. Lett. 92, 113602 (2004). [CrossRef] [PubMed]
4. E. Waks, B. C. Sanders, E. Diamanti, and Y. Yamamoto, “Highly nonclassical photon statistics in parametric down-conversion,” Phys. Rev. A 73, 033814 (2006). [CrossRef]
5. O. A. Ivanova, T. S. Iskhakov, A. N. Penin, and M. V. Chekhova, “Multiphoton correlations in parametric down-conversion and their measurement in the pulsed regime,” Quantum Electron. 36, 951 (2006). [CrossRef]
8. D. Lincoln, “A large statistics study of the performance and yields of generation-6 vlpcs (histe-vi),” Nucl. Instrum. Methods Phys. Res. A 453, 177–181 (2000). [CrossRef]
9. D. Achilles, C. Silberhorn, C. Sliwa, K. Banaszek, I. A. Walmsley, M. J. Fitch, B. C. Jacobs, T. B. Pittman, and J. D. Franson, “Photon-number-resolving detection using time-multiplexing,” J. Mod. Opt. 51, 1499–1515 (2004).
10. D. Rosenberg, A. E. Lita, A. J. Miller, and S. W. Nam, “Noise-free high-efficiency photon-number-resolving detectors,” Phys. Rev. A 71, 061803 (2005). [CrossRef]
11. B. E. Kardynal, Z. Yuan, and A. J. Shields, “An avalanche photodiode-based photon-number-resolving detector,” Nat. Photonics 2, 425 (2008). [CrossRef]
12. E. Dauler, A. Kerman, B. Robinson, J. Yang, B. Voronov, G. Goltsman, S. A. Hamiltom, and K. Berggren, “Photon-number-resolution with sub-30-ps timing using multi-element superconducting nanowire single photon detectors,” J. Mod. Opt. 56, 365 (2009). [CrossRef]
13. P. Kok, H. Lee, and J. P. Dowling, “Creation of large-photon-number path entanglement conditioned on photodetection,” Phys. Rev. A 65, 052104 (2002). [CrossRef]
14. Y. Gao, P. M. Anisimov, C. F. Wildfeuer, J. Luine, H. Lee, and J. P. Dowling, “Super-resolution at the shot-noise limit with coherent states and photon-number-resolving detectors,” J. Opt. Soc. Am. B 27, A170–A174 (2010). [CrossRef]
20. L. Mandel, “Fluctuations of photon beams: The distribution of the photo-electrons,” Proc. Phys. Soc. 74, 233 (1959). [CrossRef]
22. M. Avenhaus, H. B. Coldenstrodt-Ronge, K. Laiho, W. Mauerer, I. A. Walmsley, and C. Silberhorn, “Photon number statistics of multimode parametric down-conversion,” Phys. Rev. Lett. 101, 053601 (2008). [CrossRef] [PubMed]
23. M. Vasilyev, S.-K. Choi, P. Kumar, and G. M. D’Ariano, “Investigation of the photon statistics of parametric fluorescence in a traveling-wave parametric amplifier by means of self-homodyne tomography,” Opt. Lett. 23, 1393–1395 (1998). [CrossRef]
25. W. Wasilewski, C. Radzewicz, R. Frankowski, and K. Banaszek, “Statistics of multiphoton events in spontaneous parametric down-conversion,” Phys. Rev. A 78, 033831 (2008). [CrossRef]
26. G. Bondarenko, B. Dolgoshein, V. Golovin, A. Ilyin, R. Klanner, and E. Popova, “Limited geiger-mode silicon photodiode with very high gain,” Nucl. Phys. B. (Proc. Suppl) 61, 347 (1998). [CrossRef]
27. P. Buzhan, B. Dolgoshein, L. Filatov, A. Ilyin, V. Kaplin, A. Karakash, S. Klemin, R. Mirzoyan, A. Otte, E. Popova, V. Sosnovtsev, and M. Teshima, “Large area silicon photonmultipliers: Performance and appli,” Nucl. Instrum. Methods Phys. Res. A 567, 78 (2006). [CrossRef]
28. L. Dovrat, M. Bakstein, D. Istrati, and H. Eisenberg, “Simulations of the detection process in silicon photomultiplier detectors,” arXiv:1109.0698v1 (2011).
29. I. Afek, A. Natan, O. Ambar, and Y. Silberberg, “Quantum state measurements using multipixel photon detectors,” Phys. Rev. A 79, 043830 (2009). [CrossRef]
30. M. Akiba, K. Tsujino, K. Sato, and M. Sasaki, “Multipixel silicon avalanche photodiode with ultralow dark count rate at liquid nitrogen temperature,” Opt. Express 17, 16885–16897 (2009). [CrossRef] [PubMed]
32. H. Lee, U. Yurtsever, P. Kok, G. M. Hockney, C. Adami, S. L. Braunstein, and J. P. Dowling, “Towards photostatistics from photon-number discriminating detectors,” J. Mod. Opt. 51, 1517–1528 (2004).
33. A. Yariv, Optical Electronics in Modern Communications (Oxford University Press, 1991).
34. P. J. Mosley, J. S. Lundeen, B. J. Smith, P. Wasylczyk, A. B. U’Ren, C. Silberhorn, and I. A. Walmsley, “Heralded generation of ultrafast single photons in pure quantum states,” Phys. Rev. Lett. 100, 133601 (2008). [CrossRef] [PubMed]
35. P. Kok and S. L. Braunstein, “Postselected versus nonpostselected quantum teleportation using parametric down-conversion,” Phys. Rev. A 61, 042304 (2000). [CrossRef]