New Spin Squeezing and Other Entanglement Tests for Two Mode Systems of Identical Bosons

For any quantum state representing a physical system of identical particles, the density operator must satisfy the symmetrisation principle (SP) and for massive particles also conform to super-selection rules (SSR) that prohibit coherences between differing total particle numbers. Here we consider bi-partitite states for systems of massive bosons, where both the system and subsystems are modes (or sets of modes), particles being associated with differing mode occupancies. To define non-entangled or separable states the subsystem density operators are also required to satisfy the SP and conform to SSR, in contrast to some other approaches. Whilst in the presence of these additional constraints some of the known entanglement criteria are still valid, we present a detailed treatment to establish the validity of known criteria and derive new ones.

However, another spin squeezing entanglement test for Bose-Einstein condensates involving the variance in S z being less than the sum of the squared mean values for S x and S y divided by the boson number was based on a concept of entanglement inconsistent with the SP, and here we present a revised treatment which again leads to spin squeezing as an entanglement test. 2 Author to whom any correspondence may be addressed.

Introduction
Since the work of Einstein et al [1] on local realism, the famous cat paradox of Schrodinger [2] and the derivation of inequalites by Bell [3] and others [4], entanglement has been recognised as being one of the essential features that distinguishes quantum physics from classical physics.In macroscopic systems entanglement is of particular importance as it blurs the boundary between the classical and quantum worlds.This paper considers tests that are sufficient (though not necessary) to confirm such entanglement.
One way to detect macroscopic entanglement is by applying the so called spin squeezing inequalities [5] to a large number of particles.The importance of spin squeezing in quantum metrology was emphasised by Kitagawa and Ueda [6] and it has been demonstrated that such states beat the standard quantum limit in interferometry [7].For spin angular momentum operators Ŝx , Ŝy , Ŝz spin squeezing of Ŝx with respect to Ŝy is defined as ∆ Ŝ2 x < where ∆ Ŝ2 η is the variance of the spin operator Ŝη (η = x, y, z) Analogous definitions apply to other pairs of spin operators, or to a spin operator and any of its perpendicular components.For N spin-1/2 distinguishable particles each in the same coherent state, correlations between the spins resulted in spin squeezing in the total spin components [6], as it allows the fluctuations in one direction perpendicular to the spin component to be reduced.In such systems of distinguishable particles Sorensen et al [8,9] showed that spin squeezing required entanglement of the N particles.Here however, we consider N identical bosons occupying two modes, where the Schwinger angular momentum operators for two modes, A notable area in which the detection of macroscopic entanglement can take place is that of trapped ultra-cold gases, where entanglement may be present over micron scales.Moreover, since non-linear interactions between particles are required to generate spin squeezing [6] , ultra-cold gases are an ideal test bed for witnesses based on spin-squeezing inequalities due to the tunable atomatom interaction.In addition to achieving interferometry beyond the standard quantum limit, existing experiments have claimed to demonstrate macroscopic entanglement in ultra-cold gases via witnesses based on spin squeezing inequalities [10,11,12].
However, when detecting entanglement in systems of identical massive bosons with such inequalities care must be taken for two reasons.Firstly, entanglement is a relative concept that crucially depends on which sub-systems are being considered.A quantum state may be entangled for one choice of the subsystems but may be non-entangled if another choice of sub-systems is made.For distinguishable particles, the subsystems are usually individual particles, with internal (spin, polarization) or external (position, momentum) degrees of freedom possessing the quantum correlations.However, for identical particles the individual atoms cannot be distinguishable subsystems, and entanglement will refer to the quantum field modes that the particles may occupy [13].The modes are orthonormal single particle states, which may be localised in different spatial regions or may be delocalised where the opposite applies.The system and sub-systems are now modes, and cases with differing N are just different quantum states of the same system.Note however that a different concept of entanglement -particle entanglement -has also been applied to identical particle systems [14,15].Here we formulate spin squeezing entanglement witnesses for massive identical bosons as above in Eq.(1) via second quantization, so the symmetrization principle applies to all the identical bosons.
In the present paper we follow the approach of Werner [16] and define separable states as those which can be prepared via local operations (on the subsystems) and classical communication processes (LOCC).A relevant example of such states occurs in ultra-cold atom experiments, where mesoscopic ensembles can be prepared with definite particle numbers in optical lattice sites.This means that for bipartite systems the density operator for separable states can be written in the form where sub-systems A, B have been prepared in correlated quantum states ρ A R and ρ B with a probability P R .Entangled states are those which are not of this form.This definitition can be straightforwardly extended to multiple modes.
It is important to realise that the terms in ρsep have a physical meaning -the sub-system quantum states ρ A R and ρ B R must correspond to physical states that can actually be prepared in A, B considered as stand-alone systems.The subsystem states before the preparation begins are unimportant, apart from being a separable product of sub-system states -these might be easily preparable subsystem lowest energy pure states.This point of view regarding the nature of separable states is not new, see for example [17,18,19].However, other papers focus on the mathematical form for ρ and do not require LOCC preparation for defining what they refer to as separable but non-local states [20,21].In the latter papers the ρ A R and ρ B R are not required to represent physical states that could be prepared in the isolated sub-systems.However, for LOCC based separable states the joint probability for measurements of physical quantities Ω A , Ω B for the sub-systems resulting in eigenvalues λ Ai , λ Bj will be of the form P A,B (i, j) = R P R P A R (i)P B R (j), where Π Ai , Π Bj are the projectors onto subspaces with eigenvalues λ Ai , λ Bj , and the factors give the probabilities that the measurement results λ Ai , λ Bj occur if the sub-systems A, B are in states ρ A R and ρ B R .If the ρ A R and ρ B R did not specify possible physical states for the sub-systems, then the probabilities P A R (i) and P B R (j) would not have a physical meaning.This property of separable states in which the joint measurement probability is determined from measurement probabilities of physically possible separate sub-system states, is the key feature whose absence in non-separable states led Schrodinger to refer to such states as entangled states.Note that after preparation, a separable state may evolve into an entangled state, such as when interactions between the sub-systems take effect.
In addition, for systems of identical bosons we argue that to define the separable states (and hence the entangled states) one should take into account the presence of the particle-number superselection rule (SSR) not only at the global level -as required for any quantum state, but also crucially at the level of the local subsystems in the case of separable states.Simply put, the global SSR means that for any quantum state of massive bosons, ρ, should commute with the total particle number in the system [ρ, N ] = 0, and hence no coherences exist between Fock states with differing N .For the separable states given by Eq. ( 2) for two modes, this means in addition that for physical sub-system states, ρX R for X = A, B should also commute with the local particle numbers -thus ρX R satisfy [ρ X R , NX ] = 0. Note that both ρ and the ρX R may be mixed states, with statistical mixtures of particle number states.This additional local particle number SSR restriction leads to further tests of entanglement based on spin squeezing which are radically different to some of the well known entanglement witnesses.As we shall see below, a similar situation applies to other entanglement tests.In this letter, following an overview of the pertinant concepts, we will investigate how the spin squeezing and other inequalities are modified to detect entanglement when local particle number superselection rules are enforced for separable states.

Global and local SSR
The particle-number SSR is a fundamental constraint for systems of massive bosons.It occurs because physically realizable processes cannot create states which are coherent superpositions of different Fock states.Thus, all global operators should commute with the total particle number operator.For states that are separable, the sub-system states are also required to be physical states, so analogous constraints must also apply to the sub-systems, and hence the sub-system states also have no coherences between Fock states with differing local particle number.This feature is required for separable states in identical particle systems if -as described above -joint measurements are to be based on physical probabilities for such measurements on the sub-systems.Indeed, the SSR also limits the allowed local operations on the modes (e.g.unitary operations of the form U A ⊗ U B ) -they too have to commute with the local particle numbers.Enforcing the local SSR leads to a more restricted set of separable states than would apply if they are ignored, which can be expected to modify an entanglement witness that was based upon the non-restricted definition of separability, in which the ρ A R and ρ B R are not required to satisfy the local particle number SSR.Other authors have applied the SSR at a local level, see [20,17,22,23,24,25,26] for instance.A similar separable state was used [27] to show that the visibility of interference fringes between two spatial modes is an entanglement witness for modes of massive particles.Non local SSR conforming coherent superpositions of number states (such as Glauber coherent states) were once thought necessary for describing interferometry and coherence effects in Bose-Einstein condensates, but it is now recognised [28,29] that valid descriptions of these effects can be based on Fock states.
We will now make a brief comment about the connection between reference systems and superselection rules.It has long been known [30] that by using a suitable reference system one can, at least in principle, perform protocols (e.g.Ramsay interferometry, dense coding, Bell inequalities) whose description involves coherences between different particle numbers [31,32,24].However, this apparent inconsistency can be resolved [33].The key point is that the system states involved in a process are described differently by observers possessing different phase reference systems.If one observer describes the state as in Eq. ( 2) in terms of a reference system where the sub-system density operators σ X R do not conform to the local particle number SSR, then a second observer without access to this reference system would also describe a separable state, but now with the transformed sub-system density operators ρX R compatible with the local subsystem SSR.A more detailed discussion in terms of the U (1) symmetry group is presented in [22] (see also [34], Appendix 4).In this paper, we work from the viewpoint of the second observer.Phase reference systems such as Bose-Einstein condensates with large boson numbers -described by the second observer as being in statistical mixtures of Glauber coherent states with large fixed amplitudes and all phases having equal weight, and which satisfies the SSR for the BEC mode -are often involved in protocols such as Ramsay interferometry, dense coding, Bell inequalities.Bosons from the BEC reference interact and exchange with the primary system that is involved, for example two mode [31,32] or four mode systems [24].From the point of view of the first observer, the evolution from each large amplitude Glauber coherent state of the BEC times the initital state of the primary system can create states of the primary system that violate the SSR.However, from the point of view of the second observer, the evolved state of the primary system conforms to the super-selection rule following averaging over the BEC phases.In our approach, separable states are described by an observer with different phase reference systems with unknown phases for the different sub-systems -as is appropriate for sub-systems that are separate.On the other hand, entangled states are described by an observer with a single phase reference system (albeit one where the phase is unknown), so that entangled states are global particle number SSR compliant.This is appropriate for composite systems where the sub-systems are combined into an entangled state on which measurements for overall system operators can also be performed.In the present paper the protocol is spin squeezing and we will describe this in terms of reference systems so that the state Eq.( 2) satisfies [ρ X R , NX ] = 0 for X = A, B, which can be achieved in ultra-cold gas experiments.

Spin squeezing requires entanglement
Under the global and local particle-number superselection rule requirements for physical states of indistinguishable particles, the spin squeezing inequalities Eq. ( 1) are satisfied only if entanglement is present between the modes A and B. The proof is as follows.Using realization, ρR , of the separable state Eq. ( 2), the operators in the variance are Ŝ2 x R = Ŝ2 and Ŝx R = Ŝy R = 0, with n X R = NX R .Then, since for any operator Ω the inequality ∆ Ω2 ≥ R P R ∆ Ω2 R applies [35], the inequality ∆ Ŝ2 x

holds and likewise for ∆ Ŝ2
y .For SSR-restricted separable states, Eq. ( 2), the expectation values of the z spin component 1  2 This shows for general separable states in Eq.( 2) over two modes, A and B in which the local mode SSR applies, the spin variances in both the x and y directions satisfy

∆ Ŝ2
x ≥ and hence there is no spin squeezing of Ŝx with respect to Ŝy or vice versa It is straightforward to show that analogous results apply to the other pairs of spin components S y , S z and S z , S x or to each spin component and any perpendicular spin component.Thus we have now shown that spin squeezing in any spin component is a test of entanglement -entanglement being defined here in terms of separable states satisfying the local particle number SSR.There are, of course, entangled states which do not show spin squeezing.One example is the famous Schrodinger cat like N OON state, (|N, 0 + |0, N )/ √ 2, which describes all the particles in mode A and none in mode B superposed with no particles in mode A and all particles in mode B. For large N such states are however notoriously difficult to create in experiments.The N OON state is global particle number SSR compliant, and its density operator is the same for any observer with a different overall phase reference system.Note that for an observer with separate unrelated phase reference systems for the two modes, the state is seen as a separable state, whose density operator is a 50:50 mixture with one state having N bosons in mode A and none in mode B, the other having no bosons in mode A and N in mode B.

Other spin squeezing tests and local SSR separable states
There are several entanglement tests, based both on spin squeezing and other protocols that can be used to detect entanglement in systems, where the definition of entanglement is based on separable states which do not satisfy the physically based particle-number superselection rule.Of course any test that is derived for arbitrary ρX R must also demonstrate the present entanglement based on [ρ X R , NX ] = 0. Several such tests are discussed elsewhere [34].However, it is interesting and necessary for experiments to consider whether additional tests arise when they are re-derived for systems of identical massive bosons based on the physical definition of a separable state (2), where the local SSR is satisfied.We will focus on certain commonly used inequalities in this paper, although other tests also exist (see [36,37]).

Hillery et al 2006 spin squeezing test
A paper by Hillery and co-workers [38] derives an entanglement witness (Eq.( 4)) based upon spin squeezing, based on general (non-SSR-restricted) separable states of the EM field.It has been used to detect entanglement in e.g. the following proposals for trapped ultra-cold gases [39,40].We confirm that the main result of [38] is still valid for the SSR restricted state ( 2), but provide an example where the Hillery test based on the sum of J x , J y variances for rotated spin components fails to detect entanglement in the relevant modes, whereas a simple spin squeezing test (as in Eqs.( 3)) involving a rotated component J y is satisfied.
Hillery et al [38] show that for general separable states To obtain this result it is found that for a product state ρR that ∆S , where for non SSR compliant ρX R the terms â R and b † R are not necessarily zero.Noting that | â | 2 ≤ Na and likewise for mode b, it follows that ∆S 2 x R + ∆S 2 y R ≥ 1 2 ( N R .From ∆ Ω2 ≥ R P R ∆ Ω2 R the result in Eq.( 4) then follows.For SSR-restricted separable states however, the result ∆S N .If this applies then the state cannot be separable.
Note that we have previously shown that for local SSR compliant separable states that ∆S x + ∆S 2 y < | Ŝz | is not a valid entanglement test because there are no quantum states where this is true, as noted in [38].
Although it may be thought that because N would be more likely to demonstrate entanglement in the modes A, B than the tests obtained in this paper such as showing ∆ Ŝ2 x < 1 2 | Ŝz | or ∆ Ŝ2 y < 1 2 | Ŝz |, this is not always the case.The relative phase eigenstate [41,42] is an entangled pure state for N bosons defined by N  2 , θ p = ( , where θ p = p(2π/(N + 1)) with p = −N/2, −N/2 + 1, .., +N/2 specifies the relative phase.In terms of spin components S x , S y , S z the covariance matrix [43] for spin fluctuations is non-diagonal and spin squeezing does not occur.As the variances are such that ∆ Ŝ2 x + ∆ Ŝ2 y = ( 1 6 − π 2 64 )N 2 and which exceeds 1 2 N , the Hillery test for entanglement also fails even though the state is entangled for modes A, B. Principal spin components J x = S z , J y = (sin θ p ) S x + (cos θ p ) S y , J z = −(cos θ p ) S x + (sin θ p ) S y are obtained by a rotation, and for these the covariance matrix is diagonal with ∆ )N 2 for large N .The mean values are J x = J y = 0, J z = − π 8 N (see [42] for details).The principal spin components are related to annihilation operators c, d for new modes C, D via expressions of the form , and the relative phase state can be rewritten as a linear combination of Fock states for the new modes |N/2 − l C ⊗ |N/2 + l D with l = −N/2, ., +N/2.Note that there must be terms with differing The relative phase state is therefore an entangled state for the modes C, D. In the relative phase state J y is squeezed with respect to J x so the spin squeezing test for entanglement of modes C, D based on separable states consistent with the local particle number SSR is satisfied.However, for the Hillery test ∆ J 2 x + ∆ J 2 y ≈ 1 12 N 2 , which considerably exceeds 1  2 N .Thus the Hillery test for entanglement fails.

Sorensen et al 2001 spin squeezing test
Care must be taken when applying a spin squeezing entanglement based witness derived in an early paper by Sorensen et al [8] for systems of N distinguishable two-level particles to the situation when the particles are identical [15].The witness shows that entanglement exists when the state satisfies the following spin squeezing inequality The proof is based on writing the separable state density operator in the form ρsep = R P R ρ1 R ⊗ ρ2 R ⊗ρ 3 R ...where ρi R is the density operator for the ith distinguishable particle, whose internal states are |φ a (i) and |φ b (i) .The density matrix for ρi R is a 2 × 2 matrix and the spin operators are defined by expressions such as S x = i (|φ b (i) φ a (i)| + |φ a (i) φ b (i)|)/2 etc., which are sums over the particles.The hermiticity, positiveness, unit trace for all the density matrices lead to ξ 2 ≥ 1 for these separable states of two level distinguishable particles.This result cannot just be applied to identical bosonic particles without further development, since the Sorensen separable state density operator does not satisfy the symmetrization principle.Also Benatti et al [44] have shown that this inequality diverges for two mode states where the local particle-number superselection rule is enforced, as is easily seen by noting that S x = S y = 0 for a physical separable state, Eq.( 2) involving single modes as subsystems.Just considering N boson states for a total of two modes is not adequate.One way to revise the Sorensen et al result to an equivalent theory now based on mode entanglement involves sub-systems each consisting of two modes, which we may list as αk, where α = a, b and k lists modes with the same α.The k may correspond to spatial modes localised on different lattice sites.For identical particles occupying these modes the spin operators x etc, where the S k γ and S γ satisfy the usual commutation rules.The sub-system density operators ρk R now refer to two mode sub-systems and can be made consistent with the local particle number SSR by requiring that [ρ k R , N k a + N k b ] = 0.The symmetrization principle automatically applies for ρsep = R P R k (ρ k R ) ⊗ in this second quantization treatment.If in addition, we require that the ρ(k) R are density operators for a single boson, then the bosons in different k modes are effectively distinguishable.The proof in Sorensen et al [8] then applies, noting that for two mode subsystems S k γ R and hence S γ (γ = x, y, z) are not necessarily zero.States of this type are routinely created in experiments, for instance in [45], where the gas is frozen deep in the Mott regime with unity filling and each atom has two accessible internal states.The sites then act like distinguishable qubits.The physicality of a separable state is thus guaranteed and there may be entanglement between the modes within each subsystem.A similar proof extending the test of Sorensen et al to identical two level systems is given by Hyllus et al [15] based on a particle entanglement approach.In their approach bosons in different external modes (such as the different k) are treated as distinguishable and the symmetrization principle is ignored for such bosons.
It should be noted that for separable states where the ρX R are required to satisfy further conditions in addition to the local particle number SSR , the entanglement tests will differ from those where the additional conditions are absent.The requirement that the sub-system density operators ρk R are restricted to one boson states is an example of such an additional condition.

Non spin squeezing tests and local SSR separable states
In [38,46] it was found that for separable states based on arbitrary subsystem density operators ρX

Experimental considerations
A recent experiment [12], uses the inequality (5) to detect the entanglement in an ultra-cold gas.Here, despite their indistinguishability, they consider the particles themselves as the sub-systems and generate 'entanglement' via their internal degrees of freedom.Due to this, the inequality (5) would be valid as described above if distinguishability could be recovered.This while technically very difficult, can in principle be achieved.For instance, the gas could be frozen (without disturbing the internal states of the particles -into the Mott phase of an optical lattice so that the system behaves as a set of distinguishable qubits.Such a scheme has been considered above. On the otherhand, an earlier experiment [10] considers the entanglement of two spatial field modes.Here the application of inequality (5) should be reconsidered due to the fact that the subsystems are single modes (with no internal structure).
The experiments only involve a single test for entanglement, and it would be desirable to confirm entanglement via an independent test.

Conclusions
Only three ground breaking experiments [12,10,11] have shown spin squeezing in cold atom systems, from which the presence of entanglement is inferred.As such the field is still very much in its infancy.Our results will advance the field by allowing a greater understanding of the role of indistinguishability when detecting entanglement.A more extensive presentation of the work in this paper is in preparation [34].
) nbn is a test for entanglement.This test must also apply for separable states where the local SSR is satisfied.A particular case of the test is| âb † | 2 = | â †b | 2 > Na NbHowever, for separable states satisfying the local particle number SSR we can easily show that| âm ( b † ) n | 2 = 0. (7) Hence | âm ( b † ) n | 2 >0 is a test for entanglement based on local SSR compliant separable states.Since | âm ( b † ) n | 2 is merely required to be non-zero this test is easier to satisfy than the Hillery one based on Eq.(6).A particular case of the test is | âb † | 2 > 0. The last result also follows from S x = S y = 0 for SSR compatible states.Hence this is a simpler test for entanglement than | âb † | 2 > Na Nb and has been used instead to detect entanglement [27].