The Locations of Features in the Mass Distribution of Merging Binary Black Holes are Robust against Uncertainties in the Metallicity-dependent Cosmic Star Formation History

New observational facilities are probing astrophysical transients such as stellar explosions and gravitational-wave sources at ever-increasing redshifts, while also revealing new features in source property distributions. To interpret these observations, we need to compare them to predictions from stellar population models. Such models require the metallicity-dependent cosmic star formation history (S(Z , z )) as an input. Large uncertainties remain in the shape and evolution of this function. In this work, we propose a simple analytical function for (S(Z , z). Variations of this function can be easily interpreted because the parameters link to its shape in an intuitive way. We fit our analytical function to the star-forming gas of the cosmological TNG100 simulation and find that it is able to capture the main behavior well. As an example application, we investigate the effect of systematic variations in the (S(Z , z) parameters on the predicted mass distribution of locally merging binary black holes. Our main findings are that (i) the locations of features are remarkably robust against variations in the metallicity-dependent cosmic star formation history, and (ii) the low-mass end is least affected by these variations. This is promising as it increases our chances of constraining the physics that govern the formation of these objects (https://github.com/LiekeVanSon/SFRDfit/tree/7348a1ad0d2ed6b78c70d5100fb3cd2515493f02/).