The black hole uniqueness and the no-hair theorems imply that the quasinormal spectrum of any astrophysical black hole is determined solely by its mass and spin. The countably infinite number of quasinormal modes of a Kerr black hole are thus related to each other and any deviations from these relations provide a strong hint for physics beyond the general theory of relativity. To test the no-hair theorem using ringdown signals, it is necessary to detect at least two quasinormal modes. In particular, one can detect the fundamental mode along with a subdominant overtone or with another angular mode, depending on the mass ratio and the spins of the progenitor binary. Also in the light of the recent discovery of GW190412, studying how the mass ratio affects the prospect of black hole spectroscopy using overtones or angular modes is pertinent, and this is the major focus of our study. First, we provide ready-to-use fits for the amplitudes and phases of both the angular modes and overtones as a function of mass ratio $q\in[0,10]$. Using these fits we estimate the minimum signal-to-noise ratio for detectability, resolvability, and measurability of subdominant modes/tones. We find that performing black-hole spectroscopy with angular modes is preferable when the binary mass ratio is larger than $q\approx 1.2$ (provided that the source is not located at a particularly disfavoured inclination angle). For nonspinning, equal-mass binary black holes, the overtones seem to be the only viable option to perform a spectroscopy test of the no-hair theorem. However this would require a large ringdown signal-to-noise ratio ($\approx 100$ for a $5\%$ accuracy test with two overtones) and the inclusion of more than one overtone to reduce modelling errors, making black-hole spectroscopy with overtones impractical in the near future.

Within the framework of $f(R)=R+\alpha R^2$ gravity, we study realistic models of neutron stars, using equations of state compatible with the LIGO constraints. i.e. APR4, MPA1, SLy, and WW1. By numerically solving modified Tolman-Oppenheimer-Volkoff equations, we investigate the Mass--Radius relation in both metric and torsional $f(R)=R+\alpha R^2$ gravity models. In particular, we observe that torsion effects decrease the compactness and total mass of neutron star with respect to the General Relativity predictions, therefore mimicking the effects of a repulsive massive field. The opposite occurs in the metric theory, where mass and compactness increase with $\alpha$, thus inducing an excess of mass that overtakes the standard General Relativity limit. We also find that the sign of $\alpha$ must be reversed whether one considers the metric theory (positive) or torsion (negative) to avoid blowing up solutions. This could draw an easy test to either confirm or discard one or the other theory by determining the sign of parameter $\alpha$.

Validating the black-hole no-hair theorem with gravitational-wave observations of compact binary coalescences provides a compelling argument that the remnant object is indeed a black hole as described by the general theory of relativity. This requires performing a spectroscopic analysis of the post-merger signal and resolving the frequencies of either different angular modes or overtones (of the same angular mode). For a nearly-equal mass binary black-hole system, only the dominant angular mode ($l=m=2$) is sufficiently excited and the overtones are instrumental to perform this test. Here we investigate the robustness of modelling the post-merger signal of a binary black hole coalescence as a superposition of overtones. Further, we study the bias expected in the recovered frequencies as a function of the start time of a spectroscopic analysis and provide a computationally cheap procedure to choose it based on the interplay between the expected statistical error due to the detector noise and the systematic errors due to waveform modelling. Moreover, since the overtone frequencies are closely spaced, we find that resolving the overtones is particularly challenging and requires a loud ringdown signal. Rayleigh's resolvability criterion suggests that in an optimistic scenario a ringdown signal-to-noise ratio larger than $\sim 30$ (achievable possibly with LIGO at design sensitivity and routinely with future interferometers such as Einstein Telescope, Cosmic Explorer, and LISA) is necessary to resolve the overtone frequencies. We then conclude by discussing some conceptual issues associated with black-hole spectroscopy with overtones.

In this brief note, we clarify certain aspects related to the magnetic (i.e., odd parity or axial) tidal Love numbers of a star in general relativity. Magnetic tidal deformations of a compact star had been computed in 2009 independently by Damour and Nagar and by Binnington and Poisson. More recently, Landry and Poisson showed that the magnetic tidal Love numbers depend on the assumptions made on the fluid, in particular they are different (and of opposite sign) if the fluid is assumed to be in static equilibrium or if it is irrotational. We show that the zero-frequency limit of the Regge-Wheeler equation forces the fluid to be irrotational. For this reason, the results of Damour and Nagar are equivalent to those of Landry and Poisson for an irrotational fluid, and are expected to be the most appropriate to describe realistic configurations.

For a brief moment, a binary black hole (BBH) merger can be the most powerful astrophysical event in the visible universe. Here we present a model fit for this gravitational-wave peak luminosity of nonprecessing quasicircular BBH systems as a function of the masses and spins of the component black holes, based on numerical relativity (NR) simulations and the hierarchical fitting approach introduced by X. Jim\'enez-Forteza et al. [Phys. Rev. D 95, 064024 (2017), arXiv:1611.00332]. This fit improves over previous results in accuracy and parameter-space coverage and can be used to infer posterior distributions for the peak luminosity of future astrophysical signals like GW150914 and GW151226. The model is calibrated to the l<=6 modes of 378 nonprecessing NR simulations up to mass ratios of 18 and dimensionless spin magnitudes up to 0.995, and includes unequal-spin effects. We also constrain the fit to perturbative numerical results for large mass ratios. Studies of key contributions to the uncertainty in NR peak luminosities, such as (i) mode selection, (ii) finite resolution, (iii) finite extraction radius, and (iv) different methods for converting NR waveforms to luminosity, allow us to use NR simulations from four different codes as a homogeneous calibration set. This study of systematic fits to combined NR and large-mass-ratio data, including higher modes, also paves the way for improved inspiral-merger-ringdown waveform models.