Presently a ${>}3\sigma$ tension exists between values of the Hubble constant $H_0$ derived from analysis of fluctuations in the Cosmic Microwave Background by Planck, and local measurements of the expansion using calibrators of type Ia supernovae (SNe Ia). We perform a blinded reanalysis of Riess et al. 2011 to measure $H_0$ from low-redshift SNe Ia, calibrated by Cepheid variables and geometric distances including to NGC 4258. This paper is a demonstration of techniques to be applied to the Riess et at. 2016 data. Our end-to-end analysis starts from available CfA3 and LOSS photometry, providing an independent validation of Riess et al. 2011. We obscure the value of $H_0$ throughout our analysis and the first stage of the referee process, because calibration of SNe Ia requires a series of often subtle choices, and the potential for results to be affected by human bias is significant. Our analysis departs from that of Riess et al. 2011 by incorporating the covariance matrix method adopted in SNLS and JLA to quantify SN Ia systematics, and by including a simultaneous fit of all SN Ia and Cepheid data. We find $H_0 = 72.5 \pm 3.1$ (stat) $\pm 0.77$ (sys) km s$^{-1}$ Mpc$^{-1}$ with a three-galaxy (NGC 4258+LMC+MW) anchor. The relative uncertainties are 4.3% statistical, 1.1% systematic, and 4.4% total, larger than in Riess et al. 2011 (3.3% total) and the Efstathiou 2014 reanalysis (3.4% total). Our error budget for $H_0$ is dominated by statistical errors due to the small size of the supernova sample, whilst the systematic contribution is dominated by variation in the Cepheid fits, and for the SNe Ia, uncertainties in the host galaxy mass dependence and Malmquist bias.

We present a measurement of the volumetric rate of superluminous supernovae (SLSNe) at z~1, measured using archival data from the first four years of the Canada-France-Hawaii Telescope Supernova Legacy Survey (SNLS). We develop a method for the photometric classification of SLSNe to construct our sample. Our sample includes two previously spectroscopically-identified objects, and a further new candidate selected using our classification technique. We use the point-source recovery efficiencies from Perrett et.al. (2010) and a Monte Carlo approach to calculate the rate based on our SLSN sample. We find that the three identified SLSNe from SNLS give a rate of 91 (+76/-36) SNe/Yr/Gpc^3 at a volume-weighted redshift of z=1.13. This is equivalent to 2.2 (+1.8/-0.9) x10^-4 of the volumetric core collapse supernova rate at the same redshift. When combined with other rate measurements from the literature, we show that the rate of SLSNe increases with redshift in a manner consistent with that of the cosmic star formation history. We also estimate the rate of ultra-long gamma ray bursts (ULGRBs) based on the events discovered by the Swift satellite, and show that it is comparable to the rate of SLSNe, providing further evidence of a possible connection between these two classes of events. We also examine the host galaxies of the SLSNe discovered in SNLS, and find them to be consistent with the stellar-mass distribution of other published samples of SLSNe.

We present DES14X3taz, a new hydrogen-poor super luminous supernova (SLSN-I) discovered by the Dark Energy Survey (DES) supernova program, with additional photometric data provided by the Survey Using DECam for Superluminous Supernovae (SUDSS). Spectra obtained using OSIRIS on the Gran Telescopio CANARIAS (GTC) show DES14X3taz is a SLSN-I at z=0.608. Multi-color photometry reveals a double-peaked light curve: a blue and relatively bright initial peak that fades rapidly prior to the slower rise of the main light curve. Our multi-color photometry allows us, for the first time, to show that the initial peak cools from 22,000K to 8,000K over 15 rest-frame days, and is faster and brighter than any published core-collapse supernova, reaching 30% of the bolometric luminosity of the main peak. No physical Nickel powered model can fit this initial peak. We show that a shock-cooling model followed by a magnetar driving the second phase of the light curve can adequately explain the entire light curve of DES14X3taz. Models involving the shock-cooling of extended circumstellar material at a distance of ~400 solar radii are preferred over the cooling of shock-heated surface layers of a stellar envelope. We compare DES14X3taz to the few double-peaked SLSN-I events in the literature. Although the rise-times and characteristics of these initial peaks differ, there exists the tantalizing possibility that they can be explained by one physical interpretation.

Theoretical and observational cosmology have enjoyed a number of significant successes over the last two decades. Cosmic microwave background measurements from the Wilkinson Microwave Anisotropy Probe and Planck, together with large-scale structure and supernova (SN) searches, have put very tight constraints on cosmological parameters. Type Ia supernovae (SNIa) played a central role in the discovery of the accelerated expansion of the Universe, recognised by the Nobel Prize in Physics in 2011. The last decade has seen an enormous increase in the amount of high quality SN observations, with SN catalogues now containing hundreds of objects. This number is expected to increase to thousands in the next few years, as data from next-generation missions, such as the Dark Energy Survey and Large Synoptic Survey Telescope become available. In order to exploit the vast amount of forthcoming high quality data, it is extremely important to develop robust and efficient statistical analysis methods to answer cosmological questions, most notably determining the nature of dark energy. To address these problems my work is based on nested-sampling approaches to parameter estimation and model selection and neural networks for machine-learning. Using advanced Bayesian techniques, I constrain the properties of dark-matter haloes along the SN lines-of-sight via their weak gravitational lensing effects, develop methods for classifying SNe photometrically from their lightcurves, and present results on more general issues associated with constraining cosmological parameters and testing the consistency of different SN compilations.

It is now common practice to constrain cosmological parameters using supernovae (SNe) catalogues constructed from several different surveys. Before performing such a joint analysis, however, one should check that parameter constraints derived from the individual SNe surveys that make up the catalogue are mutually consistent. We describe a statistically-robust mutual consistency test, which we calibrate using simulations, and apply it to each pairwise combination of the surveys making up, respectively, the UNION2 catalogue and the very recent JLA compilation by Betoule et al. We find no inconsistencies in the latter case, but conclusive evidence for inconsistency between some survey pairs in the UNION2 catalogue.

A hierarchical Bayesian method is applied to the analysis of Type-Ia supernovae (SNIa) observations to constrain the properties of the dark matter haloes of galaxies along the SNIa lines-of-sight via their gravitational lensing effect. The full joint posterior distribution of the dark matter halo parameters is explored using the nested sampling algorithm {\sc MultiNest}, which also efficiently calculates the Bayesian evidence, thereby facilitating robust model comparison. We first demonstrate the capabilities of the method by applying it to realistic simulated SNIa data, based on the real 3-year data release from the Supernova Legacy Survey (SNLS3). Assuming typical values for the halo parameters in our simulations, we find that a catalogue analogous to the existing SNLS3 data set is incapable of detecting the lensing signal, but a catalogue containing approximately three times as many SNIa does produce robust and accurate parameter constraints and model selection results for two halo models: a truncated singular isothermal sphere (SIS) and a Navarro--Frenk--White (NFW) profile, thereby validating our analysis methodology. In the analysis of the real SNLS3 data, contrary to previous studies, we obtain only a very marginal detection of a lensing signal and weak constraints on the halo parameters for the truncated SIS model, although these constraints are tighter than those obtained from the equivalent simulated SNIa data set. This difference is driven by a preferred value of $\eta \approx 1$ in the assumed scaling-law $\sigma \propto L^\eta$ between velocity dispersion and luminosity, which is somewhat higher than the canonical values of $\eta = \tfrac{1}{4}$ and $\eta = \tfrac{1}{3}$ for early and late-type galaxies, respectively, and leads to a stronger lensing effect by the halo. No detection of a lensing signal is made for the NFW model.

A method is presented for automated photometric classification of supernovae (SNe) as Type-Ia or non-Ia. A two-step approach is adopted in which: (i) the SN lightcurve flux measurements in each observing filter are fitted separately; and (ii) the fitted function parameters and their associated uncertainties, along with the number of flux measurements, the maximum-likelihood value of the fit and Bayesian evidence for the model, are used as the input feature vector to a classification neural network (NN) that outputs the probability that the SN under consideration is of Type-Ia. The method is trained and tested using data released following the SuperNova Photometric Classification Challenge (SNPCC). We consider several random divisions of the data into training and testing sets: for instance, for our sample D_1 (D_4), a total of 10% (40%) of the data are involved in training the algorithm and the remainder used for blind testing of the resulting classifier; we make no selection cuts. Assigning a canonical threshold probability of p_th=0.5 on the NN output to classify a SN as Type-Ia, for the sample D_1 (D_4) we obtain a completeness of 0.78 (0.82), purity of 0.77 (0.82), and SNPCC figure-of-merit of 0.41 (0.50). Including the SN host-galaxy redshift and its uncertainty as additional inputs to the NN results in a modest 5-10% increase in these values. We find that the classification quality does not vary significantly with SN redshift. Moreover, our probabilistic classification method allows one to calculate the expected completeness, purity or other measures of classification quality as a function of the p_th, without knowing the true classes of the SNe in the testing sample. The method may thus be improved further by optimising p_th and can easily be extended to divide non-Ia SNe into their different classes.

We present a comparison of two methods for cosmological parameter inference from supernovae Ia lightcurves fitted with the SALT2 technique. The standard chi-square methodology and the recently proposed Bayesian hierarchical method (BHM) are each applied to identical sets of simulations based on the 3-year data release from the Supernova Legacy Survey (SNLS3), and also data from the Sloan Digital Sky Survey (SDSS), the Low Redshift sample and the Hubble Space Telescope (HST), assuming a concordance LCDM cosmology. For both methods, we find that the recovered values of the cosmological parameters, and the global nuisance parameters controlling the stretch and colour corrections to the supernovae lightcurves, suffer from small biasses. The magnitude of the biasses is similar in both cases, with the BHM yielding slightly more accurate results, in particular for cosmological parameters when applied to just the SNLS3 single survey data sets. Most notably, in this case, the biasses in the recovered matter density $\Omega_{\rm m,0}$ are in opposite directions for the two methods. For any given realisation of the SNLS3-type data, this can result in a $\sim 2 \sigma$ discrepancy in the estimated value of $\Omega_{\rm m,0}$ between the two methods, which we find to be the case for real SNLS3 data. As more higher and lower redshift SNIa samples are included, however, the cosmological parameter estimates of the two methods converge.