Cosmic structure sizes in generic dark energy models

The maximum allowable size of a spherical cosmic structure as a function of its mass is determined by the maximum turn around radius $R_{\rm TA,max}$, the distance from its centre where the attraction on a radial test particle due to the spherical mass is balanced with the repulsion due to the ambient dark energy. In this work, we extend the existing results in several directions. (a) We first show that for $w\neq -1$, the expression for $R_{\rm TA, max}$ found earlier using the cosmological perturbation theory, can be derived using a static geometry as well. (b) In the generic dark energy model with arbitrary time dependent state parameter $w(t)$, taking into account the effect of inhomogeneities upon the dark energy as well, where it is shown that the data constrain $w(t={\rm today})>-2.3$, and (c) in the quintessence and the generalized Chaplygin gas models, both of which are shown to predict structure sizes consistent with observations.


Introduction and overview
Since the discovery of the accelerated expansion of our current Universe at the end of last century, the dark energy, along with the cold dark matter have become the central themes in the theory and applications of modern cosmology. The simplest of the dark energy models is certainly a positive cosmological constant, Λ and the simplest model of the modern universe is Λ-Cold Dark Matter (ΛCDM). The ΛCDM has so far passed with flying colors all cosmological observation tests [1,2], starting from the redshift of type Ia Supernovae, the Hubble rate, the galaxy clustering, the microwave background radiation and so on 1 .
Since the observed value of Λ ∼ O(10 −52 m −2 ) is tiny compared to the inverse of the Planck length squared, it is natural to expect that its effect would be observable at large space-time scales only, as for instance imprinted in the light coming from some high redshift supernovae or in the data for the early universe. However, a novel potentially local check of the dark energy was proposed recently [4,5,6], pertaining the stability of large scale structures. The idea is simple, based on the observation that the maximum size of an e.g. spherical bound structure should be the distance from its center at which the attractive Newtonian force on a test mass due to the spherical overdensity is balanced with the repulsion due to the ambient dark energy. Beyond this distance, which was called "maximum turnaround radius" and labeled R TA,max a test particle cannot stay bound, but will be dragged away by the antigravity effect of the dark energy. The same idea was used in [7] to derive an upper bound on a hypothetical at that time cosmological constant on the basis of the existence of galaxies. The next step is then to use R TA,max as an observable to constrain the parameters of any cosmological model of interest, by comparing its theoretical prediction for R TA,max with the actual data.
For ΛCDM, in particular, the predicted value of R TA,max for a spherical structure equals (3M G/Λc 2 ) 1/3 , the details of which would be presented in the next sections. Fig. 1 [4] (see also references therein), shows the size versus mass of some nearby large scale structures, together with the theoretical prediction of the model. As is evident, for superclusters as large as M 10 15 M ⊙ , the prediction of ΛCDM lies very close from above, with the departure from the data being roughly only about 10%! Thus, the ΛCDM model is absolutely consistent with the stability of cosmic structures from the maximum turn around perspective 2 . It is worth noting here that the structures we are looking into are very close to us (z ≪ 1), and hence this is a local check or challenge to the dark energy -the dark energy is at work right within the structure itself. After a certain structure decouples from the expansion of the surrounding Universe, it shrinks and virializes by redistributing the kinetic energies and positions of its individual constituents [9] under the action of the attractive force due to its mass, which is now dominant. Given that the process of virialization of a structure typically enhances any non-sphericity in its initial profile, a spherical model for a structure is more appropriate when it is away from virialization. Furthermore, it was shown using the Press-Schechter mass function in [6,10] that there exists a transitional mass scale ∼ 10 13 M ⊙ , above which the structures are not virialized today. Thus, in order to challenge a dark energy model using the maximum size versus mass criterion, it is most effective to look into structures which lie above that transitional mass scale. We refer our reader to [11] for a discussion on the future evolution of the cosmic webs of such large masses and to [12,13,14] for a discussion on the possible violation of the maximum turn around bound by looking into the peculiar velocity profiles of members of a given structure. We further refer our reader to e.g. [15,16] on aspects of particle motion in the black hole Schwarzschild-de Sitter spacetime relevant to the derivation of R TA,max in ΛCDM and to [17,18,19] for effects of a positive Λ on astrophysical scenarios or phenomena such as the accretion disks.
Despite the simplicity and the most overwhelming successes of the ΛCDM model, it is well known that it suffers from the so called "fine tuning of the cosmological constant" problem, the fact that Λ, when interpreted and computed as the vacuum energy density, most outrageously disagrees with its actual observed value today [20]. Furthermore, it does not provide any insight towards an explanation of the "coincidence problem", the fact that the current numerical values of the vacuum and matter energy densities are so close to each other. Attempts to resolve these issues on the basis of the instability of de Sitter space [21,22,23,24] are not yet entirely satisfactory. Finally, in [25], a possible discrepancy between the ΛCDM model and the Baryon Acoustic Oscillation data for high redshift has been proposed. All these issues and also the so far lack of any observational evidence of a dark matter particle candidate, have triggered in recent times vigorous research for alternatives of the ΛCDM model and/or Einstein's theory of gravitation [26,27,28] (also references therein). As a result, a plethora of potential candidate theories or models have been proposed and the need to invent effective observables to distinguish them has become apparent. See also [29,30] for a recent proposal of quantum modification of the General Relativity, by treating the massive gravitons as the primary dark matter.
Since its proposal, the maximum turn around radius as a useful cosmological observable has received  considerable attention in the context of modified or alternative gravity/dark energy models [31]- [37]. The goal is of course to constrain the additional parameters of such theories. For example, for a McVittie spacetime [38,39,40], derivation of the turn around radius can be found in [31]- [33], by analyzing geodesics and also using the quasilocal mass function of [41,42]. In [34] a constraint on the matter-galileon coupling parameter was obtained for a cubic galileon model [43]. Both an upper and a lower bound was obtained for this coupling. The lower bound is new, while the upper is improved by about 50% from what was known on the basis of the Solar system [44] data. The derivation of R TA,max in the Brans-Dicke theory [45] in the presence of a positive cosmological constant Λ was carried out in [35,36]. In this theory the R TA,max is greater than that of ΛCDM. Qualitatively, this is due to the fact that the extra scalar of the Brans-Dicke model contributes to make the attractive force stronger, thus allowing for larger structures. Consequently, no constraint on its parameter ω was obtained. Two more models with the same effective behavior will be presented in Section 3.2.
There are two ways to derive the maximum turn around radius -one is to analyze geodesics in a static geometry and the other is via the study of a test fluid moving in a background cosmological McVittie spacetime [4,5]. It has been explicitly demonstrated that both these approaches yield the same result for the ΛCDM model [5], as well as for the Brans-Dicke theory with a positive Λ [36].
Let us now emphasize the new results we derive in the rest of this paper. For w = −1, the dark energy is taken to be time dependent and as we mentioned earlier, for a constant such w, the expression for R TA,max was derived earlier in [5] using cosmological scalar perturbation theory. Can we have a derivation of the same result using a static geometry as well, just like w = −1? The answer turns out to be, Yes. In Sec. 2, we consider in the context of General Relativity a static and spherically symmetric dark energy fluid with equation of state P (r) = wρ(r) with constant w. We derive a first order metric using this source and the maximum turn around radius using this metric is found to match exactly with the one derived in [5].
In many alternative dark energy models, the state parameter, w, could be time dependent. Therefore, it is important to have a derivation of R TA,max in the context of such dark energy models as well. In Sec 3.1 we shall present a derivation of the R TA,max in the cosmological scenario for a dark energy with an arbitrary, time dependent state parameter w(t) (P E (t) = w(t)ρ E (t)) and by also taking into account the effect of inhomogeneities upon the dark energy distribution. We shall further apply in Subsection 3.2 this generic result in the context of the quintessence, as well as of the generalized Chaplygin gas models, establishing their consistency with the observational data. We conclude in Section 4 with a brief discussion and outlook.
We shall use mostly positive signature (−, +, +, +) for the metric and will set c = 1 throughout.

R TA,max in the static geometry
It is well known that the positive Λ permits a static and spherically symmetric solution within the cosmological event horizon, namely, the Schwarzschild-de Sitter spacetime, where dΩ 2 = dθ 2 + sin 2 θdφ 2 is the unit sphere line element and H 0 = Λ/3. Let us consider a timelike geodesic u a (u a u a = −1) in this spacetime. The above spacetime is endowed with one timelike and three spacelike Killing vector fields generating a 2-sphere, The conserved quantities along the geodesic corresponding to each of these Killing vector fields are given by, where the 'dot' denotes differentiation with respect to the proper time along the trajectory. The above conserved quantities are respectively identified as the conserved energy and the conserved components of the orbital angular momentum of the geodesic. Expanding the on-shell condition u a u a = −1 in the background of (Eq. (1)) and using Eqs. (Eq. (3)), we havė where L 2 = L 2 1 + L 2 2 + L 2 3 = r 4 (θ 2 +φ 2 sin 2 θ) and the ± sign denotes outgoing and incoming trajectories, respectively. The maximum turn around is determined by the condition of vanishing accelerationr = 0, or, It is straightforward to convince oneself that, as expected, the consequence of angular momentum is to decrease the size of a stable structure with a given mass, because the attraction of the central mass has to counterbalance the additional centrifugal force due to the rotation. Indeed, define and write the condition (Eq. (5)) in the form The cosmic structures of interest have sizes much larger than their Schwarzschild radius and much smaller than the cosmological horizon, i.e. 2M ≪ r ≪ H −1 0 or equivalently m ≪ x ≪ m −1/2 , with m ≪ 1 for consistency. For λ = 0, (Eq. (7)) has the well known solution . Introducing a small angular momentum λ 2 = O(ǫ) the above solution changes to x = 1 − 1 3 (1 − m)λ 2 < 1, which confirms that the angular momentum diminishes the size of the cosmic structures. The above argument can be generalized numerically to all physical angular momenta.
Thus we conclude that the maximum possible size of a structure with a given mass is the maximum turn around radius for L = 0, that is In other words, this is precisely the point where the attraction due to the central mass gets balanced with the repulsion due to the dark energy, for radial timelike geodesics.
The above derivation could also be performed using the orbits of the timelike Killing vector field of (Eq. (1)). Precisely, the static observers along these timelike Killing vector field would feel a spacelike acceleration Thus the acceleration vanishes at r = R TA,max , corresponding to the maximum of the norm of the timelike Killing vector field, giving the maximum size of the structure. We also note here that for structures up to as large as M ∼ 10 16 M ⊙ and taking Λ ∼ 10 −52 m −2 , the (2M G/r + H 2 0 r 2 ) term appearing in the metric functions is much less than unity at this length scale. In other words, the maximum we are dealing with here is a very broad maximum for practical astrophysical scenarios.
Even though the above computation might seem idealistic as the framework used was static, it can be generalized, as we shall see in the next section, to include the ΛCDM cosmological perturbation theory and leads to exactly the same answer [5], as expected from the possible gauge or coordinate invariance of the maximum turn around radius [46]. We shall also further extend this result there with general dark energy with equation of state P E (t) = w(t)ρ E (t).
Staying at the constant w case, for w = −1, the dark energy is in general treated as time dependent and spatially homogeneous. For a constant such w, the maximum turn around radius was found to be (−3M G/(4π(1 + 3w)ρ E )) 1/3 [5]. For w = −1, we may identify Λ = 8πGρ E to recover the ΛCDM result. In the following, we wish to present a novel alternative derivation of this result in the framework of a static geometry, assuming accordingly that the dark energy is static as well.
The concept of static dark energy was earlier introduced in e.g. [47]- [51] to construct compact and static star-like solutions, known as the "dark energy stars", made up of a dark energy fluid. This corresponds to the assumption that in the vicinity of a structure or central mass, the dark energy becomes inhomogeneous too. Quite surprisingly, as we shall see, the answers obtained in these two ways match! However, it is important to emphasize that we are not considering any dark energy stars here, but only generic large scale structures in the ambiance of a static dark energy fluid. Certainly, as we move away from the object, the assumption of staticity would break down.
Let us then begin with the ansatz for a static spherically symmetric spacetime, We assume that the energy momentum tensor for the dark energy corresponds to that of an ideal fluid, where ρ E , P E and u a are respectively the energy density, pressure and the velocity of the fluid's world line. For a backreacting T E ab in a static geometry, u a should be parallel to the orbits of the timelike Killing vector field. If we normalize it to unity, u a u a = −1, we may take u a = f −1/2 (∂ t ) a . We assume that the parameter w in the fluid equation of state P E = wρ E , is a constant. We shall take ρ E > 0, so that w < −1/3, necessary to violate the strong energy condition and generate repulsive effects [52].
The three independent Einstein equations for this system are where a 'prime' denotes derivative with respect to the radial coordinate once and where use was made of the equation of state. Using the conservation equation ∇ a T E ab = 0, we obtain where ρ E,0 is an integration constant. The above equation can be thought of as an analogue of the time dependence of the energy density in a homogeneous cosmological spacetime, where the scale factor a(t) is equivalent to the redshift or the Tolman factor f (r) [52], for a static spacetime. We could not find an exact solution of (Eq. (10)) corresponding to the above expression for ρ E (r). However, since we are chiefly interested in the maximum turn around region, we still could find an approximate and linearized solution there as follows. Since we have seen earlier that the maximum turn around radius in a static spacetime corresponds to the maximum of the norm of the timelike Killing vector field (f ′ = 0), we may ignore the variation of ρ E around this point and write ρ E ≈ ρ E (r ∼ R TA,max ) ≈ ρ E,0 /f (R TA,max ) ≡ ρ 0 = Λ/8πG (say). Then, the solution of (Eq. (10)) for the metric is The parameter M can be regarded as the mass of a central compact object and r 0 is an arbitrary length parameter. The three other constants r H , r C and r U are the nonvanishing roots of the equation 1 − 2M G/r − Λr 2 /3 = 0. For the physically interesting case 3M G √ Λ ≤ 1, all of them are real, namely are precisely the horizon lengths of the Schwarzschild-de Sitter spacetime, while the third r U = −(r H + r C ) is clearly unphysical. However, for w = −1 they are not known to have any such physical meaning. Note that for w = −1 in Eq. (12), we recover, as expected, the Schwarzschild-de Sitter spacetime. Also, in the limit M → 0, it is easy to see that r H → 0 and r C → 3 Λ , in which case Eq. (12) becomes Finally, if we set Λ = 0 in Eq. (12), we recover the Schwarzschild spacetime. However, since the metrics Eq. (12) or Eq. (14) were derived using the a priori assumption that we are in the maximum turn around region, to be consistent, we have to linearize it now as follows. First we note that since the observed value of the dark energy is tiny, it is reasonable to assume that M G √ Λ ≪ 1.
Then it turns out that r H ≈ 2M G and r C ≈ 3 Λ in Eq. (13). Moreover, the maximum turn around length scale must be much smaller than the Hubble horizon scale ∼ O(Λ −1/2 ) and also much larger than the Schwarzschild radius of the central mass, i.e. both terms GM/r and Λr 2 /3 must be much smaller than unity at or around this length scale. Putting these all in together, we find, up to linear order the metric in Eq. (12) can be written around the maximum turn around region as, Rescaling the time coordinate t ′ = (r 0 /2M G) 3(1+w)/4 t, and dropping off the prime without any loss of generality, we obtain finally which for w = −1 leads to the linearized Schwarzschild-de Sitter spacetime with Λ identified with the cosmological constant.
As we discussed earlier, the maximum turn around radius in a static spacetime is equivalently given by the maximum of the norm of the timelike Killing vector field (i.e., ∂ r g tt = 0), giving where in the last equality we have substituted back Λ = 8πGρ E,0 /f (R TA,max ) ≡ 8πGρ 0 . This exactly matches with the result of [5], obtained by using the cosmological perturbation theory, and rules out all dark energy models with w −2.3. For the cosmological case, ρ 0 was to be understood as the dark energy density today for nearby structures. Finally, as a consistency check of Eq. (15), it is easy to see that in the maximum turn around region for structures as large as up to M ∼ 10 16 M ⊙ , both M G/r and Λr 2 terms are much smaller than unity, justifying further the simplification ρ E (r) ≈ ρ 0 we used for the derivation. 3 .

Time dependent cosmological scenario
3.1 R TA,max in models with generic state parameter w(t) We start with the ansatz for the perturbative McVittie spacetime, describing the gravitational field of a point mass sitting in the ambient spatially homogeneous and isotropic cosmological spacetime with flat spatial sections, where Ψ and Φ are the weak and linearized gravitational potentials and R = x 2 + y 2 + z 2 is the comoving radius. The above ansatz holds only for length scales much greater than the Schwarzschild radius of the central mass. In other words, it is suitable for describing large scale non-black hole structures such as the galaxies or clusters of galaxies, whose Schwarzschild radii are located much inside them, implying |Ψ|, |Φ| ≪ 1. A rather familiar example of such spacetimes would be the Schwarzschild-de Sitter spacetime written much outside the Schwarzschild radius, with a(t) = e √ Λ/3 t and Ψ(R, t) = −2GM/Ra(t) = Φ(R, t). The Christoffel symbols for the above metric are 3 As an aside, we note here that Eq. (9), Eq. (11) give us the expression for the trace of the static dark energy's energymomentum tensor, (r) Thus on the event horizon of a static black hole (f (r) = 0), the trace diverges for w < −1 which implies the divergence of the Ricci scalar too, indicating a naked curvature singularity. Since the horizon of a static black hole corresponds to the infinities of the timelike Killing coordinate (e.g. [52]), perhaps the above divergence is an analogue of the so called phantom disaster (e.g., [1]) at cosmological late times for such dark energy models.
where a 'dot' denotes differentiation once with respect to time, H =ȧ/a is the Hubble rate, the indices (i, j) run between the spatial coordinates x, y and z and we have retained terms only linear in the potentials. The nonvanishing components of the Ricci tensor are where ∂ 2 ≡ ∂ 2 x + ∂ 2 y + ∂ 2 z , whereas the components of the Einstein tensor are, We will use the above equations to do study cosmological scalar perturbation theory at the linear order with the dark energy and the cold dark matter as the sources. We shall ignore the backreaction effects due to any other matter fields such as the electromagnetic radiation. The energy momentum tensor for the dark energy is given by, where the first term on the right hand side corresponds to the homogeneous background Friedmann-Robertson-Walker (FRW) spacetime, where the indices of the Kronecker delta run over the spatial coordinates only, while the energy density and the pressure satisfy The second term on the right hand side of Eq. (22) corresponds to the inhomogeneous linear perturbation due to a central and spherical massive object, like a galaxy. Up to linear order the most general such perturbation should be written as where we have ignored second order terms, e.g. of order ∼ Ψδρ E . δT E ij is not proportional to δ ij and, in particular, it may contain off-diagonal terms. Likewise, we have the energy-momentum tensor corresponding to the cold dark matter where the first term on the right hand side corresponds to the background FRW spacetime and the second term represents the central mass M at rest with respect to the comoving frame. Since we are interested in determining the maximum size of a structure, the region of our interest is essentially outside it. Then, just outside the structure by virtue of the spherical symmetry of the problem, we may take all mass to be localized at R = 0. Finally, the third term in Eq. (25) corresponds to a backreactionless or test spherical dark matter fluid whose dynamics we wish to study, where v i = dx i /dt is the peculiar velocity of the test fluid element. Since that fluid under consideration is dark matter, it must be nonrelativistic i.e., |v i | ≪ 1 and hence we have ignored O(v 2 ) terms in the above ansatz. The negative sign in front of the second term is a mere convention -after all we may flip it to the positive sign by just changing the direction of motion of the fluid. We further assume that the weak energy condition (T tt > 0) is satisfied by all matter fields and their perturbations. This completes the specification of the sources and we shall now use them into the Einstein equations. We further refer our reader to [53,54] for some recent developments of the cosmological perturbation theory with a positive Λ, around a central point mass. However, before proceeding any further, we may simplify the Einstein tensors in Eq. (21) as follows. We note that, as we have repeatedly emphasized, we are interested here in a length scale much larger than the Schwarzschild radius and much smaller than the Hubble radius. In such scales, we may make a quasistatic approximation, in which the spatial derivatives of the potentials Ψ and Φ are much larger than their temporal derivatives. Also, if k i is a component of the spatial momenta associated with the Fourier transformed potentials, it is clear from the length scale we are interested in that |k i | ≫ H −1 (t). Putting these all in together and using Eq. (23)-Eq. (25), we obtain for the Einstein equations, for the FRW background and 1 a 2 (t) for the perturbation. It is possible to solve, at least formally if not analytically, the set of Eq. (27), Eq. (28) as follows. First, for a given dark energy model, we can solve Eq. (27) for the scale factor a(t). Coming next to the first of the above equations, integrating both sides over a constant time hypersurface, using the divergence theorem and the spherical symmetry of the problem, it is easy to obtain where δM E (R, t) is the contribution to the proper mass function due to the dark energy perturbation, Since we have already assumed that all matter fields and their perturbation satisfies the weak energy condition, we have δM E (R, t) > 0. Note that the mass function, M + δM E (R, t) could be thought of as an analogue of the Tolman-Oppenheimer-Volkoff mass function defined in a static scenario for a fluid [52], where the self gravitating fluid may now be imagined to have two components, both satisfying the positivity of the energy density.
Let us now come to the second of Eq. (28), in order to find the relation between the two potentials. Clearly Φ = Ψ only when T E ij = 0 for i = j [55]. Also, by the spherical symmetry of the problem, we have ∂ θ (Φ − Ψ) = 0 = ∂ φ (Φ − Ψ). Then, converting the indices (i, j) from Cartesian to spherical polar, we have δT E θθ = δT E φφ = δT E θφ = 0, and For the explicit form of the above equation in the case of Brans-Dicke cosmology we refer our reader to [36]. Using Eq. (30), we now obtain where we have set the integration constant to zero, as in the absence of any inhomogeneity, we must have Ψ = 1. We also have defined a total or effective mass function, Let us now come to the sector of the cold dark matter. The homogeneous part, ρ(t) satisfies the conservation equationρ (t) + 3Hρ(t) = 0 whereas the inhomogeneous backreactionless perturbation specified in Eq. (26) satisfies (the conservation equation ∇ a δT M ab = 0, in the background Eq. (18), The maximum turn around radius is defined as the point where the proper 'peculiar acceleration' (= d 2 (a(t)x i )/dt 2 ) of a fluid element vanishes which, after using the second of Eq. (34) becomes As discussed in the previous section, the largest of all the maximum turn around radii is obtained for zero orbital angular momentum. Thus, setting i = R above and using Eq. (27), Eq. (32), we obtain Writing δM := 4 3 πρ(t)R 3 a 3 (t) and Λ(t) := ρ E (t)/8πG, we find the proper maximum turn around radius where we have defined an effective mass function M eff = M tot +δM , taking into account the contribution of all sources. The above expression for the maximum turn around radius is formally similar to the result of [5] obtained for a constant equation of state parameter and without considering any dark energy perturbation, whose effects here have been dumped into the mass function M eff . Setting w(t) = −1 in Eq. (37) one recovers the ΛCDM result. Note also that in deriving the above expression, no assumption, whatsoever was made on the behaviour of w(t). Thus, from the analysis of [5] it follows that any dark energy model falling into the generic class we are considering here with w (t = t today ) −2.3 is consistent with the structure size vs mass observational data.
A few of the geometric aspects of the maximum turn around radius could be noted here. Firstly and obviously, the proper radial distance which is purely spatial, is an invariant owing to the spatial isotropy of the systems like Eq. (8), Eq. (18), [1]. In other words, for all observers obeying the isometries, R TA,max is an invariant quantity. Second, we may note here the purely geometric origin of the definition of the turn around radius in Sec. 2, from the vanishing of the acceleration of the orbits of the timelike Killing vector field. For a discussion on how such consideration in a static spacetime Eq. (8), could be linked to the general McVittie spacetime Eq. (18), we refer our reader to [36].
Before ending this section, we wish to show that Eq. (35) could in fact be related to the geodesic equation. Indeed, let us first rewrite the metric Eq. (18) in terms of the proper radial coordinate Since we are much outside the Schwarzschild radius of the structure and much inside the Hubble radius (= H −1 (t)), the metric becomes to linear order in Ψ, Φ and H(t)r where the prime denotes differentiation once with respect to the proper radial coordinate. Let us consider a test particle moving along a timelike geodesic in the above geometry, where λ is an affine parameter and we consider a radial, nonrelativistic trajectory as earlier. In this limit, λ ≈ t and hence dt/dλ ∼ 1. Thus the conservation of norm of the geodesic implies : dr/dλ ≪ 1. Then we have, Using Eq. (39) and its inverse and keeping in mind that the radial derivative of the potential dominates over the temporal one at the length scales we are interested in, we obtain Inserting into this the vanishing acceleration condition, d 2 r/dt 2 = 0, Eq. (35) is reproduced. Clearly, such equivalence between the perfect fluid and geodesic is possible due to the fact that to derive Eq. (35), we had ignored all higher order terms coming from the energy-momentum tensor of the test fluid.

Quintessence models
As an application of the above general formalism, the first example we take is the quintessence model, originally proposed in [56] to address the coincidence problem of the dark energy (see also [57] for a review). The role of the dark energy is played here by a time dependent scalar field φ 0 (t). The model is described by the action where V (φ) is some positive potential, to be chosen appropriately in order to generate the accelerated expansion and L M is the Lagrangian density corresponding to any other matter field, such as the cold dark matter. Assuming a homogeneous FRW background solution with φ = φ 0 (t), its energy density and pressure are given by respectively. There could be various choices of the potential V (φ). For example, the so called tracker model, where the parameter M has dimensions of mass and p is a positive real number. Or, the so called thawing model : with the parameters M and µ having dimensions of mass. For a review on different such choices for the potential and their cosmological implications including the perspective of particle physics, we refer our reader to [57] and references therein. For any potential V (φ) the homogeneous background is a fluid of the type studied in the previous section, with state 'parameter' given by Thus, it is clear that in order to have w(t) < 0, the potential must dominate over the kinetic term. Furthermore, for positive potential the state parameter satisfies w(t) ≥ −1. Then, according to the discussion of this section, one may immediately conclude that the generic model of this class predicts cosmic structure sizes greater than ΛCDM and, thus, is consistent with the structure size vs mass observations.

The generalized Chaplygin gas model
The other theory we wish to address here is the generalized Chaplygin gas model. The Chaplygin gas or its generalized version attempts to give a unified picture of the dark energy and the dark matter [58,59]. Let us consider the equation of state where both β and A are positive constants. The conservation equation :ρ + 3H(ρ + P ) = 0 gives where ρ m is a positive integration constant. At sufficiently late times, a(t) ≫ 1 (the current scale factor is set to unity), we have ρ ≈ A 1 β+1 and hence P ≈ −A 1 β+1 , giving the equation of state of a cosmological constant. On the other hand at sufficiently early times, a(t) ≪ 1 and we have ρ ∼ a −3 (t) and |P | ≪ ρ, i.e. the matter dominated universe.
The original Chaplygin gas model was defined specifically for β = 1. Then, it turns out that the dust dominated universe is smoothly interpolated with the dark energy dominated universe via an intermediate equation of state, P = ρ, known as the stiff matter. Such an intermediate equation of state is necessary in order to ensure that the normal matter has positive pressure. The generalized Chaplygin gas model, on the other hand, permits 0 < β ≤ 1, which we shall consider below. In this case the matter and dark energy dominated universe is interpolated via an equation of state P = βρ with β > 0 being a parameter of the theory and β = 0 corresponding to the ΛCDM.
We are interested in the current universe dominated by the dark energy. So, we may expand Eq. (47) and Eq. (48) to obtain at leading order Thus, the above two equations represent the dynamics of our universe with an effective homogeneous fluid with equation of state P (t) = −A 1 1+β + β(ρ(t) − A 1 1+β ). Solving Eqs. (Eq. (27)) for the scale factor we obtain with H 0 = 8πGA 1/(1+β) /3. For t → ∞, we recover the de Sitter space, while for β = 0 we recover the ΛCDM. In order for this model to be consistent with modern cosmology, the value of the parameter β needs to be rather small, β 0.05 [60]. Eq. (49) show that we do not need to consider the inhomogeneous dark energy perturbation in this case, and hence we have Ψ = Φ in Eq. (18). With this, we study the dynamics of a test Generalized Chaplygin gas fluid, following the method described in the preceding section. Since the parameter β is small, we take terms like βδρ or βv i to be second order only, giving the maximum turn around condition, a a R − 1 + β a 2 ∂ R Ψ = 0 The Chaplygin gas equation of motion gives, Writing 8πGA 1/(1+β) = Λ and following the procedure described in the previous section, it is straightforward to obtain to first order in the parameter β, where δM := 4 3 πR 3 a 3 (t)ρ(t) as earlier. Setting β = 0 one recovers the ΛCDM result. Thus, since β > 0 the prediction of the generalized Chaplygin gas model for the maximum turnaround radius is larger than in ΛCDM. This fact can be understood as the effect of the positive pressure fluid, that opposes the accelerated expansion and implies that the generalized Chaplygin gaz too is consistent with the data on the sizes and the masses of the cosmic structures.

Conclusions
In this work we have discussed various aspects of the maximum possible structure size, in the context of the mass versus observed sizes of some nearby large scale cosmic structures. In Sec. 2, we have given a derivation of the maximum turn around radius for a dark energy with constant equation of state parameter (P = wρ with w const.) in a purely static geometry derived at the maximum turn around length scale (Eq. (16)). This result was found earlier in [5] using cosmological perturbation theory. Even though for w = −1, the dark energy is usually taken to be time dependent and spatially homogeneous and isotropic, we have shown here that the same expression for the R TA,max holds also when we consider the dark energy to be inhomogeneous and time independent.
In Sec. 3 we have derived the expression for the R TA,max in a cosmological scenario, using the McVittie geometry in the weak field limit. The dark energy we have taken here is generic and has completely arbitrary, time dependent state 'parameter' (P (t) = w(t)ρ(t) with w(t) arbitrary). We also have taken the first order backreaction effect on the dark energy due to the central inhomogeneity. The final result obtained in Eq. (37) is formally exactly similar to that of the case of w(t) = const. and hence giving the same bound on the numerical value of the state parameter, w (t = t today ) > −2.3. We have applied this general formalism in Sec. 3.2 to two quite popular alternative dark energy/gravity models viz. the quintessence and the generalized Chaplygin gas and have shown that both these models pass the size versus mass test with flying colours.
The framework we have considered, even though it goes beyond the ΛCDM, essentially is described by the Einstein equations. It should thus be highly interesting to extend such general framework to go beyond the usual 4-dimensional Einstein equations and to address more complicated scenarios such as the brane-world [61]. We hope to return to this issue in the near future.