Equal-mass Black Holes Research Articles

Self-interacting scalar dark matter around binary black holes

Gravitational waves can provide crucial insights about the environments in which black holes live. In this work, we use numerical relativity simulations to study the behavior of self-interacting scalar (wavelike) dark matter clouds accreting onto isolated and binary black holes. We find that repulsive self-interactions smoothen the “spike” of an isolated black hole and saturate the density. Attractive self-interactions enhance the growth and result in more cuspy profiles, but can become unstable and undergo explosions akin to the superradiant bosenova that reduce the local cloud density. We quantify the impact of self-interactions on an equal-mass black hole merger by computing the dephasing of the gravitational-wave signal for a range of couplings. We find that repulsive self-interactions saturate the density of the cloud, thereby reducing the dephasing. For attractive self-interactions, the dephasing may be larger, but if these interactions dominate prior to the merger, the dark matter can undergo bosenova during the inspiral phase, disrupting the cloud and subsequently reducing the dephasing. Published by the American Physical Society 2024

Post-Minkowskian theory meets the spinning effective-one-body approach for two-body scattering

Effective-one-body (EOB) waveforms employed by the LIGO-Virgo-KAGRA Collaboration have primarily been developed by resumming the post-Newtonian expansion of the relativistic two-body problem. Given the recent significant advancements in post-Minkowskian (PM) theory and gravitational self-force formalism, there is considerable interest in creating waveform models that integrate information from various perturbative methods in innovative ways. This becomes particularly crucial when tackling the accuracy challenge posed by upcoming ground-based detectors (such as the Einstein Telescope and Cosmic Explorer) and space-based detectors (such as LISA, TianQin, or Taiji) expected to operate in the next decade. In this context, we present the derivation of the first spinning EOB (SEOB) Hamiltonian that incorporates PM results up to three-loop order: the SEOB-PM model. The model accounts for the complete hyperbolic motion, encompassing nonlocal-in-time tails. To evaluate its accuracy, we compare its predictions for the conservative scattering angle, augmented with dissipative contributions, against numerical-relativity data of nonspinning and spinning equal-mass black holes. We observe very good agreement, comparable, and in some cases slightly better to the recently proposed wEOB-potential model, of which the SEOB-PM model is a resummation around the probe limit. Indeed, in the probe limit, the SEOB-PM Hamiltonian and scattering angles reduce to the one of a test mass in Kerr spacetime. Once complemented with nonlocal-in-time contributions for bound orbits, the SEOB-PM Hamiltonian can be utilized to generate waveform models for spinning black holes on quasicircular orbits. Published by the American Physical Society 2024

Three Black Holes that are Coming Closer Together

We current the first completely relativistic lasting numerical progresses of three equal-mass black holes in a system containing of a third black hole in a close orbit about a black-hole dual. We discover that these close-three-black-hole systems have very dissimilar merger dynamics from black-hole duals. In specific, we see complex trajectories, a redeployment of energy that can communicate substantial kicks to one of the holes, characteristic waveforms, and suppression of the emitted gravitational radiation. We change two such configurations and discover very different actions. In one conformation the dual is quickly troubled and the separate holes follow complicated trajectories and coalesce with the third hole in fast succession, while in the other, the dual finishes a half-orbit before the initial merger of one of the members with the third black hole, and the subsequent two-black-hole system forms a highly indirect, well separated dual that shows no significant in spiral for (at least) the first t ∼ 1000 M of evolution.

Black hole merger simulations in wave dark matter environments

The interaction of binary black hole mergers with their environments can be studied using numerical relativity simulations. These start only a short finite time before merger, at which point appropriate initial conditions must be imposed. A key task is therefore to identify the configuration that is appropriate for the binary and its environment at this stage of the evolution. In this work we study the behavior of wave dark matter around equal mass black hole binaries, finding that there is a preferred, quasistationary profile that persists and grows over multiple orbits, in contrast to heavier mass dark matter where any overdensity tends to be dispersed by the binary motion. While different initial configurations converge to the preferred quasistationary one after several orbits, unwanted transient oscillations are generated in the process, which may have an impact on the signal in short simulation runs. We also point out that naively superimposing the matter onto a circular binary results in artificially eccentric orbits due to the matter backreaction, which is an effect of the initial conditions and not a signature of dark matter. We discuss the further work required so that comparison of waveforms obtained with environments to vacuum cases can be done in a meaningful way.

Well-Posedness of the Four-Derivative Scalar-Tensor Theory of Gravity in Singularity Avoiding Coordinates.

We show that the most general scalar-tensor theory of gravity up to four derivatives in 3+1 dimensions is well-posed in a modified version of the CCZ4 formulation of the Einstein equations in singularity-avoiding coordinates. We demonstrate the robustness of our new formulation in practice by studying equal mass black hole binary mergers for different values of the coupling constants. Although our analysis of well-posedness is restricted to cases in which the couplings are small, we find that in simulations we are able to push the couplings to larger values, so that a certain weak coupling condition is order one, without instabilities developing. Our Letter provides the means for such simulations to be undertaken by the many numerical relativity codes that rely on the moving puncture gauge to evolve black hole singularities.

Electromagnetic Signatures from Supermassive Binary Black Holes Approaching Merger

Abstract We present fully relativistic predictions for the electromagnetic emission produced by accretion disks surrounding spinning and nonspinning supermassive binary black holes on the verge of merging. We use the code Bothros to post-process data from 3D general relativistic magnetohydrodynamic simulations via ray-tracing calculations. These simulations model the dynamics of a circumbinary disk and the mini-disks that form around two equal-mass black holes orbiting each other at an initial separation of 20 gravitational radii, and evolve the system for more than 10 orbits in the inspiral regime. We model the emission as the sum of thermal blackbody radiation emitted by an optically thick accretion disk and a power-law spectrum extending to hard X-rays emitted by a hot optically thin corona. We generate time-dependent spectra, images, and light curves at various frequencies to investigate intrinsic periodic signals in the emission, as well as the effects of the black hole spin. We find that prograde black hole spin makes mini-disks brighter since the smaller innermost stable circular orbit angular momentum demands more dissipation before matter plunges to the horizon. However, compared to mini-disks in larger separation binaries with spinning black holes, our mini-disks are less luminous: unlike those systems, their mass accretion rate is lower than in the circumbinary disk, and they radiate with lower efficiency because their inflow times are shorter. Compared to a single black hole system matched in mass and accretion rate, these binaries have spectra noticeably weaker and softer in the UV. Finally, we discuss the implications of our findings for the potential observability of these systems.

Geometric horizons in binary black hole mergers

We numerically study the algebraic properties of the Weyl tensor through the merger of two non-spinning black holes (BHs). We are particularly interested in the conjecture that for such a vacuum spacetime, which is zeroth-order algebraically general, a geometric horizon (GH), on which the spacetime is algebraically special and which is identified by the vanishing of a complex scalar invariant (), characterizes a smooth foliation independent surface (horizon) associated with the BH. In the first simulation we investigate the level-0 sets of (since ) in the head-on collision of two unequal mass BHs. In the second simulation we shall investigate the level-ɛ sets of through a quasi-circular merger of two non-spinning, equal mass BHs. The numerical results, as displayed in the figures presented, provide evidence that a (unique) smooth GH can be identified throughout all stages of the binary BH merger.

On accretion discs formed in MHD simulations of black hole–neutron star mergers with accurate microphysics

ABSTRACT Remnant accretion discs formed in compact object mergers are an important ingredient in the understanding of electromagnetic afterglows of multimessenger gravitational-wave events. Due to magnetically and neutrino-driven winds, a significant fraction of the disc mass will eventually become unbound and undergo r-process nucleosynthesis. While this process has been studied in some detail, previous studies have typically used approximate initial conditions for the accretion discs, or started from purely hydrodynamical simulations. In this work, we analyse the properties of accretion discs formed from near equal-mass black hole–neutron star mergers simulated in general-relativistic magnetohydrodynamics in dynamical spacetimes with an accurate microphysical description. The post-merger systems were evolved until $120\, {\rm ms}$ for different finite-temperature equations of state and black hole spins. We present a detailed analysis of the fluid properties and of the magnetic-field topology. In particular, we provide analytic fits of the magnetic-field strength and specific entropy as a function of the rest-mass density, which can be used for the construction of equilibrium disc models. Finally, we evolve one of the systems for a total of $350\, \rm ms$ after merger and study the prospect for eventual jet launching. While our simulations do not reach this stage, we find clear evidence of continued funnel magnetization and clearing, a prerequisite for any jet-launching mechanism.

Stable circular orbits in higher-dimensional multi–black-hole spacetimes

We consider the dynamics of particles, particularly focusing on circular orbits in the higher-dimensional Majumdar-Papapetrou (MP) spacetimes with two equal mass black holes. It is widely known that in the 5D Schwarzschild-Tangherlini and Myers-Perry backgrounds, there are no stable circular orbits. In contrast, we show that in the 5D MP background, stable circular orbits can always exist when the separation of two black holes is large enough. More precisely, for a large separation, stable circular orbits exist from the vicinity of horizons to infinity; for a medium one, they appear only in a certain finite region bounded by the innermost stable circular orbit and the outermost stable circular orbit outside the horizons; for a small one, they do not appear at all. Moreover, we show that in MP spacetimes in more than 5D, they do not exist for any separations.

SAGE: finding IMBH in the black hole desert

SAGE (SagnAc interferometer for Gravitational wavE) is a project for a space observatory based on multiple 12-U CubeSats in geosynchronous equatorial orbit. The objective is a fast track mission which would fill the observational gap between LISA and ground based observatories. With albeit a lower sensitivity, it would allow early investigation of the nature and event rate of intermediate-mass black hole (IMBH) mergers, constraining our understanding of the universe formation by probing the building up of IMBH up to supermassive black holes (SMBH). Technically, the CubeSats would create a triangular Sagnac interferometer with 140.000 km roundtrip arm length, optimised to be sensitive to gravitational waves at frequencies between 10 mHz and 2 Hz. The nature of the Sagnac measurement makes it almost insensitive to position error, a feature enabling the use of spacecrafts in ballistic trajectories instead of perfect free fall. The light source and recombination units of the interferometer are based on compact fibered technologies without bulk optics. A peak sensitivity of 23 pm ()−1 is expected at 1 Hz assuming a 200 mW internal laser source and 10-centimeter diameter apertures. Because of the absence of a test mass, the main limitation would come from the non-gravitational forces applied on the spacecrafts. However, conditionally upon our ability to partially post-process the effect of solar wind and solar pressure, SAGE would allow detection of gravitational waves with strains as low as a few 10−19 within the 0.1 to 1 Hz range. Averaged over the entire sky, and including the antenna gain of the Sagnac interferometer, the SAGE observatory would sense equal mass black hole mergers in the 104 to 106 solar masses range up to a luminosity distance of 800 Mpc. Additionally, coalescence of stellar black holes (10 M) around SMBH (IMBH) forming extreme (intermediate) mass ratio inspirals could be detected within a sphere of radius 200 Mpc.

Geometric horizons in the Kastor-Traschen multi-black-hole solutions

We investigate the existence of invariantly defined quasi-local hypersurfaces in the Kastor-Traschen solution containing $N$ charge-equal-to-mass black holes. These hypersurfaces are characterized by the vanishing of particular curvature invariants, known as Cartan invariants, which are generated using the frame approach. The Cartan invariants of interest describe the expansion of the outgoing and ingoing null vectors belonging to the invariant null frame arising from the Cartan-Karlhede algorithm. We show that the evolution of the hypersurfaces surrounding the black holes depends on an upper-bound on the total mass for the case of two and three equal mass black holes. We discuss the results in the context of the geometric horizon conjectures.

Gravitational radiation driven capture in unequal mass black hole encounters

The gravitational radiation driven capture (GR capture) between unequal mass black holes without spins has been investigated with numerical relativistic simulations. We adopt the parabolic approximation which assumes that the gravitational wave radiation from a weakly hyperbolic orbit is the same as that from the parabolic orbit having the same pericenter distance. Using the radiated energies from the parabolic orbit simulations, we have obtained the critical impact parameter ($b_{\rm crit}$) for the GR capture for weakly hyperbolic orbit as a function of initial energy. The most energetic encounters occur around the boundary between the direct merging and the fly-by orbits, and can emit several percent of initial total ADM energy at the peak. When the total mass is fixed, energy and angular momentum radiated in the case of unequal mass black holes are smaller than those of equal mass black holes having the same initial orbital angular momentum for the fly-by orbits. We have compared our results with two different Post-Newtonian (PN) approximations, the exact parabolic orbit (EPO) and PN corrected orbit (PNCO). We find that the agreement between the EPO and the numerical relativity breaks down for very close encounters ($\it{e.g.}$, $b_{\rm crit} \lesssim 100$ M), and it becomes worse for higher mass ratios. For instance, the critical impact parameters can differ by more than $50\%$ from those obtained in EPO if the relative velocity at infinity $v_{\infty}$ is larger than 0.1 for the mass ratio of $m_{1}/m_{2}=16$. The PNCO gives more consistent results than EPO, but it also underestimates the critical impact parameter for the GR capture at $b_{\rm crit} \lesssim40$ M.

Puncture initial data for black-hole binaries with high spins and high boosts

We solve the Hamiltonian and momentum constraints of general relativity for two black holes with nearly extremal spins and relativistic boosts in the puncture formalism. We use a non-conformally-flat ansatz with an attenuated superposition of two Lorentz-boosted, conformally Kerr or conformally Schwarzschild 3-metrics and their corresponding extrinsic curvatures. We compare evolutions of these data with the standard Bowen-York conformally flat ansatz (technically limited to intrinsic spins $\chi=S/M^2_{\text{ADM}}=0.928$ and boosts $P/M_{\text{ADM}}=0.897$), finding, typically, an order of magnitude smaller burst of spurious radiation and agreement with inspiral and merger. As a first case study, we evolve two equal-mass black holes from rest with an initial separation of $d=12M$ and spins $\chi_i=S_i/m_i^2=0.99$, compute the waveforms produced by the collision, the energy and angular momentum radiated, and the recoil of the final remnant black hole. We find that the black-hole trajectories curve at close separations, leading to the radiation of angular momentum. We also study orbiting nonspinning and moderate-spin black-hole binaries and compare these with standard Bowen-York data. We find a substantial reduction in the nonphysical initial burst of radiation which leads to cleaner waveforms. Finally, we study the case of orbiting binary black-hole systems with spin magnitude $\chi_i=0.95$ in an aligned configuration and compare waveform and final remnant results with those of the SXS Collaboration, finding excellent agreement. This represents the first moving punctures evolution of orbiting and spinning black holes exceeding the Bowen-York limit. Finally, we study different choices of the initial lapse and lapse evolution equation in the moving punctures approach to improve the accuracy and efficiency of the simulations.

High energy collisions of black holes numerically revisited

We use fully nonlinear numerical relativity techniques to study high energy head-on collision of nonspinning, equal-mass black holes to estimate the maximum gravitational radiation emitted by these systems. Our simulations include improvements in the construction of initial data, subsequent full numerical evolutions, and the computation of waveforms at infinity. The new initial data significantly reduces the spurious radiation content, allowing for initial speeds much closer to the speed of light, i.e. $v\sim0.99c$. Using these new techniques, We estimate the maximum radiated energy from head-on collisions to be $E_{\text{max}}/M_{\text{ADM}}=0.13\pm0.01$. This value differs from the second-order perturbative $(0.164)$ and zero-frequency-limit $(0.17)$ analytic computations, but is close to those obtained by thermodynamic arguments $(0.134)$ and by previous numerical estimates $(0.14\pm0.03)$.

The detection of quasinormal mode witha/M= 0.95 would prove a sphere 99% soaking in the ergoregion of the Kerr space-time

Recent numerical relativity simulations of mergers of binary black holes suggest that the maximum final value of $a/M$ is $\sim 0.95$ for the coalescence of two equal mass black holes with aligned spins of the same magnitude $a/M=0.994$ which is close to the upper limit $a/M=0.998$ of accretion spin-up shown by Thorne. Using the WKB method, we suggest that if quasinormal modes with $a/M \sim 0.95$ are detected by the second generation gravitational wave detectors, we could confirm the strong gravity space-time based on Einstein's general relativity up to $1.33M$ which is only $\sim 1.014$ times the event horizon radius and within the ergoregion. One more message about black hole geometry is expected here. If the quasinormal mode is different from that of general relativity, we need to find the true theory of gravity which deviates from general relativity only near the black hole horizon.

Higher dimensional numerical relativity: Code comparison

The nonlinear behavior of higher dimensional black hole spacetimes is of interest in several contexts, ranging from an understanding of cosmic censorship to black hole production in high-energy collisions. However, nonlinear numerical evolutions of higher dimensional black hole spacetimes are tremendously complex, involving different diagnostic tools and "dimensional reduction methods". In this work we compare two different successful codes to evolve Einstein's equations in higher dimensions, and show that the results of such different procedures agree to numerical precision, when applied to the collision from rest of two equal-mass black holes. We calculate the total radiated energy to be E/M=9x10^{-4} in five dimensions and E/M=8.1x10^{-4} in six dimensions.

Binary black hole circular orbits computed with cocal

In this work we present our first results of binary black hole circular orbits using {\sc cocal}, the Compact Object CALculator. Using the 3+1 decomposition five equations are being solved under the assumptions of conformal flatness and maximal slicing. Excision is used and the appropriate apparent horizon boundary conditions are applied. The orbital velocity is determined by imposing a Schwarzschild behaviour at infinity. A sequence of equal mass black holes is obtained and its main physical characteristics are calculated.

Superkicks in ultrarelativistic encounters of spinning black holes

We study ultrarelativistic encounters of two spinning, equal-mass black holes through simulations in full numerical relativity. Two initial data sequences are studied in detail: one that leads to scattering and one that leads to a grazing collision and merger. In all cases, the initial black hole spins lie in the orbital plane, a configuration that leads to the so-called superkicks. In astrophysical, quasicircular inspirals, such kicks can be as large as $\ensuremath{\sim}3000\text{ }\text{ }\mathrm{km}/\mathrm{s}$; here, we find configurations that exceed $\ensuremath{\sim}15\text{ }000\text{ }\text{ }\mathrm{km}/\mathrm{s}$. We find that the maximum recoil is to a good approximation proportional to the total amount of energy radiated in gravitational waves, but largely independent of whether a merger occurs or not. This shows that the mechanism predominantly responsible for the superkick is not related to merger dynamics. Rather, a consistent explanation is that the ``bobbing'' motion of the orbit causes an asymmetric beaming of the radiation produced by the in-plane orbital motion of the binary, and the net asymmetry is balanced by a recoil. We use our results to formulate some conjectures on the ultimate kick achievable in any black hole encounter.

Simulating merging binary black holes with nearly extremal spins

Astrophysically realistic black holes may have spins that are nearly extremal (i.e., close to 1 in dimensionless units). Numerical simulations of binary black holes are important tools both for calibrating analytical templates for gravitational-wave detection and for exploring the nonlinear dynamics of curved spacetime. However, all previous simulations of binary-black-hole inspiral, merger, and ringdown have been limited by an apparently insurmountable barrier: the merging holes' spins could not exceed 0.93, which is still a long way from the maximum possible value in terms of the physical effects of the spin. In this paper, we surpass this limit for the first time, opening the way to explore numerically the behavior of merging, nearly extremal black holes. Specifically, using an improved initial-data method suitable for binary black holes with nearly extremal spins, we simulate the inspiral (through 12.5 orbits), merger and ringdown of two equal-mass black holes with equal spins of magnitude 0.95 antialigned with the orbital angular momentum.

Spinning compact binary inspiral. II. Conservative angular dynamics

We establish the evolution equations of the set of independent variables characterizing the 2PN rigorous conservative dynamics of a spinning compact binary, with the inclusion of the leading order spin-orbit, spin-spin and mass quadrupole - mass monopole effects, for generic (noncircular, nonspherical) orbits. More specifically, we give a closed system of first order ordinary differential equations for the orbital elements of the osculating ellipse and for the angles characterizing the spin orientations with respect to the osculating orbit. We also prove that (i) the relative angle of the spins stays constant for equal mass black holes, irrespective of their orientation, and (ii) the special configuration of equal mass black holes with equal, but antialigned spins, both laying in the plane of motion (leading to the largest recoil found in numerical simulations) is preserved at 2PN level of accuracy, with leading order spin-orbit, spin-spin and mass quadrupolar contributions included.