The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will soon start to scan thousands of square degrees of the northern extragalactic sky with a unique set of $56$ optical filters from a dedicated $2.55$m telescope, JST, at the Javalambre Astrophysical Observatory. Before the arrival of the final instrument (a 1.2 Gpixels, 4.2deg$^2$ field-of-view camera), the JST was equipped with an interim camera (JPAS-Pathfinder), composed of one CCD with a 0.3deg$^2$ field-of-view and resolution of 0.23 arcsec pixel$^{-1}$. To demonstrate the scientific potential of J-PAS, with the JPAS-Pathfinder camera we carried out a survey on the AEGIS field (along the Extended Groth Strip), dubbed miniJPAS. We observed a total of $\sim 1$ deg$^2$, with the $56$ J-PAS filters, which include $54$ narrow band (NB, $\rm{FWHM} \sim 145$Angstrom) and two broader filters extending to the UV and the near-infrared, complemented by the $u,g,r,i$ SDSS broad band (BB) filters. In this paper we present the miniJPAS data set, the details of the catalogues and data access, and illustrate the scientific potential of our multi-band data. The data surpass the target depths originally planned for J-PAS, reaching $\rm{mag}_{\rm {AB}}$ between $\sim 22$ and $23.5$ for the NB filters and up to $24$ for the BB filters ($5\sigma$ in a $3$~arcsec aperture). The miniJPAS primary catalogue contains more than $64,000$ sources extracted in the $r$ detection band with forced photometry in all other bands. We estimate the catalogue to be complete up to $r=23.6$ for point-like sources and up to $r=22.7$ for extended sources. Photometric redshifts reach subpercent precision for all sources up to $r=22.5$, and a precision of $\sim 0.3$% for about half of the sample. (Abridged)

We review some of our recent results about the Radial Acceleration Relation (RAR) and its interpretation as either a fundamental or an emergent law. The former interpretation is in agreement with a class of modified gravity theories that dismiss the need for dark matter in galaxies (MOND in particular). Our most recent analysis, which includes refinements on the priors and the Bayesian test for compatibility between the posteriors, confirms that the hypothesis of a fundamental RAR is rejected at more than 5$\sigma$ from the very same data that was used to infer the RAR.

The Radial Acceleration Relation (RAR) shows a strong correlation between two accelerations associated to galaxy rotation curves. The relation between these accelerations is given by a nonlinear function which depends on an acceleration scale $a_\dagger$. Some have interpreted this as an evidence for a gravity model, such as Modified Newtonian Dynamics (MOND), which posits a fundamental acceleration scale $a_0$ common to all the galaxies. However, it was later shown, using Bayesian inference, that this seems not to be the case: the $a_0$ credible intervals for individual galaxies were not found to be compatible among themselves. This type of test is a fundamental test for MOND as a theory for gravity, since it directly evaluates its basic assumption and this using the data that most favor MOND: galaxy rotation curves. Here we improve upon the previous analyses by introducing a more robust method to assess the compatibility between the credible intervals, in particular without Gaussian approximations. We directly estimate, using a Monte Carlo simulation, that the existence of a fundamental acceleration is incompatible with the data at more than $5\sigma$. We also consider quality cuts in order to show that our results are robust against outliers. In conclusion, the new analysis further supports the claim that the acceleration scale found in the RAR is an emergent quantity.

The cosmic distance ladder is the succession of techniques by which it is possible to determine distances to astronomical objects. Here, we present a new method to build the cosmic distance ladder, going from local astrophysical measurements to the CMB. Instead of relying on high-redshift cosmography in order to model the luminosity-distance relation and calibrate supernovae with BAO, we exploit directly the distance-duality relation $d_L = (1+z)^2 d_A$---valid if photon number is conserved and gravity is described by a metric theory. The advantage is that the results will not depend on the parametrization of the luminosity-distance relation at $z>0.15$: no model is adopted in order to calibrate BAO with supernovae. This method yields local measurements of the Hubble constant and deceleration parameter. Furthermore, it can directly assess the impact of BAO observations on the strong 4--5$\sigma$ tension between local and global $H_0$. Using the latest supernova, BAO and CMB observations, we found a consistently low value of $q_0$ and strong inconsistency between angular-only BAO constraints and anisotropic BAO measurements, which are, or not, in agreement with CMB depending on the kind of analysis (see Table~$4$). We conclude that, in order to understand the reasons behind the $H_0$ crisis, a first step should be clarifying the tension between angular and perpendicular anisotropic BAO as this will help understanding if new physics is required at the pre-recombination epoch or/and during the dark energy era.

We revisit the dark degeneracy that arises from the Einstein equations relating geometry to the total cosmic substratum but not resolving its individual components separately. We establish the explicit conditions for the dark degeneracy in the fluid description of the dark sector. At the background level, this degeneracy can be formally understood in terms of a unified dark sector Equation of State (EoS) that depends both on the dynamical nature of the dark energy (DE) as well as on its interaction with the pressureless dark matter. For linear perturbations, the degeneracy arises for specified DE pressure perturbations (or sound speed, equivalently) and DE anisotropic stress. Specializing to the degeneracy between non-interacting dynamical DE and interacting vacuum DE models, we perform a parameter estimation analysis for a range of dynamical DE parametrizations, where for illustration we explicitly break the degeneracy at the linear level by adopting a luminal sound speed for both scenarios. We conduct this analysis using cosmological background data alone and in combination with Planck measurements of the cosmic microwave background radiation. We find that although the overall phenomenology between the dynamical DE and interacting approaches is similar, there are some intriguing differences. In particular, there is an ambiguity in the strength of constraints on $\Omega_{m0}$ and $\sigma_8$, which are considerably weakened for interacting vacuum DE, indicating that the dark degeneracy can change the significance of tensions in cosmological parameters inferred from different data sets.

Cameron {\it et al.}~2019 (hereafter C19) recommends a more cautious and rigorous approach to statistical analysis in astronomy. We welcome this particular side of their communication as it helps stimulating the effort towards the adoption of better statistical methods in galaxy rotation curves, an effort to which we contributed with Rodrigues {\it et al.}~2018 (hereafter R18). Indeed, R18 was the first work that, in order to conclude on the universality of the acceleration scale $a_0$, studied the posterior distributions on $a_0$ of a large set of galaxies. As C19 agrees, the credible intervals were found within the Bayesian framework, that is, the marginalized posteriors on $a_0$ were found using Bayes' theorem to update the priors in light of the observational data; this process was done without introducing any approximation. Considering R18, C19 also remarks that: i) better methods to select the nuisance parameters and the corresponding priors could be used; ii) a quality cut based on $\chi^2$ values should not be used, and iii) the compatibility of the posteriors should be assessed in a more robust way. In the following, after first clarifying the context of our work, we address these criticisms.

We investigate local environmental effects from dark matter (DM) on thermonuclear supernovae (SNe Ia) using publicly available archival data of 224 low-redshift events, in an attempt to shed light on the SN Ia progenitor systems. SNe Ia are explosions of carbon-oxygen (CO) white dwarfs (WDs) that have recently been shown to explode at sub-Chandrasekhar masses; the ignition mechanism remains, however, unknown. Recently, it has been shown that both weakly interacting massive particles (WIMPs) and macroscopic DM candidates such as primordial black holes (PBHs) are capable of triggering the ignition. Here, we present a method to estimate the DM density and velocity dispersion in the vicinity of SN Ia events and nearby WDs; we argue that (i) WIMP ignition is highly unlikely, and that (ii) DM in the form of PBHs distributed according to a (quasi-) log-normal mass distribution with peak $\log_{10}(m_0/1$g$)=24.9\pm 0.9$ and width $\sigma= 3.3\pm 1.0$ is consistent with SN Ia data, the nearby population of WDs and roughly consistent with other constraints from the literature.

This white paper summarizes the activities of the Brazilian community concerning dark matter physics and highlights the importance of financial support to Brazilian groups that are deeply involved in experimental endeavours. The flagships of the Brazilian dark matter program are the Cherenkov Telescope Array, DARKSIDE, SBN and LHC experiments, but we emphasize that smaller experiments such as DAMIC and CONNIE constitute important probes to dark sectors as well and should receive special attention. Small experimental projects showing the potential to probe new regions of parameter space of dark matter models are encouraged. On the theoretical and phenomenological side, some groups are devoted to astrophysical aspects such as the dark matter density profile while others explore the signature of dark matter models at colliders, direct and indirect detection experiments. In summary, the Brazilian dark matter community that was born not long ago has grown tremendously in the past years and now plays an important role in the hunt for a dark matter particle.

The determination of the Hubble constant $H_0$ from the Cosmic Microwave Background by the Planck Collaboration [Aghanim et al. 2018] is in tension at $4.2\sigma$ with respect to the local determination of $H_0$ by the SH0ES collaboration [Reid et al. 2019]. Here, we improve upon the local determination, which fixes the deceleration parameter to the standard $\Lambda$CDM model value of $q_0=-0.55$, that is, uses information from observations beyond the local universe. First, we derive the effective calibration prior on the absolute magnitude $M_B$ of Supernovae Ia, which can be used in cosmological analyses in order to avoid the double counting of low-redshift supernovae. We find $M_B = -19.2334 \pm 0.0404$ mag. Then, we use the above $M_B$ prior in order to obtain a determination of the local $H_0$ which only uses local observations and only assumes the cosmological principle, that is, large-scale homogeneity and isotropy. This is achieved by adopting an uninformative flat prior for $q_0$ in the cosmographic expansion of the luminosity distance. We use the latest Pantheon sample and find $H_0= 75.35 \pm 1.68 \text{ km s}^{-1} {\rm Mpc}^{-1}$, which features a 2.2% uncertainty, close to the 1.9% error obtained by the SH0ES Collaboration. Our determination is at the higher tension of $4.5\sigma$ with the latest results from the Planck Collaboration that assume the $\Lambda$CDM model. Furthermore, we also constrain the deceleration parameter to $q_0= -1.08 \pm 0.29$, which disagrees with Planck at the $1.9\sigma$ level. These estimations only use supernovae in the redshift range $0.023\le z\le 0.15$.

We study the minimally and non-minimally coupled scalar field models as possible alternatives for dark energy, the mysterious energy component that is driving the accelerated expansion of the universe. After discussing the dynamics at both the background and perturbation level, we confront the two models with the latest cosmological data. After obtaining updated constraints on their parameters we perform model selection using the basic information criteria. We found that the $\Lambda$CDM model is strongly favored when the local determination of the Hubble constant is not considered and that this statement is weakened once local $H_0$ is included in the analysis. We calculate the parameter combination $S_8=\sigma_8\sqrt{\Omega_{m}/0.3}$ and show the decrement of the tension with respect to the Planck results in the case of minimally and non-minimally coupled scalar field models. Finally, for the coupling constant between DE and gravity, we obtain the constraint $\xi\simeq -0.06^{+0.19}_{-0.19}$, approaching the one from solar system tests $|\xi| \lesssim 10^{-2}$ and comparable to the conformal value $\xi=1/6$ at $1\sigma$ uncertainty.

Using almost one million galaxies from the final Data Release 12 of the SDSS's Baryon Oscillation Spectroscopic Survey, we have obtained, albeit with low significance, a first model-independent determination of the radial BAO peak with 9% error: $\Delta z_{\rm BAO}(z_{\rm eff}=0.51)= 0.0456 \pm 0.0042$. In order to obtain this measurement, the radial correlation function was computed in 7,700 angular pixels, from which mean correlation function and covariance matrix were obtained, making the analysis completely model independent. This novel method of obtaining the covariance matrix was validated via the comparison with 500 BOSS mock catalogs. This $\Delta z_{\rm BAO}$ determination can be used to constrain the background expansion of exotic models for which the assumptions adopted in the standard analysis cannot be satisfied. Future galaxy catalogs from J-PAS, DESI and Euclid are expected to significantly increase the quality and significance of model-independent determinations of the BAO peak, possibly determined at various redshift and angular positions. We stress that it is imperative to test the standard paradigm in a model-independent way in order to test its foundations, maximize the extraction of information from the data, and look for clues regarding the poorly understood dark energy and dark matter.

Despite the observational success of the standard model of cosmology, present-day observations do not tightly constrain the nature of dark matter and dark energy and modifications to the theory of general relativity. Here, we will discuss some of the ongoing and upcoming surveys that will revolutionize our understanding of the dark sector.

Since there is no known symmetry in Nature that prevents a non-minimal coupling between the dark energy (DE) and cold dark matter (CDM) components, such a possibility constitutes an alternative to standard cosmology, with its theoretical and observational consequences being of great interest. In this paper we propose a new null test on the standard evolution of the dark sector based on the time dependence of the ratio between the CDM and DE energy densities which, in the standard $\Lambda$CDM scenario, scales necessarily as $a^{-3}$. We use the latest measurements of type Ia supernovae, cosmic chronometers and angular baryonic acoustic oscillations to reconstruct the expansion history using model-independent Machine Learning techniques, namely, the Linear Model formalism and Gaussian Processes. We find that while the standard evolution is consistent with the data at $3\sigma$ level, some deviations from the $\Lambda$CDM model are found at low redshifts, which may be associated with the current tension between local and global determinations of $H_0$.

This correspondence answers the criticisms raised in McGaugh et al 2018 and Kroupa et al 2018.

The current $3.8 \sigma$ tension between local (arXiv:1804.10655) and global (arXiv:1605.02985) measurements of $H_0$ cannot be fully explained by the concordance $\Lambda$CDM model. It could be produced by unknown systematics or by physics beyond the standard model. On the other hand, it is well known that linear perturbation theory predicts a cosmic variance on the Hubble parameter $H_0$, which leads to systematic errors on its local determination. Here, we study how including in the likelihood the cosmic variance on $H_0$ affects statistical inference. In particular we consider the $\gamma$CDM, $w$CDM and $\gamma w$CDM parametric extensions of the standard model, which we constrain with the latest CMB, BAO, SNe Ia, RSD and $H_0$ data. We learn two important lessons. First, the systematic error from cosmic variance is - independently of the model - approximately $\sigma_{\text{cv}}\approx 0.88$ km s$^{-1}$ Mpc$^{-1}$ (1.2\% $H_0^{\text{loc}}$) when considering the redshift range $0.0233 \le z \le 0.15$, which is relative to the main analysis of (arXiv:1804.10655), and $\sigma_{\text{cv}}\approx 1.5$ km s$^{-1}$ Mpc$^{-1}$ (2.1\% $H_0^{\text{loc}}$) when considering the wider redshift range $0.01 \le z \le 0.15$. Although $\sigma_{\text{cv}}$ affects the total error budget on local $H_0$, it does not significantly alleviate the tension which remains at $\approx 3 \sigma$. Second, cosmic variance, besides shifting the constraints, can change the results of model selection: much of the statistical advantage of non-standard models is to alleviate the now-reduced tension. We conclude that, when constraining non-standard models it is important to include the cosmic variance on $H_0$ if one wants to use the local determination of the Hubble constant by Riess et al. (arXiv:1804.10655). Doing the contrary could potentially bias the conclusions.

The Radial Acceleration Relation confirms that a nontrivial acceleration scale can be found in the average internal dynamics of galaxies. The existence of such a scale is not obvious as far as the standard cosmological model is concerned, and it has been interpreted as a possible sign of modified gravity. The implications could be profound: it could in principle explain galactic dynamics without large amounts of yet-undetected dark matter and address issues that the standard cosmological model faces at galactic scales. Here, we consider 193 disk galaxies from the SPARC and THINGS databases and, using Bayesian inference, we show that the probability of existence of a fundamental acceleration that is common to all the galaxies is essentially zero: the $p$-value is smaller than $10^{-20}$ or, equivalently, the null hypothesis is rejected at more than 10$\sigma$. We conclude that the acceleration scale unveiled by the Radial Acceleration Relation is of emergent nature, possibly caused by a complex interplay between baryons and dark matter. In particular, the MOND theory, or any other theory that behaves like it at galactic scales, is ruled out as a fundamental theory for galaxies at more than 10$\sigma$.

We test both the FLRW geometry and $\Lambda$CDM cosmology in a model independent way by reconstructing the Hubble function $H(z)$, the comoving distance $D(z)$ and the growth of structure $f\sigma_8(z)$ using the most recent data available. We use the linear model formalism in order to optimally reconstruct the latter cosmological functions, together with their derivatives and integrals. We then evaluate four of the null tests available in literature: $Om_{1}$ by Sahni et al., $Om_{2}$ by Zunckel \& Clarkson, $Ok$ by Clarkson et al., and $ns$ by Nesseris \& Sapone. For all the four tests we find agreement, within the errors, with the standard cosmological model.

Within the standard paradigm, dark energy is taken as a homogeneous fluid that drives the accelerated expansion of the universe and does not contribute to the mass of collapsed objects such as galaxies and galaxy clusters. The abundance of galaxy clusters -- measured through a variety of channels -- has been extensively used to constrain the normalization of the power spectrum: it is an important probe as it allows us to test if the standard $\Lambda$CDM model can indeed accurately describe the evolution of structures across billions of years. It is then quite significant that the Planck satellite has detected, via the Sunyaev-Zel'dovich effect, less clusters than expected according to the primary CMB anisotropies. One of the simplest generalizations that could reconcile these observations is to consider models in which dark energy is allowed to cluster, i.e., allowing its sound speed to vary. In this case, however, the standard methods to compute the abundance of galaxy clusters need to be adapted to account for the contributions of dark energy. In particular, we examine the case of clustering dark energy -- a dark energy fluid with negligible sound speed -- with a redshift-dependent equation of state. We carefully study how the halo mass function is modified in this scenario, highlighting corrections that have not been considered before in the literature. We address modifications in the growth function, collapse threshold, virialization densities and also changes in the comoving scale of collapse and mass function normalization. Our results show that clustering dark energy can impact halo abundances at the level of 10\%--30\%, depending on the halo mass, and that cluster counts are modified by about 30\% at a redshift of unity.

We develop and apply new techniques in order to uncover galaxy rotation curves (RC) systematics. Considering that an ideal dark matter (DM) profile should yield RCs that have no bias towards any particular radius, we find that the Burkert DM profile satisfies the test, while the Navarro-Frenk-While (NFW) profile has a tendency of better fitting the region between one and two disc scale lengths than the inner disc scale length region. Our sample indicates that this behaviour happens to more than 75% of the galaxies fitted with an NFW halo. Also, this tendency does not weaken by considering "large" galaxies, for instance those with $M_*\gtrsim 10^{10} M_\odot$. Besides the tests on the homogeneity of the fits, we also use a sample of 62 galaxies of diverse types to perform tests on the quality of the overall fit of each galaxy, and to search for correlations with stellar mass, gas mass and the disc scale length. In particular, we find that only 13 galaxies are better fitted by the NFW halo; and that even for the galaxies with $M_* \gtrsim 10^{10} M_\odot$ the Burkert profile either fits as good as, or better than, the NFW profile. This result is relevant since different baryonic effects important for the smaller galaxies, like supernova feedback and dynamical friction from baryonic clumps, indicate that at such large stellar masses the NFW profile should be preferred over the Burkert profile. Hence, our results either suggest a new baryonic effect or a change of the dark matter physics.

During the process of structure formation in the universe matter is converted into radiation through a variety of processes such as light from stars, infrared radiation from cosmic dust and gravitational waves from binary black holes/neutron stars and supernova explosions. The production of this astrophysical radiation background (ARB) could affect the expansion rate of the universe and the growth of perturbations. Here, we aim at understanding to which level one can constraint the ARB using future cosmological observations. We model the energy transfer from matter to radiation through an effective interaction between matter and astrophysical radiation. Using future supernova data from LSST and growth-rate data from Euclid we find that the ARB density parameter is constrained, at the 95% confidence level, to be $\Omega_{ar_0}<0.008$. Estimates of the energy density produced by well-known astrophysical processes give roughly $\Omega_{ar_0}\sim 10^{-5}$. Therefore, we conclude that cosmological observations will only be able to constrain exotic or not-well understood sources of radiation.

The abundances of dark matter halos in the universe are described by the halo mass function (HMF). It enters most cosmological analyses and parametrizes how the linear growth of primordial perturbations is connected to these abundances. Interestingly, this connection can be made approximately cosmology independent. This made it possible to map in detail its near-universal behavior through large-scale simulations. However, such simulations may suffer from systematic effects, especially if baryonic physics is included. In this paper we ask how well observations can constrain directly the HMF. The observables we consider are galaxy cluster number counts, galaxy cluster power spectrum and lensing of type Ia supernovae. Our results show that DES is capable of putting the first meaningful constraints on the HMF, while both Euclid and J-PAS can give stronger constraints, comparable to the ones from state-of-the-art simulations. We also find that an independent measurement of cluster masses is even more important for measuring the HMF than for constraining the cosmological parameters, and can vastly improve the determination of the halo mass function. Measuring the HMF could thus be used to cross-check simulations and their implementation of baryon physics. It could even, if deviations cannot be accounted for, hint at new physics.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

We propose that dark energy in the form of a scalar field could effectively couple to dark matter inhomogeneities. Through this coupling energy could be transferred to/from the scalar field, which could possibly enter an accelerated regime. Though phenomenological, this scenario is interesting as it provides a natural trigger for the onset of the acceleration of the universe, since dark energy starts driving the expansion of the universe when matter inhomogeneities become sufficiently strong. Here we study a possible realization of this idea by coupling dark energy to dark matter via the linear growth function of matter perturbations. The numerical results show that it is indeed possible to obtain a viable cosmology with the expected series of radiation, matter and dark-energy dominated eras. In particular, the current density of dark energy is given by the value of the coupling parameters rather than by very special initial conditions for the scalar field. In other words, this model - unlike standard models of cosmic late acceleration - does not suffer from the so-called "coincidence problem" and its related fine tuning of initial conditions.

We study the validity of the Newtonian description of cosmological perturbations using the Lemaitre model, an exact spherically symmetric solution of Einstein's equation. This problem has been investigated in the past for the case of a dust fluid. Here, we extend the previous analysis to the more general case of a fluid with non-negligible pressure, and, for the numerical examples, we consider the case of radiation (P=\rho/3). We find that, even when the density contrast has a nonlinear amplitude, the Newtonian description of the cosmological perturbations using the gravitational potential \psi and the curvature potential \phi is valid as long as we consider sub-horizon inhomogeneities. However, the relation \psi+\phi={\cal O}(\phi^2), which holds for the case of a dust fluid, is not valid for a relativistic fluid and effective anisotropic stress is generated. This demonstrates the usefulness of the Lemaitre model which allows us to study in an exact nonlinear fashion the onset of anisotropic stress in fluids with non-negligible pressure. We show that this happens when the characteristic scale of the inhomogeneity is smaller than the sound horizon and that the deviation is caused by the nonlinear effect of the fluid's fast motion. We also find that \psi+\phi= \max[{\cal O}(\phi^2),{\cal O}(c_s^2\phi \, \delta)] for an inhomogeneity with density contrast \delta whose characteristic scale is smaller than the sound horizon, unless w is close to -1, where w and c_s are the equation of state parameter and the sound speed of the fluid, respectively. On the other hand, we expect \psi+\phi={\cal O}(\phi^2) to hold for an inhomogeneity whose characteristic scale is larger than the sound horizon, unless the amplitude of the inhomogeneity is large and w is close to -1.

A recently proposed technique allows one to constrain both the background and perturbation cosmological parameters through the distribution function of supernova Ia apparent magnitudes. Here we extend this technique to alternative cosmological scenarios, in which the growth of structure does not follow the $\Lambda$CDM prescription. We apply the method first to the supernova data provided by the JLA catalog combined with all the current independent redshift distortion data and with low-redshift cluster data from Chandra and show that although the supernovae alone are not very constraining, they help in reducing the confidence regions. Then we apply our method to future data from LSST and from a survey that approximates the Euclid satellite mission. In this case we show that the combined data are nicely complementary and can constrain the normalization $\sigma_8$ and the growth rate index $\gamma$ to within $0.6\%$ and $7\%$, respectively. In particular, the LSST supernova catalog is forecast to give the constraint $\gamma (\sigma_8/0.83)^{6.7} = 0.55 \pm 0.1$. We also report on constraints relative to a step-wise parametrization of the growth rate of structures. These results show that supernova lensing serves as a good cross-check on the measurement of perturbation parameters from more standard techniques.

A range of Bayesian tools has become widely used in cosmological data treatment and parameter inference (see Kunz, Bassett & Hlozek (2007), Trotta (2008), Amendola, Marra & Quartin (2013)). With increasingly big datasets and higher precision, tools that enable us to further enhance the accuracy of our measurements gain importance. Here we present an approach based on internal robustness, introduced in Amendola, Marra & Quartin (2013) and adopted in Heneka, Marra & Amendola (2014), to identify biased subsets of data and hidden correlation in a model independent way.

We argue that there is an intrinsic noise on measurements of the equation of state parameter $w=p/\rho$ from large-scale structure around us. The presence of the large-scale structure leads to an ambiguity in the definition of the background universe and thus there is a maximal precision with which we can determine the equation of state of dark energy. To study the uncertainty due to local structure, we model density perturbations stemming from a standard inflationary power spectrum by means of the exact Lema\^{i}tre-Tolman-Bondi solution of Einstein's equation, and show that the usual distribution of matter inhomogeneities in a $\Lambda$CDM cosmology causes a variation of $w$ -- as inferred from distance measures -- of several percent. As we observe only one universe, or equivalently because of the cosmic variance, this uncertainty is systematic in nature.

The Multi-coupled Dark Energy (McDE) scenario has been recently proposed as a specific example of a cosmological model characterized by a non-standard physics of the dark sector of the universe that nevertheless gives an expansion history which does not significantly differ from the one of the standard $\Lambda $CDM model. In this work, we present the first constraints on the McDE scenario obtained by comparing the predicted evolution of linear density perturbations with a large compilation of recent data sets for the growth rate $f\sigma_{8}$, including 6dFGS, LRG, BOSS, WiggleZ and VIPERS. Confirming qualitative expectations, growth rate data provide much tighter bounds on the model parameters as compared to the extremely loose bounds that can be obtained when only the background expansion history is considered. In particular, the $95\%$ confidence level on the coupling strength $|\beta |$ is reduced from $|\beta |\leq 83$ (background constraints only) to $|\beta |\leq 0.88$ (background and linear perturbation constraints). We also investigate how these constraints further improve when using data from future wide-field surveys such as supernova data from LSST and growth rate data from Euclid-type missions. In this case the $95\%$ confidence level on the coupling further reduce to $|\beta |\leq 0.85$. Such constraints are in any case still consistent with a scalar fifth-force of gravitational strength, and we foresee that tighter bounds might be possibly obtained from the investigation of nonlinear structure formation in McDE cosmologies.[Abridged]

Soon the number of type Ia supernova (SN) measurements should exceed 100,000. Understanding the effect of weak lensing by matter structures on the supernova brightness will then be more important than ever. Although SN lensing is usually seen as a source of systematic noise, we will show that it can be in fact turned into signal. More precisely, the non-Gaussianity introduced by lensing in the SN Hubble diagram dispersion depends rather sensitively on the amplitude sigma8 of the matter power spectrum. By exploiting this relation, we are able to predict constraints on sigma8 of 7% (3%) for a catalog of 100,000 (500,000) SNe of average magnitude error 0.12 without having to assume that such intrinsic dispersion is known a priori. The intrinsic dispersion has been assumed to be Gaussian; possible intrinsic non-Gaussianities in the dataset (due to the SN themselves and/or to other transients) could be potentially dealt with by means of additional nuisance parameters describing higher moments of the intrinsic dispersion distribution function. This method is independent of and complementary to the standard methods based on CMB, cosmic shear or cluster abundance observables.

The use of advanced statistical analysis tools is crucial in order to improve cosmological parameter estimates via removal of systematic errors and identification of previously unaccounted for cosmological signals. Here we demonstrate the application of a new fully-Bayesian method, the internal robustness formalism, to scan for systematics and new signals in the recent supernova Ia Union compilations. Our analysis is tailored to maximize chances of detecting the anomalous subsets by means of a variety of sorting algorithms. We analyse supernova Ia distance moduli for effects depending on angular separation, redshift, surveys and hemispherical directions. The data have proven to be robust within 2 sigma, giving an independent confirmation of successful removal of systematics-contaminated supernovae. Hints of new cosmology, as for example the anisotropies reported by Planck, do not seem to be reflected in the supernova Ia data.

We present a new programme for placing constraints on radial inhomogeneity in a dark-energy dominated universe. We introduce a new measure to quantify violations of the Copernican principle. Any violation of this principle would interfere with our interpretation of any dark-energy evolution. In particular, we find that current observations place reasonably tight constraints on possible late-time violations of the Copernican principle: the allowed area in the parameter space of amplitude and scale of a spherical inhomogeneity around the observer has to be reduced by a factor of three so as to confirm the Copernican principle. Then, by marginalizing over possible radial inhomogeneity we provide the first constraints on the cosmological constant which are free of the homogeneity prior prevalent in cosmology.

A perfect irrotational fluid with the equation of state of dust, Irrotational Dark Matter (IDM), is incapable of virializing and instead forms a cosmoskeleton of filaments with supermassive black holes at the joints. This stark difference from the standard cold dark matter (CDM) scenario arises because IDM must exhibit potential flow at all times, preventing shell-crossing from occurring. This scenario is applicable to general non-oscillating scalar-field theories with a small sound speed. Our model of combined IDM and CDM components thereby provides a solution to the problem of forming the observed billion-solar-mass black holes at redshifts of six and higher. In particular, as a result of the reduced vortical flow, the growth of the black holes is expected to be more rapid at later times as compared to the standard scenario.

With the availability of thousands of type Ia supernovae in the near future the magnitude scatter induced by lensing will become a major issue as it affects parameter estimation. Current N-body simulations are too time consuming to be integrated in the likelihood analyses used for estimating the cosmological parameters. In this paper we show that in the weak lensing regime a statistical numerical approximation produces accurate results orders of magnitude faster. We write down simple fits to the second, third and fourth central moments of the lensing magnification probability distribution as a function of redshift, of the power spectrum normalization and of the present-day matter density. We also improve upon existing models of lensing variance and show that a shifted lognormal distribution fits well the numerical one. These fits can be easily employed in cosmological likelihood analyses. Moreover, our theoretical predictions make it possible to invert the problem and begin using supernovae lensing to constrain the cosmological parameters.

The persisting consistency of ever more accurate observational data with the predictions of the standard LCDM cosmological model puts severe constraints on possible alternative scenarios, but still does not shed any light on the fundamental nature of the cosmic dark sector.As large deviations from a LCDM cosmology are ruled out by data, the path to detect possible features of alternative models goes necessarily through the definition of cosmological scenarios that leave almost unaffected the background and -- to a lesser extent -- the linear perturbations evolution of the universe. In this context,the Multi-coupled DE (McDE) model was proposed by Baldi 2012 as a particular realization of an interacting Dark Energy field characterized by an effective screening mechanism capable of suppressing the effects of the coupling at the background and linear perturbation level. In the present paper, for the first time, we challenge the McDE scenario through a direct comparison with real data, in particular with the luminosity distance of Type Ia supernovae. By studying the existence and stability conditions of the critical points of the associated background dynamical system, we select only the cosmologically consistent solutions, and confront their background expansion history with data. Confirming previous qualitative results, the McDE scenario appears to be fully consistent with the adopted sample of Type Ia supernovae, even for coupling values corresponding to an associated scalar fifth-force about four orders of magnitude stronger than standard gravity. Our analysis demonstrates the effectiveness of the McDE background screening, and shows some new non-trivial asymptotic solutions for the future evolution of the universe. Our results show how the background expansion history might be highly insensitive to the fundamental nature and to the internal complexity of the dark sector. [Abridged]

There is an approximately 9% discrepancy, corresponding to 2.4sigma, between two independent constraints on the expansion rate of the universe: one indirectly arising from the cosmic microwave background and baryon acoustic oscillations, and one more directly obtained from local measurements of the relation between redshifts and distances to sources. We argue that by taking into account the local gravitational potential at the position of the observer this tension - strengthened by the recent Planck results - is partially relieved and the concordance of the standard model of cosmology increased. We estimate that measurements of the local Hubble constant are subject to a cosmic variance of about 2.4% (limiting the local sample to redshifts z>0.010) or 1.3% (limiting it to z>0.023), a more significant correction than that taken into account already. Nonetheless, we show that one would need a very rare fluctuation to fully explain the offset in the Hubble rates. If this tension is further strengthened, a cosmology beyond the standard model may prove necessary.

We find that if we live at the center of an inhomogeneity with total density contrast of roughly 0.1, dark energy is not a cosmological constant at 95% confidence level. Observational constraints on the equation of state of dark energy, w, depend strongly on the local matter density around the observer. We model the local inhomogeneity with an exact spherically symmetric solution which features a pressureless matter component and a dark-energy fluid with constant equation of state and negligible sound speed, that reaches a homogeneous solution at finite radius. We fit this model to observations of the local expansion rate, distant supernovae and the cosmic microwave background. We conclude that the possible uncertainty from large-scale structure has to be taken into account if one wants to progress towards not just precision but also accurate cosmology.

A great deal of effort is currently being devoted to understanding, estimating and removing systematic errors in cosmological data. In the particular case of type Ia supernovae, systematics are starting to dominate the error budget. Here we propose a Bayesian tool for carrying out a systematic search for systematic contamination. This serves as an extension to the standard goodness-of-fit tests and allows not only to cross-check raw or processed data for the presence of systematics but also to pin-point the data that are most likely contaminated. We successfully test our tool with mock catalogues and conclude that the Union2.1 data do not possess a significant amount of systematics. Finally, we show that if one includes in Union2.1 the supernovae that originally failed the quality cuts, our tool signals the presence of systematics at over 3.8-sigma confidence level.

We present the exact equations governing the dynamics of a spherically-symmetric inhomogeneous model with n decoupled and non-comoving perfect fluids. Thanks to the use of physically meaningful quantities we write the set of 3+2n equations in a concise and transparent way. The n perfect fluids can have general equations of state, thus making the model extremely flexible to study a large variety of cosmological and astrophysical problems. As applications we consider a model sourced by two non-comoving dust components and a cosmological constant, and a model featuring dust and a dark energy component with negligible speed of sound.

The halo model (HM) describes the inhomogeneous universe as a collection of halos. The full nonlinear power spectrum of the universe is well approximated by the HM, whose prediction can be easily computed without lengthy numerical simulations. This makes the HM a useful tool in cosmology. Here we explore the lensing properties of the HM by use of the stochastic gravitational lensing (sGL) method. We obtain for the case of point sources exact and simple integral expressions for the expected value and variance of the lensing convergence, which encode detailed information about the internal halo properties. In particular a wide array of observational biases can be easily incorporated and the dependence of lensing on cosmology is properly taken into account. This simple setup should be useful for a quick calculation of the power spectrum and the related lensing observables, which can play an important role in the extraction of cosmological parameters from, e.g., SNe observations. Finally, we discuss the probability distribution function of the HM which encodes more information than the first two moments and can more strongly constrain the large-scale structures of the universe. To check the accuracy of our modelling we compare our predictions to the results from the Millennium Simulation.

It has been proposed that the observed dark energy can be explained away by the effect of large-scale nonlinear inhomogeneities. In the present paper we discuss how observations constrain cosmological models featuring large voids. We start by considering Copernican models, in which the observer is not occupying a special position and homogeneity is preserved on a very large scale. We show how these models, at least in their current realizations, are constrained to give small, but perhaps not negligible in certain contexts, corrections to the cosmological observables. We then examine non-Copernican models, in which the observer is close to the center of a very large void. These models can give large corrections to the observables which mimic an accelerated FLRW model. We carefully discuss the main observables and tests able to exclude them.

We revise and extend the stochastic approach to cumulative weak lensing (hereafter the sGL method) first introduced in Ref. [1]. Here we include a realistic halo mass function and density profiles to model the distribution of mass between and within galaxies, galaxy groups and galaxy clusters. We also introduce a modeling of the filamentary large-scale structures and a method to embed halos into these structures. We show that the sGL method naturally reproduces the weak lensing results for the Millennium Simulation. The strength of the sGL method is that a numerical code based on it can compute the lensing probability distribution function for a given inhomogeneous model universe in a few seconds. This makes it a useful tool to study how lensing depends on cosmological parameters and its impact on observations. The method can also be used to simulate the effect of a wide array of systematic biases on the observable PDF. As an example we show how simple selection effects may reduce the variance of observed PDF, which could possibly mask opposite effects from very large scale structures. We also show how a JDEM-like survey could constrain the lensing PDF relative to a given cosmological model. The updated turboGL code is available at turboGL.org.

We directly compare the concordance LCDM model to the inhomogeneous matter-only alternative represented by LTB void models. To achieve a "democratic" confrontation we explore LLTB models with non-vanishing cosmological constant and perform a global likelihood analysis in the parameter space of cosmological constant and void radius. In our analysis we carefully consider SNe, Hubble constant, CMB and BAO measurements, marginalizing over spectral index, age of the universe and background curvature. We find that the LCDM model is not the only possibility compatible with the observations, and that a matter-only void model is a viable alternative to the concordance model only if the BAO constraints are relaxed. Moreover, we will show that the areas of the parameter space which give a good fit to the observations are always disconnected with the result that a small local void does not significantly affect the parameter extraction for LCDM models.

We reanalyze the supernovae data from the Union Compilation including the weak lensing effects caused by inhomogeneities. We compute the lensing probability distribution function for each background solution described by the parameters Omega_M, Omega_L and w in the presence of inhomogeneities, approximately modeled with a single-mass population of halos. We then perform a likelihood analysis in the space of FLRW-parameters and compare our results with the standard approach. We find that the inclusion of lensing can move the best-fit model significantly towards the cosmic concordance of the flat LCDM model, improving the agreement with the constraints coming from the cosmic microwave background and baryon acoustic oscillations.

We study the weak gravitational lensing effects caused by a stochastic distribution of dark matter halos. We develop a simple approach to calculate the magnification probability distribution function which allows us to easily compute the magnitude bias and dispersion for an arbitrary data sample and a given universe model. As an application we consider the effects of single-mass large-scale cosmic inhomogeneities to the SNe magnitude-redshift relation, and conclude that such structures could bias the PDF enough to affect the extraction of cosmological parameters from the limited size of present-day SNe data samples. We also release turboGL, a simple and very fast (<= 1s) Mathematica code based on the method here presented.

We study the impact of cosmic inhomogeneities on the interpretation of SNe observations. We build an inhomogeneous universe model that can confront supernova data and yet is reasonably well compatible with the Copernican Principle. Our model combines a relatively small local void, that gives apparent acceleration at low redshifts, with a meatball model that gives sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different effects of voids on the data, allow the model to mimic the concordance model.

We study the impact of cosmic inhomogeneities on the interpretation of observations. We build an inhomogeneous universe model without dark energy that can confront supernova data and yet is reasonably well compatible with the Copernican Principle. Our model combines a relatively small local void, that gives apparent acceleration at low redshifts, with a meatball model that gives sizeable lensing (dimming) at high redshifts. Together these two elements, which focus on different effects of voids on the data, allow the model to mimic the concordance model.

The cosmological backreaction proposal, which attempts to account for observations without a primary dark energy source in the stress-energy tensor, has been developed and discussed by means of different approaches. Here, we focus on the concept of cosmological background solutions in order to develop a framework to study different backreaction proposals.

It has been argued that the spacetime of our universe can be accurately described by a perturbed conformal Newtonian metric, and hence even large density inhomogeneities in a dust universe can not change the observables predicted by the homogeneous dust model. In this paper we study a spherically symmetric dust model and illustrate conditions under which large spatial variations in the expansion rate can invalidate the argument.

By means of a toy Swiss-cheese cosmological model we will discuss how to set up and carry out in a physically meaningful way the idea of back-reaction, according to which dark energy could be an effective source. We will follow two distinct approaches. One is focused on how cosmological observables are affected by inhomogeneities, while the other is focused on a theoretical description of the inhomogeneous universe by means of a mean-field description.

This thesis is mainly about how to set up and carry out in a physically meaningful way the idea of back-reaction, according to which dark energy could be an effective source. There are, broadly speaking, two distinct approaches. One is focused on how cosmological observables are affected by inhomogeneities, while the other is focused on a theoretical description of the inhomogeneous universe by means of a mean-field description. Both approaches, however, share the idea of smoothing out inhomogeneities. We developed this duality in the interpretation of the back-reaction by means of toy models based on the Lemaitre-Tolman-Bondi solution of Einstein's equations. In particular we focused on voids expanding faster than the background solution.

We analyze a toy swiss-cheese cosmological model to study the averaging problem. In our model, the cheese is the EdS model and the holes are constructed from a LTB solution. We study the propagation of photons in the swiss-cheese model, and find a phenomenological homogeneous model to describe observables. Following a fitting procedure based on light-cone averages, we find that the the expansion scalar is unaffected by the inhomogeneities. This is because of spherical symmetry. However, the light-cone average of the density as a function of redshift is affected by inhomogeneities. The effect arises because, as the universe evolves, a photon spends more and more time in the (large) voids than in the (thin) high-density structures. The phenomenological homogeneous model describing the light-cone average of the density is similar to the concordance model. Although the sole source in the swiss-cheese model is matter, the phenomenological homogeneous model behaves as if it has a dark-energy component. Finally, we study how the equation of state of the phenomenological model depends on the size of the inhomogeneities, and find that the equation-of-state parameters w_0 and w_a follow a power-law dependence with a scaling exponent equal to unity. That is, the equation of state depends linearly on the distance the photon travels through voids. We conclude that within our toy model, the holes must have a present size of about 250 Mpc to be able to mimic the concordance model.

Photon geodesics are calculated in a swiss-cheese model, where the cheese is made of the usual Friedmann-Robertson-Walker solution and the holes are constructed from a Lemaitre-Tolman-Bondi solution of Einstein's equations. The observables on which we focus are the changes in the redshift, in the angular-diameter--distance relation, in the luminosity-distance--redshift relation, and in the corresponding distance modulus. We find that redshift effects are suppressed when the hole is small because of a compensation effect acting on the scale of half a hole resulting from the special case of spherical symmetry. However, we find interesting effects in the calculation of the angular distance: strong evolution of the inhomogeneities (as in the approach to caustic formation) causes the photon path to deviate from that of the FRW case. Therefore, the inhomogeneities are able to partly mimic the effects of a dark-energy component. Our results also suggest that the nonlinear effects of caustic formation in cold dark matter models may lead to interesting effects on photon trajectories.

A general model for the cosmological evolution of the fine structure constant \alpha driven by a typical Quintessence scenario is presented. We consider a coupling between the Quintessence scalar \phi and the electromagnetic kinetic term F^2, given by a general function B_F(\phi). We study the dependence of the cosmological \Delta\alpha(t) upon the functional form of B_F and discuss the constraints imposed by the data. We find that different cosmological histories for \Delta\alpha(t) are possible within the avaliable constraints. We also find that Quasar absorption spectra evidence for a time variation of alpha, if confirmed, is not incompatible with Oklo and meteorites limits.