Manipulating entanglement sudden death in two coupled two-level atoms interacting off-resonance with a radiation field: an exact treatment

We study a model of two coupled two-level atoms (qubits) interacting off-resonance (at non-zero detuning) with a single mode radiation field. This system is of special interest in the field of quantum information processing (QIP) and can be realized in electron spin states in quantum dots or Rydberg atoms in optical cavities and superconducting qubits in linear resonators. We present an exact analytical solution for the time evolution of the system starting from any initial state. Utilizing this solution, we show how the entanglement sudden death (ESD), which represents a major threat to QIP, can be efficiently controlled by tuning atom-atom coupling and non-zero detuning. We demonstrate that while one of these two system parameters may not separately affect the ESD, combining the two can be very effective, as in the case of an initial correlated Bell state. However in other cases, such as a W-like initial state, they may have a competing impacts on ESD. Moreover, their combined effect can be used to create ESD in the system as in the case of an anti-correlated initial Bell state. A clear synchronization between the population inversion collapse-revival pattern and the entanglement dynamics is observed at all system parameter combinations. Nevertheless, only for initial states that may evolve to ESD, the population inversion revival oscillations, where exchange of energy between the atoms and the field takes place, temporally coincide with the entanglement revival peaks, whereas the population collapse periods match the ESD intervals. The variation of the radiation field intensity has a clear impact on the duration of the ESD at any combination of the other system parameters.


I. INTRODUCTION
The great interest in realizing quantum information processing (QIP) systems in the last few decades, capable of performing efficient quantum simulation and quantum computing tasks [1], led to a significant progress in engineering new quantum systems that are considered very promising candidates for playing the role of a qubit. These developed artificial atomic systems (such as semiconducting quantum dots and superconducting circuits) in addition to customized natural atomic systems (such as Rydberg atoms and trapped atoms, ions and molecules), in contrary to the natural conventional atoms, enjoy a strong coupling with a similar type of system or even with a different type (when implemented in a hybrid system) [2][3][4][5], through direct or mediated interaction. A significant breakthrough in the role of artificial atomic systems in the QIP arena was achieved when superconducting qubits were successfully embedded in a superconducting microwave resonator [6][7][8][9], paving the way for the new paradigm of circuit quantum electrodynamics (cQED). Another important step came out when an architecture for quantum computation using cQED was introduced by Balis et al [10], where they proposed to use a one-dimensional (1D) transmission line resonator consisting of a full-wave section of superconducting coplanar wave guide to play the role of a cavity and a superconducting qubit as the atom. They showed that this structure can be customized to access the strong coupling regime between the field and the qubit in analogy to what can be achieved in optical cavity quantum electrodynamics (CQED). Furthermore, they demonstrated that the proposed architecture can be efficiently utilized for the coherent control, entanglement, and readout of qubits in QIP. Particularly, they emphasized the possibility of generating tunable entanglement between two qubits that are few centimeters apart within the resonator, where they introduced, in addition to the qubits-field interaction, an effective coupling between the two qubits mediated by the virtual excitation of the resonator. Soon after, a strong coupling between a single photon and a superconducting qubit using cQED was realized experimentally [11]. Several experimental works demonstrated that two close superconducting qubits can be directly coupled via local interactions [12][13][14][15][16][17][18]. Latter, the coupling of two distant superconducting qubits mediated by microwave photons confined in a superconducting transmission line was reported [19]. In fact, proposals for similar schemes in cavity QED was introduced and implemented particularly for coupling (directly or indirectly) Rydberg atoms inside optical cavities [20][21][22][23][24][25][26][27].
Currently, there is a great progress in developing superconducting multi-qubit circuits with long coherence time embedded in 1D, 2D and 3D superconducting resonators, which can perform high-fidelity quantum gates; extended review and references can be found in [5]. Remarkably, not only superconducting qubits can be integrated into superconducting resonators but also atoms (trapped and Rydberg) and spins (in Quantum dots or solid state impurities) producing hybrid circuit QED [3,5]. In these hybrid circuits, superconducting qubits can have either direct interaction with atoms (or spins) through electric or magnetic fields or indirect one mediated by the radiation field. These hybrid structures are of great interest for QIP as it combines the advantages of the two sides, insensitivity and long coherence time of atoms (spins) and rapid processing in superconducting circuits. Very recently, Nguyen et al [28] have proposed a new approach for analogue quantum simulation of spin arrays based on laser-trapped circular Rydberg atoms, benefiting from their long coherence time and insensitivity to collisions and photo-ionization, relying on the available state-ofthe-art experimental techniques. They showed that the strong coupling between the atomic dipoles can be utilized to simulate a spin-1/2 XXZ chain Hamiltonian with fully tunable nearest neighbor coupling over a wide range and studied the adiabatic time evolution of the chain. Furthermore, they suggested that this scheme can be implemented in CQED to overcome many of the challenges in Rydberg-atom CQED or even can be utilized in Hybrid cQED experiments by integrating superconducting circuits and laser-trapped Rydberg atoms.
In general, for a composite system with several interacting qubits (of same type or different) coupled to a single radiation mode, the Hamiltonian assumes the form [3,5,10,24] The first and second terms in the Hamiltonian represent the free quantized radiation field and the non-interacting qubits while the third and fourth terms represent the qubit-field and qubit-qubit interactions respectively. Ω and ω i are the frequencies of the single-mode radiation field and the ith qubit 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 spin operators representing the ith qubit. There are different regimes of the coupling strength between the radiation field and the qubit, which are mainly determined by the ratio between λ i 1 and the other energy scales of the composite system namely Ω, ω i , κ and γ i , where κ and γ i are the decay rates of the resonator (cavity) and the ith qubit respectively. When λ i 1 << Ω, ω i , κ, γ i , the system is in the weak coupling regime, whereas for κ, γ i << λ i 1 << Ω, ω i , it is in the strong coupling regime. In these regimes the rotated wave approximation (RWA) [29,30] is valid, which when applied converts the third term in the Hamiltonian into the usual Jaynes-Cummings y . The RWA fails in the ultra-strong and deep-strong regimes where λ i 1 becomes of the same order of magnitude as Ω or higher. The qubit-qubit coupling constant λ ij 2,ν varies depending on the type of qubits and the nature of the coupling (photon mediated, direct capacitive or conductive, etc.) between the qubits as well as on the field-qubit coupling strength regime but mostly is modeled as a spin-1/2 XY Z Heisenberg exchange interaction [5]. The science of quantum information processing is not only concerned with the fine preparation of such systems in a well defined state but also the controlled time evolution of them while preserving the entanglement among the different parts of the composed system within its coherence time [1]. In QIP, the classical data are mapped on the Hilbert space of the processing system, then the time evolution of the system is followed and at the end a readout measurements of the system registers are performed and the classical output is analyzed.
Quantum entanglement is considered to be the physical resource crucially needed for manipulating linear superposition of the quantum states of the different constituents of composite systems to implement the proposed schemes in QIP. However, the inevitable interaction between the quantum system and its environment leads to loss of entanglement in what is known as the decohering process. While gradual loss of entanglement that obeys the half life-time law and evolves to a state of even quite small entanglement can be treated and even reversed using different approaches such as quantum error correction [31,32], decoherence free subspace [33,34] and quantum measurement reversal [35], entanglement sudden death represents a major threat to QIP as the loss takes place very abruptly leading to a state of zero entanglement [36,37]. The entanglement dynamics in composite systems of interacting quantum systems (qubits) in presence of external magnetic fields coupled to different types of environments in absence of radiation fields have been studied intensively before, where different approaches for creating, enhancing, controlling and protecting entanglement against decoherence and dissipation were investigated and discussed [38][39][40][41][42][43][44][45][46][47][48][49][50][51][52].
The great progress in developing new systems that are promising candidates for QIP and controllable in cavity or circuit QED sparked huge interest in studying entanglement dynamics and sharing in composite systems containing two-level atoms (qubits) coupled to radiation fields. The entanglement sharing among the different constituents of a composite system containing two uncoupled atoms interacting with a radiation field was investigated and constrains on the shares were provided [53]. A big turn in the field took place when Yu and Eberly [36] showed that the entanglement between two non-interacting, initially entangled, two-level atoms vanishes within a finite period of time, where each atom was coupled to a different environment (cavity), it was called by them, for the first time, entanglement sudden death (ESD). Latter, the effect of coupling between the two atoms on the revival of the vanished entanglement was studied starting from a particular initial state, where the radiation field was treated as an environment represented by a vacuum state and the master equation of the composite system was solved [54,55]. The entanglement sudden death in two initially entangled, uncoupled, atoms under the effect of a noisy classical environment (stochastic magnetic field) was studied both collectively and separately starting from a mixed state [56]. More works were devoted to studying the ESD in systems of two uncoupled atoms, each one was in a separate independent cavity, where the effect of different initial types of Bell states was investigated [57][58][59][60]. A double JC model out of resonance with the fields was studied too, where two non-identical, uncoupled atoms were considered in two remote cavities and each atom was coupled to a single mode radiation field in its cavity [61].
It was shown that asymmetry can be an advantage for entanglement creation and evolution and the off-resonance condition may, for certain initial states, enhance entanglement transfer (between the atoms and the fields) and prevent ESD.
Recently, the effect of coupling between atoms (qubits) on the system entanglement was brought to the focus of interest, as a result of the newly engineered systems that enjoy strong interactions with each other as we explained earlier, and also due to both its central role in the system dynamics and its practical impact on QIP protocols. The Entanglement dynamics and population difference in a system of two interacting spins have been studied under the effect of coupling to Ohmic and subohmic bosonic environment [62]. It was shown that there is a bath-induced spin-spin coupling and the spin-spin entanglement dynamics can be controlled by detuning the coupling to the bosonic bath but depends critically on the initial state of the system. ESD was studied for a system of two identical interacting atoms in a double mode radiation field with frequencies ω 1 and ω 2 at resonance condition, where the atom energy gap ω 0 was such that ω 0 = ω 1 + ω 2 [63]. An analytical solution for the problem was introduced and it was shown that the time evolution of the entanglement and the ESD depends significantly not only on the initial amount of entanglement in the system but also on the type of initial state. It was demonstrated that ESD was reduced after introducing the dipole-dipole interaction between the two atoms. Also the dipole-dipole interaction was found to enhance the entanglement between the two atoms starting from a W-like initial state as long as the interaction is stronger than the atom field interaction at resonance [64]. The analytical solution for the same system was presented at off-resonance but for non-interacting atoms and it was shown that the non-zero detuning may enhance the entanglement and suppress the sudden death. An exact analytical solution for the density matrix of two identical interacting atoms coupled to a single mode radiation field at resonance was presented in [65], where the effect of the interplay between the atom-atom interaction and the coupling, at resonance, to the radiation field on both of the entanglement and purity of the system was investigated thoroughly. Recently, a non-linear model for two interacting atoms coupled to a radiation field was introduced by Sanches et al. [66]. They considered two types of interaction between the two atoms, dipole-dipole and Ising and represented the radiation field as a coherent superposition of number states. The non-linearity was introduced by introducing and multiplying photon-number-dependent function everywhere times the photon operators in the system Hamiltonian. They diagonalized the interaction Hamiltonian using a basis containing three states and solved Schrodinger equation to obtain the time evolved states at any time t in terms of the basis three-states. Although they started with a Hamiltonian for a generic system where they assumed non-identical atoms out of resonance with the field, but eventually when they came to study specific cases, they assumed identical atoms at resonance condition (zero detuning) with the field, where they provided an exact analytical solution for that case. They studied the time evolution of population inversion, purity of the atomic state and entropy of the radiation field, where they always assumed an initial state that is one of the three states in the implemented basis to be capable of performing the calculations. They carried out numerical calculations to study the entanglement between the two atoms using the concurrence function. Another work studied a non-linear model of two atoms interacting with a radiation field far from resonance where the interaction was considered intensity dependent, but the two atoms were uncoupled [67].
Very recently, schemes to avoid ESD in an evolving system of two coupled qubits exposed to a common vacuum environment using Local unitary operations was introduced, where the Lehemberg-Agarwal master equation was implemented under the Markovian approximation [68]. In a very relevant work, Gywat et al studied a system of two coupled two-level atoms (qubits) interacting with an off-resonance single mode radiation field [24]. In order to study the system dynamics, they applied a generalized Schrieffer-Wolf transformation to the system Hamiltonian and provided a perturbative analytical solution in the limit of weak interqubit coupling and an exact solution by numerical diagonalization of the Hamiltonian.
They demonstrated that the state of the two qubits can be read out using the cavity mode dispersion within the perturbative regime and studied the effect of the interqubit coupling on a cavity-mediated two-qubit gate. In another relevant work, the dynamics of entanglement in a system of two uncoupled spins (qubits) interacting with a single mode radiation field in an optical cavity was studied [69]. The time-dependent quantum correlation of Clauser-Horne-Shimony-Holt type was derived and used versus an entanglement measure, concurrence, to test the Bell inequality (BI) violation. It was shown that the interaction with the field induces decoherence and coherence revival that is characterized by the BI violation, where the important role of the field intensity in the decohering process and the inequality test was pointed out. It is essential to emphasize here that although the pioneering work of Yu and Eberly [36] introduced the ESD in the case of an open system, where the ESD took place as result of a dissipative effect (vacuum noise) and also in another work they showed a similar behavior caused by classical noise [56], the same two authors have reported the same phenomenon, ESD, taking place in a closed system. The system contained two uncoupled atoms (A, B), where each one of them is in a separate cavity with a single mode radiation field in each and the two cavities (a, b) are not coupled to each other or to any other environment [57,59]. As they explained in their work, the ESD takes place (zero value of the concurrence C A,B for finite time) in the system, although there is no interactive decoherence, by transferring the entanglement to one or more of the constituting pairs of  [60,61,63,64].
In this paper, we consider two identical coupled two-level atoms (qubits) symmetrically interacting with a single-mode quantized radiation field. We present an exact analytical solution for the the time evolution of the system at either resonance or off-resonance (nonzero detuning) interaction between the atoms and the field starting from any initial state of the composite system. The coupling between the two atoms is modeled as a spin 1/2 isotropic XY exchange interaction. As we have discussed before, this system has been treated in several works, in absence of either the coupling between the atoms or the nonzero detuning and was studied in presence of both only in a perturbative way. Our general exact analytical solution provide a mean for studying the different dynamical properties of the system while spanning the whole system parameters space taking into account the interplay among all of them without excluding any. We utilize this solution to study the entanglement dynamics of the system in general and specially in the case of ESD. Our goal is to test the effect of the coexistence and tuning of atom-atom coupling and offresonance atom-field interaction, which was not possible before, on the system dynamics and particularly the manipulation of ESD. On the other hand studying this model is not only important for its own sake as a system of coupled localized spins interacting (offresonance) with a bosonic bath but also as an enlightening step in exploring cavity (circuit) QED, with its crucial impact on QIP schemes as explained earlier. We demonstrate how these two interactions when applied separately or combined can be used to reduce, eliminate or create entanglement sudden death in the system, depending crucially on the initial state.
While one or both of the interactions may not be effective individually in treating some cases of ESD, combining them proves to be significantly different. We consider different initial states of the system of special practical interest, which contain maximum, partial or zero entanglement between the two atoms. Also, we show that there is a strong synchronization between the population inversion collapse-revival pattern and the entanglement dynamics at all system parameters combinations. However, only for initial states that may evolve to ESD, the exchange of energy between the field and the atoms enhances the entanglement between the two atoms inducing revival peaks with rapid oscillation, while the population collapse periods synchronize with that of the ESD. Also, the impact of the variation of the radiation field intensity on the ESD duration periods is investigated. This paper is organized as follows. In Sec. 2, we discuss our model. In Sec. 3, we present our exact analytical solution for the time evolution of the system. We implement our solution to study the dynamics of entanglement and atomic population inversion, starting from different initial sates, in Sec. 4. We conclude in Sec. 5.

II. THE MODEL
We consider a model 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 the same single-mode quantized radiation field with the same coupling constant λ 1 . The coupling between the two atoms is modeled as an isotropic XY exchange interaction between two spin-1/2 particles with coupling strength λ 2 . This system can be realized in either cavity or circuit QED structures as illustrated in Fig. 1. The system Hamiltonian assumes the same form as Eq. (1) except that it is for only two qubits and therefore reduces tô As in Eq. (1), the first and second terms in the Hamiltonian represent the free quantized radiation field and the non-interacting two atoms while the third and fourth terms represent the atom-field and atom-atom interactions respectively. Ω and ω • are the frequencies of the single-mode radiation field and the quantum system transition respectively,â † ,â,σ (i) ± and σ i z have the same meaning and roles as in Eq. (1). Using the Heisenberg equation of motion, which for any operatorÔ, where = 1, reads one obtains the following equations for the field and atom operators Assuming that initially the atoms are in a pure state and the field is in the coherent state, the wave function of the composite system at t = 0 can be written as where a, b, c and d, are arbitrary complex quantities that satisfy the condition 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 the two atoms are in excited state and the field has n photons, |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 and so on. The states of the quantum system satisfy the relationsσ (i) and D n+2 (t) can be obtained by solving the Schrödinger equation of the composite system, which will be discussed in the next section. Once we obtain the system wave function |ψ(t) , we can calculate the composite system density matrixρ(t) = |ψ(t) ψ(t)|. The reduced density matrix of the two atoms, ρ red (t), can be obtained by tracing out the field

III. THE ANALYTICAL SOLUTION
We devote this section to solve the Schrödinger equation of the system and provide an exact analytical expression for the time-dependent coefficients of the system wave function.
We start by rewriting the Hamiltonian (Eq. (2)) aŝ whereĤ andN where ∆ = ω • − Ω is the detuning parameter. Using Eqs. 3 and 4, one can show thatN, which represents the total number of excitations in the system, is a constant of motion, which justifies the use of its eigenstates as a basis for the expansion of the system wavefunction in Eq. (8). It is more convenient to work in the interaction picture where we defineV I = UĤ intÛ † withÛ = e iĤ•t . As a result, we obtain Now, substituting |ψ(t) and V I (t) into Schrödinger equation it yields a system of coupled differential equations which after some calculations becomes ...
with a solution where and where and Finally, the solution of the set of differential equations (16) takes the form where the initial values of the coefficients are given by As can be noticed, for Eqs. (25) to represent an acceptable physical solution, the param- and v 2 in Eqs. (22) to be appropriate for the solution whereas the others represent a nonphysical solution. In fact, each one of the two quantities v 1 and v 2 will have three, generally complex, roots. Therefore v 1 and v 2 defined by Eqs. (22) have nine possible combinations, only six of them lead to physically acceptable solution. Nevertheless, very fortunately these six combinations enable us to span the whole parameter space of the system. Finally the reduced density matrix of the two atoms defined by Eq. (9) can be obtained, utilizing that

IV. DYNAMICS OF ENTANGLEMENT AND ATOMIC POPULATION INVER-SION
In this section we implement our exact solution to study the dynamics of the bipartite entanglement between the two atoms and the atomic population inversion starting from different initial states of particular interest. For convenience, we set = 1, λ 1 = 1 and set represent the other parameters (λ 2 and ∆) in units of λ 1 . The entanglement between the two quantum system can be quantified with the help of the concurrence function C(ρ red ) as proposed by Wootters [70], which is related to the entanglement of formation E f through the formula where E is defined as h is the Shannon entropy function and the concurrence can by calculated from The ε i arranged in decreasing order are the square root of the four eigenvalues of the non-Hermitian matrix Whereρ red is the spin flipped state defined as Here ρ * red is the complex conjugate of ρ red andσ y is the Pauli spin matrix in the y direction. Both of C(ρ red ) and E f (ρ red ) go from 0 for a separable state to 1 for a maximally entangled state.
Atomic population inversion is defined as the expectation value of the operatorσ z or the difference between the probabilities of finding the atom in its excited state and ground state.
To investigate the atomic inversion we first calculate the reduced density matrix of any one of the two identical atoms, say the first,ρ 1 (t) by tracing out the other one in the two atoms reduced density matrixρ red (Eq. 27), which leads tô where Therefore, for the first atom A. Maximally entangled initial (Bell) states In Fig. 2, we explore the dynamics of entanglement and population inversion, in terms of the scaled time τ = λ 1 t, starting from a correlated initial Bell state ψ Bc = (|e 1 |e 2 + |g 1 |g 2 )/ √ 2 with the radiation field is in a coherent state. Starting from such an initial state the system shows ESD, where the entanglement changes abruptly from a non-zero to an exact zero value, which is illustrated in the different panels of the figure. In Fig. 2(a), we test the effect of the field intensity, by changing the average number of photonsn, on the entanglement dynamics and ESD time intervals for uncoupled atoms at resonance with the radiation field. As can be noticed, by increasing the intensity of the radiation field from n = 20 to 50 and then 100, the sudden death interval time increases considerably. The inset plots of Fig. 2(a)  of entanglement and population inversion is depicted in Fig. 2(b) for uncoupled atoms at resonance with the field. In the upper two panels, we compare the entanglement dynamics and the atomic population within the interval 0 ≤ τ ≤ 100, which shows three entanglement revival peaks. In the lower two panels, we compare them again but after zooming into a smaller time interval 0 ≤ τ ≤ 40 and much smaller ranges of E f and σ z , so we can focus on the first entanglement revival peak and the corresponding atomic population dynamics. As can be noticed in the upper most panel, the two-atoms start in an initial state where they are maximally entangled with each other but abruptly they lose their entanglement, showing ESD, and maintain this state for a finite period of time. However, they gain entanglement back with a revival peak as shown, at τ ∼ 27, which as we pointed out in our introduction is due to the fact that the loss of entanglement is not due to a dissipative effect but a transfer of entanglement to the atom-field subsystems, which were initially disentangled.
In the second upper panel, the atomic population shows the usual collapse revival pattern, although it doesn't collapse to zero but a constant value. This behavior of the entanglement and atomic population is repeated periodically. Clearly, the collapse periods coincide with the zero entanglement intervals, while the population revivals temporally coincide with the entanglement revivals. This behavior can be better recognized in the lower two panels where one can see that for the atoms entanglement to drop to zero a rapid oscillation of the atomic population takes place at the same time, indicating an exchange of energy between the atoms and the field is taking place. Then the zero entanglement state is maintained for a finite time before an entanglement revival occurs, where the entanglement increases from zero to a peak then back to zero, accompanied by a rapid oscillation of the atomic population that starts and finishes within the same time interval. Obviously, the entanglement between the two atoms is gained back from the atom-field subsystems before getting lost to them back through exchange of energy between the atoms and the field. In Fig. 2( increases. The effect of the atom-atom coupling at zero detuning is illustrated in Fig. 2(d), where a small coupling value, λ 2 = 3, has no effect on the ESD (dashed red line), while a higher value, 5, partially eliminates the ESD (dotted green line). Further increase of the coupling to λ 2 = 7 completely eliminates the ESD (dash dotted violet line). The inset plot in Fig. 2(d) shows the entanglement collapse-revival pattern over a longer time interval in presence of atomic coupling.
The combined effect of the atom-atom coupling and non-zero detuning on the entanglement and population inversion is considered in Fig. 3. While, as we have observed in Fig. 2, the non-zero detuning can neither remove nor reduce the ESD for uncoupled atoms, it may  Fig. 3(a)) and the population revival oscillations, for the coupled atoms at non-zero detuning, which again indicates that the entanglement revives from death and vanishes again as a result of the exchange of energy between the field and the atoms. However, the reduction or removal of ESD due to atom-atom coupling and non-zero detuning is not accompanied by any atomic population oscillation, which means they induce entanglement, though it is very weak, away from any energy exchange between the atoms and the field.
In Fig. 4, we consider a different maximum entanglement initial state that does not yield upon evolution an ESD, for uncoupled atoms at resonance with the field, namely the anti- correlated Bell state ψ Ba = (|g 1 |e 2 + |e 1 |g 2 )/ √ 2. Increasing the radiation field intensity fromn = 20 to 50 and then to 100 results in very similar impact on the system dynamics to the previous case, the time interval of the very small constant entanglement (rather than zero entanglement in the previous Bell state) increases with radiation intensity as shown in Fig. 4(a). The inset plots of Fig. 4(a) illustrates how the constant (mean) entanglement value decreases as the radiation intensity is increased and reaches a value as low as 10 −9 atn = 100. The synchronization between the periods of constant entanglement and the constant population is shown in Fig. 4(b), which again emphasizes that a steady behavior of the atomic population, where there is no exchange of energy between the atoms and the field leads to a quite small constant entanglement value, while the population revival oscillation boosts the entanglement considerably.
The effect of non-zero detuning on the entanglement between the uncoupled atoms is The combined effect of λ 2 and ∆ on the entanglement and the atomic population is illustrated in Fig. 5(a) and (b) respectively. In the inset plots of Fig. 5(a), one can see that setting λ 2 = 3 and ∆ = 0 raises the constant entanglement value from the order of 10 −9 (for uncoupled atoms at resonance with the field as shown in Fig. 4(c)) to another value that is 10 5 higher (solid blue line). Nevertheless, turning on detuning at ∆ = 2 leads to entanglement death with intermediate reviving peaks (dashed red line), but these peaks turn to a single narrow one with a higher maximum at ∆ = 4 (dotted green line).
However, applying a higher coupling, λ 2 = 5, acts to overcome the detuning effect and partially eliminates the entanglement death while shifting the entanglement reviving peaks to earlier times (dash dotted violet line). In Fig. 5

B. Partially entangled initial (W) state
A partially entangled initial state that yields ESD upon evolution, for uncoupled atoms at resonance with the field, is the W-like state ψ W = (|g 1 |g 2 +|g 1 |e 2 +|e 1 |g 2 )/ √ 3, which is considered in Fig. 6. In Fig. 6(a), as can be seen, increasing the radiation field intensity leads to longer ESD time intervals and reduces the entanglement oscillation. Clearly, the ESD time intervals in the current case is much smaller than the ones corresponding to the correlated Bell state. Testing the atom-atom coupling effect on the entanglement dynamics at zero detuning is depicted in Fig. 6(b), where increasing the coupling strength from 1 to 3 and then to 5, increases the sudden death time interval, which is clearly illustrated in the inset plot, and makes the entanglement reviving peaks narrower with a smaller maximum value, i.e. the atomic coupling enhances ESD. On the other hand, applying a non-zero detuning to uncoupled atoms removes the ESD partially, at ∆ = 3 (the dotted green line), or even completely, at ∆ = 5 (dash dotted violet line), as shown in Fig. 6(c). Therefore, the off-resonance interaction between the field and the atoms can be utilized to completely terminate ESD and the entanglement oscillation in this case indicates a reduction in the atom-atom entanglement (by transfer to the other subsystems) without completely vanishing before being gained back, which is repeated periodically. The combined presence of atomic coupling and detuning is considered in Fig. 6(d), in which the atomic coupling suppresses the entanglement, while the non-zero detuning acts the opposite way, eliminating the entanglement death and enhancing the entanglement peaks. As one can notice, setting λ 2 = 1 but ∆ = 5 is enough to completely remove the entanglement death (dotted green line) but increasing the coupling to λ 2 = 5 pushes the entanglement peaks down towards the zero value (dash-dotted violet line). The inset plot in Fig. 6(d) compares the long time behavior of the entanglement for two coupled atoms at resonance with the field, λ 2 = 1, ∆ = 0 (solid blue line), versus the same two coupled atoms but at non-zero detuning ∆ = 5, as can be seen the non-zero detuning removes the entanglement death and asymptotically sustains the entanglement oscillation.
In Fig. 7 we discuss the variation in dynamics of the atomic population as a result of changing the system parameters. Increasing the radiation intensity, not only increases the ESD time interval, as illustrated in Fig. 7(a), but also shifts the population mean value up towards the zero as shown in the inset plot. Figure 7 Figs. 7(b)), one can see that, at either zero or small ∆, ESD is present but its interval is interrupted with a smooth entanglement peak without a corresponding population oscillation. However, when an Entanglement oscillating revival peak appears latter, it is accompanied by an atomic population revival oscillation. On the other hand, when ESD is significantly reduced or even removed, the atomic population collapse does not correspond to either a zero or a quite small entanglement but a peak, while the population revival oscillation correspond to a minimum in the entanglement with a local rapid oscillation, in contrary to what was observed in the correlated Bell state case.
Since both of the entanglement and atomic population are very sensitive to the variation in the detuning parameter ∆, as can be noticed in the last two figures, it would be interesting to monitor their dynamic profile over a wide range of (negative and positive) values of ∆, which is illustrated in Fig. 8. As can be noticed, in Fig. 8(a), only positive values of ∆ can eliminate the ESD and enhance entanglement as ∆ increases. On the other hand, the atomic population reaches negative values for negative detuning and positive for positive ones, while the amplitude of the revival oscillation increases as ∆ increases as shown in Fig. 8(b).

C. Separable initial states
Now we turn to another type of initial states, which is completely separable, namely ψ e = |e 1 |e 2 , where both atoms are originally in their excited state and the field is in a coherent state. The time evolution of entanglement and population inversion, starting form that initial state is depicted in Fig. 9. For two uncoupled atoms at resonance with the field atn = 20, the entanglement starts from a zero value making a rapid oscillation with a very small amplitude before turning to a larger oscillation with big amplitude then turns to a rapid oscillation again at τ ≈ 20 (solid blue line), as illustrated in Fig. 9(a). But as the intensity of the radiation field is increased, ton = 50 then to 100 (dashed red and dotted green lines respectively), the amplitude of the oscillation slightly increases whereas the coupling parameter λ 2 with the two atoms are initially in a disentangled state ψ e = |e 1 |e 2 and the field is in a coherent state forn = 100 and ∆ = 5.
the frequency decreases to almost its half value. The inset plot of Fig. 9(a) shows a magnified view of the minima of the entanglement oscillations at differentn values, which illustrate that none of them reaches zero value. The corresponding change in the population inversion as the field intensity increases is shown in Fig. 9(b). As can be noticed, higher intensity results in longer collapse time and a slight down shift in the constant collapse value towards the zero. In the inset plots of Fig. 9(b), we gave a close look and compare the dynamics of E f and σ z . Clearly there is a strong correlation between the two, where a minimum value of entanglement with a rapid oscillation is corresponding to a revival period of the population inversion, whereas a maximum entanglement corresponds to a collapse period, where is no exchange of energy is taking place between the field and atoms. Turning on the interaction between the atoms, λ 2 = 5, at resonance is shown in Fig. 9(c)(dashed red line). The entanglement profile is slightly different from the uncoupled case except for a big oscillation that takes place after time τ = 60. On the other hand, applying a non-zero detuning, ∆ = 5, for uncoupled atoms, the entanglement oscillation shifts slightly to the right with lower peaks but the rapid oscillation appears earlier The entanglement oscillation indicates that entanglement is transferred back and forth to the other subsystems, but it is not always accompanied by an atomic population revival oscillation, which means the entanglement sharing is not always mediated by atom-field energy exchange. In Fig. 10, we depict the time evolution of entanglement and population inversion over a wide range of λ 2 in (a) and (b) respectively for ∆ = 5. One can see how increasing λ 2 spreads out and splits the population revival oscillation which also takes place at the same time for the corresponding entanglement rapid oscillation.
Finally, we consider an interesting separable initial state, which 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. As one can see in Fig. 11(a), for uncoupled atoms at resonance with the field atn = 20 starting from ψ L (solid blue line), the entanglement of the system shows at early time a number of very narrow short spikes before completely vanishing then after a long time revives again to much shorter spikes and keeps repeating this behavior continuously. Raising the field intensity ton = 100 (dashed red line) doesn't lead to a noticeable change in the entanglement profile. These two cases are depicted at a magnified scale in the left inset plot of Fig. 11(a). Setting a non-zero detuning in the system, ∆ = 2, (dotted green line) enhances the entanglement value and reduces the entanglement death period. However, turning on the coupling between the two atoms, even at a very small strength λ 2 = 0.2 at zero detuning, completely eliminates the entanglement death, which shows a collapse-revival like behavior (dash-dotted violet line), where it doesn't collapse to zero value but a constant one (≈ 7.5 × 10 −5 ). The early time behavior of the entanglement in this two last cases is illustrated in the right inset plot of Fig. 11(a) to emphasis the big impact of the atom-atom coupling on the entanglement value.
Increasing the atomic coupling further to 2 (solid blue line) then to 5 (dashed red line), the entanglement mean value increases considerably, as illustrated in Fig. 11(b). Now combining atomic coupling and non-zero detuning, λ 2 = 5, ∆ = 2 (dotted green line) the entanglement mean value decreases slightly, but higher detuning value ∆ = 5 decreases the mean value further and makes the entanglement between the rapid oscillation periods not constant any more (dash-dotted violet line).
The effect of the radiation intensity on the atomic population is shown in Fig. 11(c).
One can see that applying higher intensity, wheren = 20, 50 and finally 100 makes the collapse (constant) value approaches the zero value, also increases the collapse period and reduces the revival oscillation amplitude. Turning on the coupling, λ 2 = 5 at zero detuning  Fig. 11, one can see that the entire system starts at a separable state (the two atoms are in a completely disentangled state multiplied by the field coherent state), but the interaction between the atoms and the field, manifested as a rapid oscillation starting at τ = 0, triggers an entanglement rapid oscillation that eventually relaxes to a constant value that depends on the coupling strength and the detuning value.
The population revival oscillation is repeated periodically and is accompanied by a rapid entanglement oscillation that doesn't lead to a new entanglement value, particularly for zero or small ∆. This means in this particular case the energy exchange between the atoms and the fields do not cause entanglement transfer between the different subsystems. In Fig. 12, we show how the entanglement dynamics is very sensitive to changes in the coupling parameter λ 2 , where the ESD can be completely eliminated by increasing the coupling strength, whereas the atomic population is more sensitive to variations in ∆, where it changes considerably from negative to positive values as the detuning parameter is varied over a wide range from -5 to 5, as illustrated in (a) and (b) respectively.

V. CONCLUSION
We studied a system of two two-level atoms interacting with a single mode radiation field.
We introduced coupling between the two atoms and considered the radiation field to be out of resonance (at non-zero detuning) with the atoms. We presented an exact analytical solution for the time evolution of the system starting from any initial state. We investigated the effect of the atom-atom coupling and the non-zero detuning separately or combined (which has not been considered before) on the atom-atom entanglement dynamics and atomic population inversion, starting from different initial states of practical interest. We showed how these two parameters can be tuned to reduce, eliminate or create entanglement sudden death (ESD) in the system, which was found to depend crucially on the initial state of the system.
Particularly, we demonstrated that while one of the two interactions or both may have either negligible or weak impact on the ESD, combining them can be very effective in certain cases.
Starting from an initial correlated Bell state, the time evolution of entanglement between the two uncoupled atoms at resonance with the radiation field was found to suffer sudden death (ESD) for repeated intervals with revival oscillatory peaks in between. Turning on atom-atom coupling may reduce or even eliminate the ESD if applied at sufficient strength.
Although the non-zero detuning on its own does not affect the ESD for uncoupled atoms, it does considerably when combined with the coupling of atoms and contributes significantly in the ESD removal. For an initial anti-correlated Bell state, the entanglement evolved to intervals of very small constant value with intermediate revival peaks with no entanglement death observed for uncoupled atoms at resonance with the field. Nevertheless, applying a small value of detuning forces the entanglement to death intervals, but as the detuning was increased narrow intermediate peaks appeared. On the other hand, the atom-atom coupling at zero-detuning enhanced the entanglement between the two atoms considerably. Combining the two effects had competing impacts on the bipartite entanglement, where the non-zero detuning tended to create entanglement death while the atom-atom coupling acted to remove it. Starting from a partially entangled (W-like state), the bipartite entanglement evolved to sudden death in a very similar pattern to the correlated Bell state case but with smaller ESD intervals. However, in contrary to the Bell state, the non-zero detuning reduced or even completely eliminated the ESD, as its value was increased, leading to large entanglement values, whereas the atom-atom coupling enhanced the ESD and suppressed the entanglement revival oscillation. These competing effects sustain when both interactions are present in the system at the same time. For disentangled (separable) initial state, where both of the two atoms are in the excited state, the system never evolve to any ESD for any combination of system parameters values. Starting from that state, the entanglement showed oscillatory behavior where the non-zero detuning for uncoupled atoms raised the minima value of the entanglement oscillation and induced rapid oscillation within theses minima. Although the atom-atom coupling at zero-detuning has a negligible effect on the entanglement pattern, it enhances the entanglement considerably when applied at non-zero detuning. Finally, we considered an interesting initial separable state, which is a linear combination of all the basis states of the system. This initial state, for uncoupled atoms at zero detuning, was found to evolve to long intervals of ESD with very small intermediate spikes. Applying a non-zero detuning has no effect on the ESD, however turning the atom-atom coupling even with a small value at zero detuning completely eliminated the ESD and lead to a collapse-revival like pattern, which did not collapse to zero but a non-zero finite value that increased as the coupling strength was increased. When a non-zero detuning was applied to coupled atoms, it reduced the constant collapse value. By monitoring the atomic population inversion dynamics corresponding to these different initial states, a strong synchronization was observed in each case between the population collapse-revival pattern and the entanglement dynamics.
For all initial states that may evolve to ESD, for all system parameter combinations, the entanglement oscillatory revival peaks were found to be induced within the same intervals of the population revival peaks (where exchange of energy between the atoms and the radiation field takes place), whereas the sudden death (or the deviation from it) synchronized with the population collapse periods. In contrary, for the initial state that never evolve to any ESD under any system parameters combination, the atomic population revivals synchronized with the rapid oscillation that takes place at the minima of the entanglement oscillatory pattern, whereas the peaks of the entanglement oscillation synchronize with the population collapse periods. This means the exchange of the energy between the atoms and the fields in this case reduces the entanglement substantially. Varying the field radiation intensity showed a big impact on the ESD time intervals at all system parameters combinations. Increasing the field intensity increased the ESD period and most of the time reduced the entanglement revival oscillation amplitude. Therefore, while all the other system parameters cannot be used to significantly modify the collapse-revival temporal pattern, the field intensity can.
The ESD behavior observed in this closed system is due to a complete entanglement transfer from the atom-atom subsystem to the atom-field subsystems in absence of any decohering effects. The system shows entanglement revival peak after a finite time, due to an en-

DISCLOSURES
The authors declare no conflicts of interest.