The stochastic gravitational-wave background in the absence of horizons

Gravitational-wave astronomy has the potential to explore one of the deepest and most puzzling aspects of Einstein’s theory: the existence of black holes. A plethora of ultracompact, horizonless objects have been proposed to arise in models inspired by quantum gravity. These objects may solve Hawking’s information-loss paradox and the singularity problem associated with black holes, while mimicking almost all of their classical properties. They are, however, generically unstable on relatively short timescales. Here, we show that this ‘ergoregion instability’ leads to a strong stochastic background of gravitational waves, at a level detectable by current and future gravitational-wave detectors. The absence of such background in the first observation run of Advanced LIGO already imposes the most stringent limits to date on black-hole alternatives, showing that certain models of ‘quantum-dressed’ stellar black holes can be at most a small percentage of the total population. The future LISA mission will allow for similar constraints on supermassive black-hole mimickers.


I. INTRODUCTION
According to General Relativity and the Standard Model of particle physics, dark and compact objects more massive than ≈ 3M must be black holes (BHs) [1].These are characterized by an event horizon, causally disconnecting the BH interior from its exterior, where observations take place.In classical gravity, BHs can form from gravitational collapse [2], providing a sound and compelling theoretical support for their existence.When quantum effects are included at semiclassical level, however, BHs are not completely "dark", but evaporate by emitting thermal black-body radiation [3].For astrophysical BHs, evaporation is negligible.Therefore, BHs are commonly accepted to exist -with masses in the ranges from ∼ 10M to ∼ 60M (and perhaps larger mass), and from ∼ 10 5 M to ∼ 10 10 M -and play a fundamental role in astronomy and astrophysics [4].
It is sometimes not fully appreciated that BHs are truly "holes" in spacetime, where time "ends" and inside which the known laws of classical physics break down [5].Furthermore, the classical concept of an event horizon seems at clash with quantum mechanics, and the very existence of BHs leads to unsolved conundra such as information loss [5].Thus, in reality, the existence of BHs is an outstanding event, for which one should provide equally impressive evidence [6][7][8].Over the last decades several alternatives and arguments have been put forward according to which -in a quantum theory -BHs would either not form at all, or would just be an ensemble of horizonless quantum states [9][10][11][12][13][14][15].
From a theoretical standpoint, BHs therefore lay at the interface between classical gravity, quantum theory and thermodynamics, and understanding their nature may provide a portal to quantum gravity or other surprises.In the formal mathematical sense, it is impossible to ever show that BHs exist, since in General Relativity their definition requires knowledge of the whole spacetime, including the future [2].However, the newborn gravitational-wave (GW) astronomy allows us to constrain alternatives to BHs to unprecedented level.GW de-tectors can rule out a wide range of models, through observations of inspiralling binaries or the relaxation of the final object forming from a merger [6][7][8][16][17][18].Here, we explore one significant effect that follows from the absence of the most salient feature of a BH, the event horizon.We will show that compact, horizonless spinning geometries would fill the universe with a background of GWs detectable by current and future instruments, through a classical process known as the "ergoregion instability" [19,20].

II. ERGOREGION INSTABILITY
In Einstein's theory, the unique globally vacuum astrophysical solution for a spinning object is the Kerr geometry.It depends on two parameters only: its mass M and angular momentum J = GM 2 χ/c, with G Newton's constant, c the speed of light and |χ| ≤ 1 a pure number.The compact, dark objects in our universe could depart from the Kerr geometry in two distinct ways.The near-horizon structure might change significantly, while retaining the horizon [21][22][23].In coalescing binaries, such effects can be probed by GW measurements of the quadrupole moment, the tidal absorption and deformability [7,8,24,25], and especially the quasinormal oscillation modes of the remnant object [26,27].Here we explore a second (and more subtle) scenario, where the geometry is nearly everywhere the same as that of a BH, but the horizon is absent.Two smoking-gun effects arise in this scenario.First, the late-time ringdown consists generically of a series of slowly damped "echoes" [6,8,25,28].Furthermore, by working as a one-way membrane, horizons act as a sink for external fluctuations, including those inside the ergoregion, where negativeenergy states are possible [20,29].Such states are typically associated with instabilities: their existence allows scattering waves to be amplified, i.e. positive-energy perturbations can be produced, which can travel out of the ergoregion.Energy conservation then requires the negative-energy states inside the ergoregion to grow.In a BH, this piling up can be avoided -80 -60 -40 -20 0 20 40  A out e i ω z +A in e -i ω z FIG. 1. Schematic potential for a non-spinning ultracompact object, as a function the tortoise radial coordinate (in practice, the coordinate time t of a photon).For object radii r0 ∼ r+, the radiation travel time from the photosphere to the surface scales approximately as t0 ∼ tH | log |.The travel time from the surface to the interior is parametrized as tinterior.
by dumping the negative energy into the horizon, thus stabilizing the object.In the absence of a horizon, instead, this process leads to an exponential cascade.As a consequence, spinning BHs are linearly stable, but any horizonless object sufficiently similar to a rotating BH is unstable [19,[30][31][32].

A. Canonical model: perfectly reflecting surface
We start by the simplest model of horizonless geometries: a compact body whose exterior is described by the Kerr metric, and with a perfectly reflective surface.This spacetime defines a natural cavity, i.e. the region between the object's surface and the potential barrier for massless particles (the "photon sphere") [see Fig. 1].In this cavity, negative-energy modes, and thus instabilities, can be excited [15,33].The dynamics is controlled by two parameters: the size of the cavity and the object's angular velocity, which determine how fast the instability grows.The cavity size can be parametrized by the light travel time t 0 (as observed at infinity) between the photon sphere and the object's surface [6][7][8]20].The timescale t 0 also defines a set of possible modes, with fundamental frequencies ω = ω R + iω I and ω R ∼ π/t 0 .The instability is controlled by the amplification factor |A| 2 of the ergoregion at this frequency [20], i.e. ω I ∼ |A| 2 /t 0 .This follows from a very generic "bounce-and-amplify" argument, which was shown to accurately describe specific models [15,20].
A test scalar field Φ in such a geometry grows exponentially with time, Φ(t) ∼ e t/τ , with τ ≡ 1/ω I .The characteristic unstable modes can be computed numerically and agree well with bounce-and-amplify estimates [15].When t 0 t H ≡ GM /c 3 , these modes are well described by [15,[32][33][34] where ∆ = 1 − χ 2 , Ω(χ) is the object's angular velocity, β for a wave with angular number l [20,35], and p is an even (odd) integer for Dirichlet (Neumann) boundary conditions at the surface.Thus, the spacetime is unstable for ω R (ω R − mΩ) < 0 (i.e., in the superradiant regime [20]), on a timescale τ ≡ 1/ω I .Extensions of such calculations for electromagnetic fields with different boundary conditions at the surface shows that the simple bounceand-amplify estimate is a good upper bound to the instability timescale [36].Moreover, the amplification factors of scalars are orders of magnitude smaller than those of GWs [20,35].To be conservative, we therefore use the scalar results to estimate the background of GWs.
Equations ( 1) and ( 2) are valid for any object able to completely reflect the incoming radiation.In particular, if the object's surface sits at a constant (Boyer-Lindquist) radius where r + = GM (1 + ∆)/c 2 is the location of the (wouldbe) event horizon in these coordinates, then the travel time reads [6,8,25] Several different arguments about the magnitude of can be made.If quantum-gravity effects become important at Planck timescales t P = G/c 5 , it is natural to set = t P /t H ∼ 10 −39 − 10 −46 for stellar-mass to supermassive dark objects.These objects were dubbed ClePhOs in Ref. [15], and are impossible to rule out in practice via electromagnetic observations [8,15].

B. Modelling the interior
The above description effectively decouples the outside geometry from the inside, and is accurate when the flux across the surface vanishes.However, some models may have important transmittance.There are thus three different scenarios that need to be discussed in the general case: i.The object does not dissipate, and the light travel time inside the object is small (i.e., t interior ∼ t H t 0 ) .This situation describes most of the known models available in the literature.In such a case, the geometric center of the star effectively works as a perfectly reflecting mirror (i.e.ingoing radiation from one side exits on the other side with negligible delay), and the previous results (1)-(4) still apply.
ii.The object does not dissipate, and t interior t 0 .For these objects, a model of the interior is necessary, because even though all ingoing radiation will exit on the other side, delays/scattering due to the propagation in the interior will be important.To describe this case in a model-independent way, we assume that (1)-( 2) continue to apply, and we attempt to constrain the parameter t 0 directly, without assuming Eq. ( 4).
iii.The object dissipates radiation in its interior.In this case, the instability may be completely quenched if the absorption rate is large [34].For highly spinning objects, this requires at least 0.4% absorption rate for scalar fields, but up to 100% absorption rate for gravitational perturbations and almost maximal spins [20,37].While these numbers reduce to 0.1% for spins χ 0.7, they are still several orders of magnitude larger than achievable with viscosity from nuclear matter [34].
Based on the above arguments, we expect the following results to cover all relevant models of BH mimickers.

C. Evolution of the instability
The evolution of the object's mass and angular momentum, under energy and angular momentum losses, can be computed within the adiabatic approximation (because τ t H ) [38].The unstable mode is simply draining energy and angular momentum from the object, which we assume to have an equation of state such that = const during the evolution.From energy and angular momentum conservation, the evolution equations for each mode read where Ė0 encodes the initial preturbation of the (unstable) system.Since the instability is exponential, the overall evolution is insensitive to the precise value of these initial conditions.The equations above are valid for a monochromatic mode in a generic stationary and axisymmetric background.
The energy flux can be written as dE/df = Ė/ ḟ , where f = ω R /(2π) is the frequency associated with the mode.
From the evolution equations, we can evaluate Ṁ and J and, in turn, ωR .To leading order in the → 0 limit, we obtain ω R ∼ mΩ, and valid for any angular numbers l and m.In the same limit, the critical value of spin above which the ergoregion instability occurs reads [34] Thus, if χ(t = 0) > χ crit , the instability removes energy and angular momentum until superradiance is saturated, i.e. χ(t τ ) → χ crit .Note that the small spin value χ crit is compatible with the low measured spins of inspiralling compact objects detected via GWs so far [39].Since we are interested gravitational perturbations, when solving the evolution equations ( 5) and computing the energy flux (6) we only consider the dominant l = m = 2 mode and neglect highermodes.Note also that the above analysis assumes that the backreaction of these fields on the geometry is negligible.Our results indicate that those are always reliable approximations.

D. GW stochastic background
A population of GW sources too far and/or weak to be detected individually may still give rise to a "stochastic" background detectable by a network of interferometers, e.g. the LIGO/Virgo network, sensitive to frequency ranging from ∼ 10 to ∼ 100 Hz [40,41]; the Pulsar-Timing-Array experiments [42][43][44][45][46], which are already constraining backgrounds at frequencies ∼ 10 −9 − 10 −6 Hz; and the future LISA constellation [47], which will be sensitive to frequencies between 10 −6 Hz and ∼ 1 Hz.The background is produced by the incoherent superposition, at the detector, of the GW signals from all the unresolved sources in the population.The background can be characterized either by (i) its (dimensionless) energy spectrum (ρ gw being the background's energy density, f o the frequency measured at the detector and ρ c the critical density of the Universe at the present time), obtained by summing the energies emitted by all the unresolved sources in a given frequency bin [48]; or (ii) directly by the characteristic strain h c (f o ) observed in the detector, which can be obtained summing in quadrature (and binning in frequency) the strain amplitudes of all the unresolved sources [49].The two quantities are related by , where H 0 ≈ 68km/(s Mpc) is the Hubble rate.While these two ways of computing the background signal are equivalent [49], we have implemented both as a consistency check of our results.This also allows us to check that the number of sources contributing in each frequency bin is typically large as long as the bin size is 0.01 dex in the LISA band (which ensures that the number of sources contributing 99% of the signal in each bin ranges from thousands to millions).In the LIGO band, sources are even more numerous: frequency bins 0.01 dex yield 10 9 -10 14 sources contributing 99% of the signal in each bin.(These are mostly extragalactic sources, as Galactic ones give a negligible contribution to the background.)This in turn implies, in particular, that the background is expected to be smooth with that frequency binning [49].When computing the stochastic background of unstable exotic compact objects, the energy flux given by Eq. ( 6) is defined in the frequency range f ∈ [f min , f max ]; f max can be computed using Eq. ( 1) for a given initial mass and spin of the compact object, and f min is computed by solving the evolution equations ( 5) from the formation redshift of the compact object to the present time.
For the astrophysical populations of isolated BHs, we adopt the same models as in [48,53].For stellar-origin BHs we account for both Galactic and extragalactic BHs that form from the core collapse of massive ( 20M ) stars, by tracking the cosmic star formation history and the metallicity evolution of the Universe [54].We assume a uniform distribution for the initial spins with χ ∈ [0, 1] as the most optimistic and χ ∈ [0, 0.5] as the most pessimistic scenario.For the massive (∼ 10 4 − 10 7 M ) and supermassive (∼ 10 8 − 10 10 M ) BHs that emit respectively in the LISA and PTA bands, we adopt the semi-analytic galaxy-formation model of Ref. [55] (with later incremental improvements described in [56][57][58]), which follows the formation of these objects from their highredshift seeds and their growth by accretion and mergers.This growth is triggered in turn by the synergic co-evolution of the The black lines are the power-law integrated curves of [50], computed using noise PSDs for LISA with one year of observation time [47], LIGO's first observing runs (O1), LIGO at design sensitivity as described in [51], and an SKA-based pulsar timing array as described in [52].By definition, ρ stoch > 2 (ρ stoch = 2) when a power-law spectrum intersects (is tangent to) a power-law integrated curve.2, for an agnostic model for the compact-object (dissipationless) interior, where the light travel time t0 between the light ring and the surface is a free parameter.
BHs with their host galaxies, of which we evolve both the various baryonic components and the dark-matter halos.This model is optimistic since it predicts a spin distribution skewed towards large spins, at least at low masses.To include astrophysical uncertainties in our computation, we also consider models in between our most optimistic and most pessimistic assumptions as described in Ref. [48] (see Section III therein).

III. RESULTS
Our main results for the GW stochastic background from exotic compact objects are shown in Fig. 2 in the frequency bands relevant for LIGO/Virgo (left panel) and for LISA/an SKA-based pulsar timing arrays (right panel).The left panel suggests that the absence of a stochastic background in LIGO O1 already rules out our canonical model even for conservative spin distributions, while LIGO at design sensitivity will be able to rule out our canonical model even in more pessimistic scenarios than those assumed here, e.g. even if all BH-like objects had initial spin χ < 0.2.Similar results apply in the LISA band, whereas the stochastic signal is too small to be detectable by pulsar timing arrays, even in the SKA era.
The level of the stochastic background shown in Fig. 2 can also be understood with an approximate analytic calculation [53].The BH-mimicker mass fraction lost to GWs due to superradiance is F sr ∼ O(1%) [48,59].Because the signal spans about a decade in frequency [c.f.Eq. ( 1)], ∆ ln f ∼ 1, and Ω GW, sr = (1/ρ c )(dρ GW /d ln f ) ∼ F sr ρ BH /ρ c , with ρ GW and ρ BH the GW and BH-mimicker energy densities.In the mass range 10 4 − 10 7 M relevant for LISA, ρ BH ∼ O(10 4 )M /Mpc 3 , which gives Ω LISA GW, sr ∼ 10 −9 .To estimate the background in the LIGO band, note that the background from ordinary BH binaries is Ω GW, bin ∼ η GW F m ρ BH /ρ c , with η GW ∼ O(1%) the GW emission efficiency for BH binaries [60], and F m ∼ O(1%) [54] the fraction of stellar-mass BHs in merging binaries.This gives Ω GW, sr /Ω GW, bin ∼ F sr /(η GW F m ) ∼ 10 2 , and because the LIGO O1 results imply Ω GW, bin 10 −9 − 10 −8 [51,61], we obtain Ω LIGO GW, sr 10 −7 − 10 −6 .As shown in Fig. 3, LIGO/Virgo and LISA are also able to place model-independent constraints on the stochastic signal from exotic compact objects.At design sensivity, LIGO/Virgo can detect or rule out any model with t 0 < 10 8 t H , whereas LISA can go as far as t 0 < 10 6.5 t H .In other words, objects across which light takes 10 8 or less dynamical timescales to travel are ruled out.Finally, we note that although LIGO/Virgo rule out a wider range of the parameter space compared to LISA, it is still interesting to consider the constraints from both detectors since they probe different BH populations.

IV. DISCUSSION
Our results suggest that the current upper limits on the stochastic background from LIGO O1 already rule out the simplest models of BH mimickers at the Planck scale, setting the strongest constraints to date on exotic alternatives to BHs.The most relevant parameter for our analysis is the light travel time within the object, t 0 .LIGO/Virgo (LISA) can potentially rule out models where t 0 < 10 8 t H (t 0 < 10 6.5 t H ).These results are not significantly dependent on the astrophysical uncertainties of the extragalactic BH distributions, and even if exotic compact objects are produced all at low spin, the order of magnitude of our constraints is unaffected.
We conclude that all present and future models of exotic compact objects should either conform with t 0 10 8 t H or represent at most a fraction X of the compact-object population, the remaining being BHs.Since Eq. ( 8) scales linearly with X, Fig. 2 implies that the O1 upper limits impose X < 50% even if all BH-like objects are formed with low spins.At design sensitivity, these constraints could improve to X < 1%.
It might be possible to evade these constraints, by incorpo-rating some (exotic and still unclear) mechanism quenching the ergoregion instability, e.g.absorption rates several orders of magnitude larger than those of neutron stars [34].While such quenching mechanism might result in thermal or quasithermal electromagnetic radiation (which can be constrained by electromagnetic observations of BH candidates [62,63]), quantum-dressed BH mimickers might evade such constraints by trapping thermal energy in their interiors for very long timescales [8,15].Our results also imply that in the simplest models of nondissipative exotic ultracompact objects, GW echoes [6,8,25] can appear in the post-merger phase only after a delay time τ echo ∼ t 0 10 8 t H ≈ 10 4 [M/(20M )]s, which is much longer than what was claimed to be present in LIGO/Virgo data [64][65][66] (see also [67]).The latter claims would not be in tension with our bounds on the GW stochastic background only if one postulates exotic objects that are dissipative enough to absorb at least O(0.1)% of gravitational radiation (which is several orders of magnitude more than what typically achievable with nuclear matter).In this case, the stochastic background from a population of "echoing" merger remnants might still be detectable [68].

FIG. 2 .
FIG. 2. Extragalactic stochastic background for the canonical model in the LIGO/Virgo (left panel), LISA and PTA bands (right panel).The blue band brackets our population models (from the most pessimistic to the most optimistic, as explained in the main text).The background depends very weakly on as long as t0 ∼ tH | log | 10 5 tH , so here we show only the case t0 ∼ tH | log 10 −40 |.The black lines are the power-law integrated curves of[50], computed using noise PSDs for LISA with one year of observation time[47], LIGO's first observing runs (O1), LIGO at design sensitivity as described in[51], and an SKA-based pulsar timing array as described in[52].By definition, ρ stoch > 2 (ρ stoch = 2) when a power-law spectrum intersects (is tangent to) a power-law integrated curve.

t 0 FIG. 3 .
FIG.3.Same as in Fig.2, for an agnostic model for the compact-object (dissipationless) interior, where the light travel time t0 between the light ring and the surface is a free parameter.