From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations

Fuente: arXiv
Saved in:
Bibliographic Details
Main Authors: Levina, Sasha, Broekgaarden, Floor, van Son, Lieke, Berti, Emanuele, Romagnolo, Amedeo, Pakmor, Ruediger, Lam, Ana
Format: Preprint
Published: 2026
Subjects:
Online Access:
Tags: Add Tag
No Tags, Be the first to tag this record!
_version_ 1866910003184533504
author Levina, Sasha
Broekgaarden, Floor
van Son, Lieke
Berti, Emanuele
Romagnolo, Amedeo
Pakmor, Ruediger
Lam, Ana
author_facet Levina, Sasha
Broekgaarden, Floor
van Son, Lieke
Berti, Emanuele
Romagnolo, Amedeo
Pakmor, Ruediger
Lam, Ana
contents Observations of binary black hole (BBH) mergers provide a unique window into the lives of massive stars across cosmic time. Connecting redshift-dependent merger properties to massive star progenitors requires accurate models of cosmic star formation and chemical enrichment histories. Analytical fits for the metallicity-specific cosmic star formation rate density S(Z, z) are commonly used as proxies for the complex underlying star formation history, yet they remain unconstrained. Using the IllustrisTNG cosmological simulations, we evaluate the accuracy of these analytical S(Z, z) prescriptions and assess how simulation resolution and volume affect the inferred S(Z, z). By coupling the simulated and analytical S(Z, z) to the population synthesis code COMPAS, we investigate the resulting BBH merger rates and mass distributions. We find that analytical S(Z, z) prescriptions can overestimate BBH merger rates at high redshift ($z \gtrsim 6$) by up to a factor of $10$-$10^4$, depending on cosmological simulation resolution, and can introduce spurious features in the BBH mass distribution. For example, they can produce an artificial feature near $8\,M_\odot$ in the primary mass distribution at $z \lesssim 2$, which is absent when using the full simulation-based S(Z, z), while simultaneously suppressing high-mass features. These discrepancies arise because simple analytical models fail to capture a high-metallicity bump and a more flattened low-metallicity tail in the simulated S(Z, z) metallicity distribution. Our results highlight the importance of accurate star formation histories for modeling BBH populations, demonstrate the limitation of widely used analytical S(Z, z) fits, and underscore the need for careful integration of cosmological simulations, analytical fits, and population synthesis when interpreting gravitational-wave observations.
format Preprint
id arxiv_https___arxiv_org_abs_2601_20202
institution arXiv
publishDate 2026
record_format arxiv
spellingShingle From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations
Levina, Sasha
Broekgaarden, Floor
van Son, Lieke
Berti, Emanuele
Romagnolo, Amedeo
Pakmor, Ruediger
Lam, Ana
High Energy Astrophysical Phenomena
Astrophysics of Galaxies
Solar and Stellar Astrophysics
Observations of binary black hole (BBH) mergers provide a unique window into the lives of massive stars across cosmic time. Connecting redshift-dependent merger properties to massive star progenitors requires accurate models of cosmic star formation and chemical enrichment histories. Analytical fits for the metallicity-specific cosmic star formation rate density S(Z, z) are commonly used as proxies for the complex underlying star formation history, yet they remain unconstrained. Using the IllustrisTNG cosmological simulations, we evaluate the accuracy of these analytical S(Z, z) prescriptions and assess how simulation resolution and volume affect the inferred S(Z, z). By coupling the simulated and analytical S(Z, z) to the population synthesis code COMPAS, we investigate the resulting BBH merger rates and mass distributions. We find that analytical S(Z, z) prescriptions can overestimate BBH merger rates at high redshift ($z \gtrsim 6$) by up to a factor of $10$-$10^4$, depending on cosmological simulation resolution, and can introduce spurious features in the BBH mass distribution. For example, they can produce an artificial feature near $8\,M_\odot$ in the primary mass distribution at $z \lesssim 2$, which is absent when using the full simulation-based S(Z, z), while simultaneously suppressing high-mass features. These discrepancies arise because simple analytical models fail to capture a high-metallicity bump and a more flattened low-metallicity tail in the simulated S(Z, z) metallicity distribution. Our results highlight the importance of accurate star formation histories for modeling BBH populations, demonstrate the limitation of widely used analytical S(Z, z) fits, and underscore the need for careful integration of cosmological simulations, analytical fits, and population synthesis when interpreting gravitational-wave observations.
title From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations
topic High Energy Astrophysical Phenomena
Astrophysics of Galaxies
Solar and Stellar Astrophysics
url https://arxiv.org/abs/2601.20202