Steady State Entanglement in Cavity QED

We investigate steady state entanglement in an open quantum system, specifically a single atom in a driven optical cavity with cavity loss and spontaneous emission. The system reaches a steady pure state when driven very weakly. Under these conditions, there is an optimal value for atom-field coupling to maximize entanglement, as larger coupling favors a loss port due to the cavity enhanced spontaneous emission. We address ways to implement measurements of entanglement witnesses and find that normalized cross-correlation functions are indicators of the entanglement in the system. The magnitude of the equal time intensity-field cross correlation between the transmitted field of the cavity and the fluorescence intensity is proportional to the concurrence for weak driving fields.


I. INTRODUCTION
The study of entanglement has emerged as a central theme of quantum physics in recent years. It is driven both by fundamental questions and by the increasing interest in applications that go beyond the limit of classical physics. Entanglement as a measurable quantity is a complicated subject, in particular when the systems have multiple components. Here we choose to study entanglement and its possible avenues of quantification in an open quantum system. This system, the canonical model of cavity QED [1], has a single atom coupled to the mode of an optical cavity with two reservoirs or avenues for extracting information: spontaneous emission and losses from the cavity.
Two particles (or systems), A and B are said to be in an entangled state if the wave function of the complete system does not factorize, that is |AB = |A |B . One consequence of this form of the wavefunction is that a measurement on system A yields information about system B without any direct interaction with system B. For systems with the same dimension, in particular, a (pure) state is said to be maximally entangled if tracing over one of the two systems, say A, leaves the other one in a totally mixed state; this means that one can gain complete knowledge of system B by performing measurements on A only.
An example that is of relevance to this work is the maximally entangled state of an atom and a field mode, |Ψ = (1/ √ 2) (|1, g + |0, e ) with the first index denoting the number of photons in the field mode and the second (e = excited, g = ground) denoting the state of the atom. A measurement of the state of the atom immediately tells us the number of photons in the field mode; or a measurement of the photon number immediately tells us the state of the atom.
The von Neumann entropy E = −tr A (ρ A log 2 ρ A ) of the reduced density matrix of system A, ρ A = tr B (ρ AB ) [2] quantifies the amount of entanglement in a given bipartite quantum system in a pure state. For mixed states, on the other hand, although it is easy enough to define what is meant by a totally unentangled state-namely, one in which it is possible to represent the density operator as an incoherent superposition of factorizable states-quantifying the amount of entanglement in a partially entangled state is not, in general, simple. The natural generalization of the pure-state measure indicated above, known as the entanglement of formation, utilizes a decomposition of the quantum state ρ = j P j |ψ j ψ j | = j P j ρ j , and then defines E = min( j P j E j ) where E j is the von Neumann entropy for the density matrix ρ j = |ψ j ψ j |, and the minimum is taken over all the possible decompositions, which is in general a very challenging task [2,3]. As a result of this, alternative measures have been proposed, such as the logarithmic negativity [4]. It is also possible that some particular measurement scheme may result in a most natural unraveling of the density operator, in the sense of the quantum trajectories approach [5] (especially for systems that are continually monitored), and in that case it may be physically meaningful to focus only on the entanglement of the (conditionally pure) states obtained via that particular unraveling.
One of the main purposes of this paper is to determine how much information about the atom-field entanglement in our canonical cavity QED system can be gleaned from the kinds of measurements represented by the traditional correlation functions of quantum optics. As we shall show below, we are actually able to avoid the difficulties for mixed-state entanglement because, in the limit we are interested in, our system is, to a good approximation, in a pure state, in spite of its being an open system interacting with two reservoirs. Fig. 1 shows a two level atom in a driven optical cavity. We consider a single-ended cavity, with the intracavity field decaying via the output mirror at rate κ. The two-level atom has a spontaneous emission rate to modes out the sides of the cavity denoted by γ, which is generally less than the free space Einstein A coefficient. The resonant coupling between the atom and the field mode is given by g = µ eg ω/2 ǫ 0 V with µ eg the electric dipole matrix element, ω the transition frequency and V the volume of the cavity mode. The driving field is taken to be a large classical field ǫ incident on the input mirror, with small transmission T in , so that the incident flux (in photon units) inside the cavity is proportional

II. CAVITY QED SYSTEM
The quantum trajectory wave function that characterizes the system under a non-Hermitian Hamiltonian is: |ψ c (t) = ∞ n C g,n (t)e −iEg,nt |g, n +C e,n (t)e −iEe,nt |e, n with collapse operators associated with photons exiting the output mirror and spontaneous emission out the side of the cavity. The indices e(g) indicate the atom in the excited (ground) state, while n is the number of photons in the mode. The energies are E e,n = E g,n+1 = ω(n + 1/2). We have the usual creation (a † ) and annihilation (a) operators for the field, and Pauli raising and lowering operators σ ± for the atom.
In the weak driving limit, the system reaches a steady-state wave function: where the A ij are known [6,7]. They are A 0,e = β A 1,e = αβq (8) The quantities p and q would be 1 for coupled harmonic oscillators. In cavity QED they differ from unity due to the non-harmonic, or saturable, nature of the atom. The squares of coefficients of single excitation A 1,g , A 0,e give the rates of detection of single photons through the output mirror or in fluorescence (steady state), while the squares of the double excitation coefficients A 1,e , A 2,g give the rates of detection of two photons either in coincidence (one through the mirror, and one in fluorescence) or both out of the mirror. The variables are The one-excitation amplitudes A 1,g and A 0,e are proportional to the driving field ǫ; the twoexcitation amplitudes A 2,g , and A 1,e are proportional to the square of the driving field, ǫ 2 . The entanglement of formation for this system is calculated from the density matrix after tracing over the field variables: The eigenvalues of this matrix are, to lowest nonvanishing order, where q is defined in Eq. (13), and we have defined The entropy E = −λ 1 log 2 λ 1 − λ 2 log 2 λ 2 is then (again to lowest leading order) where we have taken the weak field limit, ǫ being the smallest rate in the problem, so ǫ/κ ≪ 1. The approximation (22) will hold provided (ǫ/κ) 2 ≪ |ξ|.
This entropy is the same as that obtained by using the density matrix for the field alone, traced over the atomic degrees of freedon.
The concurrence, first introduced by Wooters for two qubits [3], can also be used to characterize entanglement between two quantum systems of arbitrary dimension [8,9,10,11]. The concurrence for our system is To see why |ξ| ∝ |A 1,e − A 0,e A 1,g | may be a good indication of entanglement, consider what happens if the wavefunction is a product state. We could write For weak excitations, the coefficient of the ground state of the system is D 0 C g = 1, or Then the product state is Just knowing the one excitation amplitudes does not yield any information about entanglement, as it is possible to have A 1,g = D 1 and A 0,e = C e . A 2,g gives no information about entanglement, just nonclassical effects in the field, as it only involves field excitation. For weak fields D 2 is exactly A 2,g . The entanglement shows up in the value of A 1,e ; if this value does not satisfy A 1,e = D 1 C e = A 0,e A 1,g , then it is not possible to write the state as a product state.
In the presence of a non-zero vacuum contribution (as any real quantum state will have), one can learn nothing about entanglement simply by measurement of one-excitation amplitudes or probabilities. For example, the state |0, g + α(|1, g + |0, e ) is entangled, but only if one is certain that the probability amplitudes for higher excitation are truly zero. A state of the form |0, g + α(|1, g + |0, e ) + O(ǫ 2 ) cannot be said to be entangled without information on the relative size of the probability amplitude A 1,e . Measurement of one-excitation amplitudes conditioned by a previous measurement can yield information about entanglement. This can be accomplished by utilizing cross-correlation functions. A first important conclusion out of this study is that a measure of the zero time cross correlation between the atom and the field, as well as the mean transmitted and fluorescent intensities yields a measure of entanglement in the weak field limit.

III. ENTANGLEMENT FOR WEAK EXCITATION
Equation (22)  the concurrence is: This section analyzes the sensitivity of the concurrence to the different parameters that appear in Eq. (26), while trying to give physical reasons for their influence on the entanglement. Despite the fact that the rates of decay could be the same through the two reservoirs, spontaneous emission (γ) reduces entanglement more than cavity loss (κ). This is due to the fact that a γ event (spontaneous emission) must come from the atom, while a κ event (cavity transmission) could come from either the drive or a photon emitted by the atoms into the cavity mode. A spontaneous emission event unambiguously leaves the atom in the ground state, and the system wavefunction factorizes. Fig. 2 shows a remarkable result in the entanglement of the system as a function of the three rates in the problem. There is an optimal value for the coupling constant g given a set of dissipation rates κ, γ. For many interesting cavity QED effects, stronger coupling is generally better, such as the enhancement of the spontaneous emission by a factor of 1 + 2C 1 = 1 + 2g 2 /κγ (this formula strictly holds only in the bad cavity limit κ >> g, γ).
However, here increasing the coupling of the atom and field mode eventually decreases the amount of entanglement. To explain this it is instructive to recall that the concurrence C = |2αβ(q − 1)|, where α is the mean cavity field, and β = −gα/γ is the mean atomic dipole. As the coupling g increases, for a fixed weak driving field ǫ, the intracavity field α = ǫ/(κ+2g 2 /γ) decreases. The intracavity field is the sum of the driving field in the cavity ǫ/κ, and the field radiated by the atom, (−2C 1 /(1 + 2C 1 ))ǫ/κ, the minus sign resulting from the fact that the radiated field is π out of phase with the driving field on resonance. We see that as g and C 1 increase, the intracavity field decreases. This means that the steady-state wavefunction has a larger vacuum component, and consequently less entanglement. Another way to view this is that the cavity enhancement of the spontaneous emission rate means a larger loss rate for the system as the coupling increases, which is bad for entanglement.
There are two interesting limits on this Eq. (27) for the parameter f = g/γ. If f ≫ 1, the steady state wavefunction is approximately |ψ ss = |0 (|g − f α|e ) which is a product state. Also, if f ≪ 1, the steady state wavefunction is approximately |ψ ss = |g (|0 + α|1 ) which again is a product state. To have entanglement between the atom and cavity mode, we must have the parameter f ≃ 1, so as to prepare a steady state wavefunction of the form |ψ ss = |0g + α(|1g − |0e ) = |0g + α|− , a mixture of the vacuum with a small entangled state component.
The decrease of the prefactor |αβ| is the dominant reason why the concurrence decreases with increasing g for large coupling. Close inspection of Fig. 2 also shows that there is an optimal cavity loss rate κ for entanglement for a fixed g and γ. This is a result of reaching a maximum in the population of the states different from the vacuum (Eq. (5)). Our results here are consistent with the numerical results of Nha and Carmichael [5].
When the system is driven off resonance, its response is typically characterized by transmission and fluorescent spectra [12,13]. Although these are important probes of the system, they do not, in this limit, carry information about the entanglement, since they are derived from only the one-excitation amplitudes. The concurrence as a function of the detuning of the driving laser shows that the steady state entanglement decreases typically by a factor of 1/∆ 3 for large detuning, where ∆ = (ω − ω l ) with ω the resonant frequency of the atom and cavity, and ω l the frequency of the driving probe laser. But in the case where g is larger than κ and γ, the response is maximized at the vacuum-Rabi peaks [14]. Figure 3 shows a contour plot of C for parameters in the regime of cavity QED where the two decay rates are similar: 2κ/γ = 1.0. The concurrence increases with increasing g on resonance up to a saddle point, and then decreases. However the entanglement persists for detunings on the order of g, the approximate location of the vacuum-Rabi peaks in the spectra of the system.
It would seem that by continuing to tune to a vacuum-Rabi peak as g increases, it would be possible to maintain the entanglement, but Fig. 3 shows that this is not the case. Rather, as argued (for the on-resonance case) above, the crucial parameter for maximizing entanglement is f = g/γ ∝ 1/ √ n sat , where n sat = γ 2 /8g 2 is the saturation photon number. This is the dependence on the nonlinearity of the atomic system. Recall that, if these were two driven coupled harmonic oscillators, q = 1 and there would be no entanglement. A nonlinear interaction between the two harmonic oscillators would be needed to entangle them, as in the signal and idler modes in optical parametric oscillation. This nonlinear interaction would generate two-mode squeezing, which could be measured by homodyne detection of mode A(B) conditioned on detection of a photon in mode B(A), just as squeezing in one mode can be detected via conditioned homodyne detection of a mode based on a photodetection from that mode [15,16]. The nonlinearity of the two-level atom is needed to generate two-mode squeezing and entanglement between the atom and the cavity field. Even though the driving field is weak and the atom never nears saturation, there can only be entanglement with a linear atom-field coupling if the atom has a nonlinear response, as two-level atoms do.
The concurrence shows its sensitivity to different parameters. Fig. 4 shows a contour plot of C versus g/γ and ∆/γ for a case where the cavity decay rate is larger than the spontaneous excitation. The previous section shows that the concurrence is related to the rate of single photon counts out of the cavity or in fluorescence and to the rate of coincident counts from the cavity and fluorescence. These are the quantities associated in quantum optics with correlation functions, first introduced by Glauber [17,18,19,20]. Generally these correlation functions involve comparing a field (intensity) of one mode with the field (intensity) of the same mode at a later time (or different spatial location), with some exceptions [21,22,23,24,25,26]. However, entanglement in cavity QED has two components: atom and cavity mode. It is natural to look at cross correlations between the cavity mode and the fluorescent light that falls in the mode of the detector.
Consider a general cross-correlation function for two-modes of the electromagnetic field: with f 1 and f 2 well behaved functions, in the sense of a convergent Taylor series on the Hilbert space of interest. If |ψ is a product state, the correlation function G(a, b) factorizes and then is unity. If it is not a product state, then this will manifest itself in a non-unit value for the normalized cross-correlation functions.
The simplest cross correlation function to consider is g T F (0). This could be obtained by measuring the visibility of the fringe pattern formed by interfering the transmitted and fluorescent light. For the weakly driven cavity-QED system, this is so to lowest order, there is no information in this correlation function about entanglement.
To obtain information about entanglement the correlation function has to probe the two-excitation part of the state. A possibility to do this is the intensity cross correlation: This normalized correlation function is directly related to the coefficient of double excitations (See Eqs. (5), (8), (13)). If q = 1 then g T F (0) = 1 and there is no entanglement; so a non-unit value of q indicates entanglement. Using second-order intensity correlations has been proposed in the context of entangled coherent states by Stobińska and Wódkiewicz [27].
The cross correlation function g (2) T F (0) contains information about the average photon number in coincidence with a measurement of the fluorescence relative to the average photon number in the absence of any interrogation of the fluorescence. g (2) T F (0) − 1 = q 2 − 1 is an indicator of entanglement.
A way to measure q directly utilizes a field-intensity correlation function h θ (τ ) [28], that can be implemented as a homodyne measurement conditioned on the detection of a fluorescent photon, So h T F θ=0 (0) − 1 = q − 1 is also an indicator of entanglement in this system. What makes this measurement possible experimentally is the conditioning that selects only times when there is a fluctuation and the rest of the time (when the vacuum is present) no data is collected [16]. For one mode, the homodyned transmitted field conditioned by detection of a photon from that mode, is a measure of squeezing in that mode [28]. A homodyne measurement of the transmitted field conditioned by detection of a fluorescent photon is a measure of the two-mode squeezing, with the cavity field and atomic dipole as the two components. Generally, two-mode squeezing is an indicator of entanglement between the two modes. Gea Banacloche et al. explored this correlation function in a different regime of cavity QED and found it to be a witness of the dynamics of entanglement [29]. There is a particular form of the Schwartz inequality that must be satisfied for a classical field for the specific case of the system we are considering here: (g T F (0) − 1) 2 ≤ |(g T T (0) − 1)(g Here T T and F F denote zero delay intensity correlations for the transmitted and fluorescent fields respectively. In the one-atom limit, g F F (0) = 0, and g T T (0) = q 2 p 2 , so this inequality becomes |q 2 − 1| 2 ≤ |q 2 p 2 − 1| which depends on q, but also on the parameter p (Eq. (12)), which can be varied independently. There is no one-to-one relationship between Schwarz inequality violations and entanglement (by this measure) in this particular system.

V. CONCLUSION
We find that entanglement in weakly driven cavity QED is characterized by comparison of two-excitation probability amplitudes to single excitation amplitudes, in particular the amplitude involving one excitation in each subsystem. It is necessary to have a small saturation photon number to enhance the nonlinear response which generates a larger entanglement. But this is true only to a point. We find the maximal entanglement for small κ and when g/γ is on the order of unity. This stems from the dual role of the coupling g. It couples energy into the atom, but due to cavity enhanced spontaneous emission, it can also channel energy out.
Increasing γ decreases the entanglement, and this can be explained in terms of the effect of the two decay processes on the system. If we detect a fluorescent photon we know it has come from the atom, and the atom is in the ground state. If we obtain a transmitted photon, it could have been emitted from the atom into the cavity mode, or just be a driving field photon that has passed through the cavity without interaction with the atom. It is the interference of these two indistinguishable processes that leads to nonclassical effects in the transmitted field.
We have found a variety of cross-correlation functions that are indicators, or witnesses, of entanglement in this system. One can learn nothing about the entanglement by examining only one-or two-excitation amplitudes separately. In particular we find that a measurement of two-mode squeezing, or a homodyne measurement of the transmitted field conditioned on the detection of a fluorescence photon is directly proportional to the entanglement calculated via the reduced von Neumann entropy. Further work remains to generalize this approach to situations with higher drives, but the general approach of looking at entanglement together with the specific correlation function to measure gives physical insight into this problem.
We would like to thank J. P. Clemens for fruitful discussions related to the topic of this