Search

A generalization of the Majumdar-Papapetrou multi-black hole spacetime was recently constructed by Teo and Wan [1], describing charged and spinning (extremal) balanced black holes in asymptotically flat spacetime. We explore the dynamics of null geodesics on this geometry, focusing on the two-center solution. Using the topological charge formalism, we show that various light ring arrangements arise from different choices of individual angular momenta: light rings with opposite topological charges can merge and annihilate each other, resulting in configurations with a total of 4, 6, or 8 light rings. Using backward ray-tracing, we obtained the shadow and lensing of these spacetimes. The former, in particular, closely resembles those for the double-Kerr metric., Comment: 16 pages, 15 figures

We study the bifurcation phenomena between spherical and axisymmetric bosonic stars. By numerically solving for the zero-modes of spherical bosonic stars under specific axially symmetric perturbations, we discover that excited state spherical bosonic stars bifurcate into two types of axisymmetric bosonic stars under $\ell=2$ perturbations, with matter distributions resembling chains and rings, respectively. Meanwhile, $\ell=4$ axisymmetric perturbations lead spherical scalar bosonic stars to bifurcate into a new type of axisymmetric bosonic stars, exhibiting a mixed chain-like and ring-like matter distribution, which we refer to as gyroscope-like. Additionally, for the first time, we have constructed chains of scalar bosonic stars with 7 constituents and their corresponding ring-like scalar bosonic stars. Our results provide an explanation for the bifurcations in bosonic stars from the perspective of perturbations, and by analyzing physical quantities such as quadrupoles and energy densities we systematically discuss the impact of axisymmetric perturbations on spherical bosonic stars., Comment: 21 pages, 8 figures

Exact solutions describing multiple, electrically charged black holes (BHs) in a model of non-linear electrodynamics (NLE) minimally coupled to Einstein's gravity are presented. The NLE model is ModMax theory, that has attracted much attention due to its duality and conformal invariance, features shared with standard (linear) electrodynamics. In the non-extremal case, the solution has conical singularities, similarly to the multi Reissner-Nordstr\"om solution in Einstein-Maxwell theory. In the extremal case the solution is regular on and outside the event horizon; it is isometric to the Majumdar-Papapetrou solution, albeit the individual BHs have a non-unitary charge to mass ratio, due to screening effects. Using the ModMax electromagnetic duality invariance, magnetically charged and dyonic generalizations are also obtained. Finally, we construct multi-BH solutions with a positive cosmological constant., Comment: 9 pages

In relativistic mechanics, the 4-velocity and the 4-momentum need not be parallel. This allows their norm to have a different sign. This possibility occurs in nonlinear electrodynamics (NED) models minimally coupled to Einstein's theory. Surprisingly, for a large class of NED models with a Maxwell limit, for weak fields, the causal (acausal) photons, as determined by their 4-velocity, have a spacelike (timelike) 4-momentum, leading to good tachyons and bad bradyons. Departing from weak fields, this possibility is determined solely by the concavity of the NED Lagrangian, which is consistent with the Dominant Energy Condition analysis. As a corollary, some popular regular black hole solutions sourced by NED, such as the Bardeen and Hayward solutions, are acausal., Comment: 6 pages, 1 figure, 1 table

We study stationary clouds of a gauged, complex scalar field on a magnetically (and possibly electrically as well) charged Kerr-Newman black hole (BH). The existence of a magnetic charge $Q_m$ promotes a north-south $\textit{asymmetry}$ of the scalar clouds. This breakdown of the clouds' $\mathbb{Z}_2$-symmetry carries through to the spacetime $\textit{geometry}$ for the non-linear continuation of the clouds: a family of magnetically charged (or dyonic) BHs with synchronized gauged scalar hair, which we construct. Their distinct phenomenology is illustrated by their imaging, exhibiting skewed shadows and lensing. Such hairy BHs could, in principle, result from the superradiant instability of magnetically charged Kerr-Newman BHs, unveiling a dynamical mechanism for creating north-south asymmetric BHs from standard $\mathbb{Z}_2$-symmetric electrovacuum BHs., Comment: 19 pages, 10 figures

The gravitational physics landscape is evolving rapidly, driven by our ability to study strong-field regions, in particular black holes. Black Holes Inside and Out gathered world experts to discuss the status of the field and prospects ahead. We hope that the ideas and perspectives are a source of inspiration. Structure: Black Hole Evaporation - 50 Years by William Unruh The Stability Problem for Extremal Black Holes by Mihalis Dafermos The Entropy of Black Holes by Robert M. Wald The Non-linear Regime of Gravity by Luis Lehner Black Holes Galore in D > 4 by Roberto Emparan Same as Ever: Looking for (In)variants in the Black Holes Landscape by Carlos A. R. Herdeiro Black Holes, Cauchy Horizons, and Mass Inflation by Matt Visser The Backreaction Problem for Black Holes in Semiclassical Gravity by Adrian del Rio Black Holes Beyond General Relativity by Enrico Barausse and Jutta Kunz Black Holes as Laboratories: Searching for Ultralight Fields by Richard Brito Primordial Black Holes from Inflation by Misao Sasaki Tests of General Relativity with Future Detectors by Emanuele Berti Black Holes as Laboratories: Tests of General Relativity by Ruth Gregory and Samaya Nissanke Simulating Black Hole Imposters by Frans Pretorius Black Hole Spectroscopy: Status Report by Gregorio Carullo VLBI as a Precision Strong Gravity Instrument by Paul Tiede Testing the nature of compact objects and the black hole paradigm by Mariafelicia De Laurentis and Paolo Pani Some Thoughts about Black Holes in Asymptotic Safety by Alessia Platania Black Hole Evaporation in Loop Quantum Gravity by Abhay Ashtekar How the Black Hole Puzzles are Resolved in String Theory by Samir D. Mathur Quantum Black Holes: From Regularization to Information Paradoxes by Niayesh Afshordi and Stefano Liberati, Comment: 221 pages, 21 contributions

Recently, a scalar counterpart of the Schwarzschild-Melvin Universe was reported [arXiv:2410.02851]. We show this solution is a special case of a Schwarzschild black hole/mass in a scalar multipolar Universe, that can be constructed algebraically combining known vacuum solutions. This builds on the generalized Weyl construction for scalar-vacuum, that admits two, fully decoupled, independent harmonic functions: one for the gravitational sector and another for the scalar sector. Harmonic solutions with growing multipoles lead to naked singularities (at spatial infinity), including for the scalar counterpart of the Schwarzschild-Melvin; harmonic solutions with decaying (scalar) multipoles are asymptotically flat, but have a singular horizon, in accordance with no-scalar-hair theorems. The scalar-vacuum model can be mapped to a five dimensional pure gravity construction via Kaluza-Klein oxidization and, in this way, to the corresponding $D=5$ generalized Weyl construction., Comment: 9 pages; v2 Typos corrected. Matches version to appear in Phys. Lett. B

What does a black hole look like? In 1+3 spacetime dimensions, the optical appearance of a black hole is a bidimensional region in the observer's sky often called the black hole shadow, as supported by the EHT observations. In higher dimensions this question is more subtle and observational setup dependent. Previous studies considered the shadows of higher dimensional black holes to remain bidimensional. We argue that the latter should be regarded as a tomography of a higher dimensional structure, the hypershadow, which would be the structure "seen" by higher dimensional observers. As a case study we consider the cohomogeneity-one Myers-Perry black hole in 1+4 dimensions, and compute its tridimensional hypershadow., Comment: 23 pages, 5 figures

Some of the extensions to general relativity and to the Standard Model of particle physics predict families of hypothetical compact objects, collectively known as exotic compact objects (ECOs). This category can be defined to encompass non-Kerr black holes both within and beyond general relativity, as well as horizonless compact objects such as boson stars. In order to model observational signatures and identify possible detections, it is crucial to understand the interaction between these objects and their surrounding medium, usually plasmas described by the equations of general relativistic magnetohydrodynamics (GRMHD). To this end, we review the existent literature on GRMHD simulations of accretion onto these objects. These cover a variety of objects and accretion patterns. We conclude by listing possible directions to continue exploring this relatively young field., Comment: 37 pages, 8 figures, Invited chapter for the edited book "New Frontiers in GRMHD Simulations" (Eds. C. Bambi, Y. Mizuno, S. Shashank and F. Yuan, Springer Singapore, expected in 2025). Updated acknowledgements in version 2

For gauged scalar fields minimally coupled to Einstein-Maxwell theory, the Mayo-Bekenstein no-hair theorem can be circumvented when including appropriate scalar self-interactions, allowing static, electrically charged black holes to be endowed with (Abelian) gauged scalar hair. Here we show these spherically symmetric solutions can be extended to include a magnetic charge in a model with scalar multiplets. The resulting dyonic configurations share most of the properties of the electrically charged solutions, in particular satisfying the same {\it resonance} condition, with the existence of a mass gap with respect to the bald Reissner-Nordstr\"om dyonic black holes. A distinctive feature, however, is that no solitonic limit exists for a non-zero magnetic charge., Comment: 12 pages, 5 figures

Bosonic stars represent a hypothetical exotic type of compact stellar objects that could be observed from the gravitational signal of coalescing binaries in current and future gravitational wave detectors. There are two main families of bosonic stars, which depend on the nature that governs the particles that build them: Einstein-Klein-Gordon and Proca Stars. We study the multipolar structure for both families of rotating objects, using realistic potentials with the aim of finding possible universal relations and, thus, a method that allows us to distinguish between these and other compact objects in the gravitational wave paradigm. We also show how certain relevant observables can be obtained for these hypothetical but well-motivated astrophysical objects., Comment: 23 pages, 22 figures. We made some corrections

Recently, spherical and static flat space solitons (balls) and self-gravitating, everywhere regular, asymptotically flat solitons (stars) were constructed in an Einstein-Proca-Higgs model [1], where a complex vector field gains mass by coupling to a real scalar field with a Higgs-type potential. The Proca-Higgs model serves as a UV completion of a complex Proca model with self-interactions. Here, we construct and examine the mathematical and physical properties of rotating configurations. In particular, rotation allows horizon-bearing solutions, including stationary clouds surrounding Kerr black holes and their non-linear continuation into black holes with Proca-Higgs hair., Comment: 57 pages

We investigate the static boson star solutions in the so-called Einstein-Friedberg-Lee-Sirlin (E-FLS) theory, performing a complete analysis of the solution space in this model. We study the phenomenological aspects of E-FLS stars, for instance, by investigating the timelike and null geodesics with an emphasis on the analysis of circular timelike orbits and light rings. In order to study the astrophysical signatures of such stars, their images were obtained considering them surrounded by a geometrically thin accretion disk. Our results comprise two different models of accretion disks, namely the optically thin and thick disk models. We present a selection of our findings for the astrophysical images of E-FLS stars and discuss their relevance as a possible black hole mimicker., Comment: 16 pages, 13 figures

Control algorithms have been proposed based on knowledge related to nature-inspired mechanisms, including those based on the behavior of living beings. This paper presents a review focused on major breakthroughs carried out in the scope of applied control inspired by the gravitational attraction between bodies. A control approach focused on Artificial Potential Fields was identified, as well as four optimization metaheuristics: Gravitational Search Algorithm, Black-Hole algorithm, Multi-Verse Optimizer, and Galactic Swarm Optimization. A thorough analysis of ninety-one relevant papers was carried out to highlight their performance and to identify the gravitational and attraction foundations, as well as the universe laws supporting them. Included are their standard formulations, as well as their improved, modified, hybrid, cascade, fuzzy, chaotic and adaptive versions. Moreover, this review also deeply delves into the impact of universe-inspired algorithms on control problems of dynamic systems, providing an extensive list of control-related applications, and their inherent advantages and limitations. Strong evidence suggests that gravitation-inspired and black-hole dynamic-driven algorithms can outperform other well-known algorithms in control engineering, even though they have not been designed according to realistic astrophysical phenomena and formulated according to astrophysics laws. Even so, they support future research directions towards the development of high-sophisticated control laws inspired by Newtonian/Einsteinian physics, such that effective control-astrophysics bridges can be established and applied in a wide range of applications., Comment: This is a preprint of a paper that is published open access in 'Heliyon' [https://doi.org/10.1016/j.heliyon.2024.e31771]

The Einstein field equations, or generalizations thereof, are difficult to solve analytically. On the other hand, numerical solutions of the same equations have become increasingly common, in particular concerning compact objects. Whereas analytic approximations to each individual solution within a numerical family have been proposed, proxies for whole families are missing, which can facilitate studying properties across the parameter space, data compression and a wider usage of such solutions. In this work we tackle this need, proposing a simple strategy based on two different expansions of the unknown functions in an appropriately chosen basis, to build such proxy. We use as an exploratory case-study spherical, fundamental mini-boson stars, to illustrate the feasibility of such an approach, emphasise its advantage in reducing the data size, and the challenges, say, in covering large parameter spaces., Comment: 28 pages, 10 figures

We consider a model with two real Maxwell fields (or equivalently, a complex Maxwell field) minimally coupled to Einsteins gravity with a negative cosmological constant in four spacetime dimensions. Assuming a specific harmonic dependence of the vector fields, we show the existence of asymptotically anti-de Sitter (AdS) self-gravitating boson-star-like solitonic solutions, which are static and axially symmetric. Analytical solutions are found in the test-field limit, where the Maxwell equations are solved on a fixed AdS background. The fully nonlinear solutions are constructed numerically., Comment: 13 pages, 2 figures

We study the dynamical stability of Proca-Higgs stars, in spherical symmetry. These are solutions of the Einstein-Proca-Higgs model, which features a Higgs-like field coupled to a Proca field, both of which minimally coupled to the gravitational field. The corresponding stars can be regarded as Proca stars with self-interactions, while avoiding the hyperbolicity issues of self-interacting Einstein-Proca models. We report that these configurations are stable near the Proca limit in the candidate stable branches, but exhibit instabilities in certain parts of the parameter space, even in the candidate stable branches, regaining their stability for very strong self-interactions. This shows that for these models, unlike various examples of scalar boson stars, self-interactions can deteriorate, rather than improve, the dynamical robustness of bosonic stars., Comment: 14 pages, 11 figures, updated to match the published version

We study radial perturbations of spherically symmetric spin-$0$ and spin-$1$ bosonic stars, computing numerically the squared frequency of the fundamental mode. We find that not all critical points $-$ where the Arnowitt-Deser-Misner mass attains an extremum $-$ correspond to zero modes. Thus, radial stability does not $\textit{always}$ change at such critical points. The results are in agreement with the so-called critical point method., Comment: 23 pages, 8 figures and 6 tables

Ultracompact objects (UCOs) are horizonless compact objects that present light rings (LRs) - circular photon orbits. As a result, they could be black hole mimickers. Some years ago, Cunha et al. established a theorem stating that, under general assumptions, UCOs formed from smooth, quasi-Minkowski initial data, must have at least a pair of LRs, one of which must be stable. These stable LRs are supposed to trigger a non-linear instability in spacetime, potentially weakening UCOs' ability to replicate black hole phenomenology. However, this LR theorem does not extend to wormholes, which represent topologically nontrivial spacetimes. We address the wormhole case by proving the following theorem: a stationary, axisymmetric, asymptotically flat, traversable wormhole in 1+3 dimensions, connecting two different asymptotic regions, has at least one standard LR for each rotation sense. Thus, any (such) wormhole is an UCO. By filling this gap, our results not only broaden the horizon of knowledge on UCOs but also highlight their potential to closely mimic black hole phenomenology., Comment: 9 pages, 3 figures

We present several high-accuracy surrogate models for gravitational-wave signals from equal-mass head-on mergers of Proca stars, computed through the Newman-Penrose scalar $\psi_4$. We also discuss the current state of the model extensions to mergers of Proca stars with different masses, and the particular challenges that these present. The models are divided in two main categories: two-stage and monolithic. In the two-stage models, a dimensional reduction algorithm is applied to embed the data in a reduced feature space, which is then interpolated in terms of the physical parameters. For the monolithic models, a single neural network is trained to predict the waveform from the input physical parameter. Our model displays mismatches below $10^{-3}$ with respect to the original numerical waveforms. Finally, we demonstrate the usage of our model in full Bayesian parameter inference through the accurate recovery of numerical relativity signals injected in zero-noise, together with the analysis of GW190521. For the latter, we observe excellent agreement with existing results that make use of full numerical relativity., Comment: 16 pages, 14 figures. v2: Matches published version

We analyse the lensing images by dynamically robust rotating (mini-)Proca stars surrounded by thin accretion disks. Due to their peculiar geodesic structure we show that these images exhibit striking similarities with the ones of BHs, for appropriately chosen disk intensity profile, when imposing a GRMHD-motivated emission cut off. Additionally, and unlike the non-rotating case, these similarities prevail even when considering equatorial observations. This example illustrates how a horizonless compact object without light rings, with a plausible formation mechanism and dynamically robust, could mimic detailed features of black hole imagiology., Comment: 23 pages, 11 figures + appendix

Using the Ernst formalism, a novel solution of vacuum General Relativity was recently obtained [1], describing a Schwarzschild black hole (BH) immersed in a non-asymptotically flat rotating background, dubbed swirling universe, with the peculiar property that north and south hemispheres spin in opposite directions. We investigate the null geodesic flow and, in particular, the existence of light rings in this vacuum geometry. By evaluating the total topological charge $w$, we show that there exists one unstable light ring ($w=-1$) for each rotation sense of the background. We observe that the swirling background drives the Schwarzschild BH light rings outside the equatorial plane, displaying counter-rotating motion with respect to each other, while (both) co-rotating with respect to the swirling universe. Using backwards ray-tracing, we obtain the shadow and gravitational lensing effects, revealing a novel feature for observers on the equatorial plane: the BH shadow displays an odd $\mathbb{Z}_2$ (north-south) symmetry, inherited from the same type of symmetry of the spacetime itself: a twisted shadow., Comment: 21 pages, 14 figures

It has been established that Black Hole (BH) spacetimes obeying some general set of assumptions always possess, at least, one light ring (per rotation sense) [arXiv:2003.06445]. This theorem was originally established for asymptotically flat, stationary, axial symmetric, 1+3 dimensional circular spacetimes harbouring a non-extremal and topologically spherical Killing horizon. Following the mantra that a theorem is only as strong as its assumptions in this work we extend this theorem to non topologically spherical (toroidal) BHs and to spacetimes harbouring more than one BH. As in [arXiv:2003.06445], we show that each BH still contributes with, at least, one LR (per rotation sense)., Comment: 17 pages, 9 figures

We construct a new family of rotating black holes with scalar hair and a regular horizon of spherical topology, within five dimensional ($d=5$) Einstein's gravity minimally coupled to a complex, massive scalar field doublet. These solutions represent generalizations of the Kaluza-Klein monopole found by Gross, Perry and Sorkin, with a twisted $S^1$ bundle over a four dimensional Minkowski spacetime being approached in the far field. The black holes are described by their mass, angular momentum, tension and a conserved Noether charge measuring the hairiness of the configurations. They are supported by rotation and have no static limit, while for vanishing horizon size, they reduce to boson stars. When performing a Kaluza-Klein reduction, the $d=5$ solutions yield a family of $d=4$ spherically symmetric dyonic black holes with gauged scalar hair. This provides a link between two seemingly unrelated mechanisms to endow a black hole with scalar hair: the $d=5$ synchronization condition between the scalar field frequency and the event horizon angular velocity results in the $d=4$ resonance condition between the scalar field frequency and the electrostatic chemical potential.

In linear perturbation theory, a static perturber in the vicinity of a Schwarzschild black hole (BH) enhances [suppresses] the Gauss-Bonnet (GB) curvature invariant, $\mathcal{R}_{\rm GB}$, in the high [low] tide regions. By analysing exact solutions of the vacuum Einstein field equations describing one or two BHs immersed in a multipolar gravitational field, which is locally free of pathologies, including conical singularities, we study the corresponding non-linear tides on a fiducial BH, in full General Relativity (GR). We show that the tidal field due to a far away, or close by, static BH creates high/low tides that can deviate not only quantitatively but also qualitatively from the weak field/Newtonian pattern. Remarkably, the suppression in low tide regions never makes $\mathcal{R}_{\rm GB}$ negative on the BH, even though the horizon Gaussian curvature may become negative; but $\mathcal{R}_{\rm GB}$ can vanish in a measure zero set, a feature qualitatively recovered in a Newtonian analogue model. Thus, purely gravitational, static, tidal interactions in GR, no matter how strong, cannot induce GB$^-$ scalarization. We also show that a close by BH produces noticeable asymmetric tides on another (fiducial) BH., Comment: 13 pages, 13 figures

We construct and dynamically evolve dipolar, self-interacting scalar boson stars in a model with sextic (+ quartic) self-interactions. The domain of existence of such dipolar $Q$-stars has a similar structure to that of the fundamental monopolar stars of the same model. For the latter it is structured in a Newtonian plus a relativistic branch, wherein perturbatively stable solutions exist, connected by a middle unstable branch. Our evolutions support similar dynamical properties of the dipolar $Q$-stars that: 1) in the Newtonian and relativistic branches are dynamically robust over time scales longer than those for which dipolar stars without self-interactions are seen to decay; 2) in the middle branch migrate to either the Newtonian or the relativistic branch; 3) beyond the relativistic branch decay to black holes. Overall, these results strengthen the observation, seen in other contexts, that self-interactions can mitigate dynamical instabilities of scalar boson star models., Comment: 13 pages, 14 figures; movies of the numerical simulations reported can be found in http://gravitation.web.ua.pt/index.php/node/4488

General Relativity minimally coupled to a massive, free, complex scalar field, is shown to allow asymptotically flat solutions, non-singular on and outside the event horizon, describing two spinning black holes (2sBHs) in equilibrium, with co-axial, aligned angular momenta. The 2sBHs configurations bifurcate from solutions describing dipolar spinning boson stars. The BHs emerge at equilibrium points diagnosed by a test particle analysis and illustrated by a Newtonian analogue. The individual BH "charges" are mass and angular momentum only. Equilibrium is due to the scalar environment, acting as a (compact) dipolar field, providing a lift against their mutual attraction, making the 2sBHs (h)airborne. We explore the 2sBHs domain of solutions and its main features., Comment: 5 pages, 5 figures

We study extreme mass ratio inspirals (EMRIs) due to an infalling Light Compact Object (LCO) onto a generic class of stationary and axi-symmetric massive compact objects (MCO - with or without a horizon). Using the quadrupole hybrid formalism we obtain a master formula for the evolution of the radius of the LCO and find qualitatively different behaviours depending on the geodesic structure of the MCO. We then specialize the MCO to a black hole with synchronised scalar hair (BHsSH). To allow a comparison with a highly spinning Kerr BH, we consider BHsSH with dimensionless spin, $j = 1$. This yields two distinct sequences of solutions. The first harbours Kerr-like solutions with maximal hairiness of $\sim 10\%$. The corresponding EMRIs are Kerr-like, but the cut-off frequency can be a few times smaller than in Kerr, yielding waveforms with quantitatively significant non-Kerrness. The second sequence links the extremal Kerr black hole to a mini-boson star with $j=1$. Here we observe qualitative non-Kerrness, such as the non-monotonically increase of the angular velocity and stagnation endpoints, reflecting Kerr-unlike geodesic structures., Comment: 20 pages, 15 figures

We analyze the quasinormal modes (QNMs) of a recently obtained solution of a Schwarzschild black hole (BH) with corrections motivated by Loop Quantum Gravity (LQG). This spacetime is regular everywhere and presents the global structure of a wormhole, with a minimal surface whose radius depends on a LQG parameter. We focus on the investigation of massless scalar field perturbations over the spacetime. We compute the QNMs with the WKB approximation, as well as the continued fraction method. The QNM frequency orbits, for $l=0$ and $n>0$, where $l$ and $n$ are the multipole and overtone numbers, respectively, are self-intersecting, spiraling curves in the complex plane. These orbits accumulate to a fixed complex value corresponding to the QNMs of the extremal case. We obtain that, for small values of the LQG parameter, the overall damping decreases as we increase the LQG parameter. Moreover the spectrum of the quantum corrected black hole exhibits an oscillatory pattern, which might imply in the existence of QNMs with vanishing real part. This pattern suggests that the limit $n\rightarrow \infty$ for the real part of the QNMs is not well-defined, what differs from Schwarzschild's case. We also analyze the time-domain profiles for the scalar perturbations, showing that the LQG correction does not alter the Schwarzschild power-law tail. We compute the fundamental mode from the time profile by means of the Prony method, obtaining excellent agreement with the two previously mentioned methods., Comment: 12 pages, 7 figures

We study the time evolution of spherical, excited -- with $n$ radial nodes -- scalar boson stars in General Relativity minimally coupled to a complex massive scalar field with quartic self-interactions. We report that these stars, with up to $n=10$, can be made dynamically stable, up to timescales of $t\sim\frac{10^{4}}{c\mu}$, where $\mu$ is the inverse Compton wavelength of the scalar particle, for sufficiently large values of the self-interactions coupling constant $\lambda$, which depend on $n$. We observe that the compactness of these solutions is rather insensitive to $n$, for large $\lambda$ and fixed frequency. Generically, along the branches where stability was studied, these excited boson stars are not compact enough to allow for innermost stable circular orbits or light rings. Finally, we discuss the angular velocity of particles along timelike circular orbits, suggesting an application, for solutions in the Newtonian limit, to galactic rotation curves., Comment: 16 pages, 16 figures, updated with an appendix to match published version

We show that, unlike vacuum General Relativity, Einstein-scalar theories allow balanced static, neutral, asymptotically flat, double-black hole solutions, for scalar field models minimally coupled to gravity, with appropriate self-interactions. These are scalar hairy versions of the double-Schwarzschild (or Bach-Weyl) solution, but regular on and outside the two (topologically spherical) horizons. The balancing repulsive force is provided by the scalar field. An explicit illustration is presented, using a Weyl-type construction adapted to numerical solutions, requiring no partial linearisation, or integrability structure, of the Einstein-scalar equations. Fixing the couplings of the model, the balanced configurationsform a one-parameter family of solutions, labelled by the proper distance between the black holes., Comment: 22 pages, 8 figures

We present a numerical study of relativistic Bondi-Hoyle-Lyttleton (BHL) accretion onto an asymptotically flat black hole with synchronized hair. The hair is sourced by an ultralight, complex scalar field, minimally coupled to Einstein's gravity. Our simulations consider a supersonic flow parametrized by the asymptotic values of the fluid quantities and a sample of hairy black holes with different masses, angular momenta, and amount of scalar hair. For all models, steady-state BHL accretion solutions are attained that are characterized by the presence of a shock-cone and a stagnation point downstream. For the models of the sample with the largest component of scalar field, the shock-cone envelops fully the black hole, transitioning into a bow-shock, and the stagnation points move further away downstream. Analytical expressions for the mass accretion rates are obtained after fitting the numerical results, which can be used to analyze black-hole formation scenarios in the presence of ultralight scalar fields. The formation of a shock-cone leads to regions where sound waves can be trapped and resonant oscillations excited. We measure the frequencies of such quasi-periodic oscillations and point out a possible association with quasi-periodic oscillations in the X-ray light curve of Sgr~A* and microquasars., Comment: 24 pages, 3 tables, 10 figures. Minor changes

We consider a Proca-Higgs model wherein a complex vector field gains mass via spontaneous symmetry breaking, by coupling to a real scalar field with a Higgs-type potential. This vector version of the scalar Friedberg-Lee-Sirlin model, can be considered as a UV completion of a complex Proca model with self-interactions. We study the flat spacetime and self-gravitating solitons of the model, that we dub Proca-Higgs \textit{balls} and \textit{stars} respectively, exploring the domain of solutions and describing some of their mathematical and physical properties. The stars reduce to the well-known (mini-)Proca stars in some limits. The full model evades the hyperbolicity problems of the self-interacting Proca models, offering novel possibilities for dynamical studies beyond mini-Proca stars.

We investigate the gravitational lensing by spinning Proca stars and the shadows and lensing by Kerr black holes (BHs) with synchronised Proca hair, discussing both theoretical aspects and observational constraints from the Event Horizon Telescope (EHT) M87* and Sgr A* data. On the theoretical side, this family of BHs interpolates between Kerr-like solutions -- exhibiting a similar optical appearance to that of Kerr BHs -- to very non-Kerr like solutions, exhibiting exotic features such as cuspy shadows, egg-like shadows and ghost shadows. We interpret these features in terms of the structure of the fundamental photon orbits, for which different branches exist, containing both stable and unstable orbits, with some of the latter not being shadow related. On the observational side, we show that current EHT constraints are compatible with all such BHs that could form from the growth of the superradiant instability of Kerr BHs. Unexpectedly, given the (roughly) 10% error bars in the EHT data -- and in contrast to their scalar cousin model --, some of the BHs with up to 40% of their energy in their Proca hair are compatible with the current data. We estimate the necessary resolution of future observations to better constrain this model., Comment: 36 pages, 20 figures + appendix

We present a systematic study of the dynamics and gravitational-wave emission of head-on collisions of spinning vector boson stars, known as Proca stars. To this aim we build a catalogue of about 800 numerical-relativity simulations of such systems. We find that the wave-like nature of bosonic stars has a large impact on the gravitational-wave emission. In particular, we show that the initial relative phase $\Delta \epsilon =\epsilon_1-\epsilon_2$ of the two complex fields forming the stars (or equivalently, the relative phase at merger) strongly impacts both the emitted gravitational-wave energy and the corresponding mode structure. This leads to a non-monotonic dependence of the emission on the frequency of the secondary star $\omega_2$, for fixed frequency $\omega_1$ of the primary. This phenomenology, which has not been found for the case of black-hole mergers, reflects the distinct ability of the Proca field to interact with itself in both constructive and destructive manners. We postulate this may serve as a smoking gun to shed light on the possible existence of these objects., Comment: Updated version

The motion of particles on spherical $1 + 3$ dimensional spacetimes can, under some assumptions, be described by the curves on a 2-dimensional manifold, the optical and Jacobi manifolds for null and timelike curves, respectively. In this paper we resort to auxiliary 2-dimensional metrics to study circular geodesics of generic static, spherically symmetric, and asymptotically flat $1 + 3$ dimensional spacetimes, whose functions are at least $C^2$ smooth. This is done by studying the Gaussian curvature of the bidimensional equivalent manifold as well as the geodesic curvature of circular paths on these. This study considers both null and timelike circular geodesics. The study of null geodesics through the optical manifold retrieves the known result of the number of light rings (LRs) on the spacetime outside a black hole and on spacetimes with horizonless compact objects. With an equivalent procedure we can formulate a similar theorem on the number of marginally stable timelike circular orbits (TCOs) of a given spacetime satisfying the previously mentioned assumption, Comment: 7 pages + Appendices

Ultracompact objects with light-rings (LRs) but without an event horizon could mimic black holes (BHs) in their strong gravity phenomenology. But are such objects dynamically viable? Stationary and axisymmetric ultracompact objects that can form from smooth, quasi-Minkowski initial data must have at least one stable LR, which has been argued to trigger a spacetime instability; but its development and fate have been unknown. Using fully non-linear numerical evolutions of ultracompact bosonic stars free of any other known instabilities and introducing a novel adiabatic effective potential technique, we confirm the LRs triggered instability, identifying two possible fates: migration to non-ultracompact configurations or collapse to BHs. In concrete examples we show that typical migration/collapse time scales are not larger than $\sim 10^3$ light-crossing times, unless the stable LR potential well is very shallow. Our results show that the LR instability is effective in destroying horizonless ultracompact objects that could be plausible BH imitators., Comment: 5 pages + Appendices; Videos can be found in this URL http://gravitation.web.ua.pt/node/3889; v2. Title changed per request of journal. Published in Physical Review Letters as Editor's Suggestion and Featured in Physics

A mechanism for circumventing the Mayo-Bekenstein no-hair theorem allows endowing four dimensional $(D=4)$ asymptotically flat, spherical, electro-vacuum black holes with a minimally coupled $U(1)$-gauged scalar field profile: $Q$-$hair$. The scalar field must be massive, self-interacting and obey a {\it resonance condition} at the threshold of (charged) superradiance. We establish generality for this mechanism by endowing three different types of static black objects with scalar hair, within a $D=5$ Einstein-Maxwell-gauged scalar field model: asymptotically flat black holes and black rings; and black strings which asymptote to a Kaluza-Klein vacuum. These $D=5$ $Q$-hairy black objects share many of the features of their $D=4$ counterparts. In particular, the scalar field is subject to a resonance condition and possesses a $Q$-ball type potential. For the static black ring, the charged scalar hair can balance it, yielding solutions that are singularity free on and outside the horizon., Comment: 30 pages, 8 figures

The local thermodynamic stability of a black hole (BH) in the canonical ensemble is defined by the positivity of the specific heat at constant global charges. Schwarzschild BHs in thermodynamic equilibrium with an energy reservoir are always unstable against small fluctuations of energy, whereas sufficiently near-extremal Reissner-Nordstr\"{o}m/Kerr BHs are stable. One could expect that asymptotically-flat hairy BHs branching off from such stable phases would also be, by continuity, locally thermodynamically stable for vanishingly little hair. We show this is not the case in some models, including scalarized BHs bifurcating from Reissner-Nordstr\"{o}m and spinning BHs with synchronized hair bifurcating from Kerr. Specifically, it is found that quasi-bald BHs are locally thermodynamically unstable in the canonical ensemble for all global charges and regardless of being dynamically and entropically preferred over bald ones at fixed global charges., Comment: 11 pages, 4 figures. Comments are welcome

Derrick-type virial identities, obtained via dilatation (scaling) arguments, have a variety of applications in field theories. We deconstruct such virial identities in relativistic gravity showing how they can be recast as self-evident integrals of appropriate combinations of the equations of motion. In spherical symmetry, the appropriate combination and gauge choice guarantee the geometric part can be integrated out to yield a master form of the virial identity as a non-trivial energy-momentum balance condition, valid for both asymptotically flat black holes and self-gravitating solitons, for any matter model. Specifying the matter model we recover previous results obtained via the scaling procedure. We then discuss the more general case of stationary, axi-symmetric, asymptotically flat black hole or solitonic solutions in General Relativity, for which a master form for their virial identity is proposed, in a specific gauge but regardless of the matter content. In the flat spacetime limit, the master virial identity for both the spherical and axial cases reduces to a balance condition for the principal pressures, discussed by Deser., Comment: 14 pages

We present the first systematic search for exotic compact mergers in Advanced LIGO and Virgo events. We compare the short gravitational-wave signals GW190521, GW190426$\_$190642, GW200220$\_$061928 and the trigger 200114$\_$020818 (or S200114f) to a new catalogue of 759 numerical simulations of head-on mergers of horizonless exotic compact objects known as Proca stars, interpreted as self-gravitating lumps of (fuzzy) dark matter sourced by an ultralight (vector) bosonic particle. The Proca-star merger hypothesis is strongly rejected with respect to the black hole merger one by GW190426, weakly rejected by GW200220 and weakly favoured by GW190521 and S200114f. GW190521 and GW200220 yield highly consistent boson masses of $\mu_{\rm B} = 8.69^{+0.61}_{-0.75}\times10^{-13}$ eV and $\mu_{\rm B} = 9.13^{+1.18}_{-1.30}\times10^{-13}$ eV at the $90\%$ credible level. We conduct a preliminary population study of the compact binaries behind these events. Excluding (including) S200114f as a real event, and ignoring boson-mass consistencies across events, we estimate a fraction of Proca-star mergers of $\zeta = 0.27^{+0.43}_{-0.25} \ (0.39^{+0.38}_{-0.33})$. We discuss the impact of boson-mass consistency across events in such estimates. Our results maintain GW190521 as a Proca-star merger candidate and pave the way towards population studies considering exotic compact objects., Comment: Version accepted in Phys.Rev.D. 27 pages, 16 Figures, 4 Appendixes

Detection and parameter inference of gravitational-wave signals \ncor{from compact mergers} rely on the comparison of the incoming detector strain data $d(t)$ to waveform templates for the gravitational-wave strain $h(t)$ that ultimately rely on the resolution of Einstein's equations via numerical relativity simulations. These, however, commonly output a quantity known as the Newman-Penrose scalar $\psi_4(t)$ which, under the Bondi gauge, is related to the gravitational-wave strain by $\psi_4(t)=\mathrm{d}^2h(t) / \mathrm{d}t^2$. Therefore, obtaining strain templates involves an integration process that introduces artefacts that need to be treated in a rather manual way. By taking second-order finite differences on the detector data and inferring the corresponding background noise distribution, we develop a framework to perform gravitational-wave data analysis directly using $\psi_4(t)$ templates. We first demonstrate this formalism, and the impact of integration artefacts in strain templates, through the recovery of numerically simulated signals from head-on collisions of Proca stars injected in Advanced LIGO noise. Next, we re-analyse the event GW190521 under the hypothesis of a Proca-star merger, obtaining results equivalent to those in Ref.[1], where we used the classical strain framework. We find, however, that integration errors would strongly impact our analysis if GW190521 was four times louder. Finally, we show that our framework fixes significant biases in the interpretation of the high-mass GW trigger S200114f arising from the usage of strain templates. We remove the need to obtain strain waveforms from numerical relativity simulations, avoiding the associated systematic errors., Comment: Version accepted for publication in Physical Review X. 27 pages, 11 Figures, 3 Appendixes

Because the Klein-Gordon and Proca equations involve $\hbar$, they describe quantum fields. Their solutions, however, may be treated as classical if their typical action obeys $S^{\rm typical}\gg \hbar$. This is possible due to their bosonic nature, allowing states with many particles. We show, by generic arguments, that the typical action for such bosonic stars is ${\mathcal{S}^{\rm typical}}/{\hbar}\gtrsim \left({M^{\rm max}}/{M_{\rm Pl}}\right)^2\sim 10^{76}\left({M^{\rm max}}/{M_\odot}\right)^2$, where $M^{\rm max}$ is the maximal bosonic star mass of the particular model and $M_{\rm Pl}$ is the Planck mass. Thus, for models allowing $M^{\rm max}\gg M_{\rm Pl}$ and for solutions with mass $\sim M^{\rm max}$, the classical treatment is legitimate, which includes masses in the astrophysical interesting range $\gtrsim M_\odot$., Comment: 11 pages; Essay written for the Gravity Research Foundation 2022 Awards for Essays on Gravitation and selected for Honorable Mention

The Laser Interferometer Space Antenna (LISA) has the potential to reveal wonders about the fundamental theory of nature at play in the extreme gravity regime, where the gravitational interaction is both strong and dynamical. In this white paper, the Fundamental Physics Working Group of the LISA Consortium summarizes the current topics in fundamental physics where LISA observations of GWs can be expected to provide key input. We provide the briefest of reviews to then delineate avenues for future research directions and to discuss connections between this working group, other working groups and the consortium work package teams. These connections must be developed for LISA to live up to its science potential in these areas., Comment: Accepted in: Living Reviews in Relativity

We examine the constraining power of current gravitational-wave data on scalar-tensor-Gauss-Bonnet theories that allow for the spontaneous scalarization of black holes. In the fiducial model that we consider, a slowly rotating black hole must scalarize if its size is comparable to the new length scale $\lambda$ that the theory introduces, although rapidly rotating black holes of any mass are effectively indistinguishable from their counterparts in general relativity. With this in mind, we use the gravitational-wave event GW190814$\,\unicode{x2014}\,$whose primary black hole has a spin that is bounded to be small, and whose signal shows no evidence of a scalarized primary$\,\unicode{x2014}\,$to rule out a narrow region of the parameter space. In particular, we find that values of ${\lambda \in [56, 96]~M_\odot}$ are strongly disfavored with a Bayes factor of $0.1$ or less. We also include a second event, GW151226, in our analysis to illustrate what information can be extracted when the spins of both components are poorly measured., Comment: Version accepted for publication. 8 pages + references, 4 figures. Some minor clarifications added relative to v1

To what extent are all astrophysical, dark, compact objects both black holes (BHs) and described by the Kerr geometry? We embark on the exercise of defying the universality of this remarkable idea, often called the "Kerr hypothesis". After establishing its rationale and timeliness, we define a minimal set of reasonability criteria for alternative models of dark compact objects. Then, as proof of principle, we discuss concrete, dynamically robust non-Kerr BHs and horizonless imitators, that 1) pass the basic theoretical, and in particular dynamical, tests, 2) match (some of the) state of the art astrophysical observables and 3) only emerge at some (macroscopic) scales. These examples illustrate how the universality (at all macroscopic scales) of the Kerr hypothesis can be challenged., Comment: 16 pages, 3 figures, contribution to the forthcoming book "Modified and Quantum Gravity - From theory to experimental searches on all scales"

In the presence of certain non-minimal couplings between a scalar field and the Gauss-Bonnet curvature invariant, Kerr black holes can scalarize, as long as they are spinning fast enough. This provides a distinctive violation of the Kerr hypothesis, occurring only for some high spin range. In this paper we assess if strong magnetic fields, that may exist in the vicinity of astrophysical black holes, could facilitate this distinctive effect, by bringing down the spin threshold for scalarization. This inquiry is motivated by the fact that self-gravitating magnetic fields, by themselves, can also promote "spin-induced" scalarization. Nonetheless, we show that in the \textit{vicinity of the horizon} the effect of the magnetic field $B$ on a black hole of mass $M$, up to $BM\lesssim 1$, works \textit{against} spin-induced scalarization, requiring a larger dimensionless spin $j$ from the black hole. A geometric interpretation for this result is suggested, in terms of the effects of rotation $vs.$ magnetic fields on the horizon geometry., Comment: 8 pages, 2 figures

Current gravitational wave (GW) detections rely on the existence of libraries of theoretical waveforms. Consequently, finding new physics with GWs requires libraries of non-standard models, which are computationally demanding. We discuss how deep learning frameworks can be used to generate new waveforms "learned" from a simulation dataset obtained, say, from numerical relativity simulations. Concretely, we use the WaveGAN architecture of a generative adversarial network (GAN). As a proof of concept we provide this neural network (NN) with a sample of ($>500$) waveforms from the collisions of exotic compact objects (Proca stars), obtained from numerical relativity simulations. Dividing the sample into a training and a validation set, we show that after a sufficiently large number of training epochs the NN can produce from 12\% to 25\% of the synthetic waveforms with an overlapping match of at least 95\% with the ones from the validation set. We also demonstrate that a NN can be used to predict the overlapping match score, with 90\% of accuracy, of new synthetic samples. These are encouraging results for using GANs for data augmentation and interpolation in the context of GWs, to cover the full parameter space of, say, exotic compact binaries, without the need of intensive numerical relativity simulations., Comment: 20 pages, 13 figures, 4 tables

Fully non-linear numerical evolutions of the Einstein-(multi)--Klein-Gordon equations are performed to study head-on collisions of $\ell$-boson stars. Despite being spherically symmetric, $\ell$-boson stars have a (hidden) frame of reference, used in defining their individual multipolar fields. To assess the impact of their relative orientation, we perform simulations with different angles between the axes of the two colliding stars. Additionally, two scenarios are considered for the colliding stars: that they are composites of either the same or different scalar fields. Despite some model-specific behaviours, the simulations generically indicate that: 1) the collision of two sufficiently (and equally) massive stars leads to black hole formation; 2) below a certain mass threshold the end result of the evolution is a bound state of the composite scalar fields, that neither disperses nor collapses into a black hole within the simulation time; 3) this end product (generically) deviates from spherical symmetry and the equipartition of the number of bosonic particles between the different scalar fields composing the initial boson stars is lost, albeit not dramatically. This last observation indicates, albeit without being conclusive, that the end result of these collisions belongs to the previously reported larger family of equilibrium multi-field boson stars, generically non-spherical, and of which $\ell$-boson stars are a symmetry enhanced point. We also extract and discuss the waveforms from the collisions studied., Comment: 23 pages. Updated version with broader discussion, typos corrected and new analyses and figures