We study ghost-free multimetric theories for $(N+1)$ tensor fields with a coupling to matter and maximal global symmetry group $S_N\times(Z_2)^N$. Their mass spectra contain a massless mode, the graviton, and $N$ massive spin-2 modes. One of the massive modes is distinct by being the heaviest, the remaining $(N-1)$ massive modes are simply identical copies of each other. All relevant physics can therefore be understood from the case $N=2$. Focussing on this case, we compute the full perturbative action up to cubic order and derive several features that hold to all orders in perturbation theory. The lighter massive mode does not couple to matter and neither of the massive modes decay into massless gravitons. We propose the lighter massive particle as a candidate for dark matter and investigate its phenomenology in the parameter region where the matter coupling is dominated by the massless graviton. The relic density of massive spin-2 can originate from a freeze-in mechanism or from gravitational particle production, giving rise to two different dark matter scenarios. The allowed parameter regions are very different from those in scenarios with only one massive spin-2 field and more accessible to experiments.

We provide further details on a recent proposal addressing the nature of the dark sectors in cosmology and demonstrate that all current observations related to Dark Matter can be explained by the presence of a heavy spin-2 particle. Massive spin-2 fields and their gravitational interactions are uniquely described by ghost-free bimetric theory, which is a minimal and natural extension of General Relativity. In this setup, the largeness of the physical Planck mass is naturally related to extremely weak couplings of the heavy spin-2 field to baryonic matter and therefore explains the absence of signals in experiments dedicated to Dark Matter searches. It also ensures the phenomenological viability of our model as we confirm by comparing it with cosmological and local tests of gravity. At the same time, the spin-2 field possesses standard gravitational interactions and it decays universally into all Standard Model fields but not into massless gravitons. Matching the measured DM abundance together with the requirement of stability constrains the spin-2 mass to be in the 1 to 100 TeV range.

Observational evidence for the existence of Dark Matter is limited to its gravitational effects. The extensive program for dedicated searches has yielded null results so far, challenging the most popular models. Here we propose that this is the case because the very existence of cold Dark Matter is a manifestation of gravity itself. The consistent bimetric theory of gravity, the only known ghost-free extension of General Relativity involving a massless and a massive spin-2 field, automatically contains a perfect Dark Matter candidate. We demonstrate that the massive spin-2 particle can be heavy, stable on cosmological scales, and that it interacts with matter only through a gravitational type of coupling. Remarkably, these features persist in the same region of parameter space where bimetric theory satisfies the current gravity tests. We show that the observed Dark Matter abundance can be generated via freeze-in and suggest possible particle physics and gravitational signatures of our bimetric Dark Matter model.

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 a large class of infinitesimal, but fully nonlinear in the field, transformations of the Galileon and search for extended symmetries. The transformations involve powers of the coordinates $x$ and the field $\pi$ up to any finite order $N$. Up to quadratic order the structure of these symmetry transformations is the unique generalisation of both the infinitesimal version of the standard Galileon shift symmetry as well as a recently discovered infinitesimal extension of this symmetry. The only higher-order extensions of this symmetry we recover are (`Galileon dual' versions of) symmetries of the standard kinetic term.

We obtain the general cosmological evolution equations for a classically consistent theory of bimetric gravity. Their analytic solutions are demonstrated to generically allow for a cosmic evolution starting out from a matter dominated FLRW universe while relaxing towards a de Sitter (anti-de Sitter) phase at late cosmic time. In particular, we examine a subclass of models which contain solutions that are able to reproduce the expansion history of the cosmic concordance model inspite of the nonlinear couplings of the two metrics. This is demonstrated explicitly by fitting these models to observational data from Type Ia supernovae, Cosmic Microwave Background and Baryon Acoustic Oscillations.

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.

Powerful gravitational telescopes in the form of massive galaxy clusters can be used to enhance the light collecting power over a limited field of view by about an order of magnitude in flux. This effect is exploited here to increase the depth of a survey for lensed supernovae at near-IR wavelengths. A pilot SN search program conducted with the ISAAC camera at VLT is presented. Lensed galaxies behind the massive clusters A1689, A1835 and AC114 were observed for a total of 20 hours split into 2, 3 and 4 epochs respectively, separated by approximately one month to a limiting magnitude J<24 (Vega). Image subtractions including another 20 hours worth of archival ISAAC/VLT data were used to search for transients with lightcurve properties consistent with redshifted supernovae, both in the new and reference data. The feasibility of finding lensed supernovae in our survey was investigated using synthetic lightcurves of supernovae and several models of the volumetric Type Ia and core-collapse supernova rates as a function of redshift. We also estimate the number of supernova discoveries expected from the inferred star formation rate in the observed galaxies. The methods consistently predict a Poisson mean value for the expected number of SNe in the survey between N_SN=0.8 and 1.6 for all supernova types, evenly distributed between core collapse and Type Ia SN. One transient object was found behind A1689, 0.5" from a galaxy with photometric redshift z_gal=0.6 +- 0.15. The lightcurve and colors of the transient are consistent with being a reddened Type IIP SN at z_SN=0.59. The lensing model predicts 1.4 magnitudes of magnification at the location of the transient, without which this object would not have been detected in the near-IR ground based search described in this paper (unlensed magnitude J~25). (abridged)