Next Article in Journal
Comparing Information-Theoretic Measures of Complexity in Boltzmann Machines
Next Article in Special Issue
Optomechanical Analogy for Toy Cosmology with Quantized Scale Factor
Previous Article in Journal
Non-Causal Computation
Previous Article in Special Issue
Boltzmann Sampling by Degenerate Optical Parametric Oscillator Network for Structure-Based Virtual Screening
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamics of Entanglement in Jaynes–Cummings Nodes with Nonidentical Qubit-Field Coupling Strengths

Fujian Key Laboratory of Quantum Information and Quantum Optics, College of Physics and Information Engineering, Fuzhou University, Fujian 350116, China
*
Author to whom correspondence should be addressed.
Entropy 2017, 19(7), 331; https://doi.org/10.3390/e19070331
Submission received: 2 June 2017 / Revised: 18 June 2017 / Accepted: 29 June 2017 / Published: 3 July 2017
(This article belongs to the Collection Quantum Information)

Abstract

:
How to analytically deal with the general entanglement dynamics of separate Jaynes–Cummings nodes with continuous-variable fields is still an open question, and few analytical approaches can be used to solve their general entanglement dynamics. Entanglement dynamics between two separate Jaynes–Cummings nodes are examined in this article. Both vacuum state and coherent state in the initial fields are considered through the numerical and analytical methods. The gap between two nonidentical qubit-field coupling strengths shifts the revival period and changes the revival amplitude of two-qubit entanglement. For vacuum-state fields, the maximal entanglement is fully revived after a gap-dependence period, within which the entanglement nonsmoothly decreases to zero and partly recovers without exhibiting sudden death phenomenon. For strong coherent-state fields, the two-qubit entanglement decays exponentially as the evolution time increases, exhibiting sudden death phenomenon, and the increasing gap accelerates the revival period and amplitude decay of the entanglement, where the numerical and analytical results have an excellent coincidence.

1. Introduction

Generation and preservation of the quantum entanglement between the qubit and field have played an important role for both fundamental quantum theories and experiments [1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37], and they are still a challenge for real applications [38,39,40,41,42,43,44,45,46,47,48], such as quantum key distribution and teleportation. In theory, it is feasible to generate and preserve quantum entanglement based on the well-known Jaynes–Cummings model [49,50,51,52,53,54,55,56]. There have some physical systems dealing with quantum correlation and entanglement through autocorrelation function recently [57,58,59,60], such as semiconductor microcavities [59] and optomechanics [60].
However, since there is not a general measure for multi-body quantum systems so far [61,62], only two-particle entanglement is definitely quantified and the entanglement dynamics becomes hard to analytically handle as the dimension of qubit-field systems increases. Previous study [49] has shown that the entanglement sudden death and rebirth appear in two separate nodes and each node is analytically described by the Jaynes–Cummings model Hamiltonian, where two initial fields are both in the vacuum state, which are very difficult to generate and preserve in real experiments.
Coherent-state field, as a kind of continuous-variable physical system, contains infinite eigenstate spectrums and can be efficiently generated by a classical monochromatic current [63], which is an easily feasible experiment resource and most resembles a classical electromagnetic field [64,65], resisting for decoherence induced by the environment. These excellent features have made coherent-state field be widely used in many fields of quantum information processing recently, such as quantum transport [66,67] and storage [68,69,70]. However, coherent-state field is hard to obtain the analytical time-dependent dynamics when it is coupled to the qubit due to its infinite-dimensional Hilbert space. It has been theoretically reported [69] that one-ebit entanglement reciprocation between two qubits and two coherent-state fields is feasible based on numerical results. Although the time-dependent eigenfunction of separate Jaynes–Cummings nodes with coherent-state fields is known by numerical diagonalization in a truncated Hilbert space, its analytical solutions is quite necessary for clearly capturing and experimentally controlling the fundamental entanglement physics.
As far as we know, how to analytically deal with the general entanglement dynamics of separate Jaynes–Cummings nodes with continuous-variable fields is still an open question, and few analytical approaches can be directly used to solve their general entanglement dynamics. Recent study [71,72] uses a saddle point method to show that the two-qubit entanglement dynamics of two identical Jaynes–Cummings nodes can be analytically predicted by an exponentially decaying formula when the amplitudes of the coherent-state fields are both large enough. However, due to analytical diagonalization obstacle in the asymmetric infinite-dimension Hilbert space, the entanglement dynamics of two separate Jaynes–Cummings nodes with nonidentical qubit-field coupling constants have not been extensively studied [46] but is common in real experiments.
Based on the previous study [72], we focus here on the analytic entanglement dynamics between two independent, separate standard Jaynes–Cummings models, where the node is the qubit inside the cavity field and the qubit-field coupling constants for the local subsystem are different. Both vacuum-state and coherent-state fields are considered through the numerical and analytical methods. By using the saddle point method for coherent-state fields, the numerical and analytical results have an excellent coincidence. Our method is suitable for physical systems with ignored dissipation or without dissipation, such as the cavity losses and atomic spontaneous emissions. We find that the gap between two qubit-field coupling strengths shifts the revival period and changes the revival amplitude of two-qubit entanglement. For vacuum-state fields, the maximal entanglement is fully revival after a gap-dependence period, within which the entanglement nonsmoothly decreases to zero and partly recovers without exhibiting sudden death phenomenon. For strong coherent-state fields, the two-qubit entanglement decays exponentially as the evolution time increases, exhibiting sudden death phenomenon, and the increasing gap accelerates the revival period and amplitude decay of the entanglement. Our result demonstrates that when the average photon number is large enough, the decay exponent has a quadratic dependence on the qubit-field coupling strengths, and the non-full revival period is linearly shifted by the cooperative qubit-field coupling strength.

2. Separate Jaynes–Cummings Nodes

Two separate Jaynes–Cummings nodes with nonidentical qubit-field coupling strengths are described by the Hamiltonian ( = 1 )
H = ω 0 2 σ z A + ( G A σ + A a + G A * σ A a ) + ω a a + ω 0 2 σ z B + ( G B σ + B b + G B * σ B b ) + ω b b ,
where ω 0 is the transition frequency between the high level | e x and the low level | g x of the qubit x ( x = A , B ). σ z x and σ ± x are the general Pauli matrices of the qubit x. a (a) and b (b) are the creation (annihilation) operators for two single-mode fields with angular frequency ω , respectively. G x is the coupling strength between the qubit x and its localized field, and G x * is the conjugate complex of G x . The assumption G A G B , referred to as two nonidentical qubit-field coupling constants, represents a clear distinction from two standard Jaynes–Cummings nodes with G A = G B [72], which has been predicted to exhibit the exponential decay of two-qubit entanglement and its revival timing and duration by an approximately analytic formula when the fields are nearly classical. For simplicity, we further assume that there is not any interaction between two qubits or two fields and the qubit resonantly couples to its localized field, i.e., ω 0 = ω , as illustrated in Figure 1.
Since the dynamics of single Jaynes–Cummings node has been experimentally realized for different quantum systems [2], such as cavity quantum electrodynamics (QED) and circuit QED systems, the physical realization of two Jaynes–Cummings nodes with different qubit-field coupling strengths is certainly feasible when the parameters G A and G B are not perfectly identical, which is more common than that of two standard Jaynes–Cummings nodes with identical qubit-field coupling strengths in realistic experiments. Considering the reality in the experiment, it is not necessary to require that two sites are identical and here we focus on the elementary influences caused by two nonidentical qubit-field coupling strengths for two-qubit entanglement.
To analytically explain two-qubit entanglement for a little excited and a highly excited (nearly classical) field modes, we start from a simple case where the field modes are initially in their vacuum states and the two qubits are in the maximally entangled states. Since there is not any interaction between two nodes, the time-dependent evolution of the whole Jaynes–Cummings Hamiltonian can be analytically wrote out. Eigenvector evolutions of each node are the well-known Rabi oscillations [3], described by
| ψ e , n ( t ) = e i H t | e x , n y = cos ( G x t n + 1 ) | e x , n y i sin ( G x t n + 1 ) | g x , n + 1 y
and
| ψ g , n ( t ) = e i H t | g x , n y = cos ( G x t n ) | g x , n y i sin ( G x t n ) | e x , n 1 y ,
where the field mode’s photon state is denoted as the Fock state | n y (n is a positive integer) in the field mode y ( y = a , b ) and t is the evolution time.
Our target is to analytically measure the entanglement between two nonlocal qubits obtained from the time-dependent four-body state, so we adopt the Wootters concurrence C [73] as the entanglement measure for two-particle states expressed in the standard qubit basis Γ = { | e A e B , | e A g B , | g A e B , | g A g B } , which is defined as
C = max 0 , λ 1 λ 2 λ 3 λ 4 ,
where λ 1 , λ 2 , λ 3 , and λ 4 are the eigenvalues arranged in decreasing order of the following matrix:
ξ = ρ ( σ y σ y ) ρ * ( σ y σ y ) ,
where ρ is the two-qubit reduced density matrix and σ y is the Pauli matrix. C = 0 is for two unentangled qubits and C = 1 stands for the maximally entangled state for two qubits.

2.1. Vacuum-State Field

When the field modes are initially in their vacuum states and the two qubits are initially in the maximally entangled states, i.e.,
| Ψ ( 0 ) = 1 2 ( | e A g B + | g A e B ) | 0 a , 0 b ,
the evolution of the whole qubit-field system is
| Ψ ( t ) = e i H t 2 ( | e A g B + | g A e B ) | 0 a , 0 b = 1 2 [ cos ( G A t ) | e A , g B ; 0 a , 0 b i sin ( G A t ) | g A , g B ; 1 a , 0 b + cos ( G B t ) | g A , e B ; 0 a , 0 b i sin ( G B t ) | g A , g B ; 0 a , 1 b ] .
By tracing out two field modes within Γ , the resulting two-qubit mixed state has the following X form [74]
ρ = T r a , b ( | Ψ ( t ) Ψ ( t ) | ) = 0 0 0 0 0 ρ 22 ρ 23 0 0 ρ 32 ρ 33 0 0 0 0 ρ 44 ,
where ρ 22 = 1 2 cos 2 ( G A t ) , ρ 23 = ρ 32 = 1 2 cos ( G A t ) cos ( G B t ) , ρ 33 = 1 2 cos 2 ( G B t ) , and ρ 44 = 1 2 cos 2 ( G A t ) + 1 2 sin 2 ( G B t ) . The concurrence of this mixed state is analytically found to be
C = 2 max { 0 , | ρ 23 | , ρ 22 ρ 33 } = 1 2 | cos ( G A + G B ) t + cos ( G A G B ) t | .
When G A G B , the analytical period of C in Equation (9) is
T c = π | G A 2 G B 2 | LCM | G A + G B | , | G A G B | ,
where the symbol LCM ( ) represents the least common multiple among the numbers within its bracket. While when G A = G B , the period of C is π / G B . Figure 2 plots the time-dependent evolution of two-qubit concurrence for different values of two qubit-field coupling strengths. From Figure 2, we see that the maximum of C reaches 1, meaning that the maximal entanglement fully recovers after a evolution period T c even when the qubit-field coupling strengths are different. Between the oscillation periods, the entanglement nonsmoothly decreases to zero and partly recovers without staying zero for a finite interval of time, i.e., without exhibiting sudden death [49]. This concurrence shows that when the qubit-field coupling strengths are nonidentical, i.e., G A G B , its oscillation amplitude and period nonlinearly depend on the cooperative interaction terms G A + G B and G A G B , which is very different from two identical standard Jaynes–Cummings nodes where the concurrence exhibits a standard Rabi oscillation having a fixed period without depending on the qubit-field coupling strength [49].

2.2. Coherent-State Field

When the field modes are initially in the identical coherent states and the two qubits are initially in the maximally entangled states, i.e.,
| Ψ ( 0 ) = 1 2 ( | e A g B + | g A e B ) | α a , α b ,
where the coherent states are expanded by the Fock states
| α y = n = 0 A n | n y = n = 0 e | α | 2 / 2 α n n ! | n y .
Therefore, the evolution dynamics can be expressed as
| Ψ ( t ) = e i H t 2 ( | e A g B + | g A e B ) | α a , α b = 1 2 n = 0 m = 0 A n A m K n m ,
where
K n m = i C n + 1 A S m B | e A , e B ; n a , m 1 b + C n + 1 A C m B | e A , g B ; n a , m b S n + 1 A S m B | g A , e B ; n + 1 a , m 1 b i S n + 1 A C m B | g A , g B ; n + 1 a , m b i S n A C m + 1 B | e A , e B ; n 1 a , m b S n A S m + 1 B | e A , g B ; n 1 a , m + 1 b + C n A C m + 1 B | g A , e B ; n a , m b i C n A S m + 1 B | g A , g B ; n a , m + 1 b ,
in which C n x = cos ( G x t n ) and S n x = sin ( G x t n ) . Since the qubit-field coupling strengths are nonidentical in the infinite-dimension Hilbert space, the analytical two-qubit entanglement C is hard to obtain under a general condition. However, it is highly desirable to get analytical results for understanding the entanglement’s fundamental physics. In the following, we use the analytical formula of saddle point method [72] to tackle the entanglement dynamics of Equation (13) when the coherent states are both highly excited (nearly classical).
When the average photon number n ¯ > > 1 , it is reasonable to assume that photon number in the coherent state is a Poisson distribution that tightly centers around n ¯ , and photon number without centering around n ¯ can be ignored. By approximately replacing | α y with | n ¯ y and tracing out two field modes for Equation (13), we obtain the approximate X-form reduced density matrix of two qubits within Γ
ρ ρ 11 0 0 0 0 ρ 22 ρ 23 0 0 ρ 23 * ρ 33 0 0 0 0 ρ 44 ,
where ρ 23 * is the conjugate complex of ρ 23 . The concurrence of this ρ is
C = 2 max { 0 , | ρ 23 | ρ 11 ρ 44 } .
Based on the time-dependent eigenstate evolutions in Equations (2) and (3) of each node, it is workable to analytically figure out the matrix elements ρ 23 , ρ 11 , and ρ 44 in Equation (16) by tracing out two field modes ρ = T r a , b ( | Ψ ( t ) Ψ ( t ) | ) . The analytical expressions of doubly infinite summations for ρ 23 , ρ 11 , and ρ 44 are expressed as
ρ 23 = 1 2 n = 0 m = 0 ( A n 2 A m 2 C n + 1 A C n A C m + 1 B C m B + A n A n 2 A m A m + 2 S n A S n 1 A S m + 1 B S m + 2 B A n A n 1 A m A m + 1 S n A C n + 1 A S m + 1 B C m B A n A n 1 A m A m + 1 S n A C n 1 A S m + 1 B C m + 2 B ) ,
ρ 11 = 1 2 n = 0 m = 0 [ A n A m C n + 1 A S m B 2 + A n A m S n A C m + 1 B 2 + A n A n + 1 A m A m 1 C n + 1 A S n + 1 A S m B C m B + A n A n 1 A m A m + 1 S n A C n A C m + 1 B S m + 1 B ] ,
and
ρ 44 = 1 2 n = 0 m = 0 [ A n A m S n + 1 A C m B 2 + A n A m C n A S m + 1 B 2 + A n A n + 1 A m A m 1 S n + 1 A C n + 1 A S m B C m B + A n A n 1 A m A m + 1 S n A C n A S m + 1 B C m + 1 B ] .
These summation expressions cannot be analytically completed for the general average photon number. However, when n ¯ α y 2 > > 100 , it is feasible to introduce an error-deviation order 1 / n ¯ centering around the Poisson peak n m n ¯ and the terms A n A n ± 1 A m A m ± 1 can be approximated with A n 2 A m 2 . Therefore, we obtain the following approximate expressions
ρ 23 1 2 n = 0 A n 2 C n A C n + 1 A × n = 0 A n 2 C n B C n + 1 B + n = 0 A n 2 S n A S n + 1 A × n = 0 A n 2 S n B S n + 1 B n = 0 m = 0 A n 2 A m 2 S n A C n + 1 A C m B S m + 1 B ,
ρ 11 1 2 n = 0 A n C n A 2 × n = 0 A n S n B 2 + n = 0 A n S n A 2 × n = 0 A n C n B 2 + n = 0 A n 2 S n A C n A × n = 0 A n 2 S n B C n B ,
and
ρ 44 1 2 n = 0 A n S n A 2 × n = 0 A n C n B 2 + n = 0 A n C n A 2 × n = 0 A n S n B 2 + n = 0 A n 2 S n A C n A × n = 0 A n 2 S n B C n B .
It is easy to find that the approximate expressions satisfy ρ 11 ρ 44 under the condition of large average photon numbers. Based on the two-times angle cosine formula
C n x C n + 1 x = 1 2 cos [ G x t ( n + n + 1 ) ] + cos [ G x t ( n + 1 n ) ] ,
and the second-order Taylor expansion of large n ¯
n + 1 n + 1 2 n ,
the expression C n x C n + 1 x can be compressed to
C n x C n + 1 x 1 2 cos G x t 2 n + cos 2 G x t n .
Analogously, the other useful approximations are
C n x C n x = 1 2 [ 1 + cos ( 2 G x t n ) ] ,
S n x S n x = 1 2 [ 1 cos ( 2 G x t n ) ] ,
S n x S n + 1 x 1 2 cos G x t 2 n cos 2 G x t n ,
S n x C n + 1 x 1 2 sin 2 G x t n sin G x t 2 n ,
S n + 1 x C n x 1 2 sin 2 G x t n + sin G x t 2 n .
With the above approximations, ρ 23 , ρ 11 , and ρ 44 can be further simplified as
ρ 23 1 4 { [ n = 0 A n 2 cos ( 2 G A t n ) ] [ n = 0 A n 2 cos ( 2 G B t n ) ] + [ n = 0 A n 2 cos ( G A t 2 n ) ] [ n = 0 A n 2 cos ( G B t 2 n ) ] n = 0 A n 2 [ sin ( 2 G A t n ) sin ( G A t 2 n ) ] × n = 0 A n 2 [ sin ( 2 G B t n ) sin ( G B t 2 n ) ] } ,
and
ρ 11 ρ 44 1 4 { 1 [ n = 0 A n 2 cos ( 2 G A t n ) ] [ n = 0 A n 2 cos ( 2 G B t n ) ] + [ n = 0 A n 2 sin ( 2 G A t n ) ] [ n = 0 A n 2 sin ( 2 G B t n ) ] } .
The above results Equations (31) and (32) lead to
ρ 23 ρ 11 ρ 44 1 4 { [ n = 0 A n 2 cos ( G A t 2 n ) ] [ n = 0 A n 2 cos ( G B t 2 n ) ] [ n = 0 A n 2 sin ( 2 G A t n ) ] [ n = 0 A n 2 sin ( G B t 2 n ) ] + 2 [ n = 0 A n 2 cos ( 2 G A t n ) ] [ n = 0 A n 2 cos ( 2 G B t n ) ] 2 [ n = 0 A n 2 sin ( 2 G A t n ) ] [ n = 0 A n 2 sin ( 2 G B t n ) ] + [ n = 0 A n 2 sin ( G A t 2 n ) ] [ n = 0 A n 2 sin ( 2 G B t n ) ] + [ n = 0 A n 2 sin ( G A t 2 n ) ] [ n = 0 A n 2 sin ( G B t 2 n ) ] 1 } .
To analytically calculate Equation (33), we need to rewrite the sums into integrals by considering discrete n as continuations when n ¯ is large enough. Based on the integral results demonstrated in Ref. [72]
n = 0 A n 2 cos ( G x t 2 n ) 0 A n 2 cos ( G x t 2 n ) d n e G x 2 t 2 32 α 4 cos ( G x t 2 α ) ,
n = 0 A n 2 sin ( G x t 2 n ) 0 A n 2 sin ( G x t 2 n ) d n e G x 2 t 2 32 α 4 sin ( G x t 2 α ) ,
n = 0 A n 2 cos ( 2 G x t n ) 0 A n 2 cos ( 2 G x t n ) d n e G x 2 t 2 2 cos ( 2 α G x t ) + k = 1 1 π k e ( G x t 2 k π α ) 2 1 + π 2 k 2 cos [ 2 α ( G x t 2 k π α ) ] ,
n = 0 A n 2 sin ( 2 G x t n ) 0 A n 2 sin ( 2 G x t n ) d n e G x 2 t 2 2 sin ( 2 α G x t ) + k = 1 ( 1 ) k 1 π k e ( G x t 2 k π α ) 2 1 + π 2 k 2 sin [ 2 α ( G x t 2 k π α ) ] ,
we can finally obtain the effective formula
ρ 23 ρ 11 ρ 44 1 4 { e ( G A 2 + G B 2 ) t 2 32 α 4 cos [ ( G A G B ) t 2 α ] + 2 e ( G A 2 + G B 2 ) t 2 2 cos [ 2 α ( G A + G B ) t ] + 2 e G A 2 t 2 2 k = 1 1 k π e ( G B t 2 k π α ) 2 1 + k 2 π 2 cos [ 2 α G A t + ( 1 ) k ( 2 α G B t 4 k π α 2 ) ] + 2 e G B 2 t 2 2 k = 1 1 k π e ( G A t 2 k π α ) 2 1 + k 2 π 2 cos [ 2 α G B t + ( 1 ) k ( 2 α G A t 4 k π α 2 ) ] + 2 k 1 = 1 k 2 = 1 1 π 1 k 1 k 2 e ( G A t 2 k 1 π α ) 2 ( 1 + k 1 2 π 2 ) ( G B t 2 k 2 π α ) 2 ( 1 + k 2 2 π 2 ) cos [ ( 1 ) k 1 ( 2 α G A t 4 k 1 π α 2 ) + ( 1 ) k 2 ( 2 α G B t 4 k 2 π α 2 ) ] 1 } .
It is necessary to emphasize that only the terms of G x t = 2 k π α , 2 k 1 π α , and 2 k 2 π α with the corresponding k, k 1 , and k 2 in Equation (38) give a significant contribution to the sums, and the other terms’ contribution to the sums is proportional to the exponential functions in Equation (38). Thus, the sums decay exponentially with the difference from k, k 1 , and k 2 , which are the main results for the evolution of two-qubit concurrence. This result demonstrates that, when the average photon number is large enough, the decay exponent has a quadratic dependence on G A and G B , and the non-full revival period is linearly shifted by the cooperative coupling strength G A ± G B , leading to the relative revival envelope heights
2 1 k π + 1 π 1 k 1 k 2 1 e ( G A 2 + G B 2 ) t 2 32 α 4 2 .
This concurrence evolution of two Jaynes–Cummings nodes reduces to that of two identical standard Jaynes–Cummings nodes at both the decay exponent and revival period. In Figure 3, we plot the long-time system dynamics under the coherent states with the large average photon number for analytical and numerical calculations. We find that the two-qubit entanglement exhibits sudden death phenomenon, and its peaks are not fully revived and decrease quadraticly as the gap δ between two qubit-field coupling strengths enlarges. The two-qubit entanglement decays exponentially as the evolution time increases. When G A > G B , the revival period is linearly shifted to the “left” side, i.e., becoming shorter; while when G A < G B , it is linearly shifted to the “right” side, i.e., becoming longer. Note that the numerical results in Figure 3b are not perfectly predicted by the analytical results in Figure 3a, and their main difference is the absence of Rabi-type oscillations during the revivals, i.e., the disappearance of tiny revivals in numerical results, which is not predicted by this analytical method [72]. Therefore, it is safe to say that the analytical results predict the numerical results for small gaps between two qubit-field coupling strengths in both the entanglement revival peek and period. Our result is physically important because the formula reveals analytically the relation between the entanglement evolution dynamics and two nonidentical qubit-field coupling strengths, which further clarifies the mechanism of entanglement sudden death and rebirth and provides another effective direction for entanglement control.
To quantitatively show the periodic modulations by the amplitude of coherent states and two qubit-field coupling strengths, we simulate the evolution dynamics with two other average photon numbers n ¯ = 25 and n ¯ = 36 for analytical and numerical results in Figure 4. We see that even when the average photon number n ¯ decreases to 25, the analytical results can predict well the numerical results under small ratios δ . Figure 5 further shows that, for the same α , peeks of the first revival envelope P C depend quadraticly on the ratio δ , and their periods T P C linearly decrease as the ratio δ increases, which mathematically fits
2 α G B t π + 28 ( 1 δ ) 5 .
This result demonstrates that the bigger gap between two qubit-field coupling strengths becomes, the smaller and faster entanglement revival will be, meaning that the increasing gap between two qubit-field coupling strengths accelerates the revival period and amplitude decay of two-qubit entanglement, which also can be explained from the decaying contribution of ratio factor 1 + δ 2 within the exponents and 1 + δ within the cosine functions of Equation (38). Note that system asymmetry can be caused not only by the difference of the coupling strengths, but also by the difference of the amplitudes α x and α y . Based on numerical simulations, it is easy to find that the difference of the amplitudes α x and α y causes the similar influence on the two-qubit entanglement with that caused by the difference of the coupling strengths, i.e., the gap between two amplitudes affects both the period and amplitude of the entanglement revival. However, our method is hard to generalize to analytically give out the relation formula between this gap and entanglement, and this kind of interesting formula may be found elsewhere by a more effective method.

3. Conclusions

In conclusion, we have generalized the method of Ref. [72] to the model of two separate Jaynes–Cummings nodes, where the qubit-field coupling strength for each node is different. We use both the numerical and analytical methods to cases with vacuum-state and coherent-state fields, respectively. The numerical and analytical results have an excellent coincidence in the case of coherent-state fields for small gaps between two qubit-field coupling strengths. We find that the gap shifts the revival period and changes the revival amplitude of two-qubit entanglement. For vacuum-state fields, the maximal entanglement is fully revived after a gap-dependence period, within which the entanglement nonsmoothly decreases to zero and partly recovers without exhibiting a sudden death phenomenon. For strong coherent-state fields, the two-qubit entanglement decays exponentially as the evolution time increases, exhibiting sudden death phenomenon, and the increasing gap accelerates the revival period and amplitude decay of the entanglement. Our result finally demonstrates that when the average photon number is large enough, the decay exponent has a quadratic dependence on two qubit-field coupling strengths, and the non-full revival period is linearly shifted by the cooperative qubit-field coupling strength. Potential applications of our result, such as control of entanglement through changing the system parameters, are feasible in many simple experiments based on the fundamental Jaynes–Cummings Hamiltonian. In the future, we want to study further the large frequency detuning case and try to find a possible way to recover the full two-qubit entanglement under separate continuous-variable fields. Up to now, this model does not consider any dissipation, such as the cavity losses and atomic spontaneous emissions in extreme ultraviolet (XUV) and X-rays [75,76], which is an unavoidable decoherence problem and open question for keeping long-time entanglement in real experiments, and this kind of investigation could be done elsewhere by finding another effective method in the future.

Acknowledgments

This work is supported by the National Natural Science Foundation of China under Grant No. 11374054, the Natural Science Foundation of Fujian Province under Grant No. 2017J05005, and the fund from Fuzhou University under Nos. XRC-1566 and XRC-1639.

Author Contributions

All authors contributed equally to this article. The theoretical aspect of this work does not allow to point out any substantial authorship. All authors have read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nielsen, M.A.; Chuang, I.L. Quantum Computation and Quantum Information; Cambridge University Press: Cambridge, UK, 2000. [Google Scholar]
  2. Raimond, J.M.; Brune, M.; Haroche, S. Manipulating quantum entanglement with atoms and photons in a cavity. Rev. Mod. Phys. 2001, 73, 565–582. [Google Scholar] [CrossRef]
  3. 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]
  4. Zheng, S.B.; Guo, G.C. Efficient scheme for two-atom entanglement and quantum information processing in cavity QED. Phys. Rev. Lett. 2000, 85, 2392–2395. [Google Scholar] [CrossRef] [PubMed]
  5. Van Enk, S.J.; Cirac, J.I.; Zoller, P. Ideal quantum communication over noisy channels: A quantum optical implementation. Phys. Rev. Lett. 1997, 78, 4293–4296. [Google Scholar] [CrossRef]
  6. Ekert, A.K. Quantum cryptography based on Bell’s theorem. Phys. Rev. Lett. 1991, 67, 661–663. [Google Scholar] [CrossRef] [PubMed]
  7. Buluta, I.; Ashhab, S.; Nori, F. Natural and artificial atoms for quantum computation. Rep. Prog. Phys. 2011, 74, 104401. [Google Scholar] [CrossRef]
  8. Ashhab, S.; Niskanen, A.O.; Harrabi, K.; Nakamura, Y.; Picot, T.; de Groot, P.C.; Harmans, C.J.P.M.; Mooij, J.E.; Nori, F. Interqubit coupling mediated by a high-excitation-energy quantum object. Phys. Rev. B 2008, 77, 014510. [Google Scholar] [CrossRef]
  9. Rivas, A.; Huelga, S.F.; Plenio, M.B. Quantum non-Markovianity: Characterization, quantification and detection. Rep. Prog. Phys. 2014, 77, 094001. [Google Scholar] [CrossRef] [PubMed]
  10. Mortezapour, A.; Borji, M.A.; lo Franco, R. Protecting entanglement by adjusting the velocities of moving qubits inside non-Markovian environments. Laser Phys. Lett. 2017, 14, 055201. [Google Scholar] [CrossRef]
  11. Mortezapour, A.; Borji, M.A.; lo Franco, R. Non-Markovianity and coherence of a moving qubit inside a leaky cavity. arXiv, 2017; arXiv:1705.00887. [Google Scholar]
  12. Lo Franco, R.; d’Arrigo, A.; Falci, G.; Compagno, G.; Paladino, E. Spin-echo entanglement protection from random telegraph noise. Phys. Scr. 2013, 153, 014043. [Google Scholar] [CrossRef]
  13. Aolita, L.; de Melo, F.; Davidovich, L. Open-system dynamics of entanglement: A key issues review. Rep. Prog. Phys. 2015, 78, 042001. [Google Scholar] [CrossRef] [PubMed]
  14. Lo Franco, R.; Bellomo, B.; Maniscalco, S.; Compagno, G. Dynamics of quantum correlations in two-qubit systems within non-Markovian environments. Int. J. Mod. Phys. B 2013, 27, 1345053. [Google Scholar] [CrossRef]
  15. De Vega, I.; Alonso, D. Dynamics of non-Markovian open quantum systems. Rev. Mod. Phys. 2017, 89, 015001. [Google Scholar] [CrossRef]
  16. Lo Franco, R.; Compagno, G. Overview on the phenomenon of two-qubit entanglement revivals in classical environments. arXiv, 2016; arXiv:1608.05970. [Google Scholar]
  17. D’Arrigo, A.; lo Franco, R.; Benenti, G.; Paladino, E.; Falci, G. Recovering entanglement by local operations. Ann. Phys. 2014, 350, 211–224. [Google Scholar] [CrossRef]
  18. Leggio, B.; lo Franco, R.; Soares-Pinto, D.O.; Horodecki, P.; Compagno, G. Distributed correlations and information flows within a hybrid multipartite quantum-classical system. Phys. Rev. A 2015, 92, 032311. [Google Scholar] [CrossRef]
  19. González-Gutiérrez, C.A.; Román-Ancheyta, R.; Espitia, D.; lo Franco, R. Relations between entanglement and purity in non-Markovian dynamics. Int. J. Quantum Inform. 2016, 14, 1650031. [Google Scholar] [CrossRef]
  20. Sciara, S.; lo Franco, R.; Compagno, G. Universality of Schmidt decomposition and particle identity. Sci. Rep. 2017, 7, 44675. [Google Scholar] [CrossRef] [PubMed]
  21. Bellomo, B.; lo Franco, R.; Compagno, G. N identical particles and one particle to entangle them all. arXiv, 2017; arXiv:1704.06359. [Google Scholar]
  22. Lo Franco, R.; d’Arrigo, A.; Falci, G.; Compagno, G.; Paladino, E. Preserving entanglement and nonlocality in solid-state qubits by dynamical decoupling. Phys. Rev. B 2014, 90, 054304. [Google Scholar] [CrossRef]
  23. D’Arrigo, A.; Benenti, G.; lo Franco, R.; Falci, G.; Paladino, E. Hidden entanglement, system-environment information flow and non-Markovianity. Int. J. Quantum Inform. 2014, 12, 1461005. [Google Scholar] [CrossRef]
  24. Trapani, J.; Bina, M.; Maniscalco, S.; Paris, M.G. Collapse and revival of quantum coherence for a harmonic oscillator interacting with a classical fluctuating environment. Phys. Rev. A 2015, 91, 022113. [Google Scholar] [CrossRef]
  25. D’Arrigo, A.; lo Franco, R.; Benenti, G.; Paladino, E.; Falci, G. Hidden entanglement in the presence of random telegraph dephasing noise. Phys. Scr. 2013, 153, 014014. [Google Scholar] [CrossRef]
  26. Lo Franco, R.; d’Arrigo, A.; Falci, G.; Compagno, G.; Paladino, E. Entanglement dynamics in superconducting qubits affected by local bistable impurities. Phys. Scr. 2012, 147, 014019. [Google Scholar] [CrossRef]
  27. Bellomo, B.; Compagno, G.; lo Franco, R.; Ridolfo, A.; Savasta, S. Entanglement dynamics of two independent cavity-embedded quantum dots. Phys. Scr. 2011, 143, 014004. [Google Scholar] [CrossRef]
  28. Dijkstra, A.G.; Tanimura, Y. Non-Markovian entanglement dynamics in the presence of system-bath coherence. Phys. Rev. Lett. 2010, 104, 250401. [Google Scholar] [CrossRef] [PubMed]
  29. Lo Franco, R.; Compagno, G. Quantum entanglement of identical particles by standard information-theoretic notions. Sci. Rep. 2016, 6, 20603. [Google Scholar] [CrossRef] [PubMed]
  30. Bellomo, B.; lo Franco, R.; Andersson, E.; Cresser, J.D.; Compagno, G. Dynamics of correlations due to a phase-noisy laser. Sci. Rep. 2012, 147, 014004. [Google Scholar] [CrossRef]
  31. Maniscalco, S.; Francica, F.; Zaffino, R.L.; Gullo, N.L.; Plastina, F. Protecting entanglement via the quantum Zeno effect. Phys. Rev. Lett. 2008, 100, 090503. [Google Scholar] [CrossRef] [PubMed]
  32. Bellomo, B.; lo Franco, R.; Maniscalco, S.; Compagno, G. Two-qubit entanglement dynamics for two different non-Markovian environments. Phys. Scr. 2010, T140, 014014. [Google Scholar] [CrossRef]
  33. Bellomo, B.; Compagno, G.; lo Franco, R.; Ridolfo, A.; Savasta, S. Dynamics and extraction of quantum discord in a multipartite open system. Int. J. Quantum Inform. 2011, 9, 1665–1676. [Google Scholar] [CrossRef]
  34. Lo Franco, R. Nonlocality threshold for entanglement under general dephasing evolutions: A case study. Quantum Inform. Process. 2016, 15, 2393–2404. [Google Scholar] [CrossRef]
  35. Silva, I.A.; Souza, A.M.; Bromley, T.R.; Cianciaruso, M.; Marx, R.; Sarthour, R.S.; Oliveira, I.S.; lo Franco, R.; Glaser, S.J.; Deazevedo, E.R.; et al. Observation of time-invariant coherence in a nuclear magnetic resonance quantum simulator. Phys. Rev. Lett. 2016, 117, 160402. [Google Scholar] [CrossRef] [PubMed]
  36. Orieux, A.; D’Arrigo, A.; Ferranti, G.; lo Franco, R.; Benenti, G.; Paladino, E.; Falci, G.; Sciarrino, F.; Mataloni, P. Experimental on-demand recovery of entanglement by local operations within non-Markovian dynamics. Sci. Rep. 2015, 5, 8575. [Google Scholar] [CrossRef] [PubMed]
  37. Man, Z.X.; Xia, Y.J.; lo Franco, R. Cavity-based architecture to preserve quantum coherence and entanglement. Sci. Rep. 2015, 5, 13843. [Google Scholar] [CrossRef] [PubMed]
  38. Rendell, R.W.; Rajagopal, A.K. Revivals and entanglement from initially entangled mixed states of a damped Jaynes–Cummings model. Phys. Rev. A 2003, 67, 062110. [Google Scholar] [CrossRef]
  39. Rempe, G.; Walther, H.; Klein, N. Observation of quantum collapse and revival in a one-atom maser. Phys. Rev. Lett. 1987, 58, 353–356. [Google Scholar] [CrossRef] [PubMed]
  40. Brune, M.; Schmidt-Kaler, F.; Maali, A.; Dreyer, J.; Hagley, E.; Raimond, J.M.; Haroche, S. Quantum Rabi oscillation: A direct test of field quantization in a cavity. Phys. Rev. Lett. 1996, 76, 1800–1803. [Google Scholar] [CrossRef] [PubMed]
  41. Meekhof, D.M.; Monroe, C.; King, B.E.; Itano, W.M.; Wineland, D.J. Generation of nonclassical motional states of a trapped atom. Phys. Rev. Lett. 1996, 76, 1796–1799. [Google Scholar] [CrossRef] [PubMed]
  42. Boca, A.; Miller, R.; Birnbaum, K.M.; Boozer, A.D.; McKeever, J.; Kimble, H.J. Observation of the vacuum Rabi spectrum for one trapped atom. Phys. Rev. Lett. 2004, 93, 233603. [Google Scholar] [CrossRef] [PubMed]
  43. Gea-Banacloche, J. Collapse and revival of the state vector in the Jaynes–Cummings model: An example of state preparation by a quantum apparatus. Phys. Rev. Lett. 1990, 65, 3385–3388. [Google Scholar] [CrossRef] [PubMed]
  44. Phoenix, S.J.D.; Knight, P.L. Establishment of an entangled atom-field state in the Jaynes–Cummings model. Phys. Rev. A 1991, 44, 6023–6029. [Google Scholar] [CrossRef] [PubMed]
  45. Sainz, I.; Björk, G. Quantum error correction may delay, but also cause, entanglement sudden death. Phys. Rev. A 2008, 77, 052307. [Google Scholar] [CrossRef]
  46. 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]
  47. Agarwal, G.S. Vacuum-field Rabi splittings in microwave absorption by Rydberg atoms in a cavity. Phys. Rev. Lett. 1984, 53, 1732–1734. [Google Scholar] [CrossRef]
  48. Shen, L.T.; Chen, R.X.; Wu, H.Z.; Yang, Z.B. Distributed manipulation of two-qubit entanglement with coupled continuous variables. J. Opt. Soc. Am. B 2015, 32, 297–302. [Google Scholar] [CrossRef]
  49. Yu, T.; Eberly, J.H. Sudden death of entanglement. Science 2009, 323, 598–601. [Google Scholar] [CrossRef] [PubMed]
  50. Kraus, B.; Cirac, J.I. Discrete entanglement distribution with squeezed light. Phys. Rev. Lett. 2004, 92, 013602. [Google Scholar] [CrossRef] [PubMed]
  51. Sornborger, A.T.; Cleland, A.N.; Geller, M.R. Superconducting phase qubit coupled to a nanomechanical resonator: Beyond the rotating-wave approximation. Phys. Rev. A 2004, 70, 052315. [Google Scholar] [CrossRef]
  52. Paternostro, M.; Son, W.; Kim, M.S.; Falci, G.; Palma, G.M. Dynamical entanglement transfer for quantum-information networks. Phys. Rev. A 2004, 70, 022320. [Google Scholar] [CrossRef]
  53. Paternostro, M.; Son, W.; Kim, M.S. Complete conditions for entanglement transfer. Phys. Rev. Lett. 2004, 92, 197901. [Google Scholar] [CrossRef] [PubMed]
  54. Zhou, L.; Yang, G.H. Entanglement reciprocation between atomic qubits and an entangled coherent state. J. Phys. B 2006, 39, 5143–5150. [Google Scholar] [CrossRef]
  55. Hu, Y.H.; Tan, Y.G. Delaying disappearance of entanglement by the linear modulation of atom-field coupling. Phys. Scr. 2014, 89, 075103. [Google Scholar] [CrossRef]
  56. Xu, J.S.; Li, C.F.; Gong, M.; Zou, X.B.; Shi, C.H.; Chen, G.; Guo, G.C. Experimental demonstration of photonic entanglement collapse and revival. Phys. Rev. Lett. 2010, 104, 100502. [Google Scholar] [CrossRef] [PubMed]
  57. Eleuch, H.; Bennaceur, R. An optical soliton pair among absorbing three-level atoms. J. Opt. A 2003, 5, 528–533. [Google Scholar] [CrossRef]
  58. Eleuch, H. Autocorrelation function of microcavity-emitting field in the linear regime. Eur. Phys. J. D 2008, 48, 139–143. [Google Scholar] [CrossRef]
  59. Eleuch, H. Entanglement and autocorrelation function in semiconductor microcavities. Int. J. Mod. Phys. B 2010, 24, 5653–5662. [Google Scholar] [CrossRef]
  60. Sete, E.A.; Eleuch, H.; Ooi, C.R. Light-to-matter entanglement transfer in optomechanics. J. Opt. Soc. Am. B 2014, 31, 2821–2828. [Google Scholar] [CrossRef]
  61. Rafsanjani, S.M.H.; Eberly, J.H. Coherent control of multipartite entanglement. Phys. Rev. A 2015, 91, 012313. [Google Scholar] [CrossRef]
  62. Lawrence, J. Mermin inequalities for perfect correlations in many-qutrit systems. Phys. Rev. A 2017, 95, 042123. [Google Scholar] [CrossRef]
  63. Scully, M.O.; Zubairy, M.S. Quantum Optics; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  64. Braunstein, S.L.; van Loock, P. Quantum information with continuous variables. Rev. Mod. Phys. 2005, 77, 513–577. [Google Scholar] [CrossRef]
  65. Adesso, G.; Illuminati, F. Entanglement in continuous-variable systems: Recent advances and current perspectives. J. Phys. A 2007, 40, 7821–7880. [Google Scholar] [CrossRef]
  66. Casagrande, F.; Lulli, A.; Paris, M.G. Improving the entanglement transfer from continuous-variable systems to localized qubits using non-Gaussian states. Phys. Rev. A 2007, 75, 032336. [Google Scholar] [CrossRef]
  67. Blanco, P.; Mundarain, D. Faithful entanglement transference from qubits to continuous variable systems. J. Phys. B 2011, 44, 105501. [Google Scholar] [CrossRef]
  68. Ballester, D. Entanglement accumulation, retrieval, and concentration in cavity QED. Phys. Rev. A 2009, 79, 062317. [Google Scholar] [CrossRef]
  69. Lee, J.; Paternostro, M.; Kim, M.S.; Bose, S. Entanglement reciprocation between qubits and continuous variables. Phys. Rev. Lett. 2006, 96, 080501. [Google Scholar] [CrossRef] [PubMed]
  70. Paternostro, M.; Kim, M.S.; Palma, G.M. Accumulation of entanglement in a continuous variable memory. Phys. Rev. Lett. 2007, 98, 140504. [Google Scholar] [CrossRef] [PubMed]
  71. Yönaç, M.; Eberly, J.H. Qubit entanglement driven by remote optical fields. Opt. Lett. 2008, 33, 270–272. [Google Scholar] [CrossRef] [PubMed]
  72. Yönaç, M.; Eberly, J.H. Coherent-state control of noninteracting quantum entanglement. Phys. Rev. A 2010, 82, 022321. [Google Scholar] [CrossRef]
  73. Wootters, W.K. Entanglement of formation of an arbitrary state of two qubits. Phys. Rev. Lett. 1998, 80, 2245. [Google Scholar] [CrossRef]
  74. Yu, T.; Eberly, J.H. Evolution from entanglement to decoherence of bipartite mixed X states. Quantum Inf. Comput. 2007, 7, 459–468. [Google Scholar]
  75. Sete, E.A.; Svidzinsky, A.A.; Eleuch, H.; Yang, Z.; Nevels, R.D.; Scully, M.O. Correlated spontaneous emission on the Danube. J. Mod. Opt. 2010, 57, 1311–1330. [Google Scholar] [CrossRef]
  76. Sete, E.A.; Svidzinsky, A.A.; Rostovtsev, Y.V.; Eleuch, H.; Jha, P.K.; Suckewer, S.; Scully, M.O. Using Quantum Coherence to Generate Gain in the XUV and X-ray: Gain-Swept Superradiance and Lasing without Inversion. IEEE J. Sel. Top. Quantum Electron. 2012, 18, 541–553. [Google Scholar] [CrossRef]
Figure 1. Sketch illustrating two separate Jaynes–Cummings nodes. Qubits A and B are independently placed in Nodes 1 and 2, resonantly coupling to the fields a and b, respectively. There is not any interaction between A and B or between a and b.
Figure 1. Sketch illustrating two separate Jaynes–Cummings nodes. Qubits A and B are independently placed in Nodes 1 and 2, resonantly coupling to the fields a and b, respectively. There is not any interaction between A and B or between a and b.
Entropy 19 00331 g001
Figure 2. Time-dependent evolution of C for: (a) G A = 0.9 G B ; (b) G A = G B ; and (c) G A = 1.1 G B . The corresponding periods are 10 π / G B , π / G B , and 10 π / G B , respectively.
Figure 2. Time-dependent evolution of C for: (a) G A = 0.9 G B ; (b) G A = G B ; and (c) G A = 1.1 G B . The corresponding periods are 10 π / G B , π / G B , and 10 π / G B , respectively.
Entropy 19 00331 g002
Figure 3. Two-qubit concurrence as a function of the evolution time under different ratios δ = G A / G B with n ¯ = 100 : (a) analytical results and (b) numerical results, where analytical results are based on the X-form ρ in Equation (15) and numerical results are based on the original ρ = T r a , b ( | Ψ ( t ) Ψ ( t ) | ) .
Figure 3. Two-qubit concurrence as a function of the evolution time under different ratios δ = G A / G B with n ¯ = 100 : (a) analytical results and (b) numerical results, where analytical results are based on the X-form ρ in Equation (15) and numerical results are based on the original ρ = T r a , b ( | Ψ ( t ) Ψ ( t ) | ) .
Entropy 19 00331 g003
Figure 4. Two-qubit concurrence evolves with time under different ratios δ = G A / G B for analytical results ((a) n ¯ = 25 and (b) n ¯ = 36 ) and numerical results ((c) n ¯ = 25 and (d) n ¯ = 36 ).
Figure 4. Two-qubit concurrence evolves with time under different ratios δ = G A / G B for analytical results ((a) n ¯ = 25 and (b) n ¯ = 36 ) and numerical results ((c) n ¯ = 25 and (d) n ¯ = 36 ).
Entropy 19 00331 g004
Figure 5. Characteristics of the first revival envelope versus the ratio δ for different α : (a) peek P C and (b) period T P C .
Figure 5. Characteristics of the first revival envelope versus the ratio δ for different α : (a) peek P C and (b) period T P C .
Entropy 19 00331 g005

Share and Cite

MDPI and ACS Style

Shen, L.-T.; Shi, Z.-C.; Wu, H.-Z.; Yang, Z.-B. Dynamics of Entanglement in Jaynes–Cummings Nodes with Nonidentical Qubit-Field Coupling Strengths. Entropy 2017, 19, 331. https://doi.org/10.3390/e19070331

AMA Style

Shen L-T, Shi Z-C, Wu H-Z, Yang Z-B. Dynamics of Entanglement in Jaynes–Cummings Nodes with Nonidentical Qubit-Field Coupling Strengths. Entropy. 2017; 19(7):331. https://doi.org/10.3390/e19070331

Chicago/Turabian Style

Shen, Li-Tuo, Zhi-Cheng Shi, Huai-Zhi Wu, and Zhen-Biao Yang. 2017. "Dynamics of Entanglement in Jaynes–Cummings Nodes with Nonidentical Qubit-Field Coupling Strengths" Entropy 19, no. 7: 331. https://doi.org/10.3390/e19070331

APA Style

Shen, L. -T., Shi, Z. -C., Wu, H. -Z., & Yang, Z. -B. (2017). Dynamics of Entanglement in Jaynes–Cummings Nodes with Nonidentical Qubit-Field Coupling Strengths. Entropy, 19(7), 331. https://doi.org/10.3390/e19070331

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