RECORD
From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations
High Energy Astrophysical Phenomena · Submitted 28 Jan 2026 · Published version
abstract compiled · 2026-10-12 · no language-model review
The authors' main claim 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.
The compiled record
The paper's own sentences, arranged. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
Open the compiled record
The abstract, sentence by sentence
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. 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\gtrsim6$) 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\lesssim2$, 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.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | Observations of binary black hole (BBH) mergers provide a unique window into the lives of massive stars across cosmic time. | background 1.00 | 0.04 | 0.02 | asserted 1.00 |
| s1 | Connecting redshift-dependent merger properties to massive star progenitors requires accurate models of cosmic star formation and chemical enrichment histories. | background 0.98 | 0.05 | 0.03 | asserted 0.99 |
| s2 | 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. | gap 0.69 | 0.04 | 0.04 | asserted 1.00 |
| s3 | Using the IllustrisTNG cosmological simulations, we evaluate the accuracy of these analytical S(Z, z) prescriptions. | method 0.78 | 0.90 | 0.02 | asserted 0.90 |
| s4 | 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. | method 0.95 | 0.90 | 0.03 | asserted 0.82 |
| s5 | We find that analytical S(Z, z) prescriptions can overestimate BBH merger rates at high redshift ($z\gtrsim6$) by up to a factor of $10$-$10^4$, depending on cosmological simulation resolution, and can introduce spurious features in the BBH mass distribution. | result 1.00 | 0.97 | 0.72 | asserted 0.98 |
| s6 | For example, they can produce an artificial feature near $8\,M_\odot$ in the primary mass distribution at $z\lesssim2$, which is absent when using the full simulation-based S(Z, z), while simultaneously suppressing high-mass features. | result 1.00 | 0.92 | 0.10 | asserted 0.72 |
| s7 | 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. | interpretation 0.92 | 0.83 | 0.12 | asserted 1.00 |
| s8 | 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. | interpretation 0.99 | 0.58 | 0.94 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v2
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
The paper
Original abstract
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. 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\gtrsim6$) 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\lesssim2$, 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.
Human review
How the literature uses it
7 arXiv papers cite this in 8 sentences
Papers citing it · 4 sentences without a confident label
Uses its data or values 1 paper
-
The compilation includes the latest LVK constraints on the BBH, BHNS, and BNS merger rates from [ref] , together with simulated rates from a broad range of isolated binary population-synthesis studies [ref] .
Lower Your Rates: On Claims of a Binary Black Hole Merger-Rate Crisis · Method · 0.98
Cites it as support 1 paper
-
Furthermore, strong models for the progenitor formation rate and delay time distribution, even if physically motivated, may fundamentally mismatch the observed redshift evolution of the BBH merger rate [ref] .
Forbidden Formation Histories: The Binary Black Hole Merger Rate Disfavors Long Delay Times · Introduction · 0.96
Cites it as background 1 paper
-
From a theoretical perspective, the local BBH merger rate density has been estimated in numerous studies by combining a 'stellar and binary-evolution term' – representing the outcomes of binary population synthesis codes – with a metallicity-dependent star formation rate density SFRD$(z, Z)$, derived either from cosmological simulations [ref] , or from empirical relations [ref] .
Can current models predict the local black hole merger rate? · Introduction · 0.94
Discussion
Know something about this paper? Leave a comment, or highlight any passage above to comment on it.