redback.transient_models.kilonova_models.power_law_stratified_kilonova

redback.transient_models.kilonova_models.power_law_stratified_kilonova(time, redshift, mej, voffset, v0, alpha, kappa_offset, kappa_0, zeta, beta, **kwargs)[source]

Assumes a power law distribution of ejecta velocities and a power law distribution of kappa corresponding to the ejecta velocity layers for a total of 10 “shells” and calculates the kilonova light curve, using kilonova heating rate. Velocity distribution is flipped so that faster material is the outermost layer as expected for homologous expansion.

Must be used with a constraint to avoid prior draws that predict nonsensical luminosities. Or a sensible prior.

Parameters:

  • time – time in days in observer frame

  • redshift – redshift

  • mej – total ejecta mass in solar masses

  • voffset – minimum ejecta velocity in units of c

  • v0 – velocity normalization in units of c of the power law

  • alpha – power-law index of the velocity distribution i.e., vel = (xs/v0)^-alpha + voffset.

  • kappa_offset – minimum kappa value

  • kappa_0 – kappa normalization

  • zeta – power law index of the kappa distribution i.e., kappa = (xs/kappa_0)^-zeta + kappa_offset

  • beta – power law index of density profile

  • kwargs – Additional keyword arguments

  • frequency – Required if output_format is ‘flux_density’. frequency to calculate - Must be same length as time array or a single number).

  • bands – Required if output_format is ‘magnitude’ or ‘flux’.

  • output_format – ‘flux_density’, ‘magnitude’, ‘spectra’, ‘flux’, ‘sncosmo_source’

  • lambda_array – Optional argument to set your desired wavelength array (in Angstroms) to evaluate the SED on.

  • cosmology – Cosmology to use for luminosity distance calculation. Defaults to Planck18. Must be a astropy.cosmology object.

Returns:

set by output format - ‘flux_density’, ‘magnitude’, ‘spectra’, ‘flux’, ‘sncosmo_source’