# Plug-&-play generation of non-Gaussian states of light at a telecom wavelength

Mohamed Melalkia, Tecla Gabbrielli, Antoine Petitjean, Léandre Brunel, Alessandro Zavatta, Sébastien Tanzilli, Jean Etesse, Virginia D'auria

► To cite this version:

Mohamed Melalkia, Tecla Gabbrielli, Antoine Petitjean, Léandre Brunel, Alessandro Zavatta, et al. Plug-&-play generation of non-Gaussian states of light at a telecom wavelength. Optics Express, 2022, 30 (25), pp.45195. 10.1364/OE.465980. hal-03699974

HAL Id: hal-03699974

https://hal.science/hal-03699974v1

Submitted on 16 Oct 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Plug-&-play generation of non-Gaussian states of light at a telecom wavelength

MOHAMED F. MELALKIA,¹,† TECLA GABBRIELLI,²,³,† ANTOINE PETITJEAN,¹ LÉANDRE BRUNEL,¹ ALESSANDRO ZAVATTA,²,³ SÉBASTIEN TANZILLI,¹ JEAN ETESSE,¹ AND VIRGINIA D'AURIA¹,⁴,*

¹Université Côte d'Azur, CNRS, Institut de Physique de Nice, Parc Valrose, 06108 Nice Cedex 2, France

²European Laboratory for Non-linear Spectroscopy (LENS), University of Florence, Via Nello Carrara 1, 50019 Sesto Fiorentino (FI), Italy

³Istituto Nazionale di Ottica (CNR-INO), CNR, Largo Enrico Fermi 6, 50125 Firenze, Italy

⁴Institut Universitaire de France, France

*virginia.dauria@univ-cotedazur.fr

**Abstract:** This work marks an important progress towards practical quantum optical technologies in the continuous variable regime, as it shows the feasibility of experiments where non-Gaussian state generation entirely relies on plug-&-play components from guided-wave optics technologies. This strategy is demonstrated experimentally with the heralded preparation of low amplitude Schrödinger cat states based on single-photon subtraction from a squeezed vacuum. All stages of the experiment are based on off-the-shelf fiber components. This leads to a stable, compact, and easily re-configurable realization, fully compatible with existing fibre networks and, more in general, with future out-of-the-laboratory applications.

## 1. Introduction

Optical quantum technologies in continuous variable (CV) regime find applications in quantum metrology [1,2], quantum communication [3-5], and quantum computing [6,7]. At the same time, many CV quantum information protocols requires the ability of generating non-Gaussian quantum states, such as Fock [8] and Schrödinger cat-like states [9], that are essential to fault-tolerant optical quantum information processing [10] or long distance quantum communication [11,12]. A common strategy to benefit from these resources is to implement conditional state preparation, often starting from a (Gaussian) squeezed state [13]. This idea has been demonstrated in multiple CV experimental realizations [8,9,14], including the recent trend of employing a nonlinear waveguide as squeezing source [15]. Reported realizations, however, rely entirely, or largely, on bulk-optics setups that are generally hardly scalable and require careful spatial alignments [13].

This work addresses the need for compact and stable realizations that can be easily assembled out-of-the-laboratory, as demanded to push current CV quantum optical experiments towards practical demonstrations. More specifically, it shows the feasibility of non-Gaussian state generation setups fully employing plug-&-play components from guided-wave technologies. Strategically, we applied this concept to single-photon subtraction from a squeezed vacuum state. Such a scheme, commonly used for the preparation of small amplitude Schrödinger cat preparation [13], can be considered as an archetypal configuration, gathering Gaussian and non-Gaussian resources and lying at the heart of multiple complex architectures [13]. Squeezed states at a telecom wavelength are generated by single-pass spontaneous parametric down conversion (SPDC) in a commercial periodically poled lithium niobate ridge waveguide (PPLN/RW), whose input and output are fiber-pigtailed. Both photon subtraction operation and heralded state characterization via homodyne detection (HD) are then implemented using standard fiber components from classical telecom technology. The clear non-Gaussian behaviour of the measured heralded state proves the pertinence of our approach. This result opens the way to robust CV optical realizations compatible with existing telecommunication fiber networks and with considerably reduced resource overhead, spatial alignment, and mode matching issues.
## 2. Experimental setup

<figure><img src="image_2.png" /><figcaption>Fig. 1. Experimental setup. A fiber-coupled continuous wave (CW) laser delivers two outputs at $\lambda = 1560.44$ nm and $\lambda = 780.22$ nm. The $\lambda = 780.22$ nm output passes through a polarization controller (PC) and a variable fiber optical attenuator (VOA) to be used as pump source for the PPLN/RW generating the single-mode squeezed vacuum. Downstream of the PPLN/RW, squeezed light is sent to a 95:5 fiber beam-splitter. The 5% output undergoes a spectral filtering stage (SFS) and is sent to a superconducting nanowire single-photon detector (SNSPD), whose signal heralds the non-Gaussian state generation. The 95% output is sent to a homodyne detection (HD) made of a 50:50 fiber beam-splitter and a commercial balanced detector (BD). The HD local oscillator is given by the laser output at $\lambda = 1560.44$ nm, passing through a fiber phase shifter (PS) and a polarization controller (PC).</figcaption></figure>

Fig. 1 shows the detailed experimental setup. The master source is a continuous wave (CW) fiber-coupled laser (Muquans) delivering outputs at $\lambda = 1560.44$ nm and its second harmonic (SH) at half of the wavelength ($\lambda = 780.22$ nm). The SH output is sent to fiber polarization controller (PC) and an optical variable attenuator (VOA) that allow adjusting its polarization and intensity. This set of plug-&-play telecom components is then directly connected to the PPLN/RW (NTT Electronics corp). Single-mode squeezed vacuum state at $\lambda = 1560.44$ nm is generated via type-0 frequency degenerate SPDC and collected by the optical fiber at the PPLN/RW output with an efficiency $\eta_c$. Propagation losses inside the nonlinear waveguide lead to a reduced PPLN/RW transmission efficiency $\eta_{pt}$. As largely documented in the literature [13], heralded single-photon subtraction can be implemented by sending squeezed light to a low reflectivity beam-splitter followed by a single-photon detection stage. In our experiment, a 95:5 fiber beam-splitter (f-BS) splits squeezed light into two paths: the 5% output is sent to the heralding path, whereas the 95% output, carrying the heralded non-Gaussian state, goes to the HD for the measurement of its CV properties. Following an approach that we already validated for squeezing experiments [16], the HD is realized by means of a 50:50 f-BS whose outputs are connected to a commercial balanced detection (Insight) showing a bandwidth of 300 MHz. The HD local oscillator (LO) is provided by the fundamental laser output ($\lambda = 1560.44$ nm), whose phase, $\theta$, is controlled via a fiber phase shifter. The overall detection efficiency is $\eta_{HD} = \eta_t \cdot \eta_{pd} \cdot \eta_{el} = 0.72$, where $\eta_t = 0.94$ includes measured losses in the f-BSs, $\eta_{pd} = 0.80$ is the photodiode efficiency as given by the manufacturer, and $\eta_{el} = 0.96$ accounts to the residual HD electronic noise. The HD output
is monitored using an oscilloscope (Teledyne LeCroy) for data post-processing. The CV data acquisition on the oscilloscope is triggered by photon counting events from the heralding path. Light on the heralding path undergoes a cascade of optical spectral filters allowing to comply with the strong mismatch between the HD electronic bandwidth (300 MHz) and the broad squeezing optical emission bandwidth (a few THz). Indeed, photon-counting events of wavelengths far away from the CW LO could trigger the CV detection of signals whose non-Gaussian features are out-of-the-reach of HD electronics. This would weave data corresponding to the desired non-Gaussian states with those relative to squeezed vacuum, eventually resulting in the measurement of a mixed state. Such an effect further combines with the fact that single-pass squeezing can be not pure over its whole bandwidth [17]. Spectral filtering in the heralding path allows circumventing these imperfections [17, 18]. In our realization, the 5% of tapped light from the squeezed state is filtered using a commercial 500 MHz FWHM fiber Bragg grating filter (FBG from AOS GmbH) and an in-house fiber-coupled Fabry-Pérot cavity with a target bandwidth of 10 MHz. Note that, as an alternative, plug-&-play optical filters with bandwidths of a few tens of MHz are commercially available. Filtered light is then sent to a polarization controller and to a superconducting nanowire single-photon detector (SNSPD from IDQuantique). The heralding rate is of 3 kHz and the dark count rate is about 80 Hz, giving a signal to noise ratio, SNR, of 15.7 dB.

### 3. Data analysis

The properties of the heralded non-Gaussian state can be reconstructed by performing a quantum state tomography, starting from the state quadratures measured by the HD [19]. A series of 43000 homodyne traces, $x_i(t)$, of a duration of 5 µs each at a sampling rate of 5 Gs/s are recorded using the oscilloscope. Note that the LO phase is spanned with a speed of approximately $9\pi/s$ and, therefore, it can be considered as a constant over the time needed to acquire a homodyne trace. For each acquisition $x_i(t)$, the actual $\theta_i$ value is inferred from the phase dependency of the measured squeezed vacuum variance.

As discussed in [17], in the CW regime, the quadrature $X_i$ of the heralded non-Gaussian state can be obtained from raw homodyne data as [20]:

$$X_i = \int f(t)x_i(t)dt, \quad (1)$$

where $f(t)$ is the temporal mode of the heralded non-Gaussian state, as determined by the filtering stages on the heralding path. When the squeezing bandwidth is much broader than the bandwidth of the filtering cavity, $\gamma/\pi$, the theoretical $f(t)$ is expected to be $\sim e^{\gamma t}u(-t)$, $u(t)$ being the Heaviside step function. Experimentally, $f(t)$ is obtained as in [21] by computing the autocorrelation function of the heralded homodyne traces, $\langle x(t)x(t') \rangle$, and by numerically performing its single value decomposition: the eigenfunction corresponding to the highest eigenvalue gives the optimal $f(t)$. Fig. 2 shows both the reconstructed temporal mode $f(t)$ (blue) and theoretical one (red); a perfect agreement between the two curves can be observed and gives a filtering bandwidth of $\gamma/\pi = 9.75$ MHz at FWHM, in agreement with the filtering stage characteristics.

By exploiting Eq. (1), the obtained $f(t)$ is used to extract the quadratures of the squeezed state and of the heralded non-Gaussian state out of the raw homodyne data. More in details, with reference to Fig. 2, $x_i(t)$ data are post-processed within the time window $\approx [-1.25 \mu s, -1 \mu s]$ to extract the heralded non-Gaussian state quadratures. Data outside this time window are used to extract the squeezed state quadrature.
<figure><img src="image_3.png" /><figcaption>Fig. 2. The temporal mode function corresponding to the highest eigenvalue (blue curve) and the fit with the theoretical model (red curve) as functions of time. The cascade of the Fabry-Pérot cavity and the fiber Bragg grating filter in the heralding path gives a filtering bandwidth of 9.75 MHz at FWHM.</figcaption></figure>

4. Experimental results and discussion

The tomographic analysis allows measuring, for the non-heralded state, a squeezing of −1.80±0.05 dB and an antisqueezing of 3.36±0.05 dB. This corresponds to an initial squeezing of −5.39±0.05 dB when corrected for the total losses experienced by squeezed light, $\eta_{tot} = \eta_{wg} \cdot \eta_s \cdot \eta_{HD} = 0.48$. In the expression of $\eta_{tot}$, $\eta_s = 0.96$ is the transmission of the subtracting 95:5 f-BS and $\eta_{wg}$ accounts for both the propagation losses and the guide-to-fibre coupling occurring at the SPDC stage ($\eta_{wg} = \eta_{pl} \cdot \eta_c$). From the measured $\eta_s$ and $\eta_{HD}$, we deduce $\eta_{wg} = 0.69$.

The tomography of the heralded non-Gaussian state is performed with the maximum-likelihood algorithm with 200 iterations [22]. Fig. 3-(a) shows the obtained Wigner function of the generated non-Gaussian state. As can be seen, the generated state has a highly non-Gaussian Wigner function with a minimum at phase origin $W_{NG}(0,0) = 0.016 \pm 0.004$. Its profile perfectly matches the one predicted by theory ($W_{Th}(0,0) = 0.015$) following an approach similar to the one in [23, 24] and taking into account all the experimental imperfections. This result represents the very first demonstration of non-Gaussian states in CV regime obtained using a fully guided-wave approach and it attests its impact towards compact, simple, and plug-&-play setups involving both Gaussian and non-Gaussian operations.

<figure><img src="image_4.png" /><figcaption>Fig. 3. (a) The Wigner function $W_{NG}(x, p)$ of the heralded non-Gaussian state and (b) the Wigner function $W_{Corr}(x, p)$ of the state corrected by the homodyne detection efficiency $\eta_{HD} = 0.72$.</figcaption></figure>

In this regard, we observe that without modifying $\eta_{wg}$, which is in agreement with realistic
performances in such fibre-coupled ridge waveguides, our results could be straightforwardly
enhanced by simply replacing our balanced photodetection with one with higher efficiency. By
correcting the data by the HD efficiency ($\eta_{HD} = 0.72$) [19], the non-Gaussian state obtained just
after the 95:5 f-BS exhibits a negative Wigner function as expected for Schrödinger cat states, the
corresponding value at the phase origin being $W_{Corr}(0,0) = -0.065 \pm 0.004$ (see Fig. 3-(b)).
We point out that residual losses in the fiber components and low HD efficiency are not intrinsic
to the followed guided-wave approach but only to our specific realization. In particular, the lack
of efficient photodiodes, that represent the most critical loss factor in the HD, affects in the same
way both guided-wave and bulk realizations.

As concluding remark, we stress that a further improvement of the Wigner negativity could
be obtained by improving $\eta_{wg}$. The combined effect of propagation and coupling losses at
the nonlinear waveguide leads to a degradation of the purity of the squeezed state on which
the subtraction protocol is applied. This can be circumvented by adopting a full integrated
setup including squeezed state generation, photon subtraction and homodyne detection, in a
configuration similar to the one shown for Gaussian state generation, manipulation, and detection
on lithium niobate circuits [25–27]. Such a strategy would push the propagation losses down to
0.03 dB/cm and the coupling losses to 0, thus giving for a typical chip length of 5 cm, a total
$\eta_{wg} = 0.97$. In this condition, the Wigner function negativity would go to its maximal value,
$W_{Corr}(0,0) = -0.22$.

5. Conclusion

We have demonstrated that a plug-&-play, fully guided-wave, approach to CV experiments is
definitely reachable with off-the-shelves standard components and devices. This concept has been
pushed, for the first time, to the generation of non-Gaussian states of light, showing the possibility
of implementing complex setups free from any bulk optics stage. The obtained realization enables
reconfigurable and scalable systems compatible with integrated photonics [25–27], adapted to
ambitious CV quantum communication protocols. These ingredients are essential to scale-up
realizations of increasingly complex optical quantum systems capable of enabling disruptive but
realistic quantum information protocols.

**Funding.** This work has been conducted within the framework of the project OPTIMAL granted by the European Union by means of the "Fond Européen de développement regional" (FEDER). The authors also acknowledge financial support from the Agence Nationale de la Recherche (ANR) through the project HyLight (ANR-17-CE30-0006-01) and the French government through its program "Investments for the Future" under the Université Côte d'Azur UCA-JEDI project managed by the ANR (grant agreement ANR-15-IDEX-01).

**Acknowledgments.** Virginia D'Auria acknowledges financial support from the Institut Universitaire de France (IUF).

**Disclosures.** The authors declare no conflicts of interest.

**Data availability.** Data underlying the results presented in this paper are not publicly available at this time but may be obtained from the authors upon reasonable request.

†These authors contributed equally to this work.

**References**

1. R. Schnabel, "Squeezed states of light and their applications in laser interferometers," Phys. Reports **684**, 1–51 (2017).

2. Y. Wang and K. Fang, "Continuous-variable graph states for quantum metrology," Phys. Rev. A **102**, 052601 (2020).

3. S. L. Braunstein and P. Van Loock, "Quantum information with continuous variables," Rev. modern physics **77**, 513 (2005).
4. J. Dias, M. S. Winnel, N. Hosseinidehaj, and T. C. Ralph, "Quantum repeater for continuous-variable entanglement distribution," Phys. Rev. A 102, 052425 (2020).

5. T. Eberle, V. Händchen, J. Duhme, T. Franz, F. Furrer, R. Schnabel, and R. F. Werner, "Gaussian entanglement for quantum key distribution from a single-mode squeezing source," New J. Phys. 15, 053049 (2013).

6. W. Asavanant and A. Furusawa, "Optical quantum computers: A route to practical continuous variable quantum information processing," (2022).

7. U. L. Andersen, G. Leuchs, and C. Silberhorn, "Continuous-variable quantum information processing," Laser & Photonics Rev. 4, 337-354 (2010).

8. O. Morin, V. d'Auria, C. Fabre, and J. Laurat, "High-fidelity single-photon source based on a type ii optical parametric oscillator," Opt. letters 37, 3738-3740 (2012).

9. A. Ourjoumtsev, H. Jeong, R. Tualle-Brouri, and P. Grangier, "Generation of optical 'schrödinger cats' from photon number states," Nature 448, 784-786 (2007).

10. J. Niset, J. Fiurášek, and N. J. Cerf, "No-go theorem for gaussian quantum error correction," Phys. review letters 102, 120501 (2009).

11. 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, 178-181 (2010).

12. A. Tipsmark, J. S. Neergaard-Nielsen, and U. L. Andersen, "Displacement-enhanced entanglement distillation of single-mode-squeezed entangled states," Opt. Express 21, 6670-6680 (2013).

13. A. Lvovsky, P. Grangier, A. Ourjoumtsev, V. Parigi, M. Sasaki, and R. Tualle-Brouri, "Production and applications of non-gaussian quantum states of light," arXiv preprint arXiv:2006.16985 (2020).

14. W. Asavanant, K. Nakashima, Y. Shiozawa, J.-I. Yoshikawa, and A. Furusawa, "Generation of highly pure schrödinger's cat states and real-time quadrature measurements via optical filtering," Opt. Express 25, 32227-32242 (2017).

15. K. Takase, A. Kawasaki, B. K. Jeong, M. Endo, T. Kashiwazaki, T. Kazama, K. Enbutsu, K. Watanabe, T. Umeki, S. Miki, H. Terai, M. Yabuno, F. China, W. Asavanant, J. Ichi Yoshikawa, and A. Furusawa, "Generation of schrödinger cat states with wigner negativity using a continuous-wave low-loss waveguide optical parametric amplifier," Opt. Express 30, 14161-14171 (2022).

16. F. Kaiser, B. Fedrici, A. Zavatta, V. d'Auria, and S. Tanzilli, "A fully guided-wave squeezing experiment for fiber quantum networks," Optica 3, 362-365 (2016).

17. J.-i. Yoshikawa, W. Asavanant, and A. Furusawa, "Purification of photon subtraction from continuous squeezed light by filtering," Phys. Rev. A 96, 052304 (2017).

18. M. Melalkia, L. Brunel, S. Tanzilli, J. Etesse, and V. D'Auria, "Theoretical framework for photon subtraction with non-mode-selective resources," Phys. Rev. A 105, 013720 (2022).

19. U. Leonhardt, *Measuring the quantum state of light*, vol. 22 (Cambridge university press, 1997).

20. A. I. Lvovsky, "Squeezed light," Photonics: Sci. Foundations, Technol. Appl. 1, 121-163 (2015).

21. O. Morin, C. Fabre, and J. Laurat, "Experimentally accessing the optimal temporal mode of traveling quantum light states," Phys. review letters 111, 213602 (2013).

22. A. I. Lvovsky and M. G. Raymer, "Continuous-variable optical quantum-state tomography," Rev. modern physics 81, 299 (2009).

23. A. Ourjoumtsev, "Étude théorique et expérimentale de superpositions quantiques cohérentes et d'états intriqués non-gaussiens de la lumière," Ph.D. thesis, Université Paris Sud-Paris XI (2007).

24. O. Morin, "Non-gaussian states and measurements for quantum information," Ph.D. thesis, Université Pierre et Marie Curie-Paris VI (2013).

25. F. Mondain, T. Lunghi, A. Zavatta, E. Gouzien, F. Doutre, M. De Micheli, S. Tanzilli, and V. D'Auria, "Chip-based squeezing at a telecom wavelength," Photonics Res. 7, A36-A39 (2019).

26. F. Lenzini, J. Janousek, O. Thearle, M. Villa, B. Haylock, S. Kasture, L. Cui, H.-P. Phan, D. V. Dao, H. Yonezawa et al., "Integrated photonic platform for quantum information with continuous variables," Sci. advances 4, eaat9331 (2018).

27. Y. Zhao, Y. Okawachi, J. K. Jang, X. Ji, M. Lipson, and A. L. Gaeta, "Near-degenerate quadrature-squeezed vacuum generation on a silicon-nitride chip," Phys. Rev. Lett. 124, 193601 (2020).