ABSTRACT Many astrophysical environments, from star clusters and globular clusters to the discs of active galactic nuclei, are characterized by frequent interactions between stars and the compact objects that they leave behind. Here, using a suite of 3D hydrodynamics simulations, we explore the outcome of close interactions between $$1\, \mathrm{M}_{\odot }$$ stars and binary black holes (BBHs) in the gravitational wave regime, resulting in a tidal disruption event (TDE) or a pure scattering, focusing on the accretion rates, the back reaction on the BH binary orbital parameters, and the increase in the binary BH effective spin. We find that TDEs can make a significant impact on the binary orbit, which is often different from that of a pure scattering. Binaries experiencing a prograde (retrograde) TDE tend to be widened (hardened) by up to $$\simeq 20{{\ \rm per\ cent}}$$. Initially circular binaries become more eccentric by $$\lesssim 10{{\ \rm per\ cent}}$$ by a prograde or retrograde TDE, whereas the eccentricity of initially eccentric binaries increases (decreases) by a retrograde (prograde) TDE by $$\lesssim 5{{\ \rm per\ cent}}$$. Overall, a single TDE can generally result in changes of the gravitational-wave-driven merger time-scale by order unity. The accretion rates of both black holes are very highly super-Eddington, showing modulations (preferentially for retrograde TDEs) on a time-scale of the orbital period, which can be a characteristic feature of BBH-driven TDEs. Prograde TDEs result in the effective spin parameter χ to vary by ≲0.02, while χ ≳ −0.005 for retrograde TDEs.
more »
« less
Impact of eccentricity and mean anomaly in numerical relativity mergers
Abstract Accurate modelling of black hole (BH) binaries is critical to achieve the science goals of gravitational-wave detectors. Modelling such configurations relies strongly on calibration to numerical-relativity (NR) simulations. Binaries on quasi-circular orbits have been widely explored in NR, however, coverage of the broader nine-dimensional parameter space, including orbital eccentricity, remains sparse. This article develops a new procedure to control orbital eccentricity of binary BH (BBH) simulations that enables choosing initial data parameters with precise control over eccentricity and mean anomaly of the subsequent evolution, as well as the coalescence time. We then calculate several sequences of NR simulations that nearly uniformly cover the two-dimensional eccentricity–mean anomaly space for equal mass, non-spinning BBHs. We demonstrate that, for fixed eccentricity, many quantities related to the merger dynamics of BBHs show an oscillatory dependence on mean anomaly. The amplitude of these oscillations scales nearly linearly with the eccentricity of the system. We find that for the eccentricities explored in this work, the magnitude of deviations in various quantities such as the merger amplitude and peak luminosity can approach ∼5% of their quasi-circular value. We use our findings to explain eccentric phenomena reported in other studies. We also show that methods for estimating the remnant mass employed in the effective-one-body approach exhibit similar deviations, roughly matching the amplitude of the oscillations we find in NR simulations. This work is an important step towards a complete description of eccentric BBH mergers, and demonstrates the importance of considering the entire two-dimensional parameter subspace related to eccentricity.
more »
« less
- PAR ID:
- 10632257
- Publisher / Repository:
- Classical and Quantum Gravity
- Date Published:
- Journal Name:
- Classical and Quantum Gravity
- Volume:
- 42
- Issue:
- 13
- ISSN:
- 0264-9381
- Page Range / eLocation ID:
- 135011
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
The measurement of orbital eccentricity in gravitational-wave (GW) signals will provide unique insights into the astrophysical origin of binary systems, while ignoring eccentricity in waveform models could introduce significant biases in parameter estimation and tests of general relativity. Upcoming LIGO-Virgo-KAGRA observing runs are expected to detect a subpopulation of eccentric signals, making it vital to develop accurate waveform models for eccentric orbits. Here, employing recent analytical results through the third post-Newtonian order, we develop v5: a new time-domain, effective-one-body, multipolar waveform model for eccentric binary black holes with spins aligned (or antialigned) with the orbital angular momentum. Besides the dominant (2, 2) mode, the model includes the (2, 1), (3, 3), (3, 2), (4, 4), and (4, 3) modes. We validate the model’s accuracy by computing its unfaithfulness against 99 (28 public and 71 private) eccentric numerical-relativity (NR) simulations, produced by the Simulating eXtreme Spacetimes Collaboration. Importantly, for NR waveforms with initial GW eccentricities below 0.5, the maximum (2, 2)-mode unfaithfulness across the total mass range is consistently below or close to 1%, with a median value of , reflecting an accuracy improvement of approximately an order of magnitude compared to the previous-generation v4 and the state-of-the-art esumalí eccentric model. In the quasi-circular-orbit limit, v5 is in excellent agreement with the highly accurate v5 model. The accuracy, robustness, and speed of v5 make it suitable for data analysis and astrophysical studies. We demonstrate this by performing a set of recovery studies of synthetic NR-signal injections, and parameter-estimation analyses of the events GW150914 and GW190521, which we find to have no eccentricity signatures.more » « less
-
ABSTRACT Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code AREPO. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary’s semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly one out of four encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary’s merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disc with accretion rates exceeding the Eddington limit. If the debris disc cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger.more » « less
-
Abstract Despite the growing number of binary black hole coalescences confidently observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include the effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that have already been identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total source-frame massM> 70M⊙) binaries covering eccentricities up to 0.3 at 15 Hz emitted gravitational-wave frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place a conservative upper limit for the merger rate density of high-mass binaries with eccentricities 0 <e≤ 0.3 at 16.9 Gpc−3yr−1at the 90% confidence level.more » « less
-
Abstract Gravitational-wave detectors are now making it possible to investigate how the merger rate of binary black holes (BBHs) evolves with redshift. In this study, we examine whether the BBH merger rate of isolated binaries deviates from a scaled star formation rate density (SFRD)—a frequently used model in state-of-the-art research. To address this question, we conduct population synthesis simulations using COMPAS with a grid of stellar evolution models, calculate their cosmological merger rates, and compare them to a scaled SFRD. We find that our simulated rates deviate by factors up to 3.5 atz∼ 0 and 5 atz∼ 9 due to two main phenomena: (i) the formation efficiency of BBHs is an order of magnitude higher at low metallicities than at solar metallicity, and (ii) BBHs experience a wide range of delays (from a few megayears to many gigayears) between formation and merger. The deviations are similar when comparing to a delayed SFRD, and even larger (up to ∼10×) when comparing to SFRD-based models scaled to the local merger rate. Interestingly, our simulations find that the BBH delay time distribution is redshift dependent, increasing the complexity of the redshift distribution of mergers. We find similar results for simulated merger rates of black hole–neutron stars (BHNSs) and binary neutron stars (BNSs). We conclude that the rate of BBH, BHNS, and BNS mergers from the isolated channel can significantly deviate from a scaled SFRD, and that future measurements of the merger rate will provide insights into the formation pathways of gravitational-wave sources.more » « less
An official website of the United States government

