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 4 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".

The first direct observation of a binary neutron star (BNS) merger was a watershed moment in multi-messenger astronomy. However, gravitational waves from GW170817 have only been observed prior to the BNS merger, but electromagnetic observations all follow the merger event. While post-merger gravitational wave signal in general relativity is too faint (given current detector sensitivities), here we present the first tentative detection of post-merger gravitational wave "echoes" from a highly spinning "black hole" remnant. The echoes may be expected in different models of quantum black holes that replace event horizons by exotic Planck-scale structure and tentative evidence for them has been found in binary black hole merger events. The fact that the echo frequency is suppressed by $\log M$ (in Planck units) puts it squarely in the LIGO sensitivity window, allowing us to build an optimal model-agnostic search strategy via cross-correlating the two detectors in frequency/time. We find a tentative detection of echoes at $f_{\rm echo} \simeq 72$ Hz, around 1.0 sec after the BNS merger, consistent with a 2.6-2.7 $M_\odot$ "black hole" remnant with dimensionless spin $0.84-0.87$. Accounting for all the "look-elsewhere" effects, we find a significance of $4.2 \sigma$, or a false alarm probability of $1.6\times 10^{-5}$, 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. If confirmed, this finding will have significant consequences for both physics of quantum black holes and astrophysics of binary neutron star mergers [Note added: This result is independently confirmed by arXiv:1901.04138, who use the electromagnetic observations to infer $t_{\rm coll}=0.98^{+0.31}_{-0.26}$ sec for black hole formation].

Gravitational wave echoes provide our most direct and surprising observational window into quantum nature of black holes. Three years ago, the first search for echoes from Planck-scale modifications of general relativity near black hole event horizons led to tentative evidence at false detection probability of 1\% arXiv:1612.00266 . The study introduced a naive phenomenological model and used the public data release by the Advanced LIGO gravitational wave observatory for the first observing run O1 (GW150914, GW151226, and LVT151012, now GW151012). Here, we provide a status update on various observational searches for echoes by independent groups, and argue that they can all be consistent if echoes are most prominent at lower frequencies and/or in binary mergers of more extreme mass ratio. We also point out that the only reported "detection" of echoes (with $>4\sigma$ confidence) at 1.0 second after the binary neutron star merger GW170817 arXiv:1803.10454 is coincident with the formation time of the black hole inferred from electromagnetic observations.

In a recent publication (1612.00266), we demonstrated that the events in the first observing run of the Advanced LIGO gravitational wave observatory (aLIGO O1) showed tentative evidence for repeating "echoes from the abyss" caused by Planck-scale structure near black hole horizons. By considering a phenomenological echo model, we showed that the pure noise hypothesis is disfavored with a p-value of 1%, i.e. higher amplitude for echoes than those in aLIGO O1 events are only recovered in 1% of random noise realizations. A recent preprint by Westerweck, et al. (1712.09966) provides a careful re-evaluation of our analysis which claims "a reduced statistical significance ... entirely consistent with noise". It is a mystery to us why the authors make such a statement, while they also find a p-value of 2 $\pm$ 1% (given the Poisson error in their estimate) for the same model and dataset. This is p-erfectly consistent with our results, which would be commonly considered as disfavoring the null hypothesis, or "moderate to significant" evidence for "echoes". Westerweck, et al. also point to diversity of the observed echo properties as evidence for statistical fluke, but such a diversity is neither unique nor surprising for complex physical phenomena in nature.

In classical General Relativity (GR), an observer falling into an astrophysical black hole is not expected to experience anything dramatic as she crosses the event horizon. However, tentative resolutions to problems in quantum gravity, such as the cosmological constant problem, or the black hole information paradox, invoke significant departures from classicality in the vicinity of the horizon. It was recently pointed out that such near-horizon structures can lead to late-time echoes in the black hole merger gravitational wave signals that are otherwise indistinguishable from GR. We search for observational signatures of these echoes in the gravitational wave data released by advanced Laser Interferometer Gravitational-Wave Observatory (LIGO), following the three black hole merger events GW150914, GW151226, and LVT151012. In particular, we look for repeating damped echoes with time-delays of $8 M \log M$ (+spin corrections, in Planck units), corresponding to Planck-scale departures from GR near their respective horizons. Accounting for the "look elsewhere" effect due to uncertainty in the echo template, we find tentative evidence for Planck-scale structure near black hole horizons at false detection probability of $1\%$ (corresponding to $2.5\sigma$ significance level). Future observations from interferometric detectors at higher sensitivity, along with more physical echo templates, will be able to confirm (or rule out) this finding, providing possible empirical evidence for alternatives to classical black holes, such as in ${\it firewall}$ or ${\it fuzzball}$ paradigms.

In a recent paper (arXiv:1612.00266), we reported the results of the first search for echoes from Planck-scale modifications of general relativity near black hole event horizons using the public data release by the Advanced LIGO gravitational wave observatory. While we found tentative evidence (at $\simeq 3 \sigma$ level) for the presence of these echoes, our statistical methodology was challenged by Ashton, et al. (arXiv:1612.05625), just in time for the holidays! In this short note, we briefly address these criticisms, arguing that they either do not affect our conclusion or change its significance by $\lesssim 0.3\sigma$. The real test will be whether our finding can be reproduced by independent groups using independent methodologies (and ultimately more data).