Quantum speedup, non-Markovianity and formation of bound state.

In this paper, we investigate the relationship between the quantum speedup, non-Markovianity and formation of a system-environment bound state. Previous results show a monotonic relation between these three such that providing bound states with more negative energy can lead to a higher degree of non-Markovianity, and hence to a greater speed of quantum evolution. By studying dynamics of a dissipative two-level system or a V-type three-level system, when similar and additional systems are present, we reveal that the quantum speedup is exclusively related to the formation of the system-environment bound state, while the non-Markovian effect of the system dynamics is neither necessary nor sufficient to speed up the quantum evolution. On the other hand, it is shown that only the formation of the system-environment bound state plays a decisive role in the acceleration of the quantum evolution.

Introduction

The quantum speed limit (QSL) time[1,2] is defined as the minimal time a quantum system needs to evolve from an initial state to a target state, which can be used to characterize the maximal speed of evolution of a quantum system. It arises in tremendous areas of quantum physics and quantum information, such as nonequilibrium thermodynamics[3], quantum metrology[4], quantum optimal control[5], quantum computation[6,7], quantum communication[3,8], and has been studied extensively due to its fundamental importance. An important question is how to exploit the available resources to achieve the highest speed of quantum evolution where the problem of determining it is relevant to deriving physical limits in a variety of contexts. In closed systems, entanglement has been taken as a resource in the speedup of quantum evolution[9-11]. Since any real system is coupled to its environment in the practical scenarios, a considerable effort has been devoted to study controlling speedup in open quantum systems[12-14]. In this regard, Deffner and Lutz showed that the dynamical acceleration in the damped Jaynes-Cummings model for a two-level system can be induced by the non-Markovian effect, and therefore results in a smaller QSL time[14]. This phenomenon has been verified in different settings[15-18] and confirmed by the experiment in cavity quantum electrodynamics (QED) systems[19], where an atomic beam is treated as a controllable environment for the cavity field system. Recently, the authors of ref. [20] showed that both the non-Markovianity and the quantum speedup are attributed to the formation of a special eigenstate, known as a bound state, which has negative eigenenergy and is separated from the other ones by an energy gap[21-24]. It was demonstrated that if the bound state in a two-level system coupled to an environment with Ohmic spectrum is established, the evolution of the system becomes non-Markovian, and this results in quantum speedup. They found a monotonic relation between them such that providing bound states with more negative energy can result in a higher degree of non-Markovianity, and this also leads to bigger speed of quantum evolution. Despite the growing body of literature on this subject, the analysis has almost exclusively focused on the simpler settings and to our knowledge, the mechanism of speedup has not been investigated in the more complex settings.

The aim of the present study is to investigate the relationship between the quantum speedup, non-Markovianity and formation of a system-environment bound state in the more complex settings. To this aim, we study the mechanism of quantum speedup for a two-level system and a V-type three-level system in the presence of similar and additional systems. We illustrated that there is a monotonic relation between the quantum speedup and the formation of a system-environment bound state. In particular, we show that the bound states with more negative energy provided by the additional systems can speed up the quantum evolution. However, it is figure out that there is not such monotonic relation between the quantum speedup and the non-Markovianity. We demonstrate that although non-Markovianity governs the quantum speedup in the absence of additional systems, it is neither necessary nor sufficient to speed up the quantum evolution when the environment becomes more complex in the presence of additional systems. On the other hand, our results suggest that only the formation of the system-environment bound state is the essential reason for the quantum speedup.

Results

The condition for the existence of a bound state in a system including N atoms

This section is devoted to a system containing N non-interacting atoms coupled to a reservoir, as depicted in Fig. 1. As will be discussed in the following cases, we check the energy distribution of the Hamiltonian of the mentioned system by considering the atoms as two-level and three-level systems. Then, it will be shown that how the system-environment bound state, as a special eigenstate with eigenvalue residing in the band gap of the energy spectrum, can be formed in these models.

Figure 1

The figure corresponds to seven atoms (N = 7) immersed in a common reservoir. The central atom (the orange circle) is considered as our main concern of the single-atom system and the remainder ones are considered as the additional atoms (the blue circles).

Case I:

Two-level systems. By considering the atoms as two-level systems, the Hamiltonian of the whole system can be written as ()where the first term is the Hamiltonian of N two-level systems, the second term describes the reservoir Hamiltonian and the last term is the system-reservoir interaction Hamiltonian. Also, stands for the raising (lowering) operator of the lth atom with transition frequency . () creates (annihilates) the kth field mode with frequency and the coupling strength between the lth atom and the kth field mode is represented by .

The following eigenvalue equation gives the energy spectrum of the total Hamiltonian

Considering the single excitation subspace, the commutation of Hamiltonian with the number operator of excitations (i.e., ) results in the state of the whole system as followswhere means that except for the lth atom which is excited, all other atoms are in their corresponding ground states. and hold for ground and excited states, respectively, and is defined as . The vacuum state of the reservoir is denoted by , and is the kth field mode with just one excitation. Now we substitute Eqs (1) and (3) into Eq. (2) and obtain the following set of equations

If one obtain from the first equation and substitute it into the rest of the equations, the following N integro-differential equations are obtained

After summing up the Eq. (5) and eliminating from both sides, we can obtainwithwhere depending on the functionality of the spectral density of the reservoir and the number of atoms , we will have different results. In a quantum many-body system, the bound state is defined as an eigenstate with real (discrete) eigenvalue. So, it can be claimed that the system possesses a bound state, when Eq. (6) has a real root. It can be easily seen that decreases monotonically with the increase of in the regime of . As a result, when the condition is satisfied, the function on the right-hand side of Eq. (6) always has one and only one intersection with . This root just corresponds to the eigenvalue of the formed bound state in the Hilbert space of the whole system (N two-level system plus its common reservoir). Moreover, we can observe that is divergent in the positive energy range, i.e. , and thus Eq. (7) does not have a real root to make it well defined. Therefore, no real root can be found in this range to support the existence of a further bound state.

Case II:

Three-level systems. After studying the energy spectrum of the system including N two-level systems in a common reservoir, we investigate this problem by considering three-level systems in V configuration. By considering two excited states for each atom, on the other hand, each of them spontaneously decays into the ground state in such a way that the respective dipole moments of transitions may have interaction with each other through spontaneously generated interference (SGI). Also, for a given V-type atom, the decaying frequencies of excited states and into the ground state are denoted by and , respectively. The Hamiltonian for this described system can be written as

The raising and lowering of the mth excited state of the lth atom to the ground state are carried out by . Also, () annihilates (creates) the kth field mode with the frequency . Furthermore, it is assumed that , as the coupling strength between the mth excited state and the kth field mode, is identical for all atoms.

Now, we want to study the energy spectrum of N three-level systems in a common reservoir. For this purpose, the following eigenvalue equation is solvedwhereis the state of the whole system in the single excitation subspace, because the total Hamiltonian conserves the number of excitations in the system, i.e. . Also, means that all atoms are in the ground states. Moreover, and means that except for the lth atom, all other atoms are in the ground state. It is also considered that and are respectively in the first and second upper excited levels. We denote is the vacuum state of the reservoir and the state of it with only one excitation in the kth field mode. For simplicity, we assume the case of the degenerate excited levels where . Then, by substituting Eqs (8) and (10) into Eq. (9), we obtain the following differential equations

We derive from the first equation and substitute it into the rest ones, so we have

By taking a Lorentzian spectral density as given in Eq. (40) and under the assumption considered in Methods, we have where is the angle between two dipole moments of the two possible transitions in each atom (see Methods). Thus, by summing up the above equations in Eq. (12) and eliminating from both sides, a compact form of Eq. (12) can be obtained aswhere

It is clear that the functionality of the spectral density of the reservoir , the number of atoms , and the SGI parameter are effective factors for the solutions of Eq. (13). Nevertheless, as previously discussed in case I, the bound state is appears if Eq. (13) has at least one real solution for negative energies, i.e., . Otherwise, the bound state does not formed.

QSL time and non-Markovianity

Here, we aim to evaluate QSL time bound and non-Markovianity of a single atom (two-level or three-level) in the presence of () additional atoms where all of the atoms are contained in a common reservoir. In this regard, for discussed cases I and II, the first atom (l = 1) is taken as the main concern of the our single system and () remainder atoms are taken as additional atoms.

According to ref.[14], a unified expression for the QSL time is provided in open quantum systems aswhere , is the Bures angle between initial pure state and the target state governed by the time-dependent non-unitary equation , andwith , which denotes the Schatten p-norm and are the singular values of X. The intrinsic speed of a dynamical evolution is likely to evaluate by applying the recent achieved bound for the QSL time and a given actual driving time . For the case of , there is no potential capacity for further speedup and speedup does not appear. But, when we have , it indicates the potential capacity for quantum dynamical speedup.

We first consider the exactly solvable model for a two-level system in the presence of (N − 1) additional atoms which are immersed in a common reservoir (see Methods). For simplicity and without loss of generality, we suppose that our considered system (the 1th atom) is initially in the excited state and (N − 1) additional atoms are in the ground state and the environment is in a vacuum state[14]. Thus, the evolved state at time t can be obtained as , where the explicit form of is calculated by taking a Lorentzian spectral density for the reservoir in the Methods section. So, the QSL time bound for the considered atom can be evaluated as[17]

Moreover, according to ref.[25] the non-Markovian dynamics of the considered two-level atom can be studied through the notion of the reverse flow of information from the reservoir back to the system. Here, we apply this measure of non-Markovianity as followswhereand implies the trace distance of a pair of states. The trace norm of operator is defined as . Hence, stands for the changing rate of the trace distance . For our two-level system, by choosing the optimal pair of initial states obtained in ref.[26], it is easy to check that the trace distance of the evolved states can be written as . Then, one can easily verify that[17]which connects to Eq. (17) as

It is clear that, the QSL time is related to the non-Markovianity within the driving time and the atomic excited population .

Next, the QSL time bound for a V-type three-level atom in the presence of () atoms within a common reservoir will be obtained. To this aim, we consider the exactly solvable model for a V-type three-level atom in the presence of () additional atoms which are immersed in a common reservoir (see Methods) and assume the initial condition as and for , which implies that . Under these assumptions, it can be easily checked that , so the following inequality holds . Also, it is clear that . Consequently, can be evaluated as

Moreover, the measure of non-Markovianty for the considered V-type three-level atom can be obtained by considering the optimized pair of initial states proposed by ref.[27]. In this regard, we use the definition of trace distance, , where , and are the eigenvalues of . Thus, the measure of non-Markovianty in Eq. (18) is reduced to

Here, in contrast to the previous case, we emphasize that no explicit relationship can be established between the non-Markovianity and the QSL time bound of the considered V-type three-level atom.

Discussion

In the following, we illustrate the influence of non-Markovianity and formation of a bound state on the QSL time for our considered model in the last section. In Fig. 2, we present the ratio with the behaviors of non-Markovianity of a two-level atom in terms of the coupling strength and for different numbers of additional atoms in the reservoir. Here, the ratio shows the potential capacity for quantum dynamical speedup. As we can see in Fig. 2, in the critical point from no-speedup to speedup of quantum evolution, the Markovian environment becomes non-Markovian. Also, it is surprisingly seen that although increasing the number of additional atoms leads to more decrement of the QSL time, the non-Markovianity of the system does not always increase by inserting the additional atoms into the reservoir. This means that, non-Markovianity which is induced by the memory effect of environment cannot always induce dynamical acceleration in our considered model. To justify this, according to Eq. (21), it should be noted that when and provided that , no speedup can be observed in the quantum evolution, i.e. . Moreover, in the absence of additional atoms (), and by choosing the actual driving time large enough, no population will be trapped in the excited state of the considered two-level system, i.e. . Therefore, the QSL time will be inversely proportional to the non-Markovianity as , and as a result, the unique reason for speeding up of quantum evolution is the non-Markovianity and it induces dynamical acceleration which is matched with results in ref.[14]. However, by inserting the additional atoms into the reservoir (), the excited-state population approaches to a non-zero steady value in a long enough time, i.e. (see Eq. (32)). Therefore, the competition between the population, , and non-Markovianity, , takes responsibility for the intrinsic speedup of quantum evolution. On the other hand, we can claim that the main reason for the speedup is not solely due to the non-Markovianity. One question naturally arises about an essential reflection to the quantum speedup. To answer this question, we further investigate the speedup from the perspective of formation of a system-environment bound state and demonstrate that how the quantum speedup is exclusively related to the formation of the system-environment bound state. The system-environment bound state is the stationary state of the whole system in which the inhibition of spontaneous emission is the result of the formation of it[28,29]. To this aim, the ratio with the negative energy spectrum of the total Hamiltonian in Eq. (1), has been sketched in Fig. 3. It is obvious that the transition from no-speedup to speedup of quantum evolution is occurred faster by increasing the number of additional systems. Also, the critical points for forming the bound state match well with the ones for presenting quantum speedup. It is worth noting that providing bound states with more negative energy in the system-environment spectrum via inserting more additional atoms into the reservoir, can lead to greater speed of quantum evolution. Thus, we can conclude that the essential reason for the acceleration of evolution is related to the formation of a bound state.

Figure 2

The QSL time (blue solid line) and non-Markovianity (red dashed line) for a two level system in terms of the coupling strength and with different numbers of atoms in the reservoir as (a) N = 1, (b) N = 3, (c) N = 8 and (d) N = 30. The values of the used parameters are (in units of ) and (in units of ).

Figure 3

The QSL time (blue solid line) for a two level system and the energy of the formed bound state (green dashed line) in terms of the coupling strength and with different numbers of atoms in the reservoir as (a) N = 1, (b) N = 3, (c) N = 8 and (d) N = 30. The values of the used parameters are (in units of ) and (in units of ).

For more investigation of this interesting finding, we proceed further by considering a V-type three-level atom in the presence of similar and additional atoms. In this regard, the ratio with the non-Markovianity of the V-type three-level atom has been presented in Fig. 4 and with the negative energy spectrum of the total Hamiltonian in Eq. (8), has been presented in Fig. 5. In addition to the previous observations that are also seen for the model considered here, we find that the speedup phenomenon is occurred faster by increasing the SGI parameter. Also, it is clear that increasing the SGI parameter leads to improvement of degree of boundedness, which in turns causes the speedup of the evolution. Consequently, these results reveal that the appearance of system-reservoir bound state with a higher degree of boundedness ensures the enhancement in the intrinsic speed of evolution irrespective of non-Markovianity of the system.

Figure 4

The QSL time (blue solid line) and non-Markovianity (red dashed line) for a V-type three level system in terms of the coupling strength and with different numbers of atoms in the reservoir as (a) N = 1, (b) N = 3, (c) N = 8 and (d) N = 30. The values of the used parameters are (in units of ) and (in units of ). Also, lines with the square marks are plotted for θ = 1 and without the marks are plotted for θ = 0.

Figure 5

The QSL time (blue solid line) and non-Markovianity (red dashed line) for a V-type three level system in terms of the coupling strength and with different numbers of atoms in the reservoir as (a) N = 1, (b) N = 3, (c) N = 8 and (d) N = 30. The values of the used parameters are (in units of ) and (in units of ). Also, lines with the square marks are plotted for θ = 1 and without the marks are plotted for θ = 0.

Conclusion

In conclusion, the mechanism of quantum speedup for a two-level system and a V-type three-level system has been studied in the presence of additional systems from two perspectives: non-Markovianity and formation of a system-environment bound state. We have demonstrated that although non-Markovianity governs the quantum speedup in the absence of additional systems, it is neither necessary nor sufficient to speed up the quantum evolution when the environment becomes complex in the presence of additional systems. Also, it has been revealed that the bound state of the whole system plays a deterministic role in quantum speedup of the considered systems. In other words, providing bound states with more negative energy can lead to greater speed of quantum evolution. Our results may open new perspectives for acceleration of evolution through engineering the formation of a bound state.

Methods

Dynamics of a N two-level system immersed in a common reservoir

Here, we want to derive the dynamics of an open two-level atomic system in the presence of () additional systems which are immersed in a common reservoir[30]. In this regard, we obtain the exact dynamics of a system consisting of non-interacting two-level atoms in a common reservoir as considered in case I. For simplicity, we work in the interaction picturewhere the Hamiltonian in this picture is given byandwith and . Substituting Eqs (25) and (26) into Eq. (24) gives the following differential equations aswith . It is clear that , so . Then, integrating the last equation of Eq. (27) and substituting it into the others giveswhere the integration of with respect to leads to the correlation function as

It is common to use the Lorentzian function as the spectral density of the reservoir. So, we take it as followswhere is the central frequency of the reservoir, the parameter λ defines the spectral width of the coupling and γ0 is the coupling strength. Therefore, by considering the spectral density given by Eq. (30) and using the Laplace transform, the exact solutions of the probability amplitudes can be obtained as[30]where . Also, by considering the 1th atom as our main concern of system () and assuming the initial condition as for , Eq. (31) is reduced to

Moreover, tracing over the reservoir and the () additional atoms results in an explicit form for the reduced density operator of the 1st atom in the basis of as follows

Dynamics of N V-type three-level system immersed in a common reservoir

In the following, we derive the dynamics of an open V-type three-level atomic system in the presence of () additional systems[31]. To this aim, we obtain the exact dynamics of a system consisting of non-interacting V-type three-level atoms in a common reservoir as considered in case II. According to Schrödinger equation in the interaction picture, we havewhereandwith , and . It can be easily seen that , then , while by substituting Eqs (35) and (36) into Eq. (34), we have the following differential equations

Since it is important for us to obtain the coefficients , so we integrate the last differential equation in (37) and substitute it in the rest of differential equations. As a result, we obtain the following set of closed integro-differential equations

It should be noted that the kernels in Eq. (38) can be expressed in terms of spectral density as followsand is chosen as Lorentzian distributionwhere and λ are respectively the center frequency of the structured reservoir and the spectral width of the coupling. Also, show the relaxation rates of the two upper excited levels. And, , provided that , stand for the SGI between existing transitions in each atom. It is worth to note that the parameter depends on the relative angle between two dipole moment elements of the mentioned transitions. While , indicating no SGI between the two transitions, reflects the perpendicularity of the dipole moments of transitions to each other, on the other hand, reflects that the two dipole moments are parallel, verifying the strongest SGI between the transitions. The Laplace transform of Eq. (38) results in

Here, , is the Laplace transform of . According to Eq. (41), it is obvious that the following relation between the coefficients is found

Now, writing the coefficients () in terms of and inserting them into the Eq. (41), leads to the following equation

Here, we turn our attention to a special case in which the degeneration of two upper atomic states, and the resonance between the atomic transitions and the central frequency of the reservoir are considered, i.e., . Under this consideration, we can assume that and , so the kernels in Eq. (39) takes the following simple formwhere . In the following, we define the new coefficients , and rewrite Eq. (43) in terms of them astherefore, after taking inverse Laplace transform we getwhere and . Finally, it should be noted that we take the first atom () of the studied system as our main concern of the single three-level system and the () remaining atoms as the additional atoms. Hence, tracing over the reservoir and the () additional atoms leads to the explicit form of the reduced density operator of the first atom in the basis of as follows