When Peccei-Quinn (PQ) symmetry breaking happens after inflation, the axion field takes random values in causally disconnected regions. This leads to fluctuations of order one in the axion energy density around the QCD epoch. These over-densities eventually decouple from the Hubble expansion and form so-called miniclusters. We present a semi-analytical method to calculate the average axion energy density, as well as the power spectrum, from the re-alignment mechanism in this scenario. Furthermore, we develop a modified Press & Schechter approach, suitable to describe the collapse of non-linear density fluctuations during radiation domination, which is relevant for the formation of axion miniclusters. It allows us to calculate the double differential distribution of gravitationally collapsed miniclusters as a function of their mass and size. For instance, assuming a PQ scale of $10^{11}$ GeV, minicluster masses range from about $5 \times 10^{-16}$ to $3 \times 10^{-13}$ solar masses and have sizes from about $4\times 10^4$ to $7\times 10^5$ km at the time they start to collapse.

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.

In this paper we study the phenomenology of stars and galaxies in massive bigravity. We give parameter conditions for the existence of viable star solutions when the radius of the star is much smaller than the Compton wavelength of the graviton. If these parameter conditions are not met, we constrain the ratio between the coupling constants of the two metrics, in order to give viable conditions for e.g. neutron stars. For galaxies, we put constraints on both the Compton wavelength of the graviton and the conformal factor and coupling constants of the two metrics. The relationship between black holes and stars, and whether the former can be formed from the latter, is discussed. We argue that the different asymptotic structure of stars and black holes makes it unlikely that black holes form from the gravitational collapse of stars in massive bigravity.

We investigate the observational consequences of scalar instabilities in bimetric theory, under the assumption that the Vainshtein mechanism restores general relativity within a certain distance from gravitational sources. We argue that early time instabilities have a negligible impact on observed structures. Assuming that the instabilities affect sub-horizon density fluctuations, we constrain the redshift, z_i, below which instabilities are ruled out. For the "minimal" beta_1-model, observational constraints are close to the theoretical expectations of z_i = 0.5, potentially allowing the model to be ruled in or out with a more detailed study, possibly including secondary cosmic microwave background constraints.

We study the integrated Sachs-Wolfe (ISW) effect in ghost-free, massive bigravity. We focus on the infinite-branch bigravity (IBB) model which exhibits viable cosmic expansion histories and stable linear perturbations, while the cosmological constant is set to zero and the late-time accelerated expansion of the Universe is due solely to the gravitational interaction terms. The ISW contribution to the CMB auto-correlation power spectrum is predicted, as well as the cross-correlation between the CMB temperature anisotropies and large-scale structure. We use ISW amplitudes as observed by the WMAP 9-year temperature data together with galaxy and AGN data provided by the WISE mission in order to compare the theoretical predictions to the observations. The ISW amplitudes in IBB are found to be larger than the corresponding ones in the standard LCDM model by roughly a factor of 1.5, but are still consistent with the observations.

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.

We study the cosmic expansion history of massive bigravity with a viable matter coupling which treats both metrics on equal footing. We derive the Friedmann equation for the effective metric through which matter couples to the two metrics, and study its solutions. For certain parameter choices, the background cosmology is identical to that of LCDM. More general parameters yield dynamical dark energy, which can still be in agreement with observations of the expansion history. We study specific parameter choices of interest, including minimal models, maximally-symmetric models, and a candidate partially-massless theory.

We derive dynamical and gravitational lensing properties of local sources in the Hassan-Rosen bimetric gravity theory. Observations of elliptical galaxies rule out values of the effective length-scale of the theory, in units of the Hubble radius, in the interval 10^-6 < lambda_g/r_H < 10^-3, unless the proportionality constant between the metrics at the background level is far from unity, in which case general relativity is effectively restored for local sources. In order to have background solutions resembling the concordance cosmological model, without fine-tuning of the parameters of the model, we are restricted to the upper interval, or lambda_g/r_H ~ 1. Except for a limited range of parameter values, the Hassan-Rosen theory is thus consistent with the observed lensing and dynamical properties of elliptical galaxies.

We perform cosmological perturbation theory in Hassan-Rosen bimetric gravity for general homogeneous and isotropic backgrounds. In the de Sitter approximation, we obtain decoupled sets of massless and massive scalar gravitational fluctuations. Matter perturbations then evolve like in Einstein gravity. We perturb the future de Sitter regime by the ratio of matter to dark energy, producing quasi-de Sitter space. In this more general setting the massive and massless fluctuations mix. We argue that in the quasi-de Sitter regime, the growth of structure in bimetric gravity differs from that of Einstein gravity.

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.

We investigate constraints on dark energy fluctuations using type Ia supernovae. If dark energy is not in the form of a cosmological constant, that is if the equation of state is not equal to -1, we expect not only temporal, but also spatial variations in the energy density. Such fluctuations would cause local variations in the universal expansion rate and directional dependences in the redshift-distance relation. We present a scheme for relating a power spectrum of dark energy fluctuations to an angular covariance function of standard candle magnitude fluctuations. The predictions for a phenomenological model of dark energy fluctuations are compared to observational data in the form of the measured angular covariance of Hubble diagram magnitude residuals for type Ia supernovae in the Union2 compilation. The observational result is consistent with zero dark energy fluctuations. However, due to the limitations in statistics, current data still allow for quite general dark energy fluctuations as long as they are in the linear regime.

Accretion of dark energy onto black holes will take place when dark energy is not a cosmological constant. It has been proposed that the time evolution of the mass of the black holes in binary systems due to dark energy accretion could be detectable by gravitational radiation. This would make it possible to use observations of black hole binaries to measure local dark energy properties, e.g., to determine the sign of 1+w where w is the dark energy equation of state. In this Letter we show that such measurements are unfeasible due to the low accretion rates.