Accurate models for recoil velocity distribution in black hole mergers with comparable to extreme mass-ratios and their astrophysical implications
read the original abstract
Modeling the remnant recoil velocity (kick) distribution from binary black hole mergers is crucial for understanding hierarchical mergers in active galactic nuclei or globular clusters. Existing analytic models often show large discrepancies with numerical relativity (NR) data, while data-driven models are limited to mass ratios of q<=8 (aligned spins) and q<=4 (precessing spins) and break down when extrapolated outside their training ranges. Using ~5000 of NR simulations from the SXS and RIT catalogs up to q=128 and ~100 black hole perturbation theory simulations up to q=200, we present two classes of models: (i) gwModel_kick_q200 (gwModel_kick_q200_GPR), an analytic (Gaussian process regression) model for aligned-spin binaries. (ii) gwModel_kick_prec_flow, a normalizing-flow model for kick distribution from precessing binaries with isotropic spins. Our approach combines analytic insights from post-Newtonian theory with data-driven techniques to ensure correct limiting behavior and high accuracy across parameter space. Both gwModel_kick_q200 and gwModel_kick_prec_flow are valid from comparable to extreme mass ratios. Extensive validation shows all three models outperform existing ones within their respective domains. Finally, using both back-of-the-envelope estimates and 1404 detailed star cluster simulations incorporating our kick models, we find that the black hole retention probability in low mass globular clusters can vary noticeably when the gwModel_kick_prec_flow model is employed. The models are publicly available through the gwModels package.
This paper has not been read by Pith yet.
Forward citations
Cited by 5 Pith papers
-
Predicting intermediate-mass black hole formation in star clusters with machine learning
Machine learning regressors trained on Rapster simulations forecast that globular clusters rarely host black holes above 100 solar masses while a few nuclear star clusters may exceed this threshold.
-
gwBenchmarks: Stress-Testing LLM Agents on High-Precision Gravitational Wave Astronomy
LLM coding agents cannot reach the 10^{-4} relative accuracy required for gravitational wave modeling tasks and show systematic failures including metric misuse and result fabrication.
-
Milky Way Globular Clusters: Nurseries for Dynamically-Formed Binary Black Holes
A coupled galaxy-cluster model shows that Milky Way globular clusters, many from satellite galaxies, form dynamically interacting black hole binaries whose merger rates increase with redshift up to z=5.
-
Recoil kicks from binary black hole mergers in GWTC catalogs: implications for retention and hierarchical mergers
Recoil kicks are inferred for GWTC-4 binary black hole events with values up to nearly 1000 km/s for some, yielding retention probabilities of 1-5% in globular clusters and 70-100% in elliptical galaxies.
-
Remnant recoil and host environments of GWTC-4.0 binary black-hole mergers
Five O4 gravitational-wave events prefer dynamical cluster formation; typical recoil kicks eject remnants from globular clusters but allow possible retention in nuclear star clusters, disfavoring efficient hierarchica...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.