Adds a new model for compact binary merger rates to the gravitational wave(s) module, from Artale et al (2020).
There are three models, depending cumulatively in each of stellar mass, star formation rate and metallicity. Each model is a simple logarithmic model with coefficients for the log of each independent variable and a constant.
The values of the parameters are fit to objects in the EAGLE simulation suite, with independent fits on simulation output at four redshifts: 0.1, 1, 2, 6. At the moment a simple linear interpolation between these points is performed for the coefficients as a function of redshift.
ToDo
[x] Finish implementing and documenting SFR and metallicity-dependent models
[ ] Expand tests
[x] Decide on consistent treatment of units
[ ] Improve treatment of redshift dependence (function decorators?)
Description
Adds a new model for compact binary merger rates to the gravitational wave(s) module, from Artale et al (2020).
There are three models, depending cumulatively in each of stellar mass, star formation rate and metallicity. Each model is a simple logarithmic model with coefficients for the log of each independent variable and a constant.
The values of the parameters are fit to objects in the EAGLE simulation suite, with independent fits on simulation output at four redshifts: 0.1, 1, 2, 6. At the moment a simple linear interpolation between these points is performed for the coefficients as a function of redshift.
ToDo
Checklist