Cosmological tests of the gravastar hypothesis
Abstract
Gravitational vacuum stars (gravastars) have become a viable theoretical alternative to the black hole (BH) endstage of stellar evolution. These objects gravitate in vacuum like BHs, yet have no event horizons. In this paper, we present tests of the gravastar hypothesis within flat Friedmann cosmology. Such tests are complementary to optical and gravitational wave merger signatures, which have uncertainties dominated by the poorly constrained gravastar crust. We motivate our analysis directly from the action principle, and show that a population of gravastars must induce a timedependent dark energy (DE). The possibility of such a contribution has been overlooked due to a subtlety in the de facto definition of the isotropic and homogeneous stress. Conservation of stressenergy directly relates the time evolution of this DE contribution to the measured BH comoving mass function and a gravastar population crust parameter . We show that a population of gravastars formed between redshift readily produces the presentday DE density over a large range of , providing a natural resolution to the coincidence problem. We compute an effective DE equation of state as a function of presentepoch BH population observables and . Using a BH population model developed from the cosmic star formation history, we obtain and consistent with Planck bestfit values. In summary, the gravastar hypothesis leads to an unexpected correlation between the BH population and the magnitude and timeevolution of DE.
0000000269170214]K. A. S. Croker
1 Introduction
The recent direct observation of gravitational radiation from ultracompact, massive, object mergers provides definitive proof of the existence of objects with exterior geometries consistent with the classical black hole (BH) solutions. Unfortunately, the various boundaries and interiors of the classical BH geometries are severely pathological. In fact, the predicted physical curvature singularities and closed timelike curves had already motivated the exploration of BH “mimickers.” Such solutions to classical GR appear as classical BH geometries, as perceived by exterior vacuum observers, but replace the interior with another GR solution. Gliner (1966) is usually recognized as the first researcher to suggest that “ vacuum,” localized regions of vacuum with large cosmological constant, could be relevant during gravitational collapse.
A simplified scenario, which motivates more physically plausible constructions, is the following exact, spherically symmetric, GR solution
(1) 
with
(2) 
This solution is a Schwarzchild BH exterior for , and a vacuum with cosmological constant interior to . This solution is static, consistent with the intuition that material under tension stays bound. Note that the interior energy density correctly integrates to . This means that a vacuum observer at infinity will perceive a point mass . The interior region contains no physical singularity, but the infinite pressure gradient at in Eqn. (1) suggests consideration of slightly more sophisticated GR solutions.
Dymnikova (1992, Eqn. (14)) produced a “Glump” model without an infinite pressure gradient. About a decade later, Mazur & Mottola (2004) and Chapline (2003) proposed additional realizations of what is now called a gravitational vacuum star (gravastar). The essential feature of viable gravastars is a “crust” region of finite thickness slightly beyond , which surrounds the de Sitter core. These objects have no trapped surface: each interior is in causal contact with the outside universe. In fact, this ability to resolve the BH information paradox sparked renewed interest in BH mimickers during the “firewall” debates initiated by Almheiri et al. (2013). Theoretical checks on gravastar stability to perturbations (e.g. Visser & Wiltshire, 2004; Lobo, 2006; DeBenedictis et al., 2006) and stability to rotation (e.g. Chirenti & Rezzolla, 2008; Uchikata & Yoshida, 2015; Maggio et al., 2017) have found large ranges of viable parameters describing the crust. Phenomenological checks from distinct merger ringdown signatures, such as an altered quasinormal mode spectrum (e.g. Chirenti & Rezzolla, 2007) have been computed. In addition, possible optical signatures from the crust itself (e.g. Broderick & Narayan, 2007), accretion disks (e.g. Harko et al., 2009), and even direct lensing (Sakai et al., 2014) have been investigated.
The first direct detection of gravitational radiation from massive compact object mergers has greatly stimulated critical analysis of the gravastar scenario. Bhagwat et al. (2016, 2017) show that the existing aLIGO interferometer and the proposed Cosmic Explorer and Einstein telescopes can resolve the higher harmonics, which could probe the existence of gravastars. On the other hand, Chirenti & Rezzolla (2016) have already claimed, using the ringdown, that GW150914 likely did not result in a gravastar final state. Yunes et al. (2016) criticize the Chirenti & Rezzolla (2016) claim as premature, yet argue that ringdown decay faster than the light crossing time of the remnant poses the more severe challenge.
In all phenomenology to date, however, gravastar signatures depend critically on the unknown properties of the crust. In this paper, we address this problem and develop complementary observational signatures dominated instead by the known properties of the core. Before proceeding, we first comment on consistency with Birkhoff’s theorem, which gravitationally decouples the interior properties of a localized object from its surroundings. A necessary condition for Birkhoff’s theorem is vacuum boundary conditions. Thus, we will consider gravastars within Friedmann cosmology, which is nowhere vacuum.
Consistent with the theoretical motivation to resolve the BH information paradox, we will replace all BHs with gravastars. The following falsifiable experimental prediction for the physical dark energy (DE) density
(3)  
(4) 
will then be obtained. Here is a small parameter which measures the deviation of the gravastar contribution from pure de Sitter due to crusts, is a cutoff before which there are no gravastars, and is the comoving BH density, which has just begun to be constrained by aLIGO. As detailed by Dwyer et al. (2015), the proposed Cosmic Explorer interferometer will be able to definitively constrain with events per year, and out to redshift . The schematic gravastar given in Eqn. (1) has no crust, so and the prediction becomes entirely parameterfree
(5) 
The rest of this paper is organized as follows. In §2, we motivate a departure from the de facto Friedmann source (DFFS). We then define notation that simplifies our calculations and clarifies the physical observables. In §3, we construct the appropriate cosmological source from an assumed population of gravastars. In §3.1, we produce the fundamental quantitative prediction relating the DE density to the observed BH population and the gravastar crust parameter . In §4, we proceed chronologically through the history of the Universe and consider the influence of a gravastar population at three different epochs. In §4.1, we consider the formation of primordial gravastars and essentially exclude their production. In §4.2, we demonstrate that stellar processes between redshift can produce enough gravastars to account for the presentday DE density. In §4.3, we analyze the gravastarinduced DE at late times () in the dark fluid framework described in Ade et al. (2016). Here, we produce predictions for the DE equation of state parameters and in terms of presentday BH population observables.
A comprehensive set of appendices presents material complementary to the phenomenological results of this paper. Appendix A highlights nontrivial assumptions fundamental to the DFFS. Appendix B clarifies the spatiallyaveraged nature of the Friedmann source by deriving Friedmann’s equations directly from the principle of stationary action. Appendix C establishes that the positive pressure contributions of typical astrophysical systems can be ignored. Appendix D discusses how causality is maintained, given a nontrivial time dependence in the background cosmological source. Appendix E develops a simple model for the BH population in terms of the stellar population. Appendix F gives an upper bound for the pressures of large, virialized systems like clusters.
Except within the appendices, we take the time unit as the reciprocal presentday Hubble constant , the density unit as the critical density today , and set the speed of light . For efficiency, we will often describe quantities with scale factor or redshift dependence as timedependent. We will also freely use either scale factor or redshift , depending on which produces the more compact quantitative expression. Though we replace all BHs with gravastars, for consistency with existing astrophysical literature, we will often refer to these objects as BHs.
2 Isotropic and homogeneous stress in Friedmann cosmology
Given the symmetries of the RW metric, RW observers cannot distinguish points: in a homogeneous and isotropic universe, every point is observationally identical. In other words, a RW universe does not contain any notion of interior or exterior. Perhaps surprisingly, the de facto Friedmann source (DFFS) implicitly violates homogeneity and isotropy. As detailed in Appendix A, the DFFS assumes the ability to define a notion of interior and exterior, distinct for every mass distribution. This assumption is used to reduce spatially extended systems to effective point masses. These point masses are then considered to contribute only as a pressurefree “dust.”
Consistency, however, requires that the source to Einstein’s equations exhibit the same symmetries as the metric. This requirement holds at each order in a perturbative expansion of Einstein’s equations. In particular, the spatially uniform background source in Friedmann cosmology cannot contain any implicit notion of interior or exterior. In Appendix B, we construct a consistent Friedmann source (CFS) by application of the principle of stationary action to an inhomogeneous fluid ansatz. The resulting zeroorder source simply becomes a flat spatial average.
In contrast to the DFFS, the spatial average present within the CFS necessarily incorporates all pressures: cluster interiors, stellar interiors and, should they exist, gravastar interiors. In Appendix C, we show that accounting for positive astrophysical pressure contributions makes no observable changes. The effect of any negative pressure contributions, however, can critically influence the background expansion. An averaged gravastar interior contribution is both negative and strong, with . This feature of gravastars ultimately leads to the unexpected and quantitative relation between the expansion rate and the BH population.
2.1 Twocomponent approximation
Let barred variables represent physical, as opposed to comoving, background quantities. We investigate the CFS with two contributions. The first contribution, the contribution, is an approximation to the influence of gravitationally bound systems with pressure
(6) 
It is an approximation because we do not attempt to track all pressures within the universe. For the primary discussion of this work, the contribution will be exclusively from gravastars. For the discussion in Appendix C, justifying why this is an excellent approximation, the contribution will be a severe upperbound on positive pressure contributions from typical astrophysical systems.
As is well known, the constituents of tightly bound gravitational systems decouple from the expansion. The number density of these systems then dilutes with the physical volume. This behavior motivates the following definitions
(7)  
(8) 
allowing us to to rewrite Eqn. (6) in an entirely equivalent way
(9) 
The explicit appearance of factors and timedependent comoving quantities greatly simplifies many computations and will help to emphasize the difference between the DFFS and the CFS.
The second contribution, the contribution, comes from collisionless matter
(10) 
and is identical to the analogous term of the DFFS whenever . This derivative statement is an example of the utility of timedependent comoving densities, when highlighting the differences between the DFFS and the CFS.
In the gravastar scenario, effectively collisionless matter is processed into localized regions of de Sitter space. We thus populate by depleting an appropriate fraction of the collisionless contribution
(11) 
Here is the observed matter fraction today. The depletion fraction is determined through astrophysical observation.
The CFS does not conserve particle number by construction, but this is not required by Einstein’s equations (Weinberg, 1972, §2.10). Covariant conservation of energy and momentum
(12) 
however, necessarily populates both and in a manner consistent with Einstein’s equations. Thus, timedependence in comoving , , and does not violate thermodynamics or energy conservation. For a discussion of causality in the context of a nontrivially timedependent background, we direct the reader to Appendix D.
3 Gravastar contribution to the Friedmann equations
A dynamically stable gravastar, with physically plausible structure, is more sophisticated than Eqns. (1) and (2). Both Dymnikova (1992) and Mazur & Mottola (2004) introduce a transition region, which ultimately encloses a dark energy interior. The possible nature of this “crust” has been studied extensively by Visser & Wiltshire (2004); MartinMoruno et al. (2012) and others. There is consensus that substantial freedom exists in the construction of gravastar models. This freedom, while constrained, permits a range of crust thicknesses and equations of state. We will consider the following structure

: de Sitter interior

: crust region

: vacuum exterior
where is a gravastar mass and is a crust thickness. Note that we have simplified the model by placing the inner radius at the Schwarzchild radius. In any physical gravastar, this inner radius must be slightly beyond the Schwarzchild radius: this prevents the formation of a trapped surface. For our purposes, this nuance will not matter. The cosmological contribution from a plausible gravastar population thus consists of two components
(13)  
(14) 
Note that , , and are not independent degrees of freedom because each gravastar has a modeldependent internal structure.
We may define an equation of state for the gravastar source
(15) 
Mazur & Mottola (2004) assert that most of the energy density is expected to reside within each core and that within each crust. It is thus reasonable to expand the aggregate quantities and discard higherorder terms
(16)  
(17) 
This expression successfully isolates ignorance of the crust to a small dimensionless parameter, which we now define as
(18) 
To avoid confusion, we emphasize now that is not the Dark Energy equation of state as constrained by Planck. The relation between these two quantities will be explored thoroughly in §4.3.
3.1 Relation of to the BH population
Let be the comoving coordinate BH mass density. The formation of gravastars, and any subsequent accretion, depletes the baryon population. This decreases by as per Eqn (11). Consistent with the results of Appendix C, to very good approximation, we may regard all (nongravastar) positivepressure contributions as collisionless within the Friedmann source. Thus, the timevariation of comes entirely from the timevariation of . The conservation Eqn. (12) then becomes
(19) 
Switching from coordinate time to scale factor gives
(20) 
This equation can be separated with an integrating factor
(21)  
(22) 
where we have defined the proportionality to be unity and is a cutoff below which there are no gravastars. If we omit the cutoff parameter for , then we implicitly assume . The resulting energy density is
(23) 
where the lowerbound does not depend on because vanishes by definition before . When convenient, we will study the phenomenology of gravastars with the single parameter held fixed. This is consistent with other phenomenological studies of gravastars. In practice, we will find that timevariation of does not affect the essential cosmological signatures of a gravastar population.
3.1.1 Physical interpretation of DE
To develop some intuition for Eqn. (23), consider a sequence of instantaneous conversions
(24) 
where is the comoving density of baryons instantaneously converted at . Substitution into Eqn. (23) gives
(25) 
Explicitly, at the instant of the th conversion , we have the following comoving density of DE
(26) 
Thus, at the instant of conversion, agrees with a spatial average over the newly formed gravastars’ interiors. In other words, a quantity of baryonic mass has been converted to an equal quantity of DE. Beyond the instant of conversion, however, the contribution due to gravastars dilutes more slowly than the physical volume expansion. For example, in the de Sitter limit of , we find that and
(27) 
In this case, the physical baryon density becomes “frozen in” at the time of conversion. To maintain locality in general, one might conclude that each gravastars’ local mass must increase such that
(28) 
but this would violate the Strong Equivalence Principle (e.g. Will, 1993). Though such behavior is conceivable within scalartensor theories, we will remain focused on gravastars within GR.
Evidently, we cannot make any quantitative statements about the spatial location of the subsequent energy averaged to produce . This is by choice, as our analysis is designed to minimize the impact of the unknown local details that produce the gravastar. Since we have accomplished this through covariant conservation of stressenergy, we can proceed with confidence that our predictions are entirely consistent within GR. Strictly speaking, however, Eqn. (17) must be interpreted as a motivation for , and nothing more.
4 Observational implications
To begin constraining the gravastar hypothesis, we determine when gravastar production is observationally viable. Our discussion will proceed chronologically through the history of the Universe. This will inform our subsequent analysis: at late times will depend on an orderofmagnitude understanding of its earlytime behavior.
4.1 Constraint on primordial pBH as a function of temperature
The usual constraint on the fraction of critical density capable of collapse into a primordial BH is expressed in terms of the mass of the target hole: (e.g. Carr, 2003). The mass of the hole is related to via
(29) 
where is the radiation density parameter today. Thus, we may construct an analogous constraint , where is some temperature during the radiationdominated epoch. We have switched to temperature using , where K is the CMB temperature today.
Assume that all primordial BHs are produced at a single moment in time
(30) 
Here is the comoving pBH density at . For simplicity, consider fixed in time. In this setting, Eqn. (21) becomes
(31) 
Inserting this and Eqn. (30) into Eqn. (23) we find
(32) 
and evaluating at the present day gives
(33) 
The presentday dark energy density gives a conservative upper bound because we are neglecting accretion. According to Carr (2003, Eqn. 7),
(34) 
which we can relate to at the epoch of equality . With these relations, Eqn. (33) becomes
(35) 
Rearranging the previous expression for we find
(36) 
which is presented in Figure 1.
Of course we do not claim that detection of even a single gravastar could enforce the constraints shown in Figure 1. The essential point, however, is that primordial production of gravastars is heavily disfavored. For small , the amplification in energy density from the term is extreme at primordial epochs.
4.2 Gravastars as the only source of dark energy
The prediction given in Eqn. (23) tightly correlates the dark energy density to the matter density. This suggests a natural resolution to the coincidence problem: gravastars are responsible for all of the cosmological dark energy. This possibility is plausible because the in Eqn. (23) strongly amplifies the effect of conversion at early times. Unfortunately, until gravitational wave observatories begin to constrain , a precise discussion must remain focused on the latetime changes in . To frame these changes in an appropriate context though, it is useful to understand the general features of any gravastar formation history that could suffice to resolve the coincidence problem. Since this is a very important consequence of a gravastar population, we will approach the question in two complementary ways.
4.2.1 Technique I: Instantaneous formation epoch from presentday BH density
First, we will determine an instantaneous formation epoch, given the presentday BH density. In other words, we will determine an instantaneous formation time for all BH density such that the correct is obtained. When gravitational wave observatories have constrained this density directly, this estimate can be checked anew. Until then, we will use the stellar population based BH model developed in Appendix E.
To proceed with an order of magnitude estimate, we again assume that all gravastars are produced at a single moment in time
(37) 
We take the total cosmological density in gravastars to be the presentday values
(38) 
Here, rapid and delayed refer to the stellar collapse model as detailed by Fryer et al. (2012). These values are likely slight underestimates due to accretion, which is neglected in Appendix E. Note that these values are consistent with the fraction of stars that collapse to BH, as estimated by Brown & Bethe (1994). Substituting the approximation Eqn. (37) into Eqn. (23) we find
(39) 
where we have again fixed , and used that for any fixed . The instantaneous formation epoch is displayed in Figure 2.
4.2.2 Technique II: Constant formation from onset, rate consistent with ray opacity
Parameter  Value (natural units)  Reference 

Inoue et al. (2014, Fig. 2)  
Inoue et al. (2014, Fig. 2)  
Abel et al. (2002)  
Table 5 
Note. – data are used to constrain an approximate gravastar formation epoch, which suffices to produce all presentday dark energy density by some cutoff redshift . Variables defined in the text. The results of this analysis can be found in Figure 3.
The stellar model used in §4.2.1 may suffer from systematics associated with interpolating the collapse remnant distributions between only two metallicities. To estimate the effect of this systematic, we will use unrelated astrophysical constraints determined from ray opacity reported by Inoue et al. (2014). We will assume that a constant fraction of stellar density collapses to BH between stellar onset up until some cutoff time . We will determine the value of such that the produced BH density yields the observed today. We will further assume that the stellar density formation rate during is constant.
Since we are considering formation over an extended period of time, for simplicity, we continue to restrict our attention to fixed in time. Under these assumptions, we find that
(40) 
Here we have introduced a constant collapse fraction and a factor to account for accretion processes. We will take . This is reasonable, given that Li et al. (2007) have shown that accretion can be very efficient at primordial times. With the constant star formation rate assumption, we may use the chain rule to write
(41) 
We present reasonable data for the above parameters in Table 1, where the formation rates listed are upper limits. Given the matterdominated Hubble factor
(42) 
we may integrate Eqn. (41) and solve for the cutoff time. The result is
(43) 
The behavior of Eqn. (43), converted to redshift, is shown in Figure 3.
Note that we have taken the more conservative bound established at . If we take , there is no viable space. This is consistent with the stellar model of Appendix E, where the primordial collapse fraction is . It is not surprising that the collapse remnant model’s primordial is larger than the value of Brown & Bethe (1994), because their analysis depends on iron cores, which are not present in Population III stars. Note that the fraction of matter density converted to gravastar material under our assumptions,
(44) 
is consistent with the presentday used in §4.2.1.
4.2.3 Discussion
In this section, we have considered whether gravastars alone could account for the presentday observed DE density. This is plausible because the DE induced by gravastar formation is amplified by . Using two complementary techniques, we have estimated an approximate formation epoch for gravastars, which suffices to produce the presentday observed . Encouragingly, both techniques agree that sufficient production is viable over a large range of . This resolves the coincidence problem (e.g. Amendola & Tsujikawa, 2010, §6.4) of a narrow permissible value for . Furthermore, both techniques converge toward and produce excluded regions for gravastar
(45)  
(46) 
Above these thresholds, resolution of the coincidence problem by gravastars is disfavored. As we will soon seen, this exclusion is consistent with those based on latetime behavior and existing Planck constraint.
This result is very useful for many reasons. Most importantly, the viable region lies squarely in the middle of the epoch of primordial star formation. This suggests that a naturally emerging population of gravastars can produce the correct DE density. In addition, the density converted to induce this effect is of . This is well within uncertainties on the Planck bestfit value for . Since , production of a sufficient population can occur without breaking agreement with precision CMB astronomy. In other words, CMB anisotropies are not altered at any level. Finally, gravastar formation at can establish a presentday valued DE density after star formation has begun. Thus, initial conditions, and therefore results, of precision body simulations are not altered.
In any scenario where , the produced physical DE density will diminish with time . Relative to the required present day value of then, the physical DE density is larger in the past. This behavior is not relevant for the present discussion, as it represents a minor increase in any physical DE density due to gravastars with viable . During the dark ages, the matter density (dominated by dark matter) completely overwhelms any DE. This effect, however, is tracked in §4.3.
For the subsequent discussion, we interpret these estimates in the following way: if gravastars are to resolve the coincidence problem, we expect that most of the dark energy density will be established by stellar and accretion physics taking place during .
4.3 Gravastars analyzed in the dark fluid framework
Dark Energy is most commonly studied through constraint of three parameters: an energy density , and the linear Taylor expansion, about , of its equation of state parameter . This analysis of possible deviations from CDM comes from the anticipated dynamics of a minimally coupled scalar field: the dynamics of an internally conserved fluid.
In the gravastar scenario, stellar collapse rates, subsequent accretion, and evolution in the stellar remnant distribution all induce changes in the DE density . We can, however, choose to interpret this timedependence of as the dynamics of an internally conserved fluid. This leads to an effective equation of state parameter as follows. Define through
(47) 
Subtracting Eqn. (19) from Eqn. (47), switching to scale factor, substituting Eqn. (23), and solving for gives
(48) 
where prime denotes derivative with respect to . Using the definition of in Eqn. (21) and linearity of the integral, we may reexpress this relative to the presentday value of
(49) 
This form is advantageous as it samples only values very near to the present epoch, where constraint and systematics are wellcharacterized. The firstorder Taylor expansion in of Eqn. (49) does not contain at all
(50)  
(51) 
The coefficient Eqns. (50) and (51), combined with astrophysically measured values and uncertainties in BH population parameters, will give a region which can be immediately compared against any Planckstyle constraint diagram.
The results of §4.2 suggest that most of should be already produced by . In this case, we may approximate Eqn. (49) as
(52) 
We have shown that knowledge of the black hole mass function and the presentday dark energy density completely determines . In this way, gravitational wave observatory measurements of the BH population make definitive cosmological predictions. We already have a wealth of cosmological observations, however, with upcoming experiments (e.g. Levi et al., 2013, DESI) set to constrain . We may thus make predictions for gravitational wave observatories by inverting the above procedure. Solving for by differentiating Eqn. (49) and reintegrating, we find
(53) 
Note that we have removed entirely through use of Eqn. (21). Given the particularly simple assumed by Planck
(54) 
we may express Eqn. (53) in closed form
(55) 
4.3.1 Estimation of from the comoving stellar density
In the gravastar scenario, the physical origin of DE is completely different than the internal dynamics of a scalar field. The linear model Eqn. (54), based on the dark fluid assumption, may not be the most suitable for characterization of a latetime gravastar DE contribution. In this section, we use a BH population model, developed from the stellar population in Appendix E, to investigate in the gravastar scenario. These results will allow us to estimate when Eqns. (54) and (55) wellapproximate the timeevolution of the comoving BH density.
The BH population model of Appendix E does not account for accretion, so we must remain outside of any epoch where BH accretion is significant. Consistent with the simulations of Li et al. (2007), we consider only . To begin evaluation of Eqn. (49), we require , which requires . For fixed , the required integrations become trivial. We find
(56) 
and predict for Eqn. (52)
(57) 
This relation is displayed in Figure 4. The range corresponds to “rapid” and “delayed” models of stellar collapse, as determined by Fryer et al. (2012). The collapse model determines the distribution of remnant masses, which evolves in time with metallicity. Given , Eqn. (57) provides a single parameter fit, with subpercent precision, viable for . It is clear that the Planck linear ansatz is only useful for
(58) 
with departures growing significantly worse as .
4.3.2 Planck and DES constraint of at the present epoch
Physically, must be satisfied wherever ansatz Eqn. (54) is valid. From Eqn. (53), this constraint becomes an upper bound
(59) 
Regardless of any timeevolution in , we may use Eqns. (50) and (51) to constrain now. Eliminating from these equations, we find a permissible line through space
(60) 
We will use the BH population model, developed from the stellar population in Appendix E, to estimate the derivatives of . From Eqn. (E17), we numerically find that
(61)  
(62) 
which will produce a permissible band in space. For , this band is displayed in Figures 5 and 6. Evidently, Planck data disfavor gravastars with large
(63) 
Within this range, however, gravastars are consistent with Planck best fit constraints.
By inspection of Eqn. (60), simply translates the constraint region vertically. If we are to remain consistent with for , we must have
(64) 
This bound, and the positivity bounds given in Eqn. (59) are displayed in Figure 5 for a variety of .
Recent results constraining the CDM model from the Dark Energy Survey (DES) (Abbott et al., 2017) cannot be immediately applied to constrain . This is because induced by a gravastar population changes in time even if remains fixed for all time. Since we predict only small changes in , however, it is reasonable to expect that the CDM model will approximate the gravastar scenario at late times. Indeed, their reported value of (Abbott et al., 2017, Eqn. VII.5) becomes the following constraint on
(65) 
which is consistent with the constraints reported in §4.2.3.
5 Conclusion
Gravitational vacuum stars have become a viable and popular theoretical alternative to the pathological classical black hole (BH) solutions of GR. These objects appear as BHs to exterior vacuum observers, but contain de Sitter interiors beneath a thin crust. This crust is located above the classical Schwarzchild horizon, placing the entire gravastar in causal contact with the exterior universe. Many existing studies have focused on the observational consequences of the crust region for optical and gravitational wave signatures. These studies are hindered by systematics relating to the crust, which is only loosely constrained theoretically. In this paper, we have developed complementary constraints on the gravastar scenario based on the wellestablished properties of their de Sitter interiors.
We place our gravastars within a flat Friedmann cosmology, which is nowhere vacuum, and show from the action principle that the zeroorder Friedmann source must contain an averaged term sensitive to the interiors. This is consistent with Birkhoff’s theorem, which does not apply without vacuum boundaries. Through conservation of stressenergy, this term induces a timedependent Dark Energy (DE) density. This density is directly correlated to the evolution of nonlinear structure via star formation and subsequent collapse. The gravastar crust produces a small deviation from a pure de Sitter () equation of state. This deviation becomes the single parameter characterizing the gravastar population.
We replace all black holes with gravastars and consider the cosmological effects of their subsequent DE contribution at three epochs. During the primordial epoch (), we find that the fraction of matter collapsing into a primordial gravastar population with is constrained between orders of magnitude more than existing primordial BH constraints. During the dark ages we show, via two approaches, that existing astrophysical data support formation of a population of gravastars that can account for all of the presentday DE density. We show that any gravastar population with crust parameter can resolve the coincidence problem. During latetimes (), we precisely interpret the gravastar scenario in the usual language of a timevarying dark fluid. Using a BH population model built from the cosmic star formation history and stellar collapse simulations, we predict timevariation in the magnitude of that tracks star formation. We demonstrate complete consistency with Planck, given a gravastar population with . Further, we predict very little timevariation in at late times, consistent with the recent results of the Dark Energy Survey.
We make definitive predictions for both the gravitational astronomy community and dark energy surveys in the form of unexpected quantitative correlations between the timeevolution of the DE density and the BH population. In summary, the cosmological consequences of a gravastar population are unambiguous, readily testable, and already resolve many outstanding observational questions, without requiring any ad hoc departure from GR.
All code for generating the presented data and its visualizations is released publicly (oprcp4).
oprcp4, scipy (Oliphant, 2007), GNU Maxima, gnuplot
Appendix A Assumptions of the de facto Friedmann source at latetimes
In this appendix, we highlight nontrivial assumptions that enter the de facto Friedmann source (DFFS) during matter domination. In the DFFS, one imposes an additional hypothesis on the physical background matter density at latetimes (Peebles, 1980, §9, §97)
(A1) 
To understand its origin, consider the conservation of stressenergy statement for a single component perfect fluid
(A2) 
Here is the Hubble parameter (defined in §C.1) and is the background physical pressure. Equation (A1) then follows if one defines . Note that this is an assumption beyond Einstein’s equations with the RW metric. The origin of this assumption is microphysical. For , mean free paths of all particles (even hypothetical dark matter) are very large; particles are not scattering into or out of comoving volumes. This fixes the comoving number density of particles. Since the rest masses of elementary particles do not change, it follows that must be zero.
Given the lack of a formal mathematical definition for the DFFS, we consider the operational definition of the perturbations to , given by Peebles (1980, §81). He states that “we can imagine that observers spread through the universe and moving with the matter keep a record of the local density as a function of proper time, . As the observers come within the horizon, their records can be acquired and compared.” We now paraphrase this operational procedure:

Each observer measures and broadcasts (e.g. via light signals with presumably the same frequency) their own local density .

Each observer receives the others’ broadcasts and averages them with his own to produce a background .

Each observer computes his density perturbation as .
The same procedure may be performed with the local pressure , which will be nonzero in general. Between , averaged pressures remain negligible compared to averaged energy densities. Below , however, collapsed structures begin to form.
In order to justify continued use of the homogeneous and isotropic fluid ansatz for , one redefines Peebles’ observers. One regards an observer as reporting on a very large volume containing many gravitationally bound systems. Next, one invokes any number of formal embedding solutions to GR (e.g. Einstein & Straus, 1945), which glue a Schwarzchild metric into a RW metric. The existence of such formal solutions motivates defining a spatially isolated system to contribute only an active gravitational mass to the averaged energy density . In this setting, the only pressures that could contribute to would be from large structures, like rich clusters. As we show in Appendix F, however, upper bound pressures for these systems are . Thus, it may seem completely justifiable to continue to regard .
It must be emphasized, however, that formal embedding solutions are logically distinct from the perturbative treatment. The perturbative treatment assumes instead the exact RW background. Suppose we were unable to redefine our observers, and consider a Peebles’ type observer within the core of a gravastar. Here the pressure is equal in magnitude, but opposite in sign, to the energy density. The DFFS discards these contributions. In other words, the DFFS defines for because this is consistent with Newtonian intuition and static strong gravity in asymptotically flat space.
Appendix B The spatially averaged nature of the Friedmann source
In this Appendix, we rederive Friedmann’s equations directly from the EinsteinHilbert action. The action is a manifestly global quantity, from which local equations of motion are constructed. Since the action is an integral statement over the spacetime manifold, it will allow us to cleanly reconcile the isotropy and homogeneity of the RW metric ansatz with the anisotropic distribution of actual stressenergy.
Our approach will be to rewrite Friedmann’s isotropic and homogeneous model as a scalar theory in flatspace. The result will be an unambiguous interpretation of the Friedmann source as the flatspace, comoving coordinate average of the trace of the stresstensor. This wellmotivated notion of “locality” within the Friedmann model is lost if one begins with Einstein’s equations under the RW ansatz, simply because one must invoke some ad hoc procedure (e.g. Appendix A) to produce an isotropic and homogeneous source from the actual stressenergy.
The principle of stationary action requires that the total action be stable, at first order, to variations in the dynamical degrees of freedom
(B1) 
where represents the matter action and represents the gravitational action. We will consider the canonical definition of the stressenergy tensor given by Weinberg (1972, §12.2)
(B2) 
Here denotes the smooth manifold associated with the metric (O’neill, 1983, §3.2), is the variation of the inverse metric, and is the metric determinant. Recall that an overset bar denotes a physical quantity, as opposed to a comoving one. For the gravitational action , consider the usual EinsteinHilbert action
(B3) 
Here is the Ricci scalar determined from the LeviCivita connection compatible with and we are using MKS units with .
In comoving coordinates, the RW metric takes the form of Eqn. (C7). Defining the conformal time through
(B4) 
Eqn. (C7) becomes
(B5) 
Note that Eqn. (B5) takes the form of a conformal transformation of flat spacetime
(B6) 
where is the Minkowski metric of special relativity. We may thus express in the conformally flat frame (Wald, 2010, §D)
(B7)  
(B8) 
with regular partial derivatives. Similarly, the metric determinant in the conformally flat frame becomes
(B9) 
The gravitational action then becomes
(B10) 
and its variation with respect to the single physical degree of freedom is
(B11) 
where we have discarded the boundary term.
To find , we start with the most general fluid stresstensor, which applies to any type of matter and coordinate choice. It is given by Hu (2004, Eqn. 6) as
(B12) 
where and are futuredirected timelike vector fields. We have eliminated the purely formal separation of background from perturbation for clarity. Now, we may lower an index in the usual way
(B13) 
and use the chain rule to write
(B14) 
Substitution of Eqns. (B6), (B9), (B12), and (B13) into Eqn. (B2) then gives
(B15) 
We may compress this further by defining the following notation
(B16) 
to arrive at the desired expression
(B17) 
We may now make explicit our claim that the Friedmann source is a spatialslice average. Under the constraints of isotropy and homogeneity, can only be a function of . This is because the variation must be consistent with the constraints imposed on the system under study (e.g. Lanczos, 2012). We may thus expand the iterated integral of Eqn. (B17)
(B18)  
(B19) 
where is some fiducial spatial volume. Consider a typical arbitrary variation without symmetry requirements, i.e. . Since the variation is arbitrary, one is free to consider variations with compact spatiotemporal support. In such a case, the integrals appearing in the action and its variation can always be regarded as finite. In the RW setting, however, we are not free to enforce compact spatial support for . This means that is infinite, and convergence issues might arise from an iterated integral representation (i.e. Fubini’s theorem). This can be remedied by regarding as resulting from an arbitrary spatial cutoff for the integration domain. One then verifies that the result is independent of the cutoff.
Indeed, substitution of Eqns. (B11) and (B19) into Eqn. (B1) gives
(B20) 
where the fiducial spatial volume successfully divides off. The resulting field equation
(B21) 
agrees with the standard Eqns. (C8) and (C9) summed and shifted to conformal time. Note that, in the standard scenario, all are equal, and so is the usual isotropic pressure. The advantage of working directly from the action is that the spatial average of the source is made manifest. This is the main result of this appendix.
b.1 Discussion
The spatial average of Eqn. (B21) is entirely consistent with the operational definition given by Peebles (1980, §81). The temporary appearance of the spatial volume is interesting. Einstein’s equations are “local” in the sense of “localized in space.” This leads to the intuition that all dynamical fields are positiondependent. In the context of cosmology, this is reinforced from the Newtonian approximation, where defines a coordinate change (e.g. Peebles, 1980), which can vary from point to point. Together, these predispositions have sometimes (e.g. Rácz et al., 2017) led to the perception that there is some sort of local scale factor , which the Friedmann model somehow digests into .
One can decompose Eqn. (B17) into a sum of open, disjoint, spacelike regions and average over these
(B22) 
Here, the cutoff procedure discussed previously is equivalent to requiring . Until is large enough for a given epoch, these averages will not be equal and so the equations of motion will not be welldefined. This clarifies as the lower bound lengthscale implicit to the isotropic perfect fluid model (Landau & Lifshitz, 1959, §1). The novelty within cosmology is that this lengthscale is, in general and practice, timedependent. Technically, one should not assume the same perfect fluid model throughout the entire history of the universe.
Appendix C Consistency of DFFSCDM and CFSCdm
In this Appendix, we demonstrate that CDM with the CFS remains observationally viable. As established in §4.2 and discussed in Appendix A, there is no need to track pressures until the nonlinear growth regime. Before this regime, therefore, the DFFS and the CFS are equivalent. Once into the nonlinear regime, predictions are not made with linear perturbation theory. Instead, body simulations with initial conditions determined by linear perturbation theory are required. We thus show that the standard Newtonian behaviour on an expanding background is observationally unchanged. To do this, we show that the “Hubble drag” and density contrast are observationally equivalent in the DFFS and the CFS.
Our approach in this section will be to conservatively bound the influence of pressures from typical systems on both cosmological and peculiar motions. We will quantify the influence on cosmological motion with the fractional difference in the Hubble rate
(C1) 
We will quantify the influence on peculiar motions with the difference in total comoving energy densities
(C2) 
Here is defined in Eqn. (C10) as the total density appearing in the Friedmann energy Eqn. (C8). We will use DFFSCDM (with purely collisionless matter) as the reference expansion, and consider only Planck bestfit cosmological parameters for and (Planck Collaboration et al., 2016).
To build an extremely conservative upper bound on the effect of pressures neglected in the DFFS, we consider the densest stars, the neutron stars. While the neutron star equation of state is an active area of research, example stable polytrope models (e.g. Hansen et al., 2012) give a local core pressure of . For simplicity, we will regard the equation of state parameter as constant during the following analysis.
We define a depletion fraction by dividing the cosmic stellar density, given in derivative form in Eqn. (E28), by . We then numerically integrate Eqn. (C14). We choose extremely conservative (i.e. large) values for the equation of state parameter
(C3) 
The results for the expansion history are shown in the bottom panel of Figure 7. It is clear that keeping positive pressures in the CFS cannot observably alter cosmological motions.
From Eqn. (C1) and Figure 7, the expansion loses energy relative to the DFFS with fixed comoving densities. From conservation of energy and momentum given in Eqn. (12),
(C4) 
it is clear that can only remain constant at fixed when . Therefore, a small amount of energy density is subsequently lost after the conversion of collisionless stress to stress under compression . This loss, relative to the DFFS, is shown in the top panel of Fig. 7. To understand the influence of this loss, recall that peculiar motions are sourced by the energy density contrast (Peebles, 1980, §7)
(C5) 
Here is the total local density, is the Newtonian gravitational potential, and is the (background) comoving energy density defined in Eqn. (C10). In bound systems at late times,
(C6) 
It is immediately clear that can have no observable effect on the local dynamics of any contemporary system.
In summary, we have demonstrated that even conservative (i.e. large) estimates of pressure contributions to the CFS produce no observable discrepancies from the pressurefree DFFS. This is true both for cosmological motions, as characterized by , and for peculiar motions as determined by the density contrast. These results establish that the CFS is observationally consistent with the DFFS under the assumptions of CDM.
c.1 Friedmann’s equations
Here we derive the form of Friedmann’s equations used in our numerical analysis. Units will remain explicit until Eqn. (C14), which is written with the units used elsewhere in this paper. We begin with the spatially flat RobertsonWalker line element
(C7) 
with dimensionless. Substitution of Eqns. (9), (10), and (C7) into Einstein’s equations with cosmological constant term gives the Friedmann equations
(C8)  
(C9) 
where dot denotes derivative with respect to coordinate time. Here and . Define the total comoving energy density , effective equation of state , and scale velocity by
(C10)  
(C11)  
(C12) 
Substitution of Eqns. (C10) and (C12) into Eqn. (C8) gives an algebraic constraint
(C13) 
Substitution of Eqns. (11), (C11), and (C13) into Eqn. (C9) gives a single dynamical equation, which we express in terms of