Limits on manipulating conditional photon statistics via interference of weak lasers

TL;DR

This paper theoretically investigates the limits of manipulating photon statistics through interference of weak lasers, showing that while super-Poissonian statistics are achievable, sub-Poissonian anti-bunching cannot be obtained with this method.

Contribution

It clarifies the fundamental limits of using interference of weak lasers for controlling photon statistics, correcting misconceptions about achieving anti-bunching.

Findings

Conditional photon statistics can be super-Poissonian.

Photon anti-bunching cannot be achieved with interference of weak lasers.

Proper accounting of higher-order photon terms is crucial.

Abstract

Photon anti-bunching, measured via the Hanbury-Brown-Twiss experiment, is one of the key signatures of quantum light and is tied to sub-Poissonian photon number statistics. Recently, it has been reported that photon anti-bunching or conditional sub-Poissonian photon number statistics can be obtained via second-order interference of mutually incoherent weak lasers and heralding based on photon counting. Here, we report theoretical analysis on the limits of manipulating conditional photon statistics via interference of weak lasers. It is shown that conditional photon number statistics can become super-Poissonian in such a scheme. We, however, demonstrate explicitly that it cannot become sub-Poissonian, i.e., photon anti-bunching cannot be obtained in such a scheme. We point out that incorrect results can be obtained if one does not properly account for seemingly negligible higher-order photon number expansions of the coherent state.

Equations

$$|\Psi {⟩}_{\mathrm{in}}=e^{-\frac{\alpha }{2}}\sum _{\lambda =0}^{\infty }\frac{{\alpha }^{\frac{\lambda }{2}}}{\sqrt{\lambda !}}|\lambda ⟩,$$

$$|\Psi ⟩$$

$$|\Psi ({\omega }_a,{\omega }_b)⟩$$

$$\rho =\int d{\omega }_{a}d{\omega }_b\mathcal{F}({\omega }_a)\mathcal{F}({\omega }_b)|\Psi ({\omega }_a,{\omega }_b)⟩⟨\Psi ({\omega }_a,{\omega }_b)|,$$

$$\rho$$

$${\hat{a}}_H^{†}({\omega }_a)$$

$${\hat{b}}_V^{†}({\omega }_b)$$

$${\rho }_{cd}$$

$$R_{cd}(\tau )=\frac{\text{Tr}\left[{\rho }_{cd}{E}_c^{(-)}(\tau )E_d^{(-)}(0)E_d^{(+)}(0)E_c^{(+)}(\tau )\right]}{\text{Tr}\left[{\rho }_{cd}{E}_c^{(-)}(\tau )E_c^{(+)}(\tau )\right]\text{Tr}\left[{\rho }_{cd}{E}_d^{(-)}(0)E_d^{(+)}(0)\right]},$$

$$g_C^{(2)}(\tau ,{\tau }_c)=\frac{{\sum }_{p=1}^{\infty }{I}_d^{(p)}{\sum }_{p=1}^{\infty }{I}_{def}^{(p)}(\tau ,{\tau }_c)}{{\sum }_{p=1}^{\infty }{I}_{de}^{(p)}(\tau ){\sum }_{p=1}^{\infty }{I}_{df}^{(p)}(\tau +{\tau }_c)}.$$

$$I_d^{(p)}=\text{Tr}\left[{\rho }_{cd}\{E_d^{(-)}(0){\}}^p\{E_d^{(+)}(0){\}}^p\right].$$

$$I_{de}^{(p)}(\tau )$$

$$I_{df}^{(p)}(\tau )$$

$$I_{def}^{(p)}(\tau ,{\tau }_c)$$

$$g_{C,p=3}^{(2)}(\tau ,0)=\mp \frac{\left(129{\alpha }^2+176\alpha +128\right)\left(2e^{\frac{{\sigma }^2{\tau }^2}{2}}\mp 3e^{{\sigma }^2{\tau }^2}\right)}{4{\left((11\alpha +8)e^{\frac{{\sigma }^2{\tau }^2}{2}}\mp 6\alpha \mp 4\right)}^2},$$

$$g_{C,p=4}^{(2)}(\tau ,0)=\frac{2\left(5047{\alpha }^3+9288{\alpha }^2+12672\alpha +9216\right)e^{-2{\sigma }^2{\tau }^2}\left(\alpha e^{{\sigma }^2{\tau }^2}-16(3\alpha +2)e^{\frac{5{\sigma }^2{\tau }^2}{2}}+6(11\alpha +8)e^{3{\sigma }^2{\tau }^2}\right)}{9{\left(\left(258{\alpha }^2+352\alpha +256\right)e^{\frac{{\sigma }^2{\tau }^2}{2}}-155{\alpha }^2-192\alpha -128\right)}^2}.$$

$$g_{C,p=5}^{(2)}(\tau ,0)=\frac{9}{128}$$

$$g_{C,p=6}^{(2)}(\tau ,0)=\frac{64}{25}$$

$$g_{C,p=7}^{(2)}(\tau ,0)=\frac{25}{72}$$

$$g_{C,p=8}^{(2)}(\tau ,0)=\frac{36}{49}$$

$$g_{C,p=9}^{(2)}(\tau ,0)=\frac{49}{2048}$$

$$g_{C,p=10}^{(2)}(\tau ,0)=\frac{1024}{81}$$

Full text

I Introduction

The normalized second-order correlation function $g^{(2)}(\tau)$ plays an important role in quantum optics by enabling to distinguish different kinds of light and is often measured via the Hanbury-Brown-Twiss experiment hbt ; glauber63 . Most importantly, it allows to distinguish different kinds of light. The chaotic light and the coherent light are associated with $g^{(2)}(0)\geq 1$ and they can be described well with the classical electromagnetic theory of light. Any light which exhibits $g^{(2)}(0)<1$, i.e., photon anti-bunching, is considered quantum as quantization of the electromagnetic field or the concept of photons is essential to describe its behaviors. As the coherent state exhibits Poissonian photon number statistics, light which exhibits photon anti-bunching is tied to sub-Poissonian photon number statistics loudon . There are indeed a wide variety of quantum light, such as, the Fock states hong86 ; takeuchi04 ; baek08 , multi-photon entangled states yao12 ; jeong16 ; zhang16 ; wang16 ; chen17 , squeezed states breitenbach97 ; zhang15 , macroscopic superposition states sander92 ; ourj07 , etc., and applications of quantum light includes quantum communication gisin07 , quantum computing kok07 , quantum metrology giov11 , etc.

A particular quantum state of light may be post-selected yao12 ; jeong16 ; zhang16 ; wang16 ; chen17 ; jeong13 ; lee14 or heralded takeuchi04 ; baek08 ; barz10 ; kim11 ; ra15 ; ra16 . In the heralding scheme, the statistical properties of the quantum state is conditioned by the heralding signal. One of the earliest experiments on a localized single-photon state relied on the heralding signal from a single-photon detection event of a two-photon state of spontaneous parametric down-conversion hong86 . In this example, the heralding signal causes the conditional photon number statistics to be sub-Poissonian, exhibiting photon anti-bunching, whereas unheralded photon number statistics would be that of the chaotic light, i.e., super-Poissonian exhibiting photon bunching. In fact, quantum state heralding is a powerful tool in preparing a complex quantum state of light and heralding schemes have been shown to generate a variety of non-classical light states, including various entangled states barz10 ; kim11 ; ra15 ; ra16 , photon added/subtracted states bellini , etc. The idea of heralding has also been expanded to heralding quantum processes such as, quantum storage of light tanji09 , quantum gates and operations pittman03 ; xiang10 ; okamoto11 , etc.

Recently, it has been reported that photon anti-bunching or conditional sub-Poissonian photon number statistics can be obtained via second-order interference of mutually incoherent weak lasers and heralding based on photon counting silva15 ; silva16 ; amaral16 . Even with weak lasers at the single-photon regime, it is well-known that linear optical elements do not change the photon number statistics and effects such as interference and phase randomization do not produce quantum light, although they may produce chaotic light which exhibits photon bunching or super-bunching cho10 ; hong12 ; pandey14 ; choi17 . Then, the question becomes whether conventional single-photon detectors can indeed herald a non-classical light state exhibiting photon anti-bunching from the classical input light. In particular, what is the limit on manipulating conditional photon number statistics via interference of weak lasers and heralding based on photon counting?

In this paper, we report the theoretical analysis on the limits of manipulating conditional photon statistics via interference of weak lasers and heralding based on photon counting. It is shown that conditional photon number statistics can become super-Poissonian in such a scheme. We demonstrate explicitly however that, contrary to Ref. silva15 ; silva16 ; amaral16 , it cannot become sub-Poissonian, i.e., photon anti-bunching cannot be obtained in such a scheme. Theoretical and numerical analyses show that such incorrect results can be obtained if one does not properly account for seemingly negligible higher-order photon number expansions of the coherent state even at the single-photon regime.

II Experimental scheme

Consider first the well-known Shih-Alley/Hong-Ou-Mandel experiment in which two single-photons impinge on a symmetric beam splitter via the two different input ports hom ; shih . As single-photons have no definite phase, no first-order interference is formed at the output of the beam splitter. If the two photons are distinguishable, they exit the beam splitter randomly. When coincidence counts are measured between the two detectors placed at each output port of the beam splitter, random coincidence events are measured. However, if the two single-photons are made to be indistinguishable, second-order quantum interference causes the photons to coalescence. As a result, the two photons are always found together at the same output port of the beam splitter, causing the change of the photon number statistics. In this case, null coincidence counts are measured due to the second-order quantum interference. The visibility, defined as the random coincidence counts subtracted by the coincidence counts due to quantum interference normalized to the random coincidence counts, in this case can reach the maximum value of one.

The experimental setup to manipulate conditional photon statistics via second-order interference of mutually incoherent weak lasers and heralding based on photon counting is inspired by the above mentioned Shih-Alley/Hong-Ou-Mandel experiment. Instead of single-photon states at the input of a symmetric beam splitter, we have two mutually incoherent weak lasers at the two input ports, see Fig. 1. The setup is similar to the experiments in Ref. choi17 ; yongsu ; yongsu2 ; silva15 ; silva16 ; amaral16 but we make use of the additional polarization degree of freedom to more clearly present the idea. Two beams of lasers $a$ and $b$ are prepared by splitting a laser beam with a beam splitter (BS). The mutual phase coherence between the two beams $a$ and $b$ is removed by using the acousto-optic modulators (AOM1 and AOM2) that are driven by independent RF sources yongsu ; yongsu2 . This is essential to ensure that there is no first-order interference between the two beams $a$ and $b$ as we overlap the two beams at the second BS. This configuration corresponds to the Shih-Alley/Hong-Ou-Mandel experiment but with two classical light beams hom ; shih ; yongsu ; yongsu2 .

The second-order intensity correlation between the two output ports of the BS, $c$ and $d$, can be manipulated by changing the interference condition. Figure 1(a) shows a typical experimental setup to measure the second-order intensity correlation $R_{cd}(\tau)$ between the modes $c$ and $d$. Here, $\tau$ denotes the time delay between two single-photon detection events at $D_{c}$ and $D_{d}$. The interference condition can be easily changed by the polarization states of beams. In order to implement various interference conditions, the polarization states of the laser beams at $a$ and $b$ are set to be horizontal ($|H\rangle_{a}$) and vertical ($|V\rangle_{b}$), respectively. Then, the interference condition can be manipulated by changing the combinations of the polarization projection bases at the output modes $c$ and $d$ by using the polarizers (Pol). When the projection basis $\{|D\rangle_{c},|D\rangle_{d}\}$, where $|D\rangle=\frac{1}{\sqrt{2}}(|H\rangle+|V\rangle)$, is chosen for the polarizers, the second-order intensity correlation $R_{cd}(0)$ reaches the minimum value. For the projection basis $\{|A\rangle_{c},|D\rangle_{d}\}$, where $|A\rangle=\frac{1}{\sqrt{2}}(|H\rangle-|V\rangle)$, the second-order intensity correlation $R_{cd}(0)$ reaches the maximum value. We thus consider these bases in our analysis as they offer the theoretical maximum visibility of the second-order interference choi17 ; yongsu ; yongsu2 .

Since the value of the second-order intensity correlation $R_{cd}(\tau>\tau_{\rm coh})=1$, where $\tau_{\rm coh}$ is the coherence time of the input light, if the input light is two mutually incoherent lasers, $R_{cd}(0)=0.5$ for the projection basis $\{|D\rangle_{c},|D\rangle_{d}\}$. Thus, the visibility of the second-order classical interference in the Shih-Alley/Hong-Ou-Mandel setup is limited by $V\leq 0.5$ choi17 ; yongsu ; yongsu2 . This result nevertheless signals that photons are weakly bunched. On the other hand, if the measurement basis $\{|A\rangle_{c},|D\rangle_{d}\}$ is chosen, $R_{cd}(0)=1.5$ indicating that the photons tend to distribute themselves in different spatial modes. This simple argument allows us to ask what the limit would be for manipulating conditional photon statistics in a particular output mode $c$ by using photon counting at the other output mode $d$ as the heralding signal. In particular, we are interested in the conditional second-order intensity correlation function for mode $c$, $g^{(2)}_{C}(\tau,\tau_{c})$, heralded by the photon counting signal at mode $d$. The relevant experimental setup to measure $g^{(2)}_{C}(\tau,\tau_{c})$ is shown in Fig. 1(b). Conditioned on the photon counting event at mode $d$, we use the Hanbury-Brown–Twiss setup to measure $g^{(2)}_{C}(\tau,\tau_{c})$, the conditional second-order intensity correlation function at mode $c$.

III Conditional photon statistics

In this section, we show the state evolution following the experimental setup in Fig. 1. The input to the first BS is a weak laser and its quantum state can be written in the Fock basis as,

[TABLE]

where $|\lambda\rangle$ denotes the $\lambda$-photon Fock state. The symmetric 50/50 beam splitter splits the incoming beam into two, each beam having the mean photon number corresponding to $\alpha/2$. The quantum state of light at the output of the first BS is thus given by

[TABLE]

where $|m\rangle_{a}$ and $|n\rangle_{b}$ represent the $m$ and $n$ photon Fock states, respectively. The photon creation operators at modes $a$ and $b$ are denoted as $\hat{a}^{\dagger}$ and $\hat{b}^{\dagger}$, respectively.

At each output of the first BS, an acousto-optic modulator (AOM) is placed. The AOMs are driven by independent and unsynchronized RF sources, making the incoming light beams mutually incoherent. Therefore, after the AOMs, the state of each frequency mode can be represented as

[TABLE]

where $\omega_{a}(\omega_{b})$ and $\gamma_{\omega_{a}}(\gamma_{\omega_{b}})$ are the frequency mode and the phase given by the AOM1 (AOM2) at mode $a~{}(b)$, respectively. The state of light beams at the input mode $a$ and $b$ to the second BS thus can be written as,

[TABLE]

where $\mathcal{F}(\omega_{a})$ and $\mathcal{F}(\omega_{b})$ are the frequency spectra of the light given by AOM1 and AOM2, respectively. The quantum description of the light beams at the input modes $a$ and $b$ to the second BS can then be fully written as,

[TABLE]

Note that the photons in mode $a$ and $b$ are horizontally and vertically polarized, respectively. Since the BS input modes $a$ and $b$ are related to the output modes $c$ and $d$ according to the following relation,

[TABLE]

the state of light at the output modes $c$ and $d$ of the second BS is given by,

[TABLE]

where $\{C.C.\}$ denotes the complex conjugate.

Let us first consider the second-order intensity cross correlation, $R_{cd}(\tau)$, between modes $c$ and $d$, defined as

[TABLE]

where $E^{(+)}_{c}(t)=\frac{1}{\sqrt{2\pi}}\int\hat{c}(\omega)e^{-i\omega t}d\omega$ is the field operator at detector $D_{c}$ and $\hat{c}(\omega)$ is the photon annihilation operator in mode $c$. As mentioned in section II, polarization projection precedes photon detection and to ensure maximum interference visibility, the polarizer angles are set at $|D\rangle_{c}$ or $|A\rangle_{c}$. The polarizer is accounted for in the above equation by defining $\hat{c}(\omega)\equiv\frac{1}{\sqrt{2}}(\hat{c}_{H}(\omega)+\hat{c}_{V}(\omega))$ for the polarizer angle setting at $|D\rangle_{c}$. For polarization projection at $|A\rangle_{c}$, $\hat{c}(\omega)\equiv\frac{1}{\sqrt{2}}(\hat{c}_{H}(\omega)-\hat{c}_{V}(\omega))$. $E^{(+)}_{d}(t)$ and $\hat{d}(\omega)$ are defined similarly. In our analysis, we consider the polarization projection bases $\{|D\rangle_{c}$, $|D\rangle_{d}\}$ and $\{|A\rangle_{c}$, $|D\rangle_{d}\}$ as they offer the theoretical maximum visibility of the second-order interference choi17 ; yongsu ; yongsu2 . Note that Eq. (8) can be measured by using the experimental setup shown in Fig. 1(a).

Let us now consider conditional photon statistics at mode $c$, heralded by photon detection event at mode $d$. The relevant experimental setup is shown in Fig. 1(b). To investigate the conditional second-order correlation $g_{C}^{(2)}(\tau,\tau_{c})$, an Hanbury-Brown–Twiss setup is placed in mode $c$ and is triggered by the photon counting signal from detector $D_{d}$. Assuming the single-photon detectors cannot resolve the number of photons, the second-order correlation at mode $c$ heralded by photon detection at mode $d$ can be written as,

[TABLE]

Here, $I_{d}^{(p)}$ is the heralding probability caused by a $p$-photon detection at detector $D_{d}$. Since the photon counting detectors are not photon number resolving, a heralding signal may be caused either by a single-photon or by multiple photons. The heralding probability is given as,

[TABLE]

The two-fold coincidence probability $I_{de}^{(p)}$ ($I_{df}^{(p)}$) describes the coincidence events between the detectors $D_{d}$ and $D_{e}$ ($D_{d}$ and $D_{f}$), assuming that $p$-photon detection event occurs at detector $D_{d}$. $I_{de}^{(p)}$ and $I_{df}^{(p)}$ are given by,

[TABLE]

Here, $E^{(+)}_{e}(t)=\frac{1}{\sqrt{2\pi}}\int\hat{e}(\omega)e^{-i\omega t}d\omega$ and $E^{(+)}_{f}(t)=\frac{1}{\sqrt{2\pi}}\int\hat{f}(\omega)e^{-i\omega t}d\omega$ are the field operators at detector $D_{e}$ and $D_{f}$, respectively, with the annihilation operators $\hat{e}(\omega)=\frac{1}{\sqrt{2}}\hat{c}(\omega)$ and $\hat{f}(\omega)=\frac{i}{\sqrt{2}}\hat{c}(\omega)$.

The conditional coincidence count rate between $D_{e}$ and $D_{f}$, heralded by the photon detection event at $D_{d}$, is given by the triple coincidence probability,

[TABLE]

IV Analysis

Even at the single-photon regime of $|\alpha|^{2}$, the fact that the laser follows the Poisson photon number statistics makes the analytical analysis difficult. Nevertheless, the conditional second-order correlation $g_{C}^{(2)}(\tau,\tau_{c})$ in Eq. (9) can be obtained analytically if we approximate the initial input state of Eq. (1) up to a finite photon number $p$. In the following analysis, we consider the analytic forms of $g_{C}^{(2)}(\tau,\tau_{c})$ truncated at the $p$-photon Fock state term. We consider the case of $\tau_{c}=0$ because this is where the quantum nature of the heralded photon states can be best tested. For instance, the conditional second order correlation function truncated at $p=3$ for an arbitrary $\alpha$, is given as,

[TABLE]

where $\sigma$ denotes the bandwidth of the Gaussian-shaped frequency spectra of $\mathcal{F}(\omega_{a})$ and $\mathcal{F}(\omega_{b})$ in Eq. (7). Here, the $\mp$ sign is related to the polarization projection measurement $\{|D\rangle_{c}$, $|D\rangle_{d}\}$ and $\{|A\rangle_{c}$, $|D\rangle_{d}\}$. Note that, while truncating at a low photon number $p$ makes the analytic forms simpler, it could significantly alter the conditional second order correlation $g_{C}^{(2)}(\tau,0)$. The analytic forms of $g_{C}^{(2)}(\tau,0)$ up to $p=10$ are given in Appendix.

Figure 2 shows the numerical simulation of the intensity cross-correlation $R_{cd}(\tau)$ and the conditional second-order correlation $g_{C}^{(2)}(\tau,0)$ for weak coherent states with $|\alpha|=0.1$ and $|\alpha|=1.2$. To account for sufficiently large Fock state contributions, truncation is made at $p=10$. Also, we assumed that the frequency spectra $\mathcal{F}(\omega_{a})$ and $\mathcal{F}(\omega_{b})$ in Eq. (7) to be Gaussian with the bandwidth of $15/2\pi$ MHz in full width at half maximum (FWHM), corresponding to the coherence time $\tau_{\rm coh}\approx 260$ ns. The blue solid lines and red dashed lines correspond to projection measurement bases $\{|D\rangle_{c},|D\rangle_{d}\}$, and $\{|A\rangle_{c},|D\rangle_{d}\}$, respectively. The $R_{cd}(\tau)$ plots in Fig. 2(a) and Fig. 2(b) show the typical Shih-Alley/Hong-Ou-Mandel like peaks and dips with the limited visibility of $V=1/2$ hom ; shih ; yongsu ; yongsu2 . The conditional second order correlation $g_{C}^{(2)}(\tau,0)$ is shown in Fig. 2(c) and Fig. 2(d). While it is evident that the conditional photon statistics can be manipulated, even becoming super-bunched, i.e., $g_{C}^{(2)}(\tau,0)>2$ as in Fig. 2(d), $g_{C}^{(2)}(\tau,0)$ is never below 1, meaning that the heralded photon states always remain classical.

In Fig. 3, we show the conditional second order correlation $g_{C}^{(2)}(\tau,\tau_{c})$ for $|\alpha|=0.1$ and $|\alpha|=1.2$ under different $p$-photon Fock state truncation. By looking at the asymptotic behaviors of the second order correlation function, we can figure out whether truncation at the particular $p$-photon Fock state can be justified. First, the $g_{C}^{(2)}(\tau,\tau_{c})$ values are calculated for the condition $\tau=0$ and $\tau_{c}=0$, see Fig. 3(a) and Fig. 3(b). As noticed in Fig. 2, different polarization projections $\{|D\rangle_{c}$, $|D\rangle_{d}\}$ and $\{|A\rangle_{c}$, $|D\rangle_{d}\}$ result in different conditional photon number statistics.

Now, when $\tau$ is larger than the coherence time of the light, $\tau>\tau_{\rm coh}\approx 260$ ns, the conditional second order correlations are the same regardless of the polarization projection choices. Figure 3(c) shows $g_{C}^{(2)}(\tau,0)$ at $\tau=500$ ns. The conditional second order correlation $g_{C}^{(2)}(\tau>\tau_{\rm coh},0)$ starts out showing photon anti-bunching when the Fock state contributions are truncated at a low photon number $p$. However, as more and more $p$-photon Fock state components are taken into consideration, it reaches the asymptotic value of 1.5, which corresponds to the case when there is no heralding signal pandey14 .

What we find in Fig. 3 is that even for weak coherent state at the single-photon regime, $|\alpha|=0.1$, asymptotic behaviors are observed at relatively large $p=4$. For $|\alpha|=1.2$, asymptotic behaviors are not reached until $p=9$, meaning that truncation below $p=9$ would result incorrect results. If Fock state truncation is made before reaching the asymptotic value corresponding to a particular $\alpha$, it looks as though conditional photon anti-bunching were possible silva15 ; silva16 ; amaral16 . Such conditional photon anti-bunching from classical light, however, is purely due to improper handling of Fock state truncation. It is necessary to properly account for even seemingly negligible higher-order Fock state components of the coherent state.

Figure 4 shows the conditional second order correlation $g^{(2)}_{C}(\tau,\tau_{c})$, truncated at the $p=10$ photon Fock state, as functions of both $\tau$ and $\tau_{c}$. The results with the $\{|D\rangle_{c}$, $|D\rangle_{d}\}$ projection are presented in Fig. 4(a) and (b), while Fig. 4(c) and (d) show the results with the $\{|A\rangle_{c}$, $|D\rangle_{d}\}$ projection. The black solid line in each plot corresponds to the case of $\tau_{c}=0$, which is depicted in Fig. 2. It is clear that photon antibunching cannot be achieved by heralding if the input light is classical.

We also note that, when $\tau>\tau_{\rm coh}$, Figs. 2, 3, 4 show that $g_{C}^{(2)}(\tau,0)=1.5$ which is larger than the value for the coherent state. This increased value of the second order correlation is due to the fact that the light in mode $c$ is the result of interference between two mutually incoherent laser beams pandey14 .

V Conclusion

Generation of nonclassical light states is important in quantum optics and inexpensive and simple methods of generating various nonclassical light states would significantly advance the state of the art in experimental quantum information research. In a recent attempt to tackle such a challenging problem, it has been reported that photon anti-bunching or conditional sub-Poissonian photon number statistics can be obtained via second-order interference of mutually incoherent weak lasers and heralding based on photon counting silva15 ; silva16 ; amaral16 . Here, we have carried out detailed theoretical and numerical analyses on the limit of manipulating conditional photon number statistics via interference of weak lasers and heralding based on photon counting. We find that conditional photon number statistics can become super-Poissonian in such a scheme. We demonstrate explicitly however that it cannot become sub-Poissonian, i.e., photon anti-bunching cannot be obtained in such a scheme. Theoretical and numerical analyses show that such incorrect results can be obtained if one does not properly account for seemingly negligible higher-order photon number expansions of the coherent state even at the single-photon regime.

Additionally, our work shows that a light beam having super-Poissonian photon number statistics and photon super-bunching may be easily generated by interfering multiple mutually incoherent laser beams. Such a simple scheme for generating super-bunched light may be of use in optical super-resolution and ghost imaging/interference experiments cao08 ; zhou10 .

Acknowledgement

This work was supported by the National Research Foundation (Grant No. 2016R1A2A1A05005202 and No. 2016R1A4A1008978), the ICT R$\&$D program of MSIP/IITP (B0101-16-1355), and the KIST institutional program (Project No. 2E26680-16-P024, 2E27231, 2V05340).

Appendix A Solutions for $g^{(2)}_{C}(\tau,\tau_{c})$

Here, we show the analytical forms of $g^{(2)}_{C}(\tau,0)$ in Eq. (9), truncated at the $p$-photon Fock state, for the projection basis $\{|D\rangle_{c},|D\rangle_{d}\}$. The cases up to $p=10$ are shown here.

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

[TABLE]

For the $\{|A\rangle_{c},|D\rangle_{d}\}$ projection, the results are similar to the ones shown above with some signs flipped so that the ‘dip’ becomes ’peak’.