In the last decade, astronomers have found a new type of supernova called `superluminous supernovae' (SLSNe) due to their high peak luminosity and long light-curves. These hydrogen-free explosions (SLSNe-I) can be seen to z~4 and therefore, offer the possibility of probing the distant Universe. We aim to investigate the possibility of detecting SLSNe-I using ESA's Euclid satellite, scheduled for launch in 2020. In particular, we study the Euclid Deep Survey (EDS) which will provide a unique combination of area, depth and cadence over the mission. We estimated the redshift distribution of Euclid SLSNe-I using the latest information on their rates and spectral energy distribution, as well as known Euclid instrument and survey parameters, including the cadence and depth of the EDS. We also applied a standardization method to the peak magnitudes to create a simulated Hubble diagram to explore possible cosmological constraints. We show that Euclid should detect approximately 140 high-quality SLSNe-I to z ~ 3.5 over the first five years of the mission (with an additional 70 if we lower our photometric classification criteria). This sample could revolutionize the study of SLSNe-I at z>1 and open up their use as probes of star-formation rates, galaxy populations, the interstellar and intergalactic medium. In addition, a sample of such SLSNe-I could improve constraints on a time-dependent dark energy equation-of-state, namely w(a), when combined with local SLSNe-I and the expected SN Ia sample from the Dark Energy Survey. We show that Euclid will observe hundreds of SLSNe-I for free. These luminous transients will be in the Euclid data-stream and we should prepare now to identify them as they offer a new probe of the high-redshift Universe for both astrophysics and cosmology.

We consider the dispersion on the supernova distance-redshift relation due to peculiar velocities and gravitational lensing, and the sensitivity of these effects to the amplitude of the matter power spectrum. We use the MeMo lensing likelihood developed by Quartin, Marra & Amendola (2014), which accounts for the characteristic non-Gaussian distribution caused by lensing magnification with measurements of the first four central moments of the distribution of magnitudes. We build on the MeMo likelihood by including the effects of peculiar velocities directly into the model for the moments. In order to measure the moments from sparse numbers of supernovae, we take a new approach using Kernel Density Estimation to estimate the underlying probability density function of the magnitude residuals. We also describe a bootstrap re-sampling approach to estimate the data covariance matrix. We then apply the method to the Joint Light-curve Analysis (JLA) supernova catalogue. When we impose only that the intrinsic dispersion in magnitudes is independent of redshift, we find $\sigma_8=0.44^{+0.63}_{-0.44}$ at the one standard deviation level, although we note that in tests on simulations, this model tends to overestimate the magnitude of the intrinsic dispersion, and underestimate $\sigma_8$. We note that the degeneracy between intrinsic dispersion and the effects of $\sigma_8$ is more pronounced when lensing and velocity effects are considered simultaneously, due to a cancellation of redshift dependence when both effects are included. Keeping the model of the intrinsic dispersion fixed as a Gaussian distribution of width 0.14 mag, we find $\sigma_8 = 1.07^{+0.50}_{-0.76}$.

We study the feasibility of detecting weak lensing spatial correlations between Supernova (SN) Type Ia magnitudes with present (Dark Energy Survey, DES) and future (Large Synoptic Survey Telescope, LSST) surveys. We investigate the angular auto-correlation function of SN magnitudes (once the background cosmology has been subtracted) and cross-correlation with galaxy catalogues. We examine both analytical and numerical predictions, the latter using simulated galaxy catalogues from the MICE Grand Challenge Simulation. We predict that we will be unable to detect the SN auto-correlation in DES, while it should be detectable with the LSST SN deep fields (15,000 SNe on 70 deg^2) at ~6sigma level of confidence (assuming 0.15 magnitudes of intrinsic dispersion). The SN-galaxy cross-correlation function will deliver much higher signal-to-noise, being detectable in both surveys with an integrated signal-to-noise of ~100 (up to 30 arcmin separations). We predict joint constraints on the matter density parameter (Omega_m) and the clustering amplitude (sigma_8) by fitting the auto-correlation function of our mock LSST deep fields. When assuming a Gaussian prior for Omega_m, we can achieve a 25% measurement of sigma_8 from just these LSST supernovae (assuming 0.15 magnitudes of intrinsic dispersion). These constraints will improve significantly if the intrinsic dispersion of SNe Ia can be reduced.

We predict cosmological constraints for forthcoming surveys using Superluminous Supernovae (SLSNe) as standardisable candles. Due to their high peak luminosity, these events can be observed to high redshift (z~3), opening up new possibilities to probe the Universe in the deceleration epoch. We describe our methodology for creating mock Hubble diagrams for the Dark Energy Survey (DES), the "Search Using DECam for Superluminous Supernovae" (SUDSS) and a sample of SLSNe possible from the Large Synoptic Survey Telescope (LSST), exploring a range of standardisation values for SLSNe. We include uncertainties due to gravitational lensing and marginalise over possible uncertainties in the magnitude scale of the observations (e.g. uncertain absolute peak magnitude, calibration errors). We find that the addition of only ~100 SLSNe from SUDSS to 3800 Type Ia Supernovae (SNe Ia) from DES can improve the constraints on w and Omega_m by at least 20% (assuming a flat wCDM universe). Moreover, the combination of DES SNe Ia and 10,000 LSST-like SLSNe can measure Omega_m and w to 2% and 4% respectively. The real power of SLSNe becomes evident when we consider possible temporal variations in w(a), giving possible uncertainties of only 2%, 5% and 14% on Omega_m, w_0 and w_a respectively, from the combination of DES SNe Ia, LSST-like SLSNe and Planck. These errors are competitive with predicted Euclid constraints, indicating a future role for SLSNe for probing the high redshift Universe.

We assume the DE state equations w(a) = w_0+w_a(a_p-a), and study the dependence of the constraints on w_0 and w_a coefficients on the pivoting redshift 1+z_p=1/a_p. Coefficients are fitted to data including WMAP7, SNIa (Union 2.1), BAO's (including WiggleZ and SDSS results) and H_0 constraints. The fitting algorithm is CosmoMC. We find specific differences between the cases when neutrino mass is allowed or disregarded. More in detail: i) The z_p value yielding uncorrelated constraints on w_0 and w_a is different in the two cases, holding ~0.25 and ~0.35, respectively. (ii) If we consider the intervals allowed to w_0, we find that they shift when z_p increases, in opposite directions for vanishing or allowed neutrino mass. This leads to no overlap between 1sigma intervals already at z_p >~0.4. (iii) The known effect that a more negative state parameter is required to allow for neutrino mass displays its effects on w_a, rather than on w_0. (iv) The w_0-w_a constraints found by using any pivot z_p can be translated into constraints holding at a specific z_p value (0 or the z_p where errors are uncorrelated). When we do so, error ellipses exhibit a satisfactory overlap.