redback.transient_models.kilonova_models.nicholl_bns

redback.transient_models.kilonova_models.nicholl_bns(time, redshift, mass_1, mass_2, lambda_s, kappa_red, kappa_blue, mtov, epsilon, alpha, cos_theta, cos_theta_cocoon, temperature_floor_1, temperature_floor_2, temperature_floor_3, **kwargs)[source]

Kilonova model from Nicholl et al. 2021, inclides three kilonova components + shock heating from cocoon + disk winds from remnant

Parameters:

  • time – time in days in observer frame

  • redshift – redshift

  • mass_1 – Mass of primary in solar masses

  • mass_2 – Mass of secondary in solar masses

  • lambda_s – Symmetric tidal deformability i.e, lambda_s = lambda_1 + lambda_2

  • kappa_red – opacity of the red ejecta

  • kappa_blue – opacity of the blue ejecta

  • mtov – Tolman Oppenheimer-Volkoff mass in solar masses

  • epsilon – fraction of disk that gets unbound/ejected

  • alpha – Enhancement of blue ejecta by NS surface winds if mtotal < prompt collapse, can turn off by setting alpha=1

  • cos_theta – Viewing angle of observer

  • cos_theta_cocoon – Opening angle of shocked cocoon

  • temperature_floor_1 – Temperature floor of first (blue) component

  • temperature_floor_2 – Temperature floor of second (purple) component

  • temperature_floor_3 – Temperature floor of third (red) component

  • kwargs – Additional keyword arguments

  • dynamical_ejecta_error – Error in dynamical ejecta mass, default is 1 i.e., no error in fitting formula

  • disk_ejecta_error – Error in disk ejecta mass, default is 1 i.e., no error in fitting formula

  • shocked_fraction – Fraction of ejecta that is shocked by jet, default is 0.2 i.e., 20% of blue ejecta is shocked. Use 0. if you want to turn off cocoon emission.

  • nn – ejecta power law density profile, default is 1.

  • tshock – time for shock in source frame in seconds, default is 1.7s (see Nicholl et al. 2021)

  • 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.

  • dense_resolution – resolution of the grid that the model is actually evaluated on, default is 300

  • 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’