We study novel solitonic solutions to Einstein-Klein-Gordon theory in the presence of a periodic scalar potential arising in models of axion-like particles. The potential depends on two parameters: the mass of the scalar field $m_a$ and the decay constant $f_a$; the standard case of the QCD axion is recovered when $m_a\propto1/f_a$. When $f_a\to\infty$ the solutions reduce to the standard case of "mini" boson stars supported by a massive free scalar field. As the energy scale $f_a$ of the scalar self-interactions decreases we unveil several novel features of the solution: new stability branches emerge at high density, giving rise to very compact, radially stable, boson stars. Some of the most compact configurations acquire a photon sphere. When $f_a$ is at the GUT scale, a boson star made of QCD axions can have a mass up to ten solar masses and would be more compact than a neutron star. Gravitational-wave searches for these exotic compact objects might provide indirect evidence for ultralight axion-like particles in a region not excluded by the black-hole superradiant instability.

In viable models of minicharged dark matter, astrophysical black holes might be charged under a hidden $U(1)$ symmetry and are formally described by the same Kerr-Newman solution of Einstein-Maxwell theory. These objects are unique probes of minicharged dark matter and dark photons. We show that the recent gravitational-wave detection of a binary black-hole coalescence by aLIGO provides various observational bounds on the black hole's charge, regardless of its nature. The pre-merger inspiral phase can be used to constrain the dipolar emission of (ordinary and dark) photons, whereas the detection of the quasinormal modes set an upper limit on the final black hole's charge. By using a toy model of a point charge plunging into a Reissner-Nordstrom black hole, we also show that in dynamical processes the (hidden) electromagnetic quasinormal modes of the final object are excited to considerable amplitude in the gravitational-wave spectrum only when the black hole is nearly extremal. The coalescence produces a burst of low-frequency dark photons which might provide a possible electromagnetic counterpart to black-hole mergers in these scenarios.

Linear perturbations of spherically symmetric spacetimes in general relativity are described by radial wave equations, with potentials that depend on the spin of the perturbing field. In previous work we studied the quasinormal mode spectrum of spacetimes for which the radial potentials are slightly modified from their general relativistic form, writing generic small modifications as a power-series expansion in the radial coordinate. We assumed that the perturbations in the quasinormal frequencies are linear in some perturbative parameter, and that there is no coupling between the perturbation equations. In general, matter fields and modifications to the gravitational field equations lead to coupled wave equations. Here we extend our previous analysis in two important ways: we study second-order corrections in the perturbative parameter, and we address the more complex (and realistic) case of coupled wave equations. We highlight the special nature of coupling-induced corrections when two of the wave equations have degenerate spectra, and we provide a ready-to-use recipe to compute quasinormal modes. We illustrate the power of our parametrization by applying it to various examples, including dynamical Chern-Simons gravity, Horndeski gravity and an effective field theory-inspired model.

Using numerical methods, we investigate the absorption properties of a family of nonsingular solutions {which arise in different metric-affine theories, such as quadratic and Born-Infeld gravity.} These solutions continuously interpolate between Schwarzschild black holes and naked solitons with wormhole topology. The resulting spectrum is characterized by a series of quasibound states excitations, associated with the existence of a stable photonsphere.

Scalar-tensor theories of gravity where a new scalar degree of freedom couples to the Gauss-Bonnet invariant can exhibit the phenomenon of spontaneous black hole scalarization. These theories admit both the classic black holes predicted by general relativity as well as novel hairy black hole solutions. The stability of the hairy black holes is strongly dependent on the precise form of the scalar-gravity coupling. A radial stability investigation revealed that all scalarized black hole solutions are unstable when the coupling between the scalar field and the Gauss-Bonnet invariant is quadratic in the scalar, whereas stable solutions exist for exponential couplings. Here we elucidate this behavior. We demonstrate that, while the quadratic term controls the onset of the tachyonic instability that gives rise to the black hole hair, the higher-order coupling terms control the nonlinearities that quench that instability, and hence also control the stability of the hairy black hole solutions.

It has recently been shown that nontrivial couplings between a scalar and the Gauss-Bonnet invariant can give rise to black hole spontaneous scalarization. Theories that exhibit this phenomenon are among the leading candidates for testing gravity with upcoming black hole observations. All models considered so far have focused on specific forms for the coupling, neglecting scalar self-interactions. In this work, we take the first steps towards placing this phenomenon on a more robust theoretical footing by considering the leading-order scalar self-interactions as well as the scalar-Gauss-Bonnet coupling. Our approach is consistent with the principles of effective field theory and yields the simplest and most natural model. We find that a mass term for the scalar alters the threshold for the onset of scalarization, and we study the mass range over which scalarized black hole solutions exist. We also demonstrate that the quartic self-coupling is sufficient to produce scalarized solutions that are stable against radial perturbations, without the need to resort to higher-order terms in the Gauss-Bonnet coupling function. Our model therefore represents a canonical model that can be studied further, with the ultimate aim of developing falsifiable tests of black hole scalarization.

Gravitational waves emitted by distorted black holes---such as those arising from the coalescence of two neutron stars or black holes---carry not only information about the corresponding spacetime but also about the underlying theory of gravity. Although general relativity remains the simplest, most elegant and viable theory of gravitation, there are generic and robust arguments indicating that it is not the ultimate description of the gravitational universe. Here we focus on a particularly appealing extension of general relativity, which corrects Einstein's theory through the addition of terms which are second order in curvature: the topological Gauss-Bonnet invariant coupled to a dilaton. We study gravitational-wave emission from black holes in this theory, and {\bf(i)} find strong evidence that black holes are linearly (mode) stable against both axial and polar perturbations; {\bf(ii)} discuss how the quasinormal modes of black holes can be excited during collisions involving black holes, and finally {\bf(iii)} show that future ringdown detections with large signal-to-noise ratio would improve current constraints on the coupling parameter of the theory.

Gravitational waves from binary coalescences provide one of the cleanest signatures of the nature of compact objects. It has been recently argued that the post-merger ringdown waveform of exotic ultracompact objects is initially identical to that of a black-hole, and that putative corrections at the horizon scale will appear as secondary pulses after the main burst of radiation. Here we extend this analysis in three important directions: (i) we show that this result applies to a large class of exotic compact objects with a photon sphere for generic orbits in the test-particle limit; (ii) we investigate the late-time ringdown in more detail, showing that it is universally characterized by a modulated and distorted train of "echoes" of the modes of vibration associated with the photon sphere; (iii) we study for the first time equal-mass, head-on collisions of two ultracompact boson stars and compare their gravitational-wave signal to that produced by a pair of black-holes. If the initial objects are compact enough as to mimic a binary black-hole collision up to the merger, the final object exceeds the maximum mass for boson stars and collapses to a black-hole. This suggests that - in some configurations - the coalescence of compact boson stars might be almost indistinguishable from that of black-holes. On the other hand, generic configurations display peculiar signatures that can be searched for in gravitational-wave data as smoking guns of exotic compact objects.

We provide a detailed analysis of how bosonic dark matter "condensates" interact with compact stars, extending significantly the results of a recent Letter. We focus on bosonic fields with mass $m_B$, such as axions, axion-like candidates and hidden photons. Self-gravitating bosonic fields generically form "breathing" configurations, where both the spacetime geometry and the field oscillate, and can interact and cluster at the center of stars. We construct stellar configurations formed by a perfect fluid and a bosonic condensate, and which may describe the late stages of dark-matter accretion onto stars, in dark matter-rich environments. These composite stars oscillate at a frequency which is a multiple of $f=2.5\times 10^{14}\,\left(m_{B}c^2/eV\right)\,{\rm Hz}$. Using perturbative analysis and Numerical Relativity techniques, we show that these stars are generically stable, and we provide criteria for instability. Our results also indicate that the growth of the dark matter core is halted close to the Chandrasekhar limit. We thus dispel a myth concerning dark matter accretion by stars: dark mater accretion does not necessarily lead to the destruction of the star, nor to collapse to a black hole. Finally, we argue that stars with long-lived bosonic cores may also develop in other theories with effective mass couplings, such as (massless) scalar-tensor theories.

Some models (such as the Skyrme model, a low-energy effective field theory for QCD) suggest that the high-density matter prevailing in neutron star interiors may be significantly anisotropic. Anisotropy is known to affect the bulk properties of nonrotating neutron stars in General Relativity. In this paper we study the effects of anisotropy on slowly rotating stars in General Relativity. We also consider one of the most popular extensions of Einstein's theory, namely scalar-tensor theories allowing for spontaneous scalarization (a phase transition similar to spontaneous magnetization in ferromagnetic materials). Anisotropy affects the moment of inertia of neutron stars (a quantity that could potentially be measured in binary pulsar systems) in both theories. We find that the effects of scalarization increase (decrease) when the tangential pressure is bigger (smaller) than the radial pressure, and we present a simple criterion to determine the onset of scalarization by linearizing the scalar-field equation. Our calculations suggest that binary pulsar observations may constrain the degree of anisotropy or even, more optimistically, provide evidence for anisotropy in neutron star cores.

Ultracompact objects are self-gravitating systems with a light ring. It was recently suggested that fluctuations in the background of these objects are extremely long-lived and might turn unstable at the nonlinear level, if the object is not endowed with a horizon. If correct, this result has important consequences: objects with a light ring are black holes. In other words, the nonlinear instability of ultracompact stars would provide a strong argument in favor of the "black hole hypothesis," once electromagnetic or gravitational-wave observations confirm the existence of light rings. Here we explore in some depth the mode structure of ultracompact stars, in particular constant-density stars and gravastars. We show that the existence of very long-lived modes -- localized near a second, stable null geodesic -- is a generic feature of gravitational perturbations of such configurations. Already at the linear level, such modes become unstable if the object rotates sufficiently fast to develop an ergoregion. Finally, we conjecture that the long-lived modes become unstable under fragmentation via a Dyson-Chandrasekhar-Fermi mechanism at the nonlinear level. Depending on the structure of the star, it is also possible that nonlinearities lead to the formation of small black holes close to the stable light ring. Our results suggest that the mere observation of a light ring is a strong evidence for the existence of black holes.

The nature and properties of dark matter (DM) are both outstanding issues in physics. Besides clustering in halos, the universal character of gravity implies that self-gravitating compact DM configurations might be spread throughout the universe. The astrophysical signature of these objects may be used to probe fundamental particle physics, or even to provide an alternative description of compact objects in active galactic nuclei. Here we discuss the most promising dissection tool of these configurations: the inspiral of a compact stellar-size object and consequent gravitational-wave emission. The inward motion of this "test probe" encodes unique information about the nature of the central, supermassive DM configuration. When the probe travels through some compact DM profile we show that, within a Newtonian approximation, the quasi-adiabatic evolution of the inspiral is mainly driven by DM accretion into the small compact object and by dynamical friction, rather than by gravitational-wave radiation-reaction. These effects circularize the orbits and leave a peculiar imprint on the gravitational waves emitted at late time. When accretion dominates, the frequency and the amplitude of the gravitational-wave signal produced during the latest stages of the inspiral are nearly constant. In the exterior region we study a relativistic model in which the inspiral is driven by the emission of gravitational and scalar waves. Resonances in the energy flux appear whenever the orbital frequency matches the mass of the DM particle and they correspond to the excitation of the central object's quasinormal frequencies. Unexpectedly, these resonances can lead to large dephasing with respect to standard inspiral templates, to such an extent as to prevent detection with matched filtering techniques. We discuss some observational consequences of these effects for gravitational-wave detection.

We present, in closed analytic form, a general stationary, slowly rotating black hole, which is solution to a large class of alternative theories of gravity in four dimensions. In these theories, the Einstein-Hilbert action is supplemented by all possible quadratic, algebraic curvature invariants coupled to a scalar field. The solution is found as a deformation of the Schwarzschild metric in General Relativity. We explicitly derive the changes to the orbital frequency at the innermost stable circular orbit and at the light ring in closed form. These results could be useful when comparing General Relativity against alternative theories by (say) measurements of X-ray emission in accretion disks, or by stellar motion around supermassive black holes. When gravitational-wave astronomy comes into force, strong constraints on the coupling parameters can in principle be made.