1 | Indonesian Journal of Science & Technology, 6 Issue April 2021 Hal 183-192 Indonesian Journal of Science &Volume Technology 6 . 1, . Indonesian Journal of Science & Technology Journal homepage: http://ejournal. edu/index. php/ijost/ Decay of Entanglement of Correlated Qubits Through Bosonic Fields Deniz Tyrkpenye Department of Electrical Engineering, stanbul Technical University, 34467 stanbul. Turkey Correspondence E-mail: dturkpence@itu. ABSTRACT Distribution of entangled parties with the longest time possible is of importance to quantum communication. Therefore analysing the decay character of entanglement of correlated qubits in the presence of reservoir effects is of significance to the quantum-based technologies. This study covers the analysis of the temporal entanglement decay of two maximally entangled qubits against different reservoir types and system parameters. It's shown how varying the coupling type of the system to the environment affects the lifetime of entanglement. It's also found that in the presence of quantum interaction between entangled qubits, it's possible to enlarge the entanglement lifetime depending on the initialization of entanglement. Model parameters used in the numerical calculations and the results are general enough to be applied in any specific quantum-based experimental task. A 2021 Tim Pengembang Jurnal UPI ARTICLE INFO Article History: Submitted/Received 11 Aug 2020 First revised 18 Oct 2020 Accepted 20 Des 2021 First available online 04 Feb 2021 Publication date 01 Apr 2021 ____________________ Keywords: Quantum entanglement. Quantum decoherence. Open quantum system. Tyrkpenye. Decay of Entanglement of Correlated Qubits Through Bosonic Field. INTRODUCTION Quantum based technologies exploit quantum resources which have no classical Entanglement (Horodecki, et , 2. is an indispensable quantum resource with many useful applications such as quantum communication (Bennett, et al. Kimble, 2. , secure communication (Lo, et al. , 2. , or dense coding (Wang, et , 2. The major limitation of quantumbased technologies is that the quantum properties are fragile and become classical This quantum to classical Much attention has been paid for suppressing decoherence and extending quantum properties to useful timescales such as dynamical control methods . ynamical decouplin. to decouple the system from the environment (Suter and Alvarez, 2. However, these methods are costly as they require engineered decoupling Alongside these methods, the entangled systems could be analysed by the system parameters or the coupling mechanisms to the environment. By these analyses, the quantum systems could be prepared by appropriate initializations in which the quantum properties last longer. In this study, a maximally entangled twoqubit system in which one of the qubits is in contact with the reservoir . was analysed under different coupling schemes or system parameters. More specifically, the entanglement dynamics of two maximally entangled qubits explored under different coupling and initialization schemes. One of the qubits was assumed to be perfectly isolated from the environment while the other one is coupled to a bosonic reservoir. It's shown that a specific coupling type is detrimental for the lifetime of entanglement. On the other hand, it's also shown that the initialization of entanglement could extend the entanglement lifetime if there exists a quantum interaction between the qubits. The current-state-of the-art on the subject mostly show up in the field of quantum The idea 'quantum internet' coined as an application of quantum communication channels as networks (Wolfgang Dyr et al. , 2. In this approach, the communicating parties can harness the communication technology. The motivation behind such a network is performing communications much more secure than the classical algorithms which in principle, can be The parties needing such secure communication can be potentially, military users, the banks or their customers. The fragility of quantum correlations becomes an advantage, as an eavesdropping intervention will disrupt quantum communication. Therefore the security of such a communication is under the guarantee of physics laws. The current state-of-the-art allows for 96 km long quantum communication through optical fibers (Wengerowsky et al. , 2. the experiment, the Authors reported a successful distribution of entangled photon pairs travelling between Malta and Sicily over submarine optical fiber. In this experimental work, the authors show that secure quantum communication is possible over fiber links based on polarization-entangled flying By this result, it is shown that key distribution by quantum protocols could be possible by using polarization entangled Undoubtedly. O 100 km range secure communication is not sufficient for practical The secure communication distance can be enlarged by quantum repeaters in which the protocols are substantially different from the classical repeaters. The distribution of entanglement with the quantum repeater was introduced by Briegel DOI: https://doi. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 | 185 | Indonesian Journal of Science & Technology. Volume 6 Issue 1. April 2021 Hal 183-192 et al. (Briegel et al. , 1. in order to reach longer distances and currently applied as quantum interfaces between flying and stationary qubits (Hammerer et al. , 2. Although a comprehensive experiment covering hundreds of km has not yet been conducted, it has been theoretically calculated that a quantum communication of up to 1000 km can be achieved when a distance of 40 km is placed between the nodes of the quantum repeater (Munro et , 2. Finally, the most challenging task to establish quantum communication between the ground and a satellite was achieved in 2015 (Vallone, et al. , 2. ItAos shown that by using a free space channel, longer distance quantum communication is possible more efficient than optical fibers. Moreover, a Chinese group achieved a quantum communication between 1203 km distance ground stations using a satellite-based entanglement distribution (Yin et al. , 2. METHODS As pointed in the previous section the temporal evolution of entanglement studied with a very simple quantum model composed of two entangled qubits in which one of them coupled to a bosonic reservoir. Figure 1 illustrates the system with two In both schemes, the qubits are initially entangled, denoted as AoAo between the qubits. But in the second scheme, . ottom of the left panel of Figure . , the qubits have a quantum interaction, denoted as AoIiAo, alongside the initial quantum The system plus reservoir could be represented by a Hamiltonian. yaycIycI = 2 . uayc1 yuayc2 ) ycaA yca yaya where yuaycycn are Pauli-Z operators of the qubits, yca A and yca are respectively, the bosonic creation and annihilation operators, and Ea and are the characteristic system and reservoir frequencies. The first two terms of the equation belong to the system and the reservoir free terms and yaya denotes the The Hamiltonian yaya could be defined as yaya = ya. ua 1yuaOe2 Ea. ) yci. ua 2 yuaOe2 ). caA yc. where yuaOe are the Pauli raising, lowering Here, the first term of Equation 2 describes the interaction between the system qubits with an interaction strength ya and the second term describes the coupling between the second qubit and the bosonic reservoir with a coupling strength yci. These interaction and coupling strengths will be used to determine different schemes of the study as control parameters. For instance, when ya = 0, the system describes a case in which there is no interaction between qubits . the other hand, when ya O 0, the system evolution describes the entanglement decay of the interacting qubits. Likewise, the coupling of one of the qubits to the bosonic reservoir will also be evaluated under two different parameter regimes. To this end, the second term of Equation 2 could be rewritten yayaA = yci. ua 2 yuaOe2 ). caA yc. = yci. ua 2 yca yuaOe2 ycaA ) yciA. ua 2 ycaA yuaOe2 yc. Here, when yci = yciA, one recovers Equation 2. This is called the Rabi model for an interaction between a two-level quantum system and a bosonic mode. On the other hand, one could take yciA = 0 as an approximation provided that the coupling is much smaller yci O than the environment This approximation is called the rotating wave approximation (RWA) and the model with RWA is called the Jaynes-Cummings model. Therefore, in this study, one could switch from the Rabi model to the Jaynes-Cummings model simply by taking yciA = 0 in the numerical simulations. Tyrkpenye. Decay of Entanglement of Correlated Qubits Through Bosonic Field. 186 Figure 1. Two schemes of the entangled qubit system . eft pane. and the dynamics of the interaction of the second qubit to the bosonic reservoir . ight pane. Open system dynamics simulated by adopting a model based on repeated interactions process. This model is known as a collisional model (Bruneau et al. , 2. and illustrated in the right panel of Figure 1. According to this approach, the reservoir is represented by a series of identical and noninteracting reservoir units represented by a density matrix yuoycI . The subsystem coupled to the reservoir . ere, the second qubi. interacts with the reservoir unit in a finite time yua. This interaction process is unitary and after time yua the reservoir unit discards and the subsystem and interacts with the forthcoming unit. Each interaction process could be defined by a quantum dynamical 0 ) A ycIycI . = ycOycIycI . uoycIycI ycOycIycI where yuoycIycI = yuoycI . O yuoycI is the two qubit system plus reservoir density matrix and ycOycIycI = expA[OeycnyaycIycI yu. propagator covering both the system and the reservoir unit Hamiltonian yaycIycI . After nycEa interaction . r collisio. the quantum state of the system could be obtained by tracing out the reservoir units in each step as yuoycIycu = Trycu . cOycIycIycu A Tr1. cOycIycI1 . uoycI . yuoycI1 )ycOycIycI ] O A A O yuoycIycu ycOycIycI ycu Here. Trycn is the partial trace operation over the iycEa reservoir unit. The described process above, yields a standard collisional model in which the open system dynamics is Markovian. In this study, the spontaneous emission of the qubits neglected to focus on the effects of the environment type, coupling type to the environment and the initialization of the qubits to the entanglement decay As pointed in the previous section, the objective of this manuscript is to analyse the entanglement dynamics of two entangled qubits under different system and environment parameters. Concurrence (Hill and Wootters, 1. was adopted to quantify entanglement. According to this quantification the concurrence yaycuyca is defined as yaycuyca . uU] = maxA. , yuI1 Oe yuI2 Oe yuI3 Oe yuI4 } . where yuU is the density operator of a two qubit quantum system and . uIycn } are the square roots of the eigenvalues belonging to yuUyuUE where yuUE = . uayc O yuayc )yuU O. uayc O yuayc ). Here, yuU O is the complex conjugation of the two-qubit density matrix yuU. RESULTS AND DISCUSSION In this section the numerical results The time-dependent DOI: https://doi. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 | 187 | Indonesian Journal of Science & Technology. Volume 6 Issue 1. April 2021 Hal 183-192 entanglement evolution is represented in terms of concurrence against the number of collisions . As defined earlier, the system is composed of the initially entangled qubits in which one of them is coupled to a bosonic Two types of initializations is taken into account for the entangled qubits as | and | where | = [|IcIc |IeI. and | = Oo2 [|IcIe |IeI. where Oo |Ic and |Ie denote the orthogonal basis states. On the other hand, the dynamics will be examined against two different reservoir states. yuoycI = . temperature-dependent thermal Gibbs state defined by yuoycI = yce OeyuyaycI ycs where ycs = Tr. ce OeyuyaycI ] is the partition function, yu is the inverse temperature of the reservoir defined by the Hamiltonian yuoycI = AycaA yca. Temperature dependence of the reservoir will be defined by the excitation Ea ycuI = . co ycN) Oe . Oe1 where ycoyaA is the Boltzmann constant. Ea is the reduced Planck constant and ycN is the is the reservoir temperature. The fundamental physical constants are taken to be Ea = 1, ycoyaA = 1 The coupling rates are yci = A0. 1, 0. 15, 0. for all subplots . On the other hand in . , . yciA = 0 and in . , . yci = yciA O 0. The reservoir units are in the yuoycI = . vacuum state for . and in the Gibbs state for . In each collision yua = 0. 4 and there is no interaction between the qubits . a = . for all . - . The first numerical results reported As mentioned above, in all simulations, the qubits are maximally entangled initially and one of the qubits coupled to a bosonic reservoir. In Figure 2, the reservoir is in a vacuum state that is the reservoir is a zero-temperature reservoir. One observes an exponential decay in the presence of a vacuum reservoir in Figure 2 . corresponding to . ci O 0, yciA = . a JaynesCummings type emitter-bosonic mode It's typical to observe, numerically, an exponential decay in the presence of a zero temperature bath neglecting the spontaneous emission. However, when the system-reservoir interaction is assumed to be strong enough to use the Rabi model . ciA O . one observes the entanglement decay at finite time in Figure 2 . This phenomenon is called Aoentanglement sudden death' (Yu and Eberly, 2. and in general, occurs in the presence of thermal Therefore, we obtain the first notable result tells us that the Rabi type coupling is more detrimental for entanglement even in the presence of vacuum noise. On the other hand, in Figures 2 . - . , the entangled qubits placed in a thermal reservoir and the dynamics compared for yciA = 0 and yci = yciA O As obvious in Figures 2 . - . , entanglement vanishes at finite time. But for the Rabi type coupling . ci = yciA O . corresponding to Figure 2 . , entanglement lifetime is much smaller. Therefore. Rabi type coupling to the reservoir reveals its detrimental nature in the presence of the thermal noise stronger. Here, one should also note that there is no interaction . a = . between the qubits during the evolution and the initialization of entanglement has no effect on the dynamics depicted in Figures 2 . - . The initial entangled state are | for . , . , . and | for . , . , . The coupling rates are g = 0. 1, 0. 15, 0. 2 for all subplots . - . On the other hand, in . - . AyciA = 0 and in . , . yci = yciA O 0 Moreover, the qubits are interacting and ya = 0. 05 for . The reservoir units are in the yuoycI = . vacuum state for . and in the Gibbs state for . , . , . Thermal Tyrkpenye. Decay of Entanglement of Correlated Qubits Through Bosonic Fields | 188 excitation number is ycuI = 0. - . In each collision yua = 0. Next, there is an Aoalways onAo interaction between the entangled qubits . a O . in the calculations in Figure 3. Comparing Figures 3 . - . , entangled qubits placed in a vacuum reservoir for Jaynes-Cummings type . ciA = . More specifically, the qubit initial states are | for . and | for . As obvious in the plots Figures 3 . , the entanglement lifetime has been enlarged That is, as the second significant result in this study, initial preparation of entanglement could affect the entanglement dynamics only when there exists an interaction . a O . between the This result applied to the case where the reservoir is thermal and the coupling to the reservoir is yciA = 0 in Figures 3 . - . Figure 2. Decay dynamics of entanglement for a pair of entangled qubits in the presence of a bosonic reservoir in different states depending on different coupling rates yci against number of collisions . DOI: https://doi. org/10. 17509/ijost. p- ISSN 2528-1410 e- ISSN 2527-8045 | 189 | Indonesian Journal of Science & Technology. Volume 6 Issue 1. April 2021 Hal 183-192 Figure 3. Decay dynamics of entanglement for a pair of entangled qubits in the presence of a bosonic reservoir in different states depending on different coupling rates yci against number of collisions . As a typical result for a thermal reservoir, one observes entanglement sudden death in these plots. But it is clear that even under these circumstances a change in the entanglement initialization could enlarge the entanglement lifetime as seen in Figures 3 . Finally, the effect of the change of entanglement initialization for the Rabi-type interaction in the presence of a thermal reservoir for ya O 0 was examined in Figures 3 . - . As is clear in the plots, entanglement initialization has no effect for the Rabi-type interaction between the system and the reservoir. DISCUSSION The obtained results could serve for different tasks of quantum information. For instance, the obtained results about the coupling type of the entangled qubits to the reservoir could be useful for the flying qubits (Guo et al. , 2. carrying quantum On the other hand, the second significant result of the current study asserts Tyrkpenye. Decay of Entanglement of Correlated Qubits Through Bosonic Fields | 190 that the entanglement lifetime could last longer for the interacting qubits may not be suitable for flying qubits as the interaction of flying qubits is still a challenge. However, benefiting the obtained results of the current study for quantum networks (Reiserer and Rempe, 2. composed of stationary qubits would be another scheme. By the recent developments on artificial intelligence and machine learning, quantum versions of networks became crucial to the field (Levine et al. , 2019. Tyrkpenye et al,. Tyrkpenye, 2. Therefore, the obtained results could be used for both quantum communication and quantum network architectures. The suggestions can be applied in the stage of the preparation of the entangled photon pairs. The results of the article are based on the assumption that the interaction of flying qubits with the environment can be The control of inter-qubit interaction was also assumed. That can be achieved through a dynamical control method referred to as 'dynamical decoupling' as referenced in the introduction This method was applied entanglement in a quantum network (Humphreys et al. , 2. By these methods, as suggested in the current paper, letting or avoiding the qubit-qubit interaction in the creation stage of the flying qubits can extend the lifetime of quantum correlations. Moreover, instead of flying qubits, the superconducting circuits where quantum networks can be built upon them. As a recent study reports (Campagne-Ibarcq et al. , 2. the entanglement distribution between the non-interacting nodes of a quantum network becomes possible by the application of entangling gates between the stationary and flying qubits acting as quantum buses. In this architecture, the dynamical decoupling methods enabling to switch on or off the inter-qubit or qubit-reservoir coupling could be applied properly. CONCLUSION This study deals with a simple quantum system representing an entangled and open two-level quantum systems with different system and reservoir parameters. It's found that the quantum interactions could be beneficial for the entangled systems for specific initialization schemes. On the other hand, for the non-interacting quantum systems, the coupling type to the reservoir should be controlled as different types of couplings effect the entanglement lifetime More particularly. Rabi-type couplings to the bosonic reservoirs are found to be detrimental for the entanglement lifetime where the performance could be improved by letting the quantum The methods and the results could apply to any physical system covering the non-classical quantum resources beneficial for quantum technologies. AUTHORSAo NOTE The authors declare that there is no conflict of interest regarding the publication of this article. Authors confirmed that the paper was free of plagiarism. REFERENCES