Asymptotic entanglement sudden death in two atoms with dipole-dipole and Ising interactions coupled to a radiation field at non-zero detuning

We investigate the time evolution and asymptotic behavior of a system of two two-level atoms (qubits) interacting off-resonance with a single mode radiation field. The two atoms are coupled to each other through dipole-dipole as well as Ising interactions. An exact analytic solution for the system dynamics that spans the entire phase space is provided. We focus on initial states that cause the system to evolve to entanglement sudden death (ESD) between the two atoms. We find that combining the Ising and dipole-dipole interactions is very powerful in controlling the entanglement dynamics and ESD compared with either one of them separately. Their effects on eliminating ESD may add up constructively or destructively depending on the type of Ising interaction (Ferromagnetic or anti-Ferromagnetic), the detuning parameter value, and the initial state of the system. The asymptotic behavior of the ESD is found to depend substantially on the initial state of the system, where ESD can be entirely eliminated by tuning the system parameters except in the case of an initial correlated Bell state. Interestingly, the entanglement, atomic population and quantum correlation between the two atoms and the field synchronize and reach asymptotically quasi-steady dynamic states. Each one of them ends up as a continuous irregular oscillation, where the collapse periods vanish, with a limited amplitude and an approximately constant mean value that depend on the initial state and the system parameters choice. This indicates an asymptotic continuous exchange of energy (and strong quantum correlation) between the atoms and the field takes place, accompanied by diminished ESD for these chosen setups of the system. This system can be realized in spin states of quantum dots or Rydberg atoms in optical cavities, and superconducting or hybrid qubits in linear resonators.


I. INTRODUCTION
Studying the quantum phenomena in systems of atoms coupled to radiation fields has been in the center of interest in physics since the Jaynes-Cummings model was introduced in 1963 [1]. This interest was boosted by the development of new quantum structures that are considered very promising for serving as the building units in quantum information processing (QIP) systems [2], while at the same time are highly interacting with and controllable by radiation fields, which include artificial atomic systems such as semiconducting quantum dots and superconducting circuits in addition to the customized atomic systems such as Rydberg atoms and trapped atoms, ions and molecules [3][4][5][6]. Interaction between natural regular atoms, in cavity quantum electrodynamics (CQED), used to be ignored as a result of their considerably small magnitude. However, the newly developed quantum systems are characterized by strong interaction with the same type of system or even other types in hybrid quantum composite systems. Embedding superconducting qubits in a superconducting microwave resonator was a huge step toward utilizing these systems in QIP [7][8][9] and establishing the new field of circuit quantum electrodynamics (cQED), where close and distant superconducting qubits can be coupled through local interactions or microwave photons [10][11][12][13][14]. Also, other systems can be embedded in superconducting resonators, through hybrid circuit structures, such as spins in quantum dots and solid state impurities [4,6]. Similar arrangements were implemented for coupled Rydberg atoms in CQED [15][16][17][18]. Laser-trapped circular Rydberg atoms were used for analogue quantum simulation of spin arrays [19], where the strong interaction between the atoms is used to simulate an XXZ spin chain Hamiltonian, which amounts for dipole-dipole and Ising interactions between the spins. It was suggested that such a scheme can be also utilized in hybrid structures in cQED, where Rydberg atoms can be integrated into superconducting circuits.
Spin-spin interaction simulation and modeling in other systems such as optical lattices [20], trapped ions [21], coupled molecules [22,23] and microcavities [24] were performed. These important developments led to a growing interest in studying decoherence and entanglement dynamics, and population inversion in systems of atoms (qubits), with either dipole-dipole or Ising interaction, or both, in the presence of radiation fields, keeping an eye on their QIP implementations. Coupled spin (qubit) systems in absence of radiation fields but in presence of magnetic fields and different types of environments have been studied intensively [25][26][27][28][29][30][31][32][33][34][35][36].
One of the main obstacles toward realizing reliable quantum computing systems in general, and particularly using these newly customized quantum systems, is entanglement sudden death (ESD). ESD is observed when the entanglement loss in the system takes place very rapidly leading to a disentangled state. It was first identified, discussed and named by Yu and Eberly [37,38], when they studied the entanglement between two uncoupled, but initially entangled, atoms in two separate cavities. They also showed that the same phenomenon may take place under the effect of a noisy classical environment on two uncoupled atoms [39]. Several other works have studied ESD in systems of two non-interacting atoms in remote cavities as well [40][41][42][43]. The effect of an out of resonance radiation field on systems of identical and non-identical, non-interacting, atoms was studied, where it was shown that the non-zero detuning can be an advantage for preserving entanglement [44,45].
The Bell inequality was tested using a system of two uncoupled qubits interacting with a radiation field in an optical cavity [46]. Although systems of interacting atoms have been considered as well, but they were only at resonance with the field to avoid the mathematical difficulty caused by the off-resonance condition [47][48][49][50]. Particularly, ESD was studied in a system of two coupled identical atoms interacting at resonance with a double mode radiation field, where the effects of the coupling as well as the initial state of the system on the system dynamics were investigated [51]. Entanglement and purity in a system of two interaction atoms coupled to a radiation field at resonance were investigated [52], where the effect of the interplay between the atom-atom and the atom-field field couplings on the system was studied thoroughly. A system of two coupled atoms interacting off-resonance with a radiation field was studied [16], where the system was represented by XXZ model, considering dipole-dipole and Ising interactions at the same time. They provided an analytic perturbative solution for the system dynamics assuming a weak atom-atom coupling. They proposed an experimental realization for the system using spins states in quantum dots in CQED and superconducting qubits in cQED. Another work considered two identical atoms that are coupled to each other through dipole-dipole and Ising interactions while coupled to a radiation field at resonance, where an analytic solution was provided [53]. To study the entanglement in the system numerical calculations were performed. Recently, a system of two coupled qubits interacting with a common environment was investigated, where schemes to avoid ESD using Local unitary operations were provided [54]. Very recently, the entanglement dynamics of a pair of well-separated Rydberg pairs driven by a common 4 laser field while interacting via both intra-pair and inter-pair van der Waals potentials was investigated [55]. It showed in-phase (anti-phase) beating dynamics that depends on the inter-pair potentials and the field detuning.
As can be noticed, in the previous works, the dynamics of the system of atoms coupled to a radiation field were studied intensively but under certain restrictions, due to the mathematical formidability of the problem. It was considered either uncoupled atoms interacting off-resonance with the field or coupled atoms (with a dipole, Ising, or both interactions) at resonance with the field. In a recent work [56], we studied a system of two coupled atoms (qubits) with dipole-dipole interaction in the presence of an off-resonance radiation field.
We presented an analytic solution for the time evolution of the system and showed how the combined effect of the dipole coupling and the non-zero detuning can be utilized to control the ESD in the system.
In this paper, we study a system of two coupled atoms interacting off-resonance with a single-mode radiation field. We provide an exact analytic study of the system dynamics while considering simultaneously the effect of the two types of interactions, dipole-dipole and Ising, which can be modeled as a Heisenberg spin 1/2 XXZ interaction. This system is important for its own sake as a model of two coupled spins interacting with a bosonic bath, as well as for its impact on the field of cavity (circuit) QED and its implementations in QIP. It can be realized in electron spin states in quantum dots or Rydberg atoms in optical cavities as well as in superconducting or hybrid qubits in linear resonators. This work represents a crucial completion of our previous one that is vitally needed to provide a full understanding of the system dynamics and asymptotic behavior. We emphasize, using our results, the significant impact of the Ising interaction, on its own or when combined with the dipole one, on the system dynamics. Particularly, we show how the type of interaction (Ferromagnetic vs. anti-Ferromagnetic) may lead to considerably different effects on the system dynamics, especially the ESD, depending on the initial state and the detuning value. Moreover, we investigate the asymptotic behavior of the system dynamics and particularly the ESD, under the interplay of the different system parameters, which have not been addressed before in the literature or in our previous work. Most of the previous works focused on the early dynamics of the system and the treatment of the ESD at that stage. We show that the system reaches asymptotically a quasi-steady dynamic state, where the ESD can be eliminated by utilizing the interplay between the Ising and the dipole-dipole interactions as well as the detuning parameter for all initial states except the maximally entangled correlated Bell state. Also, we demonstrate how the entanglement, population inversion and quantum correlation between the two atoms ensemble and the radiation field synchronize asymptotically with remarkable profiles, where the collapse periods disappear indicating a continuous exchange of energy (and quantum correlation) between the atoms and the field, accompanied by diminished ESD for chosen setups of the system. This paper is organized as follows. In Sec. 2, we discuss our model and present the exact analytic solution. In sec. 3, we study the dynamics of entanglement and atomic population inversion, starting from different initial states. The quantum correlation between the two atoms ensemble and the radiation field is presented in sec. 4. We conclude in Sec. 5.

II. THE MODEL AND ITS EXACT SOLUTION
We consider a system of two identical atoms (qubits), each one of them is characterized by two levels: ground |g i and excited |e i , where i = 1, 2 corresponding to the first and second atoms respectively. The two atoms are coupled to each other through dipole-dipole and Ising interactions, which are modeled as an XXZ exchange interactions between two spin-1/2 particles with XX coupling strength λ 2 and Z coupling strength J. They are coupled to the same single-mode quantized radiation field with the same coupling constant λ 1 . The Hamiltonian of the system is given bŷ (1) The first and second terms in the Hamiltonian represent the free quantized radiation field and the non-interacting two atoms respectively. The third, fourth, and fifth terms represent the atom-field, dipole-dipole, and Ising interactions respectively. Ω and ω • are the frequencies of the single-mode radiation field and the atomic system transition respectively,â † andâ are creation and annihilation operators of the radiation field which satisfy the usual commutation relation [â,â † ] = 1 andσ (i) ν , where (ν = x, y, z, ±), are the usual Pauli, raising and lowering spin operators representing the ith qubit.
The atoms are assumed to be initially in a pure state and the field is in a coherent state, therefore the initial wave function of the composite system becomes where a, b, c and d, are arbitrary complex quantities that satisfy the condition |a| 2 + |b| 2 + |c| 2 + |d| 2 = 1, and |α is the coherent state defined as where |α| 2 = n is the mean photon number and |n are the photon number states, which satisfy the relations:â † |n = √ n + 1|n + 1 andâ|n + 1 = √ n + 1|n . The wave function at any time t latter can be written as where |e 1 , e 2 , n is the state in which both of the two atoms are in the excited state and the field has n photons, while |e 1 , g 2 , n + 1 is the state in which the first one is in the excited state and the second is in the ground state and the field has n + 1 photons, etc.
Rewriting the Hamiltonian (Eq. (1)) aŝ whereĤ where ∆ = ω • − Ω is the detuning parameter andN =â †â + 1 z , is the total number of excitations in the system, which is a constant of motion. It is more convenient to work in the interaction picture, where we defineV I =ÛĤ intÛ † andÛ = e iĤ•t . This yieldŝ Now substituting |ψ(t) and V I (t) into Schrödinger equation we get the following system of differential equations where α = λ 1 √ n + 1 and β = λ 1 √ n + 2.
Following the same approach that we implemented in our previous work [56], and after some calculations, the solution takes the form where and where and The initial values of the coefficients are as follows Now we can construct the system wave function |ψ(t) and calculate the system density The reduced density matrix of the two atoms,ρ red (t), can be obtained by tracing out the field

SION
Utilizing the derived analytic solution, we can investigate the entanglement dynamics between the two atoms and their population inversion starting from different initial states that lead to ESD. For convenience, we set = 1, λ 1 = 1 and represent the other parameters (λ 2 , J and ∆) in units of λ 1 . The entanglement between the two quantum systems is quantified using the concurrence function C(ρ red ) [57]. The concurrence is related to the h is the Shannon entropy function h(x) = −x log 2 x−(1−x) log 2 (1−x), and the concurrence The ε i arranged in decreasing order are the square root of the four eigenvalues of the non-Hermitian matrix R ≡ ρ redρred , Whereρ red = (σ y ⊗σ y )ρ * red (σ y ⊗σ y ). Both of C(ρ red ) and E f (ρ red ) take values from 0 for a separable (disentangled) state to 1 for a maximally entangled state. On the other hand, atomic population inversion is defined as the difference between the probabilities of finding the atom in its excited state and ground state or simply the expectation value of the operator σ z . Using the reduced density matrix of any one of the two identical atomsρ 1 (t), which can be obtained by tracing out the other one inρ red (Eq. 17), we can evaluate σ z (t) as 9 A. Initial Bell state In Fig. 1, we explore the dynamics of entanglement and population inversion, in terms of the scaled time τ = λ 1 t, starting from an initial correlated Bell state ψ Bc = (|e 1 |e 2 + |g 1 |g 2 )/ √ 2 with the radiation field is in a coherent state with field intensity corresponding ton = 100. In the forthcoming discussion, we setn = 100 everywhere except when otherwise is mentioned explicitly. Starting from such an initial state the system shows ESD, where the entanglement changes abruptly from a non-zero to an exact zero value, as illustrated in the different panels of the figure. In our forthcoming discussion and in the legends of all plots, we will use the pair notation (a, b) to refer to the values of the two parameters λ 2 and J, where "a" refers to the value of λ 2 , while "b" refers to the value of J. In Fig. 1(a), we study the entanglement dynamics at 12 here). The entanglement dynamics for negative detuning, ∆ = −2, illustrated in Fig. 1(c), does not exhibit notable differences from the positive detuning case, ∆ = 2, except the peaks (green line) for (0,2) is slightly higher, this time, than that of (0,-2) (red line), and so is (2,2) compared to (2,0). It is essential to investigate the effect of the different system parameters on the asymptotic behavior of the entanglement between the two atoms, which is considered in Fig. 1(d). While, as we concluded from Fig. 1(a) It is always very insightful to monitor the system dynamics using 3-dimensional plots, which shows the behavior of the system over a wide continuous range of the parameter values.
In Fig. 2(a) and (b), we depict the time evolution of the entanglement versus the Ising coupling parameter, where −10 ≤ J ≤ 10 at early times and asymptotically respectively.
As can be noticed, increasing the Ising coupling value may remove the ESD at early times as shown in Fig. 2(a), but eventually causes the entanglement to vanish completely without revival as illustrated in Fig. 2(b), where the ESD period length increases with J until eliminating the entanglement completely.
Testing the system dynamics starting from another maximally entangled state, the Bell anti-correlated ψ Ba = (|g 1 |e 2 + |e 1 |g 2 )/ √ 2, showed slight differences but not that significant, in the behavior of the entanglement, ESD, and the atomic population compared with the correlated Bell state at the early times, however, asymptotically the ESD can be removed entirely from the system by increasing the Ising coupling value as illustrated Fig. 2(c). Interestingly, as can be observed in Fig. 2(b) and (c), the entanglement behavior is symmetric as the values of J changes from the positive to the negative values, however the symmetry center is not J = 0 but around 0.9.
In Fig. 3, we illustrate the atomic population inversion for one of the two atoms, starting from an initial correlated Bell state ψ Bc with the radiation field is in the coherent state. which agrees with our observations from the entanglement dynamics. By comparing Fig. 1 and Fig. 3, it is clear that the collapse periods match the ESD periods, while the revivals synchronize with the entanglement revivals from death, i.e. the entanglement revives from death when the atoms exchange energy with the field.
Remarkably, there is an inverted symmetry existing when one compares the atomic population in the two cases, where each collapse line at specific parameter values , say (2,-2), is far from the (0,0) line equally for ∆ = 2 and −2 but at the opposite sides. Also, as we pointed in Fig. 3(a) Fig. 4(d). In Fig. 5(a),  Fig. 5(b), where the contour plot of the asymptotic entanglement is presented.
The atomic population is illustrated in Fig. 6(a) we depict the asymptotic dynamics of the population inversion at the parameters choice (4,-4), at the detuning values 0, 2 and -2 respectively in the top, middle and bottom panels respectively. As can be seen, in all cases, the collapse periods shrinks until completely disappearing, where the dynamics profile becomes a continuous irregular oscillation, which resembles the entanglement function asymptotic behavior and provides an explanation for it, where the exchange of energy between the two atoms and the field becomes continuous with no interruption, enhancing the mediated entanglement between the two atoms.

C. Disentangled initial state
Now we turn to another different type of initial state, which is completely separable (disentangled), where the state of the two atoms is a linear combination of all the basis states, namely ψ L = (|g 1 |g 2 + |g 1 |e 2 + |e 1 |g 2 + |e 1 |e 2 )/ √ 4, while the radiation field is in a coherent state as previously described.
This interesting initial state maintains zero entanglement between the two non-interacting atoms at zero detuning forever, as shown in Fig. 7(a) (blue line). Turning on a low AF Ising coupling, (0,2), causes the entanglement to rise initially from zero reaching a low peak before eventually suffering ESD, without any future revival (Green line). However, applying the case (0,0) before the entanglement eventually vanishes again completely once and for all.
The negligible effect of changing the detuning value can be seen also in Fig. 7(d respectively. Interestingly, when the system starts from the separable state ψ L , the Ising coupling parameter has a very peculiar impact, it can be used as a switch that turns on and off the entanglement asymptotically by crossing a critical value J c , where the entanglement vanishes for coupling values higher than J c , as illustrated in Fig. 8. The behavior of the population inversion differs significantly as the detuning parameter is varied, as shown in Fig. 9. The zero detuning case is illustrated in Fig. 9(a), which exhibits a much smaller oscillation amplitudes compared with the non-zero detuning cases.
The non-zero detuning, positive or negative, causes the collapse lines to spread over a wider range compared with the zero detuning case, as illustrated in Fig. 9 The zero detuning profile, which looks different from the non-zero detuning, still shows a collapse-revival behavior, as shown in the upper panel of Fig. 9(d). In all cases, the collapse periods disappear asymptotically and the profile turns to an oscillatory one to synchronize with the entanglement asymptotic dynamics. The separable initial states |e 1 |e 2 and |g 1 |g 2 do not show ESD upon evolution and therefore are not discussed here.
Clearly, comparing the results presented in the last three subsections shows that the dynamics of the system and its asymptotic behavior vary significantly depending on the initial state of the system. In general, when a single spin-1/2 particle is exposed to an external magnetic field, the magnetic dipole moment of the particle interacts with the field and precesses around its direction with a constant angle that depends on the initial state determined by the relative orientation (initial state) of the spins, while the frequency by the interaction strength that depends on the anisotropy of the spin-spin coupling. This explains the critical impact of the initial state, as well as the relative strength of the coupling in the x, y and z-directions, on the subsequent dynamics of the system and its asymptotic behavior.

IV. QUANTUM CORRELATION BETWEEN THE TWO ATOMS AND THE RA-DIATION FIELD
It is very insightful to study the quantum correlation between the two atoms, A 1 and A 2 , ensemble and the radiation field F , where the composite system is represented as F ⊗ (A 1 A 2 ).
This correlation can be utilized to shed some light on the reported behavior of the population inversion and entanglement dynamics that we have just discussed. Rungta and d B in an overall pure state by generalizing the spin-flip operation to apply to higher dimensional systems [58]. The resulting quantity, called the I-concurrence, is given by where ν A and ν B are arbitrary scale factors that, in general, depend on the dimensions of the subsystems. The tangle τ of a bipartite system in a pure state with arbitrary subsystem dimensions reads In our case the tangle takes the form In Fig. 10, we depict the dynamics of the system tangle starting from the maximally entangled correlated Bell state. As can be noticed, at all the different detuning values 0, 2 and ior, which we did not present here. In general, the different panels of Fig. 10 demonstrate that the quantum correlation between the ensemble of the two atoms and the radiation field minimizes during the population collapse and the ESD periods and maximizes during their revivals, when the energy exchange between the atoms and the field takes place. The tangle dynamics starting from an initial W-state is illustrated in Fig. 11(a) and (b), where it shows a smaller oscillation amplitude compared with the Bell state case, at all detuning values but similarly it evolves to a quasi steady state asymptotically with an almost constant mean value. The detuning variation is not showing a noticeable effect on the tangle behavior in this case. In Fig. 11(c) and (d), the time evolution of the tangle starting from the separable initial state is explored. The tangle starts from zero as in the two previous cases but this time it shows a much weaker oscillatory variation with a smaller amplitude and a faster increasing mean value, which reaches a quasi-steady state asymptotically as well. The effect of the non-zero detuning is to shift the tangle values up and slightly split out the tangle lines at the different parameters choices from each other.

V. CONCLUSION
We studied a system of two two-level atoms (qubits) interacting off-resonance with a single-mode radiation field. We considered the two atoms to be coupled to each other through dipole-dipole and Ferromagnetic (anti-Ferromagnetic) Ising interactions. We presented an exact analytic solution for the time evolution of the system that spans its entire parameter phase space starting from any initial state. We utilized the analytic solution to study the entanglement dynamics, between the two atoms, and its asymptotic behavior, particularly when the system starts from an initial state that leads to entanglement sudden death (ESD) upon evolution. The combination of the Ising and the dipole-dipole interaction was found to be powerful in manipulating the ESD compared with either one of them separately, especially at non-zero detuning. Their combined impact was found to vary significantly depending on the type of Ising interaction (Ferromagnetic or anti-Ferromagnetic) and the initial state 26