Magnetar+Bump Modeling of Superluminous Supernovae

These data files are associated with the article "Bumpy Declining Light Curves are Common in Hydrogen-poor Superluminous Supernovae." This repository contains four directories, described below.

MOSFiT Output Files

The "mosfit" directory contains the raw output files from the Modular Open-Source Fitter for Transients (Guillochon et al. 2017). See its documentation for the output format. Keep in mind that some detections have been converted into nondetections in these files (see the paper for a full explanation), so you should not use them as a source for the photometry.

Magnetar Models

The "models" directory contains the range of models for each supernova. These are stored as ECSV files, which can be read using Astropy tables. The columns are MJD, filter, apparent magnitude (median, minimum, maximum), and absolute magnitude (median, minimum, maximum). The metadata contains the median and standard deviation of each parameter discussed in the paper, as well as the supernova name, distance modulus, and extinction.

Magnetar Parameters

The "params" directory contains the magnetar model parameters corresponding to each model realization, also stored as ECSV files. In other words, you can use these files to reproduce the posterior for each parameter. The metadata contains the units for each parameter and a LaTeX representation of each variable name.

Bump Models

The "bumps" directory contains the raw output for the bump modeling. These are stored as compressed Numpy files (.npz) and can be read using np.load. The keywords are as follows:

"chain": the flattened MCMC chain, stored as a 2D array. The 0th dimension corresponds to the number of walkers times the number of steps. The 1st dimension corresponds to the parameters: amplitude, center (MJD), and FWHM (observed days). In the case of 2 bumps, the parameters are amplitude 1, center 1, FWHM 1, amplitude 2, center 2, FWHM 2.

"pmed": the median of each parameter, followed by the median absolute deviation for each parameter (converted to a standard deviation).

"energies": the integrated energy in each model realization (in watt-days).

"energy": the median of the “energies” array.

"denergy": the median absolute deviation of the “energies” array (converted to a standard deviation).

"integral": the integral of the magnetar residuals before modeling (in watt-days). This is not used in the paper.

"t0" and "t1": the range of MJDs used in the bump fitting.

If you have any questions about this dataset, feel free to contact the author, Griffin Hosseinzadeh.