Arxiv:2010.01203V1 [Quant-Ph] 2 Oct 2020 That Acts on the Qubit

Full Thermalization of a Photonic Qubit

A. G. de Oliveira,1 R. M. Gomes,2 V. C. C. Brasil,1 N. Rubiano da Silva,1 L. C. Céleri,2, 3 and P. H. Souto Ribeiro1, ∗ 1Departamento de F´ısica, Universidade Federal de Santa Catarina, CEP 88040-900, Florianópolis, SC, Brazil 2Institute of Physics, Federal University of Goiás,74690-900, Goiânia,GO, Brazil 3Department of Physical Chemistry, University of the Basque Country UPV/EHU, Apartado 644, E-48080 Bilbao, Spain The generalized amplitude damping (GAD) quantum channel implements the interaction between a qubit and an environment with arbitrary temperature and arbitrary interaction time. Here, we implement a photonic version of the GAD for the case of infinite interaction time (full thermalization). We also show that this quantum channel works as a thermal bath with controlled temperature. Keywords: Parametric down-conversion; Quantum maps; Quantum Thermodynamics; Photonic qubits.

I. INTRODUCTION Choosing the vertical (V ) polarization state to represent the excited state with energy 2, and the horizontal (H) Photon polarization is a very useful degree of freedom polarization state to represent the ground state with en- for the realization of several tasks in Quantum Optics ergy 1, the thermal state is written as the density oper- and Quantum Information. Having the structure of a ator: two-level system, it has been widely used in experiments to test Bell inequalities [1, 2], it has inspired the creation ρˆ = P |HihH| + P |V ihV | (1) of the paradigmatic BB84 quantum cryptography proto- thermal H V e−β1 e−β2 col [3], and it was used to demonstrate quantum tele- = | ih | + | ih |, portation [4], superdense coding [5], one-way quantum Z 1 1 Z 2 2 computing [6], and several other interesting experiments. where P = (1 − P ) ∈ [0, 1], β = 1/k T , k More recently, we have witnessed the increasing inter- H V B B is the Boltzmann constant, T is the temperature and est in the field of Quantum Thermodynamics, which at- Z = Σ e−βi acts as the partition function of the sys- tempts to combine concepts from Thermodynamics and i tem. Classical and Quantum Information theories [7, 8]. Sev- The identification between the two-level system ther- eral experimental platforms have been demonstrated for mal state and the polarization mixed state in Eq. (1) testing new ideas in this field [9–13], including the pho- allows us to obtain the effective inverse temperature β as tonic polarization [14]. a function of the populations P and P : The interaction between a single qubit and a ther- H V mal bath was interpreted as a Generalized Amplitude −1 Damping (GAD) channel by S. Jevtic et al. to inves- β = ε ln (PH /PV ), (2) tigate quantum thermometers [15]. The experimental implementation of this model of interaction using pho- where ε = (2 − 1) is the qubit energy gap. tonic polarization was presented by W. K. Tham et al. [16] and L. Mancino et al. [17]. This quantum channel is often understood as the interaction between a qubit III. GENERALIZED AMPLITUDE DAMPING and an environment with arbitrary temperature. While CHANNEL the GAD also allows control of the interaction time between qubit and environment, we realize here a simplified The action of an environment at some fixed tempera- scheme where the interaction time is set as infinite. We ture acting on a qubit system is described by the action interpret it as a thermal bath with tunable temperature of the completely positive and trace preserving map E: arXiv:2010.01203v1 [quant-ph] 2 Oct 2020 that acts on the qubit. The experimental scheme is sim- 3 ple, reliable and can be used in the study of quantum X ρˆ = E ρˆ = Γˆ ρˆ Γˆ† , (3) systems undergoing thermodynamic processes and heat out in k in k engines. k=0 whereρ în andρ ôut are the input and output states of the system and the Kraus operators Γˆ are given by II. PHOTON POLARIZATION THERMAL k STATES 1 0 0 0 Γˆ = p1 − ξ √ , Γˆ = pξ √ , 0 0 1 − p 2 p 0 √ √ We can associate a thermal state, often called Gibbs 0 p 1 − p 0 Γˆ = p1 − ξ , Γˆ = pξ . (4) state, to two-level systems like the photon polarization. 1 0 0 3 0 1

where ξ ∈ [0, 0.5] and p = (1 − e−λt) ∈ [0, 1] is the de- ∗ [email protected] cay probability associated with the excited state, while λ 2 stands for the damping constant. This association is use- interferometer there is a half-wave plate that makes the ful when the two-level system is related to some physical operation |V i → |θi and |Hi → |φi. The polarization quantity that is subject to spontaneous decay. However, components are recombined in the second PBS. Since this is not the case for the photon polarization. the interferometer is unbalanced, the recombination is If we consider an initial state in the general form incoherent and the states at the outputs of the PBS are mixed states given by: ρ ρ ρˆ = 00 01 , (5) in ρ10 ρ11 1 cos2 φ 0 ρˆ = ; (10) the output state will be s 0 cos2 θ cs H/V √ 1 − pξ −√ρ11(1 − p) ρ01 1 − p ρôut = , (6) 2 ρ10 1 − p ρ11(1 − p) + pξ 1 sin θ 0 ρˆ = . (11) l c 0 sin2 φ where we have employed the normalization condition l H/V ρ00 +ρ11 = 1. Note that in the infinite time limit t → ∞, where cs and cl are normalization factors: i.e., p → 1, the output state of the channel (GADt→∞) 2 2 is given by: cs = cos φ + cos θ; (12) c = sin2 φ + sin2 θ. 1 − ξ 0 l ρˆ = , (7) out 0 ξ After the interferometer, the trace over the spatial de- grees of freedom is made by detection of both outputs and ξ = (1 + eβε)−1 is now the thermal population of with the same detector. In practice, we use the measure- the excited state, with β being the inverse temperature ment scheme described in Ref. [18]. The final state is defined in Eq. (1). In what follows, we describe the given by experimental implementation of this channel. 1 cos2 φ + sin2 θ 0 ρôut = 2 2 , (13) 2 0 cos θ + sin φ H/V IV. IMPLEMENTING THE PHOTONIC GADt→∞ CHANNEL and the inverse temperature (assuming ε = 1) is:

2 2 A possible experimental scheme for photonic qubits is cos φ + sin θ βout = ln . (14) displayed in Fig. 1. We will show that it implements cos2 θ + sin2 φ the GADt→∞. Let us consider that the input state is incident from the left and is described by the general The identification with the action of the GADt→∞ chan- density matrix of Eq. (5) in the H-V basis. The state nel is made by comparing the density matrices in Eqs. travels through the first unbalanced Mach-Zehnder inter- (7) and (13): ferometer. The path difference is set to be larger than 1 − ξ = (cos2 φ + sin2 θ)/2 ; (15) the coherence length of the input state, so that the co- 2 2 herences of the state are destroyed. Therefore, the state ξ = (cos θ + sin φ)/2. evolves to: Here, we implement a variant of this scheme employ- ing photons produced by a polarization-entangled photon source. In our scheme, the input photonic qubits ρ00 0 ρˆ1 = . (8) are prepared in thermal states with easily controllable 0 ρ11 H/V populations, without the need of the first unbalanced in- In the next step it is transmitted through a half-wave terferometer of Fig. 1. plate (HWP) that rotates the polarization states |Hi → |Di (diagonal) and |V i → |Ai (anti-diagonal), so that V. EXPERIMENT the state is transformed to: The experimental setup is shown in Fig. 2. A ρ 0 ρˆ = 00 , (9) diode laser oscillating at 405 nm is used to pump 2 0 ρ 11 D/A a two-crystal-sandwich spontaneous parametric down- conversion (SPDC) source. This configuration is fre- where the basis was changed to diagonal/anti-diagonal quently used as a source of polarization-entangled photon (D/A). pairs [19]. The idler photon is directed to a single-photon Next, the state enters the second unbalanced interfer- counting module (SPCM) and detected without measur- ometer through a polarizing beam splitter (PBS) that ing its polarization, therefore tracing out this degree of splits the polarization components. In each arm of the freedom. In this case, the detection of an idler photon 3

HWP(θ/2) M

measurement

HWP(π/8) HWP(ϕ/2) HWP detector

PBS

QWP HWP(π/4)

FIG. 1. Scheme for implementing the GADt→∞. See main text for details.

HWP(θ/2) M pump D1 laser HWP

HWP(ϕ/2) HWP(π/8) H V HWP nal PBS sig HWP(π/4) QWP CC 405 nm idler 810 nm nonlinear crystals D2

FIG. 2. Experimental setup. The two-crystal spontaneous parametric down-conversion source is pumped by a blue laser and signal and idler photons are emitted. The idler photon is directed to a single photon counting module and the signal photon is directed to two nested unbalanced Mach-Zehnder interferometers. heralds a signal photon remotely prepared in a mixed nal and idler photons. The optical devices indicated in state for the polarization. As we have seen, a mixed Fig. 2 are half-wave plates (HWP), mirrors (M), polariz- polarization state can be interpreted as a thermal state ing beam splitters (PBS), interference filters (F), SPCM for which the temperature depends only on the popula- detectors (D1 and D2) and coincidence counting elec- tions of linear vertical (V) and horizontal (H) polarization tronics (CC), while (H) and (V) indicate the pump beam states. polarization direction that interacts with each crystal. The SPDC source consists of two Beta-Barium-Borate The overall scheme works as follows. The blue laser (BBO) nonlinear crystals. Signal and idler photons are pumps the crystals, producing entangled photon pairs. collected along directions such that they have the same The idler photons are detected and their polarization de- wavelength of 810 nm and interference filters centered at grees of freedom are traced out. This prepares the signal this wavelength and having 10-nm bandwidth are placed photons remotely in a mixed state. By adjusting the in front of the SPCMs. pump half-wave plate before the crystals, it is possible to In the path of the idler beam there are two nested, un- control the populations of the idler photon density ma- balanced Mach-Zehnder interferometers working in the trix and therefore the temperature of the corresponding non-interfering regime, because the path length differ- thermal state. In our experiment we will prepare all in- ences are bigger than the coherence length of the sig- put states for the quantum channel in an initial thermal 4

trolled temperature is not obvious. Our results show that the GAD is a good option for accomplishing this task. Moreover, in the cases where the dynamics of the coherences do not play an important role, we show that the GADt→∞ works very well. For instance, for realizing heat engines using the photon polarization as the working substance, one can use the GADt→∞. For more general applications, the GAD with variable interaction time is necessary and the experimental implementation requires one additional qubit.

VII. CONCLUSION

FIG. 3. Experimental results. The effective inverse tempera- In conclusion, we present an experimental realization ture βout of the output states after interaction with the heat bath is plotted as function of θ for four values of φ, and five of a photonic qubit quantum channel that acts like input states. a heat bath with controlled temperature and infinite interaction time. We perform measurements using heralded single photons that illustrate the operation of state and in this case the first unbalanced interferometer the scheme in the context of this interpretation. We test in Fig. 1 is not necessary. The heralded signal photons five different input thermal states, including the pure prepared in thermal states are incident in the sequence states |Hi and |V i and we find that at the output the of two nested, unbalanced Mach-Zehnder interferometers temperature of the system is only determined by the that implement the GADt→∞ and polarization analysis. parameters of the heat bath. We believe that this setup can be a useful tool in the study of thermodynamic processes in the limit of small scale quantum systems. VI. RESULTS

The experimental results are displayed in Fig. 3. The ACKNOWLEDGMENTS effective inverse temperature βout is plotted as a function of θ for four settings of φ. Five different input states were prepared in thermal states with different effective inverse The authors would like to thank the Brazilian Agencies temperatures, including the pure states |Hi (β = ∞) CNPq, FAPESC, FAPEG, and the Brazilian National and |V i (β = −∞), and the measured output inverse Institute of Science and Technology of Quantum Infor- temperatures are shown. The results indicate that the mation (INCT/IQ). This study was financed in part by output temperatures are independent of the input state, the Coordena¸cãode Aperfei¸coamento de Pessoal de N´ıvel as expected for the interaction with a heat bath during a Superior - Brasil (CAPES) - Finance Code 001. LCC time tending to infinite. We note that there are negative would like to also acknowledge support from Spanish temperatures, because in this system it is very easy to MCIU/AEI/FEDER (PGC2018-095113-B-I00), Basque have population inversion, in the sense that we can easily Government IT986-16, the projects QMiCS (820505) and prepare a state with larger population in the excited state OpenSuperQ (820363) of the EU Flagship on Quantum |V i than the ground state |Hi. Technologies and the EU FET Open Grant Quromorphic Even though two-level systems have a simple de- and the U.S. Department of Energy, Office of Science, Of- scription and the photon polarization can be controlled fice of Advanced Scientific Computing Research (ASCR) and measured with relative simplicity, the experimental quantum algorithm teams program, under field work pro- configuration that implements a heat bath with con- posal number ERKJ333.

[1] N. D. Mermin, Rev. Mod. Phys. 65, 803 (1993), URL [4] D. Bouwmeester, J. Pan, K. Mattle, M. Eibl, H. Wein- https://link.aps.org/doi/10.1103/RevModPhys.65. furter, and A. Zeilinger, Nature 390, 575 (1997). 803. [5] K. Mattle, H. Weinfurter, P. Kwiat, and A. Zeilinger, [2] N. Brunner, D. Cavalcanti, S. Pironio, V. Scarani, and Phys. Rev. Lett. 76, 4656 (1996). S. Wehner, Rev. Mod. Phys. 86, 419 (2014), URL https: [6] P. Walther, K. J. Resch, T. Rudolph, E. Schenck, H. We- //link.aps.org/doi/10.1103/RevModPhys.86.419. infurter, V. Vedral, M. Aspelmeyer, and A. Zeilinger, Na- [3] N. Gisin, G. Ribordy, W. Tittel, and H. Zbinden, Rev. ture 434, 169 (2005). Mod. Phys. 74, 145 (2002), URL https://link.aps. [7] R. Uzdin, A. Levy, and R. Kosloﬀ, Phys. Rev. X org/doi/10.1103/RevModPhys.74.145. 5, 031044 (2015), URL https://link.aps.org/doi/10. 1103/PhysRevX.5.031044. 5

[8] M. Lostaglio, K. Korzekwa, D. Jennings, and Journal of Physics 49, 783 (2019), 1905.02829. T. Rudolph, Phys. Rev. X 5, 021001 (2015), URL https: [15] S. Jevtic, D. Newman, T. Rudolph, and T. M. Stace, //link.aps.org/doi/10.1103/PhysRevX.5.021001. Phys. Rev. A 91, 012331 (2015), URL https://link. [9] J. Liphardt, S. Dumont, S. B. Smith, I. J. Tinoco, and aps.org/doi/10.1103/PhysRevA.91.012331. C. Bustamante, Science 296, 1832 (2002). [16] W. Tham, H. Ferretti, A. Sadashivan, and A. Steinberg, [10] N. C. Harris, Y. Song, and C.-H. Kiang, Phys. Rev. Lett. Scientific reports 6, 1 (2016). 99, 068101 (2007). [17] L. Mancino, M. Sbroscia, I. Gianani, E. Roccia, and [11] S. Toyabe, T. Sagawa, M. Ueda, E. Muneyuki, and M. Barbieri, Physical Review Letters 118, 130502 (2017). M. Sano, Nat. Phys. 6, 988 (2010). [18] A. Salles, F. de Melo, M. P. Almeida, M. Hor-Meyll, S. P. [12] S. An, J.-N. Zhang, M. Um, D. Lv, Y. Lu, J. Zhang, Walborn, P. H. Souto Ribeiro, and L. Davidovich, Phys. Z.-Q. Yin, H. T. Quan, and K. Kim, Nat. Phys. 11, 193 Rev. A 78, 022322 (2008), URL https://link.aps.org/ (2014). doi/10.1103/PhysRevA.78.022322. [13] R. M. de Araújo,T. Häffner, R. Bernardi, D. S. Tasca, [19] P. G. Kwiat, E. Waks, A. White, I. Appelbaum, and P. H. M. P. J. Lavery, M. J. Padgett, A. Kanaan, L. C. Céleri, Eberhard, Phys. Rev. A 60, R773 (1999), URL https: and P. H. S. Ribeiro, J. Phys. Comm. 2, 035012 (2018). //link.aps.org/doi/10.1103/PhysRevA.60.R773. [14] G. L. Zanin, T. Häffner,M. A. A. Talarico, E. I. Duzzioni, P. H. S. Ribeiro, G. T. Landi, and L. C. Céleri,Brazilian