Using a Bayesian approach, we combine measurements of neutron star macroscopic observables obtained by astrophysical and gravitational observations, to derive joint constraints on the equation of state (EoS) of matter at supranuclear density. In our analysis we use two sets of data: (i) the masses and tidal deformabilities measured in the binary neutron star event GW170817, detected by LIGO and Virgo; (ii) the masses and stellar radii measured from observations of nuclear bursts in accreting low-mass X-ray binaries. Using two different parametrizations of the equation of state, we compute the posteriorprobability distributions of the EoS parameters, using which we infer the posterior distribution for the radius and the mass of the two neutron stars of GW170817. The constraints we set on the radii are tighter than previous bounds.

Supermassive binaries detectable by the future space gravitational-wave interferometer LISA might allow to distinguish black holes from ultracompact horizonless objects, even when the latter are motivated by quantum-gravity considerations. We show that a measurement of very small tidal Love numbers at the level of $10\%$ accuracy (as achievable with "golden binaries") may also allow to distinguish between different models of these exotic compact objects, even when taking into account an intrinsic uncertainty in the object radius putatively due to quantum mechanics. We argue that there is no conceptual obstacle in performing these measurements, the main challenge remains the detectability of small tidal effects and an accurate waveform modelling. Our analysis uses only coordinate-independent quantities related to the proper radial distance and the total mass of the object.

Observations of gravitational wave (GW) signals produced by coalescing binary neutron stars (NS), like the GW event GW170817, can be exploited to constrain the equation of state (EoS) of matter in the stars' inner core. The information on the internal structure and composition of the stars is encoded in their tidal Love numbers, which leave an imprint in the waveform of the GW signal emitted from the binary during the late inspiral phase. We extended the theory of tidal deformations of compact objects by computing the spin-tidal corrections that affect the dynamics and the GW emission of a binary system at the leading post-Newtonian (PN) order and to linear order in the spin. These corrections are divided into two classes: terms due to the coupling between the standard tidal Love numbers and the spins of the objects, and terms depending on the rotational tidal Love numbers. Both enter the GW phase at 6.5PN order. We analysed the impact of the spin-tidal couplings by estimating the parameter bias induced on GW170817-like events, assuming second- and third-generation ground based interferometers. If relatively high spinning ($\chi \gtrsim 0.1$) NS binaries exist in nature, these effects might be observed by the next generation of detectors. Lastly, we proved the feasibility of solving the so-called inverse stellar problem using GW detections, i.e., reconstructing the EoS from the measurement of NS masses and tidal Love numbers. Our results show that few observations of coalescing binary NS by a network of advanced detectors would allow us to put interesting constraints on the phenomenological parameters of a piecewise polytropic representation of the EoS, and to perform a model selection among the realistic EoS proposed in the literature.

The grand challenges of contemporary fundamental physics---dark matter, dark energy, vacuum energy, inflation and early universe cosmology, singularities and the hierarchy problem---all involve gravity as a key component. And of all gravitational phenomena, black holes stand out in their elegant simplicity, while harbouring some of the most remarkable predictions of General Relativity: event horizons, singularities and ergoregions. The hitherto invisible landscape of the gravitational Universe is being unveiled before our eyes: the historical direct detection of gravitational waves by the LIGO-Virgo collaboration marks the dawn of a new era of scientific exploration. Gravitational-wave astronomy will allow us to test models of black hole formation, growth and evolution, as well as models of gravitational-wave generation and propagation. It will provide evidence for event horizons and ergoregions, test the theory of General Relativity itself, and may reveal the existence of new fundamental fields. The synthesis of these results has the potential to radically reshape our understanding of the cosmos and of the laws of Nature. The purpose of this work is to present a concise, yet comprehensive overview of the state of the art in the relevant fields of research, summarize important open problems, and lay out a roadmap for future progress.

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.

We compute the spin-tidal couplings that affect the dynamics of two orbiting bodies at the leading order in the post-Newtonian (PN) framework and to linear order in the spin. These corrections belong to two classes: (i) terms arising from the coupling between the ordinary tidal terms and the point-particle terms, which depend on the standard tidal Love numbers of order $l$ and affect the gravitational-wave (GW) phase at $(2l+5/2)$PN order and (ii) terms depending on the rotational tidal Love numbers, recently introduced in previous work, that affect the GW phase at $(2l+1/2+\delta_{2l})$PN order. For circular orbits and spins orthogonal to the orbital plane, all leading-order spin-tidal terms enter the GW phase at $1.5$PN order relative to the standard, quadrupolar, tidal deformability term (and, thus, before the standard octupolar tidal deformability terms). We present the GW phase that includes all tidal terms up to $6.5$PN order and to linear order in the spin. We comment on a conceptual issue related to the inclusion of the rotational tidal Love numbers in a Lagrangian formulation and on the relevance of spin-tidal couplings for parameter estimation in coalescing neutron-star binaries and for tests of gravity.

The LIGO/Virgo collaboration has recently announced the direct detection of gravitational waves emitted in the coalescence of a neutron star binary. This discovery allows, for the first time, to set new constraints on the behavior of matter at supranuclear density, complementary with those coming from astrophysical observations in the electromagnetic band. In this paper we demonstrate the feasibility of using gravitational signals to solve the relativistic inverse stellar problem, i.e. to reconstruct the parameters of the equation of state (EoS) from measurements of the stellar mass and tidal Love number. We perform Bayesian inference of mock data, based on different models of the star internal composition, modeled through piecewise polytropes. Our analysis shows that the detection of a small number of sources by a network of advanced interferometers would allow to put accurate bounds on the EoS parameters, and to perform a model selection among the realistic equations of state proposed in the literature.

Several quantum-gravity models of compact objects predict microscopic or even Planckian corrections at the horizon scale. We explore the possibility of measuring two model-independent, smoking-gun effects of these corrections in the gravitational waveform of a compact binary, namely the absence of tidal heating and the presence of tidal deformability. For events detectable by the future space-based interferometer LISA, we show that the effect of tidal heating dominates and allows one to constrain putative corrections down to the Planck scale. The measurement of the tidal Love numbers with LISA is more challenging but, in optimistic scenarios, it allows to constrain the compactness of a supermassive exotic compact object down to the Planck scale. Our analysis suggests that highly-spinning, supermassive binaries at 1-20 Gpc provide unparalleled tests of quantum-gravity effects at the horizon scale.