arxiv: v1 [quant-ph] 27 Sep 2016

Transcription

1 Separability conditions and limit temperatures for entanglement detection in two qubit Heisenberg XY Z models N. Canosa, R. Rossignoli Departamento de Física, Universidad Nacional de La Plata, C.C.67, La Plata (900), Argentina arxiv: v [quant-ph] 7 Sep 06 We examine the entanglement of general mixed states of a two qubit Heisenberg XYZ chain in the presence of a magnetic field, and its detection by means of different criteria. Both the exact separability conditions and the weaker conditions implied by the disorder and the von Neumann entropic criteria are analyzed. The ensuing limit temperatures for entanglement in thermal states of different XY Z models are then examined and compared with the limit temperature of the symmetrybreaking solution in a mean field type approximation. The latter, though generally lower, can also be higher than the exact limit temperature for entanglement in certain cases, indicating that symmetrybreaking does not necessarily entail entanglement. The reentry of entanglement for increasing temperatures is also discussed. PACS numbers: Pacs: Ud, a, 75.0.Jm I. INTRODUCTION Entanglement is one of the most distinctive features of quantum mechanics, representing the ability of composite quantum systems to exhibit correlations which have no classical analogue. Recognized already by Schrödinger [], it has recently become the object of intensive research due to the key role it plays in the field of quantum information [ 6]. Rigorously, a mixed state ρ of a bipartite system is said to be separable or classically correlated [7] if it can be expressedas aconvexcombinationofuncorrelated densities, i.e., ρ = ν q νρ ν A ρν B, where ρν A, ρν B are mixed states of each subsystem and q ν are non-negative numbers. Otherwise, ρ is entangled or inseparable. When separable, ρ satisfies all Bell inequalities as well as other properties characteristic of classical systems. A pure state ρ = Φ Φ is separable just for tensor product states Φ = φ A φ B, but in the case of mixed states, like thermal states ρ exp[ H/T], with H the system Hamiltonian, it is in general much more difficult to determine whether ρ is separable or not. Only in special cases, like a two-qubit or qubit+qutrit system, simple necessary and sufficient conditions for separability are known [8, 9]. Moreover, the entanglement of formation of amixedstate[0]hasbeen explicitlyquantifiedonlyfora two-qubit system []. Nonetheless, it is known that any mixed state becomes separable if it is sufficiently close to the fully mixed state [, 3]. For thermal states of finite systems, this implies that a finite limit temperature for entanglement [4], T e, will always exist such that ρ becomes separable T T e. It is then interesting to analyze if it is possible to estimate this temperature with simple separability criteria, and how it is related to the criticaltemperaturet c ofthe symmetry-breaking solution in a mean field type approximation, which is the conventional starting point for describing interacting manybody systems. Such solutions (i.e., like deformed or superconducting) normally reflect the presence of strong correlations and collective behavior. Theaimofthisworkistoexaminetheseissuesinasimple yet non-trivial model where the exact entanglement conditions and quantification can be easily obtained. For this purpose, we will consider a system of two qubits interacting through a Heisenberg XY Z Hamiltonian[5] in the presence of an external magnetic field. Interest in this model stems from the potential use of Heisenberg spin chains for gate operations in solid state quantum computers [6, 7]. The pairwise entanglement of thermal states of isotropic [8, 9] and anisotropic XY [0 ] Heisenberg models have accordingly been recently studied, and several interesting features have appeared already in the two-qubit case [], like the possibility of entanglement reentry for increasing temperatures or magnetic fields. We will first review the exact separability conditions for general mixtures of the eigenstates of arbitrary XY Z Hamiltonians, examining in particular thermal states and the possibility of entanglement reentry. We will also analyze the weaker conditions provided by the disorder criterion [3], which is the strongest one based just on the spectrum of ρ and one of its reductions, and is hence more easy to implement in general than other criteria. Violation of the disorder conditions also ensures distillability [4]. These conditions are here exact in the absence of a magnetic field. Although the disorder criterion admits a generalized entropic formulation[5], it is stronger than the von Neumann entropic criterion [6], based on the same information, whose predictions will also be analyzed. The ensuing exact and approximate limit temperatures for entanglement in thermal states of different XYZ models will then be examined. Finally, we will discuss the mean field (i.e., independent qubit) approximation for thermal states, with the aim of comparing the previous limit temperatures with the corresponding mean field critical temperature T c. It will be shown, remarkably, that for T > 0, symmetrybreaking is not necessarily a signature of entanglement, so that T c may be higher than T e, although it is usually lower. The model and methods are described in sec. II, while three different examples are analyzed in detail in sec. III. Conclusions are finally drawn in IV.

2 II. FORMALISM A. Model and separability conditions We will consider a Heisenberg XYZ chain [5] for two qubits in an external magnetic field b along the z axis. Denoting with S = s A +s B the total spin of the system, the corresponding Hamiltonian can be written as H = bs z v i s A i s B i (a) i=x,y,z = H z v + (Sx +S y ) v (Sx S y ), (b) where H z = bs z v z (S z /) and v ± = (v x ± v y )/. The ferromagnetic (antiferromagnetic) case corresponds to v i 0 ( 0), and the standard XY model to v z = 0. Its normalized eigenstates H Φ j = E j Φ j are given by Φ 0,3 = + +, E 0,3 = v z ±v +, Φ, = u± ++ u, E, = v z ±, () with = v +b /v, u ± = ±b/ and ±± ± ± the separable eigenstates of S z (standard basis). Thestates Φ 0,3 aremaximallyentangled,whereas Φ, are entangled for v 0, with concurrence v / (see Appendix). They become maximally entangled for b = 0, in which case the set of states () is just the Bell basis. We will first consider general statistical mixtures of the previous eigenstates, which can be written as 3 ρ = p j Φ j Φ j (3a) j=0 = 4 + S z S z +4 where p j 0, 3 j=0 p j = and i=x,y,z s A i sb i sa i sb i, (3b) S z = b (p p ), s A z s B z = (p +p ), s A i sb i = 4 [p 3 p 0 ± v (p p )], i = x,y, (4) with O TrρO. Eqs. (3) comprise standard thermal states as well as those arising in more general statistical descriptions [7, 8], and represent the most general two-qubit state with good permutational and phase flip symmetry U = e iπsz real in the standard basis. The two-site density matrix of an N qubit XYZ chain with cyclic boundary conditions is in fact also of this form[]. Exact separability conditions. For the state (3), they can be most easily determined with the Peres criterion [8], sufficient for two qubits [9], and can be cast as v p p p 0 +p 3, (5a) p 3 p 0 [(p +p ) b (p p ) ], (5b) or, in terms of the averages Si = sa i sb i +/, as S x S y S z, S x +S y [ S z S z ], (6a) (6b) imposing bounds on the averagesof the last two terms in (b). If ρ is entangled, only one of Eqs. (5) is violated, and its concurrence is given precisely by the difference between the left and right hand sides of the broken inequality (see Appendix). The entanglement arises essentially from one of the states Φ, ( Φ 0,3 ) if (5a) [(5b)] is broken. Eqs. (5) are always satisfied if p j /4 (4 ) j, i.e., if ρ is sufficiently close to the fully mixed state. If b = 0, ρ is diagonal in the Bell basis and Eqs. (5) reduce accordingly to p j / j [6], while Eqs. (6) to S + S i for i = x,y,z, as S z = 0. Disorder and entropic separability conditions. The disorder criterion [3] states that if ρ is separable, ρ is majorized by the reduced densities ρ A,B Tr B,A ρ, which meansthat ρis more mixed(i.e., disordered) than ρ A, ρ B. In a two qubit system, this implies that the largest eigenvalue of ρ should not exceed that of ρ A and ρ B, which is in general a necessary condition that becomes sufficient when ρ is pure or diagonal in the Bell basis [3, 5]. For the state (3), ρ α = /+ S z s α z for α = A,B, and the disorder criterion leads to the inequalities p j [+ b (p p ) ], j = 0...,3, (7) which in terms of total spin averages can be recast as S x S y S z [+ S z / S z ], S x +S y S z + S z. (8a) (8b) Eqs. (7) or (8) are clearly less stringent in general than Eqs. (5) or (6), but become exact for b = 0 ( S z = 0), i.e., when ρ is diagonal in the Bell basis. The standard entropic criterion [6], based on the von Neumann entropy S (ρ) = Trρlog ρ, states that if ρ is separable, S (ρ) S (ρ α ) for α = A,B. Although exact for pure states (in which case S (ρ) = 0 and S (ρ A ) = S (ρ B ) is just the entanglement of ρ []), for mixed states it is in general weaker than the disorder criterion [5], except when both ρ and ρ α have rank two. Fig.depicts, forp = 0and b/v =, the regionswhere the state (3) is entangled and where entanglement is detected by the disorder and the standard entropic criteria. Standard thermal state and entanglement reentry. For ρ = exp[ βh]/trexp[ βh], β /T > 0, (9) i.e., p j e βej in (3a), Eqs. (5) become v eβvz sinh β cosh(βv + ), (0a) e βvz sinh βv + [+ v sinh (β )], (0b) and determine a finite limit temperature for entanglement T e, such that they are satisfied T T e. Nonetheless, entanglement, as measured by the entanglement of formation or concurrence, may not be a decreasing function of T for T < T e when the ground state is less entangled than the first excited state [8 ], and even entanglement vanishing plus reentry may occur [], as discussed below. Note that the spectrum of H, and hence Eqs. (0), do not depend on the signs of b and v ±.

3 3 p 3 0 B B p d p c A p d p FIG. : Range of values of p and p 3 where the state (3) is entangled (shaded sectors), for p = 0, b/v =. Eq. (5a) is broken in sector A while Eq. (5b) in sectors B. Entanglement is detected by the disorder criterion in the black and dark gray sectors, and by the von Neumann entropic criterion just in the black sectors. They both coincide at the border of the triangle, where ρ has rank, and are exact for p = 0, where ρ is in addition diagonal in the Bell basis. Let us consider for instance a mixture of Φ and Φ 3 (p 0 = p = 0 in (3a)), which corresponds to the outer border in fig.. This state is separable just for p = p c (+v / ) /, with Eq. (5a) [(5b)] broken for p > p c (p < p c ). Its concurrence is C(ρ) = p /p c. The state(9)willapproximatelybeofthisformforlowt ife ande 3 aresufficientlycloseandwellbelowtheremaining levels. Hence, if E < E 3, C(ρ) will initially decrease as T increases from zero, vanishing at the temperature T r = (E 3 E )/ln[ /v ], () where p = p c, but will exhibit a reentry for T > T r, with Eq. (0a) [(0b)] broken for T < T r (T r < T < T e ). Due to the remaining levels, C(ρ) will actually vanish in a small but finite temperature interval around T r (see case 3 in sec. III). When Φ becomes separable (v /b 0), p c and T r 0, whereas when it becomes maximally entangled (b/v 0), p c / and T r, so that no reentrytakes place in this limit. In constrast, if E 3 < E, p < / p c T 0 and no reentry or enhancement of C(ρ) can take place. Nor can it occur for a mixture of Φ 0 and Φ 3 or Φ and Φ, since they are separable just for equal weights, as seen from Eqs. (5). For the state (9), the disorder conditions (7) become ( b )eβvz sinh β cosh(βv + ), (a) e βvz sinh βv + cosh(β )+ b sinh β, (b) and lead to a lower limit temperature for entanglement detection, T d e T e, with T d e = T e just for b = 0. The entropic criterion leads to an even lower limit temperature T s e T d e. The reentry effect cannot be detected by the disorder (and hence by the entropic) criterion. Violation of Eqs. (7) requires p j > / for some j, so that in the thermal case just the entanglement arising from the ground state can be detected. For a mixture of Φ and Φ 3, Eqs. (7) are broken just for p > p d = ( b/ ) or p < p d = (+ b/ ) (see Fig. ), which does not allow to detect the reentrywhen E < E 3 since p d /. B. Symmetry-breaking mean field approximation The thermal state (9) represents the density operator that minimizes the free energy F(ρ) H TS(ρ) = Trρ[H +T lnρ]. (3) In a finite temperature mean field or independent qubit approximation, Eq. (3) is minimized among the subset of uncorrelated trial densities, given in this case by ρ mf = ρ A ρ B, (4) with arbitrary ρ A, ρ B, obtaining thus an upper bound to the minimum free energy. The only way such an approximation can reflect entanglement is through symmetry breaking: the optimum density that minimizes F(ρ mf ) may break some of the symmetries present in the Hamiltonian H, and become degenerate. In these cases a critical temperature T c will exist such that the optimum density becomes symmetry conserving for T T c. At T = 0, symmetry breaking implies entanglement if the ground state of H is non-degenerate, since for pure states separability corresponds to an uncorrelated density. However, this is not necessarily the case for T > 0, where symmetry-breaking just indicates, in principle, that the true thermal state is not uncorrelated. On the other hand, entanglement does not necessarily imply symmetry-breaking either, both at T = 0 or T > 0, as correlations need to be in general sufficiently strong to induce a symmetry-breaking mean field [9]. The densities ρ α, α = A,B, can be parameterized as ρ α = exp[ βλα s α ] Trexp[ βλ α s α ] = + sα s α, (5) s α = Trρ α s α = λα tanh[ β λα ]/ λ α, so that Eq. (4) corresponds to an approximate independent qubit Hamiltonian h = α λα s α. Minimization of F(ρ mf ) with respect to λ α leads then to the selfconsistent equations (see for instance [30]) λ α i = H mf s α, i = x,y,z, (6) i where H mf = Trρ mf H. A similar equation obviously holds for the n qubit case. In the case (), H mf = b S z i v i s A i sb i and Eqs. (6) become λ A,B i = bδ iz v i s B,A i. (7) Permutational symmetry will be broken if λ A λ B, and phase flip symmetry if λ α x 0 or λα y 0. The latter has to be broken in order to see any effect from the last twointeractiontermsin(b) atthemeanfieldlevel, since otherwise their mean field averages vanish. In such a case the sign of one of the λ α x (or λα y ) remains undetermined, giving rise at least to a two-fold degeneracy. For instance, in the ferromagnetic case v i 0, H mf is minimum for s A = s B and permutational symmetry needs not be broken. Hence, λ A,B = λ. Defining

4 4 v M = Max[v x,v y ], v m = Min[v x,v y ], a phase-flip symmetry breakingsolution with λ M = 0 and λ m = 0 becomes feasible and provides the lowest free energy if v M > v z and b < b c v M v z, provided 0 T < T c, with T c = v M χ/ln[ +χ χ ], χ b /b c <. (8) T c decreases as χ increases, with T c 0 for χ and T c v M( χ /3) for χ. This solution is insensitive to v m. As discussed in sec. III, T c is usually lower than T e, but can also be higher. For example, if b = 0 and v x > v y = v z > 0, T c = v x /, but the ensuing exact thermal state, diagonal in the Bell basis, is separable T > 0 (T e = 0), as the ground state is degenerate (E = E 3 v x /) and hence p j / j,t. III. EXAMPLES We now examine in detail the previous limit temperatures in three different cases. We set in what follows b 0, v ± 0, since the concurrence and limit temperatures are independent of their signs. ) v = 0, v + > 0 (XXZ model). The states Φ, are in this case separable, with = b in Eq. (). Entanglement can then only arise through the violation of Eq. (5b), i.e., Eq. (0b) in the thermal case, which is now independent of the magnetic field b. If v + > v z, the thermal state (9) will then be entangled for any b if T > 0, up to a limit temperature T e that is independent of b. However, the ground state is Φ 3 if b < b 0 v + v z and Φ if b > b 0, so that for b > b 0, ρ becomes entangled only at finite temperature T > 0, in agreement with Eq. () (T r 0 for v /b 0). On the other hand, if v + < v z, no entanglement occurs at any temperature. These features can be appreciated in Fig. for v z = 0 (XX model), where b 0 = v + and [0] T e = αv +, α = /ln[+ ].34. (9) The disorder criterion can now detect entanglement just throughtheviolationofeq.(7)forj = 0,3,i.e., Eq.(b) in the thermal case, which can occur only for b < b 0, i.e., when Φ 3 is the ground state. The entanglement arising for T > 0 when b > b 0 cannot be detected. In addition, the limit temperature Te d determined by (b) will depend on b, decreasing as b increases and vanishing for b b 0. Its behavior for v z = 0 is shown in the central panel of fig., where Te d T e [ b/( v + )] for b 0 while Te d (v + b)/ln for b b 0 = v +. Also shown is the concurrence at T = Te d (bottom panel), which is maximum at b = b c (where C[ρ(Te d )] /3) and decreasesas α ( )b/v + for b 0. The limit temperaturete s ofthe entropiccriterionis stilllower. For v z = 0andb 0,Te s 0.478v +, withc[ρ(te)] s 84. In this case the ground state critical field b 0 coincides with the mean field critical field b c. Hence, a stable symmetry breaking mean field solution is here feasible just C(ρ) Te/v+ C[ρ(TL)]. 0.9 v + /b= 0 T/v + exact disorder entropic mean field 0 4 v + /b FIG. : Top: The concurrence as a function of temperature for v = v z = 0 and indicated values of v +/b. Center: The corresponding exact limit temperature for entanglement (constant) and the limit temperatures below which entanglement is detected by the disorder and by the von Neumann entropic criterion. The critical temperature for the symmetry-breaking mean field approximation is also shown. Bottom: The concurrence at the previous limit temperatures. for b < b 0, with T c given by Eq. (8) with v M = v +. Since now [H,S z ] = 0, this deformed solution breaks the rotational invariance around the z axis and possesses, accordingly, a continuous degeneracy. As seen in Fig., for v z = 0 T c is much lower than T e, lying actually quite close to the entropic limit temperature T s e. For b 0, T c v + /, with C[ρ(T c )] 5. Note, however, that for b b 0, T c > T d e due to the logarithmic vanishing of T c in this limit, where C[ρ(T c )] /. ) v > 0, v + = 0. This case of maximum anisotropy (v x = v y ) represents, for v z = 0, the two-qubit version of the standard Lipkin model, widely employed in nuclear physics to test symmetry-breaking mean field based descriptions[9]. It describes the interplay between a single particle term bs z and a monopole interaction that induces a deformed mean field. The states Φ, are now entangled, whereas the states Φ 0,3 become degenerate. Hence, in the thermal case p 3 = p 0, and entanglement can only arise from the states Φ,, i.e., through the violation of Eq. (5a) (Eq. (0a) in the thermal case). This requires > v z, i.e., that Φ be the ground state. ThelimittemperatureT e determinedby(0a)depends now on the field b, with ρ entangled for 0 T < T e and C(ρ) a decreasing function of T, as seen in Fig. 3 for

5 5 C(ρ) Te/v C[ρ(TL)].5 0. v /b= 0 T/v exact disorder entropic mean field 0 4 v /b FIG. 3: Same quantities as in Fig. for v + = v z = 0 and different values of v /b. v z = 0. In this case, entanglement occurs b and T e = /arcsinh[ v ], (0) with T e αv [ + ( α/ )b /v ] for b 0. A remarkable feature is that as b increases, T e now increases, even though the entanglement of Φ decreases, since the energy gap between the ground and the first excited states increases. Moreover, for b, T e b/ln(b/v ), being then possible to make ρ entangled at any temperature by increasing the field. For p 0 = p 3, entanglement will be detected by the disorder criterion through the violation of Eq. (7) for j =,, i.e., of Eq. (a). For v z = 0, this will occur for any value of b but below the lower limit temperature Te d = /arcsinh[ b ]. () For b 0, T d e T e[ αb/( v )], with C[ρ(T d e )] ( )b/v. Eq. () is not a monotonous increasing function of b, being minimum at b.5v, but for b, T d e b/ln[4b /v ] T e /, becoming then also infinite in this limit. Hence, the emergence of entanglement for large fields is also detected (since it is a ground state effect) but above a higher threshold. Note also that T d e/t e / b, with C[ρ(T d e)] The limit temperature of the entropic criterion lies very close to T d e for b (as T d e/ 0 in this limit) but becomes smaller as b decreases, with T s e 0.478v for b 0. For v z = 0, a phase flip symmetry breaking mean field solution becomes here feasible only for b < b c = v. For b > b c, ground state correlations, though non-vanishing, are not strong enough to induce a symmetry-breaking mean field, so that the entanglement effect for large fields cannot be captured by the mean field. The permutationally invariant solution corresponds to λ x 0 and λ y = 0, sothat the criticaltemperature isgivenagainbyeq. (8) with v M = v. Hence, T c v / for b 0, lying again very close to T s e in this limit, while T c 0 for b b c, where C[ρ(T c )] / ) v + > 0, v > 0. This is the case with finite anisotropy γ = v /v + > 0, where entanglement vanishing plus reentry mayoccur ast increases. For v + > v z 0, the two lowest states are Φ and Φ 3, with E < E 3 for > v + v z, i.e., b > b 0 = Max[0,(v + v z ) v ]. For b above but close to b 0, Eq. (0a) [(0b)] will be broken for 0 T < T r (T + r < T < T e ), with T r < T + r. Hence, as T increases from zero, the concurrence will first decrease, vanishing for T [T r,t + r ], but will exhibit a reentry for T > T + r, vanishing finally for T T e. This behavior is depicted in Fig.4 for v z = 0 and γ = 0.7, where b 0 0.7v + and the reentry occurs for b 0 < b < b r.v +. For b close to b 0, Tr and T r + are practically coincident and equal to the value given by Eq. (), T r = ( v + )/ln[ /v ], becoming the difference exponentially small for b b 0 (T r + T r T re v+/tr ). For b > b r the reentry disappears and T e becomes the continuation of Tr, undergoing then a sharp drop at b = b r. For b, T e, as in case, while for b 0, T e 0.93v +. At fixed T < 0.93v +, entanglement vanishing plus reentry will then also occur as b increases. As discussed in sec. II, the disorder criterion cannot detect the reentry for increasing T. Instead, the limit temperature Te d vanishes for b b 0, as seen in Fig. 4, with (a) broken for b > b 0 and (b) for b < b 0. Nevertheless, Te d T e/ for b, so that the entanglement effect for large fields will be detected, whereas Te d T e for b 0, with T e Te d b for b 0. Note also that Te s lies very close to Td e for b 5v +, but becomes lower as b decreases, with Te s 0.39v + for b 0. Now d)] < 0.37 b, with C[ρ(Td e )] = for b = b 0 and C[ρ(Te)] d 0.37 and 0.3 for b b 0. This discontinuity arises from that of C(ρ) for T = 0 (where C(ρ) = at b = b c while C(ρ) and 0.7 for b b 0 respectively). For v z = 0, a stable mean field solution breaking phase flip symmetry becomes feasible only if b < b c = v x, with v x = v + +v > b 0 andt c givenbyeq.(8)withv M = v x. TheratioT c /v + isthenlargerthanincase. Forγ = 0.7, T c lies close to T e for b v +, with T c 0.85v + and C[ρ(T c )] for b 0. However, the most striking effect is that T c > T e for. b/v +.33, i.e., for b just above the reentry interval. In this region, ρ becomes separable at a low temperature, yet correlations remain strong to induce a symmetry-breaking mean field. On the other hand, for b > b c the ground state remains entangled but correlations are not strong enough to induce symmetry-breaking, as occurs in case.

6 6 C(ρ) Te/v+.5 v + /b= T/v + exact disorder entropic mean field lowest energy levels are close and entangled, implying that such solutions, normally regarded as signatures of the presence of strong correlations in the system, are not rigorous indicators of entanglement for T > 0. It is well known that in small systems, the sharp thermal mean field transitions are to be interpreted just as rough indicators of a smooth crossover between two regimes. The concept of entanglement allows, however, to formulate a crossover precisely. Finite systems regain in this sense a critical-like behavior for increasing T, becoming classically correlated (but not uncorrelated) for T T e, and with an entanglement undetectable through the eigenvalues of ρ and one of its reductions for T d e T < T e. Acknowledgments. NC and RR acknowledge support, respectively, from CONICET and CIC of Argentina. RR also acknowledges a grant from Fundación Antorchas. C[ρ(TL)] 0 v + /b FIG. 4: Same quantities as in Fig. for v + > 0 and finite anisotropy γ = v +/v = 0.7. Entanglement vanishing plus reentry occurs as T increases for 0.9 v +/b.4, as indicated by the solid lines of the central panel. In the top panel, dashed (solid) lines depict the concurrence in the interval where Eq. 5a (5b) is broken. IV. CONCLUSIONS We have examined the exact and the disorder separability conditions for general mixed states of two qubits interacting through a general XY Z Heisenberg Hamiltonian, which can be succinctly expressed in terms of total spin expectation values. The disorder conditions are exact in the absence of a magnetic field, but become weaker as the field increases and are unable to detect the reentry of entanglement for increasing temperatures in thermal states, an effect which may here arise when the ground state is less entangled than the first excited state. The von Neumann entropic criterion leads to still lower limit temperatures and is not exact even for zero field. Nonetheless, both the disorder and entropic criteria do predict the increase in the limit temperature for large fields occurring in anisotropic models. The critical temperature for the symmetry-breaking mean field solution is normally also lower than the exact limit temperature for entanglement in the examples considered and always vanishes for sufficiently large fields. However, it can also be higher, particularly when the Appendix A The concurrence of a mixed state ρ of two qubits is a measure of the entanglement of ρ, given by [] C(ρ) = Max[λ M TrR,0], (A) where λ M is the largest eigenvalue of R = [ρ / ρρ / ] / and ρ the spin-flipped density operator, given in the standard basis by ρ = (σ y σ y )ρ (σ y σ y ), with σ y the Pauli matrix. The entanglement of formation[0] is an increasing function of C(ρ) and can be obtained as E(ρ) = ν=±q ν log q ν, q ± = (± C (ρ)). Maximum entanglement corresponds to C(ρ) =, separability to C(ρ) = 0. For a pure state ρ = Φ Φ, C(ρ) = Φ Φ and E(ρ) becomes the von Neumann entropy of the subsystems [], S (ρ A ) = S (ρ B ). For the state (3), the eigenvalues of R are λ, = {[(p +p ) b (p p ) ] ± v (p p )}, and λ 0,3 = p 0,3. Hence, If λ M = λ or λ (λ 0 or λ 3 ), Eq. (A) becomes the difference between the left and right hand sides of Eq. (5a) [(5b)] when positive. The eigenvalues of the partial transpose of (3) are q, = [p 0 +p 3 ± v (p p )] and q 0,3 = {p + p ± [(p 3 p 0 ) + b (p p ) ] }, sothat theconditionsq j 0 j also lead to Eqs. (5). Only one of them, q m, is negative when ρ is entangled [3], with q m = Min[q,q ] (Min[q 0,q 3 ]) if λ M = λ or λ (λ 0 or λ 3 ). In the first case C(ρ) = q m but in the second case, C(ρ) q m unless b = 0 or p = p. [] E. Schrödinger, Naturwissenschaften 3, 807 (935); Proc. Cambridge Philos. Soc. 3, 555 (935). [] C.H. Bennett et al., Phys. Rev. Lett. 70, 895 (993).

7 7 [3] A.K. Ekert, Phys. Rev. Lett. 67, 66 (99); Nature358, 4 (99). [4] D.P. DiVincenzo, Science 70, 55 (995). [5] C.H. Bennett and D.P. DiVincenzo, Nature (London) 404, 47 (000). [6] M.A. Nielsen and I. Chuang, Quantum Computation and Quantum Information, Cambridge Univ. Press (000). [7] R.F. Werner, Phys. Rev. A 40, 477 (989). [8] A. Peres, Phys. Rev. Lett. 77, 43 (996). [9] M. Horodecki, P. Horodecki, and R. Horodecki, Phys. Lett. A 3, (996). [0] C.H. Bennett, D.P. DiVincenzo, J.A. Smolin, W.K. Wootters, Phys. Rev. A 54, 384 (996). [] S. Hill and W.K. Wootters, Phys. Rev. Lett. 78, 50 (997);W.K. Wootters, Phys. Rev. Lett. 80, 45 (998). [] K. Zyczkowski, P. Horodecki, A. Sanpera, and M. Lewenstein, Phys. Rev. A 58, 883 (998). [3] S.L. Braunstein et al., Phys. Rev. Lett. 83, 054 (999). [4] M.A. Nielsen, Ph.D. thesis, University of New Mexico, 998 (unpublished); e-print quant-ph/ [5] V.E. Korepin, N.M. Bogoliubov,and A.G. Izergin, Quantum Inverse Scattering Method and Correlation Functions, Cambridge Univ. Press, (993). [6] A. Imamoglu et al., Phys. Rev. Lett. 83, 404 (999). [7] R. Raussendorf and H.J. Briegel, Phys. Rev. Lett. 86, 588 (00). [8] M.C. Arnesen, S. Bose and V. Vedral, Phys. Rev. Lett. 87, 0790 (00). [9] X. Wang, Phys. Rev. A 66, (00). [0] X. Wang, Phys. Rev. A 64, 033 (00). [] G.L. Kamta and A.F. Starace, Phys. Rev. Lett. 88, (00). [] T.J. Osborne and M.A. Nielsen, Phys. Rev. A 66, 030 (00). [3] M.A. Nielsen and J. Kempe, Phys. Rev. Lett. 86, 584 (00). [4] T. Hiroshima, Phys. Rev. Lett. 9, (003). [5] R. Rossignoli and N. Canosa, Phys. Rev. A 66, (00). [6] R. Horodecki and M. Horodecki, Phys. Rev. A 54, 838 (996). [7] C. Tsallis, J. Stat. Phys. 5, 479 (988). [8] N. Canosa and R. Rossignoli, Phys. Rev. Lett. 88, 7040 (00). [9] P. Ring and P. Schuck, The Nuclear Many-Body Problem, Springer Verlag, New York, (980). [30] R. Rossignoli, Phys. Rev. C 45, 60 (99). [3] A. Sanpera, R. Tarrach, G. Vidal, Phys. Rev. A 58, 86 (998).