## Abstract

Hybrid quantum information processing combines the advantages of discrete and continues variable protocols by realizing protocols consisting of photon counting and homodyne measurements. However, the mode structure of pulsed sources and the properties of the detection schemes often require the use of optical filters in order to combine both detection methods in a common experiment. This limits the efficiency and the overall achievable squeezing of the experiment. In our work, we use photon subtraction to implement the distillation of pulsed squeezed states originating from a genuinely spatially and temporally single-mode parametric down-conversion source in non-linear waveguides. Due to the distillation, we witness an improvement of 0.17 dB from an initial squeezing value of −1.648 ± 0.002 dB, while achieving a purity of 0.58, and confirm the non-Gaussianity of the distilled state via the higher-order cumulants. With this, we demonstrate the source’s suitability for scalable hybrid quantum network applications with pulsed quantum light.

© 2020 Optical Society of America under the terms of the OSA Open Access Publishing Agreement

## 1. Introduction

For practical quantum information processing (QIP), e.g. long-distance quantum communication, the ability to distribute quantum states is essential. However, this ability is limited by decoherence processes during transmission and the inability to deterministically amplify quantum states. The distillation of quantum resources, such as squeezing or entanglement, provides the means to undo parts of the detrimental effects of transmission [1].

In continuous-variable quantum information processing (CV-QIP), the field quadratures of Gaussian states are used to carry the information. Gaussian resource states such as squeezed or two-mode squeezed (entangled) states are routinely generated in laboratories e.g. through parametric down-conversion (PDC) in $\chi ^{(2)}$ crystals; either embedded in an optical cavity as an optical parametric oscillator below threshold [2–5] or in a high efficiency single-pass configuration [6–8]. Pulsed sources offer distinct benefits in terms of energy efficiency and source synchronization. However, they typically emit an ensemble of spatial and temporal modes. The disadvantage of this intrinsic multi-mode nature is that the squeezing is spread over a great number of modes, limiting the overall achievable squeezing in a single one of them, and that the separate detection of different modes requires mode-resolving detection schemes, e.g homodyne detection. In homodyne detection, a signal beam is interfered with a bright reference beam, the local oscillator (LO), and then usually detected by a pair of PIN photo diodes. The detection signal of interest is then retrieved when looking at the difference current of both diodes. Due to this interference, only those components of the signal will be detected which have a non-zero overlap with the LO field.

In addition to the properties of the source that generates the Gaussian states, we have to consider that the distillation of Gaussian states is a challenging task, which requires non-Gaussian (nG) operations [9]. This is particularly true for pulsed quantum light, which typically exhibits a rich temporal mode structure [10,11]. Among various methods of implementing nG operations [12–15], hybrid discrete- and continuous-variable operations are especially interesting because they combine the distinct features of homodyne and photon counting measurements [16–24]. In particular, photon subtraction has been used in many QIP scenarios. The photon subtraction is implemented by detecting a single photon from a pulse after a weakly reflecting beam splitter. Conventional single-photon detectors, in contrast to homodyne detectors, cannot resolve the temporal mode structure of the impinging ultrafast pulsed light field. Due to this mode-insensitivity, the photons are subtracted from an incoherent mixture of all excited modes degrading the purity of the final state [25]; posing this scheme unsuitable for multi-mode states. Usually, this issue is dealt with by employing a set of spatial and spectral filters [22–24], which in return reduces the success rate of the experiment. As an alternative solution, photon counting detectors can achieve mode-resolution when combined with an additional non-linear process that is mode-selective. Such experiments have been performed but come at the cost of an additional experimental overhead [11,26,27].

In our work, we present the successful distillation of pulsed squeezed states using a periodically poled potassium titanyl phosphate (PPKTP) waveguide as engineered source of broadband squeezed states. By using a single-mode waveguiding structure and engineering the non-linear interaction inside the crystal with ultrafast pump pulses, we are able to generate spatially single-mode and spectrally separable pulsed squeezed states of light thus reducing their effective number of modes. This reduced mode-number enables us to demonstrate the efficient combination of photon counting and homodyne measurements in a single protocol while, at the same time, avoiding the use of strict spectral filters. We witness an increase in the detected squeezing due to the distillation and can furthermore verify the non-Gaussian nature of the distilled pulsed quantum state, which both match with the predictions from our theoretical model. We want to emphasize, that the preparation of pulsed quantum states is of particular importance for synchronizing multiple conditional state preparations in network applications. Our results highlight the advantages of pulsed-state source engineering, which consist of a reduction in the experimental complexity and an increase in the efficiency of the state generation, in the context of hybrid QIP.

## 2. Theoretical background

In hybrid QIP experiments, we aim to combine the features of DV and CV measurements in order to realize operations inaccessible by DV or CV methods alone. A standard hybrid QIP protocol is the distillation of two-mode squeezing via the subtraction of photons from both modes of a two-mode squeezed state. In the following, we model a pulsed two-mode squeezed state generated by our source in an idealized setting. For this, we start by considering a PDC source generating pulsed two-mode squeezed states in a single spatio-temporal mode

The distillation is then realized by splitting the two-mode squeezed state on beam splitters with identical transmissivity $T$. The reflected portion is sent to click detectors implementing the photon subtraction while the transmitted part is forwarded to the next QIP step. We want to emphasize here that the click detectors are fully, i.e. temporally and spatially, mode-resolving in the sense that they always detect the same spatio-temporal mode that the following processing step operates on. The state after the photon subtraction can then analytically be expressed as [28]

In reality however, pulsed PDC processes emit temporal multi-mode two-mode squeezed states. The broadband generation process is usually described in the Heisenberg picture as

Yet in the actual experiment, the generated quantum states still exhibit small contributions from higher-order Schmidt modes, which need to be accounted for in the squeezing distillation process, if not dealt with by using optical filters. This is due to the fact that typical click detectors behave not ideally when detecting a multi-mode light field. In our squeezing distillation experiment as shown in Fig. 1, we need to consider the possibility that the avalanche photo diodes (APDs) randomly detect photons from different Schmidt modes. In our case however, we can restrict ourselves to the first two orders of Schmidt modes as the influence of modes of order higher than the first two is vanishing.

When one of the click detectors detects a photon, it will randomly fire from a photon from one of the Schmidt modes. We are triggering our homodyne measurement on the registration of a coincidence event, consequently several different realizations of that coincidence lead to a homodyne measurement in the multi-mode case. A coincidence might stem from two photons of the Schmidt mode coinciding with the LO, a photon in the LO mode and any other mode, and two photons from any other mode that we cannot distinguish between. Therefore, we model the detected state as a mixed state consisting of the superposition of the intended distilled state $\hat {\rho }_{\mathrm {sub(2)}}$, a state with a photon subtracted in one of the two modes $\hat {\rho }_{\mathrm {sub(1)}}$, and an undistilled squeezed state $\hat {\rho }_{\mathrm {sq}}$, each weighted with a parameter $\alpha _i$

## 3. Experimental results

#### 3.1 Experimental setup

Our experimental setup is described in Fig. 1, a mode-locked femto-second laser at a repetition rate of 82 MHz is pumping a synchronously pumped optical parametric oscillator with its central wavelength at 1540 nm which is used to drive our experiment. The light is frequency doubled via a periodically poled lithium niobate crystal, generating the pump for the PDC process. The up-converted beam is then filtered via a 4f-spectrometer in order to ascertain the spectral decorrelation of the pulsed squeezed states via an optimized pump spectral width [30,35]. In KTP, a type-II phasematched PDC process can be used to meet this decorrelation condition that allows us to generate genuinely single-mode states. In order to achieve single TM emission, we choose the pump pulse bandwidth to be 2.3 nm and its central wavelength to 769.6 nm so that we achieve the generation of frequency degenerate signal and idler beams around 1540 nm.

After the waveguide, we compensate for the birefringence of PPKTP using a Michelson-like interferometer with adequately chosen arm length so that signal and idler pulses have a high temporal overlap. We realize the photon subtraction by tapping off a small portion of the two-mode squeezed pulses via a 90/10 beam splitter which allows to compromise between a high distillation gain and the data acquisition rate. After the beam splitter, we filter signal and idler pulses using a bandpass with 10 nm FWHM to suppress any background modes. This filter bandwidth is chosen such that it does not influence the spectra of signal and idler that show spectral bandwidths $\Delta \lambda _s = 4.1~\mathrm {nm}$ and $\Delta \lambda _i = 5.8~\mathrm {nm}$ respectively. Thereafter, the modes are split using a polarizing beam splitter and subsequently sent onto a pair of APDs. Coincidence signals from the APDs will then be forwarded to herald a successful photon subtraction in the signal and idler modes. Following the brighter port of the asymmetric beam splitter, signal and idler pulses are interfered with each other on the combination of a half-wave plate and another polarizing beam splitter to generate a pair of single-mode squeezed beams of which one is sent onto a homodyne receiver with adapted LO pulses in order to record the quadrature data. The mixing of both modes of a two-mode squeezed state here can be understood as the inversion of generating two-mode squeezing by interfering two single-mode squeezed beam on a symmetric beam splitter [36].

We record the measurement data with a raw distillation rate is about $10~\mathrm {s}^{-1}$, due to technical limitations of the recording hardware however, we are limited to an effective distillation rate of $250~\mathrm {minute}^{-1}$. In total, we record 250000 time traces with a span of 100 $\mathrm {\mu s}$ and a sampling rate of $10^9~\mathrm {Samples}~\mathrm {s}^{-1}$. Each trace contains a single photon subtraction event surrounded by about 8000 squeezed states not subject to the photon subtraction. We assume, that the relative phase between squeezed light and LO stays constant during this time interval. This allows us to reconstruct the relative phase for every trace by using the variance of the squeezed states as a reference and, furthermore, allows us to directly compare the improvement in the quadrature variance due to the squeezing distillation. Further details on the phase estimation routine can be found in appendix A.

#### 3.2 Characterization

We characterize the properties of our source by a number of DV measurement techniques using APDs. A Hong-Ou-Mandel interference measurement yields us a visibility of $0.75 \pm 0.01$ when interfering the signal and idler beams. The mean photon number $\langle n \rangle$ is estimated to be $0.56 \pm 0.001$ photons per pulse which is equivalent to a two-mode squeezing of 6 dB at a pump pulse energy of $E_p = 41~\mathrm {pJ}$. The heralding efficiency after the asymmetric beam splitter is $\eta = 0.002$. Finally, we estimate the number of effective time-frequency modes of our squeezed state via a marginal $g^{(2)}$ measurement [29], resulting in a value of $g^{(2)}= 1.81 \pm 0.05$ and thus corresponding to an effective mode-number of 1.23.

By seeding the PDC process with coherent pulses at 1540 nm and measuring the optical power at the detection stage, we estimate the efficiency of the optical elements of the squeezing path to be 0.87. The optimization of the homodyne receiver, mainly the mode overlap between the pulsed LO and the squeezed state, is achieved by using the visibility between the seeded PDC and an attenuated LO as figure of merit. We optimize the visibility by shaping the spectral amplitude and phase of the LO pulses utilizing a commercially available waveshaper (Finisar waveshaper 4000s) and achieve a visibility of up to 0.91 between the seeded PDC and the LO field. When taking the quantum efficiency of 0.90 of the home-built homodyne detector into account, we result at a total receiver efficiency of about 0.72.

#### 3.3 Squeezing improvement

After successfully subtracting photons from both modes, the higher photon number components of even order become more dominant for the distilled states which is equivalent to an increase in the squeezing. We verify this by comparing the squeezing and anti-squeezing values of the initial and distilled quantum states. To do so, we use the 8000 undistilled pulses surrounding each photon subtraction event to serve as a phase reference to sort the data by their relative phase. This results in a squeezing variance of $-1.648 \pm 0.002$ dB relative to the shot noise for the undistilled state which fits well to the degradation of a 6 dB squeezed state after experiencing the losses of our setup. Due to the distillation step, this variance is increased to $-1.82 \pm 0.16$ dB below shot-noise. Note that these values have been calculated from the sorted raw data. Therefore, the statistics for the distilled events is limited thus increasing the standard deviation when compared to the undistilled squeezed state for which a greater number of states has been sampled.

Reconstructing the density matrices of both states using a maximum likelihood algorithm [37], we can directly see this improvement of the quadrature variance in the values of the density matrix elements as shown in Fig. 2. The increased contribution of the higher-order photon number components can also be visualised by plotting the Wigner functions in Fig. 3. The distilled state deviates from an elliptical shape, instead it is compressed in the $\hat {X}$ quadrature, representing the higher squeezing, and strongly elongated along the $\hat {P}$ quadrature. In order to understand these findings, we simulate our experimental data with our multi-mode theoretical model discussed earlier. These results are summed up in Table 1 and show a comparison of the reconstructed initial state $\hat {\rho }_{\mathrm {sq}}$ and the distilled state $\hat {\rho }_{\mathrm {dist}}$ with the simulation results $\hat {\rho }_{\mathrm {sq,sim}}$ and $\hat {\rho }_{\mathrm {dist,sim}}$. The figure of merit to compare theory and experiment is the fidelity between the measured and calculated density matrices, which shows a high agreement between experimental data and theoretical model. We want to emphasize here that the purity $\mathrm {tr}(\hat {\rho }^2)$ of the squeezed state can solely be explained by the optical losses present in the system. As a result of this, we are also able to see the additional drop in the purity value for the distilled state resulting from the non-Gaussian distillation process.

#### 3.4 Verifying the non-Gaussianity

For our experimental realization of the distillation, no negativity of the Wigner function is expected even for the lossless case. To quantify the non-Gaussianity of the distillation process further, we investigate the higher-order cumulants $\kappa _n, n \in \{1,2,3,4 \}$ of the Wigner function’s marginal distribution as functions of the quadrature phase, where the projected quadrature is $\hat {X}_{\theta } = \cos \theta \ \hat {X} +\sin \theta \ \hat {P}$ (see Fig. 4). In terms of the central moments $\mu _{n} = \left \langle \left ( \hat {X}_{\theta } - \langle \hat {X}_{\theta } \rangle \right )^n\right \rangle$ and the mean value $m_1$, the cumulants can be expressed in a compact form as

## 4. Conclusion

In conclusion, we have demonstrated the realization of a pulsed squeezing distillation experiment using an engineered source of squeezed states. In our experiment, source engineering allows us to generate pulsed single-mode squeezed states to efficiently combine the mode-resolving homodyne detection technique and temporal-mode-insensitive photon counting measurements in a squeezing distillation protocol without the need of strong optical filtering. We demonstrate the expected increase in squeezing due to the distillation via photon subtraction and verify the non-Gaussian nature of the distilled quantum state by calculating the higher-order statistical cumulants. We provide a simple theoretical model describing the photon subtraction process in the few-mode regime and find an excellent agreement between the theoretical predictions and the experimentally found quadrature variances. Our model also shows, that the achieved improvement due to the distillation is limited due to technical reasons and the remaining mode impurity. It has been shown recently that it is possible to precisely manipulate the mode purity by using more elaborate means of spectrally shaping the pump field [39] or the utilization of a special poling pattern of the non-linear crystal [40]. It is important to note, that the generated squeezing in this type of single-pass source is only limited by the amount of available pump energy. Experiments have reported the generation of states with a mean photon number $\langle \hat {n} \rangle = 80$ which translates to −25 dB of squeezing [41]. Finally, optimizing the losses of optics and detectors in both the photon subtraction and the homodyne detection part of the setup will increase the achievable distillation yield and open up the possibility to subtract more photons per mode. Assuming perfect optics and detectors, and only considering the for our setup imperfect HOM visibility, an initially squeezed state with −3 dB of squeezing will, after subtracting two photons from each, the signal and idler mode, would exhibit a quadrature variance of −3.19 dB relative to shotnoise. All-together, the compatibility of our source with both homodyne and photon counting measurements renders it a versatile source of optical quantum states in hybrid QIP and a candidate for the realization of larger quantum networks requiring the interference of several similar sources.

## Appendix

## A. Estimating the relative phase

The main challenge in order to reconstruct the density matrix of the state is to assign the correct LO phase for each measurement. Often, this is achieved via the precise control of the relative phase between the LO and the detected state. In our case, we are unable to control the relative phase, thus we use the undistilled squeezed states surrounding each distilled pulse as a phase reference. We illustrate this situation in Fig. 5. This has the advantage that the overall experimental setup is simplified and that the whole phase assignment can be done in post-processing assuming that all phase angles are covered equally.

To ensure an even coverage of all possible phase angles, we therefore scan the LO phase using a piezo-mounted mirror with a modulation frequency in the range of $10~\mathrm {Hz}$. We choose the measurement window of our oscilloscope to encompass $100~\mathrm {\mu s}$, roughly consisting of 8000 pulses, as we assume that it is much faster compared to the phase shift introduced by the piezo modulation and the random phase fluctuations, e.g. vibrating mirrors. In principle, different approaches to assign the LO phase are possible. We decided to use a MonteCarlo-like technique based on the experimental conditions to estimate the phase of each measurement which gives us a conservative estimate for the phase of each individual time trace of our measurement. This phase assignment procedure is described by the following steps:

- 1.
**Fitting the squeezing ellipse:**Assuming an equal sampling of all LO phases, we calculate the variance distribution by simulating the measurement of the marginal distribution $\hat {X}_{\theta } = \hat {X}\cos \theta +\hat {P}\sin \theta$ of a squeezed state using $\hat {X}$ and $\hat {P}$ as fitting parameter. We then choose the combination of $\hat {X}$ and $\hat {P}$ matching the measured data best, see Fig. 6(a). - 2.
**Phase-resolved simulation of the variances:**with the fitted squeezing and anti-squeezing values, we simulate homodyne measurements matching the experimental conditions, namely the detection of 8000 subsequent pulses and the calculation of their variance and the distribution under different phase angles $\mathrm {V}(\phi )$ assuming an uniform distribution of phases. - 3.
**Inversion of the variance distribution:**the distribution found in the previous step is inverted into a variance-dependent distribution of phases $\phi (\mathrm {V})$. - 4.
**Phase assignment for the individual measurement traces:**using the measured quadrature variances of the reference states, we randomly pick a possible phase from the phase distribution $\phi (\mathrm {V})$ to assign this measurement and the distilled pulse in its center a phase relative to the LO.

## B. Loss estimation

In our experiment, a multitude of individual parameters contribute to the overall loss figure. Most of them can be optimized in future experiments which would however not fit the scope of this first demonstration of source engineering in a hybrid QIP context. Most of the the individual contributions, such as the linear losses of the optics or the non-unity quantum efficiency of the photo diodes, can be remedied by further optimization. Other contributions cannot simply optimized by using different components. The HOM visibility of our source is intrinsically limited due to the phase matching conditions of KTP in this wavelength regime. Additionally, the the spatial overlap between LO mode and waveguide modes is limited due to differing geometrical profiles of the waveguide, which is close to quadratic, and a single-mode fiber, which is circular. Table 2 lists the individual loss contributions of our setup and the theoretically possible squeezing of a squeezed state with an initial squeezing of −6 dB.

## Funding

Horizon 2020 Framework Programme (820466 (CiViQ)) (665148 (QCUMbER)).

## Disclosures

The authors declare no conflict of interests.

## References

**1. **C. H. Bennett, H. J. Bernstein, S. Popescu, and B. Schumacher, “Concentrating partial entanglement by local operations,” Phys. Rev. A **53**(4), 2046–2052 (1996). [CrossRef]

**2. **H. Vahlbruch, M. Mehmet, K. Danzmann, and R. Schnabel, “Detection of 15 db squeezed states of light and their application for the absolute calibration of photoelectric quantum efficiency,” Phys. Rev. Lett. **117**(11), 110801 (2016). [CrossRef]

**3. **N. Takanashi, W. Inokuchi, T. Serikawa, and A. Furusawa, “Generation and measurement of a squeezed vacuum up to 100 mhz at 1550 nm with a semi-monolithic optical parametric oscillator designed towards direct coupling with waveguide modules,” Opt. Express **27**(13), 18900–18909 (2019). [CrossRef]

**4. **M. Förtsch, G. Schunk, J. U. Fürst, D. Strekalov, T. Gerrits, M. J. Stevens, F. Sedlmeir, H. G. L. Schwefel, S. W. Nam, G. Leuchs, and C. Marquardt, “Highly efficient generation of single-mode photon pairs from a crystalline whispering-gallery-mode resonator source,” Phys. Rev. A **91**(2), 023812 (2015). [CrossRef]

**5. **A. Otterpohl, F. Sedlmeir, U. Vogl, T. Dirmeier, G. Shafiee, G. Schunk, D. V. Strekalov, H. G. L. Schwefel, T. Gehring, U. L. Andersen, G. Leuchs, and C. Marquardt, “Squeezed vacuum states from a whispering gallery mode resonator,” Optica **6**(11), 1375–1380 (2019). [CrossRef]

**6. **A. Eckstein, A. Christ, P. J. Mosley, and C. Silberhorn, “Highly efficient single-pass source of pulsed single-mode twin beams of light,” Phys. Rev. Lett. **106**(1), 013603 (2011). [CrossRef]

**7. **G. Harder, V. Ansari, B. Brecht, T. Dirmeier, C. Marquardt, and C. Silberhorn, “An optimized photon pair source for quantum circuits,” Opt. Express **21**(12), 13975–13985 (2013). [CrossRef]

**8. **K.-i. Yoshino, T. Aoki, and A. Furusawa, “Generation of continuous-wave broadband entangled beams using periodically poled lithium niobate waveguides,” Appl. Phys. Lett. **90**(4), 041111 (2007). [CrossRef]

**9. **J. Eisert, S. Scheel, and M. Plenio, “Distilling gaussian states with gaussian operations is impossible,” Phys. Rev. Lett. **89**(13), 137903 (2002). [CrossRef]

**10. **B. Brecht, D. V. Reddy, C. Silberhorn, and M. G. Raymer, “Photon temporal modes: A complete framework for quantum information science,” Phys. Rev. X **5**(4), 041017 (2015). [CrossRef]

**11. **Y.-S. Ra, A. Dufour, M. Walschaers, C. Jacquard, T. Michel, C. Fabre, and N. Treps, “Non-gaussian quantum states of a multimode light field,” Nat. Phys. **16**(2), 144–147 (2020). [CrossRef]

**12. **J. Heersink, C. Marquardt, R. Dong, R. Filip, S. Lorenz, G. Leuchs, and U. L. Andersen, “Distillation of squeezing from non-gaussian quantum states,” Phys. Rev. Lett. **96**(25), 253601 (2006). [CrossRef]

**13. **R. Dong, M. Lassen, J. Heersink, C. Marquardt, R. Filip, G. Leuchs, and U. L. Andersen, “Experimental entanglement distillation of mesoscopic quantum states,” Nat. Phys. **4**(12), 919–923 (2008). [CrossRef]

**14. **B. Hage, A. Samblowski, J. DiGuglielmo, A. Franzen, J. Fiurásek, and R. Schnabel, “Preparation of distilled and purified continuous-variable entangled states,” Nat. Phys. **4**(12), 915–918 (2008). [CrossRef]

**15. **A. E. Ulanov, I. A. Fedorov, A. A. Pushkina, Y. V. Kurochkin, T. C. Ralph, and A. I. Lvovsky, “Undoing the effect of loss on quantum entanglement,” Nat. Photonics **9**(11), 764–768 (2015). [CrossRef]

**16. **U. L. Andersen, J. S. Neergaard-Nielsen, P. Van Loock, and A. Furusawa, “Hybrid discrete-and continuous-variable quantum information,” Nat. Phys. **11**(9), 713–719 (2015). [CrossRef]

**17. **O. Morin, K. Huang, J. Liu, H. Le Jeannic, C. Fabre, and J. Laurat, “Remote creation of hybrid entanglement between particle-like and wave-like optical qubits,” Nat. Photonics **8**(7), 570–574 (2014). [CrossRef]

**18. **H. Jeong, A. Zavatta, M. Kang, S.-W. Lee, L. S. Costanzo, S. Grandi, T. C. Ralph, and M. Bellini, “Generation of hybrid entanglement of light,” Nat. Photonics **8**(7), 564–569 (2014). [CrossRef]

**19. **G. Masada, K. Miyata, A. Politi, T. Hashimoto, J. L. O’brien, and A. Furusawa, “Continuous-variable entanglement on a chip,” Nat. Photonics **9**(5), 316–319 (2015). [CrossRef]

**20. **A. E. Ulanov, D. Sychev, A. A. Pushkina, I. A. Fedorov, and A. Lvovsky, “Quantum teleportation between discrete and continuous encodings of an optical qubit,” Phys. Rev. Lett. **118**(16), 160501 (2017). [CrossRef]

**21. **E. Agudelo, J. Sperling, L. Costanzo, M. Bellini, A. Zavatta, and W. Vogel, “Conditional hybrid nonclassicality,” Phys. Rev. Lett. **119**(12), 120403 (2017). [CrossRef]

**22. **J. Wenger, R. Tualle-Brouri, and P. Grangier, “Non-gaussian statistics from individual pulses of squeezed light,” Phys. Rev. Lett. **92**(15), 153601 (2004). [CrossRef]

**23. **Y. Kurochkin, A. S. Prasad, and A. Lvovsky, “Distillation of the two-mode squeezed state,” Phys. Rev. Lett. **112**(7), 070402 (2014). [CrossRef]

**24. **H. Takahashi, J. S. Neergaard-Nielsen, M. Takeuchi, M. Takeoka, K. Hayasaka, A. Furusawa, and M. Sasaki, “Entanglement distillation from gaussian input states,” Nat. Photonics **4**(3), 178–181 (2010). [CrossRef]

**25. **A. Christ, C. Lupo, M. Reichelt, T. Meier, and C. Silberhorn, “Theory of filtered type-ii parametric down-conversion in the continuous-variable domain: Quantifying the impacts of filtering,” Phys. Rev. A **90**(2), 023823 (2014). [CrossRef]

**26. **Y.-S. Ra, C. Jacquard, A. Dufour, C. Fabre, and N. Treps, “Tomography of a mode-tunable coherent single-photon subtractor,” Phys. Rev. X **7**(3), 031012 (2017). [CrossRef]

**27. **V. Ansari, G. Harder, M. Allgaier, B. Brecht, and C. Silberhorn, “Temporal-mode measurement tomography of a quantum pulse gate,” Phys. Rev. A **96**(6), 063817 (2017). [CrossRef]

**28. **T. J. Bartley and I. A. Walmsley, “Directly comparing entanglement-enhancing non-gaussian operations,” New J. Phys. **17**(2), 023038 (2015). [CrossRef]

**29. **A. Christ, K. Laiho, A. Eckstein, K. N. Cassemiro, and C. Silberhorn, “Probing multimode squeezing with correlation functions,” New J. Phys. **13**(3), 033027 (2011). [CrossRef]

**30. **A. B. U’Ren, C. Silberhorn, K. Banaszek, R. E. Ian A. Walmsley, W. P. Grice, and M. G. Raymer, “Generation of pure-state single-photon wavepackets by conditional preparation based on spontaneous parametric downconversion,” LaserPhysics **15**, 146 (2005).

**31. **G. Harder, T. J. Bartley, A. E. Lita, S. W. Nam, T. Gerrits, and C. Silberhorn, “Single-mode parametric-down-conversion states with 50 photons as a source for mesoscopic quantum optics,” Phys. Rev. Lett. **116**(14), 143601 (2016). [CrossRef]

**32. **R. J. A. Francis-Jones, R. A. Hoggarth, and P. J. Mosley, “All-fiber multiplexed source of high-purity single photons,” Optica **3**(11), 1270–1273 (2016). [CrossRef]

**33. **N. Bruno, A. Martin, T. Guerreiro, B. Sanguinetti, and R. T. Thew, “Pulsed source of spectrally uncorrelated and indistinguishable photons at telecom wavelengths,” Opt. Express **22**(14), 17246–17253 (2014). [CrossRef]

**34. **F. Kaneda, K. Garay-Palmett, A. B. U’Ren, and P. G. Kwiat, “Heralded single-photon source utilizing highly nondegenerate, spectrally factorable spontaneous parametric downconversion,” Opt. Express **24**(10), 10733–10747 (2016). [CrossRef]

**35. **W. P. Grice, A. B. U’Ren, and I. A. Walmsley, “Eliminating frequency and space-time correlations in multiphoton states,” Phys. Rev. A **64**(6), 063815 (2001). [CrossRef]

**36. **C. Silberhorn, P. K. Lam, O. Weiß, F. König, N. Korolkova, and G. Leuchs, “Generation of continuous variable einstein-podolsky-rosen entanglement via the kerr nonlinearity in an optical fiber,” Phys. Rev. Lett. **86**(19), 4267–4270 (2001). [CrossRef]

**37. **J. Řeháček, Z. c. v. Hradil, E. Knill, and A. I. Lvovsky, “Diluted maximum-likelihood algorithm for quantum tomography,” Phys. Rev. A **75**(4), 042108 (2007). [CrossRef]

**38. **C. R. Müller, G. Leuchs, C. Marquardt, and U. L. Andersen, “Optimally cloned binary coherent states,” Phys. Rev. A **96**(4), 042311 (2017). [CrossRef]

**39. **V. Ansari, E. Roccia, M. Santandrea, M. Doostdar, C. Eigner, L. Padberg, I. Gianani, M. Sbroscia, J. M. Donohue, L. Mancino, M. Barbieri, and C. Silberhorn, “Heralded generation of high-purity ultrashort single photons in programmable temporal shapes,” Opt. Express **26**(3), 2764–2774 (2018). [CrossRef]

**40. **C. Chen, J. E. Heyes, K.-H. Hong, M. Y. Niu, A. E. Lita, T. Gerrits, S. W. Nam, J. H. Shapiro, and F. N. C. Wong, “Indistinguishable single-mode photons from spectrally engineered biphotons,” Opt. Express **27**(8), 11626–11634 (2019). [CrossRef]

**41. **G. Harder, T. J. Bartley, A. E. Lita, S. W. Nam, T. Gerrits, and C. Silberhorn, “Single-mode parametric-down-conversion states with 50 photons as a source for mesoscopic quantum optics,” Phys. Rev. Lett. **116**(14), 143601 (2016). [CrossRef]