Self-homodyne enabled generation of indistinguishable photons
The rapid generation of non-classical light serves as the foundation for exploring quantum optics and developing applications such as secure communication or generation of NOON-states. While strongly coupled quantum dot-photonic crystal resonator systems have great potential as non-classical light sources due to their promise of tailored output statistics, the generation of indistinguishable photons has been obscured due to the strongly dissipative nature of such systems. Here, we demonstrate that the recently discovered self-homodyne suppression technique can be used to overcome this limitation and tune the quantum statistics of transmitted light, achieving indistinguishable photon emission competitive with state-of-the-art metrics. Furthermore, our nanocavity-based platform directly lends itself to scalable on-chip architectures for quantum information.
These authors contributed equally.
©2016 Optical Society of America. One print or electronic copy may be made for personal use only. Systematic reproduction and distribution, duplication of any material in this paper for a fee or for commercial purposes, or modifications of the content of this paper are prohibited.
Understanding the interaction between light and matter is of paramount importance for exploring the peculiar properties of quantum optics and utilizing them for applications such as non-classical light generation with implications for communication, information processing and sensing Gisin et al. (2002); Gisin and Thew (2007); O’Brien (2007); Boto et al. (2000); Giovannetti et al. (2004). In the solid-state, self-assembled quantum dots (QDs) are widely used as quantum emitters due to their strong interaction with light and ability to be integrated into nanophotonic resonators for enhanced light-matter interaction. Examples of quantum optical landmark experiments with QDs are the generation of indistinguishable photons Santori et al. (2002); Gazzano et al. (2013); He et al. (2013); Wei et al. (2014); Gschrey et al. (2015); Thoma et al. (2016), entangled photon pairs Benson et al. (2000); Akopian et al. (2006); Trotta et al. (2014); Müller et al. (2014) and the observation of Mollow triplets Nick Vamivakas et al. (2009); Flagg et al. (2009).
For off-chip applications a high photon-extraction efficiency is desirable, which can be achieved by embedding the QD into a resonator with strong vertical emission. For example, QDs have been embedded in micropillar cavities Santori et al. (2002); Gazzano et al. (2013) that also Purcell enhance the emission rate for fast photon extraction. Utilizing resonant excitation of such structures, highly indistinguishable photon generation has recently been demonstrated Somaschi et al. (2016); Ding et al. (2016); Unsleber et al. (2016). Importantly, resonant excitation enables a high degree of indistinguishability due to the absence of excitation timing jitter and electric field noise that typically results from charge fluctuations in the semiconductor environment under non-resonant excitation He et al. (2013); Wei et al. (2014).
On the other hand, for on-chip applications photonic crystal resonators are promising. Their planar geometry naturally allows for coupling with on-chip structures, including waveguides and on-chip detectors Reithmaier et al. (2013, 2015) and hence holds promise to realize fully-integrated quantum optical hardware. Photonic crystal resonators provide extremely small mode volumes which enable a large enhancement of the light-matter interaction strength with embedded quantum emitters Yoshie et al. (2004); Englund et al. (2005); Hennessy et al. (2007). Importantly, in transmission geometries for on-chip resonant generation of non-classical light can not be achieved by resonantly exciting a QD weakly coupled to a resonator. Resonant excitation of a weakly coupled QD-resonator system predominantly leads to coherent scattering from the cavity. Thus, in off-chip applications with micropillars, suppressing this coherent scattering while collecting emission from the QD using cross-polarized suppression is only enabled by a careful choice of quantum emitter structure. Specifically, the structure must allow that the polarization of the coherently scattered light is not rotated while the polarization of the QD emission is rotated. This can be achieved, for example, in bi-modal cavities and charge neutral QDs that have their symmetry axis different from the cavity and laser Giesz et al. (2016) or bi-modal cavities and charged QDs. Instead in photonic crystals, direct transmission of light through a strongly coupled QD-nanocavity system can generate a range of output quantum statistics Faraon et al. (2008); Reinhard et al. (2011); Muñoz et al. (2014); Müller et al. (2015a, b). However, the strongly dissipative nature of such systems has so far obscured the generation of indistinguishable photons.
In this letter, we demonstrate that interference which is intrinsic to photonic crystal cavities can be used to overcome this strongly dissipative nature and tune quantum statistics. Specifically, we show that this recently-discovered self-homodyne suppression (SHS) effect Fischer et al. (2016) can be used to interferometrically reject the coherent scattering off a dissipative Jaynes-Cummings system and isolate the non-classical component of the emitted light. While our experimental approach is tailored to photonic crystal cavities, self-homodyne suppression as a tool for engineering quantum statistics is widely applicable to other Jaynes-Cummings systems. Here, we demonstrate the robust and ultrafast generation of highly-indistinguishable photons from strongly coupled quantum dot-photonic crystal resonator systems with state-of-the-art indistinguishability and generation rates. Additionally, this approach circumvents the temperature limit, set by phonon-dephasing, in all previous solid-state approaches Santori et al. (2002); Gschrey et al. (2015); Thoma et al. (2016) while also facilitating on-chip integration.
The system under investigation consists of a single self-assembled quantum dot (QD) strongly coupled to a photonic crystal L3 cavity (Fig. 1a). The resulting energy level structure is well described by the Jaynes-Cummings (JC) dressed states ladder. The energies of the lowest two rungs are presented in Fig. 1b as a function of the QD-cavity detuning . They form pairs of anticrossing lines, labelled UP and LP for the upper and lower polariton of rung , respectively. Experimentally this coupling can be observed in detuning-dependent cross-polarized reflectivity measurements (Fig. 1c) that reveal the clear anticrossing of the first rung Englund et al. (2007). As the polaritonic peaks transition through the avoided-crossing (at = 0), they exchange character from cavity/QD-like to QD/cavity-like. Fitting the data results in a QD-cavity coupling strength of and a cavity energy decay rate of .
Resonant generation of single photons in such systems can be achieved by photon blockade Faraon et al. (2008). Here, an excitation laser tuned in resonance with the first rung is out of resonance with higher climbs up the ladder due to the JC anharmonicity. However, these resonances have broad linewidths and hence appreciable overlap due to the highly dissipative character of semiconducting systems. Nevertheless, by detuning the QD and cavity by a few and exciting the QD-like polariton branch, high purity and efficiency single-photon generation has recently been demonstrated (blue arrow in Fig. 1b) Müller et al. (2015a). To achieve the strongest photon blockade in pulsed on-demand applications, the pulse length has to be chosen as a compromise that minimizes both re-excitation and resonance overlap Müller et al. (2015b).
We now discuss photon blockade in the context of self-homodyne suppression. Here, the suppression results from destructive interference between the light scattered from the fundamental cavity mode and the continuum above-the-light-line modes, leading to Fano-like resonances Vasco et al. (2013). This effect can be used to significantly suppress JC coherent scattering of the excitation laser in detuned QD-cavity systems and extract the incoherent spectrum Fischer et al. (2016). To reach SHS, we optimized the excitation conditions (focus and polarisation) for the suppression of coherent scattering. The pulse length is chosen to be only to minimize re-excitation. The resulting spectrum (Fig. 1d) exhibits three distinct features: emission from the resonantly excited QD-like UP1, phonon-assisted emission from the cavity-like LP1 (purple and red arrows in Fig. 1b) and self-homodyne suppressed JC coherently scattered laser light. The last one is strongest on the sides of the UP1 peak due to a spectral dependence of SHS that results from the wavelength-dependent phase shift of the JC coherently scattered light.
To investigate the single-photon generation and photon indistinguishability under resonant excitation of UP1 (Fig. 1d), we measure photon correlations between the outputs of a fiber-based Mach-Zehnder (MZ) interferometer (Fig. 2a). Here, we excite the system with double pulses that each have a pulse area of (see Supplementary information) and a time delay that matches the delay of the interferometer. First, we perform experiments without spectral filtering. The result (Fig. 2b) is a pattern of five peaks separated by and repeated with the repetition rate of the laser (). Due to the quantum character of the emission, the three centre peaks around zero time delay are attenuated Santori et al. (2002). Note that the asymmetry of the five peaks results from the imperfect reflectivity to transmittivity ratios of the second beamsplitter in the fiber-coupled implementation (see Supplementary information for details). To quantitatively analyze the data we bin the counts in a time window of about the peaks (Fig. 2b datapoints). A fit to the data (blue columns in Fig. 2c) allows for extraction of the measured degree of second-order coherence and first-order coherence between two subsequent pulses. The extracted values of the fit are and . In the literature, when analyzing the attenuation of the center peak instead of , a quantity is often stated and defined as the single-photon mode overlap. However, this parameter would only correspond to the single-photon mode overlap for pulses of perfect single-photon character (). The limited fidelity of the measurement can be understood by recalling the emission spectrum presented in Fig. 1d. The imperfect suppression of the JC coherently scattered light limits , while the phonon-assisted emission from LP1 limits due to excitation timing jitter.
To increase the fidelity of indistinguishable photon generation through photon blockade, we now employ spectral filtering. Therefore, we repeat the correlation measurement while filtering on the UP1 emission, as indicated by the red shaded region in Fig. 1d. The result of this experiment with a filter bandwidth of is presented in Fig. 2d as empty diamonds. A fit (red columns) extracts values of and . Note that this bandwidth is much larger than the linewidth of the UP1 emission. Therefore, the improvements result only from eliminating imperfect suppression of JC coherently scattered light and phonon-assisted emission from LP1 and not from filtering with a bandwidth smaller than spectral diffusion of the QD. In both cases, frequency filtered and unfiltered we confirmed the extracted values of in second-order correlation measurements and obtained similar values. Our metrics are competitive with the best values obtained from QDs so far He et al. (2013); Wei et al. (2014); Somaschi et al. (2016); Ding et al. (2016). The UP1 state lifetime in our experiment of (measured at this detuning Müller et al. (2015b)) paves the way for on-chip generation rates over an order of magnitude faster than bulk QDs and slightly faster than those in micropillar resonators Somaschi et al. (2016); Ding et al. (2016); Unsleber et al. (2016). However, micropillar resonators are optimized for photon extraction leading to higher count rates and they also do not require spectral filtering. Nevertheless, as discussed above the photonic crystal platform facilitates scalable on-chip architectures and in this platform, its near unity coupling efficiency to waveguides matters over its emission profile. Finally, the measurements presented here have been performed at a relatively high temperature of approximately 30 K, directly contrasting with the best previously reported HOM interference visibility of at such a temperature Thoma et al. (2016). This difference can be understood from interaction with the high-temperature phonon bath: In bulk, the interaction with phonons results in dephasing of the emission which reduces the first-order coherence. Meanwhile in a strongly coupled system, the interaction with phonons leads to a population transfer from UP1 to LP1 Müller et al. (2015b), spectrally removing the dephased emission from the detection channel and ensuring robust high-fidelity operation. Thus, we have investigated indistinguishable photon generation in a dephasing regime unlike all previous experiments and found that this region is highly beneficial for photon indistinguishability.
To corroborate our finding that the combination of SHS and spectral filtering results in high-fidelity generation of indistinguishable photons we performed quantum optical simulations (see Supplementary information). The simulated emission spectrum using the measured system parameters and excitation with a pulse is presented in Fig. 3a without (blue) and with optimized (black) SHS (see Supplementary information for details). Only when including SHS is the spectrum in extremely good qualitative agreement with the experimentally measured one. To demonstrate the impact of SHS on the single-photon character of the emission, we simulate as a function of the SHS strength without (Fig. 3b black) and with (Fig. 3b red) spectral filtering of the emission from UP1 . Here, the parameter denotes the intensity of the continuum-mode scattered contribution. In both cases, a clear dip in values of very close to the measured values is obtained, with much lower for the filtered system. Moreover, at the point of best suppression, the simulations reveal that the average number of photons exiting the system per pulse is unity.
To investigate photon indistinguishability, we calculate the second-order cross correlation discussed in the original Hong-Ou-Mandel (HOM) paper Hong et al. (1987). Here, correlations between the output ports of a beamsplitter are simulated while two identical systems feed the input ports. It is important to note that in contrast to the often-confused statements in the literature, in this configuration the dip in the correlation is different than the dip of the center peak of the Fig. 2a MZ implementation (see Supplementary information for details). This distinction is important because the MZ scheme is now predominantly used to experimentally characterize single-photon source indistinguishability in our and other experiments. Here, an important result of our work is that the HOM configuration shows a dip of
while the MZ zero delay dip reduces to:
Therefore, the two are equal only for and . To visualize this difference, we look at our experimentally obtained values of and discussed above. For the frequency filtered case we obtain and while for the unfiltered case we obtain and . In the filtered case, the difference is small due to the comparable values of and but in the unfiltered case the difference is significant. Nevertheless, when extracting and , it is possible to directly compare the values obtained from the two different methods, as well as simulations and experiment.
The results of the HOM simulations are presented in Figs 3c and 3d excluding and including SHS, respectively. The figures show against the filtering wavelength, as solid lines. The area below the line is decomposed into (shaded) and (white). Without SHS only a weak dip at the wavelength of LP1 is observed. However, it is barely nonclassical since at this wavelength the emission is phonon-mediated and, thus, subject to strong excitation timing jitter which is reflected in the low value of . In contrast, when including SHS, a strong dip at the filtered wavelength of UP1 is observed with values of and , in excellent agreement with the experimentally measured values.
In summary, we have demonstrated that the self-homodyne technique can be used to interferometrically tune the output quantum statistics of Jaynes-Cummings systems. In particular, we showed that strongly coupled QD-photonic crystal nanocavity systems are capable of robust and high-fidelity generation of indistinguishable photons even at elevated temperatures, by combining resonant excitation, self-homodyne suppression and spectral filtering. Having produced indistinguishable photons from a state with a lifetime of only , our results could pave the way for sources with unprecedented rates. Moreover, the short lifetime leads to a homogeneous linewidth of the emission which is much larger than that of bulk QDs. Therefore, we expect the indistinguishability to be unaffacted by spectral diffusion of the quantum emitter even without active suppression of the spectral diffusion Prechtel et al. (2013); Hansom et al. (2014); Kuhlmann et al. (2015). Specifically, in contrast to bulk QDs He et al. (2013); Thoma et al. (2016) we expect a similarly high indistinguishability for excitation pulses with a longer time delay or for measurements interfering the emission from multiple systems. Furthermore, the generation of indistinguishable photons from strongly coupled QD-photonic crystal systems enables scalable on-chip architectures. While in our case the positioning of the QD relative to the cavity has been done probabilistically, recent progress in site-selective growth of QDs Yakes et al. (2013); Unsleber et al. (2015) as well as positioning of resonators relative to QDs Somaschi et al. (2016); Gschrey et al. (2015); Unsleber et al. (2016) provides further support for the feasibility of integrated quantum photonic circuits. Finally, the demonstration of the self-homodyne technique to isolate the quantum character of resonantly scattered light paves the way for photon bundling Muñoz et al. (2014) at significantly lower nonlinearity-cavity detunings and smaller powers for high-throughput generation of other nonclassical states of light.
We gratefully acknowledge financial support from the Air Force Office of Scientific Research, MURI center for multifunctional light-matter interfaces based on atoms and solids (Grant No. FA9550-12-1-0025), support from the Army Research Office (Grant No. W911NF1310309) and support from the National Science Foundation (Division of Materials Research - Grant. No. 1503759). KM acknowledges support from the Alexander von Humboldt Foundation. KAF acknowledges support from the Lu Stanford Graduate Fellowship and the National Defense Science and Engineering Graduate Fellowship. CD acknowledges support from the Stanford Graduate Fellowship. KGL acknowledges support from the Swiss National Science Foundation. YAK acknowledges support from the Stanford Graduate Fellowship and the National Defense Science and Engineering Graduate Fellowship.
- Gisin et al. (2002) N. Gisin, G. Ribordy, W. Tittel, and H. Zbinden, Rev. Mod. Phys. 74, 145 (2002).
- Gisin and Thew (2007) N. Gisin and R. Thew, Nature Photonics 1, 165 (2007).
- O’Brien (2007) J. L. O’Brien, Science 318, 1567 (2007), http://science.sciencemag.org/content/318/5856/1567.full.pdf .
- Boto et al. (2000) A. N. Boto, P. Kok, D. S. Abrams, S. L. Braunstein, C. P. Williams, and J. P. Dowling, Phys. Rev. Lett. 85, 2733 (2000).
- Giovannetti et al. (2004) V. Giovannetti, S. Lloyd, and L. Maccone, Science 306, 1330 (2004), http://www.sciencemag.org/content/306/5700/1330.full.pdf .
- Santori et al. (2002) C. Santori, D. Fattal, J. Vučković, G. S. Solomon, and Y. Yamamoto, Nature 419, 594 (2002).
- Gazzano et al. (2013) O. Gazzano, S. Michaelis de Vasconcellos, C. Arnold, A. Nowak, E. Galopin, I. Sagnes, L. Lanco, A. Lemaître, and P. Senellart, Nature Communications , 1425– (2013).
- He et al. (2013) Y.-M. He, Y. He, Y.-J. Wei, D. Wu, M. Atature, C. Schneider, S. Hofling, M. Kamp, C.-Y. Lu, and J.-W. Pan, Nat Nano 8, 213 (2013).
- Wei et al. (2014) Y.-J. Wei, Y.-M. He, M.-C. Chen, Y.-N. Hu, Y. He, D. Wu, C. Schneider, M. Kamp, S. Höfling, C.-Y. Lu, and J.-W. Pan, Nano Letters 14, 6515 (2014), pMID: 25357153, http://dx.doi.org/10.1021/nl503081n .
- Gschrey et al. (2015) M. Gschrey, A. Thoma, P. Schnauber, M. Seifried, R. Schmidt, B. Wohlfeil, L. Kruger, J. H. Schulze, T. Heindel, S. Burger, F. Schmidt, A. Strittmatter, S. Rodt, and S. Reitzenstein, Nat Commun 6 (2015).
- Thoma et al. (2016) A. Thoma, P. Schnauber, M. Gschrey, M. Seifried, J. Wolters, J.-H. Schulze, A. Strittmatter, S. Rodt, A. Carmele, A. Knorr, T. Heindel, and S. Reitzenstein, Phys. Rev. Lett. 116, 033601 (2016).
- Benson et al. (2000) O. Benson, C. Santori, M. Pelton, and Y. Yamamoto, Phys. Rev. Lett. 84, 2513 (2000).
- Akopian et al. (2006) N. Akopian, N. H. Lindner, E. Poem, Y. Berlatzky, J. Avron, D. Gershoni, B. D. Gerardot, and P. M. Petroff, Phys. Rev. Lett. 96, 130501 (2006).
- Trotta et al. (2014) R. Trotta, J. S. Wildmann, E. Zallo, O. G. Schmidt, and A. Rastelli, Nano Letters 14, 3439 (2014), pMID: 24845369, http://dx.doi.org/10.1021/nl500968k .
- Müller et al. (2014) M. Müller, S. Bounouar, K. D. Jöns, M. Glässl, and P. Michler, Nature Photonics 8, 224 (2014).
- Nick Vamivakas et al. (2009) A. Nick Vamivakas, Y. Zhao, C.-Y. Lu, and M. Atature, Nat Phys 5, 198 (2009).
- Flagg et al. (2009) E. B. Flagg, A. Muller, J. W. Robertson, S. Founta, D. G. Deppe, M. Xiao, W. Ma, G. J. Salamo, and C. K. Shih, Nat Phys 5, 203 (2009).
- Somaschi et al. (2016) N. Somaschi, V. Giesz, L. De Santis, J. C. Loredo, M. P. Almeida, G. Hornecker, S. L. Portalupi, T. Grange, C. Antón, J. Demory, C. Gómez, I. Sagnes, N. D. Lanzillotti-Kimura, A. Lemaítre, A. Auffeves, A. G. White, L. Lanco, and P. Senellart, Nat Photon 10, 340 (2016).
- Ding et al. (2016) X. Ding, Y. He, Z.-C. Duan, N. Gregersen, M.-C. Chen, S. Unsleber, S. Maier, C. Schneider, M. Kamp, S. Höfling, C.-Y. Lu, and J.-W. Pan, Phys. Rev. Lett. 116, 020401 (2016).
- Unsleber et al. (2016) S. Unsleber, Y.-M. He, S. Gerhardt, S. Maier, C.-Y. Lu, J.-W. Pan, N. Gregersen, M. Kamp, C. Schneider, and S. Höfling, Opt. Express 24, 8539 (2016).
- Reithmaier et al. (2013) G. Reithmaier, S. Lichtmannecker, T. Reichert, P. Hasch, K. Muller, M. Bichler, R. Gross, and J. J. Finley, Sci. Rep. 3 (2013), 10.1038/srep01901.
- Reithmaier et al. (2015) G. Reithmaier, M. Kaniber, F. Flassig, S. Lichtmannecker, K. Müller, A. Andrejew, J. Vučković, R. Gross, and J. J. Finley, Nano Letters 15, 5208 (2015), pMID: 26102603, http://dx.doi.org/10.1021/acs.nanolett.5b01444 .
- Yoshie et al. (2004) T. Yoshie, A. Scherer, J. Hendrickson, G. Khitrova, H. M. Gibbs, G. Rupper, C. Ell, O. B. Shchekin, and D. G. Deppe, Nature 432, 200 (2004).
- Englund et al. (2005) D. Englund, D. Fattal, E. Waks, G. Solomon, B. Zhang, T. Nakaoka, Y. Arakawa, Y. Yamamoto, and J. Vučković, Phys. Rev. Lett. 95, 013904 (2005).
- Hennessy et al. (2007) K. Hennessy, A. Badolato, M. Winger, D. Gerace, M. Atatüre, S. Gulde, S. Fält, E. L. Hu, and A. Imamoğlu, Nature 445, 896 (2007).
- Giesz et al. (2016) V. Giesz, N. Somaschi, G. Hornecker, T. Grange, B. Reznychenko, L. De Santis, J. Demory, C. Gomez, I. Sagnes, A. Lemaitre, O. Krebs, N. D. Lanzillotti-Kimura, L. Lanco, A. Auffeves, and P. Senellart, Nat Commun 7 (2016).
- Faraon et al. (2008) A. Faraon, I. Fushman, D. Englund, N. Stoltz, P. Petroff, and J. Vuckovic, Nature Physics 4, 859 (2008).
- Reinhard et al. (2011) A. Reinhard, T. Volz, M. Winger, A. Badolato, K. J. Hennessy, E. L. Hu, and A. Imamoglu, Nature Photonics 6, 93–96 (2011).
- Muñoz et al. (2014) C. S. Muñoz, E. del Valle, A. G. Tudela, K. Müller, S. Lichtmannecker, M. Kaniber, C. Tejedor, J. J. Finley, and F. P. Laussy, Nature Photonics 8, 550 (2014).
- Müller et al. (2015a) K. Müller, A. Rundquist, K. A. Fischer, T. Sarmiento, K. G. Lagoudakis, Y. A. Kelaita, C. Sánchez Muñoz, E. del Valle, F. P. Laussy, and J. Vučković, Phys. Rev. Lett. 114, 233601 (2015a).
- Müller et al. (2015b) K. Müller, K. A. Fischer, A. Rundquist, C. Dory, K. G. Lagoudakis, T. Sarmiento, Y. A. Kelaita, V. Borish, and J. Vučković, Phys. Rev. X 5, 031006 (2015b).
- Fischer et al. (2016) K. A. Fischer, K. Müller, A. Rundquist, T. Sarmiento, A. Y. Piggott, Y. Kelaita, C. Dory, K. G. Lagoudakis, and J. Vučković, Nat Photon 10, 163 (2016).
- Englund et al. (2007) D. Englund, A. Faraon, I. Fushman, N. Stoltz, P. Petroff, and J. Vuckovic, Nature 450, 857 (2007).
- Vasco et al. (2013) J. P. Vasco, H. Vinck-Posada, P. T. Valentim, and P. S. S. Guimãraes, Opt. Express 21, 31336 (2013).
- Hong et al. (1987) C. K. Hong, Z. Y. Ou, and L. Mandel, Phys. Rev. Lett. 59, 2044 (1987).
- Prechtel et al. (2013) J. H. Prechtel, A. V. Kuhlmann, J. Houel, L. Greuter, A. Ludwig, D. Reuter, A. D. Wieck, and R. J. Warburton, Phys. Rev. X 3, 041006 (2013).
- Hansom et al. (2014) J. Hansom, C. H. H. Schulte, C. Matthiesen, M. J. Stanley, and M. Atatüre, Applied Physics Letters 105, 172107 (2014), http://dx.doi.org/10.1063/1.4901045.
- Kuhlmann et al. (2015) A. V. Kuhlmann, J. H. Prechtel, J. Houel, A. Ludwig, D. Reuter, A. D. Wieck, and R. J. Warburton, Nat Commun 6 (2015).
- Yakes et al. (2013) M. K. Yakes, L. Yang, A. S. Bracker, T. M. Sweeney, P. G. Brereton, M. Kim, C. S. Kim, P. M. Vora, D. Park, S. G. Carter, and D. Gammon, Nano Letters 13, 4870 (2013), pMID: 23987910, http://dx.doi.org/10.1021/nl402744s .
- Unsleber et al. (2015) S. Unsleber, S. Maier, D. P. S. McCutcheon, Y.-M. He, M. Dambach, M. Gschrey, N. Gregersen, J. Mørk, S. Reitzenstein, S. Höfling, C. Schneider, and M. Kamp, Optica 2, 1072 (2015).
- Akahane et al. (2005) Y. Akahane, T. Asano, B.-S. Song, and S. Noda, Opt. Express 13, 1202 (2005).
- Minkov and Savona (2014) M. Minkov and V. Savona, Scientific Reports 4 (2014), 10.1038/srep05124.
- Johansson et al. (2013) J. Johansson, P. Nation, and F. Nori, Computer Physics Communications 184, 1234 (2013).
- Alsing et al. (1992) P. Alsing, D.-S. Guo, and H. J. Carmichael, Phys. Rev. A 45, 5135 (1992).
- del Valle et al. (2012) E. del Valle, A. Gonzalez-Tudela, F. P. Laussy, C. Tejedor, and M. J. Hartmann, Phys. Rev. Lett. 109, 183601 (2012).
- Woolley et al. (2013) M. J. Woolley, C. Lang, C. Eichler, A. Wallraff, and A. Blais, New J. Phys. 15, 105025 (2013).
Sample fabrication: The MBE-grown structure consists of an thick sacrificial layer followed by a thick GaAs layer that contains a single layer of InAs QDs. Our growth conditions result in a typical QD density of . Using e-beam lithography with ZEP resist, followed by reactive ion etching and HF removal of the sacrificial layer, we define the photonic crystal cavity. The photonic crystal lattice constant was and the hole radius . The cavity fabricated is a linear three-hole defect (L3) cavity with a measured quality factor of . To improve the cavity quality factor, holes adjacent to the cavity were shifted Akahane et al. (2005); Minkov and Savona (2014).
Optical spectroscopy: All optical measurements were performed with a liquid helium flow cryostat at temperatures in the range . For excitation and detection a microscope objective with a numeric aperture of was used. Cross-polarised measurements were performed using a polarising beam splitter. To further enhance the extinction ratio, additional thin film linear polarisers were placed in the excitation/detection pathways and a single mode fibre was used to spatially filter the detection signal. Furthermore, two waveplates were placed between the beamsplitter and microscope objective: a half-wave plate to rotate the polarisation relative to the cavity and a quarter-wave plate to correct for birefringence of the optics and sample itself. We estimate our overall efficiency to be , including setup losses and detector efficiency. Excluding the detector efficiency this corresponds to a collection efficiency into the detection single mode fiber of
Correlation measurements: Second-order autocorrelation measurements were performed using a Hanbury Brown and Twiss (HBT) setup consisting of one fibre beamsplitter and two single-photon avalanche diodes. The detected photons were correlated with a PicoHarp300 time counting module. Two-photon interference measurements were performed with the same detectors and electronics as the second-order correlation measurements but with a Mach-Zehnder interferometer replacing the single beamsplitter, as discussed in the main text. For the unfiltered two-photon interference measurements the integration time was 30 minutes and the coincidence rate (in the whole histogram) was 2336/s. For the frequency filtered two-photon interference measurements the integration time was 8.6 hours and the coincidence rate (in the whole histogram) was 22/s.
Details on the simulations
The quantum-optical simulations were performed using density matrix master equations with the Quantum Optics Toolbox in Python (QuTiP)Johansson et al. (2013), where the standardAlsing et al. (1992) Jaynes-Cummings model was used as a starting point. The effects of phonons were incorporated through the addition of incoherent decay channels with rates that were previously extractedMüller et al. (2015b).
To simulate the first order spectra of our system under excitation with a single pulse, we compute the one-sided spectrum
of the free-field mode operator . Input-output theory can relate the internal cavity mode operator to the external field operator by the radiative cavity field decay rate . Hence, for a JC system in the solid state where the QD radiative decay rate plays an insignificant role compared to Müller et al. (2015b), spectral decomposition of the cavity mode operator yields the spectrum of the detected light. Therefore, we can compute an unnormalised version of this spectrum with in equation (1). We can also compute an unnormalised version of the incoherent spectrum with in equation (1). To arrive at the version measured by a spectrometer of finite bandwidth, we convolve with the spectrometer’s response function. To simulate self-homodyne suppression (SHS), we replace in equation (1). Physically, is a slightly phase- and amplitude- shifted version of the incident laser pulse (originating from the continuum-mode scattering)Fischer et al. (2016).
In order to simulate the normalised measured degree of second-order coherence, with , we calculate
under excitation by a single pulseMüller et al. (2015a). The operators indicate the time ordering required of a physical measurementdel Valle et al. (2012). We can likewise replace in equation (2) and also model SHS with the replacement of in equation (2). Despite the simplicity of equation (1), adding spectral filtering to equation (2) is analytically and numerically quite challenging. The spectral decomposition of this equation requires a fourth order integral that is often intractable even numerically. Fortunately, the newly discovered sensor formalismdel Valle et al. (2012) allows for efficient calculation of the spectrally filtered version of the measured degree of second-order coherence. Here, we coherently attach a pair of two-level sensors to the system Hamiltonian with the addition of the sensor Hamiltonian to the Jaynes-Cummings Hamiltonian:
where is the sensor frequency, the sensor annihilation operator, and the sensor coherent coupling strength. The sensor coupling is chosen small enough so that its backaction on the system is negligible, i.e. where is the fastest transition rate in the un-sensed system. Additionally, the sensor decay terms of rate are added to the total Liouvillian. Here, in order to simulate SHS, we replace
in equation (3).
To arrive at the physically measured and spectrally filtered second-order coherence functions, the total degree of second-order coherence is computed between the two sensors:
As the sensors are degenerate in every manner, the ordering of their operation is arbitrary. In our model, the sensors are used as filters while the detector is assumed to be sufficiently broad-band to integrate the correlations over our entire experimental domain. This approximation is accurate as the detector has a timing resolution of greater than compared with the system decay time of approximately . Additionally, we integrate our data over bins.
Simulated power-dependent spectra
We began our simulations by computing the incoherent power spectrum to find the pulse driving power. As all experiments were performed with a pulse driving, we wanted to ensure an accurate match with simulations. Since the Rabi oscillations are clearly present in the incoherent spectrum without having to optimise SHS, this served as a simple way to verify the pulse power. The result of such a simulation is presented in Fig. 4 for a resonant excitation of UP1 with short pulses at a QD-cavity detuning of . For this detuning, the forward (UP1 to LP1) and reverse (LP1 to UP1) phonon-assisted transfer rates are and respectively.
Note that after optimizing SHS, a nearly identical simulation of the measured spectrum is observed (see main text). For both UP1 and LP1 Rabi oscillations are observed and the pulse area for a pulse is close to where is the width parameter of the Gaussian pulse. This analytic result comes from truncating the detuned Jaynes-Cummings level structure and calculating the Raman transition frequency between the lowest-lying levelsFischer et al. (2016). From there, the pulse driving strength is given approximately by integrating the pulse area - the exact drive area is difficult to analytically calculate due to the nonlinear damping of the system.
We first review interference in a Hong-Ou-Mandel (HOM) interferometer. In such an interferometer, two identically independent sources are interfered on two detectors by the action of a single beamsplitter. This setup is shown schematically in Fig. 5. Here, we consider two strongly coupled systems each driven with the same pulse . Their outputs, the Heisenberg free-field operators and , are then fed into the beamsplitter which mixes the two according to the unitary transformation
The detection events are then correlated electronically to arrive at a temporal coincidence histogram . In terms of the instantaneous correlations aforementioned, we wish to calculate . In order to perform this calculation, we first decompose the underlying correlations of into more manageable components based on and that will allow us to average over all quickly varying phase terms. The phase terms are difficult to observe experimentally since they depend on femtosecond phase locking and are, regardless, not required for the observation of two-photon interference. It has been shown that after performing this procedureWoolley et al. (2013) one arrives at
where the correlations between and have terms only coming from a single source (the indices could be trivially switched as and are identically independent). Experimentally due to unknown detection efficiencies, our only true intensity reference is in proportion to the source correlation . Combining this information with the discussion on detector bandwidth above, we are justified to compute only the integrated and normalised correlation
referenced in the main text (equation ). Hence the observable HOM interference for indentically independent sources simply consists of an instantaneous measurement that depends on the total first- and second-order coherence of the source. Whether this interference is truly just two-photon interference therefore depends on the nature of the source output. Additionally, we note that in many real systems the first-order coherence extracted here is wildly different than the first-order coherence extracted from a Michelson interferometer due to the long-time averaging action of the Michelson interferometer.
Experimentally, producing a single source of indistinguishable photons is often quite challenging (let alone two), however we would still like to characterise the source’s instantaneous degree of first order coherence. By doing so, we hope to quantify how it would behave in a HOM experiment. It was realised that some aspect of two-photon interference can be observed for a single source at the input to a Mach-Zehnder (MZ) interferometerSantori et al. (2002). In order to observe this interference, the system must be doubly excited at a time interval that matches the temporal delay between the two paths of the MZ interferometer (shown schematically in Fig. 6). Additionally, this interval must be long enough that the emission resulting from the two excitations is identically independent.
For perfect single photon input the MZ output correlation is in fact two-photon interference and thus equivalent to HOM interference. However, in may cases the literature has incorrectly compared the two when the source has some finite probability of emitting non-single photons. Thus, we now outline the derivation of the correct measured correlation, as discussed in the main text (equation ). Consider the schematic in Fig. 6: A single strongly coupled system is driven with the coherent driving field and has field-mode operator output . The first beamsplitter simply splits in two. Here we have dropped the vacuum operator in anticipation that our detectors only measure a normally ordered moment of their input fields. Next, one path is delayed by the excitation interval . The second beamsplitter then mixes with the time delayed version of itself and the detection apparatus computes like before.
The time delay is absolutely critical: the purpose is that for time delays much less than , can interfere with itself, and therefore reproduce a correlation similar to the HOM cross-correlation. Although computation of in terms of source correlations certainly appears daunting, the independence of for times greater than dramatically simplifies the calculation of . This expansion is further simplified by consideration of the correlations centred around time delays that are integer multiples of . It is fairly trivial to show that the majority of terms for delays larger than are zero, with only a few additional phase-dependent terms arising. We note that although one might expect the correlations to be phase locked by the MZ interferometer, our integration times of many hours quickly destroy the phase interference. Since these terms once again play no role in two-photon interference, we made no effort to further stabilise our setup to observe them.
Consider the correlations about zero delay: and its time delayed version are statistically independent so equation (7) holds but with the source correlations of instead of . Because the source is doubly excited the ratio of the correlations in terms of is altered, which can easily be seen by applying the above rules. Now, we arrive at the measured degree of correlation
choosing instead to normalise to the maximum value of this expression we arrive at the equation discussed in the main text
Next, we discuss the influence of non-equal splitting ratios at the beamsplitters to arrive at the correlations comprising each of the five observed peaks. As presented in the main text, we built our Mach-Zehnder (MZ) interferometer from fiber-coupled beam splitters and polarisation maintaining fibers. While this realization of a MZ interferometer is very stable, the commercially available fiber-coupled beam splitters have an imperfect transmittivity () to reflectivity () ratio that deviates from . Therefore, we need to derive the amplitudes of the five peak pattern observed in the correlation measurements with the individual transmitivities and reflectivities of the first beamsplitter () and second beamsplitter () included. Both of our beamsplitters have a ratio of . The transmitted component at the first beamsplitter is sent to the longer arm of the interferometer. Importantly, since this delay is realised with an additional fiber there is an additional coupling loss resulting in an effective splitting ratio of the first beamsplitter of . Then, the amplitudes of the five peaks in the correlation measurements with SPCM2 as start and SPCM1 as stop are given by:
where A is the absolute amplitude. These equations can now be used to fit the measured data. We determine from fitting a number of peaks away from zero time delay using and , followed by fitting the five peaks around zero time delay. Note that when binning the data also an extremely weak dark-count noise was subtracted.
A binned version of the frequency filtered two-photon interference measurements was presented in Fig. 2d of the main text. For completeness, the raw data of this measurement is presented in Fig. 7.
Lower limit of
Finally, we discuss the non-zero value of . For strongly coupled QD-cavity systems strongly depends on the pulse length, which has to be chosen as a compromise between re-excitation during the presence of the pulse and spectral overlap of the laser with higher climbs up the JC ladder Müller et al. (2015b). However, when using SHS and frequency filtering the situation simplifies and becomes identical to that of a two-level system, where the only limitation is re-excitation which is determined by the pulse length. To visualize this dependence, we present in Fig. 8 a simulation of obtained from a two-level system under resonant excitation with pulses as a function of the pulse length. Clearly, increases with pulse length and very small values are only reached for pulses that are much shorter than the state lifetime.