The initial-final mass relation (IFMR) of white dwarfs (WDs) plays an important role in stellar evolution. To derive precise estimates of IFMRs and explore how they may vary among star clusters, we propose a Bayesian hierarchical model that pools photo- metric data from multiple star clusters. After performing a simulation study to show the benefits of the Bayesian hierarchical model, we apply this model to five star clus- ters: the Hyades, M67, NGC 188, NGC 2168, and NGC 2477, leading to reasonable and consistent estimates of IFMRs for these clusters. We illustrate how a cluster-specific analysis of NGC 188 using its own photometric data can produce an unreasonable IFMR since its WDs have a narrow range of zero-age main sequence (ZAMS) masses. However, the Bayesian hierarchical model corrects the cluster-specific analysis by bor- rowing strength from other clusters, thus generating more reliable estimates of IFMR parameters. The data analysis presents the benefits of Bayesian hierarchical modelling over conventional cluster-specific methods, which motivates us to elaborate the pow- erful statistical techniques in this article.

Well-sampled optical light curves of 50 gamma-ray bursts (GRBs) with plateau features are compiled from the literature. By empirical fitting, we obtained the parameters of the optical plateaus, such as the decay slopes ($\alpha_{\rm 1}$ and $\alpha_{\rm 2}$), the break times ($T_{\rm b}$), and the corresponding optical fluxes ($F_{\rm b}$) at the break times. The break time of optical plateaus ranges from tens of seconds to $10^6$ seconds, with a typical value about $10^4$ seconds. We have calculated the break luminosity, and it mainly ranges from $10^{44}$ erg $s^{-1}$ to $10^{47}$ erg $s^{-1}$, which is generally two or three orders of magnitude less than the corresponding break luminosity of the X-ray afterglow plateaus. We reanalyzed the optical plateaus and also found that a significantly tighter correlation exists when we added the isotropic equivalent energy of GRBs $E_{\rm \gamma,iso}$ into the $L_{\rm b,z}-T_{\rm b,z}$ relation. The best fit correlation is obtained to be $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm \gamma,iso}^{0.4}$. We next explored the possible correlations among $L_{\rm b,z}$, $T_{\rm b,z}$ and $E_{\rm p,i}$, and found there is also a tight correlation between them, which takes the form of $L_{\rm b,z}\propto T_{\rm b,z}^{-0.9}E_{\rm p,i}^{0.5}$. We argue that these two tight $L_{\rm b,z}-T_{\rm b,z}-E_{\rm \gamma,iso}$ and $L_{\rm b,z}-T_{\rm b,z}-E_{\rm p,i}$ correlations are more physical, and it may be directly related to radiation physics of GRBs. The tight correlations are possible to be used as standard candles.

In astrophysics, we often aim to estimate one or more parameters for each member object in a population and study the distribution of the fitted parameters across the population. In this paper, we develop novel methods that allow us to take advantage of existing software designed for such case-by-case analyses to simultaneously fit parameters of both the individual objects and the parameters that quantify their distribution across the population. Our methods are based on Bayesian hierarchical modelling which is known to produce parameter estimators for the individual objects that are on average closer to their true values than estimators based on case-by-case analyses. We verify this in the context of estimating ages of Galactic halo white dwarfs (WDs) via a series of simulation studies. Finally, we deploy our new techniques on optical and near-infrared photometry of ten candidate halo WDs to obtain estimates of their ages along with an estimate of the mean age of Galactic halo WDs of [11.25, 12.96] Gyr. Although this sample is small, our technique lays the ground work for large-scale studies using data from the Gaia mission.