Cavity assisted generation of sustainable macroscopic entanglement of ultracold gases

Prospects for reaching persistent entanglement between two spatially separated atomic Bose-Einstein condensates are outlined. The system set-up comprises of two condensates loaded in an optical lattice, which, in return, is confined within a high-Q optical resonator. The system is driven by an external laser that illuminates the atoms such that photons can scatter into the cavity. In the superradiant phase a cavity field is established and we show that the emerging cavity mediated interactions between the two condensates is capable of entangling them despite photon losses. This macroscopic atomic entanglement is sustained throughout the time-evolution apart from occasions of sudden deaths/births. Using an auxiliary photon mode and coupling it to a collective quadrature of the two condensates we demonstrate that the auxiliary mode's squeezing is proportional to the atomic entanglement and as such it can serve as a probe field of the macroscopic entanglement.


Introduction
After years of experimental activity in cavity QED mainly exploring fundamentals of light-matter interactions at a true quantum level [1], in the second half of the 90's, these types of experiments took a new direction; it was demonstrated that they may constitute a platform for implementing quantum information processing [2]. In [3], for example, the light-matter interaction was monitored such that the photon state became entangled with a single two-level atom and then by letting a second atom 'absorb' the light field a two-atom entangled state was obtained. In this protocol, after the first step quantum information is stored in the light field before it gets transmitted to the second atom. Such a scheme is inevitably sensitive to photon losses in the middle step, which, therefor, requires short operation times. It was, however, realized that the 'cavity mediated atom-atom interaction' could also be utilized for directly entangling two atomic qubits without ever populating the lossy photon mode [4]. In such a scenario, photons are only virtually excited, but nevertheless, integrating out the photon field results in an effective atom-atom interaction. Then, by carefully tuning the interaction time the desired entangling operation could be implemented. Somewhat surprisingly, even for a 'classical' photon field, i.e. thermal state, it is indeed possible to use it for entanglement generation between the two atoms [5], and thus one may expect that the entanglement preparation is prone to both photon losses and other possible noise sources originating from the quantized radiation field.
Following the pioneering cavity QED experiments, it has in recent years become possible to coherently couple not only a single atom but a gas of thousands of ultracold, condensed, bosonic atoms to a single quantized light mode [6]. This scales up the effective light-matter interaction strength, and the strong coupling regime has indeed been demonstrated in several groups 1 [7]. At these sub-Kelvin atomic temperatures, the atomic motion must be treated quantum mechanically -it couples to the light field via the photon recoils. This led to explorations of a new, ultracold, regime of cavity QED. Much of the focus has especially been devoted to the collective phenomenon of superradiance 2 [8], realized by coupling the light degrees-of-freedom to the atomic vibrational modes rather than to internal atomic electronic states [6]. Here, an atomic condensate is trapped inside the resonator and an external 'perpendicular' drive/pump scatters photons into the cavity by illuminating the atoms. It turns out that at a critical pump amplitude, the atoms self-organize themselves in such a way that the different scattering processes constructively interfere making it possible for the photons to enter the resonator [9]. Below the critical point, the condensate is approximately uniformly distributed and scattering processes add up destructively leading instead to a vanishing cavity field. This change in the system properties has shown to share many resemblances with the Dicke phase transition from a normal (vacuum cavity field) to a symmetry broken superradiant phase (coherent cavity light field) [10]. Properties of this transition has been the subject of several experimental activities 3 , for example the symmetry breaking has carefully been mapped out and also the critical exponents [11,14].
Like for the single/few atom cavity QED realizations, a natural question arises whether also these 'many-body' cavity QED systems can serve as toolboxes for quantum information processing. A first step along these lines would be to study possibilities of entanglement generation among the many atoms. A crucial difference compared to the traditional cavity QED schemes is that these systems are 1 For N two-level atoms equally coupled to the field, the collective coupling goes as √ N , and thus the characteristic time scale (inverse of the Rabi frequency) is greatly reduced. 2 In this setting, superradiance does not necessarily manifest itself as an enhanced radiation process [8] but instead as collective properties of the atomic condensate.
driven. Thus, one should rather consider possible persistent entanglement in the long time limits of the driven-lossy system [15]. Just like for Refs. [3][4][5], the photon mode plays the role of an ancilla system mediating atom-atom entanglement. The situation is different from those of [4,5] though since the non-vanishing field amplitude is sustained from balancing the gain from the input field to the losses due to photon decay, while in Refs. [4,5] the system is closed and the photon field initialized in some state. Furthermore, in the present setup the atomic fields will inevitably get entangled with the photon field, which causes an effective decoherence for the atomic state, i.e. the reduced state for the atoms will be mixed. In addition, photon dissipation gives rise to another decoherence channel for the atomic fields. All in all, even if in the few atom scenario quantum coherences of the radiation field is not necessary for atomic entanglement generation [5], due to the differences with the present many atomic situation it is a priori not clear whether any atom-atom entanglement may survive.
Inspired by a new system setup of the Esslinger group [16] we analyze the above problem by considering two spatially separated atomic Bose-Einstein condensates (BECs) located within a cavity and perpendicularly driven by an external laser. Entanglement properties between the two atomic condensates (as previously mentioned the atomic degrees-of-freedom are collective atomic vibrational modes) are studied. As each condensate constitutes a macroscopic mechanical oscillator, the BEC-BEC entanglement is manifested in the corresponding vibrational states. By calculating the logarithmic negativity, persistent entanglement is demonstrated even in the open system case when photon losses are taken into account. More precisely, the BEC-BEC entanglement displays an infinite series of entanglement sudden deaths/births such that on average the two condensates keep their entanglement for an infinite time. Entanglement between the two condensates is, however, only possible in the superradiant regime. Close to and at the critical point, the main part of the quantum correlations are shared between the photon field and the atoms, while deeper in the superradiant phase the BEC-BEC entanglement becomes more prominent. This is an example of entanglement sharing in a tripartite system [17]. In order to probe the entanglement between the two condensates we introduce an auxiliary photon mode and couple it to a collective quadrature of the two condensates. We demonstrate that near the critical point the squeezing of this auxiliary mode scales linearly with the the BEC-BEC entanglement and thus provides an indirect readout of the entanglement. However, losses of this additional photon mode implies that in the long time limit no quantum entanglement survives between the two condensates.

Model system
By now, achieving a coherent light-matter coupling between a single BEC and a quantized photon mode has been demonstrated in several groups [6,7,11,18,19]. The coupling between single photons to the atomic collective vibrational modes realizes a longstanding goal of optomechanics, i.e. coherent interplay between photons and a macroscopic mechanical oscillator prepared in its ground state [20]. This has spur great activity in exploring hybrid quantum systems with the objective to devise scalable quantum architectures [21]. In particular, generating nonclassical states in atom-photon coupled hybrid quantum systems has received significant theoretical and experimental interest [22,23]. Continuing this quest, we envision two spatially separated BECs confined inside an optical resonator and explore whether macroscopic entanglement between two atomic BECs can be generated via the coupling to a common photon mode.
With the particular system in mind, the model Hamiltonian is given by [9] Here, the parameter ∆ is the cavity-pump detuning, L is the cavity length, m the atomic mass, k (= (0, 0, k)) the wave number, U 0 = g 2 0 /∆ a is the effective atom-cavity coupling with g 0 the bare Rabi frequency and ∆ a the atom-pump detuning, and η = Ωg 0 /∆ a with Ω the pump amplitude. The photon creation (annihilation) operator isâ † (â), andΨ j (x) is the atomic field operator annihilating an atom in condensate j at position x. The operators obey â,â † = 1, Ψ i (x ),Ψ † j (x) = δ ij , and the remaining commutators are identically zero. In deriving (1) the excited atomic electronic state has been adiabatically eliminated under the dispersive assumption |∆|, |∆ a | g 0 , and furthermore we consider a vanishing overlap between the two condensates such that the cross terms of the atomic operators can be neglected.
As a next step, we second quantize the Hamiltonian by expanding the field operators in the two lowest vibrational modes [9,11] whereĉ (j) 0,1 annihilates an atom of vibrational mode 0, 1 in condensate j. Inserting the operators (2) into the expression (1) for the Hamiltonian gives the low energy second quantized model, now containing five independent boson modes. Taking the conservation of number of atoms of the two condensates into account, it is practical to employ the Schwinger spin-boson mapping [24] to recast the Hamiltonian in a generalized Dicke model where The standard Dicke model is lacking the last term and the boson mode couples only to a single spin S 4 . We have introduced the total number of atoms N ; the effective atom-light coupling g = √ 2N η and the Stark shift term U = N U 0 /4. With this parameter scaling, for large atom numbers every term in the Hamiltonian O(N ). This ensures that in the thermodynamic limit, N, L → ∞ and N/L =constant, the model exhibits a second order phase transition [10]. In the following we will omit the Stark shift term, i.e. setting U = 0, valid in the regime were the pump amplitude is much larger than the Rabi frequency which is justified in typical experiments [11]. In absence of any losses the couplingg c = |∆|ω R /2 separates the two phases, normal (g <g c ) from the superradiant (g >g c ) one. In the normal phase, the ground state is simply the vacuum state of all modes (Ŝ z =T z = −N/4). In the superradiant phase, however, a Z 2 parity symmetry is spontaneously broken which is for instance reflected in a 0 or π phase of the photon field [14]. The spins in this phase are tilted away from the south pole toward the 'x'-direction [25]. Classically, the phase transition manifest as a pitchfork bifurcation [26], and it has been argued that entanglement generation should be most pronounced in the symmetry broken phase [27]. The cavity induced atom-atom interaction which underpin the entanglement is transparent when the cavity field is eliminated giving rise to an effective atomic model of the Lipkin-Meshkov-Glick type [28]. Entanglement properties of such models have been studied in the Refs. [22,28]. Integrating out the photon field and imposing an adiabatic approximation in order to derive a Lipkin-Meshkov-Glick model is not always justified, for example for moderate detunings ∆ and good cavities (low photon decay rates κ). Clearly, in this case the full system has to be considered. This may be extra important when studying entanglement properties since the presence of the cavity field will itself get entangled with the atomic fields and thereby indirectly degrade atom-atom entanglement. In addition, the photons may leak out of the cavity which result in an additional decoherence mechanism for the atoms.
To take photon decay into account we follow the standard procedure by coupling the cavity mode to a set of harmonic bath oscillators and apply the Born-Markov and secular approximations that are supposed to be valid in the optical regime [29]. Doing so we end up with a Markovian master equation of the Lindblad form. Assuming a zero temperature bath the resulting open dynamics can be modeled through the following master equation [29] andρ is the full system density operator. We should keep in mind that we consider a driven system out of equilibrium. This justifies a phenomenological introduction of losses as modeled through the master equation (5). In particular, the light-matter coupled Hamiltonian (4) is written in the frame of an external drive imparting it a non-equilibrium character. For a time-independent system coupled to an external bath it is required that the steady state will obey the principles of equilibrium statistical mechanics, while here no such restriction is enforced on the dynamics which arises out of an implicit time-dependent Hamiltonian [30].

Semiclassical analysis
As we are interested in entanglement generation on a macroscopic level we assume the number of atoms N 1. This is typically the case in experimental realizations were N 10 5 [11,14,18,31]. For particle numbers of this size, mean-field approaches qualitatively captures many of the system characteristics [11]. However, such a semiclassical mean-field approach is incapable of capturing or quantifying the quantum correlations. For bosonic (spin) systems, such quantum corrections are conveniently explored in terms of the Holstein-Primakoff representation [32], from which we can extract information about the entanglement. Before analyzing the quantum features we first need to determine the mean-field solutions of the model which are the steady state solutions of the semiclassical equations of motion governed by the master equation (5).
Let us define â = α, Ŝ − = β, T − = δ, Ŝ z = wŜ, and T z = wT . The semiclassical equations of motion are given from the Heisenberg equations under the factorization assumption, i.e., This results in a finite set of coupled nonlinear equationṡ At this level we do not include any back-action of the reservoir apart from the dissipative photon decay. This is justified when the cavity mode is coupled to a heat-bath in thermal equilibrium and the associated Langevin noise operators thus satisfy â in (t) = 0 [33]. Spin conservation for each atomic BEC at the mean-field level takes the form Making use of the above equations we can find the steady state solutions for α, β, δ, wT , wŜ. We find that the solution predicts the occurrence of two distinct phases (normal and superradiant) which are separated by a critical value of the light-matter coupling g c = .Note how the critical coupling is modified by the photon decay rate κ [9,12,34]. The open nature of the system not only shifts the transition, but also alters the critical exponents [35]. For g < g c , the steady state solution is a fully inverted (trivial) state with α = β = δ = 0, wT = wŜ = −N/2. This represents the normal phase. For g > g c , however, this trivial state becomes unstable and two new stable steady states arise with This is the aforementioned pitchfork bifurcation [26]. In particular, we will show later that a consequence of the light-matter coupling exceeding a critical value (g/g c > 1) is the generation of bi-partite entanglement between the two distant BEC samples, while they remain separable in the normal phase.

Linearized light-matter interaction
As argued above, in the large particle regime the quantum fluctuations are small and may be treated in a linearized approach. Starting with the normal phase we use the Holstein-Primakoff representation [32] whereb andĉ are the two new boson modes describing the bosonic excitations around the mean-fields (13), to reformulate the Hamiltonian (4) in terms of three boson fields. By further making use of the large N limit, we linearize the resulting Hamiltonian to arrive at In the superradiant phase (g/g c > 1) the cavity and the atomic field modes are macroscopically excited and the expansion of the Hamiltonian (4) is preceded by first displacing the operatorsâ,b,ĉ around their steady state amplitudes (13). This gives the following linearized Hamiltonian where Note that in the limit g → g c ,Ĥ lin sr equalsĤ lin n . In these linearized forms, both the Hamiltonians are quadratic for which entanglement properties have been thoroughly studied [36]. For example, the logarithmic negativity N [37] used in the next section is a good measure in order to quantify the bi-partite entanglement. To further simplify the Hamiltonians (18) and (19) we define the following set of collective operatorsp In terms of these the Hamiltonian in the normal phase (18) takes the form Similarly, the Hamiltonian in the superradiant phase (19) can be expressed as where From the above expressions it is seen that the relative modeŝ decouples from the remaining two modes. So far, the open character of the system only enters in the system parameters of the above effective low energy models. We have not taken the Lindblad terms into account. Since losses only appear for the photon mode, the Lindblad term of Eq. (5) is unaltered going to the Holstein-Primakoff representation. Nevertheless, working in the collective basis (20) we should also transform the loss term; whereĤ =Ĥ lin n orĤ =Ĥ lin sr depending on the phase studied. Note that since theŝ mode is purely atomic it does not couple to the other modes even when the system interacts with the surrounding reservoir. In particular, in the normal phase of the closed system the ground state for theŝ mode is the vacuum. In the superradiant phase, however, the corresponding ground state is a squeezed vacuum [8]. Thus, if the system'sŝ mode is initialized in vacuum it is not an eigenstate and in particular it evolves accordingly A consequence of this time-dependence of theŝ mode is that correlators involving the modesb andĉ may also become time-dependent, even though the collective modesp andq may approach a stable fixed point. For instance, if theâ,b, andĉ modes are initialized in their respective vacuum states then in the steady state one gets b †b ss = ĉ †ĉ ss = 1 4 p †p +q †q +p †q +q †p ss − 1 2 Similar expressions can be obtained for other correlators containing theb andĉ operators.

Sustainable entanglement generation
It should be remarked that the master equations (23) preserve the initial Gaussian character of the modesp,q, andŝ (or, equivalently of modesâ,b, andĉ). Using standard quantum optical techniques it allows us to convert (23) into partial differential equations for the quantum characteristic function defined as χ( p , q ) = e pp † e − * pp e qq † e − * qq . We will not enter into the details of this exercise, but will refer the reader to the same approach described in [30]. As pointed out earlier, the solutions verify that theb andĉ modes remain separable throughout the normal phase. In the superradiant phase, however, the atomicb andĉ modes can become entangled 5 . We characterize the bi-partite entanglement between the various modes in terms of the logarithmic negativity N [37].
For times t = kπ/ Ω 2 − 4ζ 2 , k ∈ Z theb andĉ modes do not directly couple to theŝ mode. Bi-partite entanglement between theb andĉ modes for such a choice of t is shown in Fig. 1 (b) and (d). Also shown in Fig. 1 (a) and (c) is the bi-partite entanglement between the cavity mode and one of the atomic BECs. As can be seen from Fig. 1 (a) and (c), bi-partite entanglement develops between the cavity field and each atomic cloud. This is true in the normal and superradiant phase. However, it is only when g/g c > 1 that the entanglement between distant atomic BECs develops. Entanglement between the cavity field and each atomic cloud is maximal near the critical point g/g c = 1, which may serve as an indicator of the quantum critical point. This is indeed a universal behavior; at the critical point where the length scales (in system with short range interactions) diverge the entanglement also diverges [38]. Furthermore, as a function of the effective detuningδ = −∆ − ω R the maximum entanglement is established at the critical point and at resonanceδ = 0. It is interesting to observe that in our model, which does not belong to the class of short range interaction ones, the bi-partite entanglement reaches its maximum at the critical points only in the cavity-BEC subsystems, but not in the BEC-BEC subsystems. In fact, as bi-partite BEC-BEC entanglement builds up (increasing the coupling strength) the cavity-BEC entanglement decreases. This is the phenomenon of entanglement sharing [17], which for a tripartite system implies that simultaneous large bi-partite entanglement among all constituents is 5 We also find that theâ andb modes (or, equivalently, modesâ andĉ) do exhibit bi-partite entanglement both in the normal and superradiant phases. and (d). The fact that the cavity-BEC or the BEC-BEC entanglement roughly complement one another in terms of amplitude is an outcome of entanglement sharing [17]. We also note that atom-atom entanglement occurs only in the superradiant phase and beyond a critical value of the detuning, i.e. |∆| > ∆ c , i.e. it is not symmetric with respect to the photon-atom detuningδ = −∆ − ω R . In (a) and (b) ∆ = −2, while in (c) and (d) κ = 0.05 and for all four plots t = kπ/ Ω 2 − 4ζ 2 , k ∈ Z. The value of the loss rate κ has been taken small, but still experimentally relevant [11,14,31], in order to achieve as large entanglement as possible. The (scaled) recoil frequency ω R = 1.
not allowed. We note that for a special case of our model this effect has been demonstrated namely for the closed Tavis-Cummings model (i.e. the Dicke model with the rotating wave approximation) for two two-level atoms [39]. An additional feature evident from Fig. 1 is that BEC-BEC entanglement appears only beyond a critical value of the detuning, i.e. ∆ < ∆ c . In other words, being in the superradiant regime is not always sufficient to warrant BEC-BEC entanglement.
The results reported in Fig. 1 are obtained for a particular choice of parameter t. Experimentally it is more relevant to ask the question how the time-averaged entanglement depends on the system parameters. One may, for example, imagine a continuous measurement and then averaging the results, or making measurements at random times and then averaging. Thus, we evaluate the time-averaged logarithmic negativity as a measure of bi-partite entanglement among photon and atomic modes, i.e. we evaluate N = T 0 N(t)dt/T for large sampling times T 's. The resulting time-averaged logarithmic negativity is plotted in Fig. 2. As is evident from Fig. 2 (a), the time-averaged bi-partite entanglement between the two BEC samples is indeed very fragile to cavity damping. As illustrated in Fig. 2 (b), the time-averaged bi-partite entanglement between the two BECs is a nontrivial function of the cavity field detuning ∆ and the coupling strength g/g c . In particular, the two BECs may become entangled also for red detunings ∆ > −ω R . Naively, one could expect thatδ = −∆ − ω R = 0 should generate maximum BEC-BEC entanglement since then the effective atom-atom coupling is the largest. However, as clear from Fig. 2 (b) the picture is more complex; a strong atom-field coupling implies also strong back-action of the field onto the atoms (and also an increased influence of the photon reservoir on the atom-atom coherence), which can counteract build-up of BEC-BEC entanglement. This tradeoff between increasing the effective atom-atom interaction and decreasing the influence of the cavity field and its decoherence is also reflected in the non-monotonous behavior of the negativity as a function of g/g c . As shown in Fig. 2 (c) due to entanglement sharing constraints in a tri-partite quantum system, the BEC-BEC entanglement grows at the expense of the cavity-BEC entanglement.

Inferring entanglement
A challenging aspect of any scheme involving macroscopic entanglement generation is the actual experimental detection of entanglement. In this section we will propose a strategy to indirectly infer the degree of entanglement generated between the two atomic BECs. The scheme can be used to indirectly infer the degree of entanglement generated between the two BECs in the superradiant phase for small values of the critical coupling ratio g/g c 1.
More precisely, we suggest to make use of an additional (auxiliary) opticalŵ mode to indirectly infer the degree of entanglement between the two atomic samples. The auxiliary mode couples to the relative atomicŝ mode with a strength Ψ as H lin−aux sr = ω pp †p + ω qq †q + Ωŝ †ŝ + (Ω + ∆) 2 (p †q +q †p ) + ζ(pq +q †p † ) + ζ(ŝ 2 +ŝ †2 ) A discussion on an actual physical realization of the above Hamiltonian has been left as an open question which might be studied in a future work. It should be pointed out that the Hamiltonian (26) has a beam-splitter (passive) interaction between the relativeŝ mode and the auxiliaryŵ mode [40]. In this sense, theŵ mode initially prepared in a classical state and evolving under the Hamiltonian (26)   our hybrid tri-partite quantum system, now interacting with an auxiliaryŵ mode, is described by the following master equation Here we have included photon decay of the auxiliary mode with a loss rate γ. It should be pointed out that a damped auxiliaryŵ mode provides an indirect dissipation channel for theŝ mode. In other words, indirect damping introduced by theŵ mode can allow theŝ mode to settle to a time-independent steady state. Fig. 3 shows transient bi-partite entanglement between the atomicb andĉ modes in the presence and absence of the auxiliaryŵ mode. The time decaying envelope of the quantum entanglement between the atomic modes can be attributed to the indirect damping suffered by theŝ mode. More precisely, without losses of theŝ mode the entanglement does not reach a vanishing steady state -the sustainable quantum correlations derive from the losslessŝ mode. With losses of theŵ mode included, however, the BEC-BEC entanglement vanishes in the long time limit. The Fig. 3 also gives a clear signature of entanglement sudden death and birth [41], arising from the common coupling of the two atomic BECs to a single cavity mode. As shown in Fig. 3, entanglement sudden death and birth occur both in the absence and presence of the auxiliaryŵ mode. Similar observation of entanglement sudden death and birth has been reported in a scheme proposing entanglement generation among distant optomechanical systems [42].
In Fig. 4 (a) we plot the logarithmic negativity both with (Ψ = 0) and without (Ψ = 0) the auxiliary mode vs. small values of the critical coupling ratio g/g c 1. The linear dependence in this parameter regime suggests a simple one-to-one mapping between the two cases. Indeed, we find that for values of the critical ratio g/g c 1, the time-averaged quantum entanglement between theb andĉ modes in the absence of the auxiliaryŵ mode (N Ψ=0 ) is proportional to the same quantity in the presence of auxiliaryŵ mode (N Ψ =0 ). Thus, the effect of the auxiliary mode on the atom-atom entanglement is to just scale it. This feature is captured in Fig. 4 (b). Moreover, the single mode squeezing S Ψ =0 [8] of thê w mode is also found to behave linearly with the bi-partite entanglement between theb andĉ modes, i.e. N Ψ =0 ∝ S Ψ =0 as shown in Fig. 4 (c). The connection between entanglement and squeezing has long been known [43], and as the later can be determined by means of homodyne measurements on the auxiliarŷ w mode [8], this allows the possibility of indirectly quantifying the macroscopic entanglement between the two BECs. It should be possible to subject the degree of single mode squeezing to a post-processing procedure and one can then rescale the bi-partite entanglement (N Ψ=0 ) according to the linear behavior shown in Fig. 4 (b) and (c).

Conclusions
In this work we have explored a possibility to generate macroscopic entanglement between two BEC atomic samples. We considered a physical setting where two laser driven spatially separated BECs are placed in a high finesse optical cavity where the photons can leak out of the cavity at a rate κ. The two BEC samples interact indirectly through the same cavity field imposing an indirect atom-atom interaction which generates the entanglement. In the thermodynamic limit our physical model is known to exhibit a Dicke-type phase transition when the light-matter coupling exceeds a critical value. This phase transition separates the normal phase from a superradiant phase. Starting with a semiclassical analysis of the steady state of our driven-dissipative model, and adding a linearized treatment of quantum fluctuations using the Holstein-Primakoff representation we solve the resulting master equations in the two phases separately. We find that the two BEC samples remain separable in the normal phase but can become entangled in the superradiant phase. This generation of bi-partite entanglement between the two BECs can be seen as a signature of our hybrid system entering the superradiant phase where quantum coherence persists over macroscopic length scales. In particular, by entering the superradiant phase the cavity-BEC entanglement peaks at the critical point and almost no BEC-BEC entanglement exists at this point. Deeper in the superradiant phase the phenomenon of entanglement sharing occurs where the Figure 3. Time-evolution of the logarithmic negativity N, serving as a measure of entanglement, between the atomicb andĉ modes in the presence, Ψ = 0, (green) and absence, Ψ = 0 (blue) of the auxiliaryŵ mode. Other parameters include ∆ = −1, ω R = 1, κ = 0.05, γ = 0.05κ, and g/g c = 1.05. In the presence of the auxiliaryŵ mode the logarithmic negativity N monotonically decreases with time. The inset shows transient oscillations of N Ψ=0 and N Ψ =0 featuring multiple occurrences of entanglement sudden death and sudden birth. At these short times, the influence of the decay of theŵ mode is not apparent; we note that the characteristic time scale for this decay γ −1 = 400. Note that with γ = κ = 0.05 we would still observe entanglement for t 50 and thereby also several deaths/births of the entanglement. cavity-BEC entanglement is lowered while the BEC-BEC entanglement gets stronger. Even though the open-driven system does not approach a steady state in the general case, the time-averaged BEC-BEC entanglement is persistent. Since the two BECs are coupled to the same leaky photon mode, we also encounter an entanglement sudden death/birth evolution. We have also outlined a strategy to indirectly infer the degree of entanglement between the two atomic BECs. This scheme is no direct measure of the entanglement between the two condensates, but it should be a clear smoking gun of it. . Entanglement (logarithmic negativity) between the two atomicb andĉ modes as a function of g/g c (a) in the presence (Ψ = 0) and absence (Ψ = 0) of the auxiliaryŵ mode. The linear dependences in the two cases lead to proportionality shown in (b); N Ψ =0 vs. N Ψ=0 . Dots are numerical data, while the solid line a least square linear fit. In (c) the logarithmic negativity is displayed vs the single mode squeezing S Ψ =0 of the auxiliary mode. As is clear, increasing the atom-atom entanglement implies an increased squeezing of theŵ mode, which suggests that squeezing would be an indirect measure of entanglement in the system. Other parameters include ∆ = −1, ω R = 1, κ = 0.05, and γ = 0.05κ. Jonas Larson conceived the main idea of the project, while Chaitanya Joshi carried out most of the theoretical and numerical calculations. Both authors contributed to writing and revising the manuscript.