In 2016, the first strongly lensed Type Ia supernova, iPTF16geu at redshift $z=0.409$ with four resolved images arranged symmetrically around the lens galaxy at $z=0.2163$, was discovered. Here, refined observations of iPTF16geu, including the time delay between images, are used to decrease uncertainties in the lens model, including the the slope of the projected surface density of the lens galaxy, $\Sigma\propto r^{1-\eta}$, and to constrain the universal expansion rate $H_0$. Imaging with HST provides an upper limit on the slope $\eta$, in slight tension with the steeper density profiles indicated by imaging with Keck after iPTF16geu had faded, potentially due to dust extinction not corrected for in host galaxy imaging. Since smaller $\eta$ implies larger magnifications, we take advantage of the standard candle nature of Type Ia supernovae constraining the image magnifications, to obtain an independent constraint of the slope. We find that a smooth lens density fails to explain the iPTF16geu fluxes, regardless of the slope, and additional sub-structure lensing is needed. The total probability for the smooth halo model combined with star microlensing to explain the iPTF16geu image fluxes is maximized at $12\,\%$ for $\eta\sim 1.8$, in excellent agreement with Keck high spatial resolution data, and flatter than an isothermal halo. It also agrees perfectly with independent constraints on the slope from lens velocity dispersion measurements. Combining with the observed time delays between the images, we infer a lower bound on the Hubble constant, $H_0 \gtrsim 40\,{\rm km \ s^{-1} Mpc^{-1}}$ at $68.3\,\%$ confidence level.

Bimetric gravity can reproduce the accelerated expansion of the Universe, without a cosmological constant. However, the stability of these solutions to linear perturbations has been questioned, suggesting exponential growth of structure in this approximation. We present a simple model of structure formation, for which an analytical solution is derived. The solution is well-behaved, showing that there is no physical instability with respect to these perturbations. The model can yield a growth of structure exhibiting measurable differences from $\Lambda$CDM.

We report the results from spectroscopic observations of the multiple images of the strongly lensed Type Ia supernova (SN Ia), iPTF16geu, obtained with ground based telescopes and the Hubble Space Telescope (HST). From a single epoch of slitless spectroscopy with HST, we can resolve spectra of individual lensed supernova images for the first time. This allows us to perform an independent measurement of the time-delay between the two brightest images, $\Delta t = 1.4 \pm 5.0$ days, which is consistent with the time-delay measured from the light-curves. We also present measurements of narrow emission and absorption lines characterizing the interstellar medium in the host galaxy at z=0.4087, as well as in the foreground lensing galaxy at z=0.2163. We detect strong Na ID absorption in the host galaxy, indicating that iPTF16geu belongs to a subclass of SNe Ia displaying "anomalously" large Na ID column densities in comparison to the amount of dust extinction derived from their light curves. For the deflecting galaxy, we refine the measurement of the velocity dispersion, $\sigma = 129 \pm 4$ km/s, which significantly constrains the lens model. Since the time-delay between the SN images is negligible, we can use unresolved ground based spectroscopy, boosted by a factor ~70 from lensing magnification, to study the properties of a high-z SN Ia with unprecedented signal-to-noise ratio. The spectral properties of the supernova, such as pseudo-Equivalent widths of several absorption features and velocities of the Si II-line indicate that iPTF16geu, besides being lensed, is a normal SN Ia, indistinguishable from well-studied ones in the local universe, providing support for the use of SNe Ia in precision cosmology. We do not detect any significant deviations of the SN spectral energy distribution from microlensing of the SN photosphere by stars and compact objects in the lensing galaxy.

Numerical integration of the field equations in bimetric relativity is necessary to obtain solutions describing realistic systems. Thus, it is crucial to recast the equations as a well-posed problem. In general relativity, under certain assumptions, the covariant BSSN formulation is a strongly hyperbolic formulation of the Einstein equations, hence its Cauchy problem is well-posed. In this paper, we establish the covariant BSSN formulation of the bimetric field equations. It shares many features with the corresponding formulation in general relativity, but there are a few fundamental differences between them. Some of these differences depend on the gauge choice and alter the hyperbolic structure of the system of partial differential equations compared to general relativity. Accordingly, the strong hyperbolicity of the system cannot be claimed yet, under the same assumptions as in general relativity. In the paper, we stress the differences compared with general relativity and state the main issues that should be tackled next, to draw a road map towards numerical bimetric relativity.

We report lensing magnifications, extinction, and time-delay estimates for the first resolved, multiply-imaged Type Ia supernova iPTF16geu, at $z = 0.409$, using $Hubble\,Space\,Telescope$ ($HST$) observations in combination with supporting ground-based data. Multi-band photometry of the resolved images provides unique information about the differential dimming due to dust in the lensing galaxy. Using $HST$ and Keck AO reference images taken after the SN faded, we obtain a total lensing magnification for iPTF16geu of $\mu = 67.8^{+2.6}_{-2.9}$, accounting for extinction in the host and lensing galaxy. As expected from the symmetry of the system, we measure very short time-delays for the three fainter images with respect to the brightest one: -0.23 $\pm$ 0.99, -1.43 $\pm$ 0.74 and 1.36 $\pm$ 1.07 days. Interestingly, we find large differences between the magnifications of the four supernova images, even after accounting for uncertainties in the extinction corrections: $\Delta m_1 = -3.88^{+0.07}_{-0.06}$, $\Delta m_2 = -2.99^{+0.09}_{-0.08}$, $\Delta m_3 = -2.19^{+0.14}_{-0.15}$ and $\Delta m_4 = -2.40^{+0.14}_{-0.12}$ mag, discrepant with model predictions suggesting similar image brightnesses. A possible explanation for the large differences is gravitational lensing by substructures, micro- or millilensing, in addition to the large scale lens causing the image separations. We find that the inferred magnification is insensitive to the assumptions about the dust properties in the host and lens galaxy.

Recently, Kenna-Allison et.al. claimed that bimetric gravity cannot give rise to a viable cosmological expansion history while at the same time being compatible with local gravity tests. In this note we review that claim and combine various results from the literature to provide several simple counter examples. We conclude that the results of Kenna-Allison et.al. cannot hold in general.

The recent detection of gravitational waves (GWs) and electromagnetic (EM) waves originating from the same source marks the start of a new multi-messenger era in astronomy. The arrival time difference between the GW and EM signal can be used to constrain differences in their propagation speed, and thus gravitational theories. We study to what extent a non-zero time delay can be explained by gravitational lensing when the line of sight to the source passes near a massive object. For galaxy scale lenses, this delay becomes relevant for GWs with frequencies between $10^{-6}$ and $10^{-9}$ Hz, sourced by super massive binary black-holes. In addition to GWs detectable by Pulsar Timing Arrays (PTAs), we expect to find also a unique and recognizable EM signal. We show that the delay between the GW and EM signal can be of the order of days to months; within reach of future observations. The effect may become important in future multi-messenger astronomy probing of gravitational propagation and interactions.

Constraints on models of the late time acceleration of the universe assume the cosmological principle of homogeneity and isotropy on large scales. However, small scale inhomogeneities can alter observational and dynamical relations, affecting the inferred cosmological parameters. For precision constraints on the properties of dark energy, it is important to assess the potential systematic effects arising from these inhomogeneities. In this study, we use the Type Ia supernova magnitude-redshift relation to constrain the inhomogeneities as described by the Dyer-Roeder distance relation and the effect they have on the dark energy equation of state ($w$), together with priors derived from the most recent results of the measurements of the power spectrum of the Cosmic Microwave Background and Baryon Acoustic Oscillations. We find that the parameter describing the inhomogeneities ($\eta$) is correlated with $w$. The best fit values $w = -0.933 \pm 0.065$ and $\eta = 0.61 \pm 0.37$ are consistent with homogeneity at $< 2 \sigma$ level. Assuming homogeneity ($\eta =1$), we find $w = -0.961 \pm 0.055$, indicating only a small change in $w$. For a time-dependent dark energy equation of state, $w_0 = -0.951 \pm 0.112$ and $w_a = 0.059 \pm 0.418$, to be compared with $w_0 = -0.983 \pm 0.127$ and $w_a = 0.07 \pm 0.432$ in the homogeneous case, which is also a very small change. Current data are not sufficient to constrain the fraction of dark matter (DM) in compact objects, $f_p$ at the 95$\%$ C.L., however at 68$\%$ C.L. $f_p < 0.73$. Future supernova surveys will improve the constraints on $\eta$, and $f_p$, by a factor of $\sim$ 10.

The $\Lambda$ Cold Dark Matter model ($\Lambda$CDM) represents the current standard model in cosmology. Within this, there is a tension between the value of the Hubble constant, $H_0$, inferred from local distance indicators and the angular scale of fluctuations in the Cosmic Microwave Background (CMB). We investigate whether the tension is significant enough to warrant new physics in the form of modifying or adding energy components to the standard cosmological model. We find that late time dark energy explanations are slightly disfavoured whereas a pre-CMB decoupling extra dark energy component has a marginally positive Bayesian evidence. A constant equation of state of the additional early energy density is constrained to 0.086$^{+0.04}_{-0.03}$. Although this value deviates significantly from 1/3, valid for dark radiation, the latter is not disfavoured based on the Bayesian evidence. If the tension persists, future estimates of $H_0$ at the 1$\%$ level will be able to decisively determine which of the proposed explanations is favoured.

Recent re-calibration of the Type Ia supernova (SNe~Ia) magnitude-redshift relation combined with cosmic microwave background (CMB) and baryon acoustic oscillation (BAO) data have provided excellent constraints on the standard cosmological model. Here, we examine particular classes of alternative cosmologies, motivated by various physical mechanisms, e.g. scalar fields, modified gravity and phase transitions to test their consistency with observations of SNe~Ia and the ratio of the angular diameter distances from the CMB and BAO. Using a model selection criterion for a relative comparison of the models (the Bayes Factor), we find moderate to strong evidence that the data prefer flat $\Lambda$ CDM over models invoking a thawing behaviour of the quintessence scalar field. However, some exotic models like the growing neutrino mass cosmology and vacuum metamorphosis still present acceptable evidence values. The bimetric gravity model with only the linear interaction term can be ruled out by the combination of SNe~Ia and CMB/BAO datasets whereas the model with linear and quadratic interaction terms has a comparable evidence value to standard $\Lambda$ CDM. Thawing models are found to have significantly poorer evidence compared to flat $\Lambda$ CDM cosmology under the assumption that the CMB compressed likelihood provides an adequate description for these non-standard cosmologies. We also present estimates for constraints from future data and find that geometric probes from oncoming surveys can put severe limits on non-standard cosmological models.

Despite its continued observational successes, there is a persistent (and growing) interest in extending cosmology beyond the standard model, $\Lambda$CDM. This is motivated by a range of apparently serious theoretical issues, involving such questions as the cosmological constant problem, the particle nature of dark matter, the validity of general relativity on large scales, the existence of anomalies in the CMB and on small scales, and the predictivity and testability of the inflationary paradigm. In this paper, we summarize the current status of $\Lambda$CDM as a physical theory, and review investigations into possible alternatives along a number of different lines, with a particular focus on highlighting the most promising directions. While the fundamental problems are proving reluctant to yield, the study of alternative cosmologies has led to considerable progress, with much more to come if hopes about forthcoming high-precision observations and new theoretical ideas are fulfilled.

We investigate Lema\^itre-Tolman-Bondi (LTB) models, whose early time evolution and bang time are homogeneous and the distance - redshift relation and local Hubble parameter are inherited from the $\Lambda$CDM model. We show that the obtained LTB models and the $\Lambda$CDM model predict different relative local expansion rates and that the Hubble functions of the models diverge increasingly with redshift. The LTB models show tension between low redshift baryon acoustic oscillation and supernova observations and including Lyman-$\alpha$ forest or cosmic microwave background observations only accentuates the better fit of the $\Lambda$CDM model compared to the LTB model. The result indicates that additional degrees of freedom are needed to explain the observations, for example by renouncing spherical symmetry, homogeneous bang time, negligible effects of pressure, or the early time homogeneity assumption.

There is a no-go theorem forbidding flat and closed FLRW solutions in massive gravity on a flat reference metric, while open solutions are unstable. Recently it was shown that this no-go theorem can be overcome if at least some matter couples to a hybrid metric composed of both the dynamical and the fixed reference metric. We show that this is not compatible with the standard description of cosmological sources in terms of effective perfect fluids, and the predictions of the theory become sensitive either to the detailed field-theoretical modelling of the matter content or to the presence of additional dark degrees of freedom. This is a serious practical complication. Furthermore, we demonstrate that viable cosmological background evolution with a perfect fluid appears to require the presence of fields with highly contrived properties. This could be improved if the equivalence principle is broken by coupling only some of the fields to the composite metric, but viable self-accelerating solutions due only to the massive graviton are difficult to obtain.

Our ability to study the most remote supernova explosions, crucial for the understanding of the evolution of the high-redshift universe and its expansion rate, is limited by the light collection capabilities of telescopes. However, nature offers unique opportunities to look beyond the range within reach of our unaided instruments thanks to the light-focusing power of massive galaxy clusters. Here we report on the discovery of one of the most distant supernovae ever found, at redshift, z=1.703. Due to a lensing magnification factor of 4.3\pm0.3, we are able to measure a lightcurve of the supernova, as well as spectroscopic features of the host galaxy with a precision comparable to what will otherwise only be possible with future generation telescopes.

Massive galaxy clusters at intermediate redshifts act as gravitational lenses that can magnify supernovae (SNe) occurring in background galaxies. We assess the possibility to use lensed SNe to put constraints on the mass models of galaxy clusters and the Hubble parameter at high redshift. Due to the standard candle nature of Type Ia supernovae (SNe Ia), observational information on the lensing magnification from an intervening galaxy cluster can be used to constrain the model for the cluster mass distribution. A statistical analysis using parametric cluster models was performed to investigate the possible improvements from lensed SNe Ia for the accurately modeled galaxy cluster A1689 and the less well constrained cluster A2204. Time delay measurements obtained from SNe lensed by accurately modeled galaxy clusters can be used to measure the Hubble parameter. For a survey of A1689 we estimate the expected rate of detectable SNe Ia and of multiply imaged SNe. The velocity dispersion and core radius of the main cluster potential show strong correlations with the predicted magnifications and can therefore be constrained by observations of SNe Ia in background galaxies. This technique proves especially powerful for galaxy clusters with only few known multiple image systems. The main uncertainty for measurements of the Hubble parameter from the time delay of strongly lensed SNe is due to cluster model uncertainties. For the extremely well modeled cluster A1689, a single time delay measurement could be used to determine the Hubble parameter with a precision of ~ 10%. We conclude that observations of SNe Ia behind galaxy clusters can be used to improve the mass modeling of the large scale component of galaxy clusters and thus the distribution of dark matter. Time delays from SNe strongly lensed by accurately modeled galaxy clusters can be used to measure the Hubble constant at high redshifts.

We present an analysis of peculiar velocities and their effect on supernova cosmology. In particular, we study (a) the corrections due to our own motion, (b) the effects of correlations in peculiar velocities induced by large-scale structure, and (c) uncertainties arising from a possible local under- or over-density. For all of these effects we present a case study of their impact on the cosmology derived by the Sloan Digital Sky Survey-II Supernova Survey (SDSS-II SN Survey). Correcting supernova redshifts for the CMB dipole slightly over-corrects nearby supernovae that share some of our local motion. We show that while neglecting the CMB dipole would cause a shift in the derived equation of state of Delta w ~ 0.04 (at fixed matter density) the additional local-motion correction is currently negligible (Delta w<0.01). We use a covariance-matrix approach to statistically account for correlated peculiar velocities. This down-weights nearby supernovae and effectively acts as a graduated version of the usual sharp low-redshift cut. Neglecting coherent velocities in the current sample causes a systematic shift of ~2% in the preferred value of w and will therefore have to be considered carefully when future surveys aim for percent-level accuracy. Finally, we perform n-body simulations to estimate the likely magnitude of any local density fluctuation (monopole) and estimate the impact as a function of the low-redshift cutoff. We see that for this aspect the low-z cutoff of z=0.02 is well-justified theoretically, but that living in a putative local density fluctuation leaves an indelible imprint on the magnitude-redshift relation.

Massive galaxy clusters at intermediate redshift can magnify the flux of distant background sources by several magnitudes and we exploit this effect to search for lensed distant supernovae that may otherwise be too faint to be detected. A supernova search was conducted at near infrared wavelengths using the ISAAC instrument at the VLT. The galaxy clusters Abell 1689, Abell 1835 and AC114 were observed at multiple epochs of 2 hours of exposure time, separated by a month. Image-subtraction techniques were used to search for transient objects with light curve properties consistent with supernovae, both in our new and archival ISAAC/VLT data. The limiting magnitude of the individual epochs was estimated by adding artificial stars to the subtracted images. Most of the epochs reach 90% detection efficiency at SZ(J) ~= 23.8-24.0 mag (Vega). Two transient objects, both in archival images of Abell 1689 and AC114, were detected. The transient in AC114 coincides - within the position uncertainty - with an X-ray source and is likely to be a variable AGN at the cluster redshift. The transient in Abell 1689 was found at SZ=23.24 mag, ~0.5 arcsec away from a galaxy with photometric redshift z=0.6 +/-0.15. The light curves and the colors of the transient are consistent with a reddened Type IIP supernova at redshift z=0.59 +/- 0.05. The lensing model of Abell 1689 predicts ~1.4 mag of magnification at the position of the transient, making it the most magnified supernova ever found and only the second supernova found behind a galaxy cluster. Our pilot survey has demonstrated the feasibility to find distant gravitationally magnified supernovae behind massive galaxy clusters. One likely supernova was found behind Abell 1689, in accordance with the expectations for this survey, as shown in an accompanying analysis paper.

We use the new SNe Ia discovered by the SDSS-II Supernova Survey together with additional supernova datasets as well as observations of the cosmic microwave background and baryon acoustic oscillations to constrain cosmological models. This complements the analysis presented by Kessler et al. in that we discuss and rank a number of the most popular non-standard cosmology scenarios. When this combined data-set is analyzed using the MLCS2k2 light-curve fitter, we find that more exotic models for cosmic acceleration provide a better fit to the data than the Lambda-CDM model. For example, the flat DGP model is ranked higher by our information criteria tests than the standard model. When the dataset is instead analyzed using the SALT-II light-curve fitter, the standard cosmological constant model fares best. Our investigation also includes inhomogeneous Lemaitre-Tolman-Bondi (LTB) models. While our LTB models can be made to fit the supernova data as well as any other model, the extra parameters they require are not supported by our information criteria analysis.