Unsupervised Transient Light Curve Analysis Via Hierarchical Bayesian Inference

Nathan Sanders,Michael Betancourt,Alicia Soderberg

DOI: https://doi.org/10.1088/0004-637X/800/1/36

2014-04-14

Abstract:Historically, light curve studies of supernovae (SNe) and other transient classes have focused on individual objects with copious and high signal-to-noise observations. In the nascent era of wide field transient searches, objects with detailed observations are decreasing as a fraction of the overall known SN population, and this strategy sacrifices the majority of the information contained in the data about the underlying population of transients. A population level modeling approach, simultaneously fitting all available observations of objects in a transient sub-class of interest, fully mines the data to infer the properties of the population and avoids certain systematic biases. We present a novel hierarchical Bayesian statistical model for population level modeling of transient light curves, and discuss its implementation using an efficient Hamiltonian Monte Carlo technique. As a test case, we apply this model to the Type IIP SN sample from the Pan-STARRS1 Medium Deep Survey, consisting of 18,837 photometric observations of 76 SNe, corresponding to a joint posterior distribution with 9,176 parameters under our model. Our hierarchical model fits provide improved constraints on light curve parameters relevant to the physical properties of their progenitor stars relative to modeling individual light curves alone. Moreover, we directly evaluate the probability for occurrence rates of unseen light curve characteristics from the model hyperparameters, addressing observational biases in survey methodology. We view this modeling framework as an unsupervised machine learning technique with the ability to maximize scientific returns from data to be collected by future wide field transient searches like LSST. \smallskip

Instrumentation and Methods for Astrophysics,High Energy Astrophysical Phenomena

What problem does this paper attempt to address?

The problem that this paper attempts to solve is how to use the hierarchical Bayesian inference method to conduct unsupervised analysis of the light - curve of transient sources in wide - field transient source surveys, in order to overcome the systematic biases in traditional methods and make full use of the information in the data to infer the properties of the transient source population. Specifically:

1. **Reducing systematic biases**: Traditional methods usually focus on individual celestial bodies with abundant observational data and high signal - to - noise ratios. This method sacrifices most of the information about the overall transient sources. The method proposed in this paper fully mines data information by simultaneously fitting all available observational data in a certain transient sub - class, avoiding some systematic biases.

2. **Improving parameter constraints**: Through the hierarchical Bayesian model, the light - curve parameters related to the physical properties of supernova progenitor stars can be better constrained. Compared with modeling each light - curve separately, this method provides more stringent parameter constraints.

3. **Evaluating the incidence of unobserved light - curve features**: Directly evaluate the occurrence probability of unobserved light - curve features from the hyper - parameters of the model, thereby solving the observational bias problem in the survey method.

4. **Maximizing scientific returns**: This modeling framework is regarded as an unsupervised machine - learning technique and can maximize the scientific returns of future wide - field transient source surveys (such as LSST) data.

In summary, this paper aims to achieve population - level modeling of the light - curves of transient sources by developing a new hierarchical Bayesian statistical model combined with an efficient Hamiltonian Monte Carlo technique, thereby reducing systematic biases and at the same time improving the understanding of the physical properties of supernova progenitor stars.