The subject of this paper is a quantification of the information content of cosmological probes of the large-scale structures, specifically of temperature and polarisation anisotropies in the cosmic microwave background, CMB-lensing, weak cosmic shear and galaxy clustering, in terms of Information theory measures like information entropies. We aim to establish relationships for Gaussian likelihoods, between conventional measures of statistical uncertainties and information entropies. Furthermore, we extend these studies to the computation of (Bayesian) evidences and the power of measurement to distinguish between competing models. We investigate in detail how cosmological data decreases information entropy by reducing statistical errors and by breaking degeneracies. In addition, we work out how tensions between data sets increase information entropy and quantify this effect in three examples: the discrepancy in $\Omega_m$ and $\sigma_8$ between the CMB and weak lensing, the role of intrinsic alignments in weak lensing data when attempting the dark energy equation of state parameters, and the famous $H_0$-tension between Cepheids in the Hubble keystone project and the cosmic microwave background as observed by Planck.

This paper is a pedagogical introduction to models of gravity and how to constrain them through cosmological observations. We focus on the Horndeski scalar-tensor theory and on the quantities that can be measured with a minimum of assumptions. Alternatives or extensions of General Relativity have been proposed ever since its early years. Because of Lovelock theorem, modifying gravity in four dimensions typically means adding new degrees of freedom. The simplest way is to include a scalar field coupled to the curvature tensor terms. The most general way of doing so without incurring in the Ostrogradski instability is the Horndeski Lagrangian and its extensions. Testing gravity means therefore, in its simplest term, testing the Horndeski Lagrangian. Since local gravity experiments can always be evaded by assuming some screening mechanism or that baryons are decoupled, or even that the effects of modified gravity are visible only at early times, we need to test gravity with cosmological observations in the late universe (large-scale structure) and in the early universe (cosmic microwave background). In this work we review the basic tools to test gravity at cosmological scales, focusing on model-independent measurements.

In this work, we use recent data on the Hubble expansion rate $H(z)$, the quantity $f\sigma_8(z)$ from redshift space distortions and the statistic $E_g$ from clustering and lensing observables to constrain in a model-independent way the linear anisotropic stress parameter $\eta$. This estimate is free of assumptions about initial conditions, bias, the abundance of dark matter and the background expansion. We denote this observable estimator as $\eta_{{\rm obs}}$. If $\eta_{{\rm obs}}$ turns out to be different from unity, it would imply either a modification of gravity or a non-perfect fluid form of dark energy clustering at sub-horizon scales. Using three different methods to reconstruct the underlying model from data, we report the value of $\eta_{{\rm obs}}$ at three redshift values, $z=0.29, 0.58, 0.86$. Using the method of polynomial regression, we find $\eta_{{\rm obs}}=0.57\pm1.05$, $\eta_{{\rm obs}}=0.48\pm0.96$, and $\eta_{{\rm obs}}=-0.11\pm3.21$, respectively. Assuming a constant $\eta_{{\rm obs}}$ in this range, we find $\eta_{{\rm obs}}=0.49\pm0.69$. We consider this method as our fiducial result, for reasons clarified in the text. The other two methods give for a constant anisotropic stress $\eta_{{\rm obs}}=0.15\pm0.27$ (binning) and $\eta_{{\rm obs}}=0.53 \pm 0.19$ (Gaussian Process). We find that all three estimates are compatible with each other within their $1\sigma$ error bars. While the polynomial regression method is compatible with standard gravity, the other two methods are in tension with it.

We introduce a methodology to test models with spatial variations of the fine-structure constant $\alpha$, based on the calculation of the angular power spectrum of these measurements. This methodology enables comparisons of observations and theoretical models through their predictions on the statistics of the $\alpha$ variation. Here we apply it to the case of symmetron models. We find no indications of deviations from the standard behavior, with current data providing an upper limit to the strength of the symmetron coupling to gravity ($\log{\beta^2}<-0.9$) when this is the only free parameter, and not able to constrain the model when also the symmetry breaking scale factor $a_{SSB}$ is free to vary.

Astrophysical tests of the stability of fundamental couplings are becoming an increasingly important probe of new physics. Motivated by the recent availability of new and stronger constraints we update previous works testing the consistency of measurements of the fine-structure constant $\alpha$ and the proton-to-electron mass ratio $\mu=m_p/m_e$ (mostly obtained in the optical/ultraviolet) with combined measurements of $\alpha$, $\mu$ and the proton gyromagnetic ratio $g_p$ (mostly in the radio band). We carry out a global analysis of all available data, including the 293 archival measurements of {\it Webb et al.} and 66 more recent dedicated measurements, and constraining both time and spatial variations. While nominally the full datasets show a slight statistical preference for variations of $\alpha$ and $\mu$ (at up to two standard deviations), we also find several inconsistencies between different sub-sets, likely due to hidden systematics and implying that these statistical preferences need to be taken with caution. The statistical evidence for a spatial dipole in the values of $\alpha$ is found at the 2.3 sigma level. Forthcoming studies with facilities such as ALMA and ESPRESSO should clarify these issues.

Recent work by Webb {\it et al.} has provided indications of spatial variations of the fine-structure constant, $\alpha$, at a level of a few parts per million. Using a dataset of 293 archival measurements, they further show that a dipole provides a statistically good fit to the data, a result subsequently confirmed by other authors. Here we show that a more recent dataset of dedicated measurements further constrains these variations: although there are only 10 such measurements, their uncertainties are considerably smaller. We find that a dipolar variation is still a good fit to the combined dataset, but the amplitude of such a dipole must be somewhat smaller: $8.1\pm1.7$ ppm for the full dataset, versus $9.4\pm2.2$ ppm for the Webb {\it et al.} data alone, both at the $68.3\%$ confidence level. Constraints on the direction on the sky of such a dipole are also significantly improved. On the other hand the data can't yet discriminate between a pure spatial dipole and one with an additional redshift dependence.

Astrophysical tests of the stability of fundamental couplings, such as the fine-structure constant $\alpha$, are a powerful probe of new physics. Recently these measurements, combined with local atomic clock tests and Type Ia supernova and Hubble parameter data, were used to constrain the simplest class of dynamical dark energy models where the same degree of freedom is assumed to provide both the dark energy and (through a dimensionless coupling, $\zeta$, to the electromagnetic sector) the $\alpha$ variation. One caveat of these analyses was that it was based on fiducial models where the dark energy equation of state was described by a single parameter (effectively its present day value, $w_0$). Here we relax this assumption and study broader dark energy model classes, including the Chevallier-Polarski-Linder and Early Dark Energy parametrizations. Even in these extended cases we find that the current data constrains the coupling $\zeta$ at the $10^{-6}$ level and $w_0$ to a few percent (marginalizing over other parameters), thus confirming the robustness of earlier analyses. On the other hand, the additional parameters are typically not well constrained. We also highlight the implications of our results for constraints on violations of the Weak Equivalence Principle and improvements to be expected from forthcoming measurements with high-resolution ultra-stable spectrographs.

Astrophysical tests of the stability of fundamental couplings, such as the fine-structure constant $\alpha$, are becoming an increasingly powerful probe of new physics. Here we discuss how these measurements, combined with local atomic clock tests and Type Ia supernova and Hubble parameter data, constrain the simplest class of dynamical dark energy models where the same degree of freedom is assumed to provide both the dark energy and (through a dimensionless coupling, $\zeta$, to the electromagnetic sector) the $\alpha$ variation. Specifically, current data tightly constrains a combination of $\zeta$ and the present dark energy equation of state $w_0$. Moreover, in these models the new degree of freedom inevitably couples to nucleons (through the $\alpha$ dependence of their masses) and leads to violations of the Weak Equivalence Principle. We obtain indirect bounds on the E\"otv\"os parameter $\eta$ that are typically stronger than the current direct ones. We discuss the model-dependence of our results and briefly comment on how the forthcoming generation of high-resolution ultra-stable spectrographs will enable significantly tighter constraints.

We use astrophysical and atomic clock tests of the stability of the fine-structure constant $\alpha$, together with Type Ia supernova and Hubble parameter data, to constrain the simplest class of dynamical dark energy models where the same degree of freedom is assumed to provide both the dark energy and (through a dimensionless coupling, $\zeta$, to the electromagnetic sector) the $\alpha$ variation. We show how current data tightly constrains a combination of $\zeta$ and the dark energy equation of state $w_0$. At the $95\%$ confidence level and marginalizing over $w_0$ we find $|\zeta|<5\times10^{-6}$, with the atomic clock tests dominating the constraints. The forthcoming generation of high-resolution ultra-stable spectrographs will enable significantly tighter constraints.