there is evidence that there is a correspondence@xcite between gravitational physics in anti  de sitter space  time and particular types of conformal field theory in one dimension fewer .
this duality is a form of `` holography ''  @xcite and a part of the correspondence operates by identifying the field theory as living on the boundary of anti  de sitter ( ads ) space  time .    to be more precise , ads@xmath0@xmath1 is the space  time of interest , and there is some @xmath2dimensional theory of gravity compactified on it .
the manifold @xmath3 can be an @xmath4sphere , @xmath5 .
the corresponding field theory is an @xmath6dimensional conformal field theory living on a space with the topology of the boundary of ads@xmath0 .
the isometries of the manifold @xmath3 appear as global symmetries of the field theory : r  symmetries if the theory is supersymmetric .
this particular form of duality between gravity and field theory is certainly intriguing . the large @xmath7 limit ( where  @xmath7 is the rank of the @xmath8 gauge group for the four dimensional yang
mills field theory , with appropriate generalizations for other dimensions ) of the field theory at strong t hooft coupling corresponds to classical supergravity . as pointed out in ref .
@xcite , following the observations in ref .
@xcite , the old program of semi  classical quantum gravity finds a new lease on life in this setting , as computations such as those performed with gravitational instantons ( at least in ads ) should have natural field theory interpretations .    in this paper , we study the thermal properties of einstein  maxwell ads ( emads ) charged black holes , and find behaviour consistent with field theory interpretations .
we do this for arbitrary dimensions ( greater than 3 see section  [ conclude ] for comments on @xmath9 ) and determine the thermal phase structure of the corresponding field theories .
the cases of ads@xmath10 , ads@xmath11 and ads@xmath12 are particularly interesting of course , as they correspond to the theories found on the world volumes of m2 , d3 , and m5branes , respectively .
the d3brane case is @xmath13 , @xmath14 supersymmetric @xmath8 yang  mills theory , while the others are exotic superconformal field theories  @xcite .
we remark on the field theory interpretation of our new results in the light of holography .
this paper is also of relevance beyond mere considerations of holography .
some of the black hole solutions and their properties ( thermodynamic or otherwise ) are presented here for the first time .
in particular , the lagrangian action calculations and subsequent determination of the phase structure are presented in their entirety here .    in section  [ sec : emads ] , we present an ansatz for obtaining the einstein  maxwell truncation of gauged ads supergravity with appropriate compactifications of @xmath15 supergravity on @xmath16 , and @xmath17 type iib supergravity on @xmath18 . in the planar or infinite  volume limit , the charged black holes in einstein  maxwell  anti
desitter correspond to the near horizon limits of rotating m2 and d3branes . in section  [ sec : holes ] , we display the solutions and note some of their properties .
the computation of the action of the solutions using a euclidean section is performed in section  [ sec : action ] , and their thermodynamic properties are uncovered in section  [ sec : thermo ] .    as the einstein  maxwell  anti
desitter truncation is naturally associated to rotating branes , ( at least in the case of em  ads@xmath10 and em
ads@xmath11 , see section  [ sec : emads ] ) it is very natural to suppose that there is an associated dual field theory , arising on the world - volume of some branes
. these would be the familiar conformal field theories the @xmath13 , @xmath14 yang  mills theory ( for coincident d3branes ) and the conformal field theory on the world  volume of coincident m2branes .
the case of em  ads@xmath12 ( _ i.e. , _ without additional scalars ) is not related to a rotating
brane truncation of the ads@xmath19 gauged supergravity ( because @xmath20 is even dimensional ) and so we can not declare that the dual field theory is the theory on the world
volume of a rotating m5brane .
however , we regard ads holography as a phenomenon which exists independently of string and m
 theory contexts@xcite .
hence , in other dimensions beyond @xmath13 and 5 , we expect that there is a dual theory .
in particular , for em
the dual field theory is probably a close cousin of the m5brane theory .
the dual field theories have their supersymmetry ( if they had any to start with ) broken due to coupling to a global background @xmath22 current ( as well as turning on a non  zero temperature ) .
the cft is in a thermal ensemble for which a certain @xmath22 charge density has also been `` turned on '' . in the ensemble ,
the expectation value of this charge breaks the global @xmath23 r  symmetry of the cft . on the ads side ,
the electromagnetic charge carried by the black holes is in the same @xmath22 of corresponding @xmath23 gauge group .
we find very interesting phase structures at intermediate temperatures ( in finite field theory volume ) , as a result of studying two complementary thermodynamic ensembles : we study thermodynamic ensembles with fixed background potential in which case the background is ads with a constant fixed potential and we also study a fixed localized charge ensemble , for which the background is an extremal black hole with that charge .    in all cases , at sufficiently high temperature the physics
is dominated by highly non  extreme black holes , and we therefore recover the `` unconfined '' behaviour characteristic of the associated field theories@xcite .
the finite horizon size of the black holes controls the behaviour of the expectation value of spatial wilson lines accordingly , yielding the area law behaviour , as follows from ref.@xcite .    at intermediate temperatures , in the fixed charge ensemble
, the presence of charge allows a new branch of black hole solutions to modify the qualitative phase structure in the low charge regime , resulting in a very interesting phase structure about which we will have more to say later in this section .    intriguingly , as there is an extremal but non  supersymmetric black hole with non
 zero entropy even at zero temperature , we must conclude something interesting about the field theory in the presence of the global background @xmath22 current : there must still be at @xmath24 a large number of states ( with the given charge ) available to the field theory in order to generate this entropy .
for the case where we hold the potential ( _ i.e. , _ _ not _ the charge ) fixed , we do not expect that this is the ground state , because the extremal black hole can decay into kaluza
klein particles , leaving ads .
this is because the extremal black hole is not supersymmetric .
this subtlety does not arise in the standard gibbons  hawking calculus of the thermodynamics of black holes which we use here because the calculations are not sensitive to the ability of the black holes to emit charged particles .    that the extremal black hole can decay by emitting charged kaluza
klein particles here follows from the fact that the charge descends from rotation in higher dimensions .
there are well  known classical processes for reducing the rotation of objects like black holes by scattering@xcite , and therefore in the context of quantum field theory , one has the analogous processes of emission in superradiant modes@xcite .
the same superradiant emission was considered in the context of charged black holes in ref.@xcite .
thus one should expect the extremal black hole in the emads truncation to decay _ via _ such supperradiant emission .
of course , the usual thermal hawking radiation may also tend to discharge nonextremal black holes@xcite . in the fixed potential ensemble , as the charge of the black hole is allowed to fluctuate while it is in contact with the thermal reservoir , superradiant and hawking emission processes can occur to reduce the charge of the black hole , allowing it to decay back to ads ( plus charge ) .
however , in the fixed charge thermodynamic ensemble ( with varying potential ) , the extremal black hole is expected to be the long  lived state at zero temperature .    translating the formula for the entropy to the field theory we find , for example , that the four dimensional yang  mills theory ( in the presence of the global background @xmath22 current )
has a zero  temperature entropy which goes like @xmath25 for large black holes , where @xmath26 measures the total charge in units of the minimal charge of kaluza
klein excitations ( _ i.e. , _
@xmath27 ) , and is proportional to the volume , @xmath28 , of the field theory . notice that the result for the four dimensional field theory is consistent with confinement at @xmath24 , as the result is independent of @xmath7 .
confinement also follows from the fact that at @xmath24 , the euclidean section of the solution has no bolt , and therefore temporal wilson lines will always be homotopic to zero , and therefore have zero expectation value . meanwhile
, spatial wilson lines can not interact with the horizon to produce an area law dependence , because at extremality the horizon recedes infinitely far away down a bertotti  robinson throat .
the phase structure which we obtain in each thermodynamic ensemble is summarized in figure  [ fig : phases ] .
0.5 cm 2.0 cm    the astute reader will recognize the figure on the right as the classic phase diagram of the liquid  gas system ! to translate , our @xmath29 is like the temperature @xmath30 of the fluid while @xmath31 is like the pressure @xmath32 .
the non  extreme black holes of type ( 1 ) ( `` small '' ) and ( 3 ) ( `` large '' ) ( see section  [ sec : action ] and [ sec : thermo ] for explanation ) are like the liquid phase and the gaseous phase , respectively . the critical line ( `` vapour pressure curve '' )
represents the place at which a first order phase transition between the liquid and gas occurs .
as is well known , there is a critical temperature at which the vapour pressure curve terminates , representing the fact that above a critical temperature , one can convert a liquid to a gas continuously .
this translates here into a critical charge above which the two types of black hole can be continuously converted into one another with no discontinuity in their size .    that this system ( first modeled by van der waals@xcite , with a crucial modification by maxwell ) appears in this ads black hole thermodynamics is fascinating , and
would not have been possible ( at least in this way ) without the presence of the extra branches of solutions which appear when there is negative cosmological constant .
we discuss this further in sections  [ sec : thermo ] and  [ sec : action ] .
further fascination may be found in the fact that the explicit shape of the free energy surface ( as a function of @xmath29 and @xmath30 ) is that of the classic `` swallowtail '' catastrophe , familiar from the study of bifurcations@xcite .
the control surface of the `` cusp '' catastrophe also appears , which ( of course ) follows from the well known fact that it is the shape of the van der waals equation of state , viewed as a surface in @xmath33 space .    that these shapes appear in this context
suggests that there is some exciting universality to be explored here : catastrophe theory is largely a classification of the possible distinct types of bifurcation shapes that can occur in a wide variety of complex systems .
this classification ( which , for the common `` elementary '' cases is of a
e type ) is equivalent to the ( perhaps more familiar ) classfication of singularities@xcite . it is of considerable interest to discover just what circumstances might give rise to the other members of the classification .
recalling that this all translates _ via _ holography into properties of a dual field theory , we would learn a great deal about universal phase structures which can occur there also .
physics near the horizon of supergravity branes can be described in terms of spontaneous compactification of supergravity . in the case of
non  dilatonic branes which will be the focus of the paper when the compactification takes place on a _ round _ @xmath4sphere the low energy degrees of freedom are described by an effective theory of einstein gravity with a negative cosmological constant coupled to @xmath23 gauge fields . the schwarzschild  anti  desitter black hole solutions of this theory have been used in the context of the ads / cft correspondence to infer thermal properties of the dual field theories  @xcite .
a natural extension of this program is to study ads black holes which are charged under a subgroup of the gauge symmetry of the gauged supergravity .
solutions of einstein  maxwell  anti
desitter in some dimensions are known , but in the context of string / m  theory , it is also interesting to determine how to make a truncation of the type  iib supergravity , or of 11 dimensional supergravity , which gives the emads effective action . in other words , we must make certain higher  dimensional choices which will result in the removal of the generic coupling of the @xmath34 term to scalars resulting from the kaluza  klein reduction .    amusingly , one simple way to introduce ( gauge ) charge on the black holes is by simply spinning or twisting the transverse ( angular ) sphere that becomes the compact space .
decoupling of the scalars is accomplished by choosing the spins in a maximally symmetric way . to be concrete , take ten dimensional iib supergravity , with the metric ansatz @xmath35 , \label{ansatz}\ ] ] where @xmath36 is a five - dimensional metric , @xmath37 , the variables @xmath38 are direction cosines on @xmath18 ( and therefore are not independent , @xmath39 we follow the notation of@xcite ) , and the @xmath40 are rotation angles on @xmath18 .
the ansatz for the rr 5-form field strength has `` electric '' components @xmath41 while the dual `` magnetic '' components are given by @xmath42 . in eqn .
( [ forms ] ) , @xmath43 is the volume form on the reduced five - dimensional space , and @xmath44 denotes hodge duality on this space .
the parameter @xmath45 measures the size of the @xmath18 and is given by the flux of the 5form field across the @xmath18 .
notice that a component @xmath46 in the time direction is interpreted as rotation of the @xmath18 in its three independent rotation planes , in equal amounts .
components in the spatial direction would instead be ` twists . ' for the sake of brevity , and since in this paper we will be mainly considering @xmath46 components , one can not define magnetic ( vector ) charges on the black holes . ] , we will refer collectively to them as `` rotations '' .    with this ansatz
( [ ansatz ] ) , the effective action in the five non  compact dimensions becomes @xmath47.\ ] ] this is precisely the einstein
maxwell  anti
desitter ( emads ) effective action we seek , with a chern  simons term .
the latter is indeed required by supersymmetry in @xmath48 five dimensional gauged supergravity  @xcite , whose bosonic sector is precisely described by the action above .
note that the gauge coupling is proportional to @xmath49 .
the ads@xmath50 gauged supergravity theory in five dimensions has an @xmath51 gauge symmetry , associated with the group of isometries of @xmath18 .
this is the r  symmetry group of the dual four dimensional @xmath14 superconformal yang
mills field theory living on the d3branes from which this near  horizon geometry arose .
the above spinning compactification corresponds to introducing rotation in the diagonal @xmath22 of the maximal abelian subgroup @xmath52 .
correspondingly , there must be a dual field theory to the emads truncation , which is simply the field theory on the world
volume of the rotating brane . from the field theory point of view
, the rotation corresponds to considering states or ensembles in which the dual global @xmath22 current ( a subgroup of the @xmath51 r  symmetry group ) has a nonvanishing expectation value . studying emads gravity and its solutions
will therefore be equivalent to studying properties of the conformal field theory in the presence of this background current vector fields , each associated with the three different independent rotations of @xmath18 , and two scalars that , roughly , measure the relative sizes of the distortions of the @xmath18 caused by rotation . for simplicity
, we will restrict ourselves to the case where all three rotations have the same magnitude , since it is only in this case that the scalars decouple and we find emads gravity .
this framework provides the cleanest interpretation in terms of the dual cft , since the number of spin parameters or charges precisely matches the number of field theory operators which are `` excited '' .
see refs .
@xcite for discussion of more general actions and solutions related to this . ] .    a similar construction can be obtained by starting from eleven dimensional supergravity .
the compactification in this case is equivalent to focusing on the near horizon region of m2branes . in this case , take @xmath53,\ ] ] leading to the ads@xmath10 theory with a maxwell term @xmath54.\ ] ] the reduction ansatz for the 4-form field strength is @xmath55 where @xmath56 is the volume form on the reduced four  dimensional space , and @xmath57 denotes hodge duality on this space .
chern  simons terms are absent in four dimensions .
appropriate inclusion of fermions leads to four dimensional @xmath48 gauged supergravity .
the more general @xmath58 theory with four independent gauge fields ( _ i.e. _ , four different rotation parameters ) 3 scalars and @xmath59 supersymmetry , as well as its black hole solutions , have been recently studied in ref .
@xcite .
we note here that there is no analogous construction for the ads@xmath19 gauged supergravity theory .
this is because  @xmath20 is even dimensional and therefore we can not have a symmetric split between @xmath22 rotations , as @xmath60 does not have an even torus for its cartan subalgebra .
this means that we can not relate the physics of the black hole solutions ( which we write later ) of the em
ads@xmath12 system to the physics of rotating m5branes of eleven dimensional supergravity .
nevertheless , as ads holography is a phenomenon which is expected to exist independently of string or m  theory realizations , we expect that the physics does have a holographic interpretation in terms of a field theory closely related to that which resides on m5brane world  volumes .
the black hole solutions of the above supergravity theories in @xmath61 were originally studied in the past in refs .
@xcite and  @xcite  more recent investigations appear in refs .  @xcite .
as we have seen in the previous section , such theories can be regarded as compactifications of the type iib and @xmath15 supergravities , where the gauge symmetry groups of the gauged supergravities are broken by a specific choice of rotation planes in the transverse compact spheres . given these considerations ,
it is natural to study the reissner
nordstrom  anti
desitter ( rnads ) black holes within the context of the ads / cft correspondence .    even if the bosonic einstein  maxwell  anti
 desitter ( emads ) theories admit supersymmetric extensions only in certain dimensions , it is easy and convenient to perform the analysis of their black hole solutions for arbitrary dimension . for space  time dimension @xmath62
, the action can be written as so as to absorb the prefactors in the action . ]
@xmath63 , \label{actionjackson}\ ] ] with @xmath64 the cosmological constant associated with the characteristic length scale @xmath45 .
then the metric on rnads may be written in static coordinates as @xmath65    where @xmath66 is the metric on the round unit @xmath67sphere , and the function @xmath68 takes the form @xmath69    here , @xmath4 is related to the adm mass of the hole , @xmath70 ( appropriately generalized to geometries asymptotic to ads@xcite ) , as @xmath71 where @xmath72 is the volume of the unit @xmath67sphere .
the parameter @xmath73 yields the charge @xmath74 of the ( pure electric ) gauge potential , which is @xmath75 where @xmath76 and @xmath77 is a constant ( to be fixed below ) .
if @xmath78 is the largest real positive root of @xmath68 , then in order for this rnads metric to describe a charged black hole with a non  singular horizon at @xmath79 , the latter must satisfy @xmath80 finally , we choose @xmath81 which then fixes @xmath82 .
the physical significance of the quantity @xmath77 , which plays an important role later , is that it is the electrostatic potential difference between the horizon and infinity .    if the inequality in eqn .
( [ extbound ] ) is saturated the horizon is degenerate and we get an extremal black hole .
this inequality imposes a bound on the black hole mass parameter of the form @xmath83 . in the cases where the theory admits a supersymmetric embedding one could naively expect to approach a supersymmetric state as we saturate this mass bound .
however , the bound that results from the supersymmetry algebra is instead  @xcite : @xmath84 , with the @xmath85 solution being a bps state , where the black hole can have magnetic charge @xmath86 , there is a magnetic ( or dyonic ) bps solution as well @xcite with @xmath87 , @xmath88 . ] .
now , it is easy to see that the mass of the extremal black hole , @xmath89 is , for finite @xmath45 , always strictly larger than @xmath90 and therefore the extremal solution is non  supersymmetric . on the other hand ,
for the supersymmetric solution one has @xmath91 which is strictly positive everywhere and therefore one finds a naked curvature singularity at @xmath92 .
in fact , all the solutions violating the bound  ( [ extbound ] ) are nakedly singular .    in the context of the ads
/ cft correspondence it is interesting to consider the limit where the boundary of ads@xmath0 is @xmath93 instead of @xmath94 as was the case above .
this can be regarded as an `` infinite volume limit '' , with particular relevance to the discussion of the dual field theory .
it should be noted that the existence of black hole solutions in this limit is possible only due to the presence of a negative cosmological constant .
in fact , black holes ( and other bolts ) in ads spaces with varied topologies ( even other than spherical and toroidal ) have been extensively studied in recent years  @xcite , including in the m
theory  @xcite . here
we will only focus on the planar ( toroidal ) solutions , which we will obtain by scaling the `` finite volume '' solutions above , as done in @xcite . to this effect , introduce a dimensionless parameter @xmath95 ( which we will shortly take to infinity ) and set @xmath96 while at the same time blowing up the @xmath97 as @xmath98 .
one finds , after taking @xmath99 , @xmath100 with @xmath101 for the supersymmetric solution , the scaling is as above except for the scaling of @xmath4 . to preserve supersymmetry , one must fix @xmath85 and so @xmath102 yielding @xmath103 notice that , compared with eqn .
( [ scaledu ] ) , the parameter @xmath4 is zero in this limit .    the resulting solution can be seen to be supersymmetric as well ( _ i.e. _ , the killing spinors remain finite in the limit @xmath104 , after appropriate rescaling ) , and nakedly singular . in this `` infinite volume ''
limit , the solutions asymptote to ads space with the horospheric slicing .
these planar solutions can be constructed with the appropriate decoupling limit@xcite of spinning d3 or m2branes , as mentioned previously .
we refer the reader to ref.@xcite for the details .
the study of the euclidean section ( @xmath105 ) of the solution , identifying the period , @xmath106 , of the imaginary time with inverse temperature , will define for us the grand canonical thermodynamic ensemble ( for fixed electric potential ) or the canonical ensemble ( for fixed electric charge ) .
we interpret this in terms of immersing the system into a thermal bath of quanta at temperature @xmath107 . for pure ads ,
the background consists of both charged and uncharged quanta free to fluctuate in the presence of fixed potential @xmath77 .
later , we consider the fixed @xmath29 ensemble . in that case
we localize all of the charge at a specific region and keep it fixed .
for such a background , as ads with a localized charge is _ not _ a solution of the emads equations , we use the extremal black hole solution as the background , and retain only neutral quanta in the thermal reservoir , in order to keep the charged fixed .
this makes sense , even though the extreme limit has zero temperature , since the euclidean section has no bolt and so can be assigned an arbitrary periodicity@xcite .
hence , the metrics and gauge fields can be matched in the asymptotic region .    with all of this in mind
we now turn to the action calculations .      with our conventions the full euclidean action
is given by analytically continuing eqn .
( [ actionjackson ] ) , where , as usual when the space is asymptotically ads the gibbons  hawking boundary term gives a vanishing contribution .
the boundary terms from the gauge field will vanish if we keep the potential @xmath46 fixed at infinity .
any possible chern  simons term will not contribute when we restrict ourselves to purely electric solutions . imposing the equations of motion we can eliminate the factors of @xmath108 in order to obtain the on  shell action @xmath109.\ ] ]    we obtain for the action ( subtracting the ads background while remembering to match the geometries of the background and black hole in the asymptotic region ) : @xmath110 here , @xmath106 denotes the period of the euclidean section of the black hole space  time . using the usual formula for the period , @xmath111 , a little algebra yields the explicit form @xmath112
this may be rewritten in terms of the potential as : @xmath113    note that the temperature is zero when the black hole is extremal .
this is because the horizon is degenerate there , and @xmath106 diverges , together with the fact that one can smoothly approach the extremal limit from non  zero temperature . from the form of the equation for @xmath106 , it is apparent that there are qualitatively two distinct types of behaviour , determined by whether @xmath77 is less than or greater than the critical value @xmath114 .
in particular , for @xmath115 , @xmath106 diverges ( @xmath30 vanishes ) at @xmath116 , while for @xmath117 , @xmath106 goes smoothly towards zero as @xmath118 .
it is instructive to plot the temperature as a function of horizon radius ( size of black hole ) for these two regimes ( see fig .
[ betaplot3 ] ) .    -0.5 cm    as can be seen from the figure , the regime of large potential ( _ i.e. , _
@xmath119 ) has a unique black hole radius associated with each temperature .
we will see later that this branch dominates the thermodynamics for all temperatures .
meanwhile the small potential regime has two branches of allowed black hole solutions , a branch with larger radii and one with smaller .
this is qualitatively similar to the familiar case of the uncharged schwarzschild black holes analyzed in ref .
@xcite , ( or the structure of the taub
bolts discovered in the thermodynamic studies of ref .
@xcite ) , which is the @xmath120 limit of the situation here . correspondingly , the smaller branch of holes is unstable , having negative specific heat .
they do not play any role in the physics .
( generally , the sign of the specific heat for a black hole of radius @xmath78 can be inferred from the local slope of the @xmath121 curve .
see also the discussion in section  [ cats ] . )
if we wish to consider a situation where instead of the potential at infinity , we fix the charge of the black hole , then the action ( [ actionone ] ) is not appropriate .
upon variation of the gauge field in the latter action , a boundary term results that vanishes only if we keep @xmath122 fixed .
that is , the on  shell action of the previous subsection is @xmath123 $ ] .
if , instead , we want to keep the charge fixed , then we must add a boundary term to @xmath124 @xcite , @xmath125 where @xmath126 is a radial unit vector pointing outwards ( notice that this boundary term is determined by the terms coming from the variation of the off  shell action ( [ actionjackson ] ) , and not ( [ master ] ) , which is on  gravity  shell .
this distinction is only relevant for @xmath127 ) .
then we get a thermodynamic function @xmath128 $ ] , in terms of the variables we wish to control .    to compute the action for the fixed charge ensemble , using as background the extremal black hole , we evaluate eqn .
( [ master ] ) for a black hole of mass @xmath129 ( and radius @xmath130 ) , and then subtract the contribution from the extremal background . remembering to match the geometries of the background and black hole in the asymptotic region , a straightforward calculation yields the final result @xmath131 .
\label{actiontwo}\ ] ]    the inverse temperature , @xmath106 , is given by eqn .
( [ betaform ] ) .
it is useful to plot the temperature as a function of horizon radius ( size of black hole ) for future use .
there are two basic scales in this expression for @xmath132 , set by @xmath73 and @xmath45 , and so we expect that there will be two distinct regimes which may display distinct phase structure : @xmath133 and @xmath134 . for comparison , we also show the case of @xmath135 ( see fig .
[ betaplot1 ] ) .    the critical charge @xmath136 is the value of @xmath73 at which the turning points of @xmath121 appear or disappear . with @xmath137 , the periodicity @xmath138 will have a point of inflection with respect to @xmath78 derivatives .
hence we can simultaneously satisfy @xmath139 with @xmath140 and @xmath141 .
a little algebra then yields @xmath142 therefore we have for @xmath143 , @xmath144 and for @xmath145 , @xmath146    -0.5 cm    in this case , the figures show that for small charge ( _ i.e. _ , below @xmath136 ) , there can be _ three _ branches of black hole solutions , to which we will refer later .
the middle branch is unstable ) , @xmath147 while the branch with the smallest radii is new , and will play an interesting role in the thermodynamics . for zero charge ,
we return to the familiar two branch situation of schwarzschild , while for large charge , we have a situation analogous to that seen for the large fixed potential .
this is the grand canonical ensemble , at fixed temperature and fixed potential .
the grand canonical ( gibbs ) potential is @xmath148 using the expression in eqn .
( [ actionone ] ) , we may compute the state variables of the system as follows : @xmath149 together , they indeed satisfy the first law : @xmath150 .    in order to study the phase structure and stability
, we must observe the free energy @xmath151 as a function of the temperature .
it is shown in figure  [ freeplot1 ] .    -.5
cm    the interpretation of this is as follows . at any non  zero temperature , for large potential ( @xmath152 ) the charged black hole is thermodynamically preferred , as its free energy ( relative to the background of ads with a fixed potential ) is strictly negative for all temperatures .
this behaviour differs sharply from the small potential ( @xmath117 ) situation , which is qualitatively the same as the uncharged case : in that situation , in finite volume , the free energy is positive for some range @xmath153 , and it is only above @xmath154 that the thermodynamics is dominated by schwarzschild black holes ( the larger , stable branch ) , after their free energy is negative .
( see the center graphs in figure  [ freeplot1 ] . )
so for high enough temperature in all cases the physics is dominated by non
extremal black holes . in this case , ( after converting gravitational to field theory quantities , @xmath155 and @xmath156 ; for @xmath145 , @xmath157 and @xmath158 ; for @xmath159 , @xmath160 and @xmath161 . ] ) the free energy and entropy behave at ultra  high temperature as @xmath162 where @xmath163 is the @xmath67dimensional spatial volume upon which the field theory lives .
this is the `` unconfined '' behaviour appropriate to the dual @xmath6dimensional field theory .
the function @xmath164 is @xmath165 when @xmath145 ; @xmath166 when @xmath143 ; @xmath167 when @xmath159 .
the resulting power of @xmath7 shows how the number of unconfined degrees of freedom of the theory goes with @xmath7 , by analogy with the case of @xmath145 where @xmath168 counts the dependence on the number of degrees of freedom on @xmath7 for an @xmath8 gauge theory .    at low temperatures , and for @xmath152
, we have something very new .
notice that as we go to @xmath24 , the free energy curve approaches a maximum value which is less than zero .
this implies that even at zero temperature the thermodynamic ensemble is dominated by a black hole . from the temperature curve  ( [ betaplot3 ] ) it is clear that it is the extremal black hole , with radius @xmath169 . for @xmath170 , at @xmath24
we recover ads space .
so this suggests that even at zero temperature the system prefers to be in a state with non
zero entropy ( given by the area of the black hole ) .
notice that this @xmath24 situation displays the `` confined '' behaviour characteristic of the ordinary conformally invariant zero  temperature phase , despite the presence of the black hole .
this follows from the fact that the temporal wilson lines will still have zero expectation value , as the fundamental strings which define them can not wind the horizon which has infinite period at zero temperature .
similarly , spatial wilson lines will not display the area law behaviour , because the fundamental string world
sheets can not be obstructed by the horizon , because at extremality , it is infinitely far away down a throat .
having pointed out this intriguing possible zero temperature behaviour , we expect that for the case of _ fixed potential _ considered here , this is not the complete story .
we must allow for the possibility that the extremal black hole might decay due to processes involving kaluza
klein particles charged under the @xmath22 .
( see the discussion near the end of section  i. ) this possibility can not be discounted because the extremal black hole is not supersymmetric , as pointed out before , and therefore not guaranteed to be stable by the supersymmetry algebra .
we expect that calculations which include the effects of charge emission will shift the free energy back to zero , representing the true , equilibrium situation .
alternatively if we consider the action ( [ actionjackson ] ) on its own merit outside of string or m  theory compactifications , it may be regarded as part of a theory without fundamental charged particles .
the resulting thermodynamic phase structure for the fixed potential ensemble is summarized in the left diagram of figure  [ fig : phases ] .
we have seen that we may consider a @xmath24 background containing an extremal black hole of charge @xmath29 .
let us now keep this charge fixed and allow the potential at infinity to vary .
this is the canonical ensemble , and the corresponding thermodynamic potential , the free energy , is @xmath171 the energy , entropy and electric potential are computed as @xmath172 in this case @xmath173 measures the energy above the ground state , which is the extremal black hole .
together , they satisfy the first law , which in this case should be written as @xmath174 .
the free energy as a function of temperature is shown below for the cases of small and large charge , respectively ( compare to the third graph in figure  ( [ freeplot1 ] ) for the uncharged case ) .    -0.5
cm    that there are three branches for the small charge case follows from the second graph in figure  [ betaplot1 ] , which is magnified and labeled in fig
.  [ freeplot3 ] , on the right .
from there , it is clear that for low temperature there can only be one solution ( `` branch 1 '' ) for the black hole radius .
at some temperature @xmath175 , the origin of two new branches ( `` branches 2 and 3 '' ) of solutions appears ( @xmath176 for the chosen parameters in the plot . ) .
above this temperature ( below @xmath177 ) , there are therefore three distinct branches of solution until at temperature @xmath178 , ( @xmath179 in the plot ) two of the branches ( 1 and 2 ) coalesce and disappear , leaving again only a single branch ( 3 ) , which persists for all higher temperatures .
returning to the free energy plot , the meaning is now clear .
starting to the extreme left of the plot , ( low temperature ) we see that there is a single branch of free energy , corresponding to the branch 1 solutions . at @xmath180
, branches 2 and 3 appear on the graph and separate from each other at higher temperatures . at @xmath181 , branches 1 and 2 coalesce and disappear , while branch 3 persists for all higher temperatures , continuing to the left .
so from zero temperature the negative free energy of branch  1 means that those non  extreme black holes dominate the thermodynamic ensemble . at temperature @xmath154 , ( @xmath182 in the plot )
the free energy of branch 3 is actually more negative than that of branch  1 , and so that branch of non  extremal black holes takes over the physics and continue to do so for all higher temperatures .
the situation at @xmath154 is a genuine finite temperature phase transition , of first order .
( notice from the first graph in figure  [ freeplot3 ] that the free energy is continuous , but its first derivative is discontinuous . )
this results from the jump ( along the dotted line in the final graph in figure  [ freeplot3 ] ) from branch 1 to branch 2 , from small to large @xmath78 black holes , as the temperature increases . as the entropy is proportional to @xmath183
, there is a jump in the entropy , or a release of `` latent heat '' .
as we approach the critical value , @xmath136 , of the charge representing the crossover into the large charge regime , the kink in the free energy and therefore the transition vanishes , as branches 1 and 3 merge ( and branch 2 disappears ) . the difference in horizon radii between the two branches , @xmath184 , may be thought of as an order parameter for the transition , as it vanishes above @xmath136 , where the transition goes away .    as noted before in the case of fixed potential ensemble , branches 2 and 3 are the exact analogues of the small and large schwarzschild black holes of hawking and page@xcite , or the small and large taub
bolts discovered in the thermodynamic studies of ref .  @xcite . in those papers , above a certain temperature @xmath180 , there were two allowed solutions at a given temperature , the smaller ( branch 2 ) being unstable and the larger ( branch 3 ) being stable , which persists to dominate the thermodynamics above some critical temperature @xmath154 .
the existence of a stable branch 1 , and its merger with branch @xmath165 to disappear at @xmath181 is a new feature when we add a small fixed charge to the story .
conversely , if we start from a situation where charge is present on the black hole but the cosmological constant vanishes , then we find branches 1 and 2 , and it is only when the negative cosmological constant is turned on that branch 3 appears .    for large charge , there is only a single branch allowed , ( see figure  [ freeplot3 ] , the cusps collide and disappear ) and the associated thermodynamic story is correspondingly simpler .
the free energy shows that the non  extreme charged black holes dominate from @xmath24 .    in all cases ( large or small @xmath29 ) ,
the ultra high temperature phases are dominated by a black hole and the free energy and entropy have the characteristic `` unconfined '' field theory behaviour shown in eqns .
( [ behave ] ) .
one might examine the approach to the critical point more closely .
in particular , consider the behaviour of the specific heat @xmath185 with @xmath141 , as the temperature approaches the critical value , one finds a singularity with @xmath186 .
this behavior may be contrasted with the @xmath187 singularity found in ref.@xcite .
the essential difference is of course that near the critical point we have a point of inflection with @xmath188 , while ref.@xcite considers a minimum with @xmath189 .    the evolution of the free energy of the system as a function of charge is particularly interesting as one goes from zero charge to large charge .
the single cusp of the uncharged ( schwarzschild ) system is joined by a second cusp which comes in from infinity , forming ( with the original one ) a section of the well known `` swallowtail '' shape , familiar as a bifurcation set or `` catastrophe '' in singularity or catastrophe theory .
the significance of this is discussed in the next section .
as we cross over into the large charge regime at some critical value of @xmath73 , the cusps merge and the free energy becomes a simple monotonic function . for completeness
, we include a series of plots showing this evolution .
( we do not put them on the same axes , as we did for the fixed potential case , for the sake of clarity . )
-0.5 cm 0.5 cm -0.5 cm    the resulting thermodynamic phase structure for the fixed charge ensemble is summarized in the diagram on the right in figure  [ fig : phases ] .
we can not refrain from further general comments upon the meaning and structure of the curves that we have uncovered in the previous sections .
although we plotted only the cases for the @xmath145 case , representing ads@xmath11 ( and hence four dimensional field theory ) , the same universal structures appear in the cases @xmath143 and 6 as well , giving the same pleasing phase structure for the fixed charge ensemble .    the phase structure that we uncovered for
the fixed charge ensemble should remind the reader of the classic van der waals
maxwell behaviour , modeling the liquid  gas system .
indeed , they are isomorphic .
the @xmath121 curve ( the middle graph of figure  [ betaplot1 ] ) should recall the graph of the @xmath190 van der waals equation of state , where @xmath32 ( the pressure ) is replaced here by @xmath106 and @xmath191 ( volume ) by @xmath78 .
the instability of branch 2 is then simply the familiar instability of the corresponding section of the van  der waals curve .
the jump from branch  1 to branch  3 which we deduced from the form of the free energy is the precise analogue of the maxwell construction . in the isomorphism between our parameters and those of the van der waals
maxwell system , our charge @xmath29 is equivalent to their temperature @xmath30 .
the instability of branch  2 in both languages makes intuitive sense : as one increases the pressure , the volume should _ decrease _ , and therefore the positively sloped branch is not stable .
a similar statement holds for the black holes after making the translation to the current situation : for black holes in equilibrium with the heat bath , an increase in the temperature results in an increase in the black hole radius and hence mass , for stable black holes . notice that this also follows from the first law , recalling that the entropy is a positive power of the radius .
so the positive slope branch of the @xmath121 curve is generally unstable .    in the language of catastrophe theory@xcite the study of jumps in some `` state variables '' as a result of smooth changes in `` control variables'' the physical solutions of the @xmath192 curve , viewed as a two dimensional surface in @xmath193 space ,
is the `` control surface '' of the `` cusp catastrophe '' .
the cusp shape is the union of points in the @xmath194 plane ( the control variables ) where the state variable ( the allowed value of @xmath78 ) jumps from branch  1 to branch  3 , as branch  2 is unstable . after applying the minimum free energy condition to determine the allowed branches ( the `` maxwell criterion '' )
, the cusp catastrophe appears in the @xmath195 plane , ( or equivalently the @xmath196 plane ) collapsed to the critical line ( see figure  [ fig : phases ] ) ( or `` vapour pressure curve '' ) along which the two types of black hole can coexist , and across which there is a phase transition .
the end of the line , at the critical value @xmath136 , where branch  2 disappears , is the point where the distinction between branch  1 and  3 goes away .
the order parameter , @xmath197 for this critical point is the radius difference of the branches @xmath198 . beyond the critical charge
there is no phase transition ( @xmath199 ) in going from branch  1 black holes to branch  3 by increasing the temperature .
this is of course the familiar statement that above a critical temperature , there is no phase transition in going from a gas to a liquid by increase of pressure .
intriguing is the fact that the two dimensional free energy surface @xmath200 forms the shape of the swallowtail catastrophe ( see figure  [ freeplots ] ) .
( note that for @xmath201 and @xmath202 the shape is the same . )
this naturally follows from the ability of the @xmath121 curve to produce three branches , and the resulting shape for the free energy curve is the union of three branches .    here
, the swallowtail does not have the usual interpretation as a bifurcation surface ( like the cusp does above ) but it is natural to wonder whether its appearance tells us that there is some universality at work here .
this is because the language of catastrophe theory is largely a classification of the possible distinct types of bifurcation shapes that can occur .
this classification ( which , for the common `` elementary '' cases is of a
e type ) is equivalent to the ( perhaps more familiar ) classfication of singularities@xcite .
a natural question is whether or not the inclusion of more control parameters will always result in a free energy curve of a shape ( and corresponding phase structure ) which falls into the classification . it would certainly be amusing to find yet another case of the a
e classification appearing in string and m  theory physics .
the study of the thermodynamics of black holes in einstein  maxwell  anti  desitter is highly relevant to the thermodynamics of certain superconformal field theories with a background global current switched on .
this follows from the logic of the ads / cft correspondence , and the fact that the emads system can arise as the near  horizon physics of rotating m2 and d3branes , and it should therefore be regarded as the effective theory of the strongly coupled field theory living on the rotating branes world  volume .
the phase structures of the charged black hole systems studied here , and summarized in figure  [ fig : phases ] , are markedly different from those of the uncharged systems studied before in this context  @xcite .
the addition of charge revealed a rich phase structure , with precise analogues to classic thermodynamic systems .
the physics is consistent with a dual field theory interpretation .    in all cases ,
the infinite volume limit can be found by taking the limits given in eqn .
( [ scaled ] ) .
this scaling may be applied to the expressions for the actions ( eqns .
( [ actionone ] ) and  ( [ actiontwo ] ) ) and the period ( eqns .
( [ betaform ] ) and ( [ betaformtwo ] ) ) . in all cases ,
the result is that there is only one branch of black hole solutions ( like the large charge and potential situations had in finite volume ) , and the free energy is negative definite , showing that the thermodynamics is dominated by black holes for all temperatures .
of course , this is what we should expect , from the field theory point of view .
as we commented before , the gauge field in the ads space naturally couples to a cft current @xmath203 , following the prescription of ref .  @xcite . from the asymptotic variation of the gauge field  ( [ pure ] ) or its corresponding field strength , one then has an expectation value @xmath204 .
thus one might think of the cft state as containing a plasma of ( globally ) charged quanta .
the precise nature of the cft state depends on the ensemble , which we were studying . for the case of the fixed potential ,
the dual statement is that a chemical potential conjugate to the global charge has been introduced leading to the expectation value .
the fixed charge calculations correspond to an ensemble of cft states with a fixed global charge .
thus the difference between the two calculations is analogous to that between the canonical ( fixed @xmath30 ) and microcanonical ( fixed @xmath173 ) ensembles .    in the context of d3 branes with @xmath145 , the @xmath51 gauge fields couple to the r  symmetry currents in the super
yang  mills theory .
this aspect of the duality has been used to great advantage to produce nontrivial consistency tests by comparing correlators protected by supersymmetry@xcite .
of course in the present case , with the truncation to emads theory , we are focussed on a particular diagonal @xmath22 generator of the @xmath51 symmetry .    in this context , we can translate the results of the supergravity calculations to quantitative statements about the strong coupling behavior of the super  yang  mills theory .
up to numerical factors , we have as usual@xcite : @xmath205 and @xmath206 , ( where @xmath207 is the type  iib string coupling ) as well as @xmath208 .
it remains fix how the black hole charge should be characterized in the cft .
the most natural approach is to measure the physical charge  ( [ thecharge ] ) in terms of the fundamental charge of the kaluza
klein excitations in the ads space , _ i.e. _ , with @xmath209 . in the field theory
then , @xmath210 ( where @xmath28 is the spatial volume of the field theory ) essentially counts the number of fundamentally charged quanta per unit volume in a given state . given this framework
, we can consider the field theory content of our results .
for example , one might wonder what the critical charge ( [ thecritics ] ) appearing in the fixed charge phase diagram corresponds to : @xmath211    in general , translating the entropy , mass or free energy to a field theory expression produces a complicated function of both the temperature @xmath30 and the charge @xmath26 .
one relatively simple case is the high temperature limit , where the charge essentially plays no role ( see eqn .
( [ behave ] ) ) .
another interesting case to consider is that of the extremal black holes for which @xmath24 . by demanding that @xmath212 and @xmath213 have a consistent solution
, one finds that the mass and charge parameters are related by the following expression : @xmath214 where @xmath215 and @xmath216 .
a simple case to consider is that of a large black hole with @xmath217 , for which @xmath218 .
further in this limit , one has that @xmath219 and so @xmath220 notice that implicitly here we are considering a regime where @xmath221 .
the lack of dependence of the entropy on @xmath7 is a signal of confined behaviour at zero temperature , despite the presence of the black hole
. it would certainly be interesting if this entropy result could be recovered by considering partitioning of the charge @xmath26 amongst the charged excitations of the cft .
we have left aside the case of compactification of six dimensional supergravity on @xmath222 to get ads@xmath223 . by setting the @xmath222 in rotation in its two independent rotation planes , in a symmetric fashion
, we get an electric potential in ads@xmath223 . doing so ,
notice that if we start from the solution describing a rotating six dimensional black string ( such as the one obtained from the d1d5 bound state ) , then , in the throat limit the rotation of the @xmath222 can be undone by a diffeomorphism@xcite .
in other words , the effective gauge field in three dimensions is pure gauge .
nevertheless , as shown in ref.@xcite , there do exist charged black hole solutions in emads theory in three dimensions .
these have an electric potential that diverges logarithmically at infinity , which prevents one from defining the ensemble at fixed potential . nevertheless ,
if the extremal black hole background is subtracted , then the fixed charge ensemble can be appropriately defined . for non  rotating black holes ( in ref.@xcite , the full kerr
newman solution is constructed ) there is only one branch , just like we have found for large fixed charge ( see figure  [ betaplot1 ] , left ) , with corresponding simple thermodynamic structure given by figure  [ freeplot1 ] ( left ) .
finally , it is also worth remarking that the close similarity that we have observed with familiar structures from equilibrium thermodynamics , and expectations from a dual field theory is further encouragement ( for those who need it ) that the quantum mechanics of black holes is not unlike that of other situations .
ac is supported by pembroke college , cambridge . re is supported by epsrc through grant gr / l38158 ( uk ) , and by grant upv 063.310eb187/98 ( spain ) .
support for cvj s research , and other support of this project , was provided by an nsf career grant , # phy9733173 ( uk ) .
rcm s research was supported by nserc ( canada ) and fonds fcar du qubec .
this paper is report # s damtp199929 , ehu
ft/9902 , dtp99/9 , uk/9902 and mcgill/9907 .
we would like to thank neil constable , mike crescimanno , susan gardner , harry lam , guy moore , malcolm perry , and alfred shapere for comments and useful conversations .
ac , re and rcm would like to thank the department of physics and astronomy of the university of kentucky for its hospitality , and we would all like to thank the staff of the william  t.  young library ( university of kentucky ) for the use of excellent meeting facilities and other services .
99    j. maldacena , adv .
theor . math
* 2 * ( 1998 ) 231 , hep - th/9711200 s. s. gubser , i. r. klebanov , and a. m. polyakov , phys .
b428 * ( 1998 ) 105 , hep - th/9802109 .
e. witten , adv .
* 2 * ( 1998 ) 253 , hep - th/9802150 .
r. penrose , nuovo cimento 1 ( 1969 ) 252 .
starobinsky , j.e.t.p .
37 ( 1973 ) 28 ; + w. unruh , phys .
d10 ( 1974 ) 3195 ; + n. deruelle and r. ruffini , phys .
lett . 52b ( 1974 )
437 ; _ ibid .
_ , 57b ( 1975 ) 248
. g. denardo and r. ruffini , phys .
lett . 45b ( 1973 ) 259 ; + t. damour and r. ruffini , phys .
35 ( 1975 ) 463 .
gibbons , commun .
44 ( 1975 ) 245 .
s.w . hawking ,
taylor - robinson , phys.rev .
* d55 * ( 1997 ) 7680 , hep - th/9702045 .
d.n . page , phys .
d14 ( 1976 ) 3260 .
see , for example r. m. wald , ann .
82 ( 1974 ) 548 ; + r.m .
`` general relativity '' _ , the university of chicago press , chicago ( 1984 ) ; + r.m .
`` quantum field theory in curved space  time and black hole thermodynamics '' _ , the university of chicago press , chicago ( 1994 ) ; + n.d .
birrell and p.c.w .
davies , _ `` quantum fields in curved space '' _ , cambridge university press , cambridge ( 1982 ) ; + p. r. brady , i. g. moss , r. c. myers , phys .
* 80 * ( 1998 ) 3432 , gr - qc/9801032 .
k. behrndt , a.h .
chamseddine and w.a .
sabra , phys .
lett . b442
( 1998 ) 97 , hep - th/9807187 ; + k. behrndt , m cvetic and w.a .
sabra , hep - th/9810227 .
russo , hep - th/9808117 ; + c. csaki , y. oz , j.g .
russo and j. terning , phys .
d59 ( 1999 ) 065008 , hep - th/9810186 ; + p. kraus , f. larsen and s.p .
trivedi , hep - th/9811120 ; + j.g .
russo and k. sfetsos , hep - th/9901056 .
gubser , hep - th/9810255 ; + r.g .
cai and k.s .
soh , hep - th/9812121 .
mann , class .
* 14 * , l109 ( 1997 ) gr - qc/9607071
. + d.r .
brill , j.  louko and p.  peldan , phys .
rev . * d56 * , 3600 ( 1997 ) gr - qc/9705012 .
+ l.  vanzo , phys .
* d56 * , 6475 ( 1997 ) gr - qc/9705004 .
r.  emparan , phys . lett .
* b432 * , 74 ( 1998 ) hep - th/9804031 .    s.  sethi and l.  susskind , phys
* b400 * , 265 ( 1997 ) hep - th/9702101 .
+ t.  banks and n.  seiberg , nucl .
* b497 * , 41 ( 1997 ) hep - th/9702187 .
+ e.  witten , hep - th/9507121 .
+ a.  strominger , phys . lett .
* b383 * , 44 ( 1996 ) hep - th/9512059 + for a review , see n.  seiberg , nucl .
. suppl . * 67 * , 158 ( 1997 ) hep - th/9705117 .
gibbons and r.e .
kallosh , phys .
* d51 * ( 1995 ) 2839 , hep - th/9407118 ; + s.w .  hawking , g.t .
horowitz and s.f .
ross , phys .
* d51 * , 4302 ( 1995 ) gr - qc/9409013 ; + s.w . hawking and gary t. horowitz , class .
* 13 * ( 1996 ) 1487 , gr - qc/950101 ; + c.  teitelboim , phys
* d51 * ( 1995 ) 4315 , hep - th/9410103 .


einstein , podolsky and rosen ( epr ) formulated their famous argument about the completness of quantum mechanics in the framework of continuous variable systems @xcite .
however , after that bohm gave a dichotomized version of it @xcite , the debate concerning nonlocality moved to systems described by discrete variables , leading bell to formulate his celebrated inequalities in a dichotomized fashion @xcite .
recently , the increasing importance of continuous variable systems leads many authors to explore the nonlocality issue in its original setting , where dichotomic observables to test bell s inequalities are not uniquely determined .
the attempts to translate bell s inequalities to continuous variable systems clarified the fact that crucial in a nonlocality test is the existence of a set of dichotomized bounded observables used to perform the test itself , from which the so called `` bell operator '' is derived .
the more debated question has dealt with the nonlocality of the normalized version of the original epr state , _
i.e. _ the twin - beam ( twb ) state of radiation produced by spontaneuos downconversion in a parametric amplifier @xcite .
nonlocality of the twb state was not clear for a long time . using the wigner function approach ,
bell argued that the original epr state , and as a consequence the twb too , does not exhibit nonlocality because its wigner function is positive , and therefore represents a local hidden variable description @xcite .
more recently , banaszek and wodkiewicz @xcite showed instead how to reveal nonlocality of the epr state through the measurement of displaced parity operator .
furthermore a subsequent work of chen et al .
@xcite showed that twb s violation of bell s inequalities may achieve the maximum value admitted by quantum mechanics upon a suitable choice of the measured observables .
indeed , the amount of violation crucially depends on the kind of bell operator adopted in the analysis , ranging from no violation to maximal violation for the same ( entangled ) quantum state .
systems which involves only two parties are the simplest setting where to study violation of local realism in quantum mechanics .
a more complex scenario arises if multipartite systems are considered .
studying the peculiar quantum features of these systems is worthwhile in view of their relevance in the development of quantum communication technology , _
e.g. _  to manipulate and distribute information in a quantum communication network @xcite .
although the study of multipartite nonlocality has originated without the use of inequalities @xcite , an approach to derive bell inequalities has been developed @xcite also for these systems and applied to characterize their entanglement properties @xcite .
being originally developed in the framework of discrete variables , these multiparty bell inequalities have found application also in the characterization of continuous variable systems @xcite .    the aim of this paper is to apply the various approaches hitherto developed to test nonlocality of two- and three - mode continuous variable systems .
we will consider tripartite gaussian states as well as non - gaussian bipartite states .
in the first case strong violation of bell inequalities is found , allowing the bell factor to reach values of @xmath0 , while in the second case enhancement of nonlocality is obtained in comparison with the twb case .
the paper is organized as follows . in sec .
[ s : rev ] we review the different approaches to test nonlocality in the framework of continuous variables and introduce notation that will be used throughout the paper .
the three - mode states we are interested in are introduced in sec .
[ s:3 m ] , and their violation of local realism is analyzed in [ ss : dp3 ] and [ ss : ps3 ] . in sec .
[ s:2 m ] , the tripartite states are considered as sources for conditional generation of non - gaussian bipartite states , whose nonlocal proprties are then studied in [ ss : dp2 ] , [ ss : ps2 ] and [ ss : h2 ] .
finally , the main results obtained are summarized in sec .
[ s : conclusions ] , which closes the paper with some concluding remarks .
in this section we will briefly recall the inequalities imposed by local realism in the cases of our interest .
let us start by focusing our attention on a bipartite system .
let @xmath1 and @xmath2 denote two possible outcomes of two possible measurements on the first subsystem and similarly @xmath3 and @xmath4 for the second subsystem .
the essential feature of this measurements is that they are local , dichotomic and bounded .
the bell s combination @xmath5 under the assumption of local realism gives rise to the well known bell - chsh inequality @xcite : @xmath6 where @xmath7 is the correlation function between the measurement results , _
i.e. _ , the expectation value of the products of the results of the experiments @xmath8 and @xmath9 .    in the case of a n - partite system ,
a nonlocality test is possible using the bell - klyshko inequalities @xcite which provides a generalization of inequality ( [ bi2 ] ) .
these inequalities are based on the following recursively defined linear combination : @xmath10f_{n-1}+\frac12 \left[m(\alpha_n)-m(\alpha_n')\right]f'_{n-1 } \label{fn}\;,\ ] ] where @xmath11 and @xmath12 refer to measurements on the n - party of the system , and @xmath13 denote the same expression as @xmath14 but with all the @xmath15 and @xmath16 inverted . in the case of a three - partite system ,
local realism assumption imposes the following inequality from combination ( [ fn ] ) : @xmath17 where again @xmath18 is the correlation function between the measurement results .
quantum mechanical systems can violate inequalities ( [ bi2 ] ) and([bi3 ] ) by a maximal amount given by , respectively , @xmath19 and @xmath20 ( see , _
e.g. _ , ref .
@xcite ) .
we now briefly review three different strategies to reveal quantum nonlocality in the framework of continuous variables systems . recall that in the case of a discrete bipartite system , for example a spin-@xmath21 two particle system , the local dichotomic bounded observable usually taken into account
is the spin of the particle in a fixed direction , say @xmath22 .
hence the correlation between two measurements performed over the two particles is @xmath23 , where the operator @xmath24 is decomposed on the pauli matrices base and @xmath25 are two unit vectors .
the bell operator is then given by the expression : @xmath26    consider now a n - partite continuous variable system identified by the creation operator @xmath27 and the annihilation operator @xmath28 ( @xmath29 ) with boson commutation relations associated .
following the original argument by epr it is quite natural attempting to reveal the nonlocality of this system trying to infer quadratures of one subsystem from those of the others .
from now on , we will refer to this procedure as a `` homodyne nonlocality test '' , as quadrature measurements of radiation field are performed through homodyne detection . here
we identify the quadrature @xmath30 according to the definition @xmath31 . as they are local but neither bounded nor dichotomic , quadrature observables are not immediately suitable to perform a nonlocality test based on bell s inequalities .
the procedure to make them bounded and dichotomic is quite arbitrary and consist in the assignment of two domains @xmath32 and @xmath33 to each observable @xcite .
when the result of a quadrature measurement falls in the domain @xmath34 the value @xmath35 is associated to it .
usually the choice @xmath36 is considered , though a choice suitable to the system under investigation may be preferable . considering a bipartite system we can introduce the following quantities : @xmath37 where @xmath38 is the joint probability distribution of the quadratures @xmath39 and @xmath40 .
we can now identify the homodyne correlation function @xmath41 as @xmath42 which can be straightforwardly used to construct the bell combination @xmath43 of eq.([bi2 ] ) and to perform the nonlocality test .
the main problem of pursuing such a nonlocality test is that it is not suitable in case of systems described by a positive wigner function , as the twb state of radiation defined as @xmath44 , where @xmath45 and @xmath46 is the squeezing parameter .
indeed , a positive wigner function can be interpreted as a hidden phase - space probability distribution , preventing violation of bell - chsh inequality unless the measured observables have an unbounded wigner representation , which is not the case of the dichotomized quadrature measurement described above .
considering in fact that @xmath38 can be determined as a marginal distribution from the wigner function one can write from eqs .
( [ ph ] ) and ( [ eh ] ) : @xmath47 where the integration is performed over the whole phase - space and without loss of generality we have considered @xmath48 .
eq.([ehwpos ] ) itself is indeed a local hidden variable description of the correlation function , hence obeying inequality ( [ bi2 ] ) .    in order to overcome this obstacle different strategies
have been considered by many authors , based essentially on parity measurements .
banaszek and wodkiewicz @xcite demonstrated the nonlocality of the twb considering as local observable on subsystem @xmath49 the parity operator on the state displaced by @xmath15 ( hence we will refer to this procedure as a `` displaced parity ( dp ) nonlocality test '' ) , which is dichotomic and bounded : @xmath50 in the above formula , @xmath51 , while @xmath52 and @xmath53 $ ] denote the number operator and the phase space displacement operator for the subsystem @xmath49 .
hence the correlation function reads : @xmath54 from which bell s combinations @xmath55 in eq.([bi2 ] ) and @xmath56 in eq.([bi3 ] ) can be easily reconstructed in the cases @xmath57 .
the reason why this procedure would be able to reveal nonlocality also in case of quantum states characterized by a positive wigner function is clear using the following relation : @xmath58 indeed , the analog of eq.([ehwpos ] ) is : @xmath59 being the dirac-@xmath60 distribution unbounded , ineqs .
( [ bi2 ] ) and ( [ bi3 ] ) are no more necessarily valid for @xmath55 and @xmath56 .
the maximal violation found with this procedure for a epr state is @xmath61 @xcite , still far from the maximum violation admitted by quantum mechanics .
another strategy , developed by chen et al .
@xcite , shares a similar behavior as the one described above , allowing to reveal nonlocality for quantum states with positive wigner function .
interestingly , this type of nonlocality test , which we will refer to as `` pseudospin ( ps ) nonlocality test '' , admits a maximum violation for the epr state . it can be seen as a generalization to continuous variable systems of the one introduced by gisin and peres for the case of discrete variable systems @xcite , hence , for the case of a pure bipartite system , it is equivalent to an entanglement test @xcite
. let us consider the following set of operators , known as _
pseudospins _ in view of their commutation relations , @xmath62 acting on the @xmath49-th subsystem @xmath63 where @xmath64 and @xmath65 is a unit vector associated to the angles @xmath66 and @xmath67 . in analogy to the spin-@xmath21 system and defining @xmath68 the correlation function
is simply given by : @xmath69 from which the bell combinations @xmath70 and @xmath71 are evaluated .
also different representations of the spin-@xmath21 algebra have been discussed in the recent literature @xcite . in particular in ref .  @xcite it has been pointed out that different representations lead to different expectation values of the bell operators .
hence , the violation of bell inequality for continuous variable systems turns out to depend , besides to orientational parameters , also to configurational ones . in the following sections we will also consider the pseudospin operators @xmath72 taken into account in ref .
@xcite , which have the following wigner representation : @xmath73 where @xmath74 , @xmath75 and @xmath76 stands for the `` principal value '' .
the correlation function obtained using operators @xmath77 will be indicated as @xmath78 .
in this section we will analyze the nonlocal properties of tripartite gaussian states .
in particular we will consider two classes of states , the first one proposed by van loock and braunstein @xcite , the second one proposed in @xcite .
the reason why we consider this two classes is that the first is a very natural and scalable way to produce multimode entanglement using only passive optical elements and single squeezers , while the second one is the simplest way to produce three mode entanglement using a single nonlinear optical device .
indeed , both states can be achieved experimentally @xcite . as concern the first class of states , it is generated with the aid of three single mode squeezed states combined in a `` tritter '' ( a three mode generalization of a beam - splitter ) .
the evolution is then ruled by a sequence of single and two mode quadratic hamiltonians . as a consequence ,
being generated from vacuum , the three - mode entangled state is gaussian and its wigner function is given by : @xmath79 \label{wvlb}\;,\ ] ] where @xmath80 , @xmath81 are the positions and momenta of the three modes and @xmath82 is the inverse of the covariance matrix , whose explicit expression reads : @xmath83 where @xmath84 , @xmath85 , @xmath86 and @xmath46 is the squeezing parameter ( with equal squeezing in all initial modes ) . the second class of tripartite entangled states is generated in a single non linear crystal through the following interaction hamiltonian : @xmath87 @xmath88 describes two interlinked bilinear interactions taking place among three modes of the radiation field coupled with the support of two parametric pumps .
it can be realized in @xmath89 media by a suitable configuration exposed in ref .
notice that the same dynamics can be implemented in different physical systems , including optomechanical couplers and bose - einstein condensates in the linear regime @xcite .
the effective coupling constants @xmath90 ,
@xmath91 , of the two parametric processes are proportional to the nonlinear susceptibilities and the pump intensities .
if we take the vacuum @xmath92 as the initial state , the evolved state @xmath93 belongs to the class of the coherent states of @xmath94 and it reads @xcite : @xmath95 where @xmath96 represent the average number of photons in the @xmath49-th mode and @xmath97 are phase factors .
the explicit expressions of @xmath98 are : @xmath99 ^ 2 \ ; , \nonumber \\ n_3 & = & \frac{|\gamma_1|^2}{\omega ^2 } \sin^2(\omega t ) \ ; , \label{ndit}\end{aligned}\ ] ] with @xmath100 and @xmath101 .
also for this second class , being the initial state gaussian and the hamiltonian quadratic , the evolved states will be gaussian .
the wigner function reads as follows @xcite : @xmath102 \label{wt}\;,\ ] ] where @xmath103 is the inverse of the covariance matrix , whose explicit expression is : @xmath104 where @xmath105 both classes of states are fully inseparable for any value of the coupling constants , namely can not be written as a factorized state for any grouping of the modes .
therefore they are good candidates to reveal true tripartite nonlocality .
being gaussian states , however , nonlocality can not be revealed by homodyne detection . in the following
we analyze the results for displaced parity and pseudospin tests .
let us start the study of tripartite systems nonlocality using the `` displaced parity test '' .
considering the correlation function @xmath106 given by eq .
( [ edp ] ) , the state ( [ wvlb ] ) was found in @xcite to give a maximal violation of @xmath107 in the limit of large squeezing and small displacement .
the study in @xcite however was performed for a particular choice of displacement parameters : @xmath108 and @xmath109 .
a numerical optimization of the displacement parameters lead us to identify a number of parameterizations that allow a significantly higher violation of bell s inequality . as an example , consider the one given by @xmath110 and @xmath111 from which follows that @xmath112 hence the remarkably high asymptotic value of @xmath113 is found for large @xmath46 and @xmath114 ( see fig .
[ f : b3dpvlbgen ] ) .
the importance of a suitable choice of the displacement parameters is apparent if this asymptotic value is compared to the violations obtained in the nonlocality study performed in ref .
@xcite . in that work in fact
generalizations to more than three modes of state ( [ wvlb ] ) were also considered , giving an increasing violation of bell inequality as the number of modes increases , but never founding a violation greater than @xmath115 .
determining the optimal choice of the displacement parameters for a given state is in general a challenging task . to our knowledge
indeed there exists no general prescription to find it out , and ultimately one must rely onto a numerical analysis ( see , _
nevertheless , a careful inspection of the symmetries of the state under consideration may be helpful . in order to clarify this observation
let us consider the explicit expression of the correlation function @xmath116 for state ( [ wvlb ] ) : @xmath117 \nonumber \\ -\frac{2}{3}e^{-2\,r}\left [ ( x_1+x_2+x_3)^2+(y_2-y_3)^2+(y_2-y_1)^2+(y_1-y_3)^2 \right ]    \bigg\}\end{aligned}\ ] ] hence , from eq .
( [ bi3 ] ) it follows that the bell combination @xmath56 is given by the sum of three positive and one negative term .
it is reasonable to expect that the maximal violation of nonlocality will be achieved for large @xmath46 .
we see from eq .
( [ vlbopt1 ] ) that , in this limit , all the correlation functions in @xmath56 vanish for nonzero displacements , unless the coefficients of @xmath118 are zero .
hence we impose the following system of equations , which allows for the three positive terms in @xmath56 not to vanish ( we consider @xmath119 and @xmath120 for @xmath121 ) : @xmath122 we see that the parameterization used to obtain eq .
( [ b3dpvlbgen ] ) is a solution of this system .
clearly , any other solution will give , in the limit of large @xmath46 , the same violation given by eq .
( [ b3dpvlbgen ] ) , namely @xmath123 .
in order to compare the violation of bell s inequality admitted by the state ( [ wvlb ] ) with the one that will be obtained below considering the state ( [ t ] ) , it is useful to rewrite eq.([b3dpvlbgen ] ) as a function of total mean photon number @xmath124 . given that @xmath125 and optimizing the displacement @xmath126 we obtain the result shown in fig [ f : b3dpn ] .
the asymptotic expression of the optimized displacement as a function of @xmath127 is @xmath128 , hence very small angles are required .    ) .
only values violating bell inequality ( [ bi3 ] ) are shown.,scaledwidth=50.0% ]    consider now the tripartite state generated by the hamiltonian ( [ inth ] ) .
the correlation function is now given by eq .
( [ edp ] ) through the wigner function ( [ wt ] ) .
the symmetry of the state suggests a maximum violation of bell inequality for @xmath129 ( recalling eq .
( [ ndit ] ) ) , while the fact that the separability of the state does nt depend on the phases @xmath130 and @xmath131 @xcite suggests that they are not crucial for the nonlocality test .
if we consider the same parametrization that led us to eq .
( [ b3dpvlbgen ] ) and fix @xmath132 , we find : @xmath133 from which follows an asymptotic violation of bell s inequalities of @xmath134 , for large @xmath127 and small @xmath126 .
a slightly better result is found if a parametrization , more suitable and numerically optimized for state ( [ t ] ) , is considered : @xmath135 and @xmath136 .
the bell combination @xmath56 for this choice of parameters is depicted in fig .
[ f : b3dpt ] .
we note that in this case a larger choice of angles allows the violation of bell inequality if compared with fig .
[ f : b3dpvlbgen ] .    )
( see text for details ) . only values violating bell inequality ( [ bi3 ] ) are shown.,scaledwidth=50.0% ]     for state ( [ wvlb ] ) , dotted line , and ( [ t ] ) , solid line .
the displacement parameters @xmath126 have been optimized to give the maximum value of @xmath137 for @xmath127 fixed.,scaledwidth=50.0% ]    as before , optimizing the displacement @xmath126 for each @xmath127 , it is possible to find the maximum violation of bell inequality as a function of @xmath127 .
we find that the asymptotic relation between the optimized displacement and the total photon number is now @xmath138 , confirming that not too small displacements are required .
the asymptotic violation of bell s inequality is now @xmath139 . to compare the results obtained for the two states ( [ wvlb ] ) and ( [ t ] ) we have plotted in fig . [ f : b3dpn ] the bell combination @xmath56 as a function of the mean total energy @xmath127 , while the displacement @xmath126 has been chosen in order to maximize @xmath56 at fixed energy .
notice that even if the two states have quite the same asymptotic violation , state ( [ wvlb ] ) reaches it for lower energies .
consider now the pseudospin nonlocality test .
let us calculate the expectation value of the correlation function ( [ eps ] ) for the state @xmath140 ( for simplicity we consider @xmath141 ) .
the only non vanishing contributes are given by : @xmath142 and by @xmath143 .
the correlation function then , according to eqs .
( [ ps ] ) and ( [ eps ] ) , reads as follows : @xmath144 hence , without loss of generality , we can fix for example @xmath145 and @xmath146 and look for the maximum violation of bell inequality ( [ bi3 ] ) constructed from eq .
( [ eps3ms ] ) .
notice that if the coefficients @xmath147 ( @xmath148 ) were equal to one then the maximum violation admitted , @xmath149 , should be reached .
considering eqs .
( [ c1c2c3 ] ) two limiting cases can be studied .
first , for large @xmath150 and small @xmath151 ( or vice - versa ) a numerical evaluation of the coefficients @xmath147 shows that @xmath152 , while the other two vanish .
hence , considering @xmath153 , the correlation function ( [ eps3ms ] ) reduces to that of a twb subjected to a pseudospin nonlocality test ( see eq .
( [ etwbps ] ) below ) , hence allowing an asymptotic violation of @xmath154 .
this result should be expected , since in this limiting case the state ( [ t ] ) reduces to a twb for modes @xmath155 and @xmath156 , while mode @xmath157 remains in the vacuum state and factors out .
consider now the case in which @xmath158 .
a numerical evaluation shows that the coefficients @xmath159 for large @xmath127 , hence also in this case the maximum violation can not be attained .
the asymptotic violation turns out to be @xmath160 .    as already mentioned in sec .
[ s : rev ] other representations for the pseudospin operators can be considered . using eqs .
( [ wpi ] ) and ( [ wt ] ) it is possible to calculate the correlation function @xmath161 . setting
again the azimuthal angles @xmath162 , the latter shows the same structure as @xmath163 where now the coefficient @xmath147 are replaced by : @xmath164 an appropriate choice of angles leads to a maximal violation of bell s inequality given by @xmath165 ( see fig .
[ f : b3ps ] ) , which is now reached for @xmath166 , value for which the coefficients @xmath167 are approximately near their maxima .
as already pointed out , we see that different representations of the pseudospin operators give rise to different expectation values for the bell operator .
applying now the same procedure to state ( [ wvlb ] ) we find the same structure for the correlation function @xmath168 , where the coefficients are now given by : @xmath169 after an optimization of the angles @xmath170 we obtain a maximal violation of @xmath171 ( see fig .
[ f : b3ps ] ) , for @xmath172 ( @xmath173 ) that maximizes the coefficients @xmath147 .     for states ( [ wt ] ) and ( [ wvlb ] ) , solid and dashed line respectively .
@xmath127 is the total number of photons as in fig .
[ f : b3dpn ] , scaledwidth=50.0% ]
in this section we consider the tripartite state ( [ t ] ) as a source of two - mode states .
in particular , we study the nonlocality of a two - mode non - gaussian state obtained by a conditional measurement performed on state ( [ t ] ) .
gaussian states are at the heart of quantum information processing with continuous variables .
the reason for this is that gaussian is the character of the vacuum state of quantum electrodynamics .
this observation , in combination with the fact that the quantum evolutions achievable with current technology are described by hamiltonian operators at most bilinear in the quantum fields , accounts for the fact that the states commonly produced in laboratories are gaussian .
in fact , bilinear evolutions preserve the gaussian character of the vacuum state .
in addition , it is worth noticing that the operation of tracing out a mode from a multipartite gaussian state preserves the gaussian character too , and the same observation is valid when the evolution of a state in a standard noisy channel is considered .
indeed , the only feasible way to `` degaussify '' a state is through a conditional measurement , or by statistically mixing it with another gaussian state .
the reason to study non - gaussian states is that when the gaussian character is lost , then immediately the wigner function of the state becomes negative , for pure states , hence stronger nonclassical properties should emerge .
actually , various authors have recently investigated the nonlocality properties of two - mode non - gaussian states .
in particular , a twin - beam state subjected to inconclusive photon subtraction ( ips state ) has been considered in refs .
@xcite and @xcite , while in ref .
@xcite it has been pointed out that if the entangled coherent states @xcite could be produced experimentally they would allow for the maximal violation ( _ i.e. _ , @xmath174 ) both in case of a dp test as well as a ps test .    the most natural way to obtain a non - gaussian state from a gaussian
one is by elimination of its vacuum component .
in fact , such a state is necessarily described by a negative wigner function ( in fact @xmath175 ) .
do to the structure of state ( [ t ] ) its vacuum component can be subtracted by a conditional measurement on mode @xmath157 , the same observation being valid for mode @xmath156 .
consider a photodetector able to distinguish only the presence or the lack of photons , _
i.e. _ , an on / off photodetector , and the state @xmath176 conditioned to the presence of at least one photon .
the probability operator measure ( povm ) is two - valued @xmath177 , @xmath178 , with the element associated to the `` no photons '' result given by @xmath179 where @xmath180 is the efficiency of the photodetector .
the probability of the outcome is given by @xmath181 = \frac{\eta n_3}{(1+\eta n_3 ) } \:,\end{aligned}\ ] ] while the conditional output state reads as follows @xmath182 \nonumber \\
\hspace{-50pt }    \eql \frac{1+\eta n_3}{(1+n_1)\eta n_3 } \sum_{p=1}^\infty    \left(\frac{n_3}{1+n_1}\right)^p \frac{1-(1-\eta)^p}{p!}(a^\dagger)^p    \sum_{n , n ' } \left(\frac{n_2}{1+n_1}\right)^{n+n ' } |n\,n\rangle    \langle n'\,n'| a^p \nonumber \\
\hspace{-50pt}\eql    \frac{1+\eta n_3}{(1+n_1+n_2)\eta n_3 } \sum_{p=1}^\infty    \left(\frac{n_3}{1+n_1}\right)^p \frac{1-(1-\eta)^p}{p ! }     ( a^\dagger)^p |x\rangle\langle x| a^p     \:,\end{aligned}\ ] ] where we have identified the twb with @xmath183 . to calculate the wigner function associated with state @xmath184 ,
consider that the characteristic function of the povm @xmath185 reads as follows : @xmath186(\mu)=\pi\delta^2(\mu)-\frac1\eta\exp\left[-|\mu|^2 \frac{2-\eta}{2\eta}\right]\;,\ ] ] hence the characteristic function of @xmath176 is given by : @xmath187(\lambda_1,\lambda_2)=\frac{1}{p_1 } \left\ { \chi[|t\rangle\langle t|](\lambda_1,\lambda_2,0 ) -\frac1\eta\int \frac{d^2\mu}{\pi}\chi[|t\rangle\langle t|](\lambda_1,\lambda_2,\mu ) \exp\left[-|\mu|^2\frac{2-\eta}{2\eta}\right ] \right\ } \;.\ ] ] after some algebra the wigner function associated with state @xmath184 can now be calculated .
it reads as follows : @xmath188 \nonumber \\
\lo- \frac{1}{\eta } \left(\frac{2}{\pi}\right)^2\frac{2}{\sqrt{\det d } } \exp\left [ -({\bf x},{\bf y } ) \left({\bf d}^{-1}\right ) ' \left ( \begin{array}{c } { \bf x } \\   { \bf y } \end{array } \right)\right ] \bigg\rbrace \label{w1}\;,\end{aligned}\ ] ] where , from now on , @xmath189 , @xmath190 , and @xmath191 .
in order to simplify the notation we have indicated with @xmath192 the @xmath193 matrix obtained from the @xmath194 matrix @xmath195 deleting the elements corresponding to the third mode ( 3-th row / column and 6-th row / column ) , due to the trace over the 3-th mode . of course , the easiest way to obtain a bipartite state from state ( [ t ] ) is to discard a mode , say the third one , by tracing over it .
the state @xmath196 then obtained is simply given by the following wigner function : @xmath197 \label{wtraced}\;.\ ] ] being the wigner function @xmath198 gaussian , we expect that this state will exhibit weaker nonlocality with respect to state ( [ rho1 ] ) . in the rest of the section the nonlocal properties of the usual twb state and of the states ( [ rho1 ] ) and ( [ wtraced ] ) will be compared .
notice that state ( [ wvlb ] ) can be considered as an extension to three modes of the twb .
all the three nonlocality tests introduced in sec .  [ s : rev ] will be taken into account .
we first study the violation of inequality ( [ bi2 ] ) in the case of a `` displaced parity test '' . as already mentioned , in ref .
@xcite banaszek and wodkiewicz found for the first time that the twb state exhibit a violation of local realism .
they obtained the following asymptotic violation for infinite energy : @xmath199 . generalizing their procedure this result
can be improved , yielding to a maximum asymptotic violation of @xmath200 @xcite .
we have considered the following parametrization to obtain the maximum violation for a twb : @xmath201 and @xmath202 .
the asymptotic relation between the squeezing parameter and the displacement angles is @xmath203 . using the same parametrization and considering the bell combination @xmath55 for the state @xmath196 , it turns out that the same asymptotic value of the twb is reached for large @xmath150 and small @xmath151 . in fact , as already noticed , when this limit is considered the original tripartite state ( [ t ] ) reduces to a factorized state composed by a twb and the vacuum state .
consider now the conditional state @xmath184 and again the case of large @xmath150 and small @xmath151 , say @xmath204 .
as in the tripartite case the phase coefficients @xmath130 and @xmath131 play no rule in the characterization of nonlocality . a stronger violation of bell inequality is then found and it is depicted in fig .
( [ f : b2dptwba ] ) , where the parametrization @xmath205 and @xmath206 has been adopted .
( see text for details ) .
only values violating inequality ( [ bi2 ] ) are shown.,scaledwidth=50.0% ]    indeed the asymptotic violation is higher then the previous , namely @xmath207 .
it can be found , for large @xmath150 , when @xmath208 .
a comparison with the violation of nonlocality attained with a ips state shows an identical asymptotic behavior @xcite . nevertheless , the scheme proposed here offers the advantage that the production rate of state @xmath176 [ _ i.e. _ , the conditional probability ( [ p1 ] ) ] is much greater then the production rate of ips state [ see ref .
@xcite , eq.(14 ) ] .
this is due to the fact that only a single on / off detection is required to produce @xmath176 , rather than the coincidence of two on / off detections for the case of ips state .
this could be useful from a practical viewpoint .
let us now focus on the `` pseudospin nonlocality test '' .
considering a twb state , it is known that the correlation function ( [ eps ] ) has the following expression ( setting to zero the azimuthal angles ) @xcite : @xmath209 where , denoting with @xmath127 the total photon number , @xmath210 it turns out that the violation of bell inequality in this contest increases monotonically to the maximum value of @xmath211 as the function @xmath212 goes to unity .
a straightforward calculation shows that an expression identical in form to eq .
( [ etwbps ] ) can be found both in case of state @xmath176 and @xmath196 , where the following functions @xmath213 and @xmath214 can be identified : @xmath215 in order to compare the violations in the three different cases , let us fix as in the previous subsection a small value for @xmath151 . a plot of the functions @xmath212 , @xmath213 and @xmath214 versus the total number of photons of the twb for the former and of the initial three - partite state for the latter two
is given in fig .
( [ f : b2ps ] )
. it can be seen that state @xmath176 achives large violations for smaller energies with respect to the other two states .
finally , a comparison with the violation attained with the ips state may be found in ref .
@xcite .
( solid line ) , @xmath216 ( dotted line ) and @xmath214 ( dashed line ) defined in the text ( the summation has been numerically performed for @xmath217 and @xmath218).,scaledwidth=50.0% ]      the negativity of the wigner function ( [ w1 ] ) may suggest to perform a nonlocality test based upon a homodyne detection scheme . while the positivity of a wigner function avoid to violate bell inequality ( [ bi2 ] ) with such a test , its negativity is yet not sufficient in general to ensure a violation .
quantum states with negative wigner function that does nt violate local realism with a homodyne test are given for example in refs .
@xcite . considering state ( [ rho1 ] ) it is necessary to calculate the correlation function ( [ ehwpos ] ) . substituting the wigner function ( [ w1 ] ) into eq .
( [ ehwpos ] ) and performing the integral we obtain the following result : @xmath219 \nonumber \\ \fl -\frac{1}{\eta } \left(\frac{2}{\pi}\right)^2\frac{2}{\sqrt{\det d } } \left [ \frac{2\,\pi \,\left ( -1 + n_3\,\left ( -2 + \eta \right )   \right)}{1 +       n_3\,\eta } \ ,      \arctan \left(\frac{2\,\cos\psi }        { { \sqrt{\frac{\left ( 1 + 2\,n_1 - n_3\,\eta \right ) \ ,                \left ( 1 + 2\,n_2 + n_3\,\eta \right ) } { \left ( 1 +                   n_1 \right ) \,n_2 } - 4\,\cos^2\psi}}}\right ) \right ] \bigg\ } \label{ehrho1}\;,\end{aligned}\ ] ] where @xmath220 . a plot of the correlation function ( [ ehrho1 ] ) is depicted in fig .
( [ f : e2h ] ) for unitary efficiency @xmath180 , together with the classical correlation function of two spin-@xmath21 particles @xcite ( see caption for details ) .
unfortunately , the comparison shows clearly that the correlation given by eq .
( [ ehrho1 ] ) is always lower then the classical one , hence despite the negativity of wigner function ( [ w1 ] ) no violation of bell inequality is achievable with this scheme .
particles classically correlated ( solid line ) and from eq .
( [ ehrho1 ] ) : @xmath221 ( dotted line),@xmath222 ( dot - dashed line),@xmath223 ( dashed line ) . in all cases we have fixed @xmath224 and @xmath225.,scaledwidth=50.0% ]    [ f : e2h ]
a detailed analysis of the nonlocality properties of multipartite continuous variables systems obtained by parametric optical systems has been presented , using the more recent approaches developed to this aim .
we have considered in particular two classes of tripartite gaussian states that seems promising for quantum communication purposes in order to implement multipartite quantum protocols .
the results show that for these states a nonlocality test based on displaced parity measurements is more suitable to reveal violation of local realism than one based on pseudospin operators .
this results are just the opposite of what have been obtained for the bipartite case .
notice , however , that a systematic approach to pseudospin operators for continuous variables havent been developed yet , hence we have only used the two inequivalent configurational parameterizations more suitable for calculations . for displaced parity test we obtained a remarkably high asymptotic value for the bell parameter , @xmath226 . in this case
the choice of a proper parametrization , suitable for the state under investigation , have revealed crucial .
we have also explored the possibility to enhance nonlocality in bipartite systems considering states endowed with a nonpositive wigner function .
we investigated a method to conditionally produce such a state from the tripartite systems considered above .
as expected , the bell parameter reaches a value higher than for a twb , namely @xmath227 .
also in the case of a pseudospin test an enhancement of nonlocality has been demonstrated , while a violation of local realism using a dichotomic quadrature measurement can not be achieved .
the authors are grateful to s.  olivares and m.  s.  kim for fruitful discussions .
99 a. einstein , b. podolsky and n. rosen , phys . rev . * 47 * , 777 ( 1935 ) .
d.  bohm , _ quantum theory _
( dover , nj , 1989 ) .
j.s . bell , physics ( long island city , n.y . ) * 1 * , 195 ( 1964 ) .
d.f . walls and g.j .
milburn , _ quantum optics _ ( berlin , springer , 1995 ) .
j.s . bell , _ speakable and unspeakable in quantum mechanics _ ( cambridge university press , cambridge , 1987 )
. k. banaszek and k. wodkiewicz , phys .
a * 58 * , 4345 ( 1998 ) .
chen , j.w .
pan , g. hou , and y.d .
zhang , phys .
. lett . * 88 * , 040406 ( 2002 ) .
p.  van loock and s.l .
braunstein , phys .
lett . * 84 * , 3482 ( 2000 ) .
h.  yonezawa , t.  aoki , and a.  furusawa , nature * 431 * , 430 ( 2004 ) .
greenberger , m.a .
horne , a.  shimony , and a.  zeilinger , am .  j.  phys . * 58 * , 1131 ( 1990 ) .
mermin , phys .
lett . * 65 * , 1838 ( 1990 ) ; d. n. klyshko , phys .
a * 172 * , 399 ( 1993 ) .
n. gisin and h. bechmann - pasquinucci , phys .
a * 246 * , 1 ( 1998 ) .
p. van loock and s.l .
braunstein , phys .
a * 63 * , 0222106 ( 2001 ) . z.b .  chen and y.d .
zhang , phys .
a * 65 * , 044102 ( 2002 ) .
clauser , m.a .
horne , a. shimony and r. a. holt , phys .
* 23 * , 880 ( 1969 ) .
a. gilchrist , p .
deuar , m.d .
reid , phys .
lett . * 80 * , 3169 ( 1998 ) .
h. jeong , w. son , m.s .
kim , d. ahn , and  .
brukner , phys .
a * 67 * , 012106 ( 2003 ) .
n. gisin and a. peres , phys .
a * 162 * , 15 ( 1992 ) .
l. mita , r. filip , and j. fiuraek , phys .
a * 65 * , 062315 ( 2002 ) ; l. praxmeyer , b.g .
englert , and k. wodkiewicz , eur .
j.  phys .
d * 32 * , 227 ( 2005 ) .
g. gour , f.c .
khanna , a. mann , m. revzen , phys .
a * 324 * , 415 ( 2004 ) ; m. revzen , p.a .
mello , a. mann , l.m .
johansen , phys .
a * 71 * , 022103 ( 2005 ) .
a. ferraro , m.g.a .
paris , m. bondani , a. allevi , e. puddu , and a. andreoni , j. opt .
b * 21 * , 1241 ( 2004 ) .
t. aoki , n. takei , h. yonezawa , k. wakui , t. hiraoka , a. furusawa , and p. van loock , phys .
rev . lett . * 91 * , 080404 ( 2003 ) .
m. bondani , a. allevi , e. puddu , a. andreoni , a. ferraro , and m.g.a .
paris , opt .
* 29 * , 180 ( 2004 ) .
n. piovella , m. cola , r. bonifacio , phys . rev . a * 67
* , 013817 ( 2003 ) .
s. pirandola , s. mancini , d. vitali , and p. tombesi , phys . rev .
a * 68 * , 062317 ( 2003 ) .
rodionov , a.s .
chirkin , pisma zh .
. fiz . * 79 * , 311 ( 2004 ) [ jetp lett . * 79 * , 253 ( 2004 ) ] .
puri , phys .
rev . a * 50 * , 5309 ( 1994 ) .
m. m. cola , m. g. a. paris , and n. piovella , phys .
a * 70 * , 043809 ( 2004 ) .
k.  banaszek , a.  dragan , k.  wdkiewicz , and c.  radzewicz , phys .
a * 66 * , 043803 ( 2002 ) .
h.  nha , and h.  j.  carmichael , phys .
lett . * 93 * , 020401 ( 2004 ) ; r. garca - patrn snchez _
et al . _ ,
* 93 * , 130409 ( 2004 ) .
s.  olivares , and m.  g.  a.  paris , phys . rev .
a * 70 * , 032112 ( 2004 ) .
b.  c.  sanders , phys .
a * 45 * , 6811 ( 1992 ) .
s.  olivares , m.  g.  a.  paris , and r.  bonifacio , phys .
a * 67 * , 032314 ( 2003 ) .
a.  ferraro , s.  olivares and m.  g.  a.  paris , _ gaussian states in c.v .
quantum information _
( bibliopolis , napoli , 2005 ) . w.j .
munro , phys .
rev a * 59 * , 4197 ( 1999 ) ; s. mancini , v.i .
manko , e.v .
shchukin , and p. tombesi , j. opt .
b : quantum s o * 5 * , s333 ( 2003 ) .
a. peres , _ quantum theory : concepts and methods _
( kluwer , dordrecht , 1993 ) .