Next Article in Journal
Magnetized Particle Motion around Black Holes in Conformal Gravity: Can Magnetic Interaction Mimic Spin of Black Holes?
Next Article in Special Issue
Editorial to the Special Issue “Probing New Physics with Black Holes”
Previous Article in Journal
Evolution of Quasiperiodic Structures in a Non-Ideal Hydrodynamic Description of Phase Transitions
Previous Article in Special Issue
Deformed General Relativity and Quantum Black Holes Interior
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Quantum Black Holes in the Sky

by
Jahed Abedi
1,2,†,
Niayesh Afshordi
3,4,5,*,†,
Naritaka Oshita
5,† and
Qingwen Wang
3,4,5,†
1
Max-Planck-Institut für Gravitationsphysik, D-30167 Hannover, Germany
2
Leibniz Universität Hannover, D-30167 Hannover, Germany
3
Department of Physics and Astronomy, University of Waterloo, Waterloo, ON N2L 3G1, Canada
4
Waterloo Centre for Astrophysics, University of Waterloo, Waterloo, ON N2L 3G1, Canada
5
Perimeter Institute For Theoretical Physics, 31 Caroline St N, Waterloo, ON N2L 2Y5, Canada
*
Author to whom correspondence should be addressed.
All authors have contributed equally to this work. The order of authors is alphabetical.
Universe 2020, 6(3), 43; https://doi.org/10.3390/universe6030043
Submission received: 16 January 2020 / Revised: 3 March 2020 / Accepted: 4 March 2020 / Published: 10 March 2020
(This article belongs to the Special Issue Probing New Physics with Black Holes)

Abstract

:
Black Holes are possibly the most enigmatic objects in our universe. From their detection in gravitational waves upon their mergers, to their snapshot eating at the centres of galaxies, black hole astrophysics has undergone an observational renaissance in the past four years. Nevertheless, they remain active playgrounds for strong gravity and quantum effects, where novel aspects of the elusive theory of quantum gravity may be hard at work. In this review article, we provide an overview of the strong motivations for why “Quantum Black Holes” may be radically different from their classical counterparts in Einstein’s General Relativity. We then discuss the observational signatures of quantum black holes, focusing on gravitational wave echoes as smoking guns for quantum horizons (or exotic compact objects), which have led to significant recent excitement and activity. We review the theoretical underpinning of gravitational wave echoes and critically examine the seemingly contradictory observational claims regarding their (non-)existence. Finally, we discuss the future theoretical and observational landscape for unraveling the “Quantum Black Holes in the Sky”.

1. Introduction

Black holes (BHs) are very interesting “stars” in the universe where both strong gravity and macroscopic quantum behavior are expected to coexist. Classical BHs in General Relativity (GR) have been thought to have only three hairs, i.e., mass, angular momentum, and charge, making observational predictions for BHs relatively easy [1,2] (compared to other astrophysical compact objects). For astrophysical BHs, due to the effect of ambient plasma, this charge is vanishingly small, leaving us with effectively two hairs for isolated black holes, with small accretion rates. In other words, finding conclusive deviations from standard predictions of these 2-parameter models, may be interpreted as fingerprints of a quantum theory of gravity or other possible deviations from GR. For example, the quasinormal modes (QNMs) of spinning BHs, which have been widely studied over the past few decades (a subject often referred to as BH spectroscopy), only depend on the mass and spin of the Kerr BH (e.g., [3]). The ringdown of the perturbations of the BH is regarded as a superposition of these QNMs, and thus can be used to test the accuracy of GR predictions and no-hair theorem (e.g., see [4]). As a result, precise detection of QNMs from the ringdown phase (from BH mergers or formation) in gravitational wave (GW) observations may enable us to test the classical and quantum modifications to GR (e.g., [5,6,7]).
A concrete path towards this goal is paved through the study of “GW echoes”, a smoking gun for near-horizon modifications of GR which are motivated from the resolutions of the proposed resolutions to the BH information paradox and dark energy problems [8,9]. The list of these models include wormholes [10], gravastars [11], fuzzballs [12], 2-2 holes [13], Aether Holes [9], Firewalls [8] and the Planckian correction in the dispersion relation of gravitational field [14,15].
The possibility of observing GW echoes was first proposed shortly after the first detection of GWs by LIGO [10,16,17], which led to several observational searches [18,19,20,21,22,23,24,25,26,27,28]. Tentative evidence for and/or detection of these echoes can be seen in the results reported by different groups [18,19,20,21,22,23,24,25] from O1 and O2 LIGO observations of binary BH and neutron star mergers, but the origin and the statistical significance of these signals remain controversial [23,25,26,27,28], motivating further investigation.
Given their uncertain theoretical and observational status, GW echoes are gathering much attention from those who are interested in the observational signatures of quantum gravity, and the field remains full of excitement, controversy and confusion. In this review article, we aim to bring some clarity to this situation, from its background, to its current status, and into its future outlook.
The review article is organized as follows: in the next section, we provide builds the motivation to investigate the quantum signatures from BHs. In Section 3, we discuss theoretical models of quantum BHs, starting from the BH information loss paradox, and then its proposed physical resolutions that lead to observable signatures. In Section 4, we review how to predict the GW echoes from spinning BHs based on the Chandrasekhar-Detweiler (CD) equation, and also review the Boltzmann reflectivity model [15,29] for quantum black holes. Section 5 is devoted to the echo searches, where we summarize positive, negative, and mixed reported outcomes, and attempt to provide a balanced and unified census. In Section 6, we discuss the future prospects for advancement in theoretical and observational studies of quantum black holes, while Section 7 concludes the review article.
Throughout the article, we use the following notations:
SymbolDescription
aspin parameter
a ¯ non-dimensional spin parameter ( a / ( G M ) )
cspeed of light
ħPlanck constant
k B Boltzmann constant
Ggravitational constant
M Pl Planck mass
E Pl Planck energy
l Pl Planck length
Mmass of a balck hole or exotic compact object
M solar mass ( 1.988 × 10 30 kg)
r g Schwarzschild radius
T H Hawking temperature
Furthermore, unless noted otherwise, we use the natural Planck units with ħ = c = 1 = G = 1 .

2. Invitation

2.1. Classical BHs

Schwarzschild and Kerr spacetimes are solutions of general relativity giving the spacetime configurations of BHs, which are the most dense objects in our universe. In some regions, matter accumulates and attracts more matter with its gravity which is classically always attractive. At the end, the force is so strong that even the light cannot escape from those regions, where then BHs form. The first and most important feature (the definition of BHs) is the formation of horizon. Inside the (event) horizon, all the light cones are directed into the singularity, and nothing can escape, unless it could travel faster than the speed of light. Therefore, horizons stand as the causal boundaries of BHs in Einstein’s theory of Relativity.
Realistic BHs in the sky have different hairs (mass, spin, and charge), and their dynamics share more complicated structure, thus, have different kinds of horizons. To list some of them, event horizons are defined as the boundaries where no light can escape to the infinite future. However, for a dynamically evolving BH, event horizons are teleological, i.e., we cannot predict them until we have the entire history of the spacetime. Apparent horizons, however, are predictable at a specific time without knowing the future. Any surface has two null normal vectors and if expansion of both of them are negative, the surface is called “trapped”. Apparent horizons are the outermost of all the trapped surfaces, which is why they are also known as the “marginally outer trapped surface”.
Here is a simple example to distinguish these two horizons—we start with a Schwarzschild BH at time t 1 , now the event and apparent horizons coincide at the Schwarzschild radius. We throw a spherical null shell into the BH and let it cool down at t 2 . This process is perfectly described by Vaidya metric [30]. The apparent horizon changes immediately when the shell falls into the BH but the event horizon starts to expand earlier, even before the shell reaches it. It is because that after throwing the shell, the gravity of BH increases. Thus it is harder for light to escape from the BH to infinity. In other words, particles might be doomed to fall into a singularity, even before they had a chance to meet the infalling gravitating matter that is responsible for their fate. Therefore, the event horizon is modified earlier than the apparent horizon. While this result is counter-intuitive, it is a result of the formal definition of the event horizons, which requires the information about the entire history of spacetime, in particular, the future!
Beyond the horizons, another intriguing trait of BHs is the curvature singularity, which sits at the centre of the BHs. Horizons can also be singular, but usually only coordinate singularities and (in classical General Relativity) removable by changing to a proper coordinate system. However, the singularities inside the BHs are where the general relativity breaks down and so far we do not have any good physics to describe them. We cannot chase the information lost into these singularities (using standard physics), which leads to the information paradox (more on this later).
Back in November 1784, John Michell, an English clergyman, advanced the idea that light might not be able to escape from a very massive object (at a fixed density). For example, light cannot escape from the surface of a star with the density of the sun, if it was 500 times bigger than the sun. Albert Einstein, later in 1915, developed general relativity. Soon after this, Karl Schwarzschild solved the Einstein vacuum field equation under spherical symmetry with a singular mass at the center, which was the first solution for BHs, the Schwarzschild metric.
While the 20th century saw a golden age of general relativity with blooming of dozens of different BH solutions, the existence of BHs was not directly confirmed until one century later in 2015. LIGO-Virgo collaboration reported unprecedented detection of GWs from the binary BH merger events [31,32,33,34,35,36,37,38]. Numerical relativity is consistent with LIGO data at least up to quite near the horizon range. However, the detection has not confirmed the existence of the horizons. We will discuss in this article how the detection opens a window for searching for quantum nature of the BHs beyond the general relativity.

2.1.1. Schwarzschild Spacetime

The Schwarzschild spacetime was the first exact solution in the Einstein theory of general relativity. It models a static gravitational field outside a mass which has spherical symmetry, zero charge and rotation. Karl Schwarzschild found this solution in 1915, and four months later, Johannes Droste published a more concrete study on this independently. The metric in the Schwarzschild coordinate is:
d s 2 = 1 2 M r d t 2 + 1 2 M r 1 d r 2 + r 2 d Ω 2 ,
where M is the mass of the centre object, 2 M is Schwarzschild radius and d Ω 2 = d θ 2 + sin 2 θ d ϕ 2 is the metric on a 2-sphere. The metric describes gravitational field outside any spherical object without charges. If the radius of the central object is smaller than the Schwarzschild radius, the object is then too dense to be stable, and will go through a gravitational collapse and form the Schwarzschild BH.
Later in 1923, G.D.Birkhoff proved that any spherically symmetric solution of the vacuum Einstein field equation must be static and asymptotically flat. Hence, Schwarzschild metric is the only solution in that case. For any static solution, the event horizon always coincides with the apparent horizon at r = 2 M . In general relativity, Schwarzschild metric is singular at the horizon but, as stated above, this is only a coordinate artifact. That is to say, a free falling observer feels no drama going through the horizon. It takes the observer a finite amount of proper time but infinite coordinate time. Particularly, we can remove the singularity by a proper coordinate transformation. In contrast, the origin r = 0 is the intrinsic curvature singularity. The scalar curvature is infinite and the general relativity is no longer valid at this point.

2.1.2. Kerr Spacetime

The Kerr spacetime [39], discovered by Roy Kerr, is a realistic generalization of the Schwarzschild spacetime. It describes the gravitational field of an empty spacetime outside a rotating object. The spacetime is stationary and has axial symmetry. The metric in the Boyer-Lindquist coordinate is:
d s 2 = 1 2 M r ρ 2 d t 2 + ρ 2 Δ d r 2 + ρ 2 d θ 2 + r 2 + a 2 + 2 M r a 2 ρ 2 sin 2 θ sin 2 θ d ϕ 2 4 M r a sin 2 θ ρ 2 d t d ϕ ,
= ρ 2 Δ Σ d t 2 + Σ ρ 2 sin 2 θ ( d ϕ ω d t ) 2 + ρ 2 Δ d r 2 + ρ 2 d θ 2 ,
where a = J / M , ρ 2 = r 2 + a 2 cos 2 θ , Δ = r 2 2 M r + a 2 , Σ = ( r 2 + a 2 ) 2 a 2 Δ sin 2 θ and ω = g t ϕ g ϕ ϕ = 2 M a r Σ . The Cartesian coordinates can be defined as:
x = r 2 + a 2 sin θ cos ϕ , y = r 2 + a 2 sin θ sin ϕ , z = r cos θ .
There are two singularities easily reading from the coordinate where the g r r and g t t vanish. The first one gives r ± = M ± M 2 a 2 corresponding to the horizon analog to the Schwarzschild metric. The larger root r + = M + M 2 a 2 is the event horizon, while the other root is inner apparent horizon.
The second singularity is related to an interesting effect in the Kerr spacetime called frame-dragging effect: When reaching close to the Kerr BHs, the observers even with zero angular momentum (ZAMOs) will co-rotate with the BHs because of the swirling of spacetime from the rotating body. We assume that u α is the four-velocity of ZAMOs, and from the conservation of angular momentum g ϕ t t ˙ + g ϕ ϕ ϕ ˙ = 0 , where an overdot is differentiation with respect to the proper time of the observers τ . Thus, d ϕ d t = g t ϕ g ϕ ϕ . Because of this frame-dragging effect, there is a region of spacetime where static observers cannot exist, no matter how much external force is applied. This region is known as the “ergosphere” r M + M 2 a 2 cos 2 θ . The rotation also leads to another interesting feature, called “superradiance”. That is, we can extract energy from scattering waves off the Kerr BHs. The exact formalism of superradiance is defined and discussed in Section 4.2.
Finally, the Kerr spacetime also possesses a curvature singularity at the origin ρ 2 = r 2 + a 2 cos 2 θ . However, in contrast to Schwarzschild case, this singularity can be avoided since it is a ring at r=0 and θ = π / 2 , where z=0 and x 2 + y 2 = a 2 . In principle, observers can go through the ring without hitting the singularity. However, it is widely believed that the inner horizon, r in Kerr spacetime is subject to an instability which would dim the analytic extension of Kerr metric beyond r unphysical [40].

2.1.3. Blue-Shift Near Horizon

As shown in the metric, different observers have different proper times. Hence, in the general relativity, the clocks at a gravitational field tick in a different speed in a different spacetime point. This is the blue(red)-shift effect, and it is extremely strong close to the dense object, especially near horizon.
Assuming static clocks in the Schwarzschild spacetime d s 2 = d τ 2 = ( 1 2 M / r o ) d t 2 , where τ is the proper (clock) time of an observer at distance r o . Hence, t is the proper time of an observer at infinity. The shifted wavelength λ o measured by observers at r o compared to observers at infinite is:
λ o λ = d τ d t = 1 2 M r o 1 / 2 .

2.1.4. Thermodynamics of Semi-Classical BH

Jacob Bekenstein and Stephen Hawking first proposed that the entropy of BHs is related to the area of their event horizons divided by the Planck area [41,42,43,44,45]. Furthermore, in 1974, Stephen Hawking showed that rather than being totally black, BHs emit thermal radiation at the Hawking temperature, T H = κ 2 π , where κ is the surface gravity at the horizon [46,47,48]. This then lead to the celebrated Bekenstein-Hawking entropy formula S BH = A 4 [44,45], where A is the area of the event horizon. However, the nature of microstates of BHs that are enumerated by this entropy remains so far unknown. String theory associates it with higher dimensional fuzzball solutions, as discussed later in Section 3.5. Loop quantum gravity relates the quantum geometries of the horizon to the microstates [49]. Both these approaches can give the right Bekenstein-Hawking entropy, given specific assumptions and idealizations.
Interestingly, not only the entropy exists for the BHs, but also Brandon Carter, Stephen Hawking and James Bardeen [50] discovered the four laws of BH thermality analogous to the four laws of thermodynamics. The latter is presented in the parentheses.
  • The zeroth law: A stationary BH has constant surface gravity κ . (A thermal equilibrium system has a constant temperature T H .)
  • The first law: A small change of mass for a stationary BH is related to the changes in the horizon area A, the angular momentum J, and the electric charge Q: d M = κ 8 π d A + Ω d J + Φ d Q , where Ω is the angular velocity and Φ is the electrostatic potential (Energy conservation: d E = T d S P d V μ d N ).
  • The second law: The area of event horizon A never decreases in general relativity (The entropy of isolated systems never decreases).
  • The third law: BHs with a zero surface gravity cannot be achieved (Matter in a zero temperature cannot be reached).

2.2. Membrane Paradigm

As mentioned above, in classical general relativity, freely falling observers experience no drama as they cross the event BH horizons, at least not until they reach the singularity inside the BH. However, to a distant and static observer outside a BH, any infalling objects are frozen at the horizons due to the blue-shift effect. Hence, the BH interior can be regarded as an irrelevant region for the static observers. Based on this complementary picture near horizon, in 1986, Kip S. Thorne, Richard H. Price and Douglas A. Macdonald published the idea of membrane paradigm [51]. They use a classically radiating membrane to model the BHs, which is motivated as a useful tool to study physics outside BHs without involving any obscure behavior within BH interior.
Let us introduce a spherical membrane located infinitesimally outside the Schwarzschild radius, a.k.a. the stretched horizon. When the membrane is sufficiently thin, one can use the Israel junction condition to nicely embed the membrane in the Schwarzschild spacetime. The condition is:
( K ( + ) f a b K a b ( + ) ) ( K ( ) f a b K a b ( ) ) = 8 π T a b ,
where f a b is the induced metric of the membrane, K a b ( ± ) is the extrinsic curvatures on its two sides, and T a b is its stress tensor. The infalling observer will cross the horizon and enter the BH interior without possibility of seeing the membrane. However the static observer outside the BH can remove irrelevant interior region from the remaining spacetime with a membrane. Assuming reflection symmetry K a b + = K a b , the Israel junction condition on the membrane becomes:
K ( + ) f a b K a b ( + ) = 4 π T a b ,
where the stress tensor T a b is no longer zero but has contributions from the extrinsic curvature on the membrane. Rewriting the left hand side of Equation (7)), one can obtain the following relation:
T b a = 1 4 π σ b a + δ b a θ 2 + κ ,
where σ b a is the shear, θ is the expansion, and κ is the surface gravity at the horizon. From the analogy with the energy momentum tensor of the 2-dimensional compressible fluid, one can read that the shear viscosity η and bulk viscosity ζ are given by η = 1 / ( 16 π ) and ζ = 1 / ( 16 π ) , respectively. The negativity of the bulk viscosity implies gravitational instability for the expansion or compression at the horizon. The viscosity at the membrane lead to the thermal dissipation of infalling gravitational waves into the horizon and in this sense, the causality of a classical BH is dual to the viscosity of a 2-dimensional fluid at the stretched horizon. This was the earliest example of fluid-gravity correspondence, which is now active area of research.
Modifying Einstein gravity which revises the structure of BHs can provide a modified structure of the thin-shell membrane. For example, by adapting the transport properties of the membrane fluid, we can investigate various models of quantum BHs. As an example, we provide a simple idea [15] which relates the reflectivity of the “horizon” of quantum BHs to viscosity in the context of membrane paradigm. We start by perturbing the Schwarzschild spacetime, whose metric is g μ ν Sch . Within Regge-Wheeler formalism [52], the axial axisymmetric perturbation g μ ν = g μ ν Sch ( r ) + δ g μ ν ( r , θ , t ) take teh form:
δ g t ϕ = ϵ e i ω t h 0 ( r ) y ( θ ) ,
δ g r ϕ = ϵ e i ω t h 1 ( r ) y ( θ ) ,
where other δ g μ ν components vanish, and ϵ 1 controls the order of perturbation. The membrane stands at r = r 0 + ϵ R ( t , θ ) , where r 0 is its unperturbed position. We apply the Israel junction conditions K a b K f a b = 4 π T a b to Brown-York stress tensor as defined in [53]. The indexes μ , ν run over ( t , r , θ , ϕ ) in the 4d spacetime, while a , b run over ( t , θ , ϕ ) on the 3d membrane. We further assume that T a b is the energy stress tensor of a viscous fluid:
T a b = [ ρ 0 + ϵ ρ 1 ( t , θ ) ] u a u b + [ p 0 + ϵ p 1 ( t , θ ) ζ Θ ] γ a b 2 η σ a b ,
σ a b = 1 2 ( u a ; c γ b c + u b ; c γ a c Θ γ a b ) ,
γ a b h a b + u a u b , Θ u ; a a ,
where ρ 0 and p 0 ( ρ 1 and p 1 ) are background (perturbation on) membrane density and pressure, and u a , η and ζ are fluid velocity, shear viscosity, and bulk viscosity, respectively. Plugging Equations (9)–(13) into the the Israel junction condition, we find in the zeroth order in ϵ :
ρ 0 ( r 0 ) = f ( r 0 ) 4 π r 0 ,
p 0 ( r 0 ) = f ( r 0 ) ( g ( r 0 ) + r 0 g ( r 0 ) ) 8 π r 0 g ( r 0 ) ,
where g ( r 0 ) = ( 1 2 M / r 0 ) 1 / 2 and f ( r 0 ) = 1 2 M / r 0 . Assuming u ϕ = 0 , equation of θ ϕ component gives in next order of ϵ :
ω h 1 ( r ) = 8 i π η [ h 1 ( r ) + ( r r g ) h 1 ( r ) ] .
We can further use ψ ω = 1 r 1 2 M r h 1 ( r ) and the tortoise coordinate x = r + 2 M log [ r / ( 2 M ) 1 ] to rewrite Equation (16) as:
ω ψ ω = 16 i π η ψ ω x .
For the classical BHs with a purely ingoing boundary condition ψ ω e i ω x at the horizon, Equation (17) gives η = 1 16 π , which is consistent with the standard membrane paradigm. If instead we assume there is no longer horizon but a reflective surface with ψ ω = A o u t e i ω x + A i n e i ω x , Equation (17) gives:
A out A in = 1 16 π η 1 + 16 π η e 2 i ω x .
Which relates the reflectivity of the membrane to the viscosity of the surface fluid.

2.3. Dawn of Gravitational Wave Astronomy

From 2015 onwards, the LIGO/Virgo collaboration reported unprecedented GW observations from binary BH merger events [31,32,33,34,35,36,37,38]. It is the first time that humankind can detect GWs after one century of the Einstein’s general theory of gravity. In 2017, Rainer Weiss, Kip Thorne and Barry C. Barish won the Nobel Prize in Physics “for decisive contributions to the LIGO detector and the observation of Gravitational Waves”.
The first and most prominent binary BH merger signal seen by LIGO, GW150914, matches well with predictions of numerical relativity simulations that settle into Kerr metric, but contrary to original claims, it could not confirm the existence of the event horizons [10]. However, it opened a new front to test general relativity in strong gravity regime and Kerr-like spacetimes (e.g., Quantum BHs) from modified gravity, which is the main topic of this review article.
This is the dawn of GW astronomy, and we stand at the threshold of a new age. We are detecting even more compact binary merger events with a better sensitivity from the O3 run of LIGO/Virgo. Future experiments such as Einstein Telescope, Cosmic Explorer, and LISA are expected to improve this by orders of magnitude. More studies on the echo-emission mechanism as well as observational strategies will be crucial for taking advantage of these new observations, to shed light on the nature of quantum BHs. It is our point of view that the best bet is on a sustained synergy between theory and observation, relying on well-motivated theoretical models (such as the Boltzmann reflectivity, aether holes, 2-2 holes, or fuzzballs, discussed in this review) to provide concrete templates for data analysis, which in turn could be used to pin down the correct theory underlying quantum BHs. With some luck, this has the potential to revolutionize our understanding of fundamental physics and quantum gravity.

2.4. Quantum Gravity and Equivalence Principle

The Einstein’s general theory of relativity is classical. However, in the Einstein field equation G μ ν = 8 π G T μ ν , the classical spacetime geometry is related to stress energy tensor of quantum matter. For decades, scientists tried to reconcile this inconsistency by embedding general relativity (or its generalizations) within some quantum mechanical framework, i.e., quantum gravity.
Conventional approach to quantizing Einstein gravity fails because it is not renormalizable. This implies that making predictions for observables, such as scattering cross-sections, requires knowledge of infinitely many parameters at high energies, leading to loss of predictivity. In the modern language, general relativity could at best be an effective field theory, and requires UV-completion beyond a cutoff near (or below) Planck energy (e.g., [54]).
Most proposals for this UV-completion involve replacing spacetime geometry with a more fundamental degree of freedom, such as strings (string theory) [55], discrete spins (loop quantum gravity) [56], spacetime atoms (causal sets) [57], or tetra-hydra (causal dynamical triangulation) [58]. More exotic possibilities include Asymptotic Safety [59], Quadratic Gravity [60], and Fakeon approach [61] that introduce a non-perturbative or non-traditional quantization schemes for 4d geometry. Yet another possibility is to modify the symmetry structure of General Relativity in the UV, as is proposed in Lorentz-violating (or Horava-Lifshitz) quantum gravity [62].
While proponents of these various proposals (with varying degrees of popularity) claimed limited success in empirical explanations of some natural phenomena, it should be fair to say that none can objectively pass muster of predicitivity. As such, for now, the greatest successes of these proposals remain in the realm of Mathematics.
Due to this lack of concrete predictivity, the EFT estimates (discussed above) are instead commonly used to argue that the quantum gravitational effects should only show up at Planck scale ∼ 10 35 m or 10 28 eV, which is far from anything accessible by current experiments. However, such arguments miss the possibilities of non-perturbative effects (such as phase transitions) which depend on a more comprehensive understanding of the full phase space of the specific quantum gravity proposal.
For example, it has been shown that the non-perturbative quantum gravitational effects may lead to Planck-scale modifications of the classical BH horizons [63]. Proposed models, such as gravastars [11], fuzzballs [12,64,65,66,67], aether BHs [9], and firewalls [68,69] among others [70,71,72,73] all drastically alter the standard structure of the BH stretched horizons with a non-classical surface. Soon after the first reported detection of gravitational waves, Ref. [10] discerned that Planck-length structure modification around horizons leads to a similar waveform as in classical GR, but followed by later repeating signal—echoes—in the ringdown from the reflective surface that replaces the classical horizon. This discovery equals a new road leading to Rome—quantum nature of gravity—and has sparked off a novel area of modeling and searching for signatures from Quantum BHs. The next section will discuss the quantum theories of BH models and possible road maps to probe them, inspired by the detection of binary BH merger events in gravitational waves.

3. Quantum BHs

3.1. Evaporation of BHs and the Information Paradox

It was already recognized by Stephen Hawking in the 1970s that the evaporation of a BH leads to an apparent breakdown of the unitarity of quantum mechanics. Here, we will briefly review this problem, which is known as the BH information loss paradox [74]. In the context of quantum field theory in curved spacetime, the energy flux out of a BH horizon is obtained by specifying a proper vacuum state and fixing the (classical) background spacetime. However, a radiating BH must lose its mass in time, and so fixing the background is valid only for a much shorter timescale than the evaporation timescale. One can roughly estimate the lifetime of a BH as follows: The energy expectation value of a Hawking particle is of the order of the Hawking temperature T H ( 8 π M ) 1 , which would be emitted over the timescale of t M . Then we can estimate the luminosity of the BH as:
d M d t T H M ( M 2 ) 1 ,
and this gives its lifetime t life :
t life M 3 .
To be consistent with the result of a more rigorous calculation (see e.g., [75]), we need a factor of about 10 5 in (20):
t life 10 5 M 3 10 75 M M 3 [ s ] ,
which is much much longer than the cosmic age of 4 × 10 17 [ s ] for astrophysical BHs whose mass are M . It may be true that BHs evaporate due to the Hawking radiation, at least, until reaching the Planck mass. However, the gravitational curvature near the horizon eventually reaches the Planckian scale and the classical picture of background gravitational field would break down. As such, the possibility of leaving a “remnant” after the evaporation has been discussed (see e.g., [76,77,78,79]), but the most natural possibility would be that only Hawking radiation is left after the completion of the BH evaporation.
If the Hawking evaporation just leaves the “thermal" radiation afterwards, one can immediately understand why the evaporation process is paradoxical. Let us suppose that a pure quantum state collapses into a BH and it radiates Hawking quanta until the BH evaporates. If the final state is a thermal mixed state, the evaporation is a process which transforms a pure to mixed state. Therefore, if the final state of any BH is a completely thermal state, one can say that the evaporation process is a non-unitary process. The information loss paradox can be also explained from the geometric aspect using the Penrose diagram. In quantum mechanics, the time-evolution of a quantum state is described by a unitary operator, U ^ , that maps an initial quantum state | in on a past Cauchy surface Σ i into a final quantum state | f on a future Cauchy surface Σ f . Since the unitary operator gives a reversible process, one can also obtain the initial state from the final state as:
| in = U ^ | f .
Although this is true in a flat space, the argument is very controversial in the existence of an evaporating BH. Assuming a gravitational collapse forms a horizon and singularity, then it eventually evaporates, leaving behind a thermal radiation, the Penrose diagram describing the whole process is given by Figure 1. Let us consider three quantum states: an initial quantum state | in on Σ i , an intermediate quantum state | mid on Σ m , and a final state | f on Σ f , where Σ i , Σ m , and Σ f are the Cauchy surfaces and Σ m intersects the future horizon H + and so one can split it into the exterior and interior regions as Σ m Σ ext Σ int (see Figure 1).
The final quantum state | f is determined by information on the exterior part of the intermediate Cauchy surface Σ ext rather than that on the whole intermediate Cauchy surface Σ m , which leads to the information loss paradox. To see this in more detail, let us consider an initial pure quantum state:
| in = i c i in | ψ i ,
where c i in is an initial vector in the Hilbert space. The intermediate state is still a pure state due to the unitary evolution of | in :
| mid = U ^ | in = i , j c i , j | ψ i int | ψ j ext ,
the time-evolution from Σ m to Σ f is non-unitary, provided that the final state on Σ f is obtained by the unitary evolution of the exterior intermediate state. The density matrix of the exterior intermediate state, denoted by ρ ^ ext , is obtained by tracing over all the internal basis states:
ρ ^ ext = k ψ k | int | mid mid | ψ k int = k , j , j c k , j c k , j | ψ j ext ψ j | ext .
The resulting density matrix, (25), is independent of the interior orthogonal basis | ψ j int due to the tracing operation. Therefore, the loss of the interior information results in a non-unitary evolution and an initial quantum state evolves to a mixed state after the BH evaporation.

3.2. BH Complementarity

The BH complementarity has been one of the leading proposals for the retrieval of BH information, which was first put forth by by Susskind, Thorlacius, and Uglum [80]. According to a distant observer, due to the infinite redshift at a BH horizon, the Hawking radiation involves modes of transplanckian frequency whose energy can be arbitrarily large in the vicinity of the horizon. In the BH complementarity proposal, the energetic modes form the membrane, which can absorb, thermalize, and reemit information, on the BH horizon. They argue that such a picture regarding the retrieval of BH information by the stretched horizon is consistent with the following three plausible postulates:
Postulate 1
(unitarity)According to a distant observer, the formation of a BH and the evaporation process can be described by the standard quantum theory. There exists a unitary S-matrix which describes a process from infalling matter to outgoing non-thermal radiation.
Postulate 2
(semi-classical equations)Outside the stretched horizon of a massive BH, physics can be approximately described by a set of semi-classical field equations.
Postulate 3
(degrees of freedom)For a distant observer, the number of microscopic states of a BH can be estimated by exp S ( M ) , where the exponent S ( M ) is the Beksntein-Hawking entropy.
On the other hand, it has been presumed that a freely infalling observer would not observe anything special when passing through the horizon due to the equivalence principle. In this sense, there are two totally different and seemingly inconsistent scenarios that co-exist in the BH complementarity. However, the contradiction arises only when attempting to compare the experiments performed inside and outside horizon, which might be impossible due to a backreaction of the high-energy modes near the stretched horizon [81].

3.3. Firewalls

In 2012, Almheiri, Marolf, Polchinski and Sully (AMPS) argued [69] that the Postulates 1–3 in the BH complementarity and the Equivalence principle of GR are mutually inconsistent for an old BH [82,83,84], provided that the monogamy of entanglement is satisfied. Then they argued that the “most conservative” resolution is a violation of the equivalence principle near the BH and its horizon should be replaced by high-energetic quanta, so-called “firewall”, to avoid the inconsistency. Before introducing the original firewall argument in more detail, let us review a theorem in quantum information theory, the monogamy of entanglement. Let us consider three independent quantum systems, A, B, and C. The strong subadditivity relation of entropy is given by:
S A B + S B C S B + S A B C .
If A and B is fully entangled, we have:
S A B = 0 and S A B C = S C .
Then the strong subadditivity relation reduces to:
S B + S C S B C 0 .
Since the left hand side in (28) is the mutual information of B and C, denoted by I B C , and it is a non-negative quantity, (28) reduces to
I B C = S B + S C S B C = 0 ,
which means that the quantum system B cannot fully correlate with C when B and A are fully entangled mutually. Therefore, any quantum system cannot fully entangle with other two quantum systems simultaneously. This is the monogamy of entanglement that is an essential theorem in the firewall argument.
Let us consider an old BH, whose origin is a gravitational collapse of a pure state, with early Hawking particles A, late Hawking particle B, and infalling particle inside the horizon C. In order for the final state of the BH to be pure state, A and B should be fully entangled mutually, that is a necessary condition for the Postulate 1. On the other hand, created pair particles, B and C, are also fully entangled according to the quantum field theory in classical background (Postulate 2). That is, imposing the Postulate 1 and 2 inevitably results in that B is fully and simultaneously entangled with both A and C, which obviously contradicts with the monogamy of entanglement. To avoid this contradiction, AMPS argued that there is no interior of BHs and the horizons should be replaced by energetic boundaries that the entanglement of Hawking pairs are broken. They called these boundaries “firewalls”. According to this proposal, any object falling into a BH would burn up at the firewall, which contradicts the equivalence principle (in vacuum) and replaces the BH complementarity proposal. Although there are some updates of this proposal, based on ER = EPR conjecture [8,85,86,87,88], backreaction due to gravitational schockwaves [89], and quantum decohenrence of Hawking pair due to the interior tidal force [90]), they do remain speculative, and at the level of toy models. However, on general grounds, if quantum effects lead to such an energetic wall at the stretched horizon, it could contribute to the reflectivity of BH which may be observable by merger events leading to the formation of BHs.

3.4. Gravastars

The gravitational vacuum condensate star (gravastar) was proposed as a final state of gravitational collapse by Mazur and Mottola [11]. According to the proposal, the resulting state of gravitational collapse is a cold compact object whose interior is a de Sitter condensate, which is separated from the outside black hole spacetime by a null surface. In this state, there is no singularity (with the exception of the null boundary) and no event horizon, which avoids the BH information loss paradox. Such gravitational condensation could be caused by quantum backreaction at the Schwarzschild horizon r = r g even for an arbitrarily large-mass collapsing object. One might wonder why the backreaction can lead to such a drastic effect for any mass since the tidal force which acts on an infalling test body can be arbitrarily weak for an arbitrarily large mass at the Schwarzschild radius. The argument is that considering a photon with asymptotic frequency ω near the Schwarzschild radius, the (infinite) blue-shift effect by which the local energy is enhanced as ħ ω / 1 r g / r , could lead to a drastic effect at the Schwarzschild radius. This is unavoidable since any object is immersed in quantum vacuum fluctuations and virtual particles always exist around them. From this argument, the gravitational condensation has been expected to take place at the final stage of gravitational collapse. The authors in [11] also estimate the entropy on the surface of gravastar by starting with a simplified vacuum condenstate model which consists of three different equations of state:
0 r < r 1 , ρ = p ,
r 1 < r < r 1 + δ r , ρ = p ,
r 1 + δ r < r , ρ = p = 0 ,
where r 1 is the radius of interior region and δ r is the thickness of the thin-shell of the gravastar. Then the obtained entropy of the shell was found out to be S 10 57 g k B M / M 3 / 2 , where g is a dimensionless constant. Recently, the derivation of gravastar-like configuration was performed by Carballo-Rubio [91]. He derived the semi-classical Tolman-Oppenheimer-Volkoff (TOV) equation by taking into account the polarization of quantum vacuum and solved it to obtain the exact solution of an equilibrium stellar configuration. It also has its de Sitter interior and thin-shell near the Schwarzschild radius, which is consistent with the original gravastar proposal [11].
From the observational point of view, the shadows of a gravastar was investigated in [92] where they argue the shadows of a BH and gravastar could be distinguishable. In addition, tests of gravastar with GW observations were discussed in e.g., [16,21,93].

3.5. Fuzzballs

Samir D. Mathur has proposed fuzzballs [65] as description of true microstates of the quantum BHs from string theory. A fuzzball state has the BH mass inside a horizon-sized region and a smooth (but higher-dimensional) geometry. Here are some crucial features of the conjecture:
  • Different fuzzball geometries represent different microstates of the quantum BH—fuzzball. Application the AdS/CFT duality [94] suggests that the counting of the microstates is consistent with the Bekenstein-Hawking entropy.
  • Fuzzballs do not possess horizons. Instead, they end with smooth “caps” near where the horizons would have been. Every microstate has almost the same geometry outside the would-be horizon matching the classical BH picture for the outside observers. However, the microstates differ from each other near the would-be horizons.
  • Fuzzball solves the information paradox by removing the horizon and singularity. The horizon is replaced by fuzzy matter and no longer vacuum. The particles created near the would-be horizon now have access to the information of fuzzball interior. Moreover, the higher-dimensional spacetime ends smoothly around the would-be horizon and is singularity-free. The infalling particles at the low frequencies interact with the “fuzz” for a relatively long time scale, while high frequency ones excite the microstates and lose their energy the same as in the classical BHs case. Hence, the traditional horizons only show up effectively from the point of view of an outside observers, over relatively short time scale M log ( M ) .
How do these higher dimensional “microstates” with the smooth and horizonless geometries looks like? We, for the first time, show a specific reduced 4D fuzzball solution has an associated 4D effective fluid near the would-be horizon. The anisotropic pressure of the fluid is crucial to the horizonless geometry.
Applying Kaluza-Klein reduction of non-supersymmetric microstates of the D1-D5-KK system [95]. the metric in 4D is:
d s 4 2 = f 2 A D ( d t + c 1 c 5 ω ) 2 + A D [ d r 2 Δ + d θ 2 + Δ f 2 sin 2 θ d φ 2 ]
Δ = r 2 r 0 2 , f 2 = Δ + r 0 2 n 2 sin 2 θ ,
A = f 2 + 2 p [ ( r r 0 ) + n 2 r 0 ( 1 + cos θ ) ] ,
B = f 2 + 2 r 0 ( r r 0 ) ( n 2 1 ) p r 0 ( 1 + n 2 ) [ ( r r 0 ) + n 2 r 0 ( 1 cos θ ) ] ,
C = 2 r 0 r 0 ( r + r 0 ) n ( n 2 1 ) p r 0 ( 1 + n 2 ) [ ( r r 0 ) + ( p + r 0 ) ( 1 cos θ ) ] ,
G = A f 2 C 2 B 2 , D = B c 1 2 c 5 2 f 2 ( c 1 2 s 5 2 + s 1 2 c 5 2 ) + G f 2 A s 1 2 s 5 2 ,
J 2 = r 0 3 p ( r + r 0 ) n 2 ( n 2 1 ) 2 p r 0 ( 1 + n 2 ) , ω 2 = 2 J sin 2 θ ( r r 0 ) f 2 d φ ,
where parameters c 1 , c 5 , s 1 , s 5 , r 0 , n and p are related to the mass, angular momentum and charges of the solution.

3.5.1. Asymptotic Behavior

Here, we study the asymptotic behavior of metric. As shown in Table 1, it behaves exactly like the Schwarzschild metric when setting K = 1, M = r g . They have different g t φ compared with Kerr metric. As stated, it resembles the Schwarzschild BH far away, but present different geometry close to the horizon.
K = c 1 2 c 5 2 c 1 2 s 5 2 s 1 2 c 5 2
M = p K + c 1 2 c 5 2 J 2 r 0 2 p n 2 ( n 2 1 )
N = 2 r g α r sin 2 θ
L = 2 r g α 3 cos 2 θ sin 3 θ

3.5.2. Matter Field

We can now study the effective 4d matter stress tensor from the Einstein tensor of the 4d fuzzball geometry (33). For a sample choice of parameters, outlined in Table 2, the (diagonalized) Energy-stress tensor T ν μ is:
ρ = 5981 1461 f ( r , θ ) g ( r , θ ) 0 0 0 0 P 1 = f ( r , θ ) 0 0 0 0 P 2 = f ( r , θ ) 0 0 0 0 P 3 = g ( r , θ )
where f ( r , θ ) and g ( r , θ ) are functions of coordinate r and θ . The concrete expressions depend on parameter setting. ρ , P 1 , P 2 and P 3 are the energy density and anisotropic pressure of the matter field. Their behavior near horizon is shown in Figure 2. Pressure P 1 = P 2 is an analytic result and true for any parameter setting, while relationship between the energy and pressure: ρ = 5981 1461 f ( r , θ ) g ( r , θ ) = 5981 1461 P 1 P 3 is a numerical approximation. The approximation is exact far away from the fuzzball with parameters in Table 2 and for any given θ except 0 and π . The relationship changes near the horizon shown as in Figure 3 after averaging over θ . At around r 1000 , we have radial pressure equals tangential pressure P 1 = P 3 . Similar to the metric, matter fields are singular at r = r 0 , θ = 0 and π .
The matter field has an anisotropic pressure. It is not traceless so it cannot be a simple electromagnetic field. We also checked that it cannot be a single scalar field. However, most fuzzball microstates still remain intractable, with no clear dimensional reduction or 4d geometry. To circumvent this obstacle, we proposed a “mock fuzzball” spacetime [96] which captures the horizonless feature of the model with an anisotropic fluid. This conjecture leads to an interesting application to dynamic binary quantum BH merger simulations, as discussed later in Section 6.2.

3.6. Mock Fuzzballs

Here, we introduce a “mock” fuzzball geometry, based on the motivation to build a generic and macroscopic metric which captures important, coarse-grained properties of the fuzzball, i.e., the metric has neither horizon nor singularity and spacetime ends around the stretched horizon. Figure 4 shows how the causal diagram of a Schwarzschild BH is changed to remove the horizon. We study different mock fuzzballs, and check the corresponding matter fields, a swell as potential observable effect. The simplest case of the mock fuzzball 1 is given by:
d s 2 = ( 1 2 M r + b ) d t 2 + 1 1 2 M r d r 2 + r 2 d Ω 2 , r > 2 M
where parameter 0 < b 1 is only important around 2 M . It ensures that g t t does not vanish and thus remove the horizon. The geometry resembles a traditional BH far away. Besides, this metric is only valid where r > 2 M to imitate a fuzzball metric which ends around the stretched horizon. The corresponding matter field is:
P r = T r r = b M 4 π r 2 ( 2 M + r + b r ) ,
P t = T θ θ = T φ φ = b M ( M + r + b r ) 8 π r 2 ( 2 M + r + b r ) 2 ,
other components vanish
The anisotropic behavior of the pressure near the stretched horizon is shown as in Figure 5 with b = 0.01 and M = 1 2 2. The absolute value of two pressures are equal at r = 3 M 1 + b . Tangential pressure is much larger than radial pressure near horizon, and opposite far away. It does not have any singularity at r = 2 M since all fields reach an extremum there.

3.6.1. Other Mock Fuzzballs

Besides the simplest case introduced in the last section, we can modify other terms in Schwarzschild metric to recover energy density which fuzzballs in Section 3.5 actually have. We also study charged and rotating BHs.
  • Schwarzschild metric
The simplest mock fuzzball studied above has no energy density. However, we can recover energy density by changing g θ θ and g φ φ within spherical symmetry:
d s 2 = ( 1 + b 2 M r ) d t 2 + 1 1 2 M r d r 2 + ( 2 d M + ( 1 d ) r ) 2 d Ω 2 , r > 2 M
ρ = T t t = d 4 M 2 + O ( r 2 M )
P r = T r r = 1 4 M 2 + O ( r 2 M )
P t = T θ θ = T φ φ = 1 + 2 b 2 b d 16 b M 2 + O ( r 2 M )
other components vanish
Small b and d ensures that ρ P r P t . In addition, b > 0 , d < 1 ensure finite Ricci scalar without curvature singularity at r > 2 M . Another possible modification to recover energy density is to assume that mass has a time dependence:
d s 2 = ( 1 + b 2 M ( t ) r ) d t 2 + 1 1 2 M ( t ) r d r 2 + r 2 d Ω 2 , r > 2 M
ρ = G t t = b M ( t ) 3 M ( t ) 2 2 b M ( t ) 3 r 2 M ( t ) + 1 8 M ( t ) 2 + ) ( r 2 M ( t ) ) 1 2
P r = G r r = b M ( t ) 3 M ( t ) 2 2 b M ( t ) 3 r 2 M ( t ) 1 8 M ( t ) 2 + O ( r 2 M ( t ) ) 1 2
P t = G θ θ = G φ φ = 3 M ( t ) 2 b ( r 2 M ( t ) ) 2 + M ( t ) 2 b M ( t ) m ( t ) a 2 M ( t ) ( r M ( t ) )
+ b 2 + 2 b 3 4 M ( t ) 2 8 b M ( t ) 2 + 8 b M ( t ) m ( t ) 16 b 3 M ( t ) 2 + O ( r 2 M ( t ) )
other components vanish
  • Extremal BH metric
For an extremal BH, the fuzzball has another interesting property: Proper length from somewhere near stretched horizon to “horizon” is finite, in contrast to the infinite throat in the traditional picture. Parameter c here captures the finite throat.
d s 2 = ( b + ( 1 r q r ) 2 ) d t 2 + 1 c + ( 1 r q r ) 2 d r 2 + r 2 d Ω 2 , r > 2 M
ρ = T t t = 1 c r q 2 + O ( r r q )
P r = T r r = 1 + c r q 2 + O ( r r q )
P t = T θ θ = T φ φ = c b r q 2 + O ( r r q )
other components vanish
  • Non-Extremal BH metric
d s 2 = ( a + 1 r s r + ( r q r ) 2 ) d t 2 + 1 c + 1 r s r + ( r q r ) 2 d r 2 + r 2 d Ω 2 , r > 2 M
ρ = T t t = c r 2 r q 2 r 4
P r = T r r = ( 1 + a ) c r 4 + ( 1 + a ) r 2 r q 2 c r 2 r q 2 r q 4 a r 3 r s + r r q 2 r s r 4 ( r q 2 + r ( r + a r r s ) )
P t = T θ θ = T φ φ dropped for simplicity
other components vanish

3.6.2. What Does an Infalling Observer See?

Assuming Einstein field equations, mock fuzzball geometries can only be sourced by matter fields with exotic (and anisotropic) equations of state. Considering simplest Schwarzschild mock fuzzball:
d s 2 = ( 1 2 M r + b ) d t 2 + 1 1 2 M r d r 2 + r 2 d Ω 2 ,
There is no longer vacuum outside the horizon, which may potentially lead to observable effects, depending on how strongly the “fuzz” matter can interact with the detectors. To visualize the signal, we assume a geodesic observer radially falling towards the stretched horizon with zero velocity at infinity. We calculate (see Appendix A for details) two observable scalars: energy density U and energy flux F , as seen by the observer:
U = T μ ν u μ u ν = 4 b M 2 r 2 ( b r 2 M + r ) 2 ,
F = T μ ν u μ a ν = b 2 M r 3 / 2 ( b r 2 M + r ) 2 ,
where T μ ν is from Einstein field equation of the mock fuzzball, a ν is the unit detector area vector and u μ is the four-velocity of the observer with a μ u μ = 0 . Both energy density and flux are finite and vanish when parameter b vanishes.
To visualize observable energy density and flux, we choose parameters M = 1 and b = 0.1 and compare it with the signal of Hawking radiation as Figure 6 and Figure 7. Both energy density and flux of the mock fuzzball are larger than those of Hawking radiation near the horizon, but drops faster away from the horizon. Energy density of the mock fuzzball is negative while Hawking radiation is positive. Flux of the mock fuzzball flows into BH while Hawking radiation flows into BH near horizon and changes direction away from the horizon.
Specially, energy density and flux of fuzzball are finite near horizon while Hawking radiation diverges. This is because we use Page’s approximation of the Hawking radiation [100], which assumes the Hartle-Hawking state for expectation value of stress-energy tensor:
T t r = 7.44 10 5 4 π r ( r 2 M ) 1 M 2 ,
T r r = 1 1 2 M / r α 1 ( β + 960 M 6 r 6 ) ,
T t t = ( 1 2 M / r ) 3 α 1 ( β 2112 M 6 r 6 ) ,
α = 368640 π 2 M 2 ,
β = 192 M r 5 + 80 M r 4 + 32 M r 3 + 12 M r 2 + 4 M r + 1 ,
The Hartle-Hawking state is a thermal equilibrium states of particles while the Hawking radiation is not in equilibrium. The inconsistency leads to infinite flux and energy density of Hawking radiation. More precise correction can be found in [101].

3.7. Aether Holes and Dark Energy

In 2009, Prescod-Weinstein, Afshordi, and Balogh [9] studied the spherically symmetric solutions of the Gravitational Aether proposal for solving the old cosmological constant problem [102,103]. Surprisingly, they showed that if one sets Planck-scale boundary conditions for aether near the horizons of stellar mass BHs, its pressure will match the observed pressure of dark energy at infinity.
In the Gravitational Aether proposal [102,103], the modified Einstein field equation is given by:
1 8 π G G μ ν = T μ ν 1 4 T α α g μ ν + T μ ν ,
T μ ν = p ( u μ u ν + g μ ν ) ,
where G = 4 3 G N , and then energy-momentum tensor of aether is assumed to be a perfect fluid with stress-energy tensor T μ ν without energy density. Here, quantum vacuum energy decouples from the gravity, as only the traceless part of the matter energy-momentum tensor appears on the right-hand side of the field equations. It can be shown that the Bianchi identity and energy-momentum conservation completely fix the dynamics, and thus the theory has no additional free parameters, or dynamical degrees of freedom, compared to General Relativity.
The modified Schwarzschild metric is the vacuum solution with spherical symmetry in modified equations, and identical to a traditional equations sourced by the aether perfect fluid. Far away from the would-be horizon but close enough to the origin ( 2 M r | p 0 | 1 / 2 ), the solution has the form:
d s 2 = ( 1 + 4 π p 0 r 2 ) d t 2 + d r 2 + r 2 d Ω 2
which can be compared to the de Sitter metric:
d s 2 = ( 1 8 3 π ρ Λ r 2 ) d t 2 + ( 1 8 3 π ρ Λ r 2 ) 1 d r 2 + r 2 d Ω 2
We see that assuming p 0 = 2 3 ρ Λ , the g tt ’s agree with each other. Therefore, the Newtonian observers (for 2 M r | p 0 | 1 / 2 ) will experience the same acceleration as in the de-Sitter metric with the cosmological constant. However, on larger scales, one has to take into account the effects of multiple black holes and other matter in the universe. The Planckian boundary conditions at the (would-be) horizon relates the pressure of the aether to the mass of the astrophysical BHs, p 0 M 3 [9]. In particular, the BH masses within the range 10 M 100 M , which correspond to the most astrophysical BHs in galaxies, yield aether pressures comparable to the pressure of Dark Energy, inferred from cosmic acceleration. Moreover, Ricci scalar is inversely proportional to g t t , so the event horizon where g t t = 0 has a curvature singularity, which is reminiscent of the firewall and fuzzball proposals discussed above.
In particular, the fuzzball paradigm is a good approach to remove the singularity. On the one hand, fuzzball gives an extra anisotropic matter field similar to the aether theory, which stands as a good evidence that quantum effects can modify the Einstein field equation with extra sources of 4d energy-momentum such as aether. Furthermore, fuzzball is a regular and horizonless geometry, which might indicate the singularity is removable in the full quantum picture of BHs.

3.8. 2-2 Holes

In general relativity, gravitational collapse of ordinary matter will always leads to singularities behind trapping horizons [104]. In [13], Holdom and Ren revisited this problem with the asymptotically free quadratic gravity, which could be regarded as a UV completion of general relativity [13]. The quantum quadratic gravity (QQG), whose action is given by:
S QQG = d 4 x g 1 2 M 2 R 1 2 f 2 2 C μ ν α β C μ ν α β + 1 3 f 0 2 R 2 ,
is famously known to be not only asymptotically free, but also perturbatively renormalizable [105,106,107,108]. However, it suffers from a spin-2 ghost due to the higher derivative terms, which is commonly regarded as a pathology of the theory. In [13], it is proposed that the ghost may not be problematic when M is sufficiently small, so that the poles in the perturbative propagators fall into the non-perturbative regime, and the perturbative analysis of ghosts is not reliable. Then it is conjectured [60] that the full graviton propagator in the IR, when M Λ Q Q G , the spin-2 ghost pole is absent in an analogy with the quantum chromodynamics (QCD) where the gluon propagator, describing off-shell gluons, also does not have a pole. Here Λ Q Q G is a certain critical value in QQG, analogous to confinment scale Λ Q C D in QCD. Based on this conjecture, the asymptotically free quadratic action in (81) may involve small quadratic corrections at super-Planckian scale, and so the super-Planckian gravity might be governed by the classical action:
S C Q G = 1 16 π d 4 x g M Pl 2 R α C μ ν α β C μ ν α β + β R 2 .
Since gravitational collapse would involve the super-Planckian energy scale, applying the classical action (82) to such a situation is interesting from a point of view of the quantum gravitational phenomenology. Then the authors in [13] found a solution of horizonless compact object, so-called 2-2 hole, in the classical quadratic gravity. 2-2 holes have an interior with a shrinking volume and a timelike curvature singularity at the origin. It also has a thin-shell configuration, leading to non-zero reflectivity at the would-be horizon, which may cause the emission of GW echoes [21]. Recently, 2-2 holes sourced by thermal gases were also investigated in [109,110].

3.9. Non-Violent Unitarization

A separate class of possible approaches to the BH information paradox involves a violation Postulate 2 in BH complementarity, i.e., non-locality of field equations well outside the stretched horizon, which is dubbed as “nonviolent unitarization” by Steve Giddings [111]. Such a possibility would allow for transfer of information outside horizon around the Page time (e.g., [112,113]), but could also lead to large scale observable deviations from general relativistic predictions in GW and electromagnetic signals [114]. However, it is not clear whether this non-locality is only limited to BH neighborhoods, and if not, how it could affect precision experimental/observational tests in other contexts. Moreover, in contrast to GW echoes that we shall discuss next, it is hard to provide concrete predictions for astrophysical observations in the nonviolent unitarization scenarios.

4. Gravitational Wave Echoes: Predictions

GW echoes may be one of the observable astrophysical signals, a smoking gun, so to speak, for the quantum gravitational processes near BH horizons. Several models of Exotic Compact Objects (ECOs) that we discussed above are expected to emit GW echoes. Some examples are wormholes [10], gravastars [16], and 2-2 holes [13]. Moreover, even Planckian correction in the dispersion relation of gravitational field [14,15,29] and the BH area quantization [115] (but see [116]) may also lead to echo signals. Not only the specific models to reproduce GW echoes but also comprehensive modeling of echo spectra in non-spinning case [117], in spinning case [29,118], and in a semi-analytical way [119,120] were investigated, which enable us to easily obtain echo spectra. In this section, we review the details of GW echoes by starting with the Chandrasekhar-Detweiler (CD) equation [121,122] that is a wave equation with a purely real angular momentum barrier in the Kerr spacetime. We also provide a short review of the GW ringdown signal that is followed by the GW echo, and the superradiance of spinning BHs. The superradiance with a high reflectivity at the would-be horizon may cause the ergoregion instability, which we shall also discuss separately.

4.1. On the Equations Governing the Gravitational Perturbation of Spinning BHs

The GW ringdown is one of the most important signals to probe the structure of BH since it mainly consists of discrete QNMs of BH characterized by mass and spin. In this subsection, we review that the QNMs can be obtained by looking for specific complex frequencies such that the mode functions of GWs satisfy the outgoing boundary condition. Let us start with the CD equation [121,122] that is the wave equation for a spin-s field and has a purely real angular momentum barrier:
d 2 d r 2 V i j s X l m ( r , ω ) = T ,
where T is the source term and the potential V i j with i , j = ± 1 is given by:
V i j = K 2 ( r 2 + a 2 ) 2 + ρ 4 Δ ( r 2 + a 2 ) 2 λ ( λ + 2 ) g + b i Δ b i Δ ρ 8 + ( κ i j ρ 2 Δ h ) ( κ i j ρ 2 g b i h ) ρ 4 ( g + b i Δ ) ( g b i Δ ) 2 + r Δ a m / ω ( r 2 + a 2 ) 2 ρ 2 2 Δ ( r 2 + a 2 ) d d r r Δ a m / ω ( r 2 + a 2 ) 2 ρ 2 .
The functions in (84) are defined by:
b ± 1 ± 3 ( a 2 a m / ω ) ,
κ i j j 36 M 2 2 λ ( a 2 a m / ω ) ( 5 λ + 6 ) 12 a 2 + 2 b i λ ( λ + 2 ) 1 / 2 ,
ρ 2 r 2 + a 2 a m / ω ,
Ξ i Δ 2 ρ 8 ( F + b i ) ,
Θ i j i ω + 1 F b i Δ ρ 2 d F d r κ i j ,
κ ( λ 2 ( λ + 2 ) 2 + 144 a 2 ω 2 ( m a ω ) 2 a 2 ω 2 ( 40 λ 2 48 λ ) + a ω m ( 40 λ 2 + 48 λ ) ) 1 / 2 + 12 i ω M ,
F λ ρ 4 + 3 ρ 2 ( r 2 a 2 ) 3 r 2 Δ Δ ,
g λ ρ 4 + 3 ρ 2 ( r 2 a 2 ) 3 r 2 Δ ,
h g Δ g Δ .
Equation (84) gives four potentials, ( i , j ) = ( 1 , + 1 ) , ( + 1 , 1 ) , ( + 1 , + 1 ) , and ( 1 , 1 ) . One has to use the different potentials in order to cover the whole frequency space with the CD potentials because 1 / ( g + b + 1 Δ ) and 1 / ( g + b 1 Δ ) in (84) are singular in different frequency regions (see Figure 8).
The CD equation is obtained as the generalized Darboux transformation of the Teukolsky equation [123] 3. In the asymptotic regions, r ± , the CD equation reduces to the following wave equation:
d 2 d r 2 + ω ˜ 2 s X l m = T , for r , d 2 d r 2 + ω 2 s X l m = T , for r + ,
where ω ˜ ω m Ω H , in terms horizon angular frequency Ω H a / ( 2 M r + ) , and horizon outre radius r + M + M 2 a 2 of the Kerr BH. In the following, we will omit the subscripts of l, m, and s for brevity. One can read that the homogeneous solutions of (94) are given by the superposition of ingoing and outgoing modes:
X = A e i ω ˜ r + B e i ω ˜ r for r , C e i ω r + D e i ω r for r + ,
where A, B, C, and D are arbitrary constants. The QNMs can be found by looking for the complex frequencies at which the homogeneous solution satisfies the outgoing boundary condition of A = D = 0 . This is equivalent to looking for the zero-points of the Wronskian between the two homogeneous solutions X + and X :
W BH X d X + d r X + d X d r = 2 i ω A in ( ω ) = 2 i ω ˜ B out ( ω ) ,
where the two homogeneous solutions satisfy the following boundary conditions:
X e i ω ˜ r for r , A out e i ω r + A in e i ω r for r + ,
X + B in e i ω ˜ r + B out e i ω ˜ r for r , e i ω r for r + .
Recently, the QNMs of Kerr spacetime were precicely investigated in [124] by using the method developed by Mano, Suzuki, and Takasugi [125,126,127,128] that enables us to obtain the solution of the Teukolsky equation in an analytic way.

4.2. Transmission and Reflection Coefficients of the Angular Momentum Barrier

The GW echoes are results of multiple reflections in the cavity between the would-be horizon (e.g., fuzzball/firewall) and angular momentum barrier (see Figure 9) and so the amplitude of echoes is mainly determined by the reflectivities of the would-be horizon and angular momentum barrier. In this subsection, we review the calculation of the reflectivity of angular momentum barrier.
From the mode functions (97) and (98), one can obtain the energy conservation law for the incident, reflected, and transmitted waves by using another Wronskian relation:
W ˜ BH X d X d r X d X d r ,
which is constant for real frequency due to the reality of the angular momentum barrier in the CD equation. Then we obtain the following relations by using W ˜ BH ( ) = W ˜ BH ( + ) for X and X + :
1 A out A in 2 = ω ˜ ω 1 | A in | 2 ,
1 B in B out 2 = ω ω ˜ 1 | B out | 2 .
From the above relations, one can read that the energy reflectivity and transmissivity for inward incident waves:
I ref | A out / A in | 2 , I trans ( ω ˜ / ω ) | 1 / A in | 2 ,
respectively, and those for outward incident waves are given by:
I ref | B in / B out | 2 , I trans ( ω / ω ˜ ) | 1 / B out | 2 .
From (102) and (103), one can calculate the reflectivity/transmissivity of the angular momentum barrier by numerically solving the homogeneous CD equation (83). In the spinning case a ¯ > 0 , the energy reflectivity is greater than 1 for m Ω H < ω ˜ < 0 , a phenomenon that is often referred to as BH superradiance. The superradiance can be characterized by the amplification factor Z ( I ref / ) 2 1 , and when only interested in the low frequency region, one can use the analytic expression [129]:
Z 4 Q β s l k = 1 l 1 + 4 Q 2 k 2 [ ω ( r + r ) ] 2 l + 1 ,
where r is the radius of inner horizon, β s l ( l s ) ! ( l + s ) ! ( 2 l ) ! ( 2 l + 1 ) ! ! and Q r + 2 + a 2 r + r ω ˜ . To give a few examples of I ref / I ref , we numerically calculate it in the frequency range 0.001 2 M ω 2 , which is shown in Figure 10. As can be seen from Figure 10, the energy flux reflectivity exceeds 1, which means that the energy of a spinning BH is extracted by reflected radiation, within the superradiance regime. We will discuss the ergoregion instability caused by the superradiance and the reflectivity of would-be horizon in Section 4.4.

4.3. Transfer Function of Echo Spectra and Geometric Optics Approximation

When the GW echo is caused by an incident wave packet repeatedly reflected between the cavity, one can use the geometric optics approximation to predict the GW echo signal, which was first pioneered in [117]. Let us start with the calculation of the Green’s function of GW ringdwon signal G BH by using the CD equation. It satisfies:
d 2 d r 2 V G BH ( r , r ) = δ ( r r ) ,
where we omit the subscripts of V i j . Once imposing the outgoing boundary condition, the Green’s function G BH is uniquely determined as:
G BH ( r , r ) = X ( r < ) X + ( r > ) W BH ,
where r < min ( r , r ) and r > max ( r , r ) . Therefore, when there is no reflectivity at the horizon, the Fourier mode of GWs at infinity and at the horizon can be obtained as:
lim r X ( r , ω ) = X + ( r ) + d r X ( r ) T ( r ) W BH X + Z ( ω ) ,
lim r X ( r , ω ) = X ( r ) + d r X + ( r ) T ( r ) W BH X Z BH ( ω ) ,
where T ( r ) is the source for the inhomogeneous CD equation. If there is no reflection near the horizon, the relevant observable spectrum is only Z , and Z BH is irrelevant for observation. On the other hand, if reflection at the would-be horizon is caused by a certain mechanism, Z BH is also observable in addition to Z .
One can obtain echo spectra by using the geometric optics approximation, which should be reliable as long as the would-be horizon and angular momentum barrier are well separated in tortoise coordinates, r . The amplitude of the first echo, Z echo ( 1 ) , can be estimated by:
Z echo ( 1 ) T BH R e 2 i ω ˜ r 0 Z BH ( ω ) ,
and the second echo may have the amplitude of:
Z echo ( 2 ) T BH R 2 R BH e 2 × 2 i ω ˜ r 0 Z BH ( ω ) .
As such, one can obtain the amplitude of n-th echo as:
Z echo ( n ) = T BH R n ( R BH ) n 1 e 2 n i ω ˜ r 0 Z BH ( ω ) ,
where R BH B in / B out and T BH ω / | ω ˜ | B out 1 . Since only the reflectivity and transmissivity of outgoing waves are involved in the echoes, we will not use R BH and T BH in the following, and so omit the symbol →. Summing up all contributions from n = 1 to n = , one obtains:
n = 1 Z BH ( n ) = T BH R e 2 i ω ˜ r 0 1 R R BH e 2 i ω ˜ r 0 Z BH K ( ω ) Z BH .
Please note that here we assume that | R R BH | < 1 , as otherwise the infinite sum of the geometric series does not converge. Finally, the spectrum taking into account the reflection at the would-be horizon is obtained:
X ( r , ω ) = ( Z + K Z BH ) e i ω r = ( 1 + K Z BH / Z ) Z e i ω r .
When R = 0 we have K = 0 and so it reduces to (107). Once we specify a specific form of the source term T, one can obtain Z BH / Z . For example, let us assume the source term located at r = r s :
T ( r ) = S ( ω ) δ ( r r s ) ,
where S ( ω ) is a non-singular function in terms of frequency. Substituting this source term in (107) and (108), one obtains:
Z BH Z = R BH + e 2 i ω ˜ r s T BH .
Please note that this is independent of the function S ( ω ) . Therefore, we finally obtain the following transfer function:
X = e i ω r Z 1 + K Z BH Z = e i ω r Z 1 + K echo + + K echo ,
K echo + R BH R e 2 i ω ˜ r 0 1 R R BH e 2 i ω ˜ r 0 ,
K echo e 2 i ω ˜ r s R e 2 i ω ˜ r 0 1 R R BH e 2 i ω ˜ r 0 .
As discussed in [130], actually K echo + and K echo represent two different trajectories of GWs in the cavity. Here we are interested in outgoing incident waves that is related to K echo + and so in the following we discard K echo from the transfer function, which does not change the qualitative feature of resulting echo signals.
Once we determine the spectrum of injected GWs, Z , one can obtain a template of GW echoes. Using the spectrum of GW ringdown may be a good approximation to obtain a realistic template. In this case, Z is given by [131]:
Z = 2 G M D o A ˜ l m 0 e i ϕ l m 0 S l m 0 ( θ ) α + + e i ϕ l m 0 S l m 0 ( θ ) α ,
α ± Im [ ω l m 0 ] Im [ ω l m 0 ] 2 + ( ω ± ω l m 0 ) 2 ,
where D o is the distance between the GW source and observer, ω l m 0 is the most long-lived QNM, A ˜ l m 0 is the initial ringdown amplitude, ϕ l m 0 is the phase of ringdown GWs, and θ is the observation angle. The amplitude A ˜ l m 0 is proportional to ϵ rd , where ϵ rd E GW / M and E GW is the total energy of GW ringdown [131]. To give a few examples, the ringdown + echo spectra are shown in Figure 11.

4.4. Ergoregion Instability and the QNMs of Quantum BH

As pointed out in the previous subsection, one should check if | R R BH | < 1 is satisfied when calculating the transfer function in the geometric optics picture. This is physically important to understand the ergoregion instability caused by the reflection at the would-be horizon. Since the common ratio of the geometric series is R R BH , the echo amplitude may be amplified and diverges when | R R BH | > 1 . This is nothing but the ergoregion instability that prevents BHs from having high spins, in lieu of dissipation.
Furthermore, the quantum structure of horizon radically modifies the spectrum of QNMs for ECOs, in a similar fashion to trapped modes (or w-modeS) in relativistic stars (e.g., [132,133,134]). Similar to Helio- or Astroseismology, detailed measurement of this spectrum can be used to reconstruct the inner structure of ECO (assuming an effective classical wave equation) [5,6,130].
One can derive the criterion from the QNMs of ECOs. The echo QNMs ω n can be obtain by looking for the poles of the Green’s function of echo GWs. That is, one can look for the poles from the zero points of the denominator of the transfer function:
1 R R BH e 2 i ( ω n m Ω H ) r 0 = 0 ,
and we obtain:
ω n = 2 π n + ( δ + δ ) Δ t echo + m Ω H + i ln | R R BH | Δ t echo ,
where Δ t echo 2 | r 0 | , δ arg [ R ] and δ arg [ R BH ] . Then we obtain the real and imaginary parts of the echo QNMs
Re [ ω n ] 2 π n + ( δ + δ ) Δ t echo + m Ω H ,
Im [ ω n ] ln | R R BH | Δ t echo ω = Re [ ω n ] .
The positivity of the imaginary part of QNMs, which leads to the instability, is equivalent to having | R R BH | > 1 . Furthermore, we can see that the real parts of the QNM frequencies depend on the phases of R and R BH , while their imaginary part depends on their absolute values. We can also rewrite the imaginary part in terms of the amplification factor:
Im [ ω n ] ln | R | Δ t echo + ln ( 1 + Z ) 2 Δ t echo ln | R | Δ t echo + Z 2 Δ t echo .
Then we obtain the analytic form of the imaginary part of QNMs in the low-frequency regime ( M ω 1 ) :
Im [ ω n ] ln | R | Δ t echo + 2 Q Δ t echo β s l k = 1 l 1 + 4 Q 2 k 2 2 π n + ( δ + δ ) Δ t echo + Ω H ( r + r ) 2 l + 1 ,
where we used (104). This is the generalization of the analytic form of QNMs [130]. This analytic form is well consistent with numerically obtained QNMs in the low frequency region as is shown in Figure 12.

4.5. Echoes from Planckian Correction to Dispersion Relation and Boltzmann Reflectivity

The Boltzmann reflection of a BH horizon has been discussed in the context of (stimulated) Hawking radiation from the path integral approach [48], quantum tunneling approach [135,136], and Feynman propagator approach [137]. Recently, two of us studied the reflectivity of a BH for incident GWs from a Lorentz violating dispersion relation and argued that it can be approximated by a Boltzmann-like reflectivity [14]. More recently, three of us used general arguments from thermodynamic detailed balance, fluctuation-dissipation theorem, and CP-symmetry to show that the reflectivity of quantum BH horizons should be universally given by a Boltzmann factor [15,29]:
F out F in horizon = exp ħ ω ˜ k B T H
The reflection of quantum BH might be understood as Hawking radiation stimulated by enormous number of incoming gravitons, and if that is so, having the dependence of the reflectivity on the Hawking temperature T H is natural. Furthermore, one can also avoid the ergoregion instability in this model [15,29]. In this subsection, we briefly review the Boltzmann reflectivity model from both theoretical and phenomenological aspects.

4.5.1. Boltzmann Reflectivity from Dissipation

The dissipative effects at the apparent horizon were discussed from the point of view of the membrane paradigm [51,53], the fluctuating geometry around a BH [138,139], and the minimal length uncertainty principle [140]. Our approach to derive the Boltzmann reflectivity starts with a heuristic assumption to model the dissipative effects, which are expected in any thermodynamic system from fluctuation-dissipation theorem. Let us assume that the wave equation governing the perturbation of BH is given by [15]:
i γ ˜ Ω ( r ) E Pl d 2 d r 2 + d 2 d r 2 + ω ˜ 2 V ( r ) ψ ω ˜ ( r ) = 0 ,
where γ ˜ is a dimensionless dissipation parameter, Ω ( x ) | ω ˜ | / | g 00 ( x ) | is the blueshifted (or proper) frequency, and V is the angular momentum barrier. The form of the dissipation term is expected from the fluctuation-dissipation theorem near the horizon, where the Hawking radiation (quantum fluctuation/dissipation) and the incoming GWs (stimulation) are blue shifted. This dissipative modification to the dispersion relation becomes dominant only when the blueshift effect is so intense that the proper frequency is comparable to the Planck energy, Ω E Pl . Furthermore, from a phenomenological point of view, the dissipative term in (128) is similar to the viscous correction to sound wave propagation in terms of shear viscosity, ν , in Navier-Stokes equation, i ( 4 / 3 ) ν Ω 2 (e.g., [141]).
Let us solve the modified wave equation by imposing a physically reasonable boundary condition (see Figure 13):
ψ ω ˜ constant . for r ,
The constant boundary condition in the limit of r means that the energy flux carried by the ingoing GWs cannot go through the horizon, and is either absorbed or reflected. That is consistent with the BH complementarity [80] or the membrane paradigm [51,142]. Although there is no unique choice of wave equation around a Kerr BH, we here choose the CD equation that has a purely real angular momentum barrier. The modified CD equation is assumed to have the form of:
i γ ˜ | ω ˜ | δ ( r ) E Pl d 2 d r 2 + d 2 d r 2 V ψ ω ˜ = 0 ,
where δ ( r ) 1 r g / r + ( a / r ) 2 is the blue shift factor in terms of the co-rotating frame [30,143]. In the near horizon limit ( r , see below for details), the CD equation reduces to the following form in the limit of r :
i γ ˜ | ω ˜ | Q E Pl e κ + r d 2 d r 2 + d 2 d r 2 ω ˜ 2 ψ ω ˜ = 0 ,
where κ + is the surface acceleration at the outer horizon, Q is defined as:
Q exp 1 2 1 a ¯ 2 r + 2 / r g 2 + a ¯ 2 / 4 r + r g + r 2 / r g 2 + a ¯ 2 / 4 2 1 a ¯ 2 log ( 1 a ¯ 2 ) + 1 2 log 1 a ¯ 2 log r + r g ,
and r ± M ( 1 ± 1 a ¯ 2 ) . The solution of (131) which satisfies the aforementioned boundary condition is:
lim r ψ ω ˜ = 2 F 1 i ω ˜ κ + , i ω ˜ κ + , 1 , i Q E Pl e κ + r γ ˜ | ω ˜ | ,
and one can read that in the intermediate region, κ + 1 log Q E Pl / ( γ ˜ | ω ˜ | ) r ± κ + 1 , ψ ω ˜ can be expressed as the superposition of outgoing and ingoing modes:
ψ ω ˜ = e π ω ˜ / ( 2 κ + ) A + e i ω ˜ r + e π ω ˜ / ( 2 κ + ) A + e i ω ˜ r for ω ˜ > 0 , e π ω ˜ / ( 2 κ + ) A e i ω ˜ r + e π ω ˜ / ( 2 κ + ) A e i ω ˜ r for ω ˜ < 0 ,
where A ± has the form of:
A ± γ ˜ | ω ˜ | Q E Pl i ω ˜ / κ + × Γ ( 2 i ω ˜ / κ + ) Γ ( i ω ˜ / κ + ) Γ ( 1 i ω ˜ / κ + ) = e i ω ˜ r 0 × Γ ( 2 i ω ˜ / κ + ) Γ ( i ω ˜ / κ + ) Γ ( 1 i ω ˜ / κ + ) ,
and:
r 0 1 κ + ln γ ˜ | ω ˜ | / ( Q E Pl ) .
Therefore, the energy reflectivity is given by:
| R | 2 = e 2 π ω ˜ / κ + for ω ˜ > 0 , e 2 π ω ˜ / κ + for ω ˜ < 0 ,
and finally we obtain:
R = exp | ω ˜ | 2 T H + i δ wall ,
where δ wall is the phase shift at the would-be horizon and it is determined by A + or A . Equation (138) then reproduces the Boltzmann energy flux reflectivity in (127). As we noted earlier, the same result can be independently derived using thermodynamic detailed balance or CP symmetry near BH horizons.
When we further modify the dispersion relation by adding a quartic correction term:
Ω ˜ 2 = K ˜ 2 + i γ ˜ Ω ˜ K ˜ 2 C d 2 K ˜ 4 ,
where C d is a constant parameter and Ω ˜ and K ˜ are the proper frequency and proper wavenumber, respectively, the exponent of the Boltzmann factor is modified, and the analytic form can be obtain for C d γ ˜ by using the WKB approximation [14]:
| R | exp 2 + 4 C d 2 / γ ˜ 2 π ( 1 + 4 C d 2 / γ ˜ 2 ) | ω ˜ | 2 T H .
One of the essential differences between the modified dispersion relation model that could give the Boltzmann reflectivity at the would-be horizon, and the Exotic Compact Object (ECO) model is the reflection radius r 0 . In the former case, r 0 depends on the frequency of incoming GWs, ω ˜ , and so the reflection surface is not uniquely determined. This is because the reflection takes place when the frequency | ω ˜ | reaches the Planckian frequency at which the modification in the dispersion relation becomes dominant. Therefore, the reflection radius depends on the initial (asymptotic) frequency of incoming GWs (see Equation (136)). On the other hand, in the ECO scenario, the reflection radius would be fixed and it would stand at ∼ a Planck proper length outside the horizon. In this case, the reflection radius is given by:
r 0 1 κ + ln M E Pl ,
which depends only on the mass of BH. For a detailed discussion of how one can observationally distinguish these scenarios, we refer the reader to [130].

4.5.2. Phenomenology of Boltzmann Reflectivity

Here we summarize some interesting phenomenological aspects of the Boltzmann reflectivity model, and refer the reader to [29] for more details.
Nearly all the previous studies of GW echoes assume a constant reflectivity model, in which the ratio of outgoing to ingoing flux at the horizon is assumed to be independent of frequency. In contrast, the echo spectrum in the Boltzmann reflectivity model can be significantly different. This difference can be seen in the sample echo spectra for both models, shown in Figure 11. As can be seen from the spectra, the echo amplitude is highly excited near m × the horizon frequency, since the Boltzmann reflectivity is sharply peaked around ω m Ω H ± T H and is exponentially suppressed outside this range. In the extremal limit a ¯ 1 , the Hawking temperature becomes zero and so the frequency range in which | R | 1 vanishes. Therefore, the peaks in echo spectrum is highly suppressed for a highly spinning BH (see Figure 14) 4.
This nature of the Boltzmann reflectivity suppresses the ergoregion instability at least up to the Thorne limit a ¯ 0.998 . In [130], a more general case is investigated, where the Hawking temperature in the Boltzmann factor is replaced by the quantum horizon temperature, T H T QH (e.g., as in Equation (140) above)
R = exp | ω ˜ | 2 T QH ,
and the ratio T H / T QH is constrained from the ergoregion instability by using | R R BH | < 1 . The constraint is T H / T QH 0.5 up to the Thorne limit [130] and so the Boltzmann reflectivity ( T H / T QH = 1 ) is safe up to a ¯ 0.998 . As an example, we show a time domain function of ringdown and echo phases with T H / T QH = 0.6 in Figure 15 by implementing the inverse Fourier transform of X = Z ( 1 + K echo + ) , where we choose Z so that it reproduces the ringdown phase [130].
Other notable phenomenological properties of quantum BHs with Boltzmann echoes are [29]:
  • The QNMs of the quantum BH are approximately those of a cavity with a complex length | r 0 | + i ( 4 T QH ) 1 .
  • For γ ˜ 1 (i.e., Planck-scale modifications), the first ∼20 echo amplitudes decay as inverse time 1 / t , and then exponentially.
  • Each QNM of the classical BH can be written as a superposition QNMs of the quantum BH for t < Δ t echo . The superposition can be approximated as a geometric series, leading to a closed-form expression for echo waveforms. In particular, the first 20 echoes have approximate temporal Lorentzian envelopes around their peaks, whose width grows linearly width echo number.
The parameter dependence of echo spectrum in the Boltzmann reflectivity model and its consistency with the tentative detection of echo in GW170817 are also investigated in [130] in more detail.

5. Gravitational Wave Echoes: Observations

For the first time in modern science history we are able to probe the smallest possible theoretical scales or highest possible theoretical energies through GW echoes. The direct observation of GWs [33] was a scientific breakthrough that has opened a vast new frontier in astronomy, providing us with possible tests of General relativity in the extreme physical conditions near the BH horizons. Motivated by the resolutions of BH information paradox that propose alternatives to BH horizons (see Section 3 above), several groups searched the LIGO/Virgo public data for GW echoes [10,16] (see Section 4 and [17] for a review), which has led to claims (and counter-claims) of tentative evidence and/or detection [18,19,20,21,22,23,24,25]. While the origins of these tentative signals remain controversial [23,24,25,26,27,28] they motivate further investigation using improved statistical and theoretical tools, and well as new observations.
In astrophysics, GW echoes from quantum BHs can be seen as a transient signal, coming from the post-coalescence phase of the binary BH merger (Figure 9) or formation of a BH (e.g., via collapse of a hypermassive neutron star; Figure 16). This section will summarize the current status of observational searches for echoes [144].
To properly model echoes, we need a full knowledge of quantum BH nonlinear dynamics, which is so far nonexistent. Therefore, any strategy to search for echoes requires parametrizing one’s ignorance, which has so far taken many shapes and form. Indeed, we need to keep a balance between having a simple tractable model (which may simply miss the real signal), or an exhaustive complex model (which may dilute a weak signal with look-elsewhere effects). Current search methods can be generally split into two: Parametrized template-based methods [18,19,23,24,145,146], and “model-agnostic” coherent methods [20,21,22,25].
Out of these, 8 studies find some observational evidence for echoes [18,19,20,21,22,23,24,25], 3 are comment notes [26,27,28], and 3 more [19,145,147] found no significant echo signals in the binary BH merger events. We can sort them into eight independent groups with 1. positive [18,19,20,21,22], 2. mixed [23,24,25], and 3. negative [19,145,147] results.
For this classification, consistent with common practice (e.g., https://en.wikipedia.org/wiki/P-value#Usage or [148]), our objective criterion for “observational evidence for echoes” is if the study reports p-values 0.05 for noise hypothesis after accounting for the look-elsewhere effects. However, if a study finds low p-values but dismisses it based on theoretical priors, we consider it as “mixed”. An executive summary can be found in Table 3 and Table 4.

5.1. Positive Results

5.1.1. Echoes from the Abyss: Echoes from binary BH mergers O1 by Abedi, Dykaar, and Afshordi (ADA)

The first search for echoes from Planck-scale modifications of general relativity near BH event horizons using the public data release by the Advanced LIGO GW observatory was developed by Abedi, Dykaar, and Afshordi (ADA) [18]. In this search, a naive phenomenological template for echoes was introduced, leading to tentative evidence at false detection probability of 1% (or ≃2.5 σ significance level5 shown in Figure 17, Figure 18 and Figure 19) for the presence of echoes [18]. This work was followed by comments, discussion, and controversy about the origin of this signal [23,24,25,26,27,28]. The ADA model was also later tested for LIGO/Virgo O2 independent events [19], which interestingly, yielded a similar percent-level p-value as O1 (see Section 5.1.3 below). The ADA search was the first phenomenological time-domain echo template search applied to real GW observations [18]. Using a standard GR inspiral-merger-ringdown template M ( t ) , a naive model including five free parameters was proposed:
For a Kerr BH with final mass M BH and dimensionless spin parameter a, the time delay of echoes from Planck-scale modifications of general relativity is: [18,20]:
Δ t echo 4 G M BH c 3 1 + 1 1 a ¯ 2 × ln M BH M planck 0.126 sec M BH 67 M 1 + 1 1 a ¯ 2 .
For the final BH (redshifted) masses and spins reported by the LIGO collaboration for each merger event, echo time delays Δ t echo and their errors constrained inside 1 σ error are as follows [18]:
Δ t echo , I ( sec ) = 0.2925 ± 0.00916 I = GW 150914 0.1778 ± 0.02789 I = GW 151012 0.1013 ± 0.01152 I = GW 151226

5.1.2. Search

In this analysis, ADA devised an echo waveform using theoretical best-fit waveform of Hanford M H , I ( t ) and Livingston M L , I ( t ) detectors (in real time series) for the BBH events, provided by the LIGO and Virgo collaborations. The search used the observed data release for the two detectors, h H , I ( t ) and h L , I ( t ) respectively, at 4096 Hz and for 32 s duration. The devised phenomenological echo waveform which was then constructed using five free parameters:
  • Δ t echo : Time-interval between successive echoes, within their 1 σ range (Equation (144)).
  • t echo : Time of arrival of the first echo, which is related to Δ t echo with corrections ± O ( 1 % ) × Δ t echo due to the non-linear dynamics of the merger.
  • t 0 : Truncation time for GR template with a smooth cut-off function,
    Θ I ( t , t 0 ) 1 2 1 + tanh 1 2 ω I ( t ) ( t t merger t 0 ) ,
    where ω I ( t ) is frequency of GR template as a function of time [150] and t merger is the time of maximum amplitude of the template. It is assumed that t 0 vary within the range t 0 ( 0.1 , 0 ) Δ t ¯ echo . Having this definition, a truncated template is introduced:
    M T , I H / L ( t , t 0 ) Θ I ( t , t 0 ) M I H / L ( t ) .
  • γ : Damping factor of successive echoes, varying between 0.1 and 0.9 .
  • A: Over-all amplitude of the echo waveform with respect to the main event. This free parameter is fitted assuming a flat prior.
The search model of echoes with all the free parameters, assuming a ( 1 ) n + 1 factor due to the phase flip at each reflection, is given by:
M T E , I H / L ( t ) A n = 0 ( 1 ) n + 1 γ n M T , I H / L ( t + t merger t echo n Δ t echo , t 0 ) .
Figure 20 shows this template using the best fit parameters within the range given above along with the main merger event GW150914.
Once this analysis has been completed for GW150914 (loudest event of O1), it has been repeated for rest of the events combined via SNR maximization:
SNR t o t a l 2 I SNR I 2 .
The proposed combination takes same γ and t 0 / Δ t ¯ echo for all events, keeping Δ t echo and A’s as free. The results are shown in Figure 17, Figure 18 and Figure 19 and Table 5 and Table 6.
  • Energy estimation [18]: Given a best-fit template for the echoes, one can provide an estimate for their total GW energy:
    E echoes I / ( M c 2 ) = 0.029 I = GW 150914 , 0.16 I = GW 151012 , 0.047 I = GW 151226 ,

5.1.3. Uchikata et al. Analysis Based on the Template of Abedi, Dykaar, and Afshordi (ADA) for O1 and O2

Uchikata et al. [19] examined GW echo signals for nine binary BH merger events observed by Advanced LIGO and Virgo during the first and second observation runs (O1 and O2 respectively). They used several models for several searches leading to positive and negative results. In this part we bring their positive results and discuss the rest in Section 5.3.3. In this search the critical p-value as 0.05 corresponding to 2 σ significance, (p-value below/above this value) indicates echo signals (are likely/unlikely) to be present in the data.
SNR is evaluated using a matched filter analysis defined using
ρ ( x | h ) = 4 Re f min f max df x ˜ ( f ) h ˜ ( f ) S n ( f )
where x ˜ ( f ) and h ˜ ( f ) are observed data and template in frequency domain respectively, and S n ( f ) is the noise power spectrum of detector. In this analysis they assume frequency band of f max = 2048 Hz, f min = 40 Hz and normalization condition of ( h | h ) = 1 [19]. They then perform an echo search by maximizing SNR, following [18]. Figure 21 presents SNR 2 for the best fit parameters of the event GW150914 (C01)6 with additional consideration of the best fit initial phase7 to the ADA [18] template.
In this part, we show the results using the same template given by ADA [18] except the cut-off parameter t 0 as described in Section 5.1.2 has been fixed to its best fit value, and set the search region of Δ t t e c h o to its 90% (rather than 68%) credible regions in ( a , M ) space. Similar to ADA, the initial phase of the template is also fixed to zero. Since variation of t 0 weakly affects SNR and has an advantage in saving computational costs, they fixed t 0 = 0.1 Δ t e c h o .
Here are the results:
  • O1 events (reanalysis of Westerweck et al. [23]): Since Uchikata et al. [19] followed the same background estimation as Westerweck et al. [23] the results are compared to their O1 results for p-values in Table 7. It is seen that they are almost consistent within the Poisson errors for all events, confirming a marginal p-value of 3–5% for ADA echoes [18].
  • O2 events: Analysis of Uchikata et al. [19] show that the six independent BBH O2 events in Table 8 have similarly small p-values for ADA echoes as O1. As shown in this table, the total p-value for the six O2 events is 0.039. Combining O2 with O1 events shown in Table 7, leads to the total p-value of 0.047.

5.1.4. Echoes from the Abyss: Binary Neutron Star Merger GW170817

A binary neutron star merger event collapsing into a black hole (Figure 16) can also enable us to test general relativity through GW echoes. Although the current LIGO/Virgo/KAGRA detector sensitivity is blind to post-merger ringdown frequency of GWs, they can be sensitive to low frequency echoes harmonics where the ringdown frequency is suppressed by l n ( M BH / M planck ) in Equation (143) [16]. In other words, since the final mass of BNS merger (2–3 M ) is much smaller than that of the binary BH mergers [18], the lowest harmonics n / Δ t echo ( n × 80 Hz) of echo chamber are shifted to the regime of LIGO sensitivity, for small n. Therefore, as first suggested by [21], an optimal model-agnostic search strategy could consist of looking for periodically spaced-harmonics in the frequency space.
Using this model-agnostic search applied to the cross-power spectrum of the two LIGO detectors, Abedi and Afshordi [20] found a tentative detection of echoes around 1.0 sec after the BNS merger, at f e c h o 72 Hz (see Figure 22, Figure 23, Figure 24 and Figure 25) using GW event data GW170817 provided by LIGO/Virgo collaboration [37,152]. As it is shown in Figure 22 and Figure 23, the main signal is also accompanied by secondary lower significance resonances at 73 Hz and t t merger = 32.9 Section It is worth noting that after this detection, Gill et al. [149] used independent Astrophysical considerations, based on GW170817 electromagnetic follow-ups, to determine that the remnant of GW170817 must have collapsed into a BH after t coll = 0.98 0.26 + 0.31 sec, which coincides with the detected GW echo signal at 1.0 second (see Figure 25). This fining of Abedi and Afshordi is consistent with a 2.6–2.7 M “BH” remnant with dimensionless spin 0.84–0.87. For this signal, considering all the “look-elsewher” effects, a significance of 4.2 σ 8 (see Figure 26), or a false alarm probability of 1.6 × 10 5 has been reported, i.e., a similar cross-correlation within the expected frequency/time window after the merger cannot be found more than 4 times in 3 days of GW data. Total energy of detected GW echoes signal using simple assumptions is around ∼ 10 2 M c 2 .

5.1.5. Model-Agnostic Search for Echoes

In this part, we describe the method that leads to the detection of [20], in some detail. The final mass of GW170817 is within ∼2–3 M which could form either a BH or a neutron star (NS). A BNS merger can end up in four possible ways: [152]:
  • A prompt collapse to a BH.
  • A formation of a BH within ≲1 s from a hypermassive NS.
  • Collapse to a BH on timescales of 10– 10 4 s from a supramassive NS.
  • A stable neutron star [37].
Abedi and Afshordi [20] considered the possibility of first and second scenario.
As it was pointed out in Section 5.1.1 the time delay for Planck scale echoes is given by [20],
Δ t echo 4 G M BH c 3 1 + 1 1 a ¯ 2 × ln M BH M planck 4.7 msec M BH 2.7 M 1 + 1 1 a ¯ 2 .
There are two natural frequencies for the waveform of the echoes: The resonance frequencies (natural harmonics) of the echo chamber (formed by the angular momentum barrier and the near-horizon quantum structure), and the BH ringdown (or classical QNM) frequencies. The high frequency harmonics that are initially excited by the merger event decay quickly, while the low frequency harmonics live for longer time [153,154,155]. While the former captures the repeat period of the echoes, the latter describes the echo internal structure. Given that the ringdown frequencies are not resolved by LIGO detector, we can roughly approximate the observable signal as a sum of Dirac delta functions, repeating with the period Δ t echo :
h ( t ) n δ D ( t n Δ t echo t 0 ) h f n δ D ( f n f echo ) .
Therefore, the method searches for coherent periodic peaks of equal amplitude in cross-power spectrum of the two detectors at integer multiples of f echo Δ t echo 1 , in following steps:
  • Fundamental frequency of echoes f echo = Δ t echo 1 within the 90% credible region range for final BH mass and spin is given by:
    63 f echo ( Hz ) 92 .
  • The prior range for echoes search is 0 < t t merger 1 sec .
  • Using amplitude spectral density (ASD) Wiener filter (rather than whiten) the data by dividing by noise variance PSD = ASD2 (rather than ASD):
    H ( t , f ) = Spectrogram IFFT FFT ( h H ( t δ t ) ) P S D H , L ( t , f ) = Spectrogram IFFT FFT ( h L ( t ) ) P S D L .
    where δ t is the time shift between detectors.
  • Cross-Correlating the Obtained Spectrograms and Sum over all the Resonance Frequencies of n × f , Since the polarizations of the LIGO detectors are opposite for GW170817, the real GW signals appears as peaks in X ( t , f ) (see Figure 22 and Figure 24).
    X ( t , f ) = n = 1 10 H ( t , n f ) × L ( t , n f ) .
The simplicity of the method (not having any arbitrary or ad-hoc cuts or parameters) along with its high significance are reasons for making this finding more interesting and reliable.

5.1.6. “GW Echoes through New Windows” by Conklin et al.

In this search, three methods (named as I, II, III) that are based on general properties of echoes has been suggested. Conklin et al. [21] mostly focused on Method II, which is based on frequency windows, while the other two methods use time windows. Window functions in these methods allow us to find quasiperiodic structures in time and/or frequency domains. Method II turns out to be the most successful one. Accordingly, in this paper we just review this method. The search methods become more optimal using correlations of data in multiple detectors. Using the suggested method and search [21] find significant evidence for GW echoes, which is shown in Table 9 (and see Figure 27 for GW170104) for both O1 and O2 LIGO/Virgo observations.
The time delay of echoes, Δ t echo which is similar to the Planckian echoes in Equation (143), is paramaterized in [21] using:
Δ t echo / M C H R = η 1 + 1 1 a ¯ 2 ln r 0 r + M
where r 0 is the location of the quantum structure outside r + . Here η = 2 corresponds to proper Planck length (Equation (143)).
Taking into account of errors in final mass, spin, and redshift it is realised that (see Figure 28) the echoes found are consistent with η = 1.7 [21,22]. Best fit properties of the peaks for O1 and O2 events also shown in Table 10 [22]. Alternatively, this result can be interpreted as the energy scale for reflection from quantum horizons to be 6 ± 2 orders of magnitude below Planck energy [90].

5.1.7. Comment on: “Gravitational Wave Echoes through New Windows by Conklin et al.”

Conklin et al. [21] remark that “We have not found signals for the two earlier events, GW150914 and GW151012, which play a significant role in ADA results in [18]”. While in their updated search [22] they indicate existence of signals for both GW150914 and GW151012, although with no p-value estimation which is crucial in this search. Therefore, it is hard to truly evaluate the significance of echoes reported in [22]. Moreover, it is not clear how much the choices made in Method II in [21] might have been affected by a posteriori statistics.

5.2. Mixed Results

5.2.1. Results of Westerweck et al. and Nielsen et al.

Westerweck et al. [23] re-analysed the same model proposed by ADA [18], using more background data and a modified procedure. They focused on the data analysis methods of ADA [18] and their significance estimation, namely the concerns presented in [26] suggesting a different significance estimate using 4096 s of LOSC data.
The results of p-value estimation and comparison with ADA results are shown in Table 11. In addition, Nielsen et al. [24] searched for echoes signals in GW data via Bayesian model selection probabilities, comparing signal and no-signal hypotheses using ADA model [18]. Accordingly, calculation of Bayes factors for the ADA model in O1 events presented in Table 12.
In the following we explain the results given in different plots:
  • In Figure 29 it is shown that depending on the overall amplitude of the injection, the signal either can be recovered or it would be difficult to recover.
  • Figure 30 shows injected and recovered values for γ having different overall amplitude A used in [18]. Although this plot shows a preference for γ = 1 , having high value of γ can be recovered easily.
  • Figure 31 shows injection of echoes signals into the Gaussian noise having different overall amplitudes A.
  • Curiously, we note that the Maximum SNR = 2.86 reported in [24] for GW150914 (Table 12), the only event with negative evidence for ADA echoes, is significantly lower than the one reported by the same group [23] and ADA [18] ( = 4.21 , see Table 6), while the numbers are consistent for other O1 events.

5.2.2. Comment on: “Low Significance of Evidence for BH Echoes in Gravitational Wave Data”

  • Comment on Search Strategy:
    As described in former part in Figure 30, Westerweck et al. [23] demonstrated the preference of γ = 1 at lower injection amplitudes for ADA model [18]. This is expected result as γ 1 extends the template to infinity. In other words, it extends the template range to infinite time, which is clearly dominated by noise. However, γ = 1 is still far from the best-fit γ = 0.9 on edge of the prior where at least 90% of energy goes to first 11 echoes. Besides, initial waveform must change significantly in subsequent echoes. Indeed, repeated template that does not damp ( γ = 1 ) is not physical. One solution to this problem might be to use a finite range of data (which is used by Westerweck et al. [23]) making the result strongly dependent on what portion of data has been taken.
    It should be also noted that using 1-sigma range for errors in Δ t echo , implies a 32% chance for the signal to be missed that causes reduced significance by diluting SNR 2 from some events.
    Westerweck et al. [23] found that the least significant event LVT151012 which is now called GW151012 has the most contribution to tentative evidence for echoes. This peculiar finding does not disfavor echoes as there is no simply reasonable justification that significance of echoes should be directly related to the significance of main event. Additionally, Wang et al. [155] showed that by changing only ± 20 % of frequency of initial condition of echoes the SNR 2 for echoes can change by 3 orders of magnitude. As BBH events have different component spins and mass ratios we might expect a significant diversity in relative echo signal amplitude for each of them. Interestingly, as will be discussed in next Section 5.2.3 mass ratio of BBH events appears to show correlation with the echo amplitude [144].
  • Comment on Abstract and Conclusion:
    The most crucial comments for Westerweck et al. [23] goes to their abstract and conclusion (also provided in [28]). Although ADA [28] strongly acknowledge the analysis by Westerweck et al. [23], which is a careful re-evaluation of ADA analysis, the Abstract/Conclusion of Westerweck et al. [23] misrepresents their finding. The most critical point of this misrepresentation is in Abstract claiming “a reduced statistical significance… entirely consistent with noise”. Contrasted to this claim in their Table I (Table 11 in this paper) they found p-value = 0.020 for the noise hypothesis, with the same model and data as in ADA (as opposed to 0.011 in ADA [18]). However, if one follows standard nomenclature (e.g., https://en.wikipedia.org/wiki/P-value#Usage or [148]), p-values < 0.05 disfavour noise hypothesis, providing “moderate to strong” Bayesian evidence for echoes, which is contrary to what they state in theirAbstract.
    To conclude, considering all the critiques of Westerweck et al. [23], we see NO evidence that their improved analysis with p-value = 0.020 ± 0.009 has reduced the significance of echoes, entirely consistent with p-value = 0.011 of ADA [18]. The fact that completely independent events of O2 also show a low p-value = 0.039 for ADA echoes (see Table 8) further boosts the statistical evidence for ADA model in LIGO/Virgo data.

5.2.3. Results of Salemi et al.

Another independent group [25] has found similar post-merger GW signals that can be attributed to GW echoes. However, the setup of this methodology, which is based on coherent WaveBurst (cWB) [156] method, was not originally developed to search for echoes. This search, which is independent of the waveform models, has been developed based on coherent excess power in events from the GWTC-1 (catalog of compact binary coalescence). Here, loose bounds on the duration and bandwidth of the signal leads to evaluation of coherent response of independent detectors.
This search has focused on detected features as deviations from GR and has presented the method to obtain their significance. It appears that from eleven events reported in the GWTC-1, two of them (GW151012 and GW151226) in Figure 32 and Figure 33 respectively, show an excess of coherent energy after the merger ( Δ t 0.2 s and 0.1 s, respectively) with p-values (0.004 and 0.03, respectively). However, [25] showed that (Figure 34) the post-merger signal from GW151012 favours different sky location than that of the main event.
In Table 13, SNR pc minR , and its upper and lower bounds are presented. In the last column of this table, estimated p-values for postmerger features are reported. Along with, 2 σ upper and lower bounds for SNR, upper and lower bounds for the p-values are reported.

5.2.4. Hint of Dependence of Significance of Echoes on Binary BH Mass Ratio, Comment on: Results of Salemi et al.

In this part, we first re-examine the interpretation of Salemi et al. [25] about the signals they found and then give a more conclusive support for echoes hypothesis.
  • Comment on: Results of Salemi et al. [25]
    Salemi et al. [25] disfavoured echoes hypothesis pointing that post-merger signal of GW151012 has arrived from a different sky location than that of the main event. However, the p-value 0.004 of this secondary signal, disfavours two signals being unrelated. We see that all the secondary (post-merger) clusters they claim as signals in Figure 32 and Figure 33 are nearly monochromatic. That means the waveform of these signals are quasi-periodic, leading to degeneracies in inferred time-delays. So looking again to the secondary signal of GW151012, the null (residual) plot in Figure 32 shows the peak of the cluster (which is mostly responsible to make a different sky localization) is at ∼130 Hz which corresponds to 7.7 m s time delay. This 7.7 m s is the same time delay of first peak and second peak in Figure 34 for post-coalescence signal (green), confirming that the monochromatic degeneracy in Figure 34 might have caused the different sky localization. Figure 35 shows better interpretation of cause of this error.
  • Dependence of Significance of Echoes on Binary BH Mass Ratio
    Consider a system of binary BHs (BBH), with progenitor masses m 1 and m 2 . It is known that these systems consist of two almost equal mass BHs m 1 m 2 . However, the inevitable diversity in the initial conditions, specifically binary mass ratio, can lead to different echo properties. Here, we review the evidence for correlation between the significance of echoes and the progenitor BBH mass ratio (first presented in [144]):
1.
We used LIGO parameter estimation samples for BBH events provided in [157]. Then obtained mass ratios by weighting all events as equal. We used full m1 vs. m2 distribution samples for “Overall_posterior”. We fit these posterior points to a straight line shown in Figure 36.
2.
We take p-values reported in [25] for each event post-coalescence signal and plot the best fit line of mass ratio vs. ( log ( p -value)) 0.5 . We used least square method [158] in fitting a straight line consisting all the posterior points provided by LIGO. Finally, we use the slope of best fit line as our primary parameter in determining significance.
3.
Finally, we obtain the significance (Figure 37) assuming that there is no relation between p-value and mass ratio by taking random p-values for BBH Catalog events within the uniform range 0 < ( log ( p value ) ) 0.5 < 2.5 ( 0.0019 < p-values <   1 ). Indeed, having no relation between p-value and mass ratio of events shall end up with zero slope in large number of random selections. Therefore, in order to find a false detection rate for correlation, number of slopes higher than the actual measured slope is calculated (see Figure 37). Accounting for the “look elsewhere” effect, we find tentative hint of mass-ratio dependence of echo significance reported in [25] at false detection probability of 1%.
Table 14 shows events and p-value of their post-coalescence signal reported in [25] versus expected Planckian echo time delays and average mass ratios. At a glance it is seen that smallest mass ratios go to smallest p-values. This can happen with 1/10 chance out of 10 BBH events. For two most significant events, being also the most extreme BBHs, the random chance becomes p-value = 1 9 × 1 10 = 0.011 which is consistent with the statistics of regression analysis in Figure 37.

5.3. Negative Results

5.3.1. Template-Based Gravitational-Wave Echoes Search Using Bayesian Model Selection by Lo et al.

Lo et al. [145] found that using a wider range of priors (listed in Table 15) compared to that of ADA model [18], including the main event template in the echo template, and keeping just three echoes, leads to lower significance on echo signals evaluated by the Bayes factor using Bayesian analysis.
They considered two hypotheses as the null hypothesis H 0 and alternative hypothesis H 1 ,
H 0 No echoes in the data d = n + h IMR H 1 There are echoes in the data d = n + h IMRE
where d and n indicate the GW data, and the instrumental noise respectively and h IMR , and h IMRE are the inspiral-merger-ringdown (IMR) gravitational-wave signal and inspiral-merger-ringdown-echo (IMRE) gravitational-wave signal respectively.
The log Bayes factor ln B in their search method gives the detection statistics which determines whether there is an IMRE signal or an IMR signal in data. From the Bayesian perspective if the log Bayes factor, is greater than 0, we can conclude that the data favor the alternative hypothesis.
The relation between p-value and null distribution of detection statistic ln B is given by
p value = Pr ( ln B ln B detected | H 0 ) = 1 ln B detected p ( ln B | H 0 ) d ln B ,
where ln B detected is the detection statistic which is obtained using a segment of data in analysis, and p ( ln B | H 0 ) is called the null distribution of ln B , i.e., the distribution of ln B assuming that H 0 is true.
Lo et al. injected an IMRE injection of template with echo parameters discussed earlier for the event GW150914 into simulated Gaussian noise. The detection statistic for this injection is given as follows,
ln B detected , Gaussian = 0.2576 < 0 .
Therefore, this finding indicates that the data slightly favor the null hypothesis from Bayesian analysis point of view. While the p-value and the corresponding statistical significance, for Gaussian noise, show that the data favors echo hypothesis,
p value = 0.01275 , statistical significance = 2.234 σ .
Table 16 shows the values of the detection statistic ln B vs. corresponding statistical significance in Gaussian and O1 backgrounds. Therefore, for a detection of gravitational-wave echoes having statistical significance 5 σ , the detection threshold would be,
ln B threshold , Gaussian = 1.9 , ln B threshold , O 1 = 5.7 ,
for Gaussian noise and O1 noise respectively.
Table 17 shows the detection statistics and the corresponding statistical significance and p-value for the O1 events. This table also shows that the ordering of the events by their statistical significance is consistent with what has been reported by Nielsen et al. [24].
ln B O 1 ( cat ) = 1.1 ,

5.3.2. Comments on Lo et al.

In their analysis, Lo et al. [145] included both the main event, as well as the ADA echo waveform in their template, but they used expanded priors in Table 15. Although expanding the ADA priors covers a larger space of possibilities, it tends to dilute marginal signals and bury them in the noise. For example, there is no good physical interpretation for repeating echoes that do not damp (with γ = 1 ), as they violate energy conservation.
For these reasons, it is not surprising that Lo et al. [145] find a smaller Bayes factors, due to their expanded priors. However, it is well-known that expanding prior into nonphysical regimes will artificially lower Bayesian evidence for any model, especially since γ = 1 is a (formal) singularity of the likelihood function.

5.3.3. Results of Uchikata et al.

As discussed in Section 5.1.3 above, Uchikata et al. [19] can approximately reproduce the evidence for ADA echoes in both O1 and O2 events. Here, we present their search using an alternative template that failed to find any evidence for echoes. To build the latter, they considered Kerr spacetime, replacing the event horizon with a reflective membrane. They then used the transmissivity of the Kerr angular momentum barrier T BH ( ω ) to filter the ADA template, which acts as a hi-pass filter (truncating the low frequency part of the ADA phenomenological waveform 147). Moreover, the overall phase shift of the waveform as a free parameter is taken into account contrary to ADA search. Using this template, they have found no significant echo signals in the binary BH merger events. The background estimation, has used the same method provided by Westerweck et al. [23].
  • The results in Table 18 gives p-values for all events. The combined p-value is well above the critical p-value 0.05. In other words, echo signals using this model do not exist in the data, or their amplitudes are too small to be detected within the current detector sensitivity.
  • Since the phase shift at the membrane (due to the reflection and boundary condition) is model dependent, it is physically reasonable to assume a total phase shift as a parameter (see Figure 21 that has used extra phase parameter for GW150914). In contrast, former studies [18,23] only considered phase inversion (or Dirichlet boundary conditions) at the reflective membrane. Therefore, the results of two cases, when the phase shift is fixed to π (result 1) and when it is a free parameter (result 2), respectively, in Table 19 are compared. In this table we see that, p-values become slightly larger (taking GW151226 as an exception) when the phase shift due to the reflection has taken as a free parameter.

5.3.4. Comment on Negative Results of Uchikata et al.

Model provided by Uchikata et al. [19] substantially truncates the low frequency part of the GR waveform (which is the basis of ADA template 147). However, one may argue that the GR waveform has already been filtered once by the transmissivity of the angular momentum barrier, T BH ( ω ) as it is what is seen by observers at infinity. In fact, [29] have shown that both GR signal (main event) and echoes can be constructed from a superposition of the QNMs of the quantum BH, which are essentially the modes trapped between the angular momentum barrier and the quantum membrane (see Section 4.4 above). Therefore, both the GR signal and the echoes pass through same barrier and are thus truncated by the same T BH ( ω ) . This implies that T BH ( ω ) cancels in the ratio of echo to main signal waveform, and in contrast to Uchikata et al. [19], no truncation is needed.
Indeed, as evidenced by their own analysis (Section 5.1.3 above), the low-frequency part of the ADA template is necessary to obtain a significant signal. This is physically justified since the high frequencies leak out of the angular momentum barrier quickly, leading to a rapid decay, while echoes can last much longer at lower frequencies [155].

5.3.5. Results of Tsang et al.

Tsang et al. [146,147] proposed a morphology-independent search method which consists of a large number of free parameters for echoes compared to ADA model [18] (9 versus 5). They search for echoes in all the significant events in (GWTC-1), and found that for all the events, the ratios of evidences for signal versus noise and signal versus glitch do not rise above their respective background. Only the smallest p-value = 3% goes to the event GW170823. Hence they found no significant evidence for echoes in GWTC-1. The results of search are given in Table 20 and Table 21.

5.3.6. Comment on: Results of Tsang et al.

Tsang et al. [147] developed a model which consists of a large number of free parameters (9, based on our correspondence with the authors). Indeed, a larger space of possibilities leads to lower significance. So it is not surprising to get large p-values out of large free-parameter space. Indeed, all the SNRs for echoes reported in ADA would be below the detection threshold reported by [146], given their large number of parameters (see [144] for more detail).

5.4. A Concordant Picture of Echoes

To have an optimal search for echoes, one may want to take the following guidelines into consideration:
  • Have a good physical model (or you will not find them!)
  • Use a simple template (avoid too many arbitrary choices)
  • Avoid a posteriori statistics (don not look at data to make your model)
Based on the positive (Section 5.1), mixed (Section 5.2), and negative (Section 5.3) results, one may offer the following general observations:
  • Coherent searches appear to give more significant evidence for echoes
  • Template searches can find evidence for echoes, if they include lower frequencies
  • Models with a large number of free parameters and/or wider priors can weaken the echoes below the detection threshold
An executive summary of these observations is shown in Table 3 and Table 4 as positive evidence (p-value 0.05 ) and failed results, respectively.
The study of post-merger GW observations with the above-mentioned motivations led to tentative signals, at varying levels of significance, by different groups. We shall outline several similarities among these findings below. However, it is also important to note that these similarities do not mean that the signals found are the same, but it does provide a preponderance of corroborating evidence for GW echoes in current observations.

5.4.1. Five Independent Groups, Five Independent Methods, Identical Results!

  • In [18] (Table 2), using the reported masses and spins of LIGO O1 events, the time delays of 0.1 s and 0.2 s for GW151226 and GW151012 were predicted respectively for Planckian echoes. These happen to be exactly the same as the times for post-merger signals found in [25].
  • Results of [18,19,20,23,24,25] all are consistent with Planckian echoes at p-values of O ( % ) .
  • Furthermore, the reconstructed detector responses for GW151226 and GW151012 [159,160] in [25] give consistent amplitudes (0.33, 0.34) × (maximum amplitude of main event) compared with [18] (Table 2). Energy reported in [18] (Appendix A) is also consistent with strength of signals found in [25]. Finally, SNR reported for GW151012 in [18] (Table 2 and Figure 6) has highest value which is also consistent with highest significance event in [25].
  • Log Base factor values in Table 2 of [24], where they found positive evidence for ADA echoes in GW151012 and GW151226 (where GW151012 is more significant) is consistent with significance of signals found in [25].
  • Also note that the echo signal of GW150914 [18] at time delay 0.3 s had narrowest time window ( ± 3 % in Table 2) and smallest energy (Table 2) compared to GW151226 and GW151012, which could explain its absence in [25], and no evidence (negative Log Base factor value) in [24] Table 2.
  • Nevertheless, the residual signal in [161] which is a supporting results for [25] is consistent with 300 m sec echo signal time delay in [18] (Table 2).
  • The percent-level evidence of Uchikata et al. [19] for ADA echoes in O1 events (shown in Table 7) are consistent with the results of other groups [18,23,24,25].
  • The results of Uchikata et al. [19] for O2 events shown in Table 8 are given ∼4% overall p-value which are small as O1 events.
  • Lo et al. [145] Table 17 by adding main event to the ADA waveform and keeping only three echoes with larger prior ranges also found similar ordering of events by their statistical significance with what reported by Nielsen et al. [24].

5.4.2. Other Findings

Here we present similarities found for echoes in binary neutron star merger GW170817:
  • We note that the time-scale of 1.0 sec after merger for collapse into BH (first reported by Abedi and Afshordi [20]) is now also independently found from purely Astrophysical considerations by [149], who found t coll = 0.98 0.26 + 0.31 s.
  • Along with echo signal found by Abedi and Afshordi [20] another group [21] claimed evidence for an echo frequency of f echo ( 0.00719 sec ) 1 = 139 Hz for GW170817, with a p-value of 1 / 300 . Noting the proximity of this value to the second harmonic of Abedi and Afshordi finding with echo frequency 2 × f echo = 144 Hz, it is feasible that the two different methods are seeing (different harmonics of) the same echo signal. However, the method applied in [21] is sub-optimal, as they whiten (rather than Wiener filter in [20]) the data, and thus could underestimate the significance of the correlation peak they found (see [28] for further discussion).
To contrast, let us point out two apparent inconsistencies:
  • The results of Conklin et al. [21,22] are finding non-Planckian echoes signal which might not be consistent with the results of other groups [18,20,23,24,25]. However, it may not be appropriate to do a one-to-one comparison of [21,22] to other studies, as the employed method is significantly different.
  • Results of Salemi et al. [25] show that using their reconstruction of cWB, the post merger signal of GW151012 seen in Figure 32 appear to come from a different sky location, compared to the main event signal. However, as we discussed above (Section 5.2.3), this might be due to time-delay degeneracy in the cWB monochromatic signals.

5.4.3. Independent Confirmations of Model Predictions

Here, we outline model predictions that have been confirmed using independent data by independent groups
  • Binary BH mergers:
    Uchikata et al. [19] used the ADA model [18] to search for echoes for both O1 (first observing run) and O2 (second observing run) while the original search of ADA only covers O1. The results for O2 (with p-value = 0.039) in Table 7 and Table 8 show similar evidence as O1 (with p-value = 0.055).
  • Binary neutron star merger:
    After detection of echoes signal with 4.2 σ significance around 1 s after BNS merger GW170817 [20] where Abedi and Afshordi claim that it has collapsed to BH at this time, Gill et al. [149] with independent Astrophysical consideration also determined that the remnant of GW170817 must have collapsed to a BH after t coll = 0.98 0.26 + 0.31 Section Error-bar for this observation compared to the detected signal of echoes by Abedi and Afshordi as a consequence of BH collapse is shown in Figure 25.

5.4.4. Concerns about ADA Searches

  • Concerns about errors in Δ t e c h o : The original ADA search [18] had used an ad-hoc method for finding symmetric 1-sigma errors for Δ t e c h o , which would miss 1 / 3 of Planckian echoes. Furthermore, the actual LIGO posteriors for these parameters (which are now publicly available, even though they were not at the time), are not Gaussian. A fully Bayesian Bayesian search, using actual mass and spin posteriors would avoid these short-comings.
  • Concerns about keeping t 0 = 0.1 and γ = 0.9 fixed: Uchikata et al. [19] fixed t 0 at its best fit value of O1 in order to search in O2. It might be a good idea to keep these parameters at their best-fit values (from O1 or O2), to make perform more efficient searches in e.g., O3 events.
  • Concerns about mass ratio dependence of echoes overall amplitude: As pointed out in Section 5.2.3 (Figure 36), it appears that the amplitudes of echo signals may depend on the BBH mass ratio, which should be taken into account if one wants to optimally combine echo signals in different events.

5.4.5. On Negative GW Echo Searches

For attempts that fail to yield any evidence for echoes [19,146], we again point out that an optimal echo search should us a simple model with minimum number of free parameters. The current positive results turn out to be weak signals with SNR 4 , which is below the threshold for those searches that consist of many free parameters ( S N R > 8 for 9 free parameters). In addition, negative results of Uchikata et al. [19] indicate that the echo signals found in [18] mostly consist of low frequency modes, which is independently confirmed in [20], and Uchikata et al’s own search in both O1 and O2, using the original ADA template.

5.4.6. Non-Gaussianity of Backgrounds

To search for signals of a given echo model, we need a proper understanding of background behaviour. Only then we might be able to determine the best statistical methodology. Since ADA [18] used different parts of LIGO data to get a combined significance, one may think about what would be the best method of combination of separate sets of data with different background behaviour. Then, we can also ask how the search changes by including three or higher number of detectors.
It was already observed that LIGO noise vary significantly and is very non-gaussian over long time-scales (see Figures 14 and 15 in [150]). This either non-stationary or non-Gaussian background makes the interpretation of p-value ambiguous, particularly in finding marginal echo signals which are often near the detection threshold. Therefore, one must examine how much this varying background affects the inferred significance of a detection. This is studied by looking at other different minute-long stretches of data within a minute of the main events [18]. As can be seen in Figure 38, the variation of p-values at the tail of the distribution is much higher than what is expected from Poisson statistics of the SNR peaks. This becomes more interesting when we see that the smallest p-value is coming from the range which is closest to the main event. This might be because the marginal LVT151012 (now called GW151012) detection is over a minimum of the LIGO (combined) detector noise.

6. Future Prospects

6.1. Towards Synergistic Statistical Methodologies

As we summarized in the previous sections, the past fours years witnessed hundreds of theoretical studies focusing on model-building for echoes, as well as dozens of observational searches and statistical methodologies. However, despite remarkable progress on both fronts, the theoretical and observational tracks largely developed independently. However, it appears that both tracks have become mature enough, so much so that the time is ripe for a synergistic convergence. For example, Bayesian methods developed in [24,145] applied to a superposition of QNMs of quantum BHs (as outlined in [29]) would put coherent methods developed by [21,28] on more sound statistical and physical footings. The analogy will be with helio- or astro-seismology, where modeling a dense spectrum of QNM frequencies can be used to infer the intrenal structure of the compact objects [130].
The real challenge will be in allowing enough freedom in our best physical models, in order to capture all the remaining theoretical uncertainties, but not any more!

6.2. Echoes in Numerical Relativity

Most studies of echoes have so far focused on the linear perturbation theory around the final BH for simplicity, but in reality the mergers start with the highly nonlinear binary BH inspiral. Hence, we need a covariant numerical implementation of binary quantum BHs within a highly nonlinear dynamical spacetime to fully address the entire dynamics, especially the initial conditions. There are several possible approaches borrowed from numerical relativity which can be modified to either include the quantum boundary condition or the full dynamics of binary quantum BHs.
For instance, the effective one body (EOB) formalism [162,163] is a concrete strategy which only needs to solve ordinary differential equations rather than to perform the costly 3d numerical relativity simulations. It uses higher-order post-Newtonian expansion in a resummed form (different from the usual the Taylor-expansion), to include the non-perturbative result using a conservative description of binary BHs dynamics, radiation-reaction and emitted GW waveform. One possible approach, that is currently underway, is to capture the nonlinear effects in echoes by modifying the boundary condition in the EOB codes to implement the quantum BH dynamics.
Another route is to directly modify numerical relativity codes that successfully produced waveforms for BBH merger events. A concrete strategy could be incorporating the mock fuzzball energy-momentum tensor (Section 3.6) as a source for Einstein equations, directly into the numerical relativity codes. If the fuzzball “fluid” manages to stay just outside the apparent BH horizons in a dynamical setting, then it can potentially generate echoes in a fully nonlinear numerical simulation of quantum BBH merger.
Recently, [164,165] presented the first numerical simulation of BBH mergers in Chern-Simon gravity. They start with the modified action and predict the dynamics order by order. It is possible that a similar iterative approach can be applied to model boundary conditions at apparent horizons, or evolution of mock fuzzballs.

6.3. Quantum Gravity, Holography, and Echoes

As we discussed in Section 4 above, any modification of event horizons that could lead to echoes should be a non-perturbative modification of general relativity, and can only be fully captured by a non-perturbative description of quantum gravity. A possible example of this is the fuzzball program in string theory (Section 3 above). However, more generally, what can non-perturbative approaches to quantum gravity tell us about BH echoes?
One of our greatest insights into the dynamics of quantum gravity has come from the Holographic Principle, that extending Bekenstein-Hawking area law for entropy of BHs [42], suggests the entire dynamics of a quantum gravitational system should be captured on its boundary. The most concrete realization of this principle was proposed by Juan Maldacena [94], in the form a conjectured duality between quantum gravity in Anti-de Sitter (AdS) spacetime and a Conformal field theory (CFT), commonly known as AdS/CFT correspondence or conjecture. It proposes that CFT in spacetime of d-1 dimension, at the asymptotic boundary of an AdS spacetime is mathematically equivalent to string theory (or quantum gravity) within the bulk AdS in d dimension. This topic has been extremely fruitful over the past two decades, offering many synergies between seemingly disparate notions in geometry and quantum information. For example, the Ryu-Takayanagi conjecture [166] relates the entanglement entropy of boundary CFTs with the areas of extremal surfaces in the bulk AdS, generalizing the notion of Bekenstein-Hawking BH enetropy to arbitrary geometries.
An intriguing connection between AdS/CFT and echoes is the appearance of echo times:
Δ t echo = t scrambling = ln ( S BH ) 2 π T H ,
as “scrambling time”, in the AdS/CFT literature [167]. Here, S BH and T H are the entropy and temperature of the BH respectively. The scrambling time refers to the time it takes to destroy quantum entanglements in a chaotic system, while BHs (and their CFT duals) are conjectured to be fast scramblers, i.e., the most efficient in destroying entanglement (e.g., [168]). Interestingly, Saraswat and Afshordi [169] recently showed that the scrambling time (computed using Ryu-Takayanagi conjecture in a dynamical setting) is identical to the Planckian echo times, for generic charged AdS BHs. Could this imply that echoes could be a generic property of (possibly a certain class of) quantum chaotic systems?
Another possible connection could come in the form of the fluid-gravity correspondence, e.g., in the context of membrane paradigm discussed in Section 2.2. For example, in [15], we argued that Boltzmann reflectivity of GW echoes, implies that viscosity of the boundary fluid should vanish at small frequencies ħ ω k T . One may also speculate that other holographic manifestations of BH echoes may appear in the Kerr/CFT conjecture [170], Braneworld BHs [171], or as Regge poles of the boundary plasma in AdS/CFT.

6.4. Einstein Telescope, Cosmic Explorer

The Einstein Telescope (ET) [172] and Cosmic Explore (CE) [173] are the third-generation ground-based GW detectors. The ET consists of three underground detectors with three arms 10 km long and CE will be realized with two arms 40 km long, which are 10 times longer than Advanced LIGO’s. These next-generation GW detectors might allow us to observe some Planckian signatures from quantum BHs such as GW echoes from merger events leading to a remnant BH. We plot the spectra of GW echoes and ringdown with the sensitivity curves of Advanced LIGO, ET, and CE in Figure 39. The detection of GW echoes with the third generation GW observatories are discussed in [119,120,174], and it may be possible to distinguish ECOs with | R | 0.3 from BHs with at 2 σ level when SNR 100 in ringdown, which would be possible for the third-generation GW detectors. The relative error on the reflectivity of would-be horizon is also investigated in [119,120,174], and the relative error for measurement of reflectivity in ground-based detectors is approximately given by
Δ R 1 R 0.5 × 8 ρ ringdown ,
where M = 30 M , ρ ringdown is the SNR in the ringdown phase, while the distance between the top of the angular momentum barrier and the would-be horizon is assumed to be longer than 50 M in the tortoise coordinate. For comparison, we note that the loudest detected BBH event, GW150914, has ρ ringdown 8 .
The detectability of GW echoes from failed supernovae, leading to the formation of BHs, with the third-generation GW observatories is also discussed in [130]. Calculating the SNR of GW spectrum consisting of echo and ringdown, ρ ringdown + echo , in the Boltzmann reflectivity model, the horizon distance D h , defined as the distance where ρ ringdown + echo = 8 , is estimated. Given the optimistic case in the Boltzmann reflectivity model, T H / T QH = e 15 ( a ¯ 1 ) , the horizon distance can be estimated as D h 10 Mpc for the Advanced LIGO at design sensitybity and D h 100 Mpc for the third-generation detectors such as ET and CE. Therefore, the authors in [130] argue that the searching for GW echoes, sourced by failed supernovae within our Galaxy and nearby galaxies, may be possible. However, in the case of T QH = T H , the horizon distance is less than or comparable with 10 Mpc and so the echo search with failed supernovae would be restricted to within the Local Group. For the comparison, the strain amplitude of GW echoes in T QH / T H = 1 and T QH / T H = e 15 ( a ¯ 1 ) are shown14 in Figure 40.

6.5. LISA

The Laser Interferometer Space Antena (LISA) is planned to be the first GW observatory in space. It will have three satellites separated by millions of kilometers and their orbits maintain near-equilateral triangular formation. LISA might enable us to reach high-precision detection of ringdown in SNR O ( 10 3 ) , which puts stronger constraints on the reflectivity of BHs [119,120].
Recently, a novel proposal to discriminate BH horizons based on the tidal heating was proposed in [175,176]. One of the main targets of the LISA mission is precision measurements of extreme-mass-ratio inspirals (EMRIs), in which the tidal heating could be important. The (partial) absorption of ECOs or BHs plays the role of dissipation at the surface, by which tides back-react on the orbital trajectory. It is argued that this tidal heating is responsible for a large dephasing between the orbits of a BH and ECO. This dephasing accumulates over the timescale of months and the accumulation speed is faster for a higher spin. The authors in [176] also found a proportionality relation between the dephasing δ ϕ and energy reflectivity | R | 2 .
To make use of this scheme to put strong constraints on the reflectivity of ECOs, one has to obtain accurate EMRI waveforms by properly taking into account the tidal heating for orbiting objects, which may decrease systematic errors in data analysis.
Not only the tidal heating, but also the tidal deformability contributes to the GW Fourier phase and it can be characterized by the tidal Love number k. The Love number of ECO of mass M may scale as 1 / | log δ | , where δ r 0 r h , with r h is the BH horizon radius of mass M and r 0 is the radius of the ECO. So the k δ relation is
δ = r h e 1 / k ,
and assuming this relation, one can infer the near-horizon structure characterized by δ from the measurement of the Love number k. For instance, if the Love number of the order of k 10 2 is measured by LISA from a supermassive BH binary signal, leading to the formation of a BH of M 10 6 M , it yields the resolution of δ l Pl .
However, the authors in [177] point out that the statistical and quantum mechanical uncertainties in measurements of near-horizon lead to some difficulty to measure δ precisely. The former one comes from the fact that the statistical uncertainty in δ is proportional to 1 / k , and the inferred value of k, where the inferred value of δ is comparable with its statistical uncertainty, is around k 0.2 . Therefore, any inferred value of δ , derived from k that is smaller than ∼0.2, would be dominated by the statistical uncertainty. The latter one comes from the uncertainty principle in quantum mechanics. Once precisely measuring δ l Pl , it may lead to the uncertainty in the mass of the ECO, which then leads to the uncertainty in the binding energy. This results in the uncertainty in the orbital and GW frequencies, which means that one cannot measure δ precisely if it is much shorter than l Pl .

6.6. Pulsar Timing Arrays

Following their first discovery in 1968 [178], over 2000 pulsars have now been detected by radio telescopes across the world. The pulsars’ intrinsic properties, as well as propagation effects in the interstellar medium, can influence the arrival times of pulsar pulses. Therefore, pulsar timing arrays (PTA) can be used as a detection tool for BH binaries [179], and thus, might be used to detect singatures of echoes from quantum BHs. In particular, millisecond pulsars stand out for their unparalleled stability (comparable to atomic clocks!) without being subject to starquakes and accretion. To give an explicit example, we show the spectrum of GW echoes predicted by the Boltzmann reflectivity model [15,29] with the sensitivity curve of International Pulsar Timing Array (IPTA) and Square Kilometre Array (SKA) (Figure 41). The lower curve in Figure 41 is for a 3 9 M BH merger at D o = 1 Gpc. Given that this mass is comparable to that of M87 supermassive BH, located at 16 Mpc, we expect 2 × 10 5 of such BHs at < Gpc. Assuming that each BH merges once every Hubble time 10 10 years, and that echoes last for 20 years (from simple mass scaling), the chances of detecting such a loud event with PTAs at any time is 0.1%. However, fainter events will be more prevalent as their number increases as SNR 3 / 2 from volume scaling. Furthermore, increase in supermassive BH merger activity observed at high redshifts shall boost this statistics.
PTAs are anticipated to detect the low frequency GW signal from supermassive BBH within the next few years [179]. We expect that the first GW detection will be a stochastic background of supermassive BH binaries. With any luck, this shall lead to new insights into the nature of quantum BHs and gravity.

7. Final Word

In this review article, we provided a comprehensive overview of the theoretical motivations for why quantum black holes in our universe may have different observable properties, in contrast to their classical counterparts in Einstein’s theory of general relativity. The most prominent and potentially observable smoking gun for these quantum black holes comes in the form of gravitational wave echoes, which have been the subject of intense theoretical and observational scrutiny over the past few years. We provided a concise account of theoretical predictions, as well as the exciting and confusing state of observational searches for echoes in LIGO/Virgo observations. We closed by article by our vision of the future of “Quantum Black Holes in the Sky”, via a synergy of statistical methodology, quantum gravity, and numerical relativity, and in light of the next generation of gravitational wave observatories.
While this review article focuses on the gravitational wave echoes, as arguably the most concrete and promising signature of quantum black holes, other possible observable signatures can be (and should be) explored. For example, interactions of photons or neutrinos with near-horizon quantum structure could lead to signatures in radio images in Event Horizon Telescope observations [180], or ultra high energy neutrinos in Ice Cube observatory [181], respectively. However, these signals will be suppressed if Boltzmann reflectivity is assumed, as they have ħ ω k T H . Another alternative to echoes may come through non-localities in non-violent unitarization, which would be observable far from the horizon (see Section 3.9). However, it is arguably difficult to pin down concrete predictions in this scenario. To conclude, the world of Quantum Black Holes remains a wide open and largely uncharted territory, spanning from the dark corners of obscure mathematical structures to the nitty-gritty details of gravitational wave detector noise. It also holds the promise to crack the century-old puzzle of quantum gravity, and yet be imminently testable in the next few years. Therefore, the study of “Quantum Black Holes in the Sky” remains extremely exciting, active, and confusing, and is bound to provide us with new surprises in the new decade, and beyond.

Author Contributions

Writing—original draft: J.A., N.O., and Q.W.; writing—review & editing: N.A.; supervision: N.A. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by Max-Planck-Institut für Gravitationsphysik (supporting J.A.), University of Waterloo (supporting N.A. and Q.W.), Natural Sciences and Engineering Research Council of Canada (NSERC; supporting N.A. and Q.W.), the JSPS Overseas Research Fellowships (supporting N.O.), and the Perimeter Institute for Theoretical Physics (supporting N.A., N.O., and Q.W.). Research at Perimeter Institute is supported in part by the Government of Canada through the Department of Innovation, Science and Economic Development Canada and by the Province of Ontario through the Ministry of Colleges and Universities.

Acknowledgments

We would like to thank Michael Balogh, Ofek Birnholtz, Avery Broderick, Vitor Cardoso, Ramit Dey, Hannah Dykaar, Will East, Steve Giddings, Bob Holdom, Badri Krishnan, Lam Hui, Luis Lehner, Luis Longo, Samir Mathur, Emil Mottola, Shinji Mukohyama, Rob Myers, Ramesh Narayan, Alex Nielsen, Paolo Pani, Joe Polchinski (RIP), Chanda Prescod-Weinstein, Jing Ren, Krishan Saraswat, Rafael Sorkin, Daichi Tsuna, Yasaman Yazdi, Huan Yang, Aaron Zimmerman, and many others for discussions and/or collaborations about Quantum Black Holes over the past decade. This research has made use of data, software and/or web tools obtained from the Gravitational Wave Open Science Center (https://www.gw-openscience.org), a service of LIGO Laboratory, the LIGO Scientific Collaboration and the Virgo Collaboration. LIGO is funded by the U.S. National Science Foundation. Virgo is funded by the French Centre National de Recherche Scientifique (CNRS), the Italian Instituto Nazionale della Fisica Nucleare (INFN) and the Dutch Nikhef, with contributions by Polish and Hungarian institutes.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Energy Density and Flux Seen by a Radially Free Falling Observer in Mock Fuzzball

We start with the geodesic of a radially falling observer in the mock fuzzball.
u μ u μ = 1 = ( 1 2 M r + b ) ( d t d τ ) 2 1 1 2 M r ( d r d τ ) 2 ,
1 = ( 1 2 M r + b ) ( d t d τ ) 2 ,
( 1 2 M r + b ) d t d τ = constant .
where τ is the proper time of the observer. θ and φ vanish with a radially observer. Equation (A1) is from the line element, Equation (A2) is initial condition of Equation (A1): the observer rest at infinite. Equation (A3) is energy conservation. Equations (A1)–(A3) give four-velocity:
u μ = ( r b r 2 M + r , 2 M ( r 2 M ) r ( b r 2 M + r ) , 0 , 0 ) .
Then we consider area vector a μ = A ( u 1 , u 0 , 0 , 0 ) which is normal to four-velocity a μ u μ = 0 . We determine A = 1 2 M / r + b 1 2 M / r by normalization of the area vector a μ a μ = 1 . In the end, the energy density and flux are defined as:
U = T μ ν u μ u ν = 4 b M 2 r 2 ( b r 2 M + r ) 2 ,
F = T μ ν u μ a ν = b 2 M r 3 / 2 ( b r 2 M + r ) 2 .

References

  1. Israel, W. Event horizons in static electrovac space-times. Commun. Math. Phys. 1968, 8, 245–260. [Google Scholar] [CrossRef]
  2. Carter, B. Axisymmetric Black Hole Has Only Two Degrees of Freedom. Phys. Rev. Lett. 1971, 26, 331–333. [Google Scholar] [CrossRef]
  3. Kokkotas, K.D.; Schmidt, B.G. Quasinormal modes of stars and black holes. Living Rev. Relativ. 1999, 2, 2. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Isi, M.; Giesler, M.; Farr, W.M.; Scheel, M.A.; Teukolsky, S.A. Testing the no-hair theorem with GW150914. Phys. Rev. Lett. 2019, 123, 111102. [Google Scholar] [CrossRef] [Green Version]
  5. Völkel, S.H.; Kokkotas, K.D. Ultra Compact Stars: Reconstructing the Perturbation Potential. Class. Quantum Gravity 2017, 34, 175015. [Google Scholar] [CrossRef] [Green Version]
  6. Völkel, S.H.; Kokkotas, K.D. Wormhole Potentials and Throats from Quasi-Normal Modes. Class. Quantum Gravity 2018, 35, 105018. [Google Scholar] [CrossRef] [Green Version]
  7. Bhagwat, S.; Forteza, X.J.; Pani, P.; Ferrari, V. Ringdown overtones, black hole spectroscopy and, no-hair theorem tests. arXiv 2019, arXiv:1910.08708. [Google Scholar] [CrossRef] [Green Version]
  8. Almheiri, A.; Marolf, D.; Polchinski, J.; Stanford, D.; Sully, J. An Apologia for Firewalls. J. High Energy Phys. 2013, 18. [Google Scholar] [CrossRef] [Green Version]
  9. Prescod-Weinstein, C.; Afshordi, N.; Balogh, M.L. Stellar Black Holes and the Origin of Cosmic Acceleration. Phys. Rev. D 2009, 80, 043513. [Google Scholar] [CrossRef] [Green Version]
  10. Cardoso, V.; Franzin, E.; Pani, P. Is the gravitational-wave ringdown a probe of the event horizon? Phys. Rev. Lett. 2016, 116, 171101. [Google Scholar] [CrossRef]
  11. Mazur, P.O.; Mottola, E. Gravitational vacuum condensate stars. Proc. Nat. Acad. Sci. USA 2004, 101, 9545–9550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Lunin, O.; Mathur, S.D. AdS/CFT duality and the black hole information paradox. Nucl. Phys. 2002, B623, 342–394. [Google Scholar]
  13. Holdom, B.; Ren, J. Not quite a black hole. Phys. Rev. 2017, D95, 084034. [Google Scholar] [CrossRef] [Green Version]
  14. Oshita, N.; Afshordi, N. Probing microstructure of black hole spacetimes with gravitational wave echoes. Phys. Rev. 2019, D99, 044002. [Google Scholar] [CrossRef] [Green Version]
  15. Oshita, N.; Wang, Q.; Afshordi, N. On Reflectivity of Quantum Black Hole Horizons. arXiv 2019, arXiv:1905.00464. [Google Scholar]
  16. Cardoso, V.; Hopper, S.; Macedo, C.F.B.; Palenzuela, C.; Pani, P. Gravitational-wave signatures of exotic compact objects and of quantum corrections at the horizon scale. Phys. Rev. 2016, D94, 084031. [Google Scholar] [CrossRef] [Green Version]
  17. Cardoso, V.; Pani, P. Testing the nature of dark compact objects: A status report. Living Rev. Relativ. 2019, 22, 4. [Google Scholar] [CrossRef] [Green Version]
  18. Abedi, J.; Dykaar, H.; Afshordi, N. Echoes from the Abyss: Tentative evidence for Planck-scale structure at black hole horizons. Phys. Rev. 2017, D96, 082004. [Google Scholar]
  19. Uchikata, N.; Nakano, H.; Narikawa, T.; Sago, N.; Tagoshi, H.; Tanaka, T. Searching for black hole echoes from the LIGO-Virgo Catalog GWTC-1. arXiv 2019, arXiv:1906.00838. [Google Scholar] [CrossRef] [Green Version]
  20. Abedi, J.; Afshordi, N. Echoes from the Abyss: A highly spinning black hole remnant for the binary neutron star merger GW170817. arXiv 2018, arXiv:1803.10454. [Google Scholar] [CrossRef] [Green Version]
  21. Conklin, R.S.; Holdom, B.; Ren, J. Gravitational wave echoes through new windows. arXiv 2017, arXiv:1712.06517. [Google Scholar] [CrossRef] [Green Version]
  22. Holdom, B. Not quite black holes at LIGO. arXiv 2019, arXiv:1909.11801. [Google Scholar]
  23. Westerweck, J.; Nielsen, A.; Fischer-Birnholtz, O.; Cabero, M.; Capano, C.; Dent, T.; Krishnan, B.; Meadors, G.; Nitz, A.H. Low significance of evidence for black hole echoes in gravitational wave data. arXiv 2017, arXiv:1712.09966. [Google Scholar] [CrossRef] [Green Version]
  24. Nielsen, A.B.; Capano, C.D.; Birnholtz, O.; Westerweck, J. Parameter estimation and statistical significance of echoes following black hole signals in the first Advanced LIGO observing run. Phys. Rev. 2019, D99, 104012. [Google Scholar] [CrossRef] [Green Version]
  25. Salemi, F.; Milotti, E.; Prodi, G.A.; Vedovato, G.; Lazzaro, C.; Tiwari, S.; Vinciguerra, S.; Drago, M.; Klimenko, S. A wider look at the gravitational-wave transients from GWTC-1 using an unmodeled reconstruction method. arXiv 2019, arXiv:1905.09260. [Google Scholar] [CrossRef] [Green Version]
  26. Ashton, G.; Birnholtz, O.; Cabero, M.; Capano, C.; Dent, T.; Krishnan, B.; Meadors, G.D.; Nielsen, A.B.; Nitz, A.; Westerweck, J. Comments on: “Echoes from the abyss: Evidence for Planck-scale structure at black hole horizons”. arXiv 2016, arXiv:1612.05625. [Google Scholar]
  27. Abedi, J.; Dykaar, H.; Afshordi, N. Echoes from the Abyss: The Holiday Edition! arXiv 2017, arXiv:1701.03485. [Google Scholar]
  28. Abedi, J.; Dykaar, H.; Afshordi, N. Comment on: “Low significance of evidence for black hole echoes in gravitational wave data”. arXiv 2018, arXiv:1803.08565. [Google Scholar]
  29. Wang, Q.; Oshita, N.; Afshordi, N. Echoes from Quantum Black Holes. Phys. Rev. D 2020, 101, 024031. [Google Scholar] [CrossRef] [Green Version]
  30. Poisson, E. A Relativist’s Toolkit: The Mathematics of Black-Hole Mechanics; Cambridge University Press: Cambridge, MA, USA, 2009. [Google Scholar]
  31. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. GW150914: The Advanced LIGO Detectors in the Era of First Discoveries. Phys. Rev. Lett. 2016, 116, 131103. [Google Scholar] [CrossRef] [Green Version]
  32. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. Binary Black Hole Mergers in the first Advanced LIGO Observing Run. Phys. Rev. 2016, X6, 041015. [Google Scholar]
  33. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. Observation of Gravitational Waves from a Binary Black Hole Merger. Phys. Rev. Lett. 2016, 116, 061102. [Google Scholar] [CrossRef] [PubMed]
  34. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. GW151226: Observation of Gravitational Waves from a 22-Solar-Mass Binary Black Hole Coalescence. Phys. Rev. Lett. 2016, 116, 241103. [Google Scholar] [CrossRef] [PubMed]
  35. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. GW170104: Observation of a 50-Solar-Mass Binary Black Hole Coalescence at Redshift 0.2. Phys. Rev. Lett. 2017, 118, 221101. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. GW170814: A Three-Detector Observation of Gravitational Waves from a Binary Black Hole Coalescence. Phys. Rev. Lett. 2017, 119, 141101. [Google Scholar] [CrossRef] [Green Version]
  37. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. GW170817: Observation of Gravitational Waves from a Binary Neutron Star Inspiral. Phys. Rev. Lett. 2017, 119, 161101. [Google Scholar] [CrossRef] [Green Version]
  38. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. GW170608: Observation of a 19-solar-mass Binary Black Hole Coalescence. Astrophys. J. 2017, 851, L35. [Google Scholar] [CrossRef]
  39. Kerr, R.P. Gravitational Field of a Spinning Mass as an Example of Algebraically Special Metrics. Phys. Rev. Lett. 1963, 11, 237–238. [Google Scholar]
  40. Poisson, E.; Israel, W. Inner-horizon instability and mass inflation in black holes. Phys. Rev. Lett. 1989, 63, 1663–1666. [Google Scholar] [CrossRef]
  41. Bekenstein, J.D. Black holes and the second law. Lett. Nuovo Cim. 1972, 4, 737–740. [Google Scholar] [CrossRef]
  42. Bekenstein, J.D. Black holes and entropy. Phys. Rev. 1973, D7, 2333–2346. [Google Scholar] [CrossRef]
  43. Bekenstein, J.D. Generalized second law of thermodynamics in black hole physics. Phys. Rev. 1974, D9, 3292–3300. [Google Scholar]
  44. Gibbons, G.W.; Hawking, S.W. Action Integrals and Partition Functions in Quantum Gravity. Phys. Rev. 1977, D15, 2752–2756. [Google Scholar] [CrossRef]
  45. Hawking, S.W. Quantum Gravity and Path Integrals. Phys. Rev. 1978, D18, 1747–1753. [Google Scholar]
  46. Hawking, S.W. Black hole explosions. Nature 1974, 248, 30–31. [Google Scholar] [CrossRef]
  47. Hawking, S.W. Particle Creation by Black Holes. Commun. Math. Phys. 1975, 43, 199–220. [Google Scholar] [CrossRef]
  48. Hartle, J.B.; Hawking, S.W. Path-integral derivation of black-hole radiance. Phys. Rev. D 1976, 13, 2188–2203. [Google Scholar] [CrossRef]
  49. Rovelli, C. Black hole entropy from loop quantum gravity. Phys. Rev. Lett. 1996, 77, 3288–3291. [Google Scholar] [CrossRef] [Green Version]
  50. Bardeen, J.M.; Carter, B.; Hawking, S.W. The Four laws of black hole mechanics. Commun. Math. Phys. 1973, 31, 161–170. [Google Scholar] [CrossRef]
  51. Thorne, K.S.; Price, R.H.; Macdonald, D.A. (Eds.) Black Holes: The Membrane Paradigm; Yale University Press: New Haven, CT, USA, 1986. [Google Scholar]
  52. Regge, T.; Wheeler, J.A. Stability of a Schwarzschild singularity. Phys. Rev. 1957, 108, 1063–1069. [Google Scholar] [CrossRef]
  53. Jacobson, T.; Mohd, A.; Sarkar, S. Membrane paradigm for Einstein-Gauss-Bonnet gravity. Phys. Rev. 2017, D95, 064036. [Google Scholar]
  54. Donoghue, J.F. General relativity as an effective field theory: The leading quantum corrections. Phys. Rev. 1994, D50, 3874–3888. [Google Scholar] [CrossRef] [Green Version]
  55. Polchinski, J. String Theory. Vol. 1: An Introduction to the Bosonic String; Cambridge Monographs on Mathematical Physics; Cambridge University Press: Cambridge, MA, USA, 2007. [Google Scholar]
  56. Ashtekar, A.; Lewandowski, J. Background independent quantum gravity: A Status report. Class. Quantum Gravity 2004, 21, R53. [Google Scholar] [CrossRef]
  57. Bombelli, L.; Lee, J.; Meyer, D.; Sorkin, R. Space-Time as a Causal Set. Phys. Rev. Lett. 1987, 59, 521–524. [Google Scholar] [CrossRef] [PubMed]
  58. Ambjorn, J.; Jurkiewicz, J.; Loll, R. Emergence of a 4-D world from causal quantum gravity. Phys. Rev. Lett. 2004, 93, 131301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Niedermaier, M.; Reuter, M. The Asymptotic Safety Scenario in Quantum Gravity. Living Rev. Relativ. 2006, 9, 5–173. [Google Scholar] [CrossRef] [Green Version]
  60. Holdom, B.; Ren, J. QCD analogy for quantum gravity. Phys. Rev. 2016, D93, 124030. [Google Scholar]
  61. Anselmi, D. On the quantum field theory of the gravitational interactions. J. High Energy Phys. 2017, 86. [Google Scholar] [CrossRef] [Green Version]
  62. Horava, P. Quantum Gravity at a Lifshitz Point. Phys. Rev. 2009, D79, 084008. [Google Scholar]
  63. Mathur, S.D. Emission rates, the correspondence principle and the information paradox. Nucl. Phys. 1998, B529, 295–320. [Google Scholar] [CrossRef] [Green Version]
  64. Lunin, O.; Mathur, S.D. Statistical interpretation of Bekenstein entropy for systems with a stretched horizon. Phys. Rev. Lett. 2002, 88, 211303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Mathur, S.D. The Fuzzball proposal for black holes: An Elementary review. Fortschritte Phys. 2005, 53, 793–827. [Google Scholar] [CrossRef] [Green Version]
  66. Mathur, S.D. Fuzzballs and the information paradox: A Summary and conjectures. arXiv 2008, arXiv:0810.4525. [Google Scholar]
  67. Mathur, S.D.; Turton, D. Comments on black holes I: The possibility of complementarity. J. High Energy Phys. 2014, 34. [Google Scholar] [CrossRef] [Green Version]
  68. Braunstein, S.L.; Pirandola, S.; Zyczkowski, K. Better Late than Never: Information Retrieval from Black Holes. Phys. Rev. Lett. 2013, 110, 101301. [Google Scholar] [CrossRef]
  69. Almheiri, A.; Marolf, D.; Polchinski, J.; Sully, J. Black Holes: Complementarity or Firewalls? J. High Energy Phys. 2013, 62. [Google Scholar] [CrossRef] [Green Version]
  70. Barceló, C.; Carballo-Rubio, R.; Garay, L.J. Where does the physics of extreme gravitational collapse reside? Universe 2016, 2, 7. [Google Scholar] [CrossRef] [Green Version]
  71. Kawai, H.; Yokokura, Y. A Model of Black Hole Evaporation and 4D Weyl Anomaly. Universe 2017, 3, 51. [Google Scholar] [CrossRef] [Green Version]
  72. Giddings, S.B. Gravitational wave tests of quantum modifications to black hole structure—With post-GW150914 update. Class. Quantum Gravity 2016, 33, 235010. [Google Scholar] [CrossRef] [Green Version]
  73. Casadio, R.; da Rocha, R. Stability of the graviton Bose–Einstein condensate in the brane-world. Phys. Lett. 2016, B763, 434–438. [Google Scholar] [CrossRef] [Green Version]
  74. Hawking, S.W. Breakdown of Predictability in Gravitational Collapse. Phys. Rev. 1976, D14, 2460–2473. [Google Scholar] [CrossRef]
  75. Fabbri, A.; Navarro-Salas, J. Modeling Black Hole Evaporation; Imperial College Press: London, UK, 2005. [Google Scholar]
  76. Aharonov, Y.; Casher, A.; Nussinov, S. The Unitarity Puzzle and Planck Mass Stable Particles. Phys. Lett. 1987, B191, 51. [Google Scholar] [CrossRef] [Green Version]
  77. Banks, T.; O’Loughlin, M.; Strominger, A. Black hole remnants and the information puzzle. Phys. Rev. 1993, D47, 4476–4482. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Giddings, S.B. Quantum mechanics of black holes. In Proceedings of the Summer School in High-Energy Physics and Cosmology, Trieste, Italy, 13 June–29 July 1994; pp. 530–574. [Google Scholar]
  79. Adler, R.J.; Chen, P.; Santiago, D.I. The Generalized uncertainty principle and black hole remnants. Gen. Relativ. Gravit. 2001, 33, 2101–2108. [Google Scholar] [CrossRef]
  80. Susskind, L.; Thorlacius, L.; Uglum, J. The Stretched horizon and black hole complementarity. Phys. Rev. 1993, D48, 3743–3761. [Google Scholar] [CrossRef] [Green Version]
  81. Susskind, L.; Thorlacius, L. Gedanken experiments involving black holes. Phys. Rev. 1994, D49, 966–974. [Google Scholar] [CrossRef] [Green Version]
  82. Page, D.N. Average entropy of a subsystem. Phys. Rev. Lett. 1993, 71, 1291–1294. [Google Scholar] [CrossRef] [Green Version]
  83. Page, D.N. Information in black hole radiation. Phys. Rev. Lett. 1993, 71, 3743–3746. [Google Scholar] [CrossRef] [Green Version]
  84. Page, D.N. Time Dependence of Hawking Radiation Entropy. J. Cosmol. Astropart. Phys. 2013, 1309, 028. [Google Scholar] [CrossRef]
  85. Papadodimas, K.; Raju, S. An Infalling Observer in AdS/CFT. J. High Energy Phys. 2013, 10, 212. [Google Scholar] [CrossRef] [Green Version]
  86. Maldacena, J.; Susskind, L. Cool horizons for entangled black holes. Fortsch. Phys. 2013, 61, 781–811. [Google Scholar] [CrossRef] [Green Version]
  87. Susskind, L. New Concepts for Old Black Holes. arXiv 2013, arXiv:1311.3335. [Google Scholar]
  88. Bousso, R. Complementarity Is Not Enough. Phys. Rev. 2013, D87, 124023. [Google Scholar] [CrossRef] [Green Version]
  89. Yoshida, B. Firewalls vs. Scrambling. J. High Energy Phys. 2019, 10, 132. [Google Scholar] [CrossRef] [Green Version]
  90. Oshita, N. Resolution to the firewall paradox: The black hole information paradox and highly squeezed interior quantum fluctuations. Class. Quantum Gravity 2017, 34, 195002. [Google Scholar] [CrossRef] [Green Version]
  91. Carballo-Rubio, R. Stellar equilibrium in semiclassical gravity. Phys. Rev. Lett. 2018, 120, 061102. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  92. Sakai, N.; Saida, H.; Tamaki, T. Gravastar Shadows. Phys. Rev. 2014, D90, 104013. [Google Scholar]
  93. Pani, P.; Berti, E.; Cardoso, V.; Chen, Y.; Norte, R. Gravitational wave signatures of the absence of an event horizon. I. Nonradial oscillations of a thin-shell gravastar. Phys. Rev. 2009, D80, 124047. [Google Scholar] [CrossRef] [Green Version]
  94. Maldacena, J.M. The Large N limit of superconformal field theories and supergravity. Int. J. Theor. Phys. 1999, 38, 1113–1133. [Google Scholar] [CrossRef] [Green Version]
  95. Giusto, S.; Saxena, A.; Ross, S.F. Non-supersymmetric microstates of the D1-D5-KK system. J. High Energy Phys. 2007, 65. [Google Scholar] [CrossRef] [Green Version]
  96. Wang, Q. Physical Structure of “Horizons” in (Mock) Fuzzballs. Master’s Thesis, Perimeter Institute and University of Waterloo, Waterloo, ON, Canada, 2016. [Google Scholar]
  97. Damour, T.; Solodukhin, S.N. Wormholes as black hole foils. Phys. Rev. 2007, D76, 024016. [Google Scholar] [CrossRef] [Green Version]
  98. Ovalle, J.; Casadio, R.; da Rocha, R.; Sotomayor, A. Anisotropic solutions by gravitational decoupling. Eur. Phys. J. 2018, C78, 122. [Google Scholar]
  99. Ovalle, J.; Casadio, R.; Rocha, R.D.; Sotomayor, A.; Stuchlik, Z. Black holes by gravitational decoupling. Eur. Phys. J. 2018, C78, 960. [Google Scholar]
  100. Page, D.N. Thermal stress tensors in static Einstein spaces. Phys. Rev. D 1982, 25, 1499. [Google Scholar] [CrossRef]
  101. Parker, L.; Toms, D. Quantum Field Theory in Curved Spacetime: Quantized Fields And Gravity; Cambridge University Press: Cambridge, MA, USA, 2009. [Google Scholar]
  102. Afshordi, N. Gravitational Aether and the thermodynamic solution to the cosmological constant problem. arXiv 2008, arXiv:0807.2639. [Google Scholar]
  103. Aslanbeigi, S.; Robbers, G.; Foster, B.Z.; Kohri, K.; Afshordi, N. Phenomenology of Gravitational Aether as a solution to the Old Cosmological Constant Problem. Phys. Rev. 2011, D84, 103522. [Google Scholar]
  104. Penrose, R. Gravitational Collapse and Space-Time Singularities. Phys. Rev. Lett. 1965, 14, 57–59. [Google Scholar] [CrossRef] [Green Version]
  105. Stelle, K.S. Renormalization of Higher Derivative Quantum Gravity. Phys. Rev. 1977, D16, 953–969. [Google Scholar] [CrossRef]
  106. Voronov, B.L. On renormalization of R2-gravitation. Yad. Fiz. (Sov. J. Nucl. Phys.) 1984, 39, 998. [Google Scholar]
  107. Fradkin, E.S.; Tseytlin, A.A. Renormalizable asymptotically free quantum theory of gravity. Nucl. Phys. 1982, B201, 469–491. [Google Scholar] [CrossRef]
  108. Avramidi, I.G.; Barvinsky, A.O. Asyptotic freedom in higher derivative quantum gravity. Phys. Lett. 1985, 159B, 269–274. [Google Scholar] [CrossRef]
  109. Holdom, B. A ghost and a naked singularity; facing our demons. Presented at the CERN Workshop “Scale Invariance in Particle Physics and Cosmology”, Geneva, Switzerland, 28 January–1 February 2019. [Google Scholar]
  110. Ren, J. Anatomy of a Burning 2-2-hole. Phys. Rev. D 2019, 100, 124012. [Google Scholar] [CrossRef] [Green Version]
  111. Giddings, S.B. Nonviolent unitarization: Basic postulates to soft quantum structure of black holes. J. High Energy Phys. 2017, 47. [Google Scholar] [CrossRef] [Green Version]
  112. Bardeen, J.M. Interpreting the semi-classical stress-energy tensor in a Schwarzschild background, implications for the information paradox. arXiv 2018, arXiv:1808.08638. [Google Scholar]
  113. Bardeen, J.M. Models for the nonsingular transition of an evaporating black hole into a white hole. arXiv 2018, arXiv:1811.06683. [Google Scholar]
  114. Giddings, S.B. Black holes in the quantum universe. Phil. Trans. Roy. Soc. Lond. 2019, A377, 20190029. [Google Scholar] [CrossRef] [Green Version]
  115. Cardoso, V.; Foit, V.F.; Kleban, M. Gravitational wave echoes from black hole area quantization. arXiv 2019, arXiv:1902.10164. [Google Scholar] [CrossRef] [Green Version]
  116. Coates, A.; Völkel, S.H.; Kokkotas,, K.D. Spectral Lines of Quantized, Spinning Black Holes and their Astrophysical Relevance. Phys. Rev. Lett. 2019, 123, 171104. [Google Scholar] [CrossRef] [Green Version]
  117. Mark, Z.; Zimmerman, A.; Du, S.M.; Chen, Y. A recipe for echoes from exotic compact objects. Phys. Rev. 2017, D96, 084002. [Google Scholar] [CrossRef] [Green Version]
  118. Conklin, R.S.; Holdom, B. Gravitational Wave “Echo” Spectra. arXiv 2019, arXiv:1905.09370. [Google Scholar] [CrossRef] [Green Version]
  119. Testa, A.; Pani, P. Analytical template for gravitational-wave echoes: Signal characterization and prospects of detection with current and future interferometers. Phys. Rev. 2018, D98, 044018. [Google Scholar] [CrossRef] [Green Version]
  120. Maggio, E.; Testa, A.; Bhagwat, S.; Pani, P. Analytical model for gravitational-wave echoes from spinning remnants. Phys. Rev. 2019, D100, 064056. [Google Scholar] [CrossRef] [Green Version]
  121. Chandrasekhar, S.; Detweiler, S.L. Equations governing gravitational perturbations of the Kerr black-hole. Proc. R. Soc. Lond. 1976, A350, 165–174. [Google Scholar]
  122. Detweiler, S.L. Resonant oscillations of a rapidly rotating black hole. Proc. Roy. Soc. Lond. 1977, A352, 381–395. [Google Scholar]
  123. Chandrasekhar, S.; Detweiler, S.L. The quasi-normal modes of the Schwarzschild black hole. Proc. Roy. Soc. Lond. 1975, A344, 441–452. [Google Scholar]
  124. Casals, M.; Longo Micchi, L.F. Spectroscopy of extremal and near-extremal Kerr black holes. Phys. Rev. 2019, D99, 084047. [Google Scholar] [CrossRef] [Green Version]
  125. Mano, S.; Suzuki, H.; Takasugi, E. Analytic Solutions of the Teukolsky Equation and Their Low Frequency Expansions. Prog. Theor. Phys. 1996, 95, 1079–1096. [Google Scholar] [CrossRef]
  126. Mano, S.; Suzuki, H.; Takasugi, E. Analytic solutions of the Regge-Wheeler equation and the postMinkowskian expansion. Prog. Theor. Phys. 1996, 96, 549–566. [Google Scholar] [CrossRef]
  127. Sasaki, M.; Tagoshi, H. Analytic black hole perturbation approach to gravitational radiation. Living Rev. Relativ. 2003, 6, 6. [Google Scholar] [CrossRef] [Green Version]
  128. Casals, M.; Zimmerman, P. Perturbations of Extremal Kerr Spacetime: Analytic Framework and Late-time Tails. arXiv 2018, arXiv:1801.05830. [Google Scholar] [CrossRef] [Green Version]
  129. Starobinskiǐ, A.A. Amplification of waves during reflection from a rotating “black hole”. Sov. J. Exp. Theor. Phys. 1973, 37, 28. [Google Scholar]
  130. Oshita, N.; Tsuna, D.; Afshordi, N. Quantum Black Hole Seismology I: Echoes, Ergospheres, and Spectra. arXiv 2020, arXiv:2001.11642. [Google Scholar]
  131. Berti, E.; Cardoso, V.; Will, C.M. On gravitational-wave spectroscopy of massive black holes with the space interferometer LISA. Phys. Rev. 2006, D73, 064030. [Google Scholar]
  132. Kokkotas, K.D. Pulsating relativistic stars. Relativistic gravitation and gravitational radiation. In Proceedings of the School of Physics, Les Houches, France, 26 September–6 October 1995; pp. 89–102. [Google Scholar]
  133. Tominaga, K.; Saijo, M.; Maeda, K.I. Gravitational waves from a test particle scattered by a neutron star: Axial mode case. Phys. Rev. 1999, D60, 024004. [Google Scholar] [CrossRef] [Green Version]
  134. Ferrari, V.; Kokkotas, K.D. Scattering of particles by neutron stars: Time evolutions for axial perturbations. Phys. Rev. D 2000, 62, 107504. [Google Scholar] [CrossRef] [Green Version]
  135. Srinivasan, K.; Padmanabhan, T. Particle production and complex path analysis. Phys. Rev. 1999, D60, 024007. [Google Scholar] [CrossRef] [Green Version]
  136. Vanzo, L.; Acquaviva, G.; Di Criscienzo, R. Tunnelling Methods and Hawking’s radiation: Achievements and prospects. Class. Quantum Gravity 2011, 28, 183001. [Google Scholar] [CrossRef]
  137. Padmanabhan, T. Thermality of the Rindler horizon: A simple derivation from the structure of the inertial propagator. Phys. Rev. 2019, D100, 045024. [Google Scholar] [CrossRef] [Green Version]
  138. Parentani, R. Quantum metric fluctuations and Hawking radiation. Phys. Rev. 2001, D63, 041503. [Google Scholar] [CrossRef] [Green Version]
  139. Barrabes, C.; Frolov, V.P.; Parentani, R. Stochastically fluctuating black hole geometry, Hawking radiation and the transPlanckian problem. Phys. Rev. 2000, D62, 044020. [Google Scholar]
  140. Brout, R.; Gabriel, C.; Lubo, M.; Spindel, P. Minimal length uncertainty principle and the transPlanckian problem of black hole physics. Phys. Rev. 1999, D59, 044005. [Google Scholar]
  141. Liberati, S.; Maccione, L. Astrophysical constraints on Planck scale dissipative phenomena. Phys. Rev. Lett. 2014, 112, 151301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  142. Damour, T. Black-hole eddy currents. Phys. Rev. D 1978, 18, 3598–3604. [Google Scholar] [CrossRef]
  143. Frolov, A.V.; Frolov, V.P. Rigidly rotating zero-angular-momentum observer surfaces in the Kerr spacetime. Phys. Rev. 2014, D90, 124010. [Google Scholar] [CrossRef] [Green Version]
  144. Abedi, J.; Afshordi, N. Echoes from the Abyss: A Status Update. arXiv 2020, arXiv:2001.00821. [Google Scholar]
  145. Lo, R.K.L.; Li, T.G.F.; Weinstein, A.J. Template-based Gravitational-Wave Echoes Search Using Bayesian Model Selection. Phys. Rev. 2019, D99, 084052. [Google Scholar] [CrossRef] [Green Version]
  146. Tsang, K.W.; Rollier, M.; Ghosh, A.; Samajdar, A.; Agathos, M.; Chatziioannou, K.; Cardoso, V.; Khanna, G.; Van Den Broeck, C. A morphology-independent data analysis method for detecting and characterizing gravitational wave echoes. Phys. Rev. 2018, D98, 024023. [Google Scholar] [CrossRef] [Green Version]
  147. Tsang, K.W.; Ghosh, A.; Samajdar, A.; Chatziioannou, K.; Mastrogiovanni, S.; Agathos, M.; Van Den Broeck, C. A morphology-independent search for gravitational wave echoes in data from the first and second observing runs of Advanced LIGO and Advanced Virgo. arXiv 2019, arXiv:1906.11168. [Google Scholar] [CrossRef] [Green Version]
  148. Goodman, S.N. Of P-Values and Bayes: A Modest Proposal. Epidemiology 2001, 12, 295. [Google Scholar] [CrossRef]
  149. Gill, R.; Nathanail, A.; Rezzolla, L. When Did the Remnant of GW170817 Collapse to a Black Hole? Astrophys. J. 2019, 876, 139. [Google Scholar] [CrossRef] [Green Version]
  150. Abbott, B.P.; Jawahar, S.; Lockerbie, N.A.; Tokmakov, K.V.; LIGO Scientific Collaboration and Virgo Collaboration. Tests of general relativity with GW150914. Phys. Rev. Lett. 2016, 116, 221101. [Google Scholar] [CrossRef] [PubMed]
  151. Gravitational Wave Open Science Center. Available online: https://www.gw-openscience.org (accessed on 10 March 2020).
  152. Abbott, B.P.; Abbott, R.; Abbott, T.D.; Acernese, F.; Ackley, K.; Adams, C.; Adams, T.; Addesso, P.; Adhikari, R.X.; LIGO Scientific Collaboration and Virgo Collaboration; et al. Search for Post-merger Gravitational Waves from the Remnant of the Binary Neutron Star Merger GW170817. Astrophys. J. 2017, 851, L16. [Google Scholar] [CrossRef]
  153. Maggio, E.; Pani, P.; Ferrari, V. Exotic Compact Objects and How to Quench their Ergoregion Instability. arXiv 2017, arXiv:1703.03696. [Google Scholar] [CrossRef] [Green Version]
  154. Bueno, P.; Cano, P.A.; Goelen, F.; Hertog, T.; Vercnocke, B. Echoes of Kerr-like wormholes. Phys. Rev. 2018, D97, 024040. [Google Scholar] [CrossRef] [Green Version]
  155. Wang, Q.; Afshordi, N. Black hole echology: The observer’s manual. Phys. Rev. 2018, D97, 124044. [Google Scholar] [CrossRef] [Green Version]
  156. Coherent WaveBurst. Available online: https://gwburst.gitlab.io (accessed on 10 March 2020).
  157. Parameter Estimation Sample Release for GWTC-1. Available online: https://dcc.ligo.org/LIGO-P1800370/public (accessed on 10 March 2020).
  158. Line of Best Fit (Least Square Method). Available online: https://www.varsitytutors.com/hotmath/hotmath_help/topics/line-of-best-fit (accessed on 10 March 2020).
  159. Reconstructed Detector Responses for GW151012. Available online: https://www.atlas.aei.uni-hannover.de/work/salemi/reports/GWOSC_LH_GW151012_R1_org/ced/ced_1128678000_1200_GWOSC_LH_GW151012_R1_org_slag0_lag0_1_job1/L1H1_1128678900.375_1128678900.375/ (accessed on 10 March 2020).
  160. Reconstructed Detector Responses for GW151226. Available online: https://www.atlas.aei.uni-hannover.de/work/salemi/reports/GWOSC_LH_GW151226_R1_org/ced/ced_1135136238_1200_GWOSC_LH_GW151226_R1_org_slag0_lag0_1_job1/L1H1_1135136350.312_1135136350.312/ (accessed on 10 March 2020).
  161. Reconstructed Detector Responses for GW150914. Available online: https://www.atlas.aei.uni-hannover.de/work/salemi/reports/GWOSC_LH_GW150914_R1_org/ced/ced_1126258800_1200_GWOSC_LH_GW150914_R1_org_slag0_lag0_1_job1/L1H1_1126259462.125_1126259462.125/ (accessed on 10 March 2020).
  162. Buonanno, A.; Damour, T. Effective one-body approach to general relativistic two-body dynamics. Phys. Rev. 1999, D59, 084006. [Google Scholar]
  163. Damour, T. Introductory lectures on the Effective One Body formalism. Int. J. Mod. Phys. 2008, A23, 1130–1148. [Google Scholar]
  164. Okounkova, M.; Stein, L.C.; Scheel, M.A.; Teukolsky, S.A. Numerical binary black hole collisions in dynamical Chern-Simons gravity. Phys. Rev. 2019, D100, 104026. [Google Scholar] [CrossRef] [Green Version]
  165. Okounkova, M.; Stein, L.C.; Moxon, J.; Scheel, M.A.; Teukolsky, S.A. Numerical relativity simulation of GW150914 beyond general relativity. arXiv 2019, arXiv:1911.02588. [Google Scholar]
  166. Ryu, S.; Takayanagi, T. Aspects of Holographic Entanglement Entropy. J. High Energy Phys. 2006, 45. [Google Scholar] [CrossRef] [Green Version]
  167. Sekino, Y.; Susskind, L. Fast Scramblers. J. High Energy Phys. 2008, 65. [Google Scholar] [CrossRef] [Green Version]
  168. Maldacena, J.; Shenker, S.H.; Stanford, D. A bound on chaos. J. High Energy Phys. 2016, 106. [Google Scholar] [CrossRef] [Green Version]
  169. Saraswat, K.; Afshordi, N. Quantum Nature of Black Holes: Fast Scrambling versus Echoes. arXiv 2019, arXiv:1906.02653. [Google Scholar]
  170. Castro, A.; Maloney, A.; Strominger, A. Hidden Conformal Symmetry of the Kerr Black Hole. Phys. Rev. 2010, D82, 024008. [Google Scholar] [CrossRef] [Green Version]
  171. Dey, R.; Chakraborty, S.; Afshordi, N. Echoes from the braneworld black holes. arXiv 2020, arXiv:2001.01301. [Google Scholar]
  172. Punturo, M.; Abernathy, M.; Acernese, F.; Allen, B.; Andersson, N.; Arun, K.; Barone, F.; Barr, B.; Barsuglia, M.; Beker, M.; et al. The Einstein Telescope: A third-generation gravitational wave observatory. Class. Quantum Gravity 2010, 27, 194002. [Google Scholar] [CrossRef]
  173. Dwyer, S.; Sigg, D.; Ballmer, S.W.; Barsotti, L.; Mavalvala, N.; Evans, M. Gravitational wave detector with cosmological reach. Phys. Rev. D 2015, 91, 082001. [Google Scholar] [CrossRef] [Green Version]
  174. Maselli, A.; Völkel, S.H.; Kokkotas, K.D. Parameter estimation of gravitational wave echoes from exotic compact objects. Phys. Rev. 2017, D96, 064045. [Google Scholar] [CrossRef] [Green Version]
  175. Datta, S.; Bose, S. Probing the nature of central objects in extreme-mass-ratio inspirals with gravitational waves. Phys. Rev. 2019, D99, 084001. [Google Scholar] [CrossRef] [Green Version]
  176. Datta, S.; Brito, R.; Bose, S.; Pani, P.; Hughes, S.A. Tidal heating as a discriminator for horizons in extreme mass ratio inspirals. arXiv 2019, arXiv:1910.07841. [Google Scholar] [CrossRef] [Green Version]
  177. Addazi, A.; Marciano, A.; Yunes, N. Can we probe Planckian corrections at the horizon scale with gravitational waves? Phys. Rev. Lett. 2019, 122, 081301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  178. Hewish, A.; Bell, S.J.; Pilkington, J.D.H.; Scott, P.F.; Collins, R.A. Observation of a rapidly pulsating radio source. Nature 1968, 217, 709–713. [Google Scholar] [CrossRef]
  179. Hobbs, G.; Archibald, A.; Arzoumanian, D.; Backer, M.; Bailes, N.D.R.; Bhat, M.; Burgay, S.; Burke-Spolaor, D.; Champio, D.; Cognard, W.; et al. The international pulsar timing array project: Using pulsars as a gravitational wave detector. Class. Quantum Gravity 2010, 27, 084013. [Google Scholar] [CrossRef]
  180. Rummel, M.; Burgess, C.P. Constraining Fundamental Physics with the Event Horizon Telescope. arXiv 2019, arXiv:2001.00041. [Google Scholar]
  181. Afshordi, N.; Yazdi, Y.K. ‘Firewall’ phenomenology with astrophysical neutrinos. Class. Quantum Gravity 2016, 33, 235017. [Google Scholar] [CrossRef] [Green Version]
1.
It turns out that this toy model spacetime coincides with the one proposed earlier by [97].
2.
We note that modified phenomenological BH spacetimes with anisotropic stresses that can remove singulairities have been previously discussed, e.g., in [98,99].
3.
Recently, it was found that the CD equation with ( i , j ) = ( + 1 , ± 1 ) is related to ( 1 , 1 ) by the Darboux transformation.
4.
Please note that in order to set the initial conditions of the QNMs of the quantum BH cavity, we choose a superposition that reproduces the time-evolution of the dominant QNM of the classical BH for t Δ t echo .
5.
In [18], 2-tailed gaussian probability assigned to significance, e.g., p-value = 68 % and 95 % correspond to 1 σ and 2 σ respectively.
6.
There are two versions of noise subtraction in LIGO open data, called C01 and C02 [151]
7.
Additional consideration of the best fit initial phase has given negative results and is discussed in negative result part below Section 5.3.3. Therefore, this plot only justifies Uchikata et al. [19] reevaluation of ADA search [18] in Figure 17, while belongs to the discussion in part Section 5.3.3.
8.
In this paper, a 1-tailed gaussian probability to assign a significance to a p-value is used, e.g., p-value = 84% and 98% correspond to 1 σ and 2 σ respectively.
9.
upper bounds are just limited by the number of trials
10.
the bandpass ranges in units of Hz are: (433, 789), (152, 733), (80, 308), (185, 794), (52, 251), (132, 351)
11.
( i j ) means echoes i through j were used
12.
the whole time range used was shifted 10 seconds later
13.
the explicit sets used: ( N E , t w / M ) = (15,40),(10,80),(15,80),(3-15,40),(5-17,40),(3-15,80)
14.
We here assume that the energy fraction of ringdown phase is ϵ rd = 2.4 × 10 6 although it highly depends on the detail of nonlinear gravitational collapse.
Figure 1. The Penrose diagram describing an evaporating BH.
Figure 1. The Penrose diagram describing an evaporating BH.
Universe 06 00043 g001
Figure 2. Near-horizon matter field of fuzzball solution from paper [95]. Shown in the figure, P 1 = P 2 . Asymptotically P 1 = P 3 .
Figure 2. Near-horizon matter field of fuzzball solution from paper [95]. Shown in the figure, P 1 = P 2 . Asymptotically P 1 = P 3 .
Universe 06 00043 g002
Figure 3. Near-horizon relationship between density and pressure of the fuzzball solution from paper [95]. The dashed line is the asymptotic behavior, and the solid line is the real behavior of P 1 P 3 ρ .
Figure 3. Near-horizon relationship between density and pressure of the fuzzball solution from paper [95]. The dashed line is the asymptotic behavior, and the solid line is the real behavior of P 1 P 3 ρ .
Universe 06 00043 g003
Figure 4. The geometry of Schwarzschild mock fuzzball. Red curve shows that the metric is modified around the horizon, and the shaded area is removed from the spacetime.
Figure 4. The geometry of Schwarzschild mock fuzzball. Red curve shows that the metric is modified around the horizon, and the shaded area is removed from the spacetime.
Universe 06 00043 g004
Figure 5. Anisotropic behavior of pressure of Schwarzschild mock fuzzball versus proper distance from horizon. The absolute value of two pressures are equal at r = 3 M 1 + a . Radial pressure (dashed line) is much larger than tangential (solid line) pressure near horizon (proper length = 0), while the trend reverses far away. The figure shows that when approaching horizon, both pressures reach their maxima, hence have no singularity at r = 2M.
Figure 5. Anisotropic behavior of pressure of Schwarzschild mock fuzzball versus proper distance from horizon. The absolute value of two pressures are equal at r = 3 M 1 + a . Radial pressure (dashed line) is much larger than tangential (solid line) pressure near horizon (proper length = 0), while the trend reverses far away. The figure shows that when approaching horizon, both pressures reach their maxima, hence have no singularity at r = 2M.
Universe 06 00043 g005
Figure 6. The energy density of mock fuzzball, compared with that of Hawking radiation as seen by a radially infalling observer along a geodesic. Here the orange curve is positive and the black curve is negative. The energy density of mock fuzzball is larger than the Hawking radiation near horizon and drops faster than the Hawking radiation away from the horizon.
Figure 6. The energy density of mock fuzzball, compared with that of Hawking radiation as seen by a radially infalling observer along a geodesic. Here the orange curve is positive and the black curve is negative. The energy density of mock fuzzball is larger than the Hawking radiation near horizon and drops faster than the Hawking radiation away from the horizon.
Universe 06 00043 g006
Figure 7. The energy flux seen by an observer falling through fuzzball geometry, compared to that of Hawking radiation. Here the orange curve is positive, while the black curves are negative. The flux of mock fuzzball is larger than the Hawking radiation near horizon and drops faster than the Hawking radiation away from the horizon.
Figure 7. The energy flux seen by an observer falling through fuzzball geometry, compared to that of Hawking radiation. Here the orange curve is positive, while the black curves are negative. The flux of mock fuzzball is larger than the Hawking radiation near horizon and drops faster than the Hawking radiation away from the horizon.
Universe 06 00043 g007
Figure 8. The CD potentials with ( i , j ) = ( 1 , + 1 ) and ( + 1 , 1 ) for a ¯ = 0.8 and = m = 2 . The potential V + 1 , 1 ( V 1 , + 1 ) is singular for r g ω = 0.5 ( 1.7 ) in this case.
Figure 8. The CD potentials with ( i , j ) = ( 1 , + 1 ) and ( + 1 , 1 ) for a ¯ = 0.8 and = m = 2 . The potential V + 1 , 1 ( V 1 , + 1 ) is singular for r g ω = 0.5 ( 1.7 ) in this case.
Universe 06 00043 g008
Figure 9. GW echoes following a BBH merger from a cavity of membrane/firewall-angular momentum barrier [18].
Figure 9. GW echoes following a BBH merger from a cavity of membrane/firewall-angular momentum barrier [18].
Universe 06 00043 g009
Figure 10. The amplification factor for a ¯ = 0.8 and = m = 2 . In the right panel, low frequency region (solid) is calculated with the potential V 1 , + 1 and the higher frequency region (dashed) is calculated from V + 1 , 1 in the CD equation.
Figure 10. The amplification factor for a ¯ = 0.8 and = m = 2 . In the right panel, low frequency region (solid) is calculated with the potential V 1 , + 1 and the higher frequency region (dashed) is calculated from V + 1 , 1 in the CD equation.
Universe 06 00043 g010
Figure 11. Spectra of echo + ringdown with a ¯ = 0.4 , = m = 2 , and M = 2.7 M , ϵ rd = 0.01 , D o = 40 Mpc, and θ = 90 . The left panel shows the spectrum in the constant reflectivity model with | R | = 0.5 and the right panel shows the spectrum in the Boltzmann reflectivity model with γ ˜ = 1 .
Figure 11. Spectra of echo + ringdown with a ¯ = 0.4 , = m = 2 , and M = 2.7 M , ϵ rd = 0.01 , D o = 40 Mpc, and θ = 90 . The left panel shows the spectrum in the constant reflectivity model with | R | = 0.5 and the right panel shows the spectrum in the Boltzmann reflectivity model with γ ˜ = 1 .
Universe 06 00043 g011
Figure 12. QNMs with | R | = 0.7 (red filled circles) and 0.9 (blue filled circles), r 0 = 40 M , a ¯ = 0.9 , and = m = 2 .
Figure 12. QNMs with | R | = 0.7 (red filled circles) and 0.9 (blue filled circles), r 0 = 40 M , a ¯ = 0.9 , and = m = 2 .
Universe 06 00043 g012
Figure 13. Mode function obtained by solving (128) with the boundary condition of ψ ω ˜ = constant . for r .
Figure 13. Mode function obtained by solving (128) with the boundary condition of ψ ω ˜ = constant . for r .
Universe 06 00043 g013
Figure 14. The energy reflection rate | R | 2 in the Boltzmann reflectivity model.
Figure 14. The energy reflection rate | R | 2 in the Boltzmann reflectivity model.
Universe 06 00043 g014
Figure 15. The time domain function with M = 2.7 M , a ¯ = 0.7 , D o = 40 Mpc, ϵ rd = 0.04 , θ = 90 , and = m = 2 in the Boltzmann reflectivity model with γ ˜ = 1 and T H / T QH = 0.6 .
Figure 15. The time domain function with M = 2.7 M , a ¯ = 0.7 , D o = 40 Mpc, ϵ rd = 0.04 , θ = 90 , and = m = 2 in the Boltzmann reflectivity model with γ ˜ = 1 and T H / T QH = 0.6 .
Universe 06 00043 g015
Figure 16. GW echoes following a collapse of binary neutron star merger event from a cavity of membrane-angular momentum barrier [20].
Figure 16. GW echoes following a collapse of binary neutron star merger event from a cavity of membrane-angular momentum barrier [20].
Universe 06 00043 g016
Figure 17. Maximized SNR2 around the expected time of merger echoes Equation (144), for the combined (top) and GW150914 (bottom) events. The significance and p-values of the peaks within the gray rectangle are specified in this plot [18].
Figure 17. Maximized SNR2 around the expected time of merger echoes Equation (144), for the combined (top) and GW150914 (bottom) events. The significance and p-values of the peaks within the gray rectangle are specified in this plot [18].
Universe 06 00043 g017
Figure 18. Same as Figure 17, but over an extended range of x = t echo t merger Δ t echo . The SNR peaks at the predicted value of 1 0.01 < x < 1 + 0.01 within gray rectangle have false detection probability of 0.11 (0.011) and significance of 1.6 σ (2.5 σ ), for GW150914 (combined events) [18] (See also [27]).
Figure 18. Same as Figure 17, but over an extended range of x = t echo t merger Δ t echo . The SNR peaks at the predicted value of 1 0.01 < x < 1 + 0.01 within gray rectangle have false detection probability of 0.11 (0.011) and significance of 1.6 σ (2.5 σ ), for GW150914 (combined events) [18] (See also [27]).
Universe 06 00043 g018
Figure 19. Average number of background peaks higher than a particular SNR-value within a time-interval 2 % × Δ t ¯ echo (gray rectangle in Figure 17 and Figure 18) for combined (left) and GW150914 (right) events [18]. The red dots show the observed SNR peak at t echo = 1.0054 Δ t echo (Figure 17 and Figure 18). The correspondence between SNR values and their significance is indicated in horizontal bar.
Figure 19. Average number of background peaks higher than a particular SNR-value within a time-interval 2 % × Δ t ¯ echo (gray rectangle in Figure 17 and Figure 18) for combined (left) and GW150914 (right) events [18]. The red dots show the observed SNR peak at t echo = 1.0054 Δ t echo (Figure 17 and Figure 18). The correspondence between SNR values and their significance is indicated in horizontal bar.
Universe 06 00043 g019
Figure 20. Original GW template of GW150914 [18], along with best fit echoes template (147).
Figure 20. Original GW template of GW150914 [18], along with best fit echoes template (147).
Universe 06 00043 g020
Figure 21. SNR 2 with respect to T = ( t echo t merger ) / Δ t echo . Solid, dashed, and dotted lines correspond to ρ 2 for combined (Hanford and Livingston), Hanford, and Livingston, respectively, for the best fit parameters of the event GW150914 [19].
Figure 21. SNR 2 with respect to T = ( t echo t merger ) / Δ t echo . Solid, dashed, and dotted lines correspond to ρ 2 for combined (Hanford and Livingston), Hanford, and Livingston, respectively, for the best fit parameters of the event GW150914 [19].
Universe 06 00043 g021
Figure 22. Time-frequency representations of X ( t , f ) (Equation (154); combining all harmonics with frequencies n × f , with n N ) around the merger for the BNS gravitational-wave event GW170817, observed through cross-correlating the two LIGO detectors [20]. The possible resonance peaks of echoes found in this plot are marked with a green squares. The color scale shows the peak at f peak = 72 ( ± 0.5 ) Hz and t t merger 1.0 s, is the highest peak in this diagram, from before and after the BNS merger (see Figure 23, Figure 24 and Figure 25). A secondary peak at the same frequency but t t merger 32.9 s is also highlighted in this plot.
Figure 22. Time-frequency representations of X ( t , f ) (Equation (154); combining all harmonics with frequencies n × f , with n N ) around the merger for the BNS gravitational-wave event GW170817, observed through cross-correlating the two LIGO detectors [20]. The possible resonance peaks of echoes found in this plot are marked with a green squares. The color scale shows the peak at f peak = 72 ( ± 0.5 ) Hz and t t merger 1.0 s, is the highest peak in this diagram, from before and after the BNS merger (see Figure 23, Figure 24 and Figure 25). A secondary peak at the same frequency but t t merger 32.9 s is also highlighted in this plot.
Universe 06 00043 g022
Figure 23. A 3d representation of Figure 22 within echo search frequency window f = 63 92 Hz [20]. This plot shows the tentative detection of echoes at f peak = 72 ( ± 0.5 ) Hz and t t merger 1.0 s clearly stands above noise.
Figure 23. A 3d representation of Figure 22 within echo search frequency window f = 63 92 Hz [20]. This plot shows the tentative detection of echoes at f peak = 72 ( ± 0.5 ) Hz and t t merger 1.0 s clearly stands above noise.
Universe 06 00043 g023
Figure 24. Amplitude-frequency plot of X ( t , f ) (Equation (154); combining all harmonics with frequencies n × f , with n N ) for the first peak (red) at 1.0 sec after the merger for the BNS merger gravitational-wave event GW170817, observed by x-correlating the LIGO detectors [20]. The same amplitude-frequency plot for a random time in data (yellow), is also shown for comparison. Solid area between 63 Hz and 92 Hz was the search frequency prior range for Planckian echoes.
Figure 24. Amplitude-frequency plot of X ( t , f ) (Equation (154); combining all harmonics with frequencies n × f , with n N ) for the first peak (red) at 1.0 sec after the merger for the BNS merger gravitational-wave event GW170817, observed by x-correlating the LIGO detectors [20]. The same amplitude-frequency plot for a random time in data (yellow), is also shown for comparison. Solid area between 63 Hz and 92 Hz was the search frequency prior range for Planckian echoes.
Universe 06 00043 g024
Figure 25. Amplitude-time plot of first echo peak at 1.0 sec after the merger at frequency of 72 Hz [20]. After this detection Gill et al. [149] with independent Astrophysical considerations also determined that the remnant of GW170817 must have collapsed to a BH after t coll = 0.98 0.26 + 0.31 Section Error-bar (in blue) is the time of collapse considering this independent observation in [149] compared to the detected signal of echoes which is also as a consequence of BH collapse. The shaded region is 0–1 s prior range after the merger.
Figure 25. Amplitude-time plot of first echo peak at 1.0 sec after the merger at frequency of 72 Hz [20]. After this detection Gill et al. [149] with independent Astrophysical considerations also determined that the remnant of GW170817 must have collapsed to a BH after t coll = 0.98 0.26 + 0.31 Section Error-bar (in blue) is the time of collapse considering this independent observation in [149] compared to the detected signal of echoes which is also as a consequence of BH collapse. The shaded region is 0–1 s prior range after the merger.
Universe 06 00043 g025
Figure 26. Average number of background peaks higher than a particular -X(t,f) within a frequency-intervals of 63–92 Hz and time-intervals of 1 Section the observed X ( t peak , f peak ) peak at 1.0 s after the merger is marked by red square. The horizontal bar shows the relation between X ( t , f ) values and their significance [20].
Figure 26. Average number of background peaks higher than a particular -X(t,f) within a frequency-intervals of 63–92 Hz and time-intervals of 1 Section the observed X ( t peak , f peak ) peak at 1.0 s after the merger is marked by red square. The horizontal bar shows the relation between X ( t , f ) values and their significance [20].
Universe 06 00043 g026
Figure 27. Correlation vs. echo time delay Δ t echo vs. N E for GW170104 using method II [21]. Here N E is the number of frequency steps between spikes.
Figure 27. Correlation vs. echo time delay Δ t echo vs. N E for GW170104 using method II [21]. Here N E is the number of frequency steps between spikes.
Universe 06 00043 g027
Figure 28. Determination of η for the events in O1 and O2 [22].
Figure 28. Determination of η for the events in O1 and O2 [22].
Universe 06 00043 g028
Figure 29. This plot shows whether it is possible to recover the potential signals with a variety of amplitudes in [23]. Here it can be seen that amplitudes less than A=0.1 in ADA [18] are difficult to be identified in data, while amplitude twice this value would be clearly identifiable.
Figure 29. This plot shows whether it is possible to recover the potential signals with a variety of amplitudes in [23]. Here it can be seen that amplitudes less than A=0.1 in ADA [18] are difficult to be identified in data, while amplitude twice this value would be clearly identifiable.
Universe 06 00043 g029
Figure 30. This plot shows injected and recovered values for γ . The diagonal line is accurate recovery [23]. The preference for γ = 1 (dashed line) at lower injection amplitudes can be clearly seen.
Figure 30. This plot shows injected and recovered values for γ . The diagonal line is accurate recovery [23]. The preference for γ = 1 (dashed line) at lower injection amplitudes can be clearly seen.
Universe 06 00043 g030
Figure 31. In this plot echoes signals are injected into Gaussian noise and are analysed for different overall amplitudes A [23]. Injection amplitudes above 10 22 of the peak strain, can be accurately recovered. For the injections with lower amplitudes no matter what value they take, recovered amplitudes are likely to be around 10 22 . Horizontal lines show the amplitudes reported in [18]. The injections are made with γ = 0.8 .
Figure 31. In this plot echoes signals are injected into Gaussian noise and are analysed for different overall amplitudes A [23]. Injection amplitudes above 10 22 of the peak strain, can be accurately recovered. For the injections with lower amplitudes no matter what value they take, recovered amplitudes are likely to be around 10 22 . Horizontal lines show the amplitudes reported in [18]. The injections are made with γ = 0.8 .
Universe 06 00043 g031
Figure 32. Reconstruction of cWB for the event GW151012 via color-coded time-frequency maps [25]. The upper plot shows the squared coherent network SNR and the plot bellow shows the normalized residual noise energy. The residual plot is given after the reconstructed signal was subtracted from the data. In this plots a secondary cluster occurring 200 ms after the merger (consistent with echo times predicted and seen by ADA, Equation (144)). The dashed vertical lines denote coalescence time for GW151012 (the network has used the Livingston detector time as a reference).
Figure 32. Reconstruction of cWB for the event GW151012 via color-coded time-frequency maps [25]. The upper plot shows the squared coherent network SNR and the plot bellow shows the normalized residual noise energy. The residual plot is given after the reconstructed signal was subtracted from the data. In this plots a secondary cluster occurring 200 ms after the merger (consistent with echo times predicted and seen by ADA, Equation (144)). The dashed vertical lines denote coalescence time for GW151012 (the network has used the Livingston detector time as a reference).
Universe 06 00043 g032
Figure 33. Reconstruction of cWB for the event GW151226, as color-coded time-frequency maps [25]. The upper plot shows the squared coherent network SNR and the plot bellow shows the normalized residual noise energy. The residual plot is given after the reconstructed signal was subtracted from the data. In this plots a secondary cluster occurring 100 ms after the merger (consistent with echo times predicted and seen by ADA, Equation (144)). The dashed vertical lines denote coalescence time for GW151226 (the network has used the Livingston detector time as a reference).
Figure 33. Reconstruction of cWB for the event GW151226, as color-coded time-frequency maps [25]. The upper plot shows the squared coherent network SNR and the plot bellow shows the normalized residual noise energy. The residual plot is given after the reconstructed signal was subtracted from the data. In this plots a secondary cluster occurring 100 ms after the merger (consistent with echo times predicted and seen by ADA, Equation (144)). The dashed vertical lines denote coalescence time for GW151226 (the network has used the Livingston detector time as a reference).
Universe 06 00043 g033
Figure 34. Maximum posteriori probability for time delay between Hanford (H) and Livingston (L) in line-of-sight frame for the main event GW151012 (blue contour) and the secondary signal (green contour) [25].
Figure 34. Maximum posteriori probability for time delay between Hanford (H) and Livingston (L) in line-of-sight frame for the main event GW151012 (blue contour) and the secondary signal (green contour) [25].
Universe 06 00043 g034
Figure 35. This plot shows how residual signal because of wavelength degeneracy can cause 7.7 m·s = 1 / ( 130 Hz ) shift in maximum poster probability for time delay. The upper plot is residual plot in Figure 32 and the plot below is the poster probability in Figure 34 for the event GW151012 [25].
Figure 35. This plot shows how residual signal because of wavelength degeneracy can cause 7.7 m·s = 1 / ( 130 Hz ) shift in maximum poster probability for time delay. The upper plot is residual plot in Figure 32 and the plot below is the poster probability in Figure 34 for the event GW151012 [25].
Universe 06 00043 g035
Figure 36. Plot of mass ratio dependence of p-values in [25]. Vertical lines are error bars for 50% credible region and central points are best value of mass ratio obtained from posteriors distribution. Because of relation of p-value to error function erf(SNR) we took roughly SNR log ( p value ) as horizontal axis [144].
Figure 36. Plot of mass ratio dependence of p-values in [25]. Vertical lines are error bars for 50% credible region and central points are best value of mass ratio obtained from posteriors distribution. Because of relation of p-value to error function erf(SNR) we took roughly SNR log ( p value ) as horizontal axis [144].
Universe 06 00043 g036
Figure 37. Histogram of slopes for uniform random choices of 0 < ( log ( p value ) ) 0.5 < 2.5 . We see that only 1.3% of these random realizations, the slope can exceed the observed value [144].
Figure 37. Histogram of slopes for uniform random choices of 0 < ( log ( p value ) ) 0.5 < 2.5 . We see that only 1.3% of these random realizations, the slope can exceed the observed value [144].
Universe 06 00043 g037
Figure 38. p-value distribution for combined events of different stretches of data within 1 min of the main events. Surprisingly, the blue line which is closest to the main event, and has used to define p-value in [18] (Figure 19), happens to give the smallest p-value. The shaded region represents the Poisson error range for blue histogram. This shows that the variation in p-values is clearly much larger. This behaviour is interpreted as non-gaussianity and/or non-stationarity of the LIGO noise. In this plot the y-axis on the left (right) shows p-value (number of higher peaks) within the mentioned range of data. In each histogram the total number of “peaks” is (38 − 9)/0.02 = 1450.
Figure 38. p-value distribution for combined events of different stretches of data within 1 min of the main events. Surprisingly, the blue line which is closest to the main event, and has used to define p-value in [18] (Figure 19), happens to give the smallest p-value. The shaded region represents the Poisson error range for blue histogram. This shows that the variation in p-values is clearly much larger. This behaviour is interpreted as non-gaussianity and/or non-stationarity of the LIGO noise. In this plot the y-axis on the left (right) shows p-value (number of higher peaks) within the mentioned range of data. In each histogram the total number of “peaks” is (38 − 9)/0.02 = 1450.
Universe 06 00043 g038
Figure 39. Spectra of ringdown and echo phases with the reflectivity of | R | = 0.99 , 0.9 , 0.6 , and 0.3 . We set D o = 40 Mpc, a ¯ = 0.1 , = m = 2 , M = 4 M , θ = 20 , and ϵ rd = 0.4 % .
Figure 39. Spectra of ringdown and echo phases with the reflectivity of | R | = 0.99 , 0.9 , 0.6 , and 0.3 . We set D o = 40 Mpc, a ¯ = 0.1 , = m = 2 , M = 4 M , θ = 20 , and ϵ rd = 0.4 % .
Universe 06 00043 g039
Figure 40. Spectra of ringdown and echo phases in the Boltzmann reflectivity model with a ¯ = 0.1 , ϵ rd = 2.4 × 10 6 , M = 2.4 M , θ = 90 , and D o = 1 Mpc. Here we also assume γ = 10 10 , T H / T QH = 1 (left) and T H / T QH = 1.37 × 10 6 (right).
Figure 40. Spectra of ringdown and echo phases in the Boltzmann reflectivity model with a ¯ = 0.1 , ϵ rd = 2.4 × 10 6 , M = 2.4 M , θ = 90 , and D o = 1 Mpc. Here we also assume γ = 10 10 , T H / T QH = 1 (left) and T H / T QH = 1.37 × 10 6 (right).
Universe 06 00043 g040
Figure 41. Spectra of GW echoes in the Boltzmann reflectivity model with a ¯ = 0.6 , = m = 2 , D o = 1 Gpc, and γ = 1 . The gray line shows the case of M = 3 × 10 9 M , ϵ rd = 0.02 , T H / T QH = 0.1 and the black line shows one for M = 8 × 10 9 M , ϵ rd = 0.04 , T H / T QH = 0.05 . We also plot the PSD for the IPTA (blue) and SKA (red).
Figure 41. Spectra of GW echoes in the Boltzmann reflectivity model with a ¯ = 0.6 , = m = 2 , D o = 1 Gpc, and γ = 1 . The gray line shows the case of M = 3 × 10 9 M , ϵ rd = 0.02 , T H / T QH = 0.1 and the black line shows one for M = 8 × 10 9 M , ϵ rd = 0.04 , T H / T QH = 0.05 . We also plot the PSD for the IPTA (blue) and SKA (red).
Universe 06 00043 g041
Table 1. Asymptotical behavior of metric.
Table 1. Asymptotical behavior of metric.
MetricFuzzballKerr BHSch. BH
g t t K + M K 3 1 r + O ( 1 r 2 ) 1 + r g 1 r + O ( 1 r 2 ) 1 + r g 1 r + O ( 1 r 2 )
g r r K + M K 1 r + O ( 1 r 2 ) 1 + r g 1 r + O ( 1 r 2 ) 1 + r g 1 r + O ( 1 r 2 )
g θ θ K r 2 + M K r + O ( 1 ) r 2 + O ( 1 ) r 2 + O ( 1 )
g φ φ K r 2 sin 2 θ + M K r + O ( 1 ) r 2 sin 2 θ + O ( 1 ) r 2 sin 2 θ + O ( 1 )
g t φ 0 + 4 c 1 c 5 J K 1 r + O ( 1 r 2 ) N r + L 1 r + O ( 1 r 2 ) 0
Table 2. Parameters Setting.
Table 2. Parameters Setting.
r 0 pm c 1 c 5 s 1 s 5
142211 1 2
Table 3. Table of positive results (p-value 0.05 ) by different groups ( The p-value for Nielsen et al. [24] is a rough estimate, based on the log-Bayes = 1.66 ).
Table 3. Table of positive results (p-value 0.05 ) by different groups ( The p-value for Nielsen et al. [24] is a rough estimate, based on the log-Bayes = 1.66 ).
AuthorsMethodDatap-Value
1Abedi, Dykaar, Afshordi (ADA) 2017 [18]ADA templateO11.1%
2Conklin, Holdom, Ren 2018 [21]spectral combO1+O20.2–0.8%
3Westerweck, et al. 2018 [23]ADA templateO12.0%
4Nielsen, et al. 2019 [24]ADA + BayesGW151012, GW1512262%
5Uchikata, et al. 2019 [19]ADA templateO15.5%
6Uchikata, et al. 2019 [19]ADA templateO23.9%
7Salemi, et al. 2019 [25]coherent WaveBurstGW151012, GW1512260.4%, 3%
8Abedi, Afshordi 2019 [20]spectral combBNS0.0016%
9Gill, Nathanail, Rezolla 2019 [149]Astro ModellingBNS EM t coll = t echo
Table 4. Table of failed searches and their possible caveat. “Infinite” prior refers to models that include unphysical non-decaying echoes. As discussed in the text, inclusion of such templates skews any statistical inference.
Table 4. Table of failed searches and their possible caveat. “Infinite” prior refers to models that include unphysical non-decaying echoes. As discussed in the text, inclusion of such templates skews any statistical inference.
AuthorsMethodDataPossible Caveat
1Westerweck, et al. 2018 [23]ADA templateO1“Infinite” prior
2Nielsen, et al. 2019 [24]ADA+BayesGW150914mass-ratio dependence,
inconsistent Max SNR
3Uchikata, et al. 2019 [19]ADA, hi-passO1,O2no low-frequencies
4Salemi, et al. 2019 [25]coherent WaveBurstO1,O2mass-ratio dependence,
only 1st echo
5Lo, et al. 2019 [145]ADA+BayesO1“Infinite” prior
6Tsang, et al. 2019 [146]BayesWaveO1+O2needs very loud echoes
(9 free parameters)
Table 5. Values of best fit echo parameters of the model Equation (147) of the highest SNR peak near the predicted Δ t echo (gray rectangle in Figure 17), and their significance [18].
Table 5. Values of best fit echo parameters of the model Equation (147) of the highest SNR peak near the predicted Δ t echo (gray rectangle in Figure 17), and their significance [18].
RangeGW150914Combined
( t echo t merger ) / Δ t echo (0.99, 1.01)1.00541.0054
γ (0.1, 0.9)0.890.9
t 0 / Δ t ¯ echo (−0.1,0)−0.084−0.1
Amplitude 0.09920.124
SNR max 4.216.96
p-value 0.11 0.011
significance 1.6 σ 2.5 σ
Table 6. Comparing the expected theoretical values of echo time delays Δ t echo ’s of each merger event (Equation (144)), to their best combined fit within the 1 σ credible region, and the contribution of each event to the combined SNR for the echoes (Equation (148)) [18].
Table 6. Comparing the expected theoretical values of echo time delays Δ t echo ’s of each merger event (Equation (144)), to their best combined fit within the 1 σ credible region, and the contribution of each event to the combined SNR for the echoes (Equation (148)) [18].
GW150914GW151012GW151226
Δ t echo , pred (sec)0.29250.17780.1013
± 0.00916± 0.02789± 0.01152
Δ t echo , best (sec)0.300680.190430.09758
| A best , I | 0.0910.340.33
SNR best , I 4.134.523.83
Table 7. p-values along with Poisson errors for O1 events [19].
Table 7. p-values along with Poisson errors for O1 events [19].
EventWesterweck et al. [23]Uchikata et al. [19]
GW150914 0.238 ± 0.043 0.157 ± 0.035
GW151012 0.063 ± 0.022 0.047 ± 0.019
GW151226 0.476 ± 0.061 0.598 ± 0.069
Total 0.032 ± 0.016 0.055 ± 0.021
Table 8. p-values for O2 events [19].
Table 8. p-values for O2 events [19].
EventUchikata et al. [19]
GW1701040.071
GW1706080.079
GW1707290.567
GW1708140.024
GW1708180.929
GW1708230.055
Total0.039
Table 9. The best-fit Δ t echo , p-value, bandpass and window parameters for the six signals [21]. Here N E is the number of frequency steps between spikes which has been chosen using injection and echoes model properties in [21].
Table 9. The best-fit Δ t echo , p-value, bandpass and window parameters for the six signals [21]. Here N E is the number of frequency steps between spikes which has been chosen using injection and echoes model properties in [21].
Event (Method)Best-Fit Δ t echo (s)p-ValueBandpass ( f min , f max ) Δ t echo Window Parameters for Average
GW151226 (I)0.0786<0.00139(34,62)10 N E = (1–29), (5-29), (9-29)11
GW151226 (II)0.07910.0076(12,58) N E = (260,270)
GW170104 (II)0.201<0.0018(16,62) N E = (100,125,150,175,200)
GW170608 (II)0.0756<0.004(14,60) N E = (140,200,260)
GW170814 (II)0.2310.04(12,58) N E = (170,190)12
GW170814 (III)0.2280.0077(30,80) N E = 10 17 , t w = 40 , 80 13
Table 10. Echoes best fit time delays and corresponding N E using method II for O1 and O2 events [22].
Table 10. Echoes best fit time delays and corresponding N E using method II for O1 and O2 events [22].
Event (Method)Best-Fit Δ t echo (s) N E Δ t echo / M
GW150914 (II)0.251200806
GW151012 (II)0.145160826
GW151226 (II)0.0791783270
GW170104 (II)0.201150831
GW170608 (II)0.0756200862
GW170729 (II)0.489180/1701240
GW170809 (II)0.235170845
GW170814 (II)0.231200878
GW170817 (II)0.00719250663
GW170818 (II)0.275140933
Table 11. Comparison of p-values obtained in [18] and using larger portion of data (4096 s of LOSC data) [23]. This data is divided into segments of 16 or 32 s length. Here different combinations of the events are considered, denoted as (GW150914, GW151012, GW151226, GW170104) ( 1 , 2 , 3 , 4 ) . Having the original priors, the Poisson errors (as suggested in [28]): for GW150914 our p-values are 0.199 ± 0.028 ( 0.238 ± 0.043 ), and for (1,2,3) our p-values are 0.02 ± 0.009 ( 0.032 ± 0.016 ). The Poisson errors for the full combination (1,2,3,4) with original priors, are 0.044 ± 0.013 ( 0.032 ± 0.016 ). The comparison of p-values using widened priors are also shown in this table.
Table 11. Comparison of p-values obtained in [18] and using larger portion of data (4096 s of LOSC data) [23]. This data is divided into segments of 16 or 32 s length. Here different combinations of the events are considered, denoted as (GW150914, GW151012, GW151226, GW170104) ( 1 , 2 , 3 , 4 ) . Having the original priors, the Poisson errors (as suggested in [28]): for GW150914 our p-values are 0.199 ± 0.028 ( 0.238 ± 0.043 ), and for (1,2,3) our p-values are 0.02 ± 0.009 ( 0.032 ± 0.016 ). The Poisson errors for the full combination (1,2,3,4) with original priors, are 0.044 ± 0.013 ( 0.032 ± 0.016 ). The comparison of p-values using widened priors are also shown in this table.
Event [18]Original 16 s (32 s)Widened Priors 16s (32s)
GW1509140.110.199 (0.238)0.705 (0.365)
GW151012-0.056 (0.063)0.124
GW151226-0.414 (0.476)0.837
GW170104-0.7250.757
(1,2)-0.0040.36
(1,3)-0.1590.801
(1,2,3)0.0110.020 (0.032)0.18 (0.144)
(1,3,4)-0.199 (0.072)0.9 (0.32)
(1,2,3,4)-0.044 (0.032)0.368 (0.112)
Table 12. Results of Bayes factor [24] using ADA model. Gaussian noise hypothesis is preferred for negative values of Log Bayes factor. Echoes hypothesis is preferred for positive values of Log Bayes factor. Log Bayes values of < 1 are “not worth more than a bare mention” [24]. However, we note that the Max SNR reported here by [24] for GW150914 is inconsistent with the one reported by the same group and ADA previously [18,23] (see Table 6).
Table 12. Results of Bayes factor [24] using ADA model. Gaussian noise hypothesis is preferred for negative values of Log Bayes factor. Echoes hypothesis is preferred for positive values of Log Bayes factor. Log Bayes values of < 1 are “not worth more than a bare mention” [24]. However, we note that the Max SNR reported here by [24] for GW150914 is inconsistent with the one reported by the same group and ADA previously [18,23] (see Table 6).
EventLog Bayes FactorMax SNR
GW150914−1.80562.86
GW1510121.24995.5741
GW1512260.41864.07
Table 13. p-values for post-coalescence deviations from GR obtained by cWB for the eleven GW events from GWTC-1 [25]. Post-coalescence SNR, SNR p c m i n R ; and estimated probability such that SNR p c m i n R produced by a noise fluctuation. Here P s u p and P i n f are refering to the probability of SNR p c i n f and SNR p c s u p , respectively.
Table 13. p-values for post-coalescence deviations from GR obtained by cWB for the eleven GW events from GWTC-1 [25]. Post-coalescence SNR, SNR p c m i n R ; and estimated probability such that SNR p c m i n R produced by a noise fluctuation. Here P s u p and P i n f are refering to the probability of SNR p c i n f and SNR p c s u p , respectively.
EventSourceSNRSNR pc minR SNR pc inf SNR pc sup p-Value pc P inf P sup
GW150914BBH25.2 5.72 5.64 6.92 0.94 ± 0.02 0.71 0.95
GW151012BBH10.5 6.60 6.26 6.54 0.0037 ± 0.0014 0.0042 0.0068
GW151226BBH11.9 4.40 4.36 4.41 0.025 ± 0.005 0.02 0.03
GW170104BBH13.0 5.29 3.95 5.30 0.07 ± 0.01 0.07 0.31
GW170608BBH14.1 1.69 1.64 1.75 0.51 ± 0.02 0.49 0.54
GW170729BBH10.2 4.81 3.43 4.86 0.09 ± 0.01 0.08 0.35
GW170809BBH11.9 3.89 3.88 4.71 0.28 ± 0.01 0.11 0.28
GW170814BBH17.2 5.98 5.94 6.02 0.10 ± 0.01 0.09 0.11
GW170817BNS29.3 0.21 0.21 0.21 0.55 ± 0.01 0.55 0.56
GW170818BBH8.6 1.97 1.76 2.04 0.87 ± 0.02 0.86 0.91
GW170823BBH10.8 3.11 2.69 3.54 0.60 ± 0.02 0.44 0.74
Table 14. Events, p-values reported in [25], and average mass ratio of event. Highlighted rows are the most significant post-coalescence signals reported in [25].
Table 14. Events, p-values reported in [25], and average mass ratio of event. Highlighted rows are the most significant post-coalescence signals reported in [25].
Eventp-Value ± 2 σ Average Mass Ratio
GW150914 0.94 ± 0.02 0.86
GW151012 0.0037 ± 0.0014 0.58
GW151226 0.025 ± 0.005 0.56
GW170104 0.07 ± 0.01 0.65
GW170608 0.51 ± 0.02 0.68
GW170729 0.09 ± 0.01 0.68
GW170809 0.28 ± 0.01 0.68
GW170814 0.10 ± 0.01 0.82
GW170818 0.87 ± 0.02 0.75
GW170823 0.60 ± 0.02 0.74
Table 15. Prior range proposed by Lo et al. [145] of the echo parameters of ADA model [18].
Table 15. Prior range proposed by Lo et al. [145] of the echo parameters of ADA model [18].
ParameterPrior Range
A[0.0, 1.0]
γ [0.0, 1.0]
t 0 (s)[−0.1,0.01]
t echo (s)[0.05, 0.5]
Δ t echo (s)[0.05, 0.5]
Table 16. Detection statistic ln B vs. its corresponding statistical significances shown for both Gaussian and O1 backgrounds [145].
Table 16. Detection statistic ln B vs. its corresponding statistical significances shown for both Gaussian and O1 backgrounds [145].
Statistical SignificanceDetection Statistic (Gaussian Noise)Detection Statistic (O1 Noise)
1 σ −0.90.1
2 σ −0.41.5
3 σ 1.14.0
4 σ 1.55.4
5 σ 1.95.7
Table 17. The detection statistic and its corresponding statistical significance and p-value for O1 events [145]. The ordering of events by their statistical significance is consistent with what reported by Nielsen et al. [24].
Table 17. The detection statistic and its corresponding statistical significance and p-value for O1 events [145]. The ordering of events by their statistical significance is consistent with what reported by Nielsen et al. [24].
EventDetection Statisticp-ValueStatistical Significance ( σ )
GW150914−1.30.806 < 1
GW1510120.40.08731.4
GW151226−0.20.254 < 1
Table 18. Obtained p-values for each event along with total p-value [19]. A hyphen means that 4096-s of data are not available. Data version.
Table 18. Obtained p-values for each event along with total p-value [19]. A hyphen means that 4096-s of data are not available. Data version.
EventC01C02
GW1509140.9920.984
GW1510120.6460.882
GW1512260.276-
GW1701040.7170.677
GW170608-0.488
GW170729-0.575
GW170814-0.472
GW170818-0.976
GW170823-0.315
Total0.9760.921
Table 19. p-value for each event and total p-value [19]. Result 1 is the case when the phase shift is fixed to π , and result 2 is the case when the total phase shift is also a parameter.
Table 19. p-value for each event and total p-value [19]. Result 1 is the case when the phase shift is fixed to π , and result 2 is the case when the total phase shift is also a parameter.
EventResult 1Result 2
GW1509140.6380.992
GW1510120.4170.646
GW1512260.9530.276
GW1701040.2130.717
Total0.5280.976
Table 20. Log Bayes factors for signal versus noise and signal versus glitch, and the corresponding p-values, for events seen in two detectors of GWTC-1 [147]. The bottom row shows the combined p-values for all these events together.
Table 20. Log Bayes factors for signal versus noise and signal versus glitch, and the corresponding p-values, for events seen in two detectors of GWTC-1 [147]. The bottom row shows the combined p-values for all these events together.
EventLog B S / N p S / N Log B S / G p S / N
GW1509142.320.262.950.43
GW151012−0.590.700.350.88
GW151226−0.670.722.480.53
GW1701041.090.443.800.28
GW170608−0.900.750.900.82
GW1708236.110.035.290.11
Combined-0.34-0.57
Table 21. Same as Table 20, while for the events that three detectors has been involved [147]. In the case of GW170817, in order to cover 1.0 sec after the merger for echoes found by Abedi and Afshordi [20], additional search as first echo being from this time has been set. For this particular event, latter prior choice has been taken for combined p-values.
Table 21. Same as Table 20, while for the events that three detectors has been involved [147]. In the case of GW170817, in order to cover 1.0 sec after the merger for echoes found by Abedi and Afshordi [20], additional search as first echo being from this time has been set. For this particular event, latter prior choice has been taken for combined p-values.
EventLog B S / N p S / N Log B S / G p S / N
GW1707294.240.675.640.62
GW1708099.050.3112.690.09
GW1708148.750.338.540.34
GW17081711.050.1910.300.20
GW170817 + 1 s 6.190.529.390.27
GW17081810.390.239.360.27
Combined-0.47-0.22

Share and Cite

MDPI and ACS Style

Abedi, J.; Afshordi, N.; Oshita, N.; Wang, Q. Quantum Black Holes in the Sky. Universe 2020, 6, 43. https://doi.org/10.3390/universe6030043

AMA Style

Abedi J, Afshordi N, Oshita N, Wang Q. Quantum Black Holes in the Sky. Universe. 2020; 6(3):43. https://doi.org/10.3390/universe6030043

Chicago/Turabian Style

Abedi, Jahed, Niayesh Afshordi, Naritaka Oshita, and Qingwen Wang. 2020. "Quantum Black Holes in the Sky" Universe 6, no. 3: 43. https://doi.org/10.3390/universe6030043

APA Style

Abedi, J., Afshordi, N., Oshita, N., & Wang, Q. (2020). Quantum Black Holes in the Sky. Universe, 6(3), 43. https://doi.org/10.3390/universe6030043

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop