Photon statistics and bunching of a chaotic semiconductor laser
Yanqiang Guo Chunsheng Peng Yulin Ji Pu Li Yuanyuan Guo Xiaomin Guo






1Key Laboratory of Advanced Transducers and Intelligent Control System, Ministry of Education and Shanxi Province, Taiyuan University of Technology, Taiyuan 030024, China

2Institute of Optoelectronic Engineering, College of Physics and Optoelectronics, Taiyuan University of Technology, Taiyuan 030024, China

*guoxiaomin@tyut.edu.cn 







The photon statistics and bunching of a semiconductor laser with external optical feedback are investigated experimentally and theoretically. In a chaotic regime, the photon number distribution is measured and undergoes a transition from Bose-Einstein distribution to Poisson distribution with increasing the mean photon number. The second order degree of coherence decreases gradually from 2 to 1. Based on Hanbury Brown-Twiss scheme, pronounced photon bunching is observed experimentally for various injection currents and feedback strengths, which indicates the randomness of the associated emission light. Near-threshold injection currents and strong feedback strengths modify exactly the laser performance to be more bunched. The macroscopic chaotic dynamics is confirmed simultaneously by high-speed analog detection. The theoretical results qualitatively agree with the experimental results. It is potentially useful to extract randomness and achieve desired entropy source for random number generator and imaging science by quantifying the control parameters.


(140.1540) Chaos; (140.5960) Semiconductor lasers; (190.3100) Instabilities and chaos; (270.5290) Photon statistics. 



99
Soriano13 M. C. Soriano, J. García-Ojalvo, C. R. Mirasso, and I. Fischer, “Complex photonics: Dynamics and applications of delay-coupled semiconductors lasers,” 85(1), 421–470 (2013).

Sciamanna15 M. Sciamanna and K. A. Shore, “Physics and applications of laser diode chaos,” Nat. Photonics 9(3), 151–162 (2015).

Broom69 R. F. Broom, “Self modulation at gigahertz frequencies of a diode laser coupled to an external cavity,” Electron. Lett. 5(23), 571–572 (1969).

Lang80 R. Lang and K. Kobayashi, “External optical feedback effects on semiconductor injection laser properties,” 16(3), 347–355 (1980).

Torcini06 A. Torcini, S. Barland, G. Giacomelli, and F. Marin, “Low-frequency fluctuations in vertical cavity lasers: experiments versus Lang-Kobayashi dynamics,” 74(6), 063801 (2006).

Argyris05 A. Argyris, D. Syvridis, L. Larger, V. Annovazzi-Lodi, P. Colet, I. Fischer, J. García-Ojalvo, C. R. Mirasso, L. Pesquera, and K. A. Shore, “Chaos-based communications at high bit rates using commercial fibre-optic links,” 438(7066), 343–346 (2005).

Yoshimura12 K. Yoshimura, J. Muramatsu, P. Davis, T. Harayama, H. Okumura, S. Morikatsu, H. Aida, and A. Uchida, “Secure key distribution using correlated randomness in lasers driven by common random light,” 108(7), 070602 (2012).

Uchida08 A. Uchida, K. Amano, M. Inoue, K. Hirano, S. Naito, H. Someya, I. Oowada, T. Kurashige, M. Shiki, S. Yoshimori, K. Yoshimura, and P. Davis, “Fast physical random bit generation with chaotic semiconductor lasers,” Nat. Photonics 2(12), 728–732 (2008).

Reidler09 I. Reidler, Y. Aviad, M. Rosenbluh, and I. Kanter, “Ultrahigh-speed random number generation based on a chaotic semiconductor laser,” 103(2), 024102 (2009).

Wang13 A. B. Wang, P. Li, J. G. Zhang, J. Z. Zhang, L. Li, and Y. C. Wang, “4.5 Gbps high-speed real-time physical random bit generator,” 21(17), 20452–20462 (2013).

Lin04 F. Y. Lin and J. M. Liu, “Chaotic radar using nonlinear laser dynamics,” 40(6), 815–820 (2004).

Wang08 Y. C. Wang, B. J. Wang, and A. B. Wang, “Chaotic correlation optical time domain reflectometer utilizing laser diode,” IEEE Photon. Technol. Lett. 20(19), 1636–1638 (2008).

Brunner13 D. Brunner, M. C. Soriano, C. R. Mirasso, and I. Fischer, “Parallel photonic information processing at gigabyte per second data rates using transient states,” Nat. Communications 4(4), 1364 (2013).

Argyris10 A. Argyris, S. Deligiannidis, E. Pikasis, A. Bogris, and D. Syvridis, “Implementation of 140Gb/s true random bit generator based on a chaotic photonic integrated circuit,” 18(18), 18763–18768 (2010).

Oliver11 N. Oliver, M. C. Soriano, D. W. Sukow, and I. Fischer, “Dynamics of a semiconductor laser with polarization-rotated feedback and its utilization for random bit generation,” 36(23), 4632–4634 (2011).

Kanter10 I. Kanter, Y. Aviad, I. Reidler, E. Cohen, and M. Rosenbluh, “An optical ultrafast random bit generator,” Nat. Photonics 4(1), 58–61 (2010).

Locquet10 A. Locquet, “Chaos-Based Secure Optical Communications Using Semiconductor Lasers,” in Handbook of Information and Communication Security, M. Stamp and P. Stavroulakis, eds. (Springer, 2010), pp. 451–478.

Heiligenthal11 S. Heiligenthal, T. Dahms, S. Yanchuk, T. Jüngling, V. Flunkert, I. Kanter, E. Schöll, and W. Kinzel, “Strong and Weak Chaos in Nonlinear Networks with Time-Delayed Couplings,” 107(23), 234102 (2011).

Toomey14 J. P. Toomey and D. M. Kane, “Mapping the dynamic complexity of a semiconductor laser with optical feedback using permutation entropy,” 22(2), 1713–1725 (2014).

Porte14 X. Porte, M. C. Soriano, and I. Fischer, “Similarity properties in the dynamics of delayed-feedback semiconductor lasers,” 89(2), 023822 (2014).

Heil99 T. Heil, I. Fischer, W. Elsässer, J. Mulet, and C. R. Mirasso, “Statistical properties of low-frequency fluctuations during single-mode operation in distributed-feedback lasers: experiments and modeling,” 24(18), 1275–1277 (1999).

Sukow99 D. W. Sukow, T. Heil, I. Fischer, A. Gavrielides, A. Hohl-AbiChedia, and W. Elsäβer, “Picosecond intensity statistics of semiconductor lasers operating in the low-frequency fluctuation regime,” 60(1), 667–673 (1999).

Sacher89 J. Sacher, W. Elsässer, and E. O. Göbel, “Intermittency in the coherence collapse of a semiconductor laser with external feedback,” 63(20), 2224–2227 (1989).

Sano94 T. Sano, “Antimode dynamics and chaotic itinerancy in the coherence collapse of semiconductor lasers with optical feedback,” 50(3), 2719–2726 (1994).

Fischer96 I. Fischer, G. H. M. van Tartwijk, A. M. Levine, W. Elsässer, E. Göbel, and D. Lenstra, “Fast pulsing and chaotic itinerancy with a drift in the coherence collapse of semiconductor lasers,” 76(2), 220–223 (1996).

Li14 N. Li, B. Kim, A. Locquet, D. Choi, W. Pan, and D. S. Citrin, “Statistics of the optical intensity of a chaotic external-cavity DFB laser,” 39(20), 5949–5952 (2014).

Albert11 F. Albert, C. Hopfmann, S. Reitzenstein, C. Schneider, S. Höfling, L. Worschech, M. Kamp, W. Kinzel, A. Forchel, and I. Kanter, “Observing chaos for quantum-dot microlasers with external feedback,” Nat. Communications 2(1), 366 (2011).

Schulze14 F. Schulze, B. Lingnau, S. M. Hein, A. Carmele, E. Schöll, K. Lüdge, and A. Knorr, “Feedback-induced steady-state light bunching above the lasing threshold,” 89(4), 041801(R) (2014).

Hanbury56 R. Hanbury Brown and R. Q. Twiss, “Correlation between photons in two coherent beams of light,” 177(4497), 27–29 (1956).

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

Mandel95 L. Mandel and E. Wolf, Optical Coherence and Quantum Optics (Cambridge, 1995).

Arecchi65 F. T. Arecchi, “Measurement of the statistical distribution of Gaussian and laser sources,” 15(24), 912–916 (1965).

Kimble77 H. J. Kimble, M. Dagenais, and L. Mandel, “Photon antibunching in resonance fluorescence,” 39(11), 691–695 (1977).

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

Banaszek09 K. Banaszek, R. Demkowicz-Dobrzański, and I. A. Walmsley, “Quantum states made to measure,” Nat. Photonics 3(12), 673–676 (2009).

Davidovich96 L. Davidovich, “Sub-Poissonian processes in quantum optics,” 68(1), 127–173 (1996).

Schultheiss16 V. H. Schultheiss, S. Batz, and U. Peschel, “Hanbury Brown and Twiss measurements in curved space,” Nat. Photonics 10(2), 106-110 (2016).

Ryczkowski16 P. Ryczkowski, M. Barbier, A. T. Friberg, J. M. Dudley, G. Genty, “Ghost imaging in the time domain,” Nat. Photonics 10(3), 167–170 (2016).

Loaiza16 O. S. Magaña-Loaiza, M. Mirhosseini, R. M. Cross, S. M. Hashemi Rafsanjani, and R. W. Boyd, “Hanbury Brown and Twiss interferometry with twisted light,” Sci. Adv. 2(4), e1501143 (2016).

Guo12 Y. Q. Guo, G. Li, Y. F. Zhang, P. F. Zhang, J. M. Wang, T. C. Zhang, “Efficient fluorescence detection of a single neutral atom with low background in a microscopic optical dipole trap,” Sci. China: Phys. Mech. Astron. 55(9), 1523–1528 (2012).

Redding12 B. Redding, M. A. Choma, and H. Cao, “Speckle-free laser imaging using random laser illumination,” Nat. Photonics 6(6), 355–359 (2012).

Hayenga16 W. E. Hayenga, H. Garcia-Gracia, H. Hodaei, C. Reimer, R. Morandotti, P. LiKamWa, and M. Khajavikhan, “Second-order coherence properties of metallic nanolasers,” Optica 3(11), 1187–1193 (2016).

Hagerstrom15 A. M. Hagerstrom, T. E. Murphy, and R. Roy, “Harvesting entropy and quantifying the transition from noise to chaos in a photon-counting feedback loop,” PNAS 112(30), 9258–9263 (2015).

Li07 Y. Li, G. Li, Y. C. Zhang, X. Y. Wang, J. Zhang, J. M. Wang, and T. C. Zhang, “Effects of counting rate and resolution time on a measurement of the intensity correlation function,” 76(1), 013829 (2007).

Kondakci15 H. E. Kondakci, A. F. Abouraddy, and B. E. A. Saleh, “A photonic thermalization gap in disordered lattices,” Nat. Physics 11(11), 930–935 (2015).






§ INTRODUCTION


Semiconductor lasers (SCL) with external optical feedback (EOF) provide a rich platform to study the nonlinear effects and complex photonics dynamics <cit.>. EOF is the most prominent configuration for achieving laser chaos. A common aspect of such systems is the presence of chaotic behaviour that consists of very short and random spiking of the laser intensity. The chaotic behaviour in these systems is extremely sensitive to time-delayed returning field of their own emission and was identified as early as the late 1960s <cit.>. To verify the observed chaotic dynamics, the Lang-Kobayashi (LK) model <cit.> is employed to describe the system qualitatively, and even quantitatively in some cases <cit.>. With the growing understanding of the chaotic dynamics and its control, the optical-feedback-induced laser chaos has gained more and more attention for chaos-based secure communications <cit.>, chaos key distribution <cit.>, physical random number generation (RNG) <cit.>, chaotic optical sensing <cit.>, and bioinspired information processing <cit.>. Thus, it is indispensable to understand and elucidate the underlying mechanisms, the nonlinear dynamics under different conditions, and how to distinguish and identify different chaotic properties. Intensity statistics and autocorrelation (AC) associated to the dynamics of SCL with EOF are important and useful to characterize chaotic processes. The intensity statistics is closely related to the maximum extractible rate of randomness <cit.> and post-processing techniques <cit.> in RNG application. In optical chaos communications, the AC can be used as a good indicator of a bandwidth of chaotic laser limiting the speed of modulating chaotic carrier <cit.>. Additionally, the AC function can present useful information about the weak-strong chaos transition in the laser system with EOF <cit.>. Previous research has been mainly devoted to elucidating the intensity statistics of chaotic laser that operates in the low-frequency fluctuations (LFF) regime <cit.>. A typical characteristic of LFF is a sudden power drop followed by a gradual power recovery and the frequency of LFF is distinctly low compared to the laser's intrinsic relaxation oscillation frequency <cit.>. This phenomenon occurs near the laser threshold for low-to-moderate feedback and has been explained as chaotic itinerancy with a drift <cit.>. The chaotic itinerancy predicts fast chaotic pulsations of the output intensity. Since the pulsing behavior skews the intensity distribution, the probability density distributions of the laser intensity becomes extremely asymmetric <cit.>. As the pump current and feedback strength are increased, the laser linewidth is broadened greatly from a few MHz to tens of GHz and the LFF typically evolves into the laser's coherence collapse (CC) <cit.>. The research on the intensity statistics of high-dimensional chaotic waveforms in the CC regime is sparse, in contrast to that of LFF dynamics <cit.>. However, it is still unclear whether the intensity statistics of chaotic laser would change when varying the injection current and feedback strength. Currently, there is a significant discrepancy between experimental and theoretical probability density distributions <cit.>. Recent researches reveal that the photon statistics is a candidate to indicate random-intensity fluctuations associated with the emission of chaotic laser <cit.>. Photon number distribution and photon correlations are more sensitive to control parameters compared to probability density distributions of laser intensity <cit.>. Photon statistical investigation of chaotic laser is also rare in the continuous transition region between the LFF regime and the fully developed CC chaotic regime.

 The photon statistics and correlations are pivotal and fundamental in characterizing the quantum statistics of light sources. The pioneering experiments conducted by Hanbury Brown and Twiss (HBT) were a landmark linking the temporal and spatial second-order degree of coherence g^(2)(τ ,x) of a thermal source <cit.>. The g^(2)(τ ,x) was formalized by Glauber in the 1960s <cit.> and described by correlation function within statistical optics <cit.>, which has information on both the nature and dynamics of the photon emission process of the underlying field. Subsequently, the observation of g^(2) effect is harnessed in numerous experiments and applications, such as measuring photon bunching and photon number distribution <cit.>, characterizing nonclassicality of light field <cit.>, and so on. With the rapid development of photon counting technology, single photon counting detection, as the most sensitive and very widespread method of optical measurement, has been expanded towards quantum information <cit.> and precision measurement applications <cit.>. Up until now, g^(2)(τ) already allows one to categorize different light fields and carry a great deal of information on the photon statistical properties of a light field. In general, g^(2)(0) (i.e., at zero time difference) >1 or Bose-Einstein distribution, characteristic of an incoherent or chaotic light, while g^(2)(0)=1 or Poissonian distribution, characteristic of a coherent light or a classical stable field, and g^(2)(0)<1 or sub-Poissonian distribution, characteristic of a quantum light emission <cit.>. In this regard, g^(2)(τ ) is fundamentally different from the first-order degree of coherence, which serves as a description of phase and cannot sufficiently extract quantum properties of light field. Thus, HBT scheme combined with single-photon detection have boosted the study of spatial interference <cit.>, ghost imaging <cit.>, azimuthal HBT effect <cit.>, deterministic manipulation and detection of single-photon source <cit.>, etc. Moreover, the random or chaotic properties of light is an essential element of the HBT effect and have been applied in full-field imaging science <cit.> and nanophotonics <cit.>. However, quantum statistics of chaotic laser is unexplored. More recent experiments of photon statistics in chaotic regime aim to bridge the gap between chaotic laser and quantum optics <cit.>. A comprehensive understanding of photon statistics and correlation of chaotic laser remains an open question.

In this paper, we present the experimental and numerical results of photon statistics and coherence of chaotic laser based on single photon counting technique. The chaotic laser operating in the mW (high-gain) regime, consists of SCL with time-delayed optical feedback. In order to compare macroscopic dynamics and photon statistical transition, we simultaneously measure the optical signals by photon-counting method and high-speed analog detection. Implementing HBT interferometry, we demonstrate that the control parameters of chaotic laser, i.e., feedback strength and injection current, can substantially affect the photon correlation. The photon number distribution of the chaotic laser undergoes a transition from Bose-Einstein distribution to Poissonian distribution as the mean photon number increases. It shows a good agreement between experimental photon number distribution and theoretical fitting without post-processing. In the chaotic regime, the laser emission shows bunching effect, and the control parameters modify exactly the laser performance to be more bunched. Lower injection current and stronger feedback strength are beneficial for enhancing bunching effect. The macroscopic chaotic dynamics is confirmed experimentally. The phenomenon is also verified numerically using the LK equations where the results are in good agreement with the experimental data. This demonstration well reveals photon statistics and coherence of chaotic laser and provide a necessary and better understanding of chaotic process. In this sense, it contributes to the research of chaos with quantum optics.



§ EXPERIMENTAL SETUP AND THEORETICAL MODEL


The experimental setup is schematically illustrated in Fig. 1. The chaotic laser is composed of a distributed feedback laser diode (DFB-LD, WTD LDM5S752) subject to external optical feedback, operating at 1.55 μm and with threshold current J_th of 10.5 mA. The temperature and current of DFB-LD are stabilized to an accuracy of 0.01 K and 0.1 mA, respectively. The DFB-LD connects to an 80:20 optical coupler (OC1) whose principal output passes through a variable optical attenuator (VOA1) onto a optical circulator, forming a fiber-based feedback cavity with a time delay of 99.85 ns. The VOA1 and polarization controller (PC) are used to accurately define the feedback conditions. The optical feedback drives the DFB-LD into high-frequency chaotic oscillation. The output of chaotic laser is split via 50:50 optical coupler (OC2). At one port photon statistics and coherence are measured based on a HBT scheme: the laser emission is coupled out and collimated by triplet lens coupler. The freespace beam goes through a filter and is divided intensity-equally into two parts. The transmitted and reflected photons eventually arrive at dual-channel single photon detector (DSPD, Aurea Technology LYNXEA-NIR-M2-SM-01). The quantum efficiency of DSPD is 25% at 1550 nm and the counting rate is controlled below 0.4 Mcounts/s by a variable optical attenuator (VOA2). The DSPD outputs are then connected to a time-correlated coincidence unit and the time difference between the two signals originating from two emitted photons is repeatedly measured with 60 ps resolution. The other port of OC2 is used for detection of intensity and frequency dynamics by two 50 GHz photodetectors (PD, Finisar XPDV2120RA). All data are recorded simultaneously using a 26.5 GHz RF spectrum analyzer (SA, Agilent N9020A, 3 MHz RBW, 3 KHz VBW) and a 40 Gsamples/s real-time oscilloscope (OSC, Lecroy, LabMaster10-36Zi) with 36 GHz bandwidth.



In the high-gain regime, chaotic dynamics of a single-mode semiconductor laser with time delayed optical feedback can be well modeled by the LK equations, which are capable of confirming and reproducing the observations <cit.>. The LK equations for the three real variables of the optical field amplitude or the photon flux density E(t), the field phase φ (t), and the carrier density N(t) are described as follows:

    Ė(t)=1/2[G(t)-τ _p^-1)E(t)+κ E(t-τ _ext)cos
    [ϕ (t)],


    φ̇(t)=α/2[G(t)-τ _p^-1]-κ E(t-τ
    _ext)E(t)^-1sin [ϕ (t)],


    Ṅ(t)=J/e-N(t)/τ _N-G(t)| E(t)|
    ^2,


    ϕ (t)=ωτ +φ (t)-φ (t-τ _ext),

where G(t)=G_N[N(t)-N_0]/(1+ε| E(t)| ^2) is the nonlinear optical gain (with G_N being the gain coefficient and ε the saturation coefficient), N_0 is the carrier density at transparency, κ =(1-r_in^2)r_0/(r_inτ _in) is the optical feedback strength, r_in is the reflectivity of the internal cavity, r_0 is the reflectivity of the external mirror, τ _in is the optical round-trip time in internal cavity, τ _p is the photon lifetime, τ _N is the carrier lifetime, α is the linewidth-enhancement factor, τ _ext is the feedback delay time of the external cavity, ω =2π c/λ is the angular optical frequency, c is the speed of light, and λ is the optical wavelength, and J=ρ J_th is the injection current density (with ρ being the pump factor). The laser threshold current results from Eq. (3) as J_th=(1/τ_N)[N_0+1/(G_Nτ _p)]≈ 10.5 mA, corresponding to J_th of the DFB laser. The following values for the above parameters are fixed at α =5, τ _p=2.5 ps, τ _N=2.3 ns, G_N=2.56× 10^-8 ps^-1, N_0=1.35× 10^8, τ _ext=99.85 ns, λ =1.55 μ m, ε =5× 10^-7. We numerically integrated the Eqs. (1)–(3) with a fourth order Runge-Kutta routine by employing time steps of h=2× 10^-12 s.

For the initial conditions, Eqs. (1)–(3) are solved many times (typically 10^5), each solution being referred to as a “realization” of the laser operation. The second-order degree of coherence g^(2)(τ ) of optical field can be calculated as

    g^(2)(τ )=⟨ E^∗(t)E^∗(t+τ )E(t+τ
    )E(t)⟩/⟨ E^∗(t)E(t)⟩ ^2=
    ∫_tI(t)I(t+τ )dt/(∫_tI(t)dt)^2=⟨
    n_ph(t)n_ph(t+τ )⟩/⟨ n_ph(t)⟩
    ^2,

where I(t)=| E^∗(t)E(t)| denotes the laser intensity, τ is the time difference between the two photon detection events, n_ph(t) is the photon number and ⟨·⟩ designates the statistical averaging that is done over a large ensemble of different realizations of the laser field. For coherent light, the photon number distribution P(n_ph) obeys Poisson distribution P(n_ph)=⟨ n_ph⟩ ^n_phe^-⟨ n_ph⟩/n_ph!, where ⟨ n_ph⟩ =∑_n_phn_phP(n_ph) is the mean photon number. For chaotic light, the photon number distribution P(n_ph) obeys Bose-Einstein distribution P(n_ph)=⟨ n_ph⟩ ^n_ph/(1+⟨ n_ph⟩ )^n_ph+1. Figure 2(a) shows the calculated emission intensity for J=1.5J_th, κ =35 ns^-1, κ =50 ns^-1, κ =65 ns^-1, and the parameter values mentioned above. For κ =35 ns^-1, κ =50 ns^-1, and κ =65 ns^-1 the model predicts chaotic intensity fluctuations. The related photon number distributions are depicted in Fig. 2(b). The histograms from the LK model are better fitted by Bose-Einstein distribution (blue solid curve), whereas the Poisson distribution (black dotted curve) shows relatively poor fit. The chaotic intensity fluctuations and the spiked emission of photons result in bunching effect with g^(2)(0)>1. As the optical feedback strength κ increases the bunching effect becomes strong, and the associated g^(2)(τ ) results are shown in Fig. 2(c).

For κ =35 ns^-1, the central bunching maximum is observed with 1<g^(2)(0)<2 that is associated with sub-chaotic light. For κ =50 ns^-1, the bunching maximum is obtained with a magnitude g^(2)(0)∼ 2, which indicates that the light approaches fully chaotic limit (g^(2)=2), or chaotic limit. For strong optical feedback strength [κ =65 ns^-1 in Fig. 2(c)], the pronounced bunching maximum is observed with a magnitude g^(2)(0) substantially exceeding 2 and the light becomes super chaotic, which is more chaotic than the fully chaotic light mentioned above. Owing to the laser's intrinsic relaxation oscillation, slight oscillations of g^(2)(τ ) appear around the central bunching area.



§ RESULTS


We start by evaluating the photon statistics of chaotic laser. Figures 3(a) and 3(b) show the measured time traces of the chaotic laser and its associated RF spectrum operating at J=1.5J_th and experimental feedback strength η =12.8%. In our experiment, η =25% of the optical power is fed back for κ =20 ns^-1 and we keep it below. Comparing with the noise floor, the output of the laser shows high bandwidth and large amplitude fluctuation, which indicate that the laser is operating in a chaotic region. In Fig. 3(a), the gray line (red for online version) is noise floor from the stable operation laser. According to the 80% bandwidth definition, the bandwidth of the chaotic laser is about 9.85 GHz, as shown in Fig. 3(b).



Under this condition, the measured photon number distribution P(n_ph) and second-order degree of coherence g^(2)(0) of chaotic laser are illustrated in Fig. 4. The incident photon number is controlled by finely adjusting the programmable VOA2. With an increase of the mean photon number, the photon number distribution changes continuously from Bose-Einstein distribution to Poisson distribution and the g^(2)(0) decreases gradually from 2 to 1. The maximum value of measured g^(2)(0) is 2.02 and the associated photon number distribution is almost identical to the Bose-Einstein distribution of the same mean photon number (⟨ n_ph⟩ =0.69). As we increase the mean photon number, the photon number distribution of chaotic laser gradually deviates from the Bose-Einstein distribution. When the mean photon number ⟨ n_ph⟩ reaches 1.8, the measured photon number distribution is between the Bose-Einstein distribution and Poisson distribution and the g^(2)(0) becomes 1.21. As the mean photon number is increased further, the photon number distribution gets closer to the Poisson distribution and the g^(2)(0) is reduced. Eventually, when the mean photon number  ⟨ n_ph⟩ is 2.61, the photon number distribution is nearly identical to the Poisson distribution and the corresponding g^(2)(0) is 1.03. Due to the fact that the multiphoton events increase, the photon arriving interval becomes shorter compared to the sampling time. From another point of view, the sampling time gets longer than the coherence time of the emitted photons. Thus, the measured photon number distribution gradually approaches Poisson distribution. Without any post-processing, the experimental results are in good agreement with the theoretical fitting. The accuracy is improved compared to the measurement of macroscopic intensity statistics <cit.>. The measured photon statistics of the light field also can be profoundly affected by experimental conditions. In this work, we choose the proper experimental conditions (e.g., below 1.62 for the mean photon number and 80 Counts/s for the photon count rate) and the results can reflect correctly the photon statistics of light field <cit.>.



Despite its utility, the photon number distribution P(n_ph) is a single description and does not constitute a complete representation of the light. A unique description requires knowledge of the photon correlation g^(2), which we have shown in Fig. 5 for various injection current J and feedback strength κ in the chaotic regime. The g^(2)(0) also characterizes the randomness of the associated photon number n_ph and is a measurement of photon bunching, which refers to the tendency of photons to arrive together. As the g^(2)(0) increases, the photons become more bunched and chaotic <cit.>. Figure 5 shows theoretical 5(a) and experimental 5(b) g^(2)(0) as functions of κ (η ) for five different values of J in chaotic regime: 1.07J_th, 1.12J_th, 1.2J_th, 1.5J_th, and 2.0J_th. At all the set J values, the g^(2)(0) gradually grows as κ (η ) increases. The κ (η ) affects little on the threshold J_th, and there is a negligible influence on photon correlation measurement. The experimental results shows good agreement with the theory. Additionally, the bunching effect become weak as the inject current J increases. In Figs. 5(c) and 5(d), the theory and experiment for g^(2)(0) at various J/J_th and κ (η ) are shown. At all κ (η ) shown, the g^(2)(0) decreases monotonically as the current J increases and the theory is in good agreement with the experiment. For larger J, more multiphoton events result in the g^(2)(0) falls faster, leading to approach the stable laser output. Near-threshold injection currents and strong feedback strengths are beneficial for enhancing bunching effect. The pronounced photon bunching is helpful to boost the image contrast in ghost imaging. Moreover, the bunched chaotic laser with low coherence is well suited for speckle-free full-field imaging <cit.>. The results also provide a basic picture and understanding of the higher order correlation of chaotic light in the single photon counting regime.



In what follows, we confirm the chaotic dynamics to help interpret the near-threshold bunching effect. Figure 6 depicts the measured temporal waveform and RF spectrum for three values of η =3.1%, η =12.5%, and η =25%, when J=1.07J_th [Figs. 6(a1)-6(a4)], J=1.2J_th [Figs. 6(b1)-6(b4)], and J=1.5J_th [Figs. 6(c1)-6(c4)]. At all the three J values, the amplitude fluctuation and bandwidth of chaotic laser increase as the feedback strength η grows. From Fig. 6(a4) to Fig. 6(c4), the bandwidth of the laser is broadened from 2.65 GHz to 9.9 GHz as the injection current and feedback strength increase. In the near-threshold regime, wider bandwidth and larger amplitude fluctuation contribute to obtaining stronger bunching effect when the injection current is fixed. While the feedback strength is fixed, it's an opposite case that stronger bunching effect happens on the bandwidth and amplitude fluctuation drop.



In addition, the external-cavity delay provide useful information on the chaotic dynamics and indicate the weak-strong chaos transition <cit.>. The maximum height h of normalized autocorrelation function near the first delay echo is inversely proportional to the chaotic strength. Figure 7 illustrates theory 7(a) and experiment 7(b) for h at various J and κ (η ). At all J shown, the h follows a nonmonotonic dependence on the feedback strength κ (η ), revealing a dip for intermediate κ (η ). The theory is in good agreement with experiment. In the week feedback regime, the chaotic strength and bunching effect are both enhanced as κ (η ) increases. On the other side (strong feedback), the chaotic laser becomes more bunched but the chaotic strength decreases when increasing κ (η ).





§ CONCLUSIONS


To conclude, the photon statistics and second-order photon correlation g^(2) are investigated theoretically and experimentally for a SCL with optical feedback and current modulation in the chaotic regime. The photon number distribution of the laser operating in a chaotic regime transitions continuously from Bose-Einstein distribution to Poisson distribution as the mean photon number increases. The theoretical results predict the pronounced photon bunching which is observed in HBT experiment. We have found a good agreement between the theoretical model and experimental observations for the g^(2)(0) as functions of J and κ (η ). To describe the photon statistics and bunching of SCLs under optical feedback in the high-gain regime, the LK model is well employed. The bunching effect is enhanced by increasing the feedback strength near the threshold current. The results also provide a better understanding of the super-chaotic bunching phenomena. The new type of photon bunching spectroscopy with the current and optical feedback detuning allow us to extract information about photon statistics and randomness. In this sense, the finding will boost the study of ghost imaging utilizing chaotic laser with tunable photon correlation and speckle-free full-field imaging with low-coherence chaotic laser, and contribute to improvements of chaos-based random number generation and secure communication with chaos synchronization.



§ FUNDING


National Natural Science Foundation of China (NSFC) (61405138, 61505136, 61775158, 61671316); Shanxi Scholarship Council of China (SXSCC) (2017-040); Natural Science Foundation of Shanxi Province (201701D221116).



§ ACKNOWLEDGMENTS


The authors thank H. Shen for helpful discussions.

