Quantum Correlations Between Identical and Unidentical Atoms in a Dissipative Environment
We have studied the dynamics of quantum correlations such as entanglement, Bell-nonlocality and quantum discord between identical as well as unidentical atoms interacting with a single-mode cavity field and subject to cavity decay. The effect of single atom detuning, cavity decay rate and initial preparation of the atoms on the corresponding correlation measures have been investigated. It is found that even under strong dissipation, time evolution can create high quantum discord while entanglement and Bell nonlocality stay zero for an initially separable state. Quantum discord increases while entanglement decreases in a certain time period under dissipation for the initial state that both atoms are in the excited state if the qubits are identical. For some type of initial states, cavity decay is shown to drive the system to a stationary state with high entanglement and quantum discord.
pacs:03.65.Ud, 03.65.Yz, 03.67.-a
Entanglement, nonlocality and quantum discord (QD) are all different but somehow related aspects of the quantumness of the correlated composite quantum systems. Entanglement is a kind of quantum correlation and determines whether the given state is separable or not . On the other hand, the quantum systems may contain quantum correlations other than entanglement. Quantum discord, which is defined as the difference between the quantum versions of two analogous classically equivalent definitions of mutual information, captures all types of nonclassical correlations including entanglement . Although entanglement and QD are equivalent for pure states , QD is found to be nonzero for some separable states [3, 4]. Such states were shown to be useful in the deterministic quantum computations with one pure qubit context [5, 6]. Therefore QD is believed a new resource for quantum computation. The nonlocality as measured by the violation of Bell’s inequalities is also a signature of the inseparability of a quantum state and signifies entangled states whose correlations cannot be reproduced by a classical local hidden-variable model . However it must be emphasized here that the presence of entanglement does not always imply violation of Bell’s inequalities, while the violation of Bell inequalities implies entanglement. Such entangled states violating Bell’s inequalities play a central role in some applications in quantum information science , such as to guarantee the safety of device-independent key distrubition protocols in quantum cryptography [9, 10].
The relations among different measures of quantum correlations and their dynamical behaviour have been an active area of research with an aim to understand the fundamental question of what quantumness is as well as to characterize these correlations as useful resources for various quantum tasks [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]. One of the most important obstacle to realizing quantum operations is the interaction between the qubits and the environment which tends to wipe out the quantumness of the system . Characterization of the effect of environment on the system such as decoherence and dissipation is important. The dynamics of quantum correlations such as entanglement, quantum discord and Bell nonlocality for a two-qubit system subject to various environments have been studied by many groups recently [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]. Entanglement sudden death (ESD) is one of the new concepts that arise within the context of these studies . ESD refers to death of entanglement between two qubits in a finite time while the coherence of the single qubit decays only exponentially. In contrary to entanglement, QD is more robust in a decohering environment, actually it was shown that almost all states have non-zero quantum discord . QD was shown to decay exponentially in Markovian environments, even at finite temperatures [14, 16] and can become zero momentarily under non-Markovian dynamics . QD was also experimentally investigated in an all-optical set up recently . Furthermore, QD and entanglement were analyzed for a Heisenberg spin chain with quantum phase transitions and QD was found to detect the critical points of the transition at finite temperatures while entanglement could not . Very recently, Sun et al. investigated the induced quantum correlations between two qubits coupled to a common environment modelled by Ising spin chain in a transverse field and during the time evolution, non-zero quantum discord is found to be created for specific initial states whereas entanglement cannot be induced . In Ref.  quantum discord is found to be totally unaffected by non-dissipative noise for a long time interval under certain conditions, while entanglement decays to zero quickly. These remarkable properties of quantum discord makes it more practical than the entanglement to characterize quantum correlations.
As mentioned, the nonlocality is curicial for some applications in quantum information. Therefore, many theoretical studies have been devoted to the comparision of entanglement and the violations of Bell inequalities [29, 30, 31, 32, 33, 34, 35, 36]. It was demonstrated by several authors that the survival time for entanglement is much larger than that of the Bell-inequality violation [29, 30, 31, 35]. In Ref. , the nonlocal entanglement identified by the violations of a Bell inequality is found to be preserved during the time evolution of a system consisting of two qubits which are embedded independently in zero-temperature bosonic reservoirs and are initially in an entangled mixed state. Bellomo et al. showed that Bell inequality might not be violated for a state with high entanglement for a two-qubit system subject to amplitude damping . Moreover, the nonlocal entanglement between two qubits subject to classical dephasing environment modelled as Ornstein-Uhlenbeck noise has been investigated in Refs. [30, 31], very recently, and it was demonstrated that the strong non-Markovian effect can protect the nonlocal entanglement [30, 31]. Additionally, the nonlocal entanglement induced by the qubit-qubit interaction is found to be protected against the sudden death despite the strong influence of the environmental noise for a certain class of initial states regardless of whether the environment is Markovian or non-Markovian .
Although the dynamical relation between the entanglement and quantum discord as well as entanglement and Bell nonlocality for a bipartite system subject to different environments have been investigated by a large number of groups as mentioned above, further investigations would help in clarifying the relation of these three correlation measures. Along these lines, in the present paper, we have investigated the dynamics of quantum correlations such as quantum discord, entanglement and Bell nonlocality between two qubits in a lossy cavity. Very recently, entanglement, Bell nonlocality and quantum discord have been studied in cavity QED [15, 16, 18, 20, 22, 24, 26, 32, 34, 35, 40]. Most of them assume identical atoms, or the atom-cavity resonant case, or only one excitation. Here, we have investigated the effect of atom-cavity detuning and cavity decay rate on the corresponding correlation measures for a system which is initially prepared in separable or maximally entangled state. The most important results of the present study are: depending on the initial preparation of the system and the type of the detunings between the frequency of the cavity mode and the transition frequency of the qubits, the induced quantum discord is found to increase for a long time under cavity decay loss mechanism for certain conditions, while there is no increase in entanglement and Bell nonlocality. In some situations, the quantum discord is induced even while the atoms remain separable at all times. Also, the preservation of induced entanglement and quantum discord was observed under strong dissipation for certain type of initial states.
The paper is organized as follows. In Sec. 2, the model is introduced and the exact differential equations governing the dynamics of two atoms and single cavity mode are derived. The correlation measures; concurrence, CHSH-Bell inequality and quantum discord are also introduced, briefly. In Secs. 3.1 and 3.2, the effect of cavity damping on the corresponding correlation measures are revealed for different detunings for the atoms initially prepared in separable or maximally entangled state. The summary of the important results is given in Sec. 4.
2 The Model and Correlation Measures
Here we consider two two-level atoms with excited states and ground states interacting with a single-mode cavity field of frequency . The Hamiltonian of this system is given as :
where and are the transition frequencies for the atoms and , and are the annihilation and creation operators for the field and and are the coupling constants to the cavity mode for the atoms and , respectively.
When maximum two-photon is considered to be present in the cavity (i.e., the total excitation is equal to two) and if there is no dissipation in the cavity modes, the closed form of the total state vector at any time may be written as [41, 42]:
where are the time-dependent probability amplitudes. For example, represents that the atoms are in their ground states while there are two photons inside the cavity. For this case, the dynamics of the atoms and single mode cavity field can be solved by using time-dependent Schrödinger equation and recently many studies have been reported for this case [41, 42, 43]. On the other hand, when cavity decay is taken into account, since the cavity decay only changes the number of photons inside the cavity, the atoms and the cavity mode can be found in any of the states [40, 41, 44]: . Then the density matrix of the system of two atoms interacting with a single mode cavity field will evolve according to the master equation [40, 41, 45]:
where and are the spin inversion operators. Here the first term in Eq. (3) is generated by the Hamiltonian given in Eq. (1) and represents the coherent evolution, while the other term denotes the incoherent dynamics; is the rate of loss of photons from the cavity mode. By setting , , , , , , , as a basis and considering the Eqs. (1) and (3), the differential equations governing the dynamics of two atoms and the single mode cavity field can be obtained easily. In our analysis we will restrict ourselves for the initial states: and , where . One should note for these initial states the diagonal elements and some of the off-diagonal elements form a closed evolution. They can be listed as:
where and are called the detunings and describes the frequency difference between the energy levels of atoms and the cavity mode, and . We will consider two particular detuning types: for the atoms will be called identical since the transition frequencies of the atoms are the same and the atoms will be called unidentical for case for which the individual atom transition frequencies are different . In practice, the detunings might be slightly different due to broadening, in this work we will consider relatively large values of . One should note that it is quite difficult to obtain an analytic solution for this set of coupled differential equations, therefore we will solve these equations numerically for the considered initial states.
The density matrix of the two atoms is obtained by taking a partial trace over the cavity degrees of freedom: and used to calculate the quantum correlations as measured by CHSH-Bell inequality, concurrence and quantum discord between two atoms interacting with a single mode cavity field. For the considered initial states whose dynamics is given by Eq. (2), the reduced density matrix of the atoms in the two-qubit standart basis can be calculated as:
where and . It should be noted that the density matrix in Eq. (5) has an X-form and the dynamics preserve its form. The considered correlation measures for this type of states can be calculated easily [3, 11, 23, 24, 32]. The concurrence  as an entanglement measure is given by :
where and .
QD is defined as :
where measures the total correlation between the atoms; is the von-Neumann entropy and is the reduced density matrix obtained by tracing over the subsystem . The other quantity is the amount of classical correlations between the atoms defined as the maximum information one can get about the atom by performing measurements on the atom , or vice versa. Its calculation requires maximization over one of the subsystem and generally obtaining an analytic result is not an easy task . Recently, for the density matrix having X-form the explicit expressions of QD and classical correlation are reported [3, 23, 24]. We have used the expressions given in Ref. . The classical correlation for the density matrix Eq. (5) is given as:
where . And QD is given as:
where , with being the eigenvalues of , , , and . We have also calculated QD numerically by considering projective measurements [2, 14] and observed that the results obtained by using Eq. (10) and the numerical results agree perfectly.
As is well known, concurrence is equal to 1 (0) for a maximally entangled (separable) state. Although QD is equal to 1 for a maximally entangled state, it may or may not be equal to zero for a separable state because it was shown that even some separable states can carry non-zero quantum discord [3, 4]. For pure states QD is found to be equal to the entanglement of formation . On the other hand when the CHSH-Bell inequality is violated which signifies that the correlations cannot be accessible by any classical local model .
3.1 Effects of Cavity Damping-Separable States
In this section, we will examine the effect of cavity decay rate, , on the quantum correlations as measured by concurrence, CHSH-Bell inequality and quantum discord. To do this we will solve the differential equations in Eq. (2) numerically for the initial states: and , where and we will set (symmetric coupling), for identical atoms and for unidentical atoms. One should note that the detunings are chosen somewhat large because at higher detunings the qubit-cavity energy exchange is low [41, 42, 48] and a high-degree of cavity-induced quantum correlations can be created .
We first consider the dynamics of concurrence and quantum discord for initially separable state and display and versus dimensionless time, in Figs. 1 and 2 for unidentical and identical atoms, respectively. For this initial state neither non-zero concurrence nor violation in Bell inequality ever occurs for unidentical qubits, so their dynamics are not plotted in Fig. 1. However, this does not indicate the loss of quantum correlations; QD induced by the atom-field interaction oscillates nearly between 0 and 0.18 for (Fig. 1(a)). A similar result is found by Sun et al., in Ref.  where the authors studied the entanglement and quantum discord between two qubits which are independently coupled to a common environment modelled as an Ising spin chain. It can be seen from the analysis of Figs. 1(b)-1(e) that QD is created independently of the cavity decay rate, but its maximum value depends on the decay rate inversely and the effect of cavity decay on the dynamics of QD is to create a damped oscillatory behaviour; damping in QD is proportional to the cavity decay rate. Also, in order to understand the nature of non-zero quantum discord, we have plotted the total correlation, , and classical correlation, , for in the inset of Fig. 1(c). It should be noted that for this initial state, the total as well as classical correlations are also induced and is always greater than as a consequence non-zero quantum discord is present. This inset also demonstrates that and are also damped oscillatorily as . Correlations in the system are gradually lost to the reservoir because of the cavity decay as found in the system studied in Ref. .
We have calculated the dynamics of CHSH inequality, concurrence and quantum discord for the same initial state for the identical atoms () case to understand the effect of the type of detuning for the same leaky cavity and plot the resulting and as a function of in Figs. 2(a)-2(e) and insets of Figs. 2(a)-2(c), respectively. As can be seen from the insets of Figs. 2(a)-2(c), non-zero entanglement, although small, is created which decays to zero as cavity decay rate is increased and the induced entanglement has no dynamics for higher values of ( and ). Bell inequality is found to be not violated for any value of , which is expected because entanglement is low and generally Bell inequality violation is shown to be a hallmark of high [31, 34, 35], so is not displayed in Fig. 2. One should note the difference in the dynamics of QD in the leaky cavity for unidentical and identical atom cases (Figs. 1 and 2). While QD has an oscillatory decay type dynamics for the unidentical atoms (Fig. 1), it increases until the time and then decays for the identical atoms. Also, for the same decay parameters, QD is found to have a much longer lifetime for the identical atoms. The higher leakage in the cavity decreases and the lifetime of QD after the time . On the other hand, has no appreciable effect on the maximum value of . Also note that the fast oscillations in QD are gradually smoothed for larger values of and longer times as can be seen from Figs. 2(b)-2(e).
Now we consider another initial separable state of the form which has one excitation in the field and the other excitation in one of the atoms to see the effect of initial partition of the excitation on the dynamics of the quantum correlations. The time dependent of CHSH-inequality, concurrence and quantum discord are displayed in the first, second and third columns of Fig. 3, respectively. In these figures, insets show the dynamics for unidentical atoms () while main subfigures are for the identical atoms (). Inspection of subfigures demonstrates the pronounced effect of qubits being identical or unidentical. For identical atoms, the atom-field interaction induces high values of and and Bell inequality is also violated at the high values of entanglement for a high-quality cavity (Fig. 3(a), 3(b) and 3(c)). In the case of unidentical atoms, and are also created, but their maximum are small; 0.018 and 0.13, respectively, and is not violated. When unidentical atoms interact with a low-quality cavity, entanglement can cease to exist very quickly and QD shows a fast oscillatory decaying behaviour (see the insets in Figs. 3(e)-3(r)) which is very similar to what we have found for initial state. On the other hand, this situation is completely and significantly different in the case of identical atoms. Although CHSH-Bell inequality goes below 2 very quickly for and no longer is violated for , the dynamics of entanglement and quantum discord are very similar that both approach to constant values in an oscillatory manner; the oscillations die down as the cavity decay rate is increased. In effect, cavity decay entangles initially separable state and the resulting entanglement as well as quantum discord are robust. In many studies, the dissipation of cavity fields were found to play a constructive role in the generation of quantum correlations as measured by concurrence and quantum discord [22, 25, 27, 28, 48, 49, 50]. This situation can be understood by looking at the steady-state of the dynamics given in Eq. (3). For , and , the equation has a solution independent of and and the atom-atom reduced density matrix is given as: . The concurrence for this state is 0.50 and quantum discord of it can be calculated from Eq. (10) to be 0.412 which are seen clearly from Figs. 3(p) and 3(r). Also, we want to emphasize that the dynamics of , and between identical atoms for the initial state is the same as that for , since we assume symmetric coupling (i.e., ) between the atoms and the single mode cavity field.
To further elucidate the asymptotic constant entanglement and quantum discord for the initial state , we display the time dependence of the density matrix elements of the atom-field system given by Eq. (2) and the purity of the atomic subsystem for the density matrix Eq. (5) for identical and unidentical detuning cases in Figs. 4(a)-4(c) for . As can be seen from Fig. 4(a), all populations, except and , become zero for a short time while goes to and go to asymptotically. The key mechanism for the creation of the entanglement and QD in this case seems to be the coherence generated by the interaction of the atom with the dissipative cavity. The absolute value of also reaches to 1/4 asymptotically as can be seen from the Fig. 4(a). The evolved state is a trapping state which has the property that the probabilities of the different decay channels interfere destructively and the state is robust against the cavity loss mechanism [51, 52]. Generation or protecting entanglement by population trapping has been discussed by a number of groups before [48, 49, 51, 52, 53]. The magnitude of has no effect on the long time limit of and , it only affects how long other population components can be nonzero. The dynamics of atomic purity (which is equal to 1 for a pure state and 1/4 for a maximally mixed state in four dimensional Hilbert space) displayed in Fig. 4(c) as the solid line for indicates that the state approaches to a partially mixed state which has high entanglement and quantum discord as indicated before in Figs. 3(p) and 3(r). For unidentical atoms, the short time dynamics is similar to the identical detuning case that all states are populated for a short time, but the asymptotic state is now as can be seen from Fig. 4(b). The atomic purity (dotted line in Fig. 4(c)) indicates that the state of the atoms firstly goes to a mixed state in a short time and then approaches to the pure state which is the reason of the loss of correlations between unidentical atoms in the asymptotic limit.
3.2 Effects of Cavity Decay-Maximally Entangled Initial State
In this subsection, we will investigate the effect of initial correlations in the atom-atom part of the state on the dynamics of quantum correlations.
We take as the initial state which is the maximally entangled Bell state for the atomic system. QD and Bell inequality violation are also maximal for this initial state. We display , and as a function of for several cavity decay rates for unidentical (Fig. 5) and identical (Fig. 6) atoms. At , all correlations show oscillatory behaviour; for the unidentical atoms , and oscillate nearly in the ranges [, 2.2], [1.0, 0.63] and [1.0, 0.58], respectively (Fig. 5). On the other hand, and for have periodic sudden deaths and births for the identical detunings of the atoms (Figs. 6(a) and 6(b)), while oscillates nearly between 0.3 and 1 (Fig. 6(c)). For small values of , the effect of cavity decay on the dynamics of correlations can be described as underdamped oscillations for both identical and unidentical detuning cases. As is increased to , the oscillation becomes overdamped and finally at large ’s () damping becomes critical, i.e. and show exponential decay for the unidentical atom cases. For the identical atoms at and concurrence has sudden death, while QD decays exponentially. One should note that Bell nonlocality is the most fragile quantum correlation considered here; it suffers sudden death in a short time as long as .
We have investigated the dynamics of quantum correlations as quantified by entanglement, Bell nonlocality and quantum discord for two qubits embedded in a common leaky cavity for three different initial states with double excitations. The qubits were considered as identical or unidentical based on the detuning of their transition frequency with respect to the frequency of the field mode. The main findings of the study can be summarized as: Qubits being identical or not have profound influence on the dynamics of the considered correlations; starting from the same initial state, the correlations show very different behaviour depending on whether the detunings of the individual qubits satisfy or . For unidentical qubits initially in their excited states and interacting with a high or low quality cavity, quantum discord is found to be induced for all considered cavity decay rates and it decays oscillatorily for a leaky cavity although the qubits remain unentangled for all times. For the same initial state and the leaky cavity, QD is found to increase in the long time for identical qubits, while entanglement is either ”never created” or decays to zero in a short time.
One of the most interesting finding is the constant high entanglement and quantum discord state which is obtained as the long-time limit of the dynamics of initial state for the identical qubit case. Here the final state is independent of the cavity decay rate (if ) and the magnitude of the detuning and has high concurrence () and quantum discord (). From an analysis of the time dependence of the density matrix elements it is seen that the considered initial state leads to a trapping state which happen to have high degree of quantum correlations. On the other hand, for the same initial state, the correlations between the unidentical atoms are found to be destroyed quickly even for a small cavity decay rate.
For the atomic Bell-type initial state we have found ESD for the identical qubits interacting with a low quality cavity, while concurrence decays exponentially for the unidentical qubits. Quantum discord is found to have asymptotical decay behaviour independent of the type of atoms for this initial state and leaky cavity. On the other hand, CHSH-Bell inequality violations survive only for a finite time for for both identical and unidentical qubits. It is interesting to note that steady-state nonzero correlation state of (or ) is not observed for the . This shows that the stationary states depend strongly on the initial conditions .
We have also demonstrated that Bell nonlocality as manifested by CHSH inequality is the most fragile correlation measure under dissipative environment and is violated only at the high values of entanglement.
-  Nielsen M A and Chuang I L, Quantum computation and quantum information (Cambridge University Press, 2000).
-  Ollivier H and Zurek W H 2001 Phys. Rev. Lett. 88 017901
-  Ali M, Rau A R P and Alber G 2010 Phys. Rev. A 81 042105
-  Henderson L and Vedral V 2001 J. Phys. A: Math. Gen. 34 6899
-  Datta A, Shaji A and Caves C M 2008 Phys. Rev. Lett. 100 050502
-  Knill E and Laflamme R 1998 Phys.Rev. Lett. 81 5672
-  Bell J S 1964 Physics 1 195
-  Acin A, Gisin N, Masanes L and Scarani V 2004 Int. J. Quantum Inf. 2 23
-  Acin A, Gisin N and Masanes L 2006 Phys. Rev. Lett. 97 120405
-  Gisin N and Thew R 2007 Nature Photon 1 165
-  Altintas F and Eryigit R 2010 Phys. Lett. A 374 4283
-  Sun Z, Lu X M and Song L 2010 J. Phys. B: At. Mol. Opt. Phys. 43 215504
-  Lu X M, Xi Z, Sun Z and Wang X 2010 Quant. Inf. Comput. 10 0994
-  Werlang T, Souza S, Fanchini F F and Boas C J V 2009 Phys. Rev. A 80 024103
-  Fanchini F F, Werlang T, Brasil C A, Arruda L G E and Caldeira A O 2010 Phys. Rev. A 81 052107
-  Wang B, Xu Z Y, Chen Z Q and Feng M 2010 Phys. Rev. A 81 014101
-  Werlang T, Trippe C, Ribeiro G A P and Rigolin G 2010 Phys. Rev. Lett. 105, 095702
-  Ge R C, Gong M, Li C F, Xu J S and Guo G C 2010 Phys. Rev. A 81 064103
-  Werlang T and Rigolin G 2010 Phys. Rev. A 81 044101
-  Pei P, Li C, Jin J S and Song H S 2011 J. Phys. B: At. Mol. Opt. Phys. 44 035501
-  Mazzola L, Piilo J, and Maniscalco S 2010 Phys. Rev. Lett. 104 200401
-  Francica F, Plastina F and Maniscalco S 2010 Phys. Rev. A 82 052118
-  Hao X, Ma C L and Sha J 2010 J. Phys. A: Math. Theor. 43 425302
-  Wang C Z, Liu C X, Nie L Y and Li J F 2011 J. Phys. B: At. Mol. Opt. Phys. 44 015503
-  Fanchini F F, Castelano L K and Caldeira A O 2010 New. J. Phys. 12 073009
-  Li G X, Yi Z and Ficek Z 2011 arXiv:1101.4983v1
-  Zhang J S, Chen L, Aty A and Chen A X 2011 arXiv:1101.5429v1
-  Yuan J B, Kuang L M and Liao J Q 2010 J. Phys. B: At. Mol. Opt. Phys. 43 165503
-  Kofman A G and Korotkov A N 2008 Phys. Rev. A 77 052329
-  Li J Q and Liang J Q 2010 Phys. Lett. A 374 1975
-  Altintas F and Eryigit R 2010 J. Phys. A: Math. Theor. 43 415306
-  Mazzola L, Bellomo B, Franco R L and Compagno G 2010 Phys. Rev. A 81 052116
-  Derkacz L and Jakobczyk L 2005 Phys. Rev. A 72 042321
-  Bellomo B, Franco R L and Compagno G 2009 Adv. Sci. Lett. 2 459
-  Bellomo B, Franco R L and Compagno G 2008 Phys. Rev. A 78 062309
-  Werner R F 1989 Phys. Rev. A 40 4277
-  Yu T and Eberly J H 2004 Phys. Rev. Lett. 93 140404
-  Ferraro A, Aolita L, Cavalcanti D, Cucchietti F M and Acin A 2010 Phys. Rev. A 81 052318
-  Xu J S, Xu X Y, Li C F, Zhang C J, Zou X B and Guo G C 2010 Nat. Commun. 1 7
-  Mazzola L, Maniscalco S, Suominen K A and Garraway B M 2009 Quantum Inf. Process. 8 577
-  Pathak P K and Agarwal G S 2004 Phys. Rev. A 70 043807
-  Kim M S and Agarwal G S 1998 Phys. Rev. A 57 3059
-  Z. Deng 1985 Opt. Commun. 54 222
-  Mahmood S and Zubairy M S 1990 J. Mod. Optic 37 243
-  Agarwal G S and Gupta S D 1990 Phys. Rev. A 42 1737
-  Wootters W K 1998 Phys. Rev. Lett. 80 2245
-  Clauser J F, Horne M A, Shimony A and Holt R A 1969 Phys. Rev. Lett. 23 880
-  Francica F, Maniscalco S, Piilo J, Plastina F, Suominen K A 2009 Phys. Rev. A 79 032310
-  Maniscalco S, Francica F, Zaffino R L, Gullo N L and Plastina F 2008 Phys. Rev. Lett. 100 090503
-  He Z, Zou J, Shao B and Kong S Y 2010 J. Phys. B: At. Mol. Opt. Phys. 43 115503
-  Garraway B M 1997 Phys. Rev. A 55 4636
-  Bellomo B, Franco R L, Maniscalco S, Compagno G 2008 Phys. Rev. A 78 060302
-  Wang F Q, Zhang Z M, Liang R S 2008 Phys. Rev. A 78 042320