Probing higher order correlations of the photon field with photon number resolving avalanche photodiodes

J. F. Dynes; Z. L. Yuan; A. W. Sharpe; O. Thomas; A. J. Shields

doi:10.1364/OE.19.013268

1. Introduction

Photon correlation functions, formalized by Glauber in the 1960s [1], have become fundamental tools in characterizing the quantum statistical properties of light. From an experimental point of view, it is the class of normally ordered correlation functions [2] corresponding to photon absorption as opposed to the anti-normally ordered class [3] that feature as most accessible. Up until now, the most significant normally ordered photon correlation function (other than first order) is the second order correlation function, first measured by Hanbury-Brown and Twiss (HBT) [4, 5]. For example, the second order correlation function makes possible photon states to be categorized as classical or quantum mechanical. Poissonian or bunched photon statistics would indicate classical behavior - or if sub-Poissonian statistics are observed then it would be concluded the photon field had a quantum character [6].

Despite the strong interest in measuring the second order correlation function, comparatively little research to date has focused on measuring higher order correlations. Partly this is due to the complexity of the task at hand; characterizing higher order correlations of a given light source can require as many detectors as the order of the correlation function being measured itself [7]. Nevertheless, from both a practical and fundamental point of view, characterization of higher order correlation functions is of considerable importance. In view of applications, for example quantum key distribution (QKD) [8], statistical knowledge of the source is vital for security. Furthermore, it has been shown for QKD using realistic single photon sources such as semiconductor quantum dots [9–11], the secure transmission distance can approach that of using a theoretical, true single photon source with knowledge of the source higher order photon correlation functions [12].

Traditionally, monitoring the photon field statistics by measuring higher order correlations is a time consuming process [7]. The time taken to acquire a large number of coincidences to achieve acceptable statistical bounds on the correlation function scales exponentially with the order of the correlation. This is usually impractical in many situations. Alternative techniques for measuring higher order correlations have been put forward such as two-photon detection [13], spatial multiplexing of superconducting nanowire detectors [14] and temporal multiplexing based on linear optics [15]. While these techniques have their own individual merits, they suffer from other significant drawbacks including electrical cross-talk, cryogenic operating temperatures and low coincidence rates. Moreover, to measure higher order correlations where the order is greater than about four, the above techniques may prove problematic due to scalability issues. We note that very recently a proposal to filter out cross-talk in multi-pixel detectors for accessing higher order correlations has been reported [16].

Measuring higher order correlations with just two detectors is a relatively unexplored area. One possible method is to make use of the detector photon number resolving capability. Most single photon detectors though can only differentiate between the zero photons and one or more photons [17]. A single or multi-photon event produces the same voltage pulse height, making photon number resolving impossible. On the other hand, photon number resolving (PNR) detectors produce voltage pulses in proportion to the number of photons. Hence, n-photon events can be precisely selected depending on the discrimination levels. In this paper, we present a proof of principle experiment for demonstrating sensitivity to higher order correlations based on high speed PNR detectors. We employ two gigahertz gated, PNR avalanche photodiodes (APDs) in a normal HBT setup [4, 5]. This setup can measure the usual normally ordered and normalized second order correlation function g(2) defined as $g (2) = \frac{< â^{† 2} â^{2} >}{< â^{†} â >^{2}}$ where â ^† and â are the usual photon creation and annihilation operators. The form of g(2) presented above is the most salient in terms of categorizing photon states. However, we demonstrate that due to the PNR capability of the detectors, two-photon, three photon and up to n-photon events can be selected using photon counting discriminators, where n is the PNR capability of the detector. Coupled with the joint photon detection of using both PNR detectors together, the setup is sensitive to higher order correlations. As we show below for a bunched photon source, the sensitivity to photon bunching is significantly improved over the standard g(2) measurement.

2. Theory

According to Glauber [1], the definition of the nth order, normally ordered, normalized correlation function is given by:

g (n) = \frac{< â^{† n} â^{n} >}{< â^{†} â >^{n}}

When n = 2, g(2) can be measured using the usual Hanbury-Brown Twiss experiment with a double single photon detector arrangement. Such a setup is depicted in Fig. 1(a) with two InGaAs APDs as the detectors. In this case, both detector discriminators are set so 0-photons are rejected while 1-photon and multi-photon events are sampled. This is readily achieved using most non-photon number resolving single photon detectors. As an example, if the source under investigation were purely thermal with a coherence time greater than the temporal resolution of the detector, a standard HBT experiment would yield g(2) = 2. Specifically, the measurement involves sampling the coincident counts at times smaller than the coherence time of the source and dividing by the accidental coincidences (obtained at times greater than the source coherence time).

Fig. 1 (a) Schematic of the experimental setup. The solid circles represent photons from the source under test. SD-APD1 & SD-APD2: self-differencing avalanche photodiodes (APDs), discr.: electrical signal discriminator.(b) Electrical setup for measuring weak avalanches. (c) The three types of photon source employed and pictorial representations of their frequency spectra: (i) Filtered multi-mode laser operating slightly above lasing threshold (FML) (the mode in red box denotes allowed mode), (ii) DFB laser operating near lasing threshold (LNT) & (iii) DFB laser operated well above threshold (LAT).

Download Full Size | PDF

If the detectors are photon number resolving, then information on higher order correlation functions can be obtained by adjusting the discrimination levels. Sampling only 2-photon events by rejecting 0, 1 and ≥ 3 photon events can be achieved with window discriminators with the window set to count only 2-photon events. We can define a quantity called the higher order coincidence, for example, γ(4):

γ (4) = \frac{< â^{† 4} â^{4} >}{< â^{† 2} â^{2} > < â^{† 2} â^{2} >}

γ(4) clearly contains information regarding higher order correlation functions as we can rewrite γ(4) in terms of the usual normalized correlation functions:

γ (4) = \frac{g (4)}{g (2) * g (2)}

In the example given above with the thermal source, the higher order coincidence rate would be higher than g(2) with γ(4) = 6.

The above analysis can be extended to higher photon numbers. We can generalize Eq. (2) as:

γ (n_{1} + n_{2}) = \frac{< â^{† (n_{1} + n_{2})} â^{(n_{1} + n_{2})} >}{< â^{† n_{2}} â^{n_{2}} > < â^{† n_{1}} â^{n_{1}} >}

Here n_i indicates that the window discrimination levels of detector i is placed around the n_i photon number states. Eq. (4) can also be re-written in terms of normalized correlation functions:

γ (n_{1} + n_{2}) = \frac{g (n_{1} + n_{2})}{g (n_{1}) * g (n_{2})}

If the photon number resolution of both detectors is n, then the higher order coincidence can be measured up to 2n. Note that this arrangement circumvents the non-unity detector efficiency of the detectors and any inherent optical losses in the system.

For simplicity and to illustrate the proof of principle experiment presented below, we use threshold discriminators rather than window discriminators. To incorporate that fact it is necessary to integrate over all photon number states greater than the lowest allowed photon number state. In this case Eq. (4) is modified to read:

γ (n_{1} + n_{2}) = \frac{Σ_{n_{1}, n_{2}}^{n_{max}} < â^{† (n_{1} + n_{2})} â^{(n_{1} + n_{2})} > η_{1}^{n_{1}} η_{2}^{n_{2}}}{Σ_{n_{1}, n_{2}}^{n_{max}} < â^{† n_{2}} â^{n_{2}} > < â^{† n_{1}} â^{n_{1}} > η_{1}^{n_{1}} η_{2}^{n_{2}}}

where n_max is the maximum photon number the detector can measure and η_i is the i system single photon detection efficiency. The photon sources appropriate for use with this technique have detected mean photon fluxes much less than 2n_max. Note that Eq. (6) approximates Eq. (4) when η_i is relatively small. This is due to the weighting of the higher order terms in the sums of Eq. (6) scaling as

η_{i}^{n_{i}} η_{j}^{n_{j}}

.

3. Photon number resolution with InGaAs avalanche photodiodes

Figure 1(a) shows the layout of the experimental setup. The photon source under test is coupled into single mode fiber and photons propagate along the fiber to a 50:50 fiber beamsplitter. Photons in each arm are detected by an InGaAs APD operating in a self-differencing mode at a clock rate of 1GHz [18]. The self-differencer permits cancelation of the strong APD capacitive response leaving behind the previously obscured weak avalanche. This is achieved by dividing the electrical APD output into two arms, delaying one arm by an integer number of clock cycles before recombining the two arms using a subtractive combiner, Fig. 1(b). After self-differencing, the avalanches from the detectors are amplified before being discriminated by ultra-fast (1ns) discriminators. The discriminated outputs are then sent to a two channel time-correlated photon counting card (correlator) for recording the photon arrival times. For measuring the detector avalanche voltage distribution, a high speed oscilloscope is employed directly on the amplified detector output. The APDs feature single photon detection efficiencies of around 15–20%.

Figure 2 shows the experimentally measured avalanche voltage distribution of a detector for two types of photon source. When excited with a 1550nm distributed feedback (DFB) type pulsed laser operating well above threshold (LAT, Fig. 1(c)), the avalanche voltage distribution shows a number of distinct peaks, black circles Fig. 2. These peaks correspond to different photon number, n in the incident light field. The avalanche voltage distribution can be satisfactorily modeled assuming the photon number distribution of the source is Poissonian with a mean detected photon flux, μ ∼ 2.6, black line in Fig. 2. Each photon peak, n ≥ 1 is assumed to be Gaussian and to reflect the statistical broadening due to avalanche noise. The widths of the photon peaks are scaled as √n relative to the 1-photon peak width. Such a technique has successfully been employed to model photon number distributions for both InGaAs and Silicon APDs [19, 20]. The areas of each photon peak are proportional to the expected Poisson photon number distribution probabilities [19]. The main features of the LAT trace (black circles) in Fig. 2 are reproduced by the black line.

Fig. 2 The experimentally measured distribution of avalanche voltages for a pulsed laser source operating well above threshold (black circles) and for a pulsed, mode filtered non-DFB laser source operating close to threshold (red circles). Also shown are the theoretical distributions of avalanches for a Poissonian source (black line) and a partially bunched source (red line). The dashed lines correspond to the cross-over avalanche voltages for the experimental and theoretical avalanche distributions, as described in the text. Inset: Depicts the theoretical avalanche distributions for the Poissonian source (black line) and the partially bunched source (red line) up to an avalanche voltage of 0.75V.

Download Full Size | PDF

The ability to detect differences in the photon statistics for dissimilar photon sources is now illustrated. The photon source was replaced with a multi-mode pulsed laser operating slightly above threshold with a 1nm bandpass filter that serves to pass a single spectral mode close to a wavelength of 1550nm (FML in Fig. 1(c)(i)). Here the average detected photon flux was similar to before with μ ∼ 2.8.

This photon source gives rise to a rather different avalanche voltage distribution, (red circles, Fig. 2). The 0-photon peak is higher than for the laser source operating well above threshold, indicating the presence of more vacuum states. Conversely there is a suppression of the 2-photon and 3-photon number states. The photon source is expected to exhibit intensity fluctuations due to operation near to lasing threshold. Such fluctuation is manifest experimentally as photon bunching [21]. Bunching can be viewed as a departure from the independent statistics of a Poissonian source. The degree of bunching is accentuated through selection of a single mode by filtering [22]. The mathematical form of the photon number distribution for this particular case is not known precisely, so we choose an analytic photon counting distribution function representing a linear superposition of a coherent state field with mean photon number μ_s and a chaotic field with mean photon number μ_n [23, 24]. The following relations fix μ_s & μ_n from experimentally measurable quantities, namely average detected photon flux, μ = μ_s + μ_n and g(2) = μ_n(μ_n + 2μ_s)/μ ² +1 [23,24]. Using μ = 2.8 and g(2) = 1.2 [25], the resulting modeled avalanche distribution is shown in Fig. 2 (red line).

For both types of photon source, the theoretical avalanche distribution overestimates large avalanche heights. The origin of this effect is unclear at present. However, it is quite likely due to the quenching by the APD series resistance causing large avalanches from photon number states with n > 3 to be reduced in size (the linear model assumes a linear avalanche voltage dependence on photon number). The quenching effect can be somewhat quantified by comparing the experimental curves crossover point with the theoretical crossover point. These occur at 0.18 and 0.26V respectively, as indicated by the dashed lines in Fig. 2. Hence there is a 30% reduction in avalanche height at an avalanche voltage of 0.26V compared to the simple linear model prediction.

4. Higher order correlations of the photon field

Having shown that the detectors are sensitive to differences in the photon number distribution, we now present experimental results that detectors can be used to measure higher order correlations of the photon field. The setup illustrated in Fig. 1(a) is used. We stress that the setup will be exactly the same as a conventional one that measures the second order correlation function g(2) if the photon number resolving function is disabled. By setting the discrimination level at 300mV, a value far greater than single photon generated avalanches, for both detectors SD-APD 1 and SD-APD 2, the measurement gives a correlation result that is different from g(2). Plotted in Fig. 3(a), is a histogram for source FML when both the discrimination levels of SD-APD 1 and SD-APD 2 are at high values. At zero time delay (Δt = 0) a prominent correlation peak is observed due to photon bunching. The correlation here is given by Eq. (4), γ. The correlation value is measured as the ratio of peak height at t=0 to the average height at Δt ≠ 0. γ is much higher than the the equivalent g(2) value of g(2) = 1.2, as shown by the histogram in the inset of Fig. 3(a).

Fig. 3 (a) Experimentally measured photon count histogram as a function of detector relative delay for source FML. SD-APD 1 and SD-APD 2 discrimination levels are both set at a high value. Inset: Second order correlation measurement for source FML with g(2) = 1.2. SD-APD 1 and SD-APD 2 discrimination levels are both set at the lowest possible discrimination level. (b) Photon correlations, γ of the three sources shown in Fig. 1(c): FML (squares), LNT (circles) & LAT (triangles). The arrows correspond to the mean positions of the photon number states assuming a avalanche voltage linear dependence on photon number. Inset: Detail of γ for source LAT. Error bars are ± 1 standard deviation of the plotted (mean) values and are plotted for all data; however only some error bars are visible due to the size of the data points.

Download Full Size | PDF

For different discrimination levels, correlation plots similar to Fig. 3(a) are collected and the correlations, γ evaluated. The detector discriminator level of detector SD-APD 1 was kept at 300mV. Figure 3(b) plots the resulting correlations, γ as a function of SD-APD 2 discrimination level for the three sources. All three sources show γ > 1 even at the lowest discrimination level used. This can be readily understood in terms of elevated photon bunching through higher order correlations. The n^th order correlation is the expectation of the joint detections of n photons correlations. For photons that are indistinguishable and statistically dependent, there is a factorial increase of the available permutations of photon amplitudes as n rises [26, 27]. At higher discrimination levels, higher photon number states feature more markedly than lower photon number states and the correlation is therefore expected to be higher than for lower photon number states. This is in direct contrast to what would be expected for the case of statistically independent photons. In this case the photon statistics are Poissonian and no bunching is possible for any order of correlation.

Photon source LAT shows the least amount of correlation; this is expected for a laser source operating well above threshold with largely Poissonian statistics [2]. The correlation for LAT is relatively flat over the entire discrimination range of SD-APD 2. γ for this source attains ∼ 1.3 at the very highest D2 discrimination level employed, Fig. 3(b), inset. Operating the same laser near to threshold (LNT) shows a marked increase in correlation which is due to increased intensity fluctuations leading to increased photon bunching. Source FML displays the highest correlation values; attaining γ ∼ 47.1 ±4.5 for these measurements. Intensity fluctuations for this source are by far the highest of all sources under test.

All three sources under study feature different photon number statistics which is incorporated in our model through the mean photon numbers μ_s & μ_n. These mean photon numbers are derived from the total average photon number μ and g(2). The second order correlation function g(2) was measured separately for each of the sources FML, LNT & LAT and the values g(2) = 1.2, 1.075 & 1.01 respectively were used in the calculations (the inset of Fig. 3(a) shows the measured g(2) photon coincidence histogram for source FML).

As remarked previously, avalanche height self-limiting due to the APD series resistance is likely to distort the linear avalanche voltage dependence on photon number for high photon number. We can incorporate this fact by increasing n_max in the model for the detectors. Figure 4 shows the simulated results based on Eq. (6) for γ plotted with up to n_max = 7 photon detection events. Each point corresponds to selecting exactly n avalanche detections based on a linear dependence of the avalanche peaks (as shown by the arrows in Fig. 3(b)). The x-axis is plotted reciprocally to emphasize the avalanche self-limiting. γ for all sources shows a similar growth and dependence to that observed experimentally, qualitatively confirming our underlying model.

Fig. 4 Simulated correlation value as a function of discriminated SD-APD 2 Avalanche Voltage (V) state for SD-APD1 discriminator level set at n = 7 photon number state. FML (squares), LNT (circles) & LAT (triangles).

Download Full Size | PDF

Finally, we remark that although the agreement between the experimental and theoretical data is good, it is only qualitative. The photon number resolving capability of the detector worsens at higher photon number due to avalanche broadening. Avalanche broadening causes adjacent photon number avalanches to overlap. This fact is not incorporated in the model for simplicity. The next step is to realize a SD-APD with narrow photon number avalanche peaks having minimal overlap; similar to the performance of visible light photon counters (VLPCs), for example [28]. This will permit efficient window discrimination between n and n ± 1 photon number avalanche peaks.

5. Conclusion

In summary we have demonstrated a unique, proof of principle technique for inferring higher order photon correlations based on high speed (GHz) photon number resolving APDs. Using three light sources with differing photon statistics we have shown that higher order correlations of these sources can be detected. Moreover, accessing these higher order correlations illustrates the sensitivity of the method when trying to distinguish between sources compared to a usual g(2) measurement. We believe the technique will be of use and importance in characterizing photon sources involved in the rapidly evolving field of quantum information science.

Acknowledgments

We thank M. B. Ward for enlightening discussions and a critical reading of the manuscript. This work was partially supported by EU FP7 ICT project “Q-Essence”.

References and links

1. R. J. Glauber, “The quantum theory of optical coherence,” Phys. Rev. 130(6), 2529–2539 (1963). [CrossRef]

2. D. F. Walls and G. J. Milburn, “Coherence properties of the electromagnetic field,” in Quantum Optics (Springer-Verlag, 2008), pp. 29–55. [CrossRef]

3. K. Usami, Y. Nambu, B. S. Shi, A. Tomita, and K. Nakamura, “Observation of antinormally ordered Hanbury Brown-Twiss correlations,” Phys. Rev. Lett. 92(11), 113601 (2004). [CrossRef] [PubMed]

4. R. Hanbury-Brown and R. Q. Twiss, “Correlation between photons in two coherent beams of light,” Nature 177(4497), 27–29 (1956). [CrossRef]

5. T. Jeltes, J. M. McNamara, W. Hogervorst, W. Vassen, V. Krachmalnicoff, M. Schellekens, A. Perrin, H. Chang, D. Boiron, A. Aspect, and C. I. Westbrook, “Comparison of the Hanbury Brown-Twiss effect for bosons and fermions,” Nature 445(7126), 402–405 (2007). [CrossRef] [PubMed]

6. H. J. Kimble, M. Dagenais, and L. Mandel, “Photon antibunching in resonance fluorescence,” Phys. Rev. Lett. 39(11), 691–695 (1977). [CrossRef]

7. G. S. Agarwal, “Field–correlation effects in multiphoton absorption processes,” Phys. Rev. A 1(5), 1445–1459 (1970). [CrossRef]

8. N. Gisin, G. Ribordy, W. Tittel, and H. Zbinden, “Quantum cryptography,” Rev. Mod. Phys. 74(1), 145–195 (2002). [CrossRef]

9. P. Michler, A. Kiraz, C. Becher, W. V. Schoenfeld, P. M. Petroff, L. Zhang, E. Hu, and A. Imamoglu, “A quantum dot single–photon turnstile device,” Science 290(5500), 2282–2285 (2000). [CrossRef] [PubMed]

10. C. Santori, M. Pelton, G. Solomon, Y. Dale, and Y. Yamamoto, “Triggered single photons from a quantum dot,” Phys. Rev. Lett. 86(8), 1502–1505 (2001). [CrossRef] [PubMed]

11. Z. Yuan, B. E. Kardynal, R. M. Stevenson, A. J. Shields, C. J. Lobo, K. Cooper, N. S. Beattie, D. A. Ritchie, and M. Pepper, “Electrically driven single–photon source,” Science 295(5552), 102–105 (2002). [CrossRef]

12. Y. Adachi, T. Yamamoto, M. Koashi, and N. Imoto, “Boosting up quantum key distribution by learning statistics of practical single-photon sources,” New J. Phys. 11(11), 113033 (2009). [CrossRef]

13. Y. Qu, S. Singh, and C. D. Cantrell, “Measurements of higher order photon bunching of light beams,” Phys. Rev. Lett. 76(8), 1236–1239 (1996). [CrossRef] [PubMed]

14. M. J. Stevens, B. Baek, E. A. Dauler, A. J. Kerman, R. J. Molnar, S. A. Hamilton, K. K. Berggren, R. P. Mirin, and S. W. Nam, “High–order temporal coherences of chaotic and laser light,” Opt. Express 18(2), 1430–1437 (2010), http://www.opticsexpress.org/abstract.cfm?URI=OPEX-18-2-1430. [CrossRef] [PubMed]

15. M. Avenhaus, K. Laiho, M. V. Chekhova, and C. Silberhorn, “Accessing higher order correlations in quantum optical states by time multiplexing,” Phys. Rev. Lett. 104(6), 063602 (2010). [CrossRef] [PubMed]

16. D. A. Kalashikov, S. H. Tan, M. V. Chekhova, and L. A. Krivitsky, “Accessing photon bunching with photon number resolving multi-pixel detector,” Opt. Express 19(10), 9352–9363 (2011), http://www.opticsinfobase.org/oe/abstract.cfm?uri=oe-19-10-9352. [CrossRef]

17. R. H. Hadfield, “Single-photon detectors for optical quantum information applications,” Nat. Photonics 3, 696–705 (2009). [CrossRef]

18. Z. L. Yuan, B. E. Kardynal, A. W. Sharpe, and A. J. Shields, “High speed single photon detection in the near infrared,” Appl. Phys. Lett. 91(4), 041114 (2007). [CrossRef]

19. B. E. Kardynal, Z. L. Yuan, and A. J. Shields, “An avalanche–photodiode–based photon–number–resolving detector,” Nat. Photonics 2(7), 425–428 (2008).

20. O. Thomas, Z. L. Yuan, J. F. Dynes, and A. J. Shields, “Efficient photon number detection with silicon avalanche photodiodes,” Appl. Phys. Lett. 97(3), 031102 (2010). [CrossRef]

21. L. Mandel and E. Wolf, “Quantum theory of photoelectric light detection,” in Optical Coherence and Quantum Optics (Cambridge University Press, 1995), pp. 683–740.

22. A. R. Dixon, J. F. Dynes, Z. L. Yuan, A. W. Sharpe, A. J. Bennett, and A. J. Shields, “Ultrashort dead time of photon-counting InGaAs avalanche photodiodes,” Appl. Phys. Lett. 94(23), 231113 (2009). [CrossRef]

23. A. W. Smith and J. A. Armstrong, “Laser photon counting distributions near threshold,” Phys. Rev. Lett. 16(25), 1169–1172 (1966). [CrossRef]

24. G. Lachs, “Theoretical aspects of mixtures of thermal and coherent radiation,” Phys. Rev. 138(4B), B1012–B1016 (1965). [CrossRef]

25. The value of g(2) ∼ 1.2 is corroborated experimentally by an independent measurement of g(2).

26. M. Aßmann, F. Veit, M. Bayer, M. van der Poel, and J. M. Hvam, “Higher-order photon bunching in a semiconductor microcavity,” Science 325(5938), 297–300 (2009). [CrossRef] [PubMed]

27. P. Meystre and M. Sargent, “Field quantization,” in Elements of Quantum Optics (Springer–Verlag, 1998), pp. 263–285.

28. E. Waks, E. Diamanti, and Y. Yamamoto, “Generation of photon number states,” New J. Phys. 8(1), 4 (2006). [CrossRef]

Probing higher order correlations of the photon field with photon number resolving avalanche photodiodes

Abstract

1. Introduction

2. Theory

3. Photon number resolution with InGaAs avalanche photodiodes

4. Higher order correlations of the photon field

5. Conclusion

Acknowledgments

References and links

Cited By

Figures (4)

Equations (6)

Optics Express