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

Introduction to the show: ident: Astrophysics Radio. Generated commentary on the latest astrophysics papers.

Vera: Today's paper: "From cosmological simulations to binary black hole mergers".

Jocelyn: The impact of using analytical star formation history models on gravitational-wave source populations is investigated by coupling cosmological simulations with population synthesis codes to assess how these analytical prescriptions affect…

Vera: First, who's behind it and why it matters.

Title and authors: Jocelyn: So they are diving into what this paper actually does. They take analytical prescriptions for S(Z, z)—which is just their way of approximating how much star formation happened at different redshifts and metallicities—and they compare that against what the simulations themselves generated.

Vera: And they found that these analytical fits can cause some pretty big changes in the predicted merger rates depending on which model you use. They're looking at how different S(Z, z) prescriptions affect the BBH merger rates.

Subrahmanyan: This is important because if we rely too much on simple analytical models without checking them against the simulation data, our predictions for when and where these mergers happen could be skewed by a lot.

Jocelyn: They specifically point out that at high redshift, around z greater than six the analytical fits give merger rates that are systematically larger than what you get directly from the TNG simulation itself <ref:2601.20202#pg1>.

Vera: How big of a difference are we talking here? The paper gives some pretty striking numbers when they compare different simulations.

The paper's summary: Subrahmanyan: The researchers summarize that they constructed their analytical model for S(Z, z) by combining the total star formation rate density, S(z), with a skewed log-normal metallicity distribution, dP/dZ. This is the specific setup they used to test things.

Jocelyn: So it’s not just one simple fit; they are using a combination of elements to build their analytical model for the star formation history that they then plug into their merger rate calculations using COMPAS.

Vera: The summary highlights that these analytical prescriptions can cause significant changes in BBH merger rates, and this change depends heavily on which specific model they chose to use for S(Z, z).

Subrahmanyan: They found that at high redshift, specifically z greater than six the analytical fits yield systematically larger BBH merger rates compared to those derived directly from the IllustrisTNG simulation <ref:2601.20202#pg1>.

Jocelyn: And it’s not just a slight difference; they found these factors can be as large as one hundred four times bigger for the TNG300-one simulation, and around ten times bigger for the TNG50-one and TNG100-one simulations.

Vera: That factor of one hundred four is a big one, and it shows that you can get really far off if you aren't careful about your star formation history approximations.

The paper's improvements: Vera: So what do the authors suggest we should do differently? They look at how to fix these issues in their analytical models for S(Z, z). The main critique is that simple analytical models fail to capture some key features of the underlying data.

Subrahmanyan: Specifically, they say that these simple fits don't manage to capture the bimodality of S(Z), which is how star formation usually shows up across different redshifts.

Jocelyn: And they also point out that these simple models fail to account for the truncation of the low-metallicity tail in their metallicity distribution, which is another important part of that history.

Vera: They suggest future analytical fits need to be much more flexible. They should allow for a secondary peak in the metallicity distribution and a flatter low-metallicity tail instead of just one simple log-normal distribution.

Subrahmanyan: That flexibility would help them better capture the structure of S(Z, z) that they see in their simulations, which is crucial because those features affect the resulting BBH mass distribution too.

Jocelyn: They also suggest improving how we optimize these fits to balance getting an accurate picture of the main peak of star formation with accurately describing that low-metallicity tail.

Conclusion: Vera: So to wrap up, the central message from this paper is a caution about using simple analytical models for star formation history when predicting black hole mergers. They show that these simple models can lead to serious overestimations of merger rates at high redshift.

Subrahmanyan: The implication is that we need more complex modeling tools for S(Z, z) if we want robust predictions about the BBH population, because those simple approximations miss important physical details about how galaxies evolve.

Jocelyn: They emphasize that caution is needed when using these simple analytical fits to represent the more complex data coming out of cosmological simulations.

Vera: And they suggest that future work should look at a direct comparison with observational S(Z, z) distributions, to see if the features where these analytical fits fail in simulations actually exist in what we observe.

Subrahmanyan: I think the paper also hints that binary evolution uncertainties often dominate over the uncertainties in S(Z, z) when it comes to how much mergers happen.

Jocelyn: And they also noted that resolution matters; lower resolution simulations show a bigger disagreement between the analytical fit and the simulation's actual star formation history.

Vera: So, this paper on "From cosmological simulations to binary black hole mergers: The impact of using analytical star formation history models on gravitational-wave source populations" is a reminder that the complexity of galaxy evolution needs complex modeling to get accurate results for these kinds of astrophysical populations.