Dynamics Of Binary Black Holes Research Articles

The effect of interstellar medium on LVK’s black holes

Abstract Gravitational radiation alone is not efficient in hardening the orbit of a wide binary black hole (BBH). By employing a toy model for the interstellar medium (ISM) surrounding BBHs, here we discuss the effect of this baryonic medium on BBH dynamics. Depending on the BBH’s mass, we show that a binary surrounded by an isotropic cold neutral medium (i.e., an asymptotic temperature T∞ ≈ 100 K) with a time-averaged particle density of $\langle n_H \rangle = \mathcal {O}(1)$ cm−3 can play a significant role in hardening the binary orbit over a $\mathcal {O}(10^9)$ yr time scale. Additionally, this causes the black hole’s mass to grow at a rate ∝m2. We thus discuss the impact of the ISM on the LIGO-Virgo-KAGRA (LVK) observables and quantify the properties of the ISM under which the latter could act as an additional important pathway for driving a subset of LVK’s BBH mergers.

Universality in binary black hole dynamics: An integrability conjecture

The waveform of a binary black hole coalescence appears to be both simple and universal. In this essay we argue that the dynamics should admit a separation into “fast and slow” degrees of freedom, such that the latter are described by an integrable system of equations, accounting for the simplicity and universality of the waveform. Given that Painlevé transcendents are a smoking gun of integrable structures, we propose the Painlevé-II transcendent as the key structural element threading a hierarchy of asymptotic models aiming at capturing different (effective) layers in the dynamics. Ward’s conjecture relating integrable and (anti-)self-dual solutions can provide the avenue to encode background binary black hole data in (nonlocal) twistor structures.

Dynamics of binary black holes in young star clusters: the impact of cluster mass and long-term evolution

ABSTRACT Dynamical interactions in dense star clusters are considered one of the most effective formation channels of binary black holes (BBHs). Here, we present direct N-body simulations of two different star cluster families: low-mass (∼500–800 M⊙) and relatively high-mass star clusters (≥5000 M⊙). We show that the formation channels of BBHs in low- and high-mass star clusters are extremely different and lead to two completely distinct populations of BBH mergers. Low-mass clusters host mainly low-mass BBHs born from binary evolution, while BBHs in high-mass clusters are relatively massive (chirp mass up to ∼100 M⊙) and driven by dynamical exchanges. Tidal disruption dramatically quenches the formation and dynamical evolution of BBHs in low-mass clusters on a very short time-scale (≲100 Myr), while BBHs in high-mass clusters undergo effective dynamical hardening until the end of our simulations (1.5 Gyr). In high-mass clusters, we find that 8 per cent of BBHs have primary mass in the pair-instability mass gap at metallicity Z = 0.002, all of them born via stellar collisions, while only one BBH with primary mass in the mass gap forms in low-mass clusters. These differences are crucial for the interpretation of the formation channels of gravitational-wave sources.

Bounding dark charges on binary black holes using gravitational waves

In models of minicharged dark matter associated with a hidden $U(1)$ symmetry, astrophysical black holes may acquire a "dark" charge, in such a way that the inspiral dynamics of binary black holes can be formally described by an Einstein-Maxwell theory. Charges enter the gravitational wave signal predominantly through a dipole term, but their effect is known to effectively first post-Newtonian order in the phase, which enables measuring the size of the charge-to-mass ratios, $|q_i/m_i|$, $i = 1,2$, of the individual black holes in a binary. We set up a Bayesian analysis to discover, or constrain, dark charges on binary black holes. After testing our framework in simulations, we apply it to selected binary black hole signals from the second Gravitational Wave Transient Catalog (GWTC-2), namely those with low masses so that most of the signal-to-noise ratio is in the inspiral regime. We find no evidence for charges on the black holes, and place typical 1-$\sigma$ bounds on the charge-to-mass ratios of $|q_i/m_i| \lesssim 0.2 - 0.3$.

Dynamics of binary black holes in low-mass young star clusters

ABSTRACT Young star clusters are dynamically active stellar systems and are a common birthplace for massive stars. Low-mass star clusters (∼300–103 M⊙) are more numerous than massive systems and are characterized by a two-body relaxation time-scale of a few Myr: the most massive stars sink to the cluster core and dynamically interact with each other even before they give birth to compact objects. Here, we explore the properties of black holes (BHs) and binary black holes (BBHs) formed in low-mass young star clusters, by means of a suite of 105 direct N-body simulations with a high original binary fraction (100 per cent for stars with mass >5 M⊙). Most BHs are ejected in the first ∼20 Myr by dynamical interactions. Dynamical exchanges are the main formation channel of BBHs, accounting for ∼40–80 per cent of all the systems. Most BBH mergers in low-mass young star clusters involve primary BHs with mass <40 M⊙ and low-mass ratios are extremely more common than in the field. Comparing our data with those of more massive star clusters (103 − 3 × 104 M⊙), we find a strong dependence of the percentage of exchanged BBHs on the mass of the host star cluster. In contrast, our results show just a mild correlation between the mass of the host star cluster and the efficiency of BBH mergers.

Integrability of eccentric, spinning black hole binaries up to second post-Newtonian order

Accurate and efficient modeling of the dynamics of binary black holes (BBHs) is crucial to their detection and parameter estimation through gravitational waves, both with LIGO/Virgo and LISA. General BBH configurations will have misaligned spins and eccentric orbits, eccentricity being particularly relevant at early times. Modeling these systems is both analytically and numerically challenging. Even though the 1.5 post-Newtonian (PN) order is Liouville integrable, numerical work has demonstrated chaos at 2PN order, which impedes the existence of an analytic solution. In this article we revisit integrability at both 1.5PN and 2PN orders. At 1.5PN, we construct four (out of five) action integrals. At 2PN, we show that the system is indeed integrable - but in a perturbative sense - by explicitly constructing five mutually-commuting constants of motion. Because of the KAM theorem, this is consistent with the past numerical demonstration of chaos. Our method extends to higher PN orders, opening the door for a fully analytical solution to the generic eccentric, spinning BBH problem.

Comparison of post-Newtonian mode amplitudes with numerical relativity simulations of binary black holes

Gravitational waves from the coalescence of two black holes carry the signature of the strong field dynamics of binary black holes. In this work we have used numerical relativity simulations and post-Newtonian theory to investigate this dynamics. Post-Newtonian theory is a low-velocity expansion that assumes the companion bodies to be point-particles, while numerical relativity treats black holes as extended objects with horizons and fully captures their dynamics. There is a priori no reason for the waveforms computed using these disparate methods to agree with each other, especially at late times when the black holes move close to the speed of light. We find, remarkably, that the leading order amplitudes in post-Newtonian theory agree well with the full general relativity solution for a large set of spherical harmonic modes, even in the most dynamical part of the binary evolution, with only some modes showing distinctly different behavior than that found by numerical relativity simulations. In particular, modes with spherical harmonic indices as well as are least modified from their dominant post-Newtonian behavior. Understanding the nature of these modes in terms of the post-Newtonian description will aid in formulating better models of the emitted waveforms in the strong field regime of the dynamics.

Revisiting the second post-Minkowskian eikonal and the dynamics of binary black holes

In this paper we study the two-body gravitational scattering of massive scalars with different masses in general spacetime dimensions. We focus on the Regge limit (eikonal regime) of the resulting scattering amplitudes and discuss how to extract the classical information representing the scattering of two black holes. We derive the leading eikonal and explicitly show the resummation of the first leading energy contribution up to second order in Newton's gravitational constant. We also calculate the subleading eikonal showing that in general spacetime dimensions it receives a non-trivial contribution from the box integral. From the eikonal we extract the two-body classical scattering angle between the two black holes up to the second post-Minkowskian order (2PM). Taking various probe-limits of the two-body scattering angles we are able to show agreement between our results and various results in the literature. We highlight that the box integral also has a log-divergent (in energy) contribution at subsubleading order which violates perturbative unitarity in the ultra-relativistic limit. We expect this term to play a role in the calculation of the eikonal at the 3PM order.

Spinning-black-hole scattering and the test-black-hole limit at second post-Minkowskian order

Recently, the gravitational scattering of two black holes (BHs) treated at the leading order in the weak-field, or post-Minkowskian (PM), approximation to General Relativity has been shown to map bijectively onto a simpler effectively one-body process: the scattering of a test BH in a stationary BH spacetime. Here, for BH spins aligned with the orbital angular momentum, we propose a simple extension of that mapping to 2PM order. We provide evidence for the validity and utility of this 2PM mapping by demonstrating its compatibility with all known analytical results for the conservative local-in-time dynamics of binary BHs in the post-Newtonian (weak-field and slow-motion) approximation and, separately, in the test-BH limit. Our result could be employed in the construction of improved effective-one-body models for the conservative dynamics of inspiraling spinning binary BHs.

Spin-multipole effects in binary black holes and the test-body limit

We discuss the effects of the black holes' spin-multipole structure in the orbital dynamics of binary black holes according to general relativity, focusing on the leading-post-Newtonian-order couplings at each order in an expansion in the black holes' spins. We first review previous widely confirmed results up through fourth order in spin, observe suggestive patterns therein, and discuss how the results can be extrapolated to all orders in spin with minimal information from the test-body limit. We then justify this extrapolation by providing a complete derivation within the post-Newtonian framework of a canonical Hamiltonian for a binary black hole, for generic orbits and spin orientations, which encompasses the leading post-Newtonian orders at all orders in spin. At the considered orders, the results reveal a precise equivalence between arbitrary-mass-ratio two-spinning-black-hole dynamics and the motion of a test black hole in a Kerr spacetime, as well as an intriguing relationship to geodesic motion in a Kerr spacetime.

Apparent ambiguities in the post-Newtonian expansion for binary systems

We discuss the source of the apparent ambiguities arising in the calculation of the dynamics of binary black holes within the Post-Newtonian framework. Divergences appear in both the near and far zone calculations, and may be of either ultraviolet (UV) or infrared (IR) nature. The effective field theory (EFT) formalism elucidates the origin of the singularities which may introduce apparent ambiguities. In particular, the only (physical) 'ambiguity parameters' that necessitate a matching calculation correspond to unknown finite size effects, which first appear at fifth Post-Newtonian (5PN) order for non-spinning bodies. We demonstrate that the ambiguities linked to IR divergences in the near zone, that plague the recent derivations of the binding energy at 4PN order, both in the ADM and 'Fokker-action' approach, can be resolved by implementing the so-called 'zero-bin' subtraction in the EFT framework. The procedure yields ambiguity-free results without the need of additional information beyond the PN expansion.

Testing general relativity using golden black-hole binaries

The coalescences of stellar-mass black-hole binaries through their inspiral, merger, and ringdown are among the most promising sources for ground-based gravitational-wave (GW) detectors. If a GW signal is observed with sufficient signal-to-noise ratio, the masses and spins of the black holes can be estimated from just the inspiral part of the signal. Using these estimates of the initial parameters of the binary, the mass and spin of the final black hole can be uniquely predicted making use of general-relativistic numerical simulations. In addition, the mass and spin of the final black hole can be independently estimated from the merger--ringdown part of the signal. If the binary black hole dynamics is correctly described by general relativity (GR), these independent estimates have to be consistent with each other. We present a Bayesian implementation of such a test of general relativity, which allows us to combine the constraints from multiple observations. Using kludge modified GR waveforms, we demonstrate that this test can detect sufficiently large deviations from GR, and outline the expected constraints from upcoming GW observations using the second-generation of ground-based GW detectors.

Merger states and final states of black hole coalescences: Comparing effective-one-body and numerical-relativity

We study to what extent the effective-one-body description of the dynamical state of a nonspinning, coalescing binary black hole (considered either at merger, or after ringdown) agrees with numerical relativity results. This comparison uses estimates of the integrated losses of energy and angular momentum during ringdown, inferred from recent numerical-relativity data. We find that the values, predicted by the effective-one-body formalism, of the energy and angular momentum of the system agree at the per mil level with their numerical-relativity counterparts, both at merger and in the final state. This gives a new confirmation of the ability of effective-one-body theory to accurately describe the dynamics of binary black holes even in the strong-gravitational-field regime. Our work also provides predictions (and analytical fits) for the final mass and the final spin of coalescing black holes for all mass ratios

Binary black hole mergers inf(R)theory

In the near future, gravitational wave detection is set to become an important observational tool for astrophysics. It will provide us with an excellent means to distinguish different gravitational theories. In effective form, many gravitational theories can be cast into an f(R) theory. In this article, we study the dynamics and gravitational waveform of an equal-mass binary black hole system in f(R) theory. We reduce the equations of motion in f(R) theory to the Einstein-Klein-Gordon coupled equations. In this form, it is straightforward to modify our existing numerical relativistic codes to simulate binary black hole mergers in f(R) theory. We considered binary black holes surrounded by a shell of scalar field. We solve the initial data numerically using the Olliptic code. The evolution part is calculated using the extended AMSSNCKU code. Both codes were updated and tested to solve the problem of binary black holes in f(R) theory. Our results show that the binary black hole dynamics in f(R) theory is more complex than in general relativity. In particular, the trajectory and merger time are strongly affected. Via the gravitational wave, it is possible to constrain the quadratic part parameter of f(R) theory in the range |a_2|<10^{11} m^2 . In principle, a gravitational wave detector can distinguish between a merger of binary black hole in f(R) theory and the respective merger in general relativity. Moreover, it is possible to use gravitational wave detection to distinguish between f(R) theory and a non self-interacting scalar field model in general relativity.

Simulating the dynamics of binary black holes in nuclear gaseous discs

AbstractWe study the pairing of massive black holes embedded in a massive circum–nuclear, rotationally supported disc, until they form a close binary. Using high resolution SPH simulations, we follow the black hole dynamics, and in particular the eccentricity evolution, as a function of the composition in stars and gas of the disc. Binary–disc interaction always leads to orbital decay and, in case of co–rotating black holes, to orbit circularization. We present also a higher resolution simulation performed using the particle–splitting technique showing that the binary orbital decay is efficient down to a separation of ~ 0.1 pc, comparable to our new resolution limit. We detail the gaseous mass profile bound to each black hole. Double nuclear activity is expected to occur on an estimated timescale of ≲ 10 Myrs.

Dynamics of precessing binary black holes using the post-Newtonian approximation

We investigate the (conservative) dynamics of binary black holes using the Hamiltonian formulation of the post-Newtonian (PN) equations of motion. The Hamiltonian we use includes spin-orbit coupling, spin-spin coupling, and mass monopole/spin-induced quadrupole interaction terms. In the case of both quasi-circular and eccentric orbits, we search for the presence of chaos (using the method of Lyapunov exponents) for a large variety of initial conditions. For quasi-circular orbits, we find no chaotic behavior for black holes with total mass 10 - 40 solar masses when initially at a separation corresponding to a Newtonian gravitational-wave frequency less than 150 Hz. Only for rather small initial radial distances, for which spin-spin induced oscillations in the radial separation are rather important, do we find chaotic solutions, and even then they are rare. Moreover, these chaotic quasi-circular orbits are of questionable astrophysical significance, since they originate from direct parametrization of the equations of motion rather than from widely separated binaries evolving to small separations under gravitational radiation reaction. In the case of highly eccentric orbits, which for ground-based interferometers are not astrophysically favored, we again find chaotic solutions, but only at pericenters so small that higher order PN corrections, especially higher spin PN corrections, should also be taken into account.

Binary black-hole dynamics at the third-and-a-half post-Newtonian order in the ADM formalism

We specialize the radiation-reaction part of the Arnowitt-Deser-Misner (ADM) Hamiltonian for many non-spinning point-like bodies, calculated by Jaranowski and Schaefer [1], to third-and-a-half post-Newtonian approximation to general relativity, to binary systems. This Hamiltonian is used for the computation of the instantaneous gravitational energy loss of a binary to 1PN reactive order. We also derive the equations of motion, which include PN reactive terms via Hamiltonian and Euler-Lagrangian approaches. The results are consistent with the expressions for reactive acceleration provided by Iyer-Will formalism in Ref. [2] in a general class of gauges.

The binary black-hole dynamics at the third post-Newtonian order in the orbital motion

The orbital dynamics of the binary point-mass systems is ambiguous at the third post-Newtonian order of approximation. The static ambiguity is known to be related to the difference between the Brill-Lindquist solution and the Misner-Lindquist solution of the time-symmetric conformally flat initial value problem for binary black holes.The kinetic ambiguity is noticed to violate in general the standard relation between the center-of-mass velocity and the total linear momentum as demanded by global Lorentz invariance.