Next Article in Journal
Reply to Pessoa, P.; Arderucio Costa, B. Comment on “Tsallis, C. Black Hole Entropy: A Closer Look. Entropy 2020, 22, 17”
Next Article in Special Issue
Local Information as an Essential Factor for Quantum Entanglement
Previous Article in Journal
Multi-Person Tracking and Crowd Behavior Detection via Particles Gradient Motion Descriptor and Improved Entropy Classifier
Previous Article in Special Issue
Mutual Information and Quantum Discord in Quantum State Discrimination with a Fixed Rate of Inconclusive Outcomes
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Asymptotic Entanglement Sudden Death in Two Atoms with Dipole–Dipole and Ising Interactions Coupled to a Radiation Field at Non-Zero Detuning

1
Department of Applied Physics and Astronomy, University of Sharjah, Sharjah 27272, United Arab Emirates
2
Department of Physics, Ain Shams University, Cairo 11566, Egypt
3
Department of Physics, Imam Muhammad Ibn Saud Islamic University, Riyadh 11432, Saudi Arabia
4
Department of Physics, University of Siegen, 57068 Siegen, Germany
5
Department of Physics, University of Hamburg, 22589 Hamburg, Germany
*
Author to whom correspondence should be addressed.
Submission received: 19 April 2021 / Revised: 11 May 2021 / Accepted: 13 May 2021 / Published: 18 May 2021
(This article belongs to the Collection Quantum Information)

Abstract

:
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.

1. 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]. Furthermore, 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 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. Furthermore, 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 Section 2, we discuss our model and present the exact analytic solution. In Section 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 Section 4. We conclude in Section 5.

2. 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 X X Z exchange interaction between two spin-1/2 particles with X X 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 by
H ^ = Ω a ^ a ^ + ω 2 i = 1 , 2 σ ^ z ( i ) + λ 1 i = 1 , 2 ( a ^ σ ^ + ( i ) + a ^ σ ^ ( i ) ) + λ 2 ( σ ^ ( 1 ) σ ^ + ( 2 ) + σ ^ + ( 1 ) σ ^ ( 2 ) ) + J σ z ( 1 ) σ z ( 2 ) .
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, a ^ and a ^ are creation and annihilation operators of the radiation field which satisfy the usual commutation relation [ a ^ , a ^ ] = 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
| ψ ( 0 ) = [ a | e 1 , e 2 + b | e 1 , g 2 + c | g 1 , e 2 + d | g 1 , g 2 ] | α ,
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
| α = n Q n | n ; Q n = α n n ! exp | α | 2 2 ,
where | α | 2 = n ¯ is the mean photon number and | n are the photon number states, which satisfy the relations: a ^ | n = n + 1 | n + 1 and a ^ | n + 1 = n + 1 | n . The wave function at any time t latter can be written as
| ψ ( t ) = n [ A n ( t ) | e 1 , e 2 , n + B n + 1 ( t ) | e 1 , g 2 , n + 1 + C n + 1 ( t ) | g 1 , e 2 , n + 1 + D n + 2 ( t ) | g 1 , g 2 , n + 2 ] ,
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 (Equation (1)) as
H ^ = H ^ + H ^ i n t ,
where
H ^ = Ω N ^ + Δ 2 i = 1 , 2 σ ^ z ( i ) ,
H ^ i n t = λ 1 i = 1 , 2 ( a ^ σ ^ + ( i ) + a ^ σ ^ ( i ) ) + λ 2 ( σ ^ ( 1 ) σ ^ + ( 2 ) + σ ^ + ( 1 ) σ ^ ( 2 ) ) + J σ z ( 1 ) σ z ( 2 ) ,
where Δ = ω Ω is the detuning parameter and N ^ = a ^ a ^ + 1 2 i = 1 , 2 σ ^ z ( i ) , 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 define V ^ I = U ^ H ^ i n t U ^ and U ^ = e i H ^ t . This yields
V ^ I ( t ) = λ 1 k = 1 , 2 ( a ^ e i Δ t σ ^ + ( k ) + a ^ e i Δ t σ ^ ( k ) ) + λ 2 ( σ ^ ( 1 ) σ ^ + ( 2 ) + σ ^ + ( 1 ) σ ^ ( 2 ) ) + J σ z ( 1 ) σ z ( 2 ) .
Now substituting | ψ ( t ) and V I ( t ) into Schrödinger equation
i t | ψ ( t ) = V ^ I ( t ) | ψ ( t ) ,
we get the following system of differential equations
i A ˙ n ( t ) = α e i Δ t ( B n + 1 ( t ) + C n + 1 ( t ) ) + J A n ( t ) , i B ˙ n + 1 ( t ) = α e i Δ t A n ( t ) + β e i Δ t D n + 2 ( t ) + λ 2 C n + 1 ( t ) J B n + 1 ( t ) , i C ˙ n + 1 ( t ) = α e i Δ t A n ( t ) + β e i Δ t D n + 2 ( t ) + λ 2 B n + 1 ( t ) J C n + 1 ( t ) , i D ˙ n + 2 ( t ) = β e i Δ t ( B n + 1 ( t ) + C n + 1 ( t ) ) + J D n + 2 ( t ) ,
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
A n ( t ) = A n ( 0 ) e i J t i α j = 1 3 [ δ j m j + i ( Δ + J ) ( e ( m j + i Δ ) t e i J t ) ] , B n + 1 ( t ) = 1 2 [ ( B n + 1 ( 0 ) C n + 1 ( 0 ) ) e i ( λ 2 + J ) t + j = 1 3 δ j e m j t ] , C n + 1 ( t ) = 1 2 [ ( C n + 1 ( 0 ) B n + 1 ( 0 ) ) e i ( λ 2 + J ) t + j = 1 3 δ j e m j t ] , D n + 2 ( t ) = D n + 2 ( 0 ) e i J t i β j = 1 3 [ δ j m j i ( J Δ ) ( e ( m j i Δ ) t e i J t ) ] ,
where
δ 1 = ( B n + 1 ( 0 ) + C n + 1 ( 0 ) ) ( δ 2 + δ 3 ) , δ 2 = 1 ( m 1 m 2 ) ( m 3 m 2 ) { 2 α A n ( 0 ) [ i ( m 1 + m 3 ) λ 2 Δ ] + 2 β D n + 2 ( 0 ) [ i ( m 1 + m 3 ) λ 2 + Δ ] + [ i ( m 1 + m 3 ) ( λ 2 J i m 1 ) 2 ( α 2 + β 2 ) ( λ 2 J ) 2 m 1 2 ] × ( B n + 1 ( 0 ) + C n + 1 ( 0 ) ) } , δ 3 = 1 ( m 1 m 3 ) ( m 2 m 3 ) { 2 α A n ( 0 ) [ i ( m 1 + m 2 ) λ 2 Δ ] + 2 β D n + 2 ( 0 ) [ i ( m 1 + m 2 ) λ 2 + Δ ] + [ i ( m 1 + m 2 ) ( λ 2 J i m 1 ) 2 ( α 2 + β 2 ) ( λ 2 J ) 2 m 1 2 ] × ( B n + 1 ( 0 ) + C n + 1 ( 0 ) ) } ,
and
m 1 = ( v 1 + v 2 ) i λ 2 + J 3 , m 2 = v 1 + v 2 2 + i 3 2 ( v 1 v 2 ) i λ 2 + J 3 , m 3 = v 1 + v 2 2 i 3 2 ( v 1 v 2 ) i λ 2 + J 3 ,
where
v 1 = [ μ 2 + ( μ 2 4 + η 3 27 ) 1 2 ] 1 3 ; v 2 = [ μ 2 ( μ 2 4 + η 3 27 ) 1 2 ] 1 3 ,
and
μ = 2 27 i [ 27 α 2 Δ + 27 β 2 Δ + 8 J 3 + 18 α 2 J + 18 β 2 J 18 Δ 2 J + λ 2 9 α 2 9 β 2 + 9 Δ 2 12 J 2 + 6 J λ 2 2 λ 2 3 ] ,
η = 2 α 2 + β 2 + Δ 2 + 1 3 ( λ 2 2 J ) 2 .
The initial values of the coefficients are as follows
A n ( 0 ) = Q n a , B n + 1 ( 0 ) = Q n + 1 b , C n + 1 ( 0 ) = Q n + 1 c , D n + 2 ( 0 ) = Q n + 2 d .
Now we can construct the system wave function | ψ ( t ) and calculate the system density matrix ρ ^ ( t ) = | ψ ( t ) ψ ( t ) | . The reduced density matrix of the two atoms, ρ ^ red ( t ) , can be obtained by tracing out the field
ρ ^ red ( t ) = Tr field ρ ^ ( t ) = l l | ψ ( t ) ψ ( t ) | l .

3. Time Evolution of Entanglement and Population Inversion

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 entanglement of formation E f through the formula E f ( ρ red ) = ε ( C ( ρ red ) ) , where ε is defined as
ε ( C ( ρ red ) ) = h 1 + 1 C 2 ( ρ red ) 2 ,
h is the Shannon entropy function h ( x ) = x log 2 x ( 1 x ) log 2 ( 1 x ) , and the concurrence can be calculated from C ( ρ red ) = max [ 0 , ε 1 ε 2 ε 3 ε 4 ] . 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 ρ ^ r e d (Equation (17)), we can evaluate σ ^ z ( t ) as
σ ^ z ( t ) = T r [ ρ ^ 1 ( t ) σ ^ z ] = n = 0 | A n ( t ) | 2 + | B n + 1 ( t ) | 2 | C n + 1 ( t ) | 2 | D n + 2 ( t ) | 2 .

3.1. Initial Bell State

In Figure 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 ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 with the radiation field is in a coherent state with field intensity corresponding to n ¯ = 100 . In the forthcoming discussion, we set n ¯ = 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 Figure 1a, we study the entanglement dynamics at resonance, i.e., zero detuning ( Δ = 0 ) , at different choices of λ 2 and J. As can be noticed, when the two parameters are (0, 0), i.e., non-interacting atoms, there is an ESD for a whole period of time from about 2 to 27 before reviving again, represented by the solid (blue) line. Turning on the Anti-Ferromagnetic (AF) Ising coupling, J > 0 , with a small value such as 1 shows a negligible effect on the ESD, but increasing that value to 2, shows a partial removal of the ESD, as shown clearly in the upper inner inset, dashed (green) line. When the Ferromagnetic (FM) Ising coupling, J < 0 , is considered instead, it shows the same exact behavior as the Ferromagnetic case for this initial state at λ 2 = 0. Turning on the dipole coupling λ 2 at J = 0 , i.e., for (2, 0), shows no effect on ESD at small values, such as 2 represented by the dot-dashed (red line), however applying higher values of λ 2 , removes the ESD partially and then eliminates it totally as it increases. Interestingly, while we could remove the ESD completely by setting the values at (4, 0), turning on AM Ising even at a small value is devastating for the entanglement, (4, 1) for instance, shows a complete ESD. However, applying an FM Ising proves to be an enhancement for the entanglement, (2, 2 ) removes ESD partially, represented by the (brown) solid line with dot marks, while (2, 4 ) (the (magenta) solid line with x marks) and (4, 4 ) (the (pink) solid line with star marks) remove the ESD completely and boost the entanglement value considerably. This means, at resonance, combining λ 2 and FM J is very useful in removing ESD and more effective than any of them separately, while combining λ 2 and AF J cancels out their individual effects and maintains the ESD.
In Figure 1b, we consider an out of resonance case, with positive detuning ( Δ = 2 ) . For non-interacting atoms, with parameter values (0,0), the ESD sustains over the entire time period (blue line) from about 2 to 27. However, turning on FM (or AF) Ising has a strong impact on the ESD compared with the Δ = 0 case, where setting the parameters at ( 0 , ± 2 ) removes completely the ESD with the AM entanglement peaks (red line) is slightly higher than that of the FM one (green line). Turning on the dipole interaction while keeping J = 0 , removes the ESD partially, for (2, 0), and for higher values of λ 2 , it completely eliminates the ESD and enhances the Entanglement. Adding the AF Ising to the dipole interaction, (2, 2), does not have a large impact on the entanglement behavior except for a slight shift (brown line). However, adding the FM Ising removes the ESD and enhances the entanglement significantly, even for parameter values as low as (2, 2 ) (magenta line). Increasing the value of (FM) J further increases the entanglement considerably (not shown here). The entanglement dynamics for negative detuning, Δ = 2 , illustrated in Figure 1c, 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 Figure 1d. While, as we concluded from Figure 1a–c, combining the dipole–dipole and FM Ising interactions is very effective in eliminating the ESD and enhancing the entanglement in the first period, it is clear from panel (d) that setting the parameter values even at (4, 4 ) can remove the ESD from the first and second periods but not the latter ones. In fact, the entanglement exhibits asymptotically an irregular oscillation interrupted with persistent ESD small periods, as shown in the upper inset of panel (d). Applying higher values of the parameters, such as (4, 8 ) (green line), may remove the ESD in more early periods, as shown in the lower inset of panel (d), however, it causes the entanglement to vanish entirely asymptotically.
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 Figure 2a,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 Figure 2a, but eventually causes the entanglement to vanish completely without revival as illustrated in Figure 2b, 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 ψ B a = ( | 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 in Figure 2c. Interestingly, as can be observed in Figure 2b,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 Figure 3, we illustrate the atomic population inversion for one of the two atoms, starting from an initial correlated Bell state ψ B c with the radiation field is in the coherent state. Under all different conditions, the atomic population shows the usual collapse-revival pattern that is a sign of atom-radiation field interaction. It collapses to a constant value that depends on the chosen set of parameters. In Figure 3a, we discuss the resonance case, where, as can be noticed, at the parameter values (0, 0), the collapse constant value (blue line) is around 2.45 × 10 3 but when J is turned on, ( 0 , ± 2 ) , the constant value shifts down (green line) and the collapse period slightly increases. On the other hand, for the values (2, 0) (which coincides with (2, 2)), the constant value (red line) gets closer to the (0, 0) case. Combining the dipole and FM Ising couplings, (2, 2 ) (purple line), shows a big shift away from the (0, 0) case, which agrees with our observations from the entanglement dynamics. By comparing Figure 1 and Figure 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.
In Figure 3b,c, we discuss the non-zero detuning cases, Δ = 2 , and 2 , respectively. 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. Furthermore, as we pointed in Figure 3a, the relative separation among the different collapse lines matches the corresponding entanglement values relative separation. Comparing the three different detuning cases, in panels (a), (b), and (c), it is clear that the off-resonance condition splits the collapse lines away from each other considerably compared with the resonance case. The asymptotic behavior of the atomic population is illustrated in Figure 3d at the parameter values (4, 4 ) (blue line) and (4, 8 ) (red line). As can be seen, the collapse periods get shorter with time before disappearing asymptotically and the general profile turns into an irregular continuous oscillation with varying reduced amplitude (as shown in the inset), which indicates a continuous energy exchange between the atoms and the field.

3.2. Partially Entangled Initial (W) State

In Figure 4, Figure 5 and Figure 6, we discuss the system dynamics starting from a partially entangled initial state that yields ESD upon evolution, namely, the W-like state ψ W = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 ) / 3 . The ESD periods that we used to observe in the maximally entangled initial states, for instance, the first period from τ 3 to 28, are interrupted, in the current case, in the middle by a single reviving peak at all choices of λ 2 and J, as shown in all panels of Figure 4.
The width of these peaks and their heights vary depending on the different system parameters. The maximum height corresponds to (0, 2), i.e., with only AF Ising coupling between the two atoms, at all Δ values, while the largest width takes place at Δ = 2 as shown in Figure 4b. The ESD can be removed partially, particularly in the first ESD period at the detuning Δ = 2 , mostly by combining the dipole and the FM Ising couplings as shown in the inner inset of Figure 4c. Interestingly, increasing the combined values of the dipole and FM Ising interactions further, (4, −4), as shown in Figure 4d, yields a complete elimination of the ESD at the different detuning values, most effectively at Δ = 2 and less effectively at Δ = 2 . The asymptotic behavior of the entanglement, illustrated in Figure 4d, shows that the entanglement collapse periods diminish with time until completely disappearing and the entanglement takes an oscillatory form without any ESD interruptions as illustrated in the inset in the bottom panel in Figure 4d. In Figure 5a, we plot the entanglement dynamics versus the Ising coupling J over a wide range from 5 to 5. As can be noticed, applying high values of Ferromagnetic Ising coupling eliminates the ESD, while having an anti-Ferromagnetic coupling leads to ESD, although it is accompanied by revivals, as shown in the figure. The same effect of the Ising coupling type and strength is confirmed in Figure 5b, where the contour plot of the asymptotic entanglement is presented.
The atomic population is illustrated in Figure 6a–c for the detuning values 0, 2 and 2 , respectively. Again, the separation of the collapse lines from the (0, 0) blue line, at all parameter values match the same one of the lines in the entanglement case, Figure 4. The collapse lines separation from each other increases from Δ = 0 to 2 and mostly at 2 . The detuning effect is clear in shifting the collapse lines upwards, for Δ = 2 , and downwards, for Δ = 2 , compared with the Δ = 0 case. Noticeably, the inverted reflection symmetry of the population lines, as we turn from Δ = 2 to Δ = 2 , that was observed in the Bell state case, in Figure 3, is not seen in the current case. In Figure 6d, 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.

3.3. 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 Figure 7a (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 a low FM Ising coupling, (0, 2 ), saves the entanglement from any sudden death and makes it maintain a non-zero value for the whole time, with a collapse-revival-like pattern (red line). Furthermore, turning on the dipole interaction in the absence of the Ising interaction, (2, 0), causes the entanglement to exhibit a very similar behavior to that of the setup (0, 2 ) but at a smaller mean value (purple line). Adding AF Ising interaction to the dipole one, (2, 2), turns out to be devastating to the entanglement, where the ESD is reached early right after the entanglement rises from zero (brown line), as illustrated in the inner inset. However, as we have observed in the previous cases, combining the dipole and FM Ising interactions is very effective in removing the ESD and boosting the entanglement, which is the case here too, where applying (2, 2 ) causes the entanglement to rise from zero to a high persisting mean value (magenta line), compared with the other choices. Remarkably, changing the detuning to a non-zero value ( ± 2 ) does not show a significant impact on the entanglement dynamics, as can be noticed in Figure 7b,c, except for quite small reviving peaks in 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 Figure 7d, where the asymptotic behavior of the entanglement shows the disappearance of the collapse periods and the emergence of a continuous irregular oscillation with a very similar profile at all the detuning values. Figure 8 shows the decisive role the Ising coupling plays when the system starts from the separable initial state ψ L , where the entanglement dynamics is tested versus the Ising coupling J. The ESD takes place at early times and asymptotically for positive values of J that are higher than approximately 0.95 , but can be completely eliminated for smaller positive values of J and all negative values.
In fact, the critical role played by the Ising coupling in controlling the system dynamics and asymptotic behavior is highlighted when one particularly compares the 3-dimensional plots in Figure 2, Figure 5 and Figure 8. They clearly show how the type of Ising coupling (Ferromagnetic or anti-Ferromagnetic) and its strength, depending on the initial state, may lead to completely different entanglement profiles. While Figure 2b illustrates how a strong Ising coupling (of any type) could be devastating for the asymptotic entanglement starting from an initial correlated Bell state, Figure 2c shows how the same coupling can considerably boost the asymptotic entanglement starting from an anti-correlated Bell state. On the other hand, Figure 5 shows how the effect of the Ising coupling type is crucial when the system starts from a partially entangled (W−like) state, where only a Ferromagnetic coupling is capable of eliminating the ESD at early times and asymptotically as illustrated in Figure 5a,b, 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 Figure 8.
The behavior of the population inversion differs significantly as the detuning parameter is varied, as shown in Figure 9. The zero detuning case is illustrated in Figure 9a, which exhibits a much smaller oscillation amplitude 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 Figure 9a,b. The positive detuning shifts the collapse line upwards, whereas the negative one shifts the lines downwards compared with the zero detuning case. Remarkably, similar to the maximally entangled initial state case and in contrary to the w-state one, the collapse–revival lines in the non-detuning cases are inverted images from each other. The asymptotic behavior of the atomic population, at (2, 2 ), is depicted in Figure 9d for the three different detuning cases. The profiles of the positive and negative detunings are very similar, as shown in the middle and bottom panels of Figure 9d, though the main values are positive and negative, respectively. 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 Figure 9d. 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 (orientation) of the spin with respect to the field. The precession (Larmor) frequency is determined by the magnetic field strength. On the other hand, when two spin-1/2 particles interact with each other, one of them could be considered as (a magnetic dipole) precessing in the field (magnetic moment) of the other one, where the dynamics of the system is 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.

4. Quantum Correlation between the Two Atoms and the Radiation 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 et al. derived an analytic form for the concurrence of a bipartite system with arbitrary dimensions d A 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
C ( ψ ) = 2 ν A ν B 1 T r ρ A 2
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
τ A , B ( ψ ) C 2 ( ψ ) = 2 ν A ν B 1 T r ρ A 2
In our case, the tangle takes the form
τ F A 1 A 2 = 2 1 T r ρ F 2 = 2 1 T r ρ A 1 A 2 2
In Figure 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 2 , shown in panels (a), (b) and (c), respectively, the tangle starts at zero but rises abruptly to a value 1 within the same time period where the entanglement and the population have their initial rapid oscillations. Then, the tangle decreases reaching a minimum in the middle of a time period that coincides with the population collapse and ESD periods. It peaks up again reaching a maximum value that is higher than the previous one, with a local minimum, within a time period that coincides with the entanglement and population revival oscillation periods. This behavior continues afterward as shown in the different panels, where the maxima and minima values increase with time and get closer while the time periods are shrinking in a similar fashion to the entanglement and population inversion profiles. The inner insets in panel (a) show that for Δ = 0 the curves coincide for ( 0 , ± 2 ) as well as for (2, 0) and (2, 2), but they split at a non-zero detuning, as illustrated in the inner insets of panels (b) and (c). The asymptotic behavior of the tangle is presented in Figure 10d, at Δ = 2 at two choices of the parameters; (4, 4 ) (blue line) and (4, 8 ) (red line). The tangle profile is oscillatory with a decreasing amplitude and increasing mean value that reach an approximately constant value asymptotically for both cases, which also agrees with the observed behavior of the entanglement and the population inversion. In fact, the other choices of the parameter values, at the different detuning values, show a similar asymptotic behavior, which we did not present here. In general, the different panels of Figure 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 Figure 11a,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 Figure 11c,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 parameter choices from each other.

5. Conclusions

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 of the system. The atomic population inversion and the quantum correlation between the two atoms ensemble and the radiation field synchronized with the entanglement dynamics, where the ESD periods coincided with the collapse ones in the population inversion and the valleys in the quantum correlations, while the revival oscillations of entanglement and atomic population matched the peaks in the quantum correlation. This means the entanglement revives from death and the quantum correlation maximizes when the two atoms exchange energy with the field. The asymptotic behavior of the system showed crucial dependence on the initial state, where the ESD was found to be always removable by tuning the system parameters except in the case of a maximally entangled initial correlated Bell state, where ESD was found to be persistent in that case for any choice of the parameters values. The entanglement, atomic population and quantum correlations were found to synchronize and reach asymptotically a quasi-steady dynamic states of continuous irregular oscillation with a limited amplitude and an approximately constant mean value, indicating a continuous exchange of energy between the two atoms ensemble and the field along with a strong correlation between the atoms and the field, and diminished ESD upon tuning the system parameters.

Author Contributions

Conceptualization, G.S. and W.A.-D.; Data curation, S.S. and H.E.; Formal analysis, G.S., W.A.-D., S.S. and H.E.; Investigation, G.S., W.A.-D., S.S. and H.E.; Methodology, G.S., W.A.-D.; Project administration, G.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by University of Sharjah, Office of Vice Chancellor of Research, grant No. 2002143094-P.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Jaynes, E.T.; Cummings, F.W. Comparison of quantum and semiclassical radiation theories with application to the beam maser. Proc. IEEE 1963, 51, 89–109. [Google Scholar] [CrossRef] [Green Version]
  2. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information, 1st ed.; Cambridge University Press: Cambridge, UK, 2010. [Google Scholar]
  3. Buluta, I.; Ashhab, S.; Nori, F. Natural and artificial atoms for quantum computation. Rep. Prog. Phys. 2011, 74, 104401. [Google Scholar] [CrossRef] [Green Version]
  4. Xiang, Z.-L.; Ashhab, S.; You, J.Q.; Nori, F. Hybrid quantum circuits: Superconducting circuits interacting with other quantum systems. Rev. Mod. Phys. 2013, 85, 623. [Google Scholar] [CrossRef] [Green Version]
  5. Lodahl, P.; Mahmoodian, S.; Stobbe, S. Interfacing single photons and single quantum dots with photonic nanostructures. Rev. Mod. Phys. 2015, 87, 347. [Google Scholar] [CrossRef] [Green Version]
  6. Wendin, G. Quantum information processing with superconducting circuits: A review. Rep. Prog. Phys. 2017, 80, 106001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Yang, C.-P.; Chu, S.-I.; Han, S. Possible realization of entanglement, logical gates, and quantum-information transfer with superconducting-quantum-interference-device qubits in cavity QED. Phys. Rev. A 2003, 67, 042311. [Google Scholar] [CrossRef] [Green Version]
  8. Yang, C.-P.; Chu, S.-I.; Han, S. Quantum information transfer and entanglement with SQUID qubits in cavity QED: A dark-state scheme with tolerance for nonuniform device parameter. Phys. Rev. Lett. 2004, 92, 117902. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Blais, A.; Huang, R.-S.; Wallraff, A.; Girvin, S.M.; Schoelkopf, R.J. Cavity quantum electrodynamics for superconducting electrical circuits: An architecture for quantum computation. Phys. Rev. A 2004, 69, 062320. [Google Scholar] [CrossRef] [Green Version]
  10. Yamamoto, T.; Pashkin, Y.A.; Astafiev, O.; Nakamura, Y.; Tsai, J.-S. Demonstration of conditional gate operation using superconducting charge qubits. Nature 2003, 425, 941–944. [Google Scholar] [CrossRef]
  11. Majer, J.B.; Paauw, F.G.; Ter Haar, A.C.J.; Harmans, C.J.P.M.; Mooij, J.E. Spectroscopy on two coupled superconducting flux qubits. Phys. Rev. Lett. 2005, 94, 090501. [Google Scholar] [CrossRef] [Green Version]
  12. Steffen, M.; Ansmann, M.; Bialczak, R.C.; Katz, N.; Lucero, E.; McDermott, R.; Neeley, M.; Weig, E.M.; Cleland, A.N.; Martinis, J.M. Measurement of the entanglement of two superconducting qubits via state tomography. Science 2006, 313, 1423–1425. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Van der Ploeg, S.H.W.; Izmalkov, A.; van den Brink, A.M.; Hübner, U.; Grajcar, M.; Il’Ichev, E.; Meyer, H.-G.; Zagoskin, A.M. Controllable coupling of superconducting flux qubits. Phys. Rev. Lett. 2007, 98, 057004. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Niskanen, A.O.; Harrabi, K.; Yoshihara, F.; Nakamura, Y.; Lloyd, S.; Tsai, J.S. Quantum coherent tunable coupling of superconducting qubits. Science 2007, 316, 723–726. [Google Scholar] [CrossRef] [PubMed]
  15. Osnaghi, S.; Bertet, P.; Auffeves, A.; Maioli, P.; Brune, M.; Raimond, J.-M.; Haroche, S. Coherent control of an atomic collision in a cavity. Phys. Rev. Lett. 2001, 87, 037902. [Google Scholar] [CrossRef] [Green Version]
  16. Gywat, O.; Meier, F.; Loss, D.; Awschalom, D.D. Dynamics of coupled qubits interacting with an off-resonant cavity. Phys. Rev. B 2006, 73, 125336. [Google Scholar] [CrossRef] [Green Version]
  17. Guerlin, C.; Brion, E.; Esslinger, T.; Mølmer, K. Cavity quantum electrodynamics with a Rydberg-blocked atomic ensemble. Phys. Rev. A 2010, 82, 053832. [Google Scholar] [CrossRef] [Green Version]
  18. Donaire, M.; Muñoz-Castañeda, J.M.; Nieto, L. Dipole-dipole interaction in cavity QED: The weak-coupling, nondegenerate regime. Phys. Rev. A 2017, 96, 042714. [Google Scholar] [CrossRef] [Green Version]
  19. Nguyen, T.L.; Raimond, J.-M.; Sayrin, C.; Cortinas, R.; Cantat-Moltrecht, T.; Assemat, F.; Dotsenko, I.; Gleyzes, S.; Haroche, S.; Roux, G.; et al. Towards quantum simulation with circular Rydberg atoms. Phys. Rev. X 2018, 8, 011032. [Google Scholar] [CrossRef] [Green Version]
  20. Sorensen, A.; Molmer, K. Spin-Spin Interaction and Spin Squeezing in an Optical Lattice. Phys. Rev. Lett. 1999, 83, 2274. [Google Scholar] [CrossRef] [Green Version]
  21. Porras, D.; Cirac, J.I. Effective Quantum Spin Systems with Trapped Ions. Phys. Rev. Lett. 2004, 92, 207901. [Google Scholar] [CrossRef] [Green Version]
  22. Cusati, T.; Napoli, A.; Messina, A. Competition between inter- and intra-molecular energy exchanges in a simple quantum model of a dimer. J. Mol. Struct. Theochem. 2006, 769, 3. [Google Scholar] [CrossRef]
  23. Napoli, A.; Messina, A.; Cusati, T.; Draganescu, G. Quantum signatures in the dynamics of two dipole-dipole interacting soft dimers. Eur. Phys. J. B 2006, 50, 419. [Google Scholar] [CrossRef]
  24. Hartmann, M.J.; Brandão, G.S.L.; Plenio, M.B. Effective Spin Systems in Coupled Microcavities. Phys. Rev. Lett. 2007, 99, 160501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Khlebnikov, S.; Sadiek, G. Decoherence by a nonlinear environment: Canonical versus microcanonical case. Phys. Rev. A 2002, 66, 032312. [Google Scholar] [CrossRef] [Green Version]
  26. Huang, Z.; Sadiek, G.; Kais, S. Time evolution of a single spin inhomogeneously coupled to an interacting spin environment. J. Chem. Phys. 2006, 124, 144513. [Google Scholar] [CrossRef] [PubMed]
  27. Abliz, A.; Gao, H.J.; Xie, X.C.; Wu, Y.S.; Liu, W.M. Entanglement control in an anisotropic two-qubit Heisenberg XYZ model with external magnetic fields. Phys. Rev. A 2006, 74, 052105. [Google Scholar] [CrossRef] [Green Version]
  28. Tsomokos, D.I.; Hartmann, M.J.; Huelga, S.F.; Plenio, M.B. Entanglement dynamics in chains of qubits with noise and disorder. New J. Phys. 2007, 9, 79. [Google Scholar] [CrossRef]
  29. Burić, N. Influence of the thermal environment on entanglement dynamics in small rings of qubits. Phys. Rev. A 2008, 77, 012321. [Google Scholar] [CrossRef] [Green Version]
  30. Dubi, Y.; Di Ventra, M. Relaxation times in an open interacting two-qubit system. Phys. Rev. A 2009, 79, 012328. [Google Scholar] [CrossRef] [Green Version]
  31. Sadiek, G.; Alkurtass, B.; Aldossary, O. Entanglement in a time-dependent coupled XY spin chain in an external magnetic field. Phys. Rev. A 2010, 82, 052337. [Google Scholar] [CrossRef] [Green Version]
  32. Xu, Q.; Sadiek, G.; Kais, S. Dynamics of entanglement in a two-dimensional spin system. Phys. Rev. A 2011, 83, 062312. [Google Scholar] [CrossRef] [Green Version]
  33. Sadiek, G.; Kais, S. Persistence of entanglement in thermal states of spin systems. J. Phys. B 2013, 46, 245501. [Google Scholar] [CrossRef] [Green Version]
  34. Duan, L.; Wang, H.; Chen, Q.-H.; Zhao, Y. Entanglement dynamics of two qubits coupled individually to Ohmic baths. J. Chem. Phys. 2013, 139, 044115. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wu, N.; Nanduri, A.; Rabitz, H. Rabi oscillations, decoherence, and disentanglement in a qubit–spin-bath system. Phys. Rev. A 2014, 89, 062105. [Google Scholar] [CrossRef] [Green Version]
  36. Sadiek, G.; Almalki, S. Entanglement dynamics in Heisenberg spin chains coupled to a dissipative environment at finite temperature. Phys. Rev. A 2016, 94, 012341. [Google Scholar] [CrossRef] [Green Version]
  37. Yu, T.; Eberly, J.H. Finite-time disentanglement via spontaneous emission. Phys. Rev. Lett. 2004, 93, 140404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Yu, T.; Eberly, J.H. Sudden death of entanglement. Science 2009, 323, 598–601. [Google Scholar] [CrossRef] [Green Version]
  39. Yu, T.; Eberly, J.H. Sudden death of entanglement: Classical noise effects. Opt. Commun. 2006, 264, 393–397. [Google Scholar] [CrossRef] [Green Version]
  40. Yönaç, M.; Yu, T.; Eberly, J.H. Sudden death of entanglement of two Jaynes–Cummings atoms. J. Phys. B 2006, 39, S621–S625. [Google Scholar] [CrossRef]
  41. Sainz, I.; Klimov, A.B.; Roa, L. Entanglement dynamics modified by an effective atomic environment. Phys. Rev. A 2006, 73, 032303. [Google Scholar] [CrossRef] [Green Version]
  42. Yönaç, M.; Yu, T.; Eberly, J.H. Pairwise concurrence dynamics: A four-qubit model. J. Phys. B 2007, 73, S45–S59. [Google Scholar] [CrossRef]
  43. Sainz, I.; Björk, G. Entanglement invariant for the double Jaynes-Cummings model. Phys. Rev. A 2007, 76, 042313. [Google Scholar] [CrossRef] [Green Version]
  44. Chan, S.; Reid, M.D.; Ficek, Z. Entanglement evolution of two remote and non-identical Jaynes–Cummings atoms. J. Phys. B 2009, 42, 065507. [Google Scholar] [CrossRef] [Green Version]
  45. Tavassoly, M.K.; Daneshmand, R.; Rustaee, N. Entanglement dynamics of linear and nonlinear interaction of two two-level atoms with a quantized Phase-damped field in the dispersive regime. Int. J. Theor. Phys. 2018, 57, 1645–1658. [Google Scholar] [CrossRef]
  46. Bai, X.-M.; Gao, C.-P.; Li, J.-Q.; Liu, N.; Liang, J.-Q. Entanglement dynamics for two spins in an optical cavity–field interaction induced decoherence and coherence revival. Opt. Express 2017, 25, 17051–17065. [Google Scholar] [CrossRef] [Green Version]
  47. Liu, R.-F.; Chen, C.-C. Role of the Bell singlet state in the suppression of disentanglement. Phys. Rev. A 2006, 74, 024102. [Google Scholar] [CrossRef] [Green Version]
  48. Ficek, Z.; Tanaś, R. Dark periods and revivals of entanglement in a two-qubit system. Phys. Rev. A 2006, 74, 024304. [Google Scholar] [CrossRef] [Green Version]
  49. Deng, T.; Yan, Y.; Chen, L.; Zhao, Y. Dynamics of the two-spin spin-boson model with a common bath. J. Chem. Phys. 2016, 144, 144102. [Google Scholar] [CrossRef]
  50. Li, C.; Xiao-Qiang, S.; Shou, Z. The influences of dipole–dipole interaction and detuning on the sudden death of entanglement between two atoms in the Tavis–Cummings model. Chin. Phys. B 2009, 18, 888–893. [Google Scholar] [CrossRef]
  51. Zhang, G.-F.; Chen, Z.-Y. The entanglement character between atoms in the non-degenerate two photons Tavis–Cummings model. Opt. Commun. 2007, 275, 274–277. [Google Scholar] [CrossRef] [Green Version]
  52. Torres, J.M.; Sadurní, E.; Seligman, T.H. Two interacting atoms in a cavity: Exact solutions, entanglement and decoherence. J. Phys. A 2010, 43, 192002. [Google Scholar] [CrossRef]
  53. de los Santos-Sánchez, O.; González-Gutiérrez, C.; Récamier, J. Nonlinear Jaynes–Cummings model for two interacting two-level atoms. J. Phys. B 2016, 49, 165503. [Google Scholar] [CrossRef] [Green Version]
  54. Chathavalappil, N.; Satyanarayana, S.V.M. Schemes to avoid entanglement sudden death of decohering two qubit system. Eur. Phys. J. D 2019, 73, 36. [Google Scholar] [CrossRef]
  55. Zhang, H.; Fan, C.; Wu, J. In-phase and anti-phase entanglement dynamics of Rydberg atomic pairs. Opt. Exp. 2020, 28, 35350–35362. [Google Scholar] [CrossRef] [PubMed]
  56. Sadiek, J.; Al-Dress, W.; Abdallah, M.S. Manipulating entanglement sudden death in two coupled two-level atoms interacting off-resonance with a radiation field: An exact treatment. Opt. Exp. 2019, 27, 33799–33825. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Wootters, W.K. Entanglement of formation of an arbitrary state of two qubits. Phys. Rev. Lett. 1998, 80, 2245–2248. [Google Scholar] [CrossRef] [Green Version]
  58. Rungta, P.; Buzek, V.; Caves, C.M.; Hillery, H.; Milburn, G.J. Universal state inversion and concurrence in arbitrary dimensions. Phys. Rev. A 2001, 64, 042315. [Google Scholar] [CrossRef] [Green Version]
Figure 1. The entanglement E f versus the scaled time τ = λ 1 t with the two atoms are initially in a correlated Bell state ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 and the field is in a coherent state with mean number of photons n ¯ = 100 , where the detuning parameter Δ takes the values (a) 0; (b) 2; (c) 2 ; (d) 2 , while the dipole–dipole and Ising coupling parameters ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 1. The entanglement E f versus the scaled time τ = λ 1 t with the two atoms are initially in a correlated Bell state ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 and the field is in a coherent state with mean number of photons n ¯ = 100 , where the detuning parameter Δ takes the values (a) 0; (b) 2; (c) 2 ; (d) 2 , while the dipole–dipole and Ising coupling parameters ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g001
Figure 2. The entanglement E f versus the scaled time τ = λ 1 t and the Ising coupling J with the two atoms are initially in the correlated Bell state ψ B c and the field is in a coherent state, where λ 2 = 2 and Δ = 2 , at early time in (a) and asymptotically (contour plot) in (b), and asymptotically (contour plot) for the initial anti-correlated Bell state ψ B a = ( | g 1 | e 2 + | e 1 | g 2 ) / 2 in (c).
Figure 2. The entanglement E f versus the scaled time τ = λ 1 t and the Ising coupling J with the two atoms are initially in the correlated Bell state ψ B c and the field is in a coherent state, where λ 2 = 2 and Δ = 2 , at early time in (a) and asymptotically (contour plot) in (b), and asymptotically (contour plot) for the initial anti-correlated Bell state ψ B a = ( | g 1 | e 2 + | e 1 | g 2 ) / 2 in (c).
Entropy 23 00629 g002
Figure 3. The population inversion σ z versus the scaled time τ = λ 1 t with the two atoms are initially in a correlated Bell state ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 2 , while ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 3. The population inversion σ z versus the scaled time τ = λ 1 t with the two atoms are initially in a correlated Bell state ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 2 , while ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g003
Figure 4. The entanglement E f versus the scaled time τ = λ 1 t with the two atoms are initially in a W-like state ψ W = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 ) / 3 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 4. The entanglement E f versus the scaled time τ = λ 1 t with the two atoms are initially in a W-like state ψ W = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 ) / 3 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g004
Figure 5. The entanglement E f versus the scaled time τ = λ 1 t and the Ising coupling J with the two atoms are initially in a W-like state ψ W and the field is in a coherent state where λ 2 = 4 and Δ = 2 , at early time in (a) and asymptotically at (b).
Figure 5. The entanglement E f versus the scaled time τ = λ 1 t and the Ising coupling J with the two atoms are initially in a W-like state ψ W and the field is in a coherent state where λ 2 = 4 and Δ = 2 , at early time in (a) and asymptotically at (b).
Entropy 23 00629 g005
Figure 6. The population inversion σ z versus the scaled time τ = λ 1 t with the two atoms are initially in a W-like state ψ W = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 ) / 3 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 6. The population inversion σ z versus the scaled time τ = λ 1 t with the two atoms are initially in a W-like state ψ W = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 ) / 3 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g006
Figure 7. The entanglement E f versus the scaled time τ = λ 1 t with the two atoms are initially in a separable state ψ L = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 + | e 1 | e 2 ) / 4 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 7. The entanglement E f versus the scaled time τ = λ 1 t with the two atoms are initially in a separable state ψ L = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 + | e 1 | e 2 ) / 4 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g007
Figure 8. The entanglement E f versus the scaled time τ = λ 1 t and the Ising coupling J with the two atoms are initially in the separable state ψ L and the field is in a coherent state, where λ 2 = 2 and Δ = 2 , at early time in (a) and asymptotically at (b).
Figure 8. The entanglement E f versus the scaled time τ = λ 1 t and the Ising coupling J with the two atoms are initially in the separable state ψ L and the field is in a coherent state, where λ 2 = 2 and Δ = 2 , at early time in (a) and asymptotically at (b).
Entropy 23 00629 g008
Figure 9. The population inversion σ z versus the scaled time τ = λ 1 t with the two atoms are initially in a separable state ψ L = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 + | e 1 | e 2 ) / 4 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 9. The population inversion σ z versus the scaled time τ = λ 1 t with the two atoms are initially in a separable state ψ L = ( | g 1 | g 2 + | g 1 | e 2 + | e 1 | g 2 + | e 1 | e 2 ) / 4 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 0 , 2 and 2 . ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g009
Figure 10. The tangle versus the scaled time τ = λ 1 t with the two atoms are initially in a correlated Bell state ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 2 , while ( λ 2 , J ) take various values as shown in the legend in each panel.
Figure 10. The tangle versus the scaled time τ = λ 1 t with the two atoms are initially in a correlated Bell state ψ B c = ( | e 1 | e 2 + | g 1 | g 2 ) / 2 and the field is in a coherent state with n ¯ = 100 , where in (a) Δ = 0 ; (b) 2; (c) 2 ; (d) 2 , while ( λ 2 , J ) take various values as shown in the legend in each panel.
Entropy 23 00629 g010
Figure 11. The tangle versus the scaled time τ = λ 1 t with the field is in a coherent state with n ¯ = 100 and the two atoms are initially in the W-like state ψ W in (a,b), and in the separable state ψ L in (c,d), at different values of the detuning as indicated in every panel. The values of the parameters ( λ 2 , J ) in panels (a,c) are as shown in the legend of Figure 10a.
Figure 11. The tangle versus the scaled time τ = λ 1 t with the field is in a coherent state with n ¯ = 100 and the two atoms are initially in the W-like state ψ W in (a,b), and in the separable state ψ L in (c,d), at different values of the detuning as indicated in every panel. The values of the parameters ( λ 2 , J ) in panels (a,c) are as shown in the legend of Figure 10a.
Entropy 23 00629 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sadiek, G.; Al-Dress, W.; Shaglel, S.; Elhag, H. Asymptotic Entanglement Sudden Death in Two Atoms with Dipole–Dipole and Ising Interactions Coupled to a Radiation Field at Non-Zero Detuning. Entropy 2021, 23, 629. https://0-doi-org.brum.beds.ac.uk/10.3390/e23050629

AMA Style

Sadiek G, Al-Dress W, Shaglel S, Elhag H. Asymptotic Entanglement Sudden Death in Two Atoms with Dipole–Dipole and Ising Interactions Coupled to a Radiation Field at Non-Zero Detuning. Entropy. 2021; 23(5):629. https://0-doi-org.brum.beds.ac.uk/10.3390/e23050629

Chicago/Turabian Style

Sadiek, Gehad, Wiam Al-Dress, Salwa Shaglel, and Hala Elhag. 2021. "Asymptotic Entanglement Sudden Death in Two Atoms with Dipole–Dipole and Ising Interactions Coupled to a Radiation Field at Non-Zero Detuning" Entropy 23, no. 5: 629. https://0-doi-org.brum.beds.ac.uk/10.3390/e23050629

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop