redback.transient_models.supernova_models.arnett_with_features

redback.transient_models.supernova_models.arnett_with_features(time, redshift, f_nickel, mej, **kwargs)[source]

A version of the arnett model where SED has time-evolving spectral features.

Parameters:

  • time – time in days

  • redshift – source redshift

  • f_nickel – fraction of nickel mass

  • mej – total ejecta mass in solar masses

  • kwargs – Must be all the kwargs required by the specific interaction_process, photosphere, sed methods used e.g., for Diffusion and TemperatureFloor: kappa, kappa_gamma, vej (km/s), temperature_floor

  • interaction_process – Default is Diffusion. Can also be None in which case the output is just the raw engine luminosity, or another interaction process.

  • photosphere – Default is TemperatureFloor. kwargs must have vej or relevant parameters if using different photosphere model

  • sed – Default is BlackbodyWithFeatures for Type Ia spectra.

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

  • feature_list – Optional list of spectral features. If None, uses default Type Ia features.

Feature Parameters (dynamically numbered): Features are defined by groups of parameters with pattern: {param}_feature_{N} where N starts from 1. All features with the same N are grouped together.

Required for each feature N: :param rest_wavelength_feature_N: Central wavelength in Angstroms :param sigma_feature_N: Gaussian width in Angstroms :param amplitude_feature_N: Amplitude (negative=absorption, positive=emission), percentage of continuum (e.g., -0.4 = 40% absorption) :param t_start_feature_N: Start time in source-frame days :param t_end_feature_N: End time in source-frame days

Optional for each feature N (smooth mode only): :param t_rise_feature_N: Rise time in source-frame days (default: 2.0) :param t_fall_feature_N: Fall time in source-frame days (default: 5.0)

General parameters: :param evolution_mode: ‘smooth’ or ‘sharp’ (default: ‘smooth’) :param use_default_features: If True and no custom features found, use defaults (default: True)

Examples:

# Single custom feature result = model(time, z, f_ni, mej,

rest_wavelength_feature_1=6355.0, sigma_feature_1=400.0, amplitude_feature_1=-0.4, t_start_feature_1=0, t_end_feature_1=30, output_format=’magnitude’, bands=’lsstg’)

# Multiple features result = model(time, z, f_ni, mej,

rest_wavelength_feature_1=6355.0, sigma_feature_1=400.0, amplitude_feature_1=-0.4, t_start_feature_1=0, t_end_feature_1=40, rest_wavelength_feature_2=3934.0, sigma_feature_2=300.0, amplitude_feature_2=-0.5, t_start_feature_2=0, t_end_feature_2=60, rest_wavelength_feature_3=8600.0, sigma_feature_3=500.0, amplitude_feature_3=-0.3, t_start_feature_3=0, t_end_feature_3=50, evolution_mode=’smooth’, output_format=’magnitude’, bands=’lsstg’)

return:

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