in a previous paper @xcite , we presented a method for reducing a general quadratic theory of gravitation to a canonical second - order form .
the quadratic action provides an example of a higher - derivative theory , where the gravitational equations of motion are higher than second order .
this means that the theory contains more degrees of freedom than just the simple massless graviton . by rewriting the action
, we showed that quadratic gravitation is classically equivalent to einstein s gravity coupled to a massive real scalar field and a massive symmetric tensor field describing a spin - two field , with a specific lagrangian .
we found that , even for simple quadratic gravity , the reduced action gave highly non - trivial potential energy and kinetic coupling terms .
this suggests that higher - derivative theories have a great deal of structure , not immediately apparent from a simple linear analysis where one expands around flat space .
one particular feature , which is the main thesis of this paper , is that the transformed theory may have a complicated vacuum structure . here , by a vacuum solution , we mean a stable solution of the second - order theory .
in this paper , we will only consider vacua in which the auxiliary fields are covariantly constant .
it is important to note that only once the vacuum of the transformed theory has been identified can the nature of the elementary field excitations be discussed .
thus , for instance , the exact mass and couplings of the excitations , as well as which excitations are potentially ghost - like , may be very different around different vacuum states .
in general , they will bear little relation to the excitations of the linearized analysis .
the vacuum structure of the quadratic theories remained comparatively simple .
the only stable vacuum was flat space with zero vacuum expectation value for both of the auxiliary fields .
the elementary excitations were a massive scalar field and a massive ghost - like spin - two field .
we would like to investigate how other higher - derivative theories can introduce a more interesting vacuum structure .
a natural generalization is to consider actions which are not quadratic but general functions of the curvature tensors .
since these tensors involve only second derivatives of the metric , the corresponding field equations can be at most fourth - order ; that is , the same order as the quadratic actions
. we can start by asking how many new degrees of freedom we expect in such theories . to give a rough count
, we turn to the cauchy problem . first recall that even higher - derivative gravity theories remain diffeomorphism invariant .
thus we are always able to transform away four components of the symmetric metric , leaving six free components .
since the field equations are fourth - order , if the cauchy problem can be solved we expect to be required to give four initial conditions for each independent component of the metric , namely the component itself and the first three time derivatives . this would imply that we have at most twelve degrees of freedom .
rewriting the theory in a second - order form , since the auxiliary fields are set equal to terms involving second derivatives of the metric , we expect that giving the initial conditions of the auxiliary fields is equivalent to fixing second and third derivatives of the metric , a total of twelve conditions .
thus we find that the auxiliary fields should carry six degrees of freedom , just as in the general quadratic theory .
this leaves a possible six further degrees of freedom in the metric .
however , since the metric in the second - order theory obeys einstein gravity , the number of degrees of freedom is reduced to the two helicity states of the usual massless graviton . by this rough argument
, we expect an action in the form of any general function of the curvature tensors to describe , at most , the propagation of a massless graviton plus an additional six degrees of freedom ; that is , the same as for the special case of the general quadratic theory discussed in our previous paper @xcite .
be this as it may , the structure of such theories is potentially much richer than in the quadratic case . as we will see , the vacuum structure is now non - trivial .
a second generalization is to consider gravity actions which have higher derivatives acting on the curvature tensors .
since these theories involve more than two derivatives on the metric , the corresponding field equations will be generically higher than fourth - order .
consequently , we now expect new degrees of freedom in addition to those of the general fourth - order theory .
again , we will show that the vacuum structure of such theories is , in general , non - trivial .
some discussion of reducing both types of generalized theories to a second - order form already exists in the literature . that the equations of motion following from actions given by a general function @xmath1 are equivalent to the equations of motion of a scalar field dilatonically coupled to gravity , was first shown by teyssandier and tourrenc  @xcite .
the equivalence at the level of the action was given by magnano _
et al . _
@xcite , who also rewrote the reduce theory in canonical form .
this latter group and jakubiec and kijowski  @xcite also rewrote actions given by a general function @xmath2 in second - order form , though without the canonical separation of the new degrees of freedom .
actions including derivatives of the scalar curvature were considered in the context of inflation by gottlber _
et al . _
@xcite , who again showed , in some special cases , the equivalence of the field equations to those describing scalar fields coupled to gravity . still at the level of the field equations , this work was extended to a general class of actions by schmidt  @xcite and wands  @xcite .
possible vacua of the @xmath1 theory are discussed , in terms of the original higher - derivative field equations , by barrow and ottewill  @xcite , and later in terms of the second - order field equations by barrow and cotsakis  @xcite .
this paper is organized as follows . in the next three section
we will consider actions which are general functions first of the curvature scalar only , then of the ricci tensor , and finally of the full riemann tensor .
rewriting these theories in a canonical second - order form , we will find a rich vacuum structure .
this will require a careful discussion of the different `` branches '' of the theory , a concept we define below
. we will show that the vacuum states are in a one - to - one correspondence with the stable , constant curvature , desitter or anti - desitter solutions of the higher - derivative theory .
furthermore , we will show that if the vacuum of the transformed theory has zero cosmological constant , then the corresponding spacetime for the original higher - derivative theory is also flat .
we will then discuss the elementary excitations around a given vacuum and argue that they remain a graviton , a scalar field and a ghost - like spin - two field , with masses and couplings depending on the particular vacuum and the functional form of the original action . in section [ fboxr ]
we will consider a class of gravity actions that are general functions of the curvature scalar @xmath3 and the derivatives @xmath0 , with @xmath4 a positive integer .
we show how to rewrite such theories in second - order form and show that they have , in general , a non - trivial vacuum structure . in particular , we present an example with two new scalar degrees of freedom , neither of which is ghost - like , coupled to einstein gravity with a stable anti - desitter space as its vacuum .
we briefly present our conclusions in section [ conclusion ] .    throughout the paper our conventions are to use a metric of signature @xmath5 and define the ricci tensor as @xmath6 .
as in the case of quadratic gravity , the higher - derivative theory described by actions which are general functions of the curvature scalar is classically completely equivalent to the canonical second - order theory of a scalar field coupled to gravity .
this equivalence was first shown at the level of the action by magnano _
et al . _
we shall derive this equivalent theory in a slightly different form , by first introducing an auxiliary field to reduce the action to second order and then making a suitable conformal transformation .
we will point out the need to consider different branches of the theory when making the reduction .
we introduce the auxiliary field in two steps .
first we write @xmath7 } , \label{f(r)s}\end{aligned}\ ] ] where @xmath8 .
the auxiliary field @xmath9 has the equation of motion @xmath10 provided we have @xmath11 , this gives @xmath12 .
substituting back into the action we return to the original higher - derivative form .
thus , the reduced action is equivalent to the original theory , but only away from the critical points defined by @xmath13
. a continuous region of @xmath9 between critical points , where @xmath14 never vanishes , will be called a branch .
typically , there will be several branches in the theory .
since @xmath9 gets set equal to @xmath3 , in terms of the original theory , the condition for a critical point can be written @xmath15 .
thus , branches in @xmath9 correspond to branches in @xmath3 in the original higher - derivative theory .    as a concrete example , let us assume that @xmath16 where @xmath17 . solving @xmath13 implies that there is a single critical point at @xmath18 .
therefore , the second - order formalism in terms of the auxiliary field @xmath9 is valid in each of the two branches @xmath19 and @xmath20 , but breaks down at the critical point @xmath18 .
since @xmath9 gets set equal to @xmath3 , these regions correspond to branches in the space of curvature given by @xmath21 and @xmath22 with a critical point at @xmath23 .
the second step , which will then allow us to rewrite the action in canonical form , is to change variables to a new field @xmath24 . in terms of this new variable ,
the action becomes @xmath25 } .\ ] ] clearly this action is well defined only for those regions where the @xmath9 action is valid ; that is , away from the critical points .
therefore , the @xmath26 formulation is only defined over ranges of @xmath9 for which @xmath11 .
note that to define the action we must be able to invert the expression @xmath27 to give @xmath28 .
locally this requires the same non - degenerate condition @xmath11 .
globally , there may still be many different roots when we solve for @xmath9 in terms of @xmath26 .
however , in any given branch of @xmath9 there is only a single root .
thus there is a valid formulation of the theory in terms of @xmath26 for each branch of @xmath9 .
it is important to note that , in each branch , the inverted function @xmath29 which enters is different .    as a concrete example , consider once again @xmath1 defined in .
it follows that @xmath30 note that @xmath31 .
this expression can be inverted to give @xmath32 clearly one must take the @xmath33 root in the @xmath20 branch and the @xmath34 root in the @xmath19 branch .
this situation is quite generic , as we will see below .    for a particular reduced theory , corresponding to a given branch of the original higher - derivative theory , we can define @xmath35 and perform the conformal transformation @xmath36 .
action then becomes , @xmath37 } .\ ] ] as promised , we find a theory of einstein gravity coupled to a scalar field with a particular potential dependent on the choice of the original function @xmath1 . as in the quadratic case
there is one subtlety in defining @xmath38 . to keep @xmath38 real
, we must have @xmath39 .
when @xmath40 we must introduce , instead , the field @xmath41 , the effect of which is to change the sign of the overall normalization of the action as compared with the form given in above .
further , at the special point @xmath42 the action can not be put in canonical form . in this paper
we will restrict our attention to the @xmath43 case only .    as a concrete example , consider @xmath1 defined in . as pointed out above , in this case @xmath31 in either branch of @xmath9 .
therefore , @xmath35 is always well defined .
we point out , however , that the range of @xmath38 is restricted to @xmath44 .
note that the scalar field kinetic energy has the usual sign and , hence , @xmath38 is not a ghost .
we can now look for the vacua in any given branch of the theory .
the vacua of the theory described by action are defined to be stable solutions of the @xmath45 and @xmath38 equations of motion given by @xmath46 -            \tfrac32{\overline g}_{\mu\nu}v(\chi ) \\
\label{canr2chieom }     { \overline \nabla}^2\chi = \frac{dv}{d\chi } \end{gathered}\ ] ] respectively , where the potential is given by @xmath47 where we recall that the form of @xmath48 depends on the branch in question .
as stated previously , we will only consider vacua satisfying the covariant constant condition @xmath49 it follows that @xmath38 is a constant and , from , that it must extremize the potential .
furthermore , it follows from that the vacuum is a space of constant curvature with @xmath50 since @xmath51 and @xmath38 is constant , these vacua also correspond to spaces of constant curvature with respect to the original metric , with @xmath52 .
since the vacuum field @xmath38 must be an extremum of the potential , we have the condition @xmath53 we started with a completely general function @xmath54 , so this condition may generically have multiple solutions , including solutions away from @xmath55 .
further , given a particular solution for @xmath38 , the value of the potential at this point need not be zero ; that is , the vacuum generally has non - zero cosmological constant . in this case , the curvature scalars @xmath3 and @xmath56 will be non - vanishing .    as a concrete example of all this
, we return to the example specified in . as discussed previously , there are two branches where the second order theory is defined ; branch a where @xmath20 and branch b where @xmath19 .
the field @xmath57 , which is given as a function of @xmath9 in , satisfies @xmath31 everywhere and , hence , @xmath58 is well defined in both branches , although its range in restricted to @xmath44 .
expression was inverted to give @xmath9 as a function of @xmath57 in .
this expression is branch dependent , being given by @xmath59 with the positive root chosen in branch a and the negative root chosen in branch b. it follows that the form of the potential energy defined in also depends on the branch chosen .
it is given by @xmath60 for branch a and b respectively .
the two potentials are plotted in figures [ fig : r3vp ] and [ fig : r3vm ] .
we see that each has a stationary point .
however , since vacua must be minima of the potential , it follows that only branch b has a stable vacuum .
this is located at @xmath61 , which is non - zero .
note that @xmath62 , so that the vacuum state has a non - vanishing negative cosmological constant .
it follows that @xmath63 and @xmath64 , which implies that the vacuum state is an anti - de sitter space both in terms of the metric @xmath45 and the metric @xmath65 .
thus we have given an example of a higher - derivative theory which has a new vacuum state away from the flat - space solution of ordinary einstein gravity .
it is important to note that none of this structure would have been evident if we had made a simple expansion of the original action around flat space in terms of @xmath66 .
in fact , keeping only the first non - trivial , quadratic terms in @xmath67 , we would not even have been aware that there was an additional scalar degree of freedom in the theory .
the lowest order term in @xmath68 is cubic in @xmath67 , and so in a quadratic expansion we would only see the usual einstein term .
the conclusion is that generic higher - derivative corrections to pure einstein gravity introduce completely new vacuum states into the theory .
flat space is no longer a unique point in field space , so that if , for instance , we wish to investigate the fundamental excitations in the theory , we must start by specifying which vacuum state we are considering .    at this point
, we would like to present another specific example that further illustrates the preceding discussion .
consider the cubic function @xmath69 with @xmath70 .
we now have @xmath71 so that the theory has a single critical point at @xmath72 .
it follows that the second - order theory is defined on two branches ; branch a where @xmath73 and branch b where @xmath74 .
on either branch we can define a new field @xmath24 , which in this case is given by @xmath75 note that , since @xmath76 , @xmath57 is always a positive real number in the range @xmath77 for both branches .
expression can be inverted to give @xmath78 clearly the @xmath79 solution is correct on branch a whereas the @xmath80 solution is to be used in branch b. since @xmath57 is always positive , the definition of @xmath81 is valid in both branches .
note that @xmath38 is then restricted to lie in the range @xmath82 .
since the expression for @xmath9 in terms of @xmath57 , and hence @xmath38 , given in is different in each branch , it follows from that the potential energy function also is different in each branch .
we have @xmath83 where @xmath84 and the @xmath85 expression is valid in branch a and the @xmath86 expression valid in branch b. the two potentials are plotted in figures [ fig : r3genvp ] and [ fig : r3genvm ] . unlike the previous example , in this case each branch contains a single , stable minimum . in branch
a , this minimum is located at @xmath55 and has vanishing cosmological constant @xmath87 .
it follows that , for this vacuum , @xmath88 and spacetime is flat with respect to both the @xmath45 and @xmath65 metrics . in branch b
, the minimum is located at @xmath89 with a non - vanishing negative cosmological constant given by @xmath90 .
for this vacuum @xmath91 and , hence , @xmath64 .
it follows that the spacetime is an anti de - sitter space with respect to both metrics @xmath45 and @xmath65 .
thus we have found an example which , aside from a conventional minimum at @xmath23 with zero cosmological constant , has an additional minimum with negative cosmological constant .
furthermore , there is an unstable maximum with positive cosmological constant .
the mass of the scalar field is different at each of the different minima . here
a quadratic expansion of the action around flat - space would have identified a new scalar degree of freedom , but would never have revealed the presence of a second stable vacuum state .
it will not have escaped the readers notice that finding the vacua in the second - order formalism requires a careful discussion of the branches in the @xmath9 variable .
the branching structure was not too difficult in the preceding examples , but as we will show below , it can be , and generally is , extremely complicated .
however , the following remarks will allow us to define a simpler procedure for determining the vacua .
recall that , for a covariantly constant scalar @xmath38 satisfying @xmath92 , the corresponding spacetime structure in terms of the metric @xmath65 is that of a space of constant curvature with @xmath93 .
provided we are not at a degenerate point , spaces of constant curvature @xmath3 and vacuum solutions are in one - to - one correspondence , since the relation @xmath93 not only means that constant @xmath38 implies constant @xmath3 , but also the reverse .
note that , as a corollary , any vacuum solution with zero cosmological constant , namely @xmath94 , must correspond to flat space in the original higher - derivative theory
. an important consequence of this result is that , since flat space is a single point in the field space of the original theory , there can be only one vacuum state with zero cosmological constant in the second - order theory .
we conclude that all vacuum solutions of the above type can be found as constant curvature solutions of the equations of motion for the higher - derivative @xmath65 theory .
all such vacua can be found directly , without reference to any branch structure .
having found a vacuum of interest , one can then proceed in reverse , introducing the second - order theory in the appropriate branch containing that vacuum .
this is often an easier procedure .
the @xmath65 equations derived from the higher - derivative action are @xmath95 if we look for solutions of constant curvature , then @xmath96 where @xmath3 is a constant .
it follows that the derivative terms in the equations of motion drop out and we are left with the simple condition , first derived by barrow and ottewill  @xcite , @xmath97 we see that this is exactly the condition we obtained for stationary points of the potential aside from a factor of @xmath98 . in the latter case
, however , the expression was taken to be a function of the scalar field @xmath38 expressed in terms of @xmath9 , and as such only valid in a branch - by - branch sense , whereas here the expression is valid for all curvature @xmath3 .
this now provides us with a procedure for finding all the vacua of the theory as well as the nature of the excitations around a given vacuum .
we start by looking for constant curvature solutions of the original higher - derivative equations of motion , solving the equation , which is valid globally . having identified these points we can then , locally around each solution , make a transformation to the reduced theory , choosing , if necessary , the appropriate branch of @xmath48 . in the transformed frame , the new scalar degree of freedom
is made explicit .
we can then address questions of local stability and the mass of the scalar field .
the only breakdown of this procedure occurs when the minimum is at a degenerate point of the transformation ; that is , a point where @xmath15 . in this case
, any discussion of the particle spectrum and local stability must be in terms of the original theory .
in general , we shall not consider such degenerate points .
as an example of this procedure , consider the function @xmath99.\ ] ] where @xmath100 is a constant .
the condition for a degenerate point is then @xmath101 this equation has solutions @xmath102 where @xmath103 is any integer .
the theory thus has an infinite number of branches given by @xmath104 . rather than treating each branch separately , we follow our procedure for finding the vacuum states of the theory by solving the condition for the constant curvature solutions of the original theory .
the condition reads @xmath105 = 0.\ ] ] the left - hand side of this expression is plotted in figure [ fig : constrex ] .
we have the solutions @xmath106 where @xmath103 is an integer and @xmath107 are the solutions of @xmath108 explicitly , the first few solutions of near @xmath109 are @xmath110 furthermore @xmath111 and as @xmath112 becomes large the solutions approach @xmath113 .
the first observation to make is that the constant curvature solutions at @xmath114 , including the point @xmath23 , all correspond to degenerate points . as such
, a second - order form of the theory does not exist at these points .
thus , although these are constant curvature solutions , we can not identify them as vacuum solutions of a canonical second - order theory .
however the solutions at @xmath115 are not at degenerate points . from our previous discussion
this implies that they must correspond to vacuum solutions of the canonical second - order theory .
nonetheless , we find that each vacuum solution lies in a different branch of the theory .
let us consider two particular solutions , @xmath116 and @xmath117 , in the second - order formulation .
for the first solution we must take the branch @xmath118 .
following our usual procedure for writing the theory in a second - order form , we first introduce the auxiliary field @xmath9 and then define @xmath119.\ ] ] however , as usual , we are required to invert this expression to give @xmath9 as a function of @xmath26 .
it follows that @xmath120 we note that @xmath26 is only defined in the range @xmath121 . to put the action in canonical second - order form
we define @xmath81 .
the potential for @xmath38 , is then given by @xmath122\ ] ] where , @xmath123 we now have that @xmath38 is restricted to the range @xmath124 .
the potential is plotted in figure [ fig : sinrvp ] .
we find that there is a single maximum of the potential at @xmath125 .
the value of the potential at this point is @xmath126 , so that @xmath127 and @xmath128
. thus we see that the maximum corresponds to the flat space solution @xmath129 as expected , and we can conclude that this vacuum is unstable .    for the solution at @xmath130 , we must take the branch @xmath131 .
the expression for the potential is the same as the expression above , except that we must now define @xmath132 the resulting potential is plotted in figure [ fig : sinrvm ] .
again @xmath38 is restricted to lie in the range @xmath124 .
we now find a stable minimum at @xmath125 .
this gives @xmath133 , @xmath134 and @xmath135 , showing that the minimum of the potential does indeed correspond to the constant curvature solution @xmath130 .
we can concluded that this is a stable vacuum state .
in fact , further analysis shows that all the vacua @xmath115 with @xmath136 are unstable maxima , while those with @xmath137 are stable minima .
this final example demonstrates that the vacuum structure of higher - derivative theories can be complicated .
we find an infinite number of stable and unstable vacua , each in a different branch of the theory .    in summary
, we have shown that theories which are general functions of the scalar curvature @xmath3 , can be reduced to a canonical second - order form describing einstein gravity coupled to a scalar field with a potential which may have a vacuum structure of arbitrary complexity .
further these vacua all correspond to constant curvature solutions of the original theory .
the advantage of discussing the theory in the second - order formalism is that it is in a canonical form where we can easily identify the excitations around the vacuum .
the disadvantage is that the problem of inverting @xmath138 means that we are forced to break the theory into regions separated by degenerate points .
for this reason it is often easier to identify vacua as constant curvature solutions of the higher - derivative theory and then make a local transformation in the region of each solution to investigate the excitations around that vacuum .
all of the non - trivial vacuum structure arises from non - linear terms in the field equations and so is missed in a quadratic expansion of the higher - derivative theory in terms of the fluctuation of the metric around flat space .
further , any analysis of the particle content of the theory , in this case the presence of a new scalar degree of freedom and identification of its mass , can only be made once the vacuum is identified .
in this and the following section we shall extend our discussion to include actions which are general functions of first the ricci tensor and then the riemann tensor .
we shall use the same techniques as in the previous section , namely introducing auxiliary fields , showing that vacuum states correspond to constant curvature solutions of the higher derivative theory , and then investigating the excitations around a given vacuum .
the reduction to a second - order form was first discussed by magnano _
et al . _
@xcite and jakubiec and kijowski  @xcite . here
, we shall use a slightly different analysis , again pointing out the need for identifying different branches of the theory and also giving a canonical separation of the new variables about any given vacuum .
we shall only give the general formulation without specific examples .
again we will find a rich structure of vacua which would be completely missed in a linearized analysis .
this will re - emphasize the need to identify the vacuum in question before investigating the masses and couplings of the elementary excitations of the theory .    starting with an action
@xmath2 we can introduce an auxiliary field and put the theory in a second - order form @xmath139              } \notag \\       & = \frac{1}{2\kappa^2 } \int { d^4x \sqrt{-g } \left [               \pi^{\mu\nu } \left(r_{\mu\nu } - x_{\mu\nu}(\pi_{{\rho}{\sigma}})\right )               + f(x_{\mu\nu}(\pi_{{\rho}{\sigma } } ) ) \right ] } .\end{aligned}\ ] ] again we introduce the auxiliary field in two steps , first writing the action in terms of the field @xmath140 which gets set to @xmath141 on solving its equation of motion .
we then define @xmath142 and invert the expression to give @xmath140 as a function of @xmath143 .
note again that the introduction of both @xmath140 and @xmath143 requires the non - degeneracy condition @xmath144 if the auxiliary field is to be properly eliminated to return to the original action .
furthermore , in inverting , it may be necessary to divide the theory into branches corresponding to different possible roots for @xmath140 .
thus for the variable @xmath143 it may be necessary to introduce a collection of auxiliary variable theories each valid for a different branch of @xmath141 .
we can ask how many degrees of freedom are there in the auxiliary field ? we have argued that we expect six .
if we consider the @xmath140 formulation for now , we see first that we can again derive a spin - two divergence condition .
since the auxiliary field equation of motion gives @xmath145 , we clearly can show , using the bianchi identity @xmath146 , that @xmath147 . constraining the components of @xmath140 ,
these four conditions imply that we do indeed have six new propagating degrees of freedom .
we would like to be able to separate these into scalar and spin - two degrees of freedom as was done in the quadratic case in @xcite .
for a general function @xmath54 , this is not possible since we are not able to obtain the necessary trace condition from the @xmath65 equation of motion .
however , as we will show below , having identified a vacuum state , we can always make the separation locally in an expansion around the vacuum solution .
what then are the vacua of the @xmath2 theory ?
since we are unable to separate the spin - two and scalar degrees of freedom in the auxiliary field , we can not put the theory in a canonical form and look to minimize the potential as we did for the @xmath1 theory . however , we recall that the vacuum states in question are none other than states where the auxiliary field is covariantly constant .
this means , since the @xmath140 equation of motion gives @xmath145 that they are states of constant tensor curvature .
we can also impose the condition that the states should have no preferred direction , so that @xmath140 is proportional to the metric @xmath65
. thus vacuum states are states of constant scalar curvature @xmath148 with @xmath3 constant .
consequently , one approach to finding vacuum solutions is simply to look for constant curvature solutions of the original theory .
the equations of motion of the original higher derivative theory are given by @xmath149 to incorporate the presence of the metric in @xmath2 in deriving the equations of motion , we consider @xmath54 as a function of the mixed index object @xmath150 , with all contractions made between raised and lowered indices so that the metric does not enter explicitly .
then in the above expression @xmath151 , so that , for instance , @xmath152 . looking for constant curvature solutions of the form @xmath96 , we obtain the condition @xmath153 where @xmath154 and @xmath155 , and we have used the fact that @xmath156 evaluated at the constant curvature solution must be proportional to @xmath157 . note that this condition has exactly the same form as the condition we obtained for actions which were general functions of the scalar curvature .    given a particular vacuum solution satisfying the condition , we would like to investigate the excitations around the vacuum state .
the natural way to do this is to consider an expansion of the action around a given vacuum @xmath158 .
expanding @xmath2 and keeping terms to quadratic order in the curvature only , we have @xmath159 } \notag \\       & = \frac{\beta}{2\kappa^2 } \int { d^4x \sqrt{-g } \left [            - \tfrac12 r_0 + r + \frac{1}{6{m_0}^2}r^2            - \frac{1}{{m_2}^2}\left(r_{\mu\nu}r^{\mu\nu }                -\tfrac13r^2\right ) \right ] } ,   \label{expandf}\end{aligned}\ ] ] here we have evaluated @xmath54 and its derivative at the vacuum @xmath160 and defined @xmath161 using the general symmetric decomposition of the second two expressions .
the variables @xmath162 , @xmath163 and @xmath164 are then given by @xmath165 where in writing the final line of we have used the constant curvature condition @xmath166 .
we see that we have succeeded in putting the action in a quadratic form , though with the addition of a cosmological constant .
the final line of is identical to the general quadratic form we discussed in our previous paper @xcite , except for the addition of a cosmological constant term and a renormalization of the gravitational coupling constant by a factor @xmath162 .
we can thus introduce auxiliary fields @xmath38 and @xmath167 exactly as we did in @xcite to give @xmath168 . \end{gathered}\ ] ] the only effect of the cosmological term is to modify the potential for the scalar field , adding a term @xmath169 .
we find that the extremum of the potential is now at @xmath170 , which , relating the auxiliary field @xmath38 back to the original curvature , gives @xmath158 , as required for a consistent expansion .
following exactly the analysis of our previous paper , the auxiliary field @xmath167 satisfies generalized divergence and trace conditions , so it does indeed describe the degrees of freedom of a spin - two field . by making a final field redefinition
, we can write the action in canonical form with explicit kinetic energy terms for the spin - two field . again following the discussion in our previous paper
, the spin - two field will have the correct pauli - fierz limit , but will unfortunately be ghost - like .
thus we have shown that locally , around any vacuum of the higher - derivative theory ( that is a solution of constant curvature ) , we can expand the theory to identify a new scalar and a new spin - two degree of freedom ( provided we are not at a degenerate point ) . the spin - two field satisfies divergence and trace conditions as before but importantly ,
we see that it remains ghost - like .
thus generalizing to @xmath2 actions fails to remove the problem of the ghost spin - two degree of freedom .
the masses of the degrees of freedom are fixed by the form of the function @xmath54 around the constant curvature solution .
it is worth noting that we could equally well have done this analysis in the second - order form , looking for solutions with constant @xmath140 proportional to @xmath65 .
then expanding in small @xmath140 about such solutions , gives a linear coupling between @xmath140 and @xmath141 and quadratic ` mass ' terms for @xmath140 . from this form
we could then extract the scalar and spin - two parts of @xmath140 . in this sense ,
the auxiliary field @xmath140 always carries six degrees of freedom , which about any given vacuum can be decomposed into a scalar field and a ghost spin - two field , with the form of the decomposition changing as we go from vacuum to vacuum .
our last generalization is to consider actions which are general functions of the riemann tensor @xmath171 . as mentioned earlier , the suggestion is that such theories have an additional six degrees of freedom , which we know in the quadratic case can be decomposed into scalar and spin - two fields .
we start the discussion by demonstrating that , as in all previous cases , we can introduce an auxiliary field to write the action in a second order form .
@xmath172 } \notag \\       & = \frac{1}{2\kappa^2 } \int { d^4x \sqrt{-g } \left [               \rho^{{\lambda}\mu\nu{\rho } } \left(r_{{\lambda}\mu\nu{\rho } }                    - x_{{\lambda}\mu\nu{\rho}}(\rho_{\alpha\beta\gamma\delta})\right )               + f(x_{{\lambda}\mu\nu{\rho}}(\rho_{\alpha\beta\gamma\delta } ) ) \right ] } .\end{aligned}\ ] ] again we introduce the auxiliary field in two stages ; first introducing @xmath173 , which gets set equal to the riemann tensor on solving its equation of motion , and then defining @xmath174 in both cases we require the non - degeneracy condition @xmath175 in order to be able to eliminate the auxiliary field and return to the original action .
when using the variable @xmath176 we may be required to break the theory into branches , introducing a collection of second - order theories , each taking a different branch when inverting to find @xmath177 in terms of @xmath176 .    turning to the vacuum states
, we again find that states with constant @xmath177 have constant riemann curvature , since @xmath177 is set equal to @xmath178 by its equation of motion .
furthermore , imposing the condition that there is no preferred direction in spacetime we require that @xmath177 is proportional to @xmath65 .
therefore , given the symmetries of @xmath177 , we have @xmath179 and , hence , we are considering solutions of constant ricci scalar curvature .
as before the easiest way to obtain such solutions is from the equations of motion of the original higher - derivative action .
again , to circumvent the problem , when deriving the equation of motion , of the metric explicitly entering the function @xmath54 , we consider @xmath54 as a function of @xmath180 with all contractions made between raised and lowered indices .
we then derive the equations of motion @xmath181 where @xmath182 . restricting to constant curvature solutions of the form , we get the familiar condition @xmath183 where now @xmath184 and @xmath185 , and we have used the fact that , by symmetry , @xmath186 is proportional to @xmath187 .    to investigate the excitations around a given vacuum we expanding the action about the constant curvature solution .
we have , keeping terms up to quadratic order in the curvature only , @xmath188 \right )               \left ( r^{{\lambda}\mu\nu{\rho } } - \tfrac{1}{12}r_0                  \left[g^{{\lambda}\nu}g^{\mu{\rho}}-g^{{\lambda}{\rho}}g^{\mu\nu}\right ] \right )             \right ] \notag \\       & = \frac{\beta}{2\kappa^2 } \int { d^4x \sqrt{-g } \left [             - \tfrac12 r_0 + r + \frac{1}{6{m_0}^2}r^2             - \frac{1}{{m_2}^2}\left(r_{\mu\nu}r^{\mu\nu }                -\tfrac13r^2\right )           \right . } \notag \\ & \qquad\qquad \left .
+ \gamma\left(r_{{\lambda}\mu\nu{\rho}}r^{{\lambda}\mu\nu{\rho } }                -4r_{\mu\nu}r^{\mu\nu}+r^2\right ) \right ] .
\label{riemannexpand}\end{aligned}\ ] ] here we evaluate @xmath54 and its derivatives at the constant curvature solution , defining , for the contraction of the derivatives with a tensor @xmath189 which has the symmetries of the riemann tensor but in otherwise arbitrary , @xmath190 with @xmath191 and @xmath192 .
we have also introduced the parameters @xmath193 and have used the constant curvature condition @xmath166 in the final line of .
the final expression in is a quadratic action with a gauss - bonnet term , a weyl - squared term and a ricci - scalar - squared term , together with a cosmological constant .
the gauss - bonnet term can be dropped classically as a total divergence , leaving the action in the same quadratic form as discussed in our previous paper @xcite . reducing
the action to a second - order form then follows exactly as in the case for @xmath2 actions , so that we obtain the same transformed action .
again we can verify that the transformed action has a vacuum solution at @xmath158 as is required for the consistency of our expansion .
we also note that here too the expansion could have been made in terms of the variable @xmath173 , which could then be decomposed into its scalar and spin - two parts , the form of the decomposition depending on which vacuum is being considered .    in conclusion ,
general actions of the form @xmath194 may have a variety of vacuum solutions , generically not apparent in a linear analysis . around any vacuum the new degrees of freedom in the theory , aside from the massless graviton , can always be separated into a scalar field and a spin - two field .
unfortunately , the spin - two field is always ghost - like .
in this section , we shall briefly discuss how to extend our analysis to actions with higher - order equations of motion .
we have seen that an action involving any function of the curvature tensors gives at most fourth - order equations of motion . for higher - order equations we need to consider actions which include some derivatives of the curvature . here
, we shall be concerned with only the simplest form , @xmath195 where we require that the function @xmath54 has been reduced , by integration by parts and dropping total derivatives , so as to minimize @xmath4 .
the equations of motion following from similar actions have been considered by schmidt  @xcite and wands  @xcite , who showed them to be equivalent to those for a set of scalar fields coupled to gravity . here
will shall show the full equivalence in a new way , working at the level of the action
.    it will be important to distinguish between possible forms of the function @xmath54 .
if we write @xmath196 , ( where in the action we have @xmath197 , @xmath198 and so on ) , we find there are two possible cases : either @xmath199 is a function of @xmath200 ( case one ) or it is not ( case two ) .
if it is not then it must be a function of @xmath201 , since otherwise the original form was reducible ; that is , the action was not written in a form minimized with respect to @xmath4 .
thus , in case two , we can always decompose @xmath54 as @xmath202 further , we see that for case one the equations of motion are @xmath203th - order , while for case two they are @xmath204th - order .
we would like to reduce the action to a canonical second - order form , by introducing auxiliary fields . given the order of the higher - derivative equations of motion we see that in case one we must introduce @xmath205 new fields , while in case two we need only @xmath206 new fields .
the procedure we shall use is essentially a generalization of ostrogradski s method for reducing a higher - order action to a first - order form @xcite , only that here we shall be reducing to a second - order form .
( the ostrogradski result is usual given as a hamiltonian , but this can always be rewritten as a helmhotz lagrangian , the analog of the form we shall use . )    we start by introducing a set of lagrange multipliers , so @xmath207.\end{gathered}\ ] ] clearly eliminating fields via the lagrange multiplier equations of motion , starting with @xmath208 and working down in @xmath4 returns one to the original action .
however , we have introduced at least one too many new fields .
it is clear that the @xmath200 equation of motion is purely algebraic and eliminating it will not introduce higher - derivatives ; that is , the action will remain second order .
the @xmath209 equation of motion reads @xmath210 where we have distinguished between the two cases discussed above , and substituted the special form of @xmath54 in case two . in case one , to form the analog of ostrogradski s lagrangian
, we solve the equation to give @xmath200 as a function of @xmath211 and @xmath208 , writing @xmath212 .
it should be noted that in general the solution is not unique , and we must divide the original theory into pieces corresponding to different branches of the solution , just as in the case of actions of the form @xmath1 discussed in the previous section . in case two
we simply substitute for @xmath208 and the special form of @xmath54 .
we get in case one @xmath213 , \end{gathered}\ ] ] the exact analog of the ostrogradski - helmhotz lagrangian , while for case two we have a slightly different form , @xmath214 .
\end{gathered}\ ] ] we see that , as expected , in case one we have a total of @xmath205 auxiliary fields , while in case two we have only @xmath206 new fields .
all that remains is to transform the action into a canonical form , growing canonical kinetic energy terms for all the new auxiliary fields .
let us concentrate on actions of case one .
terms of the form @xmath215 are easy to deal with . simply introducing a pair of new fields , @xmath216 and @xmath217 by @xmath218 we have , in the action , @xmath219 }                \notag \\        & = \int { d^4x\sqrt{-g } \left [                - \left(\nabla\chi_i\right)^2 + \left(\nabla\psi_i\right)^2                \right ] } .
\end{aligned}\ ] ] we see that the @xmath216 field has a canonical kinetic term , but the @xmath217 field has the wrong sign ; it is a ghost .
this is characteristic of higher- order theories . in reducing the theory to second - order
the new fields always enter as a pair of a ghost - like field with a ordinary field .
if we define a potential function , @xmath220 where it is understood that the right - hand side is evaluated at @xmath221 , @xmath222 , the action in case one becomes @xmath223 } .\ ] ] to complete the transformation to canonical form all that is left is to make a conformal rescaling of the metric to remove the @xmath224 coupling . as usual
we define @xmath225 and rescale to a metric @xmath36 , giving @xmath226   .\end{gathered}\ ] ] we conclude that the original case one higher - derivative gravity theory is equivalent to canonical einstein gravity coupled to @xmath205 scalar fields , @xmath227 of which , @xmath38 and @xmath216 , propagate physically and @xmath4 of which , @xmath217 , are ghost - like .
to put the case two action in canonical form is more complicated because of the @xmath228 term .
however , in principle , it is always possible to introduce a set of new fields @xmath229 which simultaneously diagonalize the kinetic terms for the @xmath230 and @xmath231 , though now the form of the transformation will depend on the function @xmath232 .
at least @xmath233 of the new fields will be ghost - like .
we can then make a conformal rescaling as in case one to put the action in the same canonical form , though the potential function will have a different form .
the conclusion is that there is a procedure for rewriting the higher - order action in a canonical second - order form .
however , @xmath4 of the new fields in case one and at least @xmath233 of the fields in case two will be ghost - like . in each case
we obtain a specific potential , and so we can again look for vacuum states as stationary points of the potential .
generically , as in the case of @xmath1 actions discussed in the previous section , the non - trivial vacua do not correspond to the flat - space solution of the original higher - derivative theory .    as an simple example of this procedure consider the function @xmath234 .
this example is of the first case since writing @xmath235 with @xmath197 and @xmath198 , we have @xmath236 , which is not independent of @xmath237 . following our general procedure ,
the first step is to introduce a set of lagrange multipliers , @xmath238 } \notag \\       & = \frac{1}{2\kappa^2 } \int { d^4x \sqrt{-g }            \left [ { \lambda}_1 + \alpha{{\lambda}_2}^2            + \mu \left(r-{\lambda}_1\right ) + \mu_1 \left(\nabla^2{\lambda}_1-{\lambda}_2\right )           \right ] } .
\end{aligned}\ ] ] as discussed above we have introduced one too many auxiliary fields .
we can eliminate @xmath237 by solving its equation of motion , which reads @xmath239 implying @xmath240 .
substituting back into the action gives @xmath241 } .\ ] ] next , to put the kinetic energy for @xmath242 and @xmath243 in canonical form , we define @xmath244 and @xmath245 , so that @xmath246 } .\ ] ] finally we make the conformal rescaling @xmath36 with @xmath225 to put the action in canonical form @xmath247 . \end{gathered}\ ] ] thus we see that the original higher - derivative gravity theory is equivalent to canonical einstein gravity coupled to three scalar fields .
one field is ghost - like , and the potential has a single , unstable stationary point at @xmath248 .    as a second example consider @xmath249 writing @xmath197 and @xmath198 , we now have @xmath250 , independent of @xmath237 , so this example clearly falls under case two . repeating our procedure , we introduce lagrange multipliers to give @xmath251 }           \notag \\       & = \frac{1}{2\kappa^2 } \int { d^4x \sqrt{-g } \left [            \alpha + \beta{\lambda}_1 + \gamma{{\lambda}_1}^2            + \epsilon{\lambda}_1{\lambda}_2 + \mu \left(r-{\lambda}_1\right )            + \mu_1 \left(\nabla^2{\lambda}_1-{\lambda}_2\right ) \right ] } .
\end{aligned}\ ] ] again we have too many new fields .
the @xmath237 equation of motion now reads @xmath252 which can not be solved for @xmath237 since the action is case two . however substituting this solution into the action eliminates both @xmath237 and @xmath243 , leaving the second - order form , @xmath253 } .\ ] ] the kinetic energy for @xmath242 is already in canonical form , but we must make a final conformal rescaling by @xmath254 to put the action in the completely canonical form @xmath255 } , \ ] ] where we have the potential , @xmath256 again the original higher - derivative gravity is shown to be equivalent to ordinary einstein gravity though now coupled to only two scalar fields . by choosing @xmath17
we can ensure that neither field is ghost - like .
note that we argued above that case two theories must have at least @xmath233 ghost - like fields .
it is thus only in this special case of @xmath257 that we are able to have all the scalar fields non - ghost - like .    as before we obtain a specific potential for the fields .
we find that @xmath258 has a single stationary point at @xmath259 , @xmath260 , provided @xmath261 . expanding around this point
we find that , if @xmath262 and @xmath263 , we have a stable minimum .
the value of the potential at the minimum is @xmath264 which is negative .
we conclude that , for the given range of @xmath265 , @xmath162 and @xmath266 , the theory has a single stable vacuum state with negative cosmological constant . from einstein
s equation we find that this state is an anti - desitter space with @xmath267 .
using the fact that @xmath52 for covariantly constant @xmath38 , we find that @xmath268 .
thus the vacuum state also corresponds to an anti - desitter space in the original higher - derivative theory .
the most important conclusion of this paper is that higher - derivative theories of gravitation generically have multiple stable vacua .
one of these may be trivial , corresponding to flat spacetime , but all the other vacua are non - trivial with the associated manifold being either desitter or anti - desitter spacetime with non - vanishing cosmological constant . while of interest from various points of view , such non - trivial vacua can not represent the universe as it is now since the radius of curvature of these solutions is of order the inverse planck mass .
thus , one might conclude that non - trivial gravitational vacua are irrelevant for particle physics theories .
however , this is not the case .
we have recently shown that if we extend the methods and results of this paper to the realm of @xmath269 supergravity , then non - trivial vacua can exist with vanishing cosmological constant @xcite .
furthermore , we find that supersymmetry is generically spontaneously broken in these vacuum states .
it follows that higher - derivative @xmath269 supergravitation could play a pivotal role in high energy physics .
this possibility is being pursued elsewhere @xcite .
this work was supported in part by doe grant no .  de - fg02 - 95er40893 and nato grand no .
crg-940784 .
a.  hindawi , b.  a. ovrut , and d.  waldram , _ soft supersymmetry breaking induced by higher - derivative supergravitation in the electroweak standard model _ ,
381b * ( 1996 ) , 154162 , hep - th/9602075 .


superclusters are the largest structures in the universe .
their sizes span the range from a few times the radius of a typical galaxy cluster ( @xmath3 3 - 4 @xmath4 ) to the lengths of the `` great walls '' of galaxies observed by the cfa redshift survey @xcite and the sloan digital sky survey @xcite .
these latter sizes of @xmath5 represent a significant fraction ( @xmath6% ) of the horizon scale .
because they are so large relative to a typical survey volume , superclusters have proven difficult to study , especially in a statistical manner .
observationally , attempts to characterize the large - scale structure using distributions of galaxies date back to @xcite , but these methods did not become powerful until the advent of galaxy cluster and galaxy redshift surveys .
@xcite used the palomar sky survey to construct a catalog of rich clusters , with which he was able to find evidence for what he called `` second - order clusters , '' i.e. clusters of clusters , which are indeed superclusters .
subsequent work included @xcite s study of the supercluster containing the rich coma cluster , and @xcite s work on the perseus supercluster . other observational approaches to supercluster identification
have been based , like that of @xcite , on cluster catalogs ( e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) or on analysis of smoothed density fields @xcite .
these have been important in characterizing the low - redshift universe , but they have not yet allowed the study of supercluster evolution from early cosmic times ( i.e. , from redshifts of order unity ) .
how do superclusters form , and how do their shapes , multiplicities , sizes , and structures evolve with time ? the answers to these fundamental questions are not yet known .
recently , there has been exciting progress on both the observational and computational fronts .
the sloan digital sky survey ( sdss , * ? ? ?
* ; * ? ? ?
* ) is collecting deep , broadband imaging and spectroscopy of a large , contiguous portion of the sky .
sdss has created the largest galaxy redshift survey catalog to date , and much work is also underway to improve photometric redshift techniques ( e.g. * ? ? ?
* ; * ? ? ?
* ) , which would permit measurements of 3-d galaxy positions for the much larger sdss imaging catalog . on the theoretical side , a concordance cosmological model has gained widespread acceptance @xcite . applying this model to cosmological dynamics codes such as that described below
, we can now reliably predict the properties and evolution of superclusters .
the time is therefore ideal for a computational analysis of large - scale structure using new , powerful simulations , the results of which can be used to determine the properties and evolution of superclusters . in addition , these results can be compared with observations , thus providing a new test of the current cosmological model
complementary to that provided by cmb , large - scale power spectra , and other studies  as well as shed light on structure formation and evolution .
we begin by describing , in  2 , the process by which clusters and in turn superclusters were identified in the simulation output .
we provide quantitative descriptions of the multiplicity and size distributions of the identified superclusters in  [ sec : mult ] and  [ sec : size ] .
we examine the radial structure of our superclusters in  [ sec : radial ] and the dimensionality of superclusters in  [ sec : shapes ] .
the evolution of the above properties with redshift is presented in  [ sec : evo ] .
results for superclusters selected from only higher - mass clusters are presented in
[ sec : him ] .
our conclusions are summarized in
[ sec : conc ] .
we use the efficient , parallel dark matter simulation code of @xcite . a large - scale
, high - resolution simulation was evolved with a particle - mesh method to compute long - range gravitational forces , and a tree code to treat high - density regions .
we use the cosmological parameters of the concordance model : @xmath7 and @xmath8 ( cf .
@xcite ; see also @xcite ) .
a total of @xmath9 particles , each assigned the mass @xmath10 , is evolved in a periodic box @xmath11 on a side , and positions are saved at the appropriate times in order to construct a light cone with vertex at one corner of the box .
clusters are then identified in two different , partially overlapping volumes : a `` low - redshift '' sphere of radius @xmath11 , corresponding to a maximum redshift of @xmath12 , and a `` high - redshift cone , '' which includes the low-@xmath13 sphere s positive octant ( @xmath14 ) and extends this octant out to @xmath15 , a distance of @xmath16 . using the cluster selection scheme described below ,
a total of 1,442,616 clusters with @xmath17 are found in the low-@xmath13 sphere , and 720,550 in the high-@xmath13 cone .
further details on this light cone simulation are given by @xcite , hereafter hbb .
clusters are identified in the simulation using the friends - of - friends ( fof ) algorithm ( see hbb ; * ? ? ?
the fof results depend on the chosen linking length @xmath18 .
a common parameterization uses the linkage parameter @xmath19 , where @xmath20 is the mean interparticle separation .
this is useful in applications where the mean density of objects , and in turn their mean separation , changes over time . holding @xmath21 constant imposes the same minimum overdensity for clusters at all redshifts .
we use the value @xmath22 ( see also hbb ) , corresponding to an overdensity of @xmath23 times the mean density .
the center of each fof cluster is taken to be the position of the most bound particle .
this center is used to compute the virial radius @xmath24 , defined as the radius enclosing the virial overdensity expected from spherical top - hat collapse ( which evolves with redshift in the @xmath0cdm cosmology ) . in this analysis
, we consider only clusters with virial mass threshold ( i.e. , the total mass within the virial radius ) @xmath17 .
this mass threshold is typical for a poor cluster of galaxies .
the mean mass of all clusters in the low-@xmath13 sphere is @xmath25 . in
[ sec : him ] we discuss superclusters selected from clusters with a higher mass threshold .
the ellipticities of the clusters identified at low and high @xmath13 are studied in detail by hbb .
relevant to our present work on superclusters , they find alignments between nearby cluster pairs out to separations of @xmath26 ; the strength of alignment increases with redshift and with decreasing cluster separation .
hbb explain these alignments as being due to clusters forming and gathering material along `` filamentary superstructures '' ; they provide evidence for overdensities along the separation vectors between pairs of aligned clusters .
as demonstrated by @xcite , @xcite and references therein , superclusters can be identified by their overdensity of clusters .
we use the fof method to identify superclusters from the cluster sample defined above (  2.1 - 2.2 ) .
the mean density of clusters in the low-@xmath13 sphere ( we shall not discuss the high-@xmath13 cone until  [ sec : evo ] ) is @xmath27 , yielding a mean intercluster separation of @xmath28 .
we therefore select linking lengths of @xmath29 and @xmath30 , corresponding to @xmath31 and @xmath32 , respectively .
all scales in this paper are in comoving coordinates .
the different linking lengths sample different portions of the supercluster geometry .
small linking lengths identify only the dense cores of superclusters , whereas larger linking lengths allow the fof algorithm to percolate to lower - density filaments of clusters .
large @xmath18-values therefore reach to the outer parts of superclusters , as well as identify looser superclusters , and often connect structures that appear discrete with smaller @xmath18 .
we use monte carlo simulations to compare the significance of the selected superclusters with that of random associations .
clusters are placed at randomly chosen positions within a sphere such that the total number density is the same as that for the simulated clusters , and the fof method is applied , as before , with @xmath29 and @xmath30 .
we find significantly fewer clusters per supercluster , and significantly fewer superclusters down to multiplicity as small as 3 , for all four linking lengths tested .
this result implies that most of the superclusters are real , even for @xmath33 .
the only exception is the case of binary clusters ( superclusters with only 2 members ) selected with @xmath33 , of which we find only @xmath34% more using the simulated clusters than we do with the monte carlo sample .
the random rate of superclusters is typically @xmath35% for @xmath36 superclusters , decreasing to @xmath37% for @xmath38 superclusters .
the monte carlo results are compared with the @xmath0cdm multiplicity functions below .
to help visualize the typical geometries of superclusters identified by the fof method , we present supercluster maps of a typical region ( a cube of side length @xmath39 , concentric with the low-@xmath13 sphere ) for @xmath40 ( figs@xmath41[fig1]-[fig2 ] ) . each point is a cluster , and a given symbol represents membership in the same supercluster .
the maps present a projection of the 3-d distribution onto the @xmath42-plane in the first plot , and onto the @xmath43-plane in the second plot .
viewing the plots in turn corresponds to a @xmath44 rotation of the viewing angle about the @xmath45-axis . for clarity
, we show only superclusters with at least 5 member clusters . members of a supercluster that extend outside the limits of our cubical region are not shown , so some superclusters appear to have fewer members than the imposed minimum . comparing these maps to analogous figures generated with @xmath46 and @xmath30 reveals , as expected ,
that increasing the linking length often causes neighboring superclusters to be connected into one .
this is especially true when @xmath18 increases from 7 to @xmath30 , resulting in many large , complicated structures with central cores joining multiple lower - density filaments .
the _ multiplicity _ of a supercluster is the number of clusters it contains , @xmath47 . in fig@xmath41[fig3 ] we present , for various linking lengths , the integrated supercluster multiplicity function  the number density ( per cubic comoving @xmath4 ) of superclusters with @xmath47 or more members , as a function of @xmath47 . also shown for comparison
are the corresponding results from random catalogs , as described in the previous section .
longer linking lengths yield more superclusters of all multiplicities , as expected .
shorter linking lengths yield curves that rise more steeply toward low multiplicities , implying that they preferentially find small superclusters with low multiplicity .
conversely , longer linking lengths find a higher percentage of larger , richer superclusters .
the shorter linking lengths naturally trace high - density regions ( e.g. , dense cores of superclusters ) , while larger linking lengths characterize lower - density regions , such as the outer parts of superclusters or diffuse superclusters .
small @xmath18 can also divide structures that are identified with larger linking lengths .
the fact that we see superclusters with @xmath48 or more clusters may seem surprising , but the total number density of such structures is only @xmath49 per @xmath50 , and it is based on a very low mass threshold of cluster members , @xmath17 .
such clusters are numerous , and any comparison with observations should of course account for this threshold . for superclusters
selected using higher - mass clusters , see  [ sec : mevo ] .
we repeated the above analysis after selecting superclusters in redshift - space coordinates ; the resulting multiplicity functions are nearly identical to the real - space functions at low multiplicities , and are shifted to slightly higher @xmath47 ( by up to 20% ) at the highest multiplicities .      the supercluster size function  i.e. , the number of superclusters above a given size  is presented in fig@xmath41[fig4 ] for different linking lengths .
the size of a supercluster is defined as the maximum distance of any member cluster from the supercluster center of mass . as with the multiplicity function ,
the size function rises more steeply for shorter linking lengths .
this implies that the longer linking lengths find a greater ratio of large superclusters to small superclusters , analogous to the result for the multiplicity function .
similar results are found using redshift - space coordinates ; for the largest linking length used , the sizes are slightly ( @xmath51 ) larger .
we find that the largest supercluster radii are @xmath52 , corresponding to diameters of @xmath53 .
this is consistent with the largest structures seen in observations ( e.g. * ? ? ?
* ; * ? ? ?
* ) , and it is also the largest distance at which hbb detect cluster alignment in the simulation .      the average multiplicity of superclusters is plotted as a function of supercluster size for different linking lengths in fig@xmath41[fig5 ] . for @xmath54 or @xmath55
 objects more aptly called binary or triple clusters than superclusters  the maximum radius is set mainly by the separation(s ) of the 2 or 3 clusters ; thus the average multiplicity changes slowly over this range of separations . beyond this regime ,
the multiplicity and size of superclusters appear to be related nearly linearly .
this result suggests that the spherically averaged density profile in superclusters is roughly _ isothermal _ ; i.e. , @xmath56 , where the cluster number density @xmath57 satisfies : @xmath58 this profile yields multiplicity @xmath59 for @xmath56 .
we can investigate the spherically averaged density profile more directly by examining the distribution of clusters within individual superclusters and determining the average number density profile for the entire sample .
to do so , we superpose the centers - of - mass ( calculated using only the mass in clusters ) of all superclusters , and scale all superclusters to the same size by dividing all cluster coordinates ( relative to the center - of - mass ) for a given supercluster by its @xmath60 .
we add up the total number of clusters inside radial bins , divided by the total number of clusters in superclusters for the given linking length .
we have included only superclusters with 5 or more members ; because the multiplicity function drops off rapidly , most of the superclusters included have only 5 or 6 members .
this fact affects the combined profile in two ways .
first , only in few cases will one of the five members happen to lie very near the supercluster s center - of - mass ; thus it is only at @xmath61 that we begin to find a significant fraction of the clusters .
second , one out of roughly five clusters in each supercluster lies at @xmath62 .
thus , for the smaller linking lengths , the integrated number fraction approaches a limit of @xmath63 ; for the larger linking lengths , which find richer superclusters , the limiting fraction is slightly higher .
the fraction of clusters within @xmath64 increases linearly with radius for the range @xmath65 .
this linear relationship indeed implies @xmath66 for the supercluster density profile in this range .
the profile steepens on larger scales .
the average mass of clusters in superclusters , @xmath67 for the @xmath40 superclusters , is somewhat larger than the average mass of clusters in the total sample , @xmath68 , indicating that superclusters contain slightly more massive clusters , on average , than the mean .
the @xmath40 superclusters with @xmath69 members have an even higher average cluster mass of @xmath70 , implying that richer superclusters contain more massive clusters than poorer superclusters .
we measure the shape , and in particular the dimensionality , of a given supercluster by fitting a 3-d ellipsoid to the distribution of its member clusters .
we use the code developed by hbb to fit ellipsoids to the dark matter particles comprising their clusters .
the code constructs , for a given supercluster , the @xmath71 matrix of second moments of cluster positions relative to the supercluster center - of - mass , i.e. @xmath72 where @xmath73 and @xmath74 are two of the position coordinates for a given cluster , and @xmath75 is the mass of the cluster ; the sum is over all member clusters .
the eigenvalues of this matrix are simply the three axis lengths of the best - fit ellipsoid for the supercluster , times a known constant .
following the notation used by hbb , we shall denote these lengths @xmath76 for the primary , secondary , tertiary axes , respectively ; i.e. , @xmath77 .
we restrict our analysis to superclusters with @xmath69 ; for poorer superclusters , the `` best - fit ellipsoid '' would not be very meaningful .
we have already examined size distributions for superclusters (  [ sec : size ] ) , so we now focus on the _ axis ratios _ , which allow us to probe the shape and dimensionality of superclusters .
fig@xmath41[fig6 ] displays the distribution of _ primary axis ratios _
@xmath78 for superclusters selected with @xmath79 .
fig@xmath41[fig7 ] shows the distribution of _ secondary axis ratios _ @xmath80 , and fig@xmath41[fig8 ] shows @xmath81 .
the mean and peak axis ratios from each curve in these figures are listed in table 1 .
the secondary axis ratios are generally smaller than the primary axis ratios .
we find that a primary axis ratio of @xmath82 is typical for superclusters selected using a broad range of linking lengths ; this corresponds to a primary ellipticity of @xmath83 .
the secondary axis ratio is more strongly dependent on the linking length used to find the superclusters . for @xmath84
, we find a peak @xmath80 of @xmath85 , corresponding to a secondary ellipticity of @xmath86 ; for @xmath33 , we find a peak of @xmath87 , or @xmath88 .
the peak @xmath81 ratios range from @xmath89 for @xmath33 to @xmath90 for @xmath84 .
these measurements allow us to construct a picture of the typical shape and dimensionality of superclusters .
the interpretation of different ellipticities can be described as follows :    1 .
a combination of high primary and low secondary axis ratios is typical of two - dimensional pancake - like structures ( two dimensions are large , and one is small ) .
a combination of low primary and high secondary axis ratios is typical of one - dimensional filamentary structures ( one dimension is large , and two are small ) .
a combination of primary and secondary axis ratios that are both significantly less than unity is typical of structures most appropriately described as triaxial ( each dimension is a different size ) .
4 .   if both axis ratios are close to unity , then the distribution is nearly spherical .
the results we find suggest that superclusters are typically triaxial structures nearing 2-d , pancake - like structures at @xmath91 .
they are neither spherical nor purely filamentary in nature . as the linking length
is increased , the resulting superclusters become more two - dimensional ; the cores of superclusters are thus more triaxial , spreading out to more pancake - like structures on larger scales .
we plot the bivariate distribution of @xmath78 and @xmath80 for @xmath40 in fig@xmath41[fig9 ] ( the distributions for @xmath92 and @xmath93 are similar ) . applying the interpretation rules listed above ,
1-d structures reside in the upper left corner of fig@xmath41[fig9 ] , 2-d in the lower right corner , and 3-d ( spherical ) in the upper right .
points in the middle or toward the lower left corner represent triaxial structures .
fig@xmath41[fig9 ] shows that a large majority of the superclusters are more nearly two - dimensional than one - dimensional .
no nearly spherical superclusters are found with this linking length , but there is a smaller population of filamentary structures .
similar distributions are found for @xmath92 and @xmath93 , though showing more flattened structures ( 2d ) as @xmath18 increases , and more triaxial structures as @xmath18 decreases .
this is the same trend depicted in fig@xmath41[fig7 ] as a shift toward lower secondary axis ratio with increasing linking length .
the shorter linking lengths also yield significantly higher fraction of filamentary structures  i.e. , superclusters with low primary axis ratio and a secondary axis ratio approaching unity .
this trend suggests that most superclusters at @xmath94 are triaxial structures with one dimension considerably smaller than the other two .
this geometry is apparent in fig@xmath41[fig9 ] , and becomes more significant when a larger fof linking length is used .
since the superclusters are not completely flat , smaller linking lengths may probe only the dense core regions , which are smaller in size than the supercluster thickness .
thus , while the thickness sets an upper bound to the tertiary axis length identified with any linking length , the higher @xmath18-values find longer primary and secondary axes by probing farther out into the pancake .
this causes larger linking lengths to find a smaller mean secondary axis ratio , and reveal the underlying 2-d nature of the large - scale superclusters .
similarly , for a given @xmath18 , the largest , highest - multiplicity superclusters ( @xmath95 ) are on average flatter ( @xmath96 for @xmath97 ) than those with @xmath98 5 - 10 ( @xmath99 ) .
this flatness is still apparent in projection : for superclusters selected ( from @xmath17 clusters ) in redshift - space , the projected axis ratios have a broad distribution with a mean @xmath100 and a peak @xmath101 , for all linking lengths @xmath1025 , 7 , and 10@xmath4 .
the triaxiality of superclusters can also be quantified using the measurement commonly employed for geometrical descriptions of elliptical galaxies , first introduced explicitly by @xcite : @xmath103 with @xmath77 , @xmath104 will take on values between 0 and 1 ( unless @xmath105 , in which case @xmath104 is undefined ) .
two - dimensional pancake structures approach @xmath106 , while filaments approach @xmath107 ( note , however , that the increase is not linear with axis ratio ) .
intermediate values of @xmath104 represent triaxial systems , with triaxiality increasing with @xmath104 . by this measure
our low-@xmath13 superclusters are triaxial , ranging from a mean of @xmath108 for @xmath33 to @xmath109 for @xmath84 .    for all linking lengths
we see a tail in the distribution that describes structures which are very nearly one - dimensional .
thus , some mass concentrations are filamentary at @xmath91 , especially in the relatively dense regions probed by @xmath84 .
indeed , we would expect structures that have collapsed along two dimensions to be of higher average density than those that have collapsed along only one dimension .
this result agrees with hbb s identification of filamentary density enhancements in the large - scale structure ; here we add the prediction of more diffuse flattened structures surrounding and possibly connecting the filaments .
we examine the time evolution of the supercluster properties studied above .
we study superclusters in two redshift slices at @xmath110 and @xmath111 , denoted @xmath112 and @xmath113 , respectively .
superclusters are identified by the fof method as before . however , since the cluster abundance decreases at higher redshift , and their mean separation thus increases , we compare superclusters found at different redshifts _ not _ with the same linking length @xmath18 , but instead with the same linkage parameter @xmath19 , where @xmath20 is the mean comoving intercluster separation . as discussed in  2.2
, this approach sets a minimum overdensity for detection of superclusters , which is constant with redshift .    this use of a constant linkage parameter accounts for the reduced density of clusters at high redshift .
however , because clusters are clustered , the density of superclusters depends on the relative rates of cluster formation in regions of high density vs. low density .
hierarchical supercluster formation , in which clusters gravitate toward each other over time , also contributes to a reduced supercluster density at early times .
we indeed find a lower density of superclusters in the high-@xmath13 cone .
the total number of superclusters decreases by a factor of a few from @xmath114 to @xmath112 , and by over an order of magnitude from @xmath112 to @xmath113 . in order to construct a sufficiently large sample ,
we consider broad ranges in redshift , especially for the @xmath113 slice .
in addition , we omit the @xmath21-value corresponding to the @xmath115 ( @xmath114 ) case from our subsequent analysis , as it results in too few superclusters at high redshift .    as stated in  2.3 , our @xmath114 linking lengths of @xmath79 correspond to @xmath116 , respectively .
we use these latter linkage parameters @xmath21 to refer to the linking lengths used in our subsequent analysis .
a given value of the @xmath21-parameter can be converted to an actual comoving linking length at a given redshift simply by multiplying by the mean intercluster separation at that redshift .
the mean cluster densities in our @xmath112 and @xmath113 shells are @xmath117 and @xmath118 , corresponding to mean cluster separations of @xmath119 and @xmath120 , respectively . therefore , @xmath121 0.23 , 0.33 , and 0.47 correspond to @xmath102 7 , 10 , and 14.3@xmath4 at @xmath122 , and 16.8 , 24 , and 34.3@xmath4 at @xmath123 .
the evolution of the multiplicity function is presented in fig@xmath41[fig10 ] .
we present the integrated multiplicity function for all three redshifts , using two different values of @xmath21 .
the number density of superclusters drops noticeably with increasing redshift .
monte carlo simulations analogous to those described in  2.3 were also run at @xmath112 and @xmath124 , with results similar to those obtained at low @xmath13 : the multiplicity functions are physically significant for all multiplicities and linking lengths , except for binary clusters in the @xmath125 case , where only @xmath126 more binary clusters are found over the monte carlo results .
the redshift evolution of the multiplicity function depends on the linking length of the superclusters .
the vertical distances between the @xmath127 curves in fig@xmath41[fig10 ] increase with increasing multiplicity , implying that the abundance of the richest ( highest-@xmath47 ) superclusters found with this linking length drops off especially swiftly with increasing @xmath13 .
however , for @xmath128 the vertical distances _ decrease _ with increasing @xmath47 , to the extent that the abundance of the richest superclusters found with this linking length is in fact roughly the same at @xmath112 as it is today .
we interpret these trends as further evidence that small linking lengths ( e.g. @xmath128 ) select the densest regions of superclusters . whereas the larger linking lengths do not find many rich superclusters at high redshift because of the lack of clusters , the smaller linking lengths find the same dense supercluster cores even at much earlier times .
these dense cores are likely the sites of the largest overdensities present at early times , and the first places where observable structures form through gravitational collapse .
thus smaller linking lengths relative to the mean intercluster separation are needed to locate the highest overdensities in the observable universe .
the lack of rich superclusters at high @xmath13 is also evidence that these structures form largely through gravitational accretion of clusters that are initially too distant to be identified as supercluster members .
supercluster size remains highly correlated with multiplicity at high @xmath13 , but with a lower @xmath129 normalization than at low @xmath13 .
the size function does not evolve significantly with redshift for the larger superclusters ( see fig@xmath41[fig11 ] ) , while , as expected , the number of cluster members decreases .      using the method described in
[ sec : shapes ] , we fit ellipsoids to the superclusters identified at @xmath112 and @xmath113 . because we only perform these fits on superclusters with @xmath69 , which are rare at high redshift , the statistical uncertainties are larger . for example , only 68 superclusters are identified at @xmath113 for @xmath128 .
we use the axis ratios of the best - fit ellipsoids to characterize the evolution of supercluster dimensionality .
tables 2 and 3 present the mean primary and secondary axis ratios , respectively , for @xmath116 at @xmath130 .
we also construct the axis ratio distributions ( similar to figs@xmath41[fig6]-[fig7 ] ) as a function of @xmath13 for @xmath131 and @xmath32 ( there are too few @xmath128 superclusters at high redshift to create meaningful distributions ) .
the results for @xmath127 are shown in figs@xmath41[fig12]-[fig13 ] .
the number of high - redshift superclusters is too low to construct bivariate distributions as in fig@xmath41[fig9 ] , but we can consider how the changes in _
mean _ ratios correspond to translations of these distributions .
the strongest evolution observed is the increase in the secondary axis ratio with increasing redshift , for @xmath125 ( fig@xmath41[fig13 ] ) .
most of the change occurs between @xmath114 and @xmath112 , with only a slight additional increase between @xmath112 and @xmath113 ( see table 3 ) .
typical secondary axis ratios increase from @xmath132 at low @xmath13 to @xmath133 at high @xmath13 .
thus the flattened structures observed with @xmath125 at @xmath114 ( see  [ sec : shapes ] ) appear less flattened at high redshifts , indicating their gravitational collapse over time .
very little evolution is evident in the secondary axis ratio for superclusters identified with @xmath131 , but the population of high - density filaments with @xmath134 evolves between @xmath112 and @xmath114 , as there are almost none of these filamentary structures at high redshift .
the @xmath135 distributions for @xmath136 are shown in figs@xmath41[fig6]-[fig7 ] ; at higher redshifts the distributions are similar to the corresponding curves in figs@xmath41[fig12]-[fig13 ] .
the primary axis ratios are slightly smaller at higher redshifts for both @xmath131 and @xmath32 .
typical values decrease from @xmath82 at @xmath114 to @xmath137 at @xmath138 1 - 2 . in a plot similar to fig@xmath41[fig9 ] , this would correspond to a leftward shift in the distribution , i.e. a shift toward greater triaxiality at high redshift , and a shift to a 2-d geometry as the systems approach @xmath114 .
quantitatively , our high-@xmath13 superclusters have mean triaxiality @xmath139 for both linking lengths , significantly higher than the mean @xmath104 for any linking length at low @xmath13 .
we conclude that superclusters had highly triaxial structure at early times , and that between a cosmological look - back time of @xmath140 gyr ( corresponding to @xmath112 for the cosmology described in  2.1 ) and the present day , many collapsed along their smallest dimensions to form the more nearly pancake - like structures observed in the simulation at @xmath114 .
in the same span of time , some structures in regions of highest density also collapsed along a second dimension to form nearly one - dimensional filaments .
how do the supercluster properties change if superclusters are selected using higher - mass clusters as their seeds ?
while the abundance of high - mass clusters is significantly lower , and the statistical uncertainties larger , the massive clusters are more easily detected in observational surveys .
we repeat the analysis described above for superclusters selected using the rarer , higher - mass clusters with @xmath141 .
this threshold corresponds to typical `` rich clusters '' ( e.g. * ? ? ?
* ; * ? ? ?
there are @xmath142 clusters above this mass threshold at @xmath91 , and @xmath143 clusters at @xmath122 .
the cluster abundance is @xmath144 at @xmath91 , decreasing to @xmath145 at @xmath122 . the intercluster mean separation is @xmath146 at @xmath91 , increasing to @xmath147 at @xmath122 .
superclusters are selected from these clusters at @xmath91 and @xmath122 using the linkage parameters @xmath136 and @xmath32 ; these correspond to @xmath148 and @xmath149 respectively at @xmath91 , and @xmath150 and @xmath151 at @xmath122 .
the large linking lengths @xmath18 reflect the large mean intercluster separation of these massive clusters .
we also include , for comparison , the smaller linking length of @xmath97 at @xmath91 .
there are fewer superclusters identified , and a smaller number of cluster members in each , as compared to those selected using the low - mass clusters .
the overall structure of the superclusters is not significantly changed , as described below .
the multiplicity function of the superclusters is presented in fig@xmath41[fig14 ] for both @xmath91 and @xmath152 .
the results show , as expected , a similar shape to the multiplicity function of superclusters selected using the more numerous low - mass clusters ( figs@xmath41[fig3 ] , [ fig10 ] , [ fig11 ] ) , but with a greatly reduced amplitude : the number of superclusters at a given multiplicity is lower , and
equivalently  the supercluster multiplicity ( richness ) is smaller for a given supercluster abundance .
this is partly due to the much smaller number of massive clusters .
the evolution of the multiplicity function is stronger than that of superclusters selected with low - mass clusters ( cf .
fig@xmath41[fig10 ] ) .
superclusters with linkage parameters of @xmath153 contain up to @xmath154 rich clusters at @xmath91 .
the supercluster multiplicity decreases sharply at higher redshifts .
the supercluster size function is presented in fig@xmath41[fig15 ] .
the maximum size of @xmath91 superclusters ranges from @xmath155 ( i.e. , twice the plotted @xmath60 radius ) for @xmath97 , to @xmath156 for @xmath157 ( @xmath136 ) , to @xmath158 for @xmath159 ( @xmath127 ) . for a given linkage parameter @xmath21 ,
the maximum supercluster size is somewhat larger than that of superclusters selected with low - mass clusters , since the linking length values @xmath18 are larger . for a given value of @xmath21 , no significant evolution
is observed in the size of the largest superclusters , reflecting the existence of these extended structures at early times .
the same result is observed for superclusters selected with the lower mass threshold clusters .
the number of cluster members within these large structures decreases with redshift , as seen above , due to the lowered cluster abundance at high @xmath13 .
the supercluster multiplicity versus size relation is similar  but reduced in amplitude , as expected  to the superclusters selected with the lower mass threshold ( fig@xmath41[fig5 ] ) . a nearly linear relation between multiplicity and size
is observed for all but the binary clusters .
the ellipticity distribution shows that superclusters exhibit a triaxial geometry on average , but with significantly larger scatter than the more numerous superclusters selected with lower - mass clusters ( figs@xmath41[fig12]-[fig13 ] ) .
the mean primary and secondary axis ratios are each @xmath160 to @xmath161 ; for example , for @xmath136 at @xmath91 , the mean values are @xmath162 and @xmath163 .
these ratios do not change significantly with either linking length or with redshift to @xmath122 .
the distribution of ellipticities around this mean is , however , quite broad .
the supercluster shapes are not significantly different from those selected using the lower - mass clusters , except for the stronger flattening observed in the latter at @xmath91 , where the mean triaxial shape has evolved slightly towards a flattened two - dimensional shape ( with @xmath164 and @xmath165 ) . the superclusters selected with high - mass clusters
are triaxial , with somewhat higher triaxiality at @xmath91 than those selected using lower - mass clusters  i.e .
, the high - mass superclusters are approaching a slightly more filamentary nature , with mean axis ratios of @xmath160 in both @xmath78 and @xmath80 ( cf .
fig@xmath41[fig9 ] ) , and a mean triaxiality of @xmath166 .
we have investigated the multiplicities , sizes , radial structures , and dimensionalities of superclusters from @xmath114 to @xmath113 in @xmath0cdm cosmology . the results
improve our understanding of the properties ( especially shapes ) and evolution of superclusters , and can also be used as predictions for testing the cosmology with future observations .    at low redshift , we predict a total number density of @xmath167 , and @xmath168 superclusters with at least five member clusters with mass @xmath17 for superclusters selected with linking length @xmath169 , and @xmath93 , respectively .
the maximum size of these superclusters ranges from @xmath170 for @xmath97 to @xmath171 for @xmath172 .
the abundance of superclusters decreases rapidly with increasing redshift .
there are many more superclusters at low @xmath13 than at high @xmath13 largely because there are many more clusters that have formed by the present time and that have gravitated to form superclusters .
the richest superclusters are almost entirely missing at high redshifts ; these structures form late in cosmic time . however , the densest `` core '' regions of superclusters , traced well by fof linking lengths @xmath173 of the mean intercluster separation , are present at @xmath122 almost as abundantly as they are today .
these are likely the sites of the highest overdensities in the early universe .
the spherically averaged density profiles of superclusters are well fit by an isothermal profile , @xmath174 , over a broad range in radius .
the profile is shallower at small radii and is steeper at large radii . a nearly linear relation exists between supercluster size and multiplicity , arising from this density profile .
we find that superclusters are triaxial in shape , especially at early times , where the mean triaxiality of our sample is @xmath139 . over time , many collapse along a single dimension to approach two - dimensional shapes on the largest scales ( i.e. , those sampled by the fof linking length @xmath33 at @xmath114 , using low - mass clusters with @xmath17 ) ; their cores remain triaxial .
quantitatively , the mean triaxiality for long linking lengths decreases to @xmath108 at @xmath114 , whereas for smaller @xmath18 ( i.e. , in the cores ) it decreases to @xmath109 . the mean primary axis ratio of superclusters at @xmath114 is @xmath175 , corresponding to a primary ellipticity @xmath176 .
typical secondary axis ratios at @xmath114 are @xmath177 and @xmath178 for @xmath84 , and @xmath179 and @xmath180 for @xmath33 .
gravitational collapse increases ellipticity , as the potential gradient is larger along the minor rather than the major axis @xcite ; collapse thus leads to lower dimensionality ( cf . * ? ? ? * for a more recent treatment ) .
a nonnegligible population of dense , filamentary superstructures is also present at @xmath114 , though not at higher redshift .
superclusters selected using higher - mass clusters , @xmath181 , show consistent results , with a reduced multiplicity function , as expected .
these superclusters are triaxial in shape , with somewhat greater mean triaxiality at @xmath91 ( @xmath182 0.8 - 0.85 ) than superclusters with lower - mass clusters
.    these properties , derived from a large - scale , high - resolution cosmological simulation , provide direct information on the size , content , shape , and evolution of superclusters , as well as testable predictions of the currently favored @xmath0cdm model with which future and currently accumulating observations can be compared .
once sdss is complete , the best statistical analyses to date of large - scale structure using observations will be possible , especially if accurate and well - behaved photo-@xmath13 algorithms are developed .
there are other exciting observing programs on the horizon .
the pan - starrs project @xcite is expected to see first light in 2006 , and will carry out multi - band imaging surveys to identify galaxy clusters , as well as gravitational weak lensing observations to map the matter distribution in the universe .
the large synoptic survey telescope @xcite , which is expected to see first light in 2012 , will provide improved weak lensing maps that will probe the dark matter distribution to the largest scales .
the analysis presented here will aid in the comparison of these major observational data with the predictions of the concordance cosmological model .
this research was supported in part by nsf grant ast-0407305 .
computations were performed on the national science foundation terascale computing system at the pittsburgh supercomputing center ; additional computational facilities at princeton were provided by nsf grant ast-0216105 .
abell , g. o. 1958 , apjs , 3 , 211 bahcall , n. a. 1988 , ara&a , 26 , 631 bahcall , n. a. & soneira , r. m. 1984 , apj , 277 , 27 bahcall , n. a. , ostriker , j. p. , perlmutter , s. , & steinhardt , p. 1999 , science , 284 , 1481 bennett , c. l. , et al .
2003 , apjs , 148 , 1 bode , p. , & ostriker , j. p. 2003 , apj , 145 , 1 budavri , t. , szalay , a. s. , connolly , a. j. , csabai , i. , & dickinson , m. e. 2000 , aj , 120 , 1588 csabai , i. , connolly , a. j. , szalay , a. s. , & budavri , t. 2000 , aj , 119 , 69 davis , m. , efstathiou , g. , frenk , c. s. , & white , s. d. m. 1985 , apj , 292 , 371 einasto , m. , einasto , j. , tago , e. , dalton , g. , & andernach , h. 1994 , mnras , 269 , 301 einasto , j. , et al .
2003a , a&a , 405 , 425 einasto , j. , et al .
2003b , a&a , 410 , 425 einasto , j. , et al .
2006 , a&a , submitted ( astro - ph/0603764 ) erdogdu , p. , et al .
2004 , , 352 , 939 geller , m. j. , & huchra , j. p. 1989 , science , 246 , 897 gott , j. r. , et al . 2005 , apj , 624 , 463 gregory , s. a. , & thompson , l. a. 1978 , apj , 222 , 784 gregory , s. a. , thompson , l. a. , & tifft , w. g. 1981 , apj , 243 , 411 hopkins , p. f. , bahcall , n. , & bode , p. 2005 , apj , 618 , 1 hubble , e. 1936 , the realm of the nebulae ( new haven : yale university press ) kaiser , n. , et al .
2002 , proc .
spie , 4836 , 154 kolokotronis , v. , basilakos , s. , & plionis , m. 2002 , mnras , 331 , 1020 lin , c. c. , mestel , l. , & shu , f. h. 1965 , apj , 142 , 1431 navarro , j. f. , frenk , c. s. , & white , s. d. 1996 , apj , 462 , 563 navarro , j. f. , frenk , c. s. , & white , s. d. 1997 , apj , 490 , 493 oort , j. h. 1983 , ara&a , 21 , 373 ostriker , j. p. , & steinhardt , p. j. 1995 , nature , 377 , 600 peacock , j. a. 1999 , cosmological physics ( cambridge , uk : university press ) spergel , d. n. , et al .
2003 , apjs , 148 , 175 spergel , d. n. , et al . 2006 , astro - ph/0603449 statler , t. s. 1994 , aj , 108 , 111 stoughton , c. , et al .
2002 , aj , 124 , 1810 tyson , j. a. 2002 , proc .
spie , 4836 , 10 york , d. g. , et al .
2000 , aj , 120 , 1579 yoshisato , a. , morikawa , m. , gouda , n. , & mouri , h. 2006 , apj , 637 , 555    rrrrrrr @xmath18 & mean @xmath78 & peak @xmath78 & mean @xmath80 & peak @xmath80 & mean @xmath81 & peak @xmath81 + 5 & 0.60 @xmath183 .009 & 0.64 & 0.50 @xmath183 .007 & 0.41 & 0.29 @xmath183 .004 & 0.23 + 7 & 0.57 @xmath183 .004 & 0.62 & 0.43 @xmath183 .003 & 0.31 & 0.23 @xmath183 .001 & 0.19 + 10 & 0.54 @xmath183 .002 & 0.64 & 0.33 @xmath183 .001 & 0.20 & 0.16 @xmath183
.001 & 0.13 +    rrrrrrr @xmath13 & @xmath184 & @xmath136 & @xmath127 + 0 & 0.60 @xmath183 .009 & 0.57 @xmath183
.004 & 0.54 @xmath183 .002 + 1 & 0.44 @xmath183 .03 & 0.47 @xmath183 .01 & 0.48 @xmath183 .01 + 2 & 0.52 @xmath183 .06 & 0.50 @xmath183 .03 & 0.50 @xmath183 .02 +    rrrrrrr @xmath13 & @xmath184 & @xmath136 & @xmath127 + 0 & 0.50 @xmath183 .007 & 0.43 @xmath183 .003 & 0.33 @xmath183 .001 + 1 & 0.36 @xmath183 .02 & 0.40 @xmath183 .01 & 0.44 @xmath183 .01 + 2 & 0.41 @xmath183 .05 & 0.41 @xmath183 .03 & 0.46 @xmath183 .01 +