Kilonova#

class uvex_transients.transients.kilonovae.Kilonova(cosmology: Cosmology | None = None)[source]#

GW170817-like kilonova population.

Modeled with KilonovaCoolingBlackbodySED (a Gaussian-rise/broken-power-law-decline light curve with a cooling blackbody photosphere, calibrated against GW170817), a constant volumetric rate of \(110\ \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}\) out to \(z=0.2\), and a 30-day duration window. This rate is Fishbach et al. 2026’s total BNS merger rate, not just their GW170817-like mass-bin sub-rate, so every BNS merger sampled here is assumed to be feasibly GW170817-like: a simplifying assumption of this simulation, not a claim made by that paper.

The \(z=0.2\) redshift limit is conservative: a source with bolometric luminosity below \(10^{43}\ \mathrm{erg\,s^{-1}}\) (well above anything reported in the literature) would still remain detectable at UVEX’s \(m<27\) band limit out to \(z=2\), assuming a flat spectrum and no K-correction. The 30-day duration is a conservative upper bound on the total light curve; the blue/early kilonova component this SED targets fades below detectability closer to \(\sim10\) days.

See also

uvex_transients.models.kilonovae.kne.KilonovaCoolingBlackbodySED: The SED this class pairs with DEFAULT_MODEL.
uvex_transients.transients.base.ExtragalacticTransient: The base class supplying rate/sampling machinery.

Methods

`compute_all_sky_yield`(duration)	Compute the expected intrinsic event count over the full sky.
`compute_all_sky_yield_ci`(duration)	Compute rate-only bounds on the expected all-sky intrinsic event count.
`event_rate`(z)	Compute the comoving event rate density at redshift(s) z, \(R(z;A)=A f(z)\).
`rate_shape`(z)	Return the (trivial, constant) rate shape of kilonovae at a given redshift.
`registry`()	dict[str, type[TransientBase]]: A copy of every concrete subclass imported so far, keyed by class name.
`sample_event_count`(solid_angle[, duration, ...])	Sample the number of events in a given solid angle and duration.
`sample_event_redshift`(n_samples, *[, rng])	Draw n_samples redshifts from the rate-weighted redshift distribution.
`sample_events_on_healpix_grid`(nside, *, t_start)	Draw a Monte Carlo realization of events on a HEALPix grid, over a time window.

Attributes

`DEFAULT_DURATION`	The duration of this transient.
`DEFAULT_Z_GRID_SIZE`
`DEFAULT_Z_LIM`
`RATE_CI`	Multiplicative (lower, upper) 90% confidence factors on rate.
`all_sky_rate`	Expected all-sky event rate per observer year.
`all_sky_rate_ci`	Rate-only bounds on `all_sky_rate`.
`cosmology`	The cosmology used for luminosity-distance/volume calculations.
`duration_limit`	Upper bound on this transient's total duration.
`effective_volume`	\(\mathcal V\), the all-sky rate-weighted comoving volume.
`integrated_rate`	Expected event rate per observer time and solid angle.
`integrated_rate_ci`	Rate-only bounds on `integrated_rate`.
`luminosity_distance_grid`	\(D_L(z)\) at each point of redshift_grid.
`rate`	The volumetric kilonova rate, constant in \(z\).
`rate_ci`	`(R_L, R_U)`, the confidence bounds on rate.
`redshift_grid`	The cached redshift grid backing integrated_rate/sampling.
`redshift_grid_size`	Number of points in redshift_grid; reassigning invalidates the rate cache.
`redshift_limit`	Upper redshift bound of redshift_grid; reassigning invalidates the rate cache.
`sed`	This instance's SED -- flux, magnitude, and spectrum evaluation all live here.