skypy.gravitational_wave.b_band_merger_rate

skypy.gravitational_wave.b_band_merger_rate(luminosity, population='NS-NS', optimism='low')

Model of Abadie et al (2010), Table III

Compact binary merger rates as a linear function of a galaxies B-band luminosity.

Parameters

luminosity(ngal,) array-like

The B-band luminosity of the galaxies to generate merger rates for, in units of solar luminosity.

population{‘NS-NS’, ‘NS-BH’, ‘BH-BH’}

Compact binary population to get rate for. ‘NS-NS’ is neutron star - neutron star ‘NS-BH’ is neutron star - black hole ‘BH-BH’ is black hole - black hole

optimism{‘low’, ‘realistic’, ‘high’}

Optimism of predicted merger rates. For ‘NS-NS’ there is an extra option ‘max’.

Returns

merger_ratearray_like

Merger rates for the galaxies in units of year^-1

References

Examples

>>> import numpy as np
>>> from skypy.gravitational_wave import b_band_merger_rate

Sample 100 luminosity values near absolute magnitude -20.5.

>>> luminosities = 10.**(-0.4*(-20.5 + np.random.randn(100)))

Generate merger rates for these luminosities.

>>> rates = b_band_merger_rate(luminosities,
...                            population='NS-NS',
...                            optimism='low')