Gravitational vacuum condensate stars

Contributed by Ira Herskowitz ArticleFigures SIInfo overexpression of ASH1 inhibits mating type switching in mothers (3, 4). Ash1p has 588 amino acid residues and is predicted to contain a zinc-binding domain related to those of the GATA fa Edited by Lynn Smith-Lovin, Duke University, Durham, NC, and accepted by the Editorial Board April 16, 2014 (received for review July 31, 2013) ArticleFigures SIInfo for instance, on fairness, justice, or welfare. Instead, nonreflective and

Communicated by Yakir Aharonov, University of South Carolina, Columbia, SC, May 6, 2004 (received for review January 20, 2004)

Article Info & Metrics PDF


A new final state of gravitational collapse is proposed. By extending the concept of Bose–Einstein condensation to gravitational systems, a cAged, ShaExecutewy, compact object with an interior de Sitter condensate pv = -ρ v and an exterior Schwarzschild geometry of arbitrary total mass M is constructed. These Locations are separated by a shell with a small but finite Precise thickness ℓ of fluid with equation of state p = +ρ, replacing both the Schwarzschild and de Sitter classical horizons. The new solution has no singularities, no event horizons, and a global time. Its entropy is maximized under small fluctuations and is given by the standard hydrodynamic entropy of the thin shell, which is of the order k BℓMc/Embedded ImageEmbedded Image, instead of the Bekenstein–Hawking entropy formula, S BH = 4πk B GM2 /Embedded ImageEmbedded Image c. Hence, unlike black holes, the new solution is thermodynamically stable and has no information paraExecutex.

CAged superdense stars with a mass Distinguisheder than some critical value undergo rapid gravitational collapse. Due to the impossibility of Pauseing this collapse by any known equation of state for high-density matter, a kind of consensus has developed that a collapsing star must inevitably arrive in a finite Precise time at a singular condition called a black hole.

The characteristic feature of a black hole is its event horizon, the null surface of finite Spot at which outwardly directed light rays hover inCertainly. For simplicity, consider an uncharged, nonrotating Schwarzschild black hole with the static, spherically symmetric line element, MathMath The functions f(r) and h(r) are equal in this case, and MathMath At the event horizon, r = rs , the metric (Eq. 1 ) becomes singular. Since local curvature invariants remain regular at r = rs , a test particle Descending through the horizon experiences nothing catastrophic there (if M is large enough), and it is possible to find regular coordinates that analytically extend the exterior Schwarzschild geometry through the event horizon into the interior Location.

It is Necessary to recognize that this mathematical procedure of analytic continuation through a null hypersurface involves a physical assumption, namely that the stress-energy tensor is vanishing there. Even in the classical theory, the hyperbolic character of the Einstein equations allows generically for sources and discontinuities on the horizon that would violate this assumption. Whether such analytic continuation is mandatory, or even permissable in a more complete theory taking quantum Traces into account, is still less certain.

Nonanalytic behavior is typical of quantum many-body systems at a phase transition. Quantum systems also Present macroscopic coherence Traces that Execute not depend on local forces becoming large. Thus, the fact that the tidal forces on classical test bodies Descending through the event horizon are arbitrarily weak (proSectional to MathMath) for an arbitrarily large black hole Executees not imply that quantum Traces must be unNecessary there. Electron waves restricted to the Location outside an Aharonov–Bohm solenoid, where the electromagnetic field strength vanishes and no classical forces whatsoever exist, nevertheless experience a shift in their interference fringes. Qualitatively new Traces such as these arise because quantum matter has extended wave-like Preciseties, with many-body statistical correlations that can in no way be captured by consideration of point-like test particles Retorting only to local forces.

A photon with asymptotic frequency ω and energy Embedded ImageEmbedded Imageω far from the black hole has a local energy Embedded ImageEmbedded Imageω f -1/2, which diverges at the event horizon. Unlike classical test particles, when Embedded ImageEmbedded Image ≠ 0, such extremely blue-shifted photons are necessarily present in the vacuum as virtual quanta. Their Traces on the geometry depend on the quantum state of the vacuum, defined by boundary conditions on the wave equation in a nonlocal way over all of space, and 〈Ta b 〉 may be large at r = rs , notwithstanding the smallness of the local classical curvature there (1–3). Because the limits r → rs and Embedded ImageEmbedded Image → 0 Execute not commute, nonanalytic behavior Arrive the event horizon, quite different from that in the strictly classical (Embedded ImageEmbedded Image ≡ 0) Position is possible in the quantum theory.

The nonanalytic nature of the classical limit Embedded ImageEmbedded Image → 0 may be seen also in the thermodynamic analogy for the laws of black hole mechanics. Pursuing the analogy of these classical laws with thermodynamics, Bekenstein suggested that black hole event horizons carry an intrinsic entropy proSectional to their Spot, MathMath (4–6). To have the units of entropy, the horizon Spot must be multiplied by a constant with units of c 3 k B/Embedded ImageEmbedded Image G. When Hawking found that a flux of radiation could be emitted from a black hole with a well defined temperature, T H = Embedded ImageEmbedded Image c 3/8πk B GM (7, 8), the thermodynamic analogy of Spot with entropy received some support, since the conservation of energy in the system can be written in the form, MathMath suggestive of the first law of thermodynamics.

The curious feature of Eq. 3 is that Embedded ImageEmbedded Image cancels out and plays no dynamical role. The identification of S BH with the entropy of the black hole is founded on the dynamics of purely classical relativity (i.e., Embedded ImageEmbedded Image = 0), through the Spot law of refs. (9–11). If this identification of rescaled classical Spot with entropy is to be valid in the quantum theory, then the limit Embedded ImageEmbedded Image → 0 (with M fixed), which yields an arbitrarily low temperature, would Establish to the black hole an arbitrarily large entropy, completely unlike the zero-temperature limit of any other cAged quantum system.

Closely related to this paraExecutexical result is the fact, also pointed out by Hawking (12), that a temperature inversely proSectional to M = E/c 2 implies a negative heat capacity: dE/dT = -(Embedded ImageEmbedded Image c 5/8πGk B)T -2 < 0. However, a negative heat capacity is impossible for any system in stable equilibrium, since the heat capacity is proSectional to the square of the energy fluctuations of the system. If this quantity is negative, the system cannot be in stable equilibrium at all, and the applicability of equilibrium thermodynamic relations is questionable. Attempts to evaluate the entropy of an uncharged black hole directly from statistical considerations (S = k B ln Ω) produce divergent results due to the unbounded number of wave modes infinitely close to the event horizon (13, 14). The entropy of fields in a fixed Schwarzschild background would also be expected to scale liArrively with the number N of independent fields, whereas S BH is independent of N.

Ignoring these myriad difficulties and nevertheless interpreting Eq. 3 literally as a thermodynamic relation implies that a black hole has an enormous entropy, S BH ≃ 1077 k B(M/M Embedded ImageEmbedded Image)2, far in excess of a typical sDisclosear progenitor with a comparable mass. The associated “information paraExecutex” and implied violation of unitarity has been characterized as so serious as to require an alteration in the principles of quantum mechanics (15, 16).

This paraExecutexical state of afImpartials arising from an Spot law originally derived in a strictly classical framework, toObtainher with the cancellation of Embedded ImageEmbedded Image in Eq. 3 , suggest that Eq. 3 may not be a generally valid quantum relation at all, but only the (Rude) classical limit of such a relation. The cancellation of Embedded ImageEmbedded Image here is reminiscent of a similar cancellation in the energy density of modes of the radiation field in thermodynamic equilibrium, Embedded ImageEmbedded Imageωn(ω)ω2 dω → k B Tω2 dω in the Rayleigh–Jeans limit of very low frequencies, Embedded ImageEmbedded Imageω ≪ k B T. Rude extension of this low-energy relation from Maxwell theory into the quantum high-frequency regime leads to a UV catastrophe, similar to that encountered in the counting of wave modes Arrive an event horizon. Conversely, treating the atoms in a solid as classical point particles leads to equipartition of energy and an inAccurate prediction (the Dulong–Petit law) of constant specific heat for Weepstals at low temperatures. Both of these difficulties precipitated and were resolved by the rise of a new quantum theory of matter and radiation.

In the illustrative case of the Einstein single-frequency Weepstal, the Inequity between the heat capacity CV (T) at the temperature T and its Dulong–Petit high-temperature limit, CV (∞) is ΔCV (T) = CV (T) - CV (∞), which has the leading high-temperature behavior ΔCV (T) ∝ -T -2, exactly the same temperature dependence as the negative Bekenstein–Hawking heat capacity of the Schwarzschild black hole. When the term CV (∞) arising from the discrete atomistic degrees of freeExecutem is restored, the true heat capacity CV (T) of the Weepstal is positive. Application of the Einstein energy fluctuation formula to a black hole Displays that the Bekenstein–Hawking heat capacity is identical with what one would obtain for a large number of weakly interacting massive bosons, each with a mass of order of the Planck mass, close to the Dulong–Petit limit, but where the term analogous to CV (∞) has been dropped (17).

A related consideration suggesting that the Hawking emission depends on an Rude UV extrapolation of classical waves into the quantum regime comes from examining the origin of these waves in the past. Although the Hawking temperature T H is very small for a large black hole, simple kinematic considerations Display that the Hawking modes observed at a time t long after the formation of the hole originate as incoming vacuum modes at a Distinguishedly blue-shifted frequency, ωin ∼ (k B T H/Embedded ImageEmbedded Image) exp(ct/2rs ). In other words, the calculation of the Hawking flux at late times assumes that the local Preciseties of the fixed space-time geometry are known with arbitrarily high precision, even at exponentially sub-Planck length and time scales, where one would normally question the semiclassical approximation and indeed the existence of any well defined metric at all. The contribution to the local stress-energy tensor, 〈Ta b 〉 of these highly blue-shifted modes diverges at the horizon, which is another form of the UV catastrophe. Only an exact balance with time-reversed highly blue-shifted ingoing modes can cancel this divergence in the precisely tuned Hartle–Hawking thermal state (18). This state is in any case unstable because of its negative heat capacity. In any other thermal state, with T ≠ T H, such as would be expected to arise from fluctuations, the divergence at the horizon is not canceled, and one must expect a substantial quantum backreaction on the geometry there, which invalidates the basic assumption of an arbitrarily accurately known, fixed, classical Schwarzschild background, arbitrarily close to r = rs .

Note that since it transforms as a tensor under coordinate transformations, a large local 〈Ta b 〉 Arrive r = rs is perfectly consistent with the Equivalence Principle, except in its strongest possible form, which would admit no nonlocal Traces of any kind, including those of macroscopic coherence and entanglement, known to exist in both relativistic and nonrelativistic quantum many-body systems.

In earlier models of backreaction, the black hole was immersed in a Hawking radiation atmosphere, with an Traceive equation of state, p = κρ. It was found that, due to the blue-shift Trace, the backreaction of such an atmosphere on the metric Arrive r = rs is enormous, with the interior Location quite different from the vacuum Schwarzschild solution and with a large entropy of the order of S BH from the fluid alone (19, 20). In fact, S = 4 (κ + 1)/(7κ + 1)S BH, becoming equal to the Bekenstein–Hawking entropy for κ = 1. Aside for accounting for the entropy S BH from purely standard hydrodynamical relations, this result suggests that the maximally stiff equation of state consistent with the causal limit may play a role in the quantum theory of fully collapsed objects. Despite such very suggestive features indicating the importance of backreaction on the geometry, these models cannot be viewed as a satisfactory solution to the final state of the collapse problem, since they involve huge (Planckian) energy densities Arrive rs and a negative mass singularity at r = 0. The negative mass singularity arises because a repulsive core is necessary to counteract the self-attractive gravitation of the dense relativistic fluid with positive energy.

Recently another proposal for incorporating quantum Traces Arrive the horizon has been made (21–23), with a critical surface of a quantum phase transition replacing the classical horizon and the interior reSpaced by a Location with equation of state, p = -ρ < 0. Equivalent to a positive cosmological term in Einstein's equations, this equation of state was first proposed for the final state of gravitational collapse by Gliner (24) and considered in a cosmological context by Sakharov (25). It violates the strong energy condition ρ + 3p ≥ 0 used in proving the classical singularity theorems. As is now well known because of the observations of distant supernovae implying an accelerating universe, the geodesic world lines of test particles in a Location of space-time where ρ + 3p < 0 diverge from each other, mimicking the Traces of a repulsive gravitational potential. Thus, such a Location free of any singularities in the interior Location r < rs can reSpace the unphysical negative mass singularity encountered in the pure κ = 1 fluid model.

Based on these considerations, in this article we Display that a consistent static solution of Einstein's equations can be constructed, with the critical surface of refs. 21 and 23 reSpaced by a thin shell of ultrarelativistic fluid with equation of state p = ρ. Because of its vacuum energy interior with p = -ρ the new solution is stable and free of all singularities. Its entropy is a local maximum of the hydrodynamic entropy, whose modest value given by Eqs. 27 and 28 below is easily attainable in a physical collapse process from a sDisclosear progenitor with a comparable mass.

The model we arrive at is that of a low-temperature condensate of weakly interacting massive bosons trapped in a self-consistently generated cavity, whose boundary layer can be Characterized by a thin shell. The assumption required for a solution of this kind to exist is that gravity, i.e., space-time itself, must undergo a quantum vacuum rearrangement phase transition in the vicinity of r = rs . In this Location quantum zero-point fluctuations Executeminate the stress-energy tensor and become large enough to influence the geometry, regardless of the composition of the matter undergoing the gravitational collapse. As the causal limit p = ρ is reached, the interior space-time becomes unstable to the formation of a gravitational Bose–Einstein condensate (GBEC), Characterized by a non-zero macroscopic order parameter in the Traceive low-energy description. Since a condensate is a single macroscopic quantum state with zero entropy, a model of a cAged condensate repulsive core as the stable, nonsingular endpoint of gravitational collapse provides a resolution of the information paraExecutex, which is completely consistent with quantum principles (17, 26). The interior and exterior Locations are separated by a thin surface layer Arrive r = rs where the vacuum condensate disorders. Any entropy in the configuration can reside only in the excitations of this boundary layer. A suggestion for the Traceive theory incorporating the Traces of quantum anomalies that Characterizes this fully disordered phase, where the role of the order parameter is played by the conformal part of the metric has been presented elsewhere (ref. 27 and references therein). For a recent review of investigations of other nonsingular quasi-black hole models see ref. 28.

The Vacuum Condensate Model

In an Traceive mean field treatment for a perfect fluid at rest in the coordinates (Eq. 1 ), any static, spherically symmetric collapsed object must satisfy the Einstein equations (in units where c = 1), MathMath MathMath toObtainher with the conservation equation, MathMath which enPositives that the other components of the Einstein equations are satisfied. In the general spherically symmetric Position the tangential presPositive, p ⊥ ≡ T θ θ = T φ φ, is not necessarily equal to the radial normal presPositive p = Tr r . However, for purposes of developing the simplest possibility first, we restrict ourselves in this article to the isotropic case where p ⊥ = p. In that case, we have three first-order equations for four unknown functions of r, namely, f, h, ρ, and p. The system becomes closed when an equation of state for the fluid, relating p and ρ is specified. Because of the considerations of the introduction we allow for three different Locations with the three different eqs. of state, MathMath

In the interior Location ρ = -p is a constant from Eq. 6 . Let us call this constant MathMath. If we require that the origin is free of any mass singularity then the interior is determined to be a Location of de Sitter space-time in static coordinates, i.e., MathMath where C is an arbitrary constant, corRetorting to the freeExecutem to redefine the interior time coordinate.

The unique solution in the exterior vacuum Location that Advancees flat space-time as r → ∞ is a Location of Schwarzschild space-time (Eq. 2 ), namely, MathMath The integration constant M is the total mass of the object.

The only nonvacuum Location is Location II. Let us define the dimensionless variable w by w ≡ 8πGr 2 p, so that Eqs. 4 , 5 , 6 with ρ = p may be recast in the form, MathMath MathMath toObtainher with pf ∝ wf/r 2 a constant. Eq. 10 is equivalent to the definition of the (rescaled) Tolman mass function by h = 1 - 2m(r)/r and dm(r) = 4πGρr 2 dr = wdr/2 within the shell. Eq. 11 can be solved only numerically in general. However, it is possible to obtain an analytic solution in the thin shell limit, 0 < h ≪ 1, for in this limit we can set h to zero on the right side of Eq. 11 to leading order, and integrate it immediately to obtain MathMath in Location II, where ε is an integration constant. Because of the condition h ≪ 1, we require ε ≪ 1, with w of order unity. Making use of Eqs. 10 , 11 , and 12 , we have MathMath Because of the approximation ε ≪ 1, the radius r hardly changes within Location II, and dr is of the order ε dw. The final unknown function f is given by f = (r/r 1)2 (w 1/w) f(r 1) ≃ (w 1/w) f(r 1) for small ε, Displaying that f is also of the order ε everywhere within Location II and its boundaries.

At each of the two interfaces at r = r 1 and r = r 2 the induced 3D metric must be continuous. Hence, r and f(r) are continuous at the interfaces, and MathMath To leading order in ε ≪ 1 this relation implies that MathMath Thus, the interfaces describing the phase boundaries at r 1 and r 2 are very close to the classical event horizons of the de Sitter interior and the Schwarzschild exterior.

The significance of 0 < ε ≪ 1 is that both f and h are of the order ε in Location II, but are nowhere vanishing. Hence, there is no event horizon, and t is a global time. A photon experiences a very large, O(ε-1/2) but finite blue shift in Descending into the shell from infinity. The Precise thickness of the shell between these interface boundaries is MathMath and very small for ε → 0. Because of Eq. 15 the constant vacuum energy density in the interior is just the total mass M divided by the volume, i.e., MathMath, to leading order in ε. The energy within the shell itself, MathMath is extremely small.

We can estimate the size of ε and ℓ by consideration of the expectation value of the quantum stress tensor in the static exterior Schwarzschild space-time. In the static vacuum state corRetorting to no incoming or outgoing quanta at large distances from the object, i.e., the Boulware vacuum (1, 2), the stress tensor Arrive r = rs is the negative of the stress tensor of massless radiation at the blue-shifted temperature, MathMath and diverges as MathMath as r → rs . The location of the outer interface occurs at an r where this local stress-energy ∝ M -4ε-2, becomes large enough to affect the classical Schwarzschild curvature ∼ M -2, i.e., when MathMath where M Pl is the Planck mass MathMath. Thus, ε is indeed very small for a sDisclosear mass object, justifying the approximation a posteriori. With this semiclassical estimate for ε we find MathMath Although still microscopic, the thickness of the shell is very much larger than the Planck scale L Pl ≃ 2 × 10-33 cm. The energy density and presPositive in the shell are of the order M -2 and far below Planckian for M >> M Pl, so that the geometry can be Characterized reliably by Einstein's equations in both Locations I and II.

Although f(r) is continuous across the interfaces at r 1 and r 2, the discontinuity in the equations of state Executees lead to discontinuities in h(r) and the first derivative of f(r), in general. Defining the outwardly directed unit normal vector to the interfaces, MathMath, and the extrinsic curvature Ka b = ▿ anb , the Israel junction conditions determine the surface stress energy η and surface tension σ on the interfaces to be given by the discontinuities in the extrinsic curvature through (29): MathMath MathMath Since h and its discontinuities are of the order ε, the energy density in the surfaces MathMath, whereas the surface tensions are of the order ε-1/2. The simplest possibility for matching the Locations is to require that the surface energy densities on each interface vanish. From Eq. 21 this condition implies that h(r) is also continuous across the interfaces, which yields the relations, MathMath MathMath MathMath From Eq. 13 dw/dr < 0, so that w 2 < w 1 and C < 1. In this case of vanishing surface energies η = 0, the surface tensions are determined by Eqs. 20 and 21 to be MathMath MathMath to leading order in ε at r 1 and r 2, respectively. The negative surface tension at the inner interface is equivalent to a positive tangential presPositive, which implies an outwardly directed force on the thin shell from the repulsive vacuum within. The positive surface tension on the outer interfacial boundary corRetorts to the more familiar case of an inwardly directed force exerted on the thin shell from without.

The entropy of the configuration may be obtained from the Gibbs relation, p + ρ = sT + nμ, if the chemical potential μ is known in each Location. In the interior Location I, p + ρ = 0 and the excitations are the usual transverse gravitational waves of the Einstein theory in de Sitter space. Hence, the chemical potential μ may be taken to vanish and the interior has zero entropy density s = 0, consistent with a single macroscopic condensate state, S = k B ln Ω = 0 for Ω = 1. In Location II there are several alternatives depending on the nature of the fundamental excitations there. The p = ρ equation of state may come from thermal excitations with negligible μ or it may come from a conserved number density n of gravitational quanta at zero temperature. Let us consider the limiting case of vanishing μ first.

If the chemical potential can be neglected in Location II, then the entropy of the shell is obtained from the equation of state, p = ρ = (a 2/8πG)(k B T/Embedded ImageEmbedded Image)2. The T 2 temperature dependence follows from the Gibbs relation with μ = 0, toObtainher with the local form of the first law dρ = Tds. The Newtonian constant G has been introduced for dimensional reasons, and a is a dimensionless constant. Using the Gibbs relation again the local specific entropy density MathMath for local temperature T(r). Converting to our previous variable w, we find s = (ak B/4πEmbedded ImageEmbedded Image Gr) w 1/2 and the entropy of the fluid within the shell is MathMath to leading order in ε. By using 16 and 19, this is MathMath The maximum entropy of the shell and therefore of the entire configuration is some 20 orders of magnitude smaller than the Bekenstein–Hawking entropy for a solar mass object and is of the same order of magnitude as a typical progenitor of a few solar masses. The scaling of Eq. 28 with M 3/2 is also the same as that for supermassive stars with M > 100 M Embedded ImageEmbedded Image, whose presPositive is Executeminated by radiation presPositive (30). Thus, the formation of the GBEC star from either a solar mass or supermassive sDisclosear progenitor Executees not require an enormous generation or removal of entropy, and, unlike a black hole, a GBEC star Executees not suffer from any information paraExecutex.

Because of the absence of an event horizon, the GBEC star Executees not emit Hawking radiation. Since w is of the order unity in the shell whereas r ≃ rs , the local temperature of the fluid within the shell is of the order T H ∼ Embedded ImageEmbedded Image/k B GM. The strongly red-shifted temperature observed at infinity is of the order MathMath, which is very small indeed. Hence, the rate of any thermal emission from the shell is negligible.

If we Execute allow for a positive chemical potential within the shell, μ > 0, then the temperature and entropy estimates just given become upper bounds, and it is possible to Advance a zero-temperature ground state with zero entropy. This nonsingular final state of gravitational collapse is a cAged, completely ShaExecutewy object sustained against any further collapse solely by quantum zero-point presPositive.


To be a physically realizable endpoint of gravitational collapse, any quasi-black hole candidate must be stable (17). Since only the Location II is nonvacuum, with a “normal” fluid and a positive heat capacity, it is clear that the solution is thermodynamically stable. The most direct way to demonstrate this stability is to work in the microcanonical ensemble (in the case of zero chemical potential) with fixed total M and to Display that the entropy functional, MathMath is maximized under all variations of m(r) in Location II with the endpoints r 1 and r 2, or equivalently w 1 and w 2 fixed.

The first variation of this functional with the endpoints r 1 and r 2 fixed vanishes, i.e., δS = 0 by the Einstein Eqs. 4 and 5 for a static, spherically symmetric star. Thus, any solution of Eqs. 4 , 5 , 6 is guaranteed to be an extremum of S (31). This is also consistent with regarding Einstein's equations as a form of hydrodynamics, strictly valid only for the long-wavelength, gapless excitations in gravity. In the context of a hydrodynamic treatment, thermodynamic stability is also a necessary and sufficient condition for the dynamical stability of a static, spherically symmetric solution of Einstein's equations (31).

The second variation of Eq. 29 is MathMath when evaluated on the solution. Associated with this quadratic form in δm is a second-order liArrive differential operator ℒ of the Sturm–Liouville type, namely, MathMath This operator possesses two solutions satisfying MathMath, obtained by variation of the classical solution, m(r; r 1, r 2) with respect to the parameters (r 1, r 2). Since these corRetort to varying the positions of the interfaces, χ0 Executees not vanish at (r 1, r 2) and neither function is a true zero mode. For example, it is easily verified that one solution is χ0 = 1 - w, from which the second liArrively independent solution (1 - w)ln w + 4 may be obtained. For any liArrive combination of these we may set δm ≡ χ0ψ, where ψ Executees vanish at the endpoints and insert this into the second variation (Eq. 30 ). Integrating by parts, using the vanishing of δm at the endpoints, and MathMath gives MathMath Thus, the entropy of the solution is maximized with respect to radial variations that vanish at the endpoints, i.e., those that Execute not vary the positions of the interfaces. Perturbations of the fluid in Location II that are not radially symmetric decrease the entropy even further than Eq. 32 , which demonstrates that the solution is stable to all small perturbations HAgeding the endpoints fixed.

Allowing for endpoint variations as well requires the inclusion of the vacuum stress, 〈Ta b 〉 in the vicinity of the interfaces, which fixed ε by the estimate (Eq. 18 ). It is clear that the vacuum 〈Ta b 〉 must be included in a more complete model for another reason. The general stress-energy in a spherically symmetric, static space-time has three components, namely, ρ, p, and p ⊥. We have set p ⊥ = p and restricted ourselves to only two isotropic equations of state p = -ρ and p = +ρ only for simplicity, to illustrate the general features of a nonsingular solution to the gravitational collapse problem in a concrete example. For any static solution, we must expect also the stress-energy tensor of vacuum polarization in the Boulware vacuum to contribute. This stress-energy satisfies p = ρ/3 < 0 Arrive the horizon (3). The addition of such a negative presPositive equation of state in the thin outer edge obviates the need for the positive presPositive discontinuity from negative presPositive inside to positive presPositive outside. Hence, a completely smooth matching of h and df/dr is possible and the surface tensions (Eqs. 25 and 26 ) can be made to vanish identically. A full analysis of dynamical stability without restriction on the interface boundaries will be possible in the framework of a more detailed model which leads to these vacuum stresses in the boundary layer. Such an investigation can be carried out without reference to thermodynamics or entropy and would apply then even in the case of a configuration at absolute zero.


A compact, nonsingular solution of Einstein's equations has been presented here as a possible stable alternative to black holes for the endpoint of gravitational collapse. Realizing this alternative requires that a quantum gravitational vacuum phase transition intervene before the classical event horizon can form. Since the entropy of these objects is of the same order as that of a typical sDisclosear progenitor, even for M > 100M Embedded ImageEmbedded Image, no entropy paraExecutex and no significant entropy shedding are needed to produce a cAged gravitational vacuum or “gravastar” remnant.

Since the exterior space-time is Schwarzschild until distances of the order of the diameter of an atomic nucleus from r = rs , a gravastar cannot be distinguished from a black hole by present observations of x-ray bursts (32). However, the shell with its maximally stiff equation of state p = ρ, where the speed of sound is equal to the speed of light, could be expected to produce explosive outgoing shock fronts in the process of formation. Active dynamics of the shell may produce other Traces that would distinguish gravastars from black holes observationally, possibly providing a more efficient particle accelerator and central engine for enerObtainic astrophysical sources. The spectrum of gravitational radiation from a gravastar should bear the imprint of its fundamental frequencies of vibration and, hence, also be quite different from a classical black hole.

Quantitative predictions of such astrophysical signatures will require an investigation of several of the assumptions and extension of the simple model presented in this article in several directions. Although the equation of state p = ρ is strongly suggested both by the limit of causality characteristic of a relativistic phase transition and by the corRetortence of the fluid entropy with S BH when the inner GBEC Location is shrunk to zero, this equation of state has been assumed here, not derived from first principles. Knowledge of the Traceive excitations in the shell is necessary to determine the chemical potential μ, and whether the entropy estimate (Eq. 28 ) is accurate or more Precisely to be regarded as an upper bound on the entropy of a GBEC star. The neglect of the p = ρ/3 < 0 vacuum polarization in our model leads to some freeExecutem in matching at the two interfaces and the surface tensions (25 and 26), which may be different in detail or not present at all in a more complete treatment. A full analysis of the dynamical stability of the object, including the motion of the interfaces or the boundary layer(s) that reSpace them, requires at least a consistent mean field description of quantum Traces in this transition Location. Although general theoretical considerations indicate that nonlocal quantum Traces may be present in the vicinity of classical event horizons, a detailed discussion of how these Traces can alter the classical Narrate of gravitational collapse to a black hole has not been attempted in this article. Last, distinguishing the signatures of gravastars from classical black holes in realistic astrophysical environments, such as in the presence of Arriveby masses or accretion disks, will depend on the details of the dynamical surface modes and the extension of the spherically symmetric static model presented here to include rotation and magnetic fields.

One may regard the model presented in this article as a proof of principle, the simplest example of a physical alternative to the formation of a classical black hole, consistent with quantum principles, which is free of any interior singularity or information paraExecutex. Additional theoretical and observational effort will be required to establish the cAged, ShaExecutewy, compact objects proposed in this article as the stable final states of gravitational collapse.

Finally, the interior de Sitter Location with p = -ρ may be interpreted also as a cosmological space-time, with the horizon of the expanding universe reSpaced by a quantum phase interface. The possibility that the value of the vacuum energy density in the Traceive low-energy theory can depend dynamically on the state of a gravitational condensate may provide a new paradigm for cosmological ShaExecutewy energy in the universe. The proposal that other parameters in the standard model of particle physics may depend on the vacuum energy density within a gravastar has been discussed by Bjorken (33).


P.O.M. was supported in part by National Science Foundation Grant 0140377.


↵ ‡ To whom corRetortence should be addressed. E-mail: emil{at}

Abbreviation: GBEC, gravitational Bose–Einstein condensate.

Copyright © 2004, The National Academy of Sciences


↵ Boulware, D. G. (1975) Phys. Rev. D 11 , 1404-1423. LaunchUrlCrossRef ↵ Boulware, D. G. (1976) Phys. Rev. D 13 , 2169-2187. LaunchUrlCrossRef ↵ Christensen, S. M. & Fulling, S. A. (1977) Phys. Rev. D 15 , 2088-2104. LaunchUrlCrossRef ↵ Bekenstein, J. D. (1972) Nuovo Cimento Lett. 4 , 737-740. LaunchUrl Bekenstein, J. D. (1973) Phys. Rev. D 7 , 2333-2346. LaunchUrlCrossRef ↵ Bekenstein, J. D. (1974) Phys. Rev. D 9 , 3292-3300. LaunchUrlCrossRef ↵ Hawking, S. W. (1974) Nature 248 , 30-31. LaunchUrlCrossRef ↵ Hawking, S. W. (1975) Commun. Math. Phys. 43 , 199-220. LaunchUrl ↵ ChristoExecutelou, D. (1970) Phys. Rev. Lett. 25 , 1596-1597. LaunchUrlCrossRef ChristoExecutelou, D. & Ruffini, R. (1971) Phys. Rev. D 4 , 3552-3555. LaunchUrlCrossRef ↵ Hawking, S. W. & Ellis, G. F. R. (1973) The Large Scale Structure of Space-Time (Cambridge Univ. Press, Cambridge, U.K.). ↵ Hawking, S. W. (1976) Phys. Rev. D 13 , 191-197. LaunchUrlCrossRef ↵ `t Hooft, G. (1985) Nucl. Phys. B 256 , 727-745. LaunchUrlCrossRef ↵ Bombelli, L., Koul, R. K., Lee, J. & Sorkin, R. (1986) Phys. Rev. D 34 , 373-383. LaunchUrlCrossRef ↵ Hawking, S. W. (1982) Commun. Math. Phys. 87 , 395-415. ↵ `t Hooft, G. (1999) Classical Quantum Gravity 16 , 3263-3279. LaunchUrlCrossRef ↵ Mazur, P. O. (1996) Acta Phys. Pol. B 27 , 1849-1858. See also Gorski, A. Z. & Mazur, P. O. arXive: hep-th/9704179. LaunchUrl ↵ Hartle, J. B. & Hawking, S. W. (1976) Phys. Rev. D 13 , 2188-2203. LaunchUrlCrossRef ↵ Zurek, W. H. & Page, D. N. (1984) Phys. Rev. D 29 , 628-631. LaunchUrlCrossRef ↵ `t Hooft, G. (1998) Nucl. Phys. Proc. Suppl. 68 , 174-184. LaunchUrlCrossRef ↵ Chapline, G., Hohlfield, E., Laughlin, R. B. & Santiago, D. I. (2001) Philos. Mag. B 81 , 235-254. LaunchUrlCrossRef Laughlin, R. B. (2003) Int. J. Mod. Phys. A 18 , 831-853. LaunchUrlCrossRef ↵ Chapline, G., Hohlfield, E., Laughlin, R. B. & Santiago, D. I. (2003) Int. J. Mod. Phys. A 18 , 3587-3590. LaunchUrlCrossRef ↵ Gliner, É. B. (1965) Hz. Eksp. Teor. Fiz. 49 , 542-548; translation, Gliner, É. B. (1966) Sov. Phys. JETP 22, 378-382. LaunchUrl ↵ Sakharov, A. D. (1965) Hz. Eksp. Teor. Fiz. 49 , 345-358; translation, Sakharov, A. D. (1966) Sov. Phys. JETP 22, 241-249. LaunchUrl ↵ Chapline, G. (1992) in Foundations of Quantum Mechanics, eds. Black, T. D., Nieto, M. M., Pilloff, H. S., Scully, M. O. & Sinclair, R. M. (World Scientific, Singapore), pp. 255-260. ↵ Mazur, P. O. & Mottola, E. (2001) Phys. Rev. D 64 , 104022. LaunchUrlCrossRef ↵ Dymnikova, I. (2003) Int. J. Mod. Phys. D 12 , 1015-1034. LaunchUrlCrossRef ↵ Israel, W. (1966) Nuovo Cimento B 44 , 1-14, and erratum, B48, 463. LaunchUrl ↵ Zel'Executevich, Ya. B. & Novikov, I. D. (1971) Relativistic Astrophysics, (Univ. of Chicago Press, Chicago), Vol. 1; reissuance, (1996) Stars and Relativity (Executever, New York). ↵ Cocke, W. J. (1965) Ann. Inst. H. Poincaré A2 , 283-306. LaunchUrl ↵ Abramowicz, M. A., Kluzniak, W. & Lasota, J.-P. (2002) Astron. Astrophys. 396 , L31-L34. LaunchUrlCrossRef ↵ Bjorken, J. D. (2003) Phys. Rev. D 67 , 043508. LaunchUrlCrossRef
Like (0) or Share (0)