We investigate the possibility of dark matter being a pure geometrical effect, rather than a particle or a compact object, by exploring a specific modified gravity model: mimetic dark matter. We present an alternative formulation of the theory, closer to the standard cosmological perturbation theory framework. We make manifest the presence of arbitrary parameters and extra functions, both at background level and at first order in perturbation theory. We present the full set of independent equations of motion for this model, and we discuss the amount of tuning needed to match predictions of the theory to actual data. By using the matter power spectrum and cosmic microwave background angular power spectra as benchmark observables, we explicitly show that since there is no natural mechanism to generate adiabatic initial conditions in this specific model, extra fine-tuning is required. We modify the publicly available Boltzmann code \texttt{CLASS} to make accurate predictions for the observables in mimetic dark matter. Our modified version of \texttt{CLASS} is available on GitHub. We have used mimetic dark matter as an illustration of how much one is allowed to change the initial conditions before contradicting observations when modifying the laws of gravity as described by General Relativity but we point out that modifying gravity without providing a natural mechanism to generate adiabatic initial conditions will always lead to highly fine-tuned models.

The concept of blind analysis, a key procedure to remove the human-based systematic error called confirmation bias, has long been an integral part of data analysis in many research areas. In cosmology, blind analysis is recently making its entrance, as the field progresses into a fully fledged high-precision science. The credibility, reliability and robustness of results from future sky-surveys will dramatically increase if the effect of confirmation bias is kept under control by using an appropriate blinding procedure. Here, we present a catalog-level blinding scheme for galaxy clustering data apt to be used in future spectroscopic galaxy surveys. We shift the individual galaxy positions along the line of sight based on 1) a geometric shift mimicking the Alcock-Paczynski effect and 2) a perturbative shift akin to redshift-space distortions. This procedure has several advantages. After combining the two steps above, it is almost impossible to accidentally unblind. The procedure induces a shift in cosmological parameters without changing the galaxies' angular positions, hence without interfering with the effects of angular systematics. Since the method is applied at catalog level, there is no need to adopt several blinding schemes tuned to different summary statistics, likelihood choices or types of analyses. By testing the method on mock catalogs and the BOSS DR12 catalog we demonstrate its performance in blinding galaxy clustering data for relevant cosmological parameters sensitive to the background expansion rate and the growth rate of structures.

In the era of precision cosmology, establishing the correct magnitude of statistical errors in cosmological parameters is of crucial importance. However, widely used approximations in galaxy surveys analyses can lead to parameter uncertainties that are grossly mis-estimated. These approximations can be introduced at three different levels: in the form of the likelihood, in the theoretical modelling of the observable and in the numerical computation of the observable. Their consequences are important both in data analysis through e.g., Markov Chain Monte Carlo parameter inference, and when survey instrument and strategy are designed and their constraining power on cosmological parameters is forecasted, for instance using Fisher matrix analyses. In this work, considering the galaxy angular power spectrum as the target observable, we report one example for each of such three categories and derive their impact on the estimated parameters uncertainties and degeneracies. We show that there are cases where these commonly used approximations lead, perhaps counter-intuitively, to unacceptably large mis-estimates of parameters errors and correlations. Furthermore, we stress how these approximations might even spoil the benefits of the nascent multi-tracer and multi-messenger cosmology. Hence we recommend that the type of analysis presented here should be repeated for every approximation adopted in survey design or data analysis, to quantify how it may affect the results. To this aim, we have developed \texttt{Multi\_CLASS}, a new extension of \texttt{CLASS} that includes the angular power spectrum for multiple (galaxy and other tracers such as gravitational waves) populations. The public release of \texttt{Multi\_CLASS} is associated with this paper.

Cosmological parameter estimation from forthcoming experiments promise to reach much greater precision than current constraints. As statistical errors shrink, the required control over systematic errors increases. Therefore, models or approximations that were sufficiently accurate so far, may introduce significant systematic biases in the parameter best-fit values and jeopardize the robustness of cosmological analyses. We present a general expression to estimate a priori the systematic error introduced in parameter inference due to the use of insufficiently good approximations in the computation of the observable of interest or the assumption of an incorrect underlying model. Although this methodology can be applied to measurements of any scientific field, we illustrate its power by studying the effect of modeling the angular galaxy power spectrum incorrectly. We also introduce Multi_CLASS, a new, public modification of the Boltzmann code CLASS, which includes the possibility to compute angular cross-power spectra for two different tracers. We find that significant biases in most of the cosmological parameters are introduced if one assumes the Limber approximation or neglects lensing magnification in modern galaxy survey analyses, and the effect is in general larger for the multi-tracer case, especially for the parameter controlling primordial non-Gaussianity of the local type, $f_{\rm NL}$.

We present a new estimator for the cross-correlation signal between line intensity maps to probe the Epoch of Reionization. The proposed estimator is the hitherto neglected antisymmetric component of the cross-correlation, under the exchange of line-of-sight positions. We consider the cross-correlation between HI and CO fluctuations, and forecast the improvement in precision on reionization parameters when the antisymmetric contribution is accounted for. As a way to break the degeneracy between astrophysics and cosmology in the intensity mapping power spectrum, we study the ratio between the antisymmetric and symmetric components. While our results depend on the highly uncertain astrophysical modelling, we show that in most standard scenarios including the antisymmetric contribution as a complementary probe can lead to a significant gain in information.

Baryon acoustic oscillations (BAO) provide a robust standard ruler, and can be used to constrain the expansion history of the Universe at low redshift. Standard BAO analyses return a model-independent measurement of the expansion rate and the angular diameter distance as function of redshift, normalized by the sound horizon at radiation drag. However, this methodology relies on anisotropic distance distortions of a fixed, pre-computed template (obtained in a given fiducial cosmology) in order to fit the observations. Therefore, it may be possible that extensions to the consensus $\Lambda$CDM add contributions to the BAO feature that cannot be captured by the template fitting. We perform mock BAO fits to power spectra computed assuming cosmological models which modify the growth of perturbations prior to recombination in order to test the robustness of the standard BAO analysis. We find no significant bias in the BAO analysis for the models under study ($\Lambda$CDM with a free effective number of relativistic species, early dark energy, and a model with interactions between neutrinos and a fraction of the dark matter), even for cases which do not provide a good fit to Planck measurements of the CMB power spectra. This result supports the use of the standard BAO analysis and its measurements to perform cosmological parameter inference and to constrain exotic models. In addition, we provide a methodology to reproduce our study for different models and surveys, as well as discuss different options to handle eventual biases in the BAO measurements.

Constraints on dark matter annihilation or decay offer unique insights into the nature of dark matter. We illustrate how surveys dedicated to detect the highly redshifted 21 cm signal from the dark ages will offer a new window into properties of particle dark matter. The 21 cm intensity mapping signal and its fluctuations are sensitive to energy injection from annihilation or decay of long-lived particles in a way that is complementary to other probes. We present forecasted constraints from forthcoming and next-generation radio surveys. We show that, while SKA might be capable of a detection for some cases, the most promising opportunity to detect dark matter in the 21 cm intensity mapping signal is with a futuristic radio array on the lunar far-side, with the potential to detect a signal many orders of magnitude weaker than current or maximal constraints from other probes.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

This paper outlines the science case for line-intensity mapping with a space-borne instrument targeting the sub-millimeter (microwaves) to the far-infrared (FIR) wavelength range. Our goal is to observe and characterize the large-scale structure in the Universe from present times to the high redshift Epoch of Reionization. This is essential to constrain the cosmology of our Universe and form a better understanding of various mechanisms that drive galaxy formation and evolution. We argue that the proposed frequency range would make it possible to probe important metal cooling lines such as [CII] up to very high redshift as well as a large number of rotational lines of the CO molecule. These can be used to trace molecular gas and dust evolution and constrain the buildup in both the cosmic star formation rate density and the cosmic infrared background (CIB). Moreover, surveys at the highest frequencies will detect FIR lines which are used as diagnostics of galaxies and AGN. Tomography of these lines over a wide redshift range will enable invaluable measurements of the cosmic expansion history at epochs inaccessible to other methods, competitive constraints on the parameters of the standard model of cosmology, and numerous tests of dark matter, dark energy, modified gravity and inflation. To reach these goals, large-scale structure must be mapped over a wide range in frequency to trace its time evolution over a reasonable fraction of the volume of the observable Universe. In addition, the surveyed area needs to be very large to beat cosmic variance and to probe the largest scales where its easier to separate the astrophysical and cosmological contributions to the observed signal. Only, a space-borne mission can properly meet these requirements.

Line-intensity mapping (LIM) of emission form star-forming galaxies can be used to measure the baryon acoustic oscillation (BAO) scale as far back as the epoch of reionization. This provides a standard cosmic ruler to constrain the expansion rate of the Universe at redshifts which cannot be directly probed otherwise. In light of growing tension between measurements of the current expansion rate using the local distance ladder and those inferred from the cosmic microwave background, extending the constraints on the expansion history to bridge between the late and early Universe is of paramount importance. Using a newly derived methodology to robustly extract cosmological information from LIM, which minimizes the inherent degeneracy with unknown astrophysics, we show that present and future experiments can gradually improve the measurement precision of the expansion rate history, ultimately reaching percent-level constraints on the BAO scale. Specifically, we provide detailed forecasts for the SPHEREx satellite, which will target the H$\alpha$ and Lyman-$\alpha$ lines, and for the ground-based COMAP instrument---as well as a future stage-3 experiment---that will target the CO rotational lines. Besides weighing in on the so-called Hubble tension, reliable LIM cosmic rulers can enable wide-ranging tests of dark matter, dark energy and modified gravity.

Line-intensity mapping (LIM) provides a promising way to probe cosmology, reionization and galaxy evolution. However, its sensitivity to cosmology and astrophysics at the same time is also a nuisance. Here we develop a comprehensive framework for modelling the LIM power spectrum, which includes redshift space distortions and the Alcock-Paczynski effect. We then identify and isolate degeneracies with astrophysics so that they can be marginalized over. We study the gains of using the multipole expansion of the anisotropic power spectrum, providing an accurate analytic expression for their covariance, and find a 10%-60% increase in the precision of the baryon acoustic oscillation scale measurements when including the hexadecapole in the analysis. We discuss different observational strategies when targeting other cosmological parameters, such as the sum of neutrino masses or primordial non-Gaussianity, finding that fewer and wider bins are typically more optimal. Overall, our formalism facilitates an optimal extraction of cosmological constraints robust to astrophysics.

The expansion of the Universe is understood to have accelerated during two epochs: in its very first moments during a period of Inflation and much more recently, at $z < 1$, when Dark Energy is hypothesized to drive cosmic acceleration. The undiscovered mechanisms behind these two epochs represent some of the most important open problems in fundamental physics. The large cosmological volume at $2 < z < 5$, together with the ability to efficiently target high-$z$ galaxies with known techniques, enables large gains in the study of Inflation and Dark Energy. A future spectroscopic survey can test the Gaussianity of the initial conditions up to a factor of ~50 better than our current bounds, crossing the crucial theoretical threshold of $\sigma(f_{NL}^{\rm local})$ of order unity that separates single field and multi-field models. Simultaneously, it can measure the fraction of Dark Energy at the percent level up to $z = 5$, thus serving as an unprecedented test of the standard model and opening up a tremendous discovery space.

HST/STIS spectra of the small clumps and filaments closest to the central star in VY CMa reveal that the very strong K I emission and TiO and VO molecular emission, long thought to form in a dusty circumstellar shell, actually originate in a few small clumps 100's of AU from the star. The K I lines are 10 to 20 times stronger in these nearest ejecta than on the star. The observations also confirm VO as a circumstellar molecule. In this letter we discuss the spectra of the features, their motions and ages, and the identification of the molecular emission. The strength of the atomic and molecular features in the small clumps present an astrophysical problem for the excitation process. We show that the clumps must have a nearly clear line of sight to the star's radiation.

Our current understanding of the Universe is established through the pristine measurements of structure in the cosmic microwave background (CMB) and the distribution and shapes of galaxies tracing the large scale structure (LSS) of the Universe. One key ingredient that underlies cosmological observables is that the field that sources the observed structure is assumed to be initially Gaussian with high precision. Nevertheless, a minimal deviation from Gaussianityis perhaps the most robust theoretical prediction of models that explain the observed Universe; itis necessarily present even in the simplest scenarios. In addition, most inflationary models produce far higher levels of non-Gaussianity. Since non-Gaussianity directly probes the dynamics in the early Universe, a detection would present a monumental discovery in cosmology, providing clues about physics at energy scales as high as the GUT scale.

The LIGO discoveries have rekindled suggestions that primordial black holes (BHs) may constitute part to all of the dark matter (DM) in the Universe. Such suggestions came from 1) the observed merger rate of the BHs, 2) their unusual masses, 3) their low/zero spins, and 4) also from the independently uncovered cosmic infrared background (CIB) fluctuations signal of high amplitude and coherence with unresolved cosmic X-ray background (CXB). Here we summarize the prospects to resolve this important issue with electromagnetic observations using the instruments and tools expected in the 2020's. These prospects appear promising to make significant, and potentially critical, advances. We demonstrate that in the next decade, new space- and ground-borne electromagnetic instruments, combined with concurrent theoretical efforts, should shed critical light on the long-considered link between primordial BHs and DM. Specifically the new data and methodologies under this program will involve: I) Probing with high precision the spatial spectrum of source-subtracted CIB with Euclid and WFIRST, and its coherence with unresolved cosmic X-ray background using eROSITA and Athena, II) Advanced searches for microlensing of Galactic stars by the intervening Galactic Halo BHs with OGLE, Gaia, LSST and WFIRST, III) Supernovae (SNe) lensing in the upcoming surveys with WFIRST, LSST and also potentially with Euclid and JWST, IV) Advanced theoretical work to understand the details of PBH accretion and evolution and their influence on cosmic microwave background (CMB) anisotropies in light of the next generation CMB experiments, V) Better new samples and theoretical understanding involving stability and properties of ultra faint dwarf galaxies, pulsar timing, and cosmological quasar lensing.

The Evolutionary Map of the Universe (EMU) is an all-sky survey in radio-continuum which uses the Australian SKA Pathfinder (ASKAP). Using galaxy angular power spectrum and the integrated Sachs-Wolfe effect, we study the potential of EMU to constrain models beyond $\Lambda$CDM (i.e., local primordial non-Gaussianity, dynamical dark energy, spatial curvature and deviations from general relativity), for different design sensitivities. We also include a multi-tracer analysis, distinguishing between star-forming galaxies and galaxies with an active galactic nucleus, to further improve EMU's potential. We find that EMU could measure the dark energy equation of state parameters around 35\% more precisely than existing constraints, and that the constraints on $f_{\rm NL}$ and modified gravity parameters will improve up to a factor $\sim2$ with respect to Planck and redshift space distortions measurements. With this work we demonstrate the promising potential of EMU to contribute to our understanding of the Universe.

Scattering and decay processes of thermal bath particles involving heavy leptons can dump hot axions in the primordial plasma around the QCD phase transition. We compute their relic density, parameterized by an effective number $\Delta N_{\rm eff}$ of additional neutrinos. For couplings allowed by current bounds, production via scattering yields $\Delta N_{\rm eff} \lesssim 0.6$ and $\Delta N_{\rm eff} \lesssim 0.2$ for the cases of muon and tau, respectively. Flavor violating tau decays to a lighter lepton plus an axion give $\Delta N_{\rm eff} \lesssim 0.3$. Such values of $\Delta N_{\rm eff}$ can alleviate the tension between the direct local measurement of the Hubble constant $H_0$ and the inferred value from observations of the Cosmic Microwave Background, assuming $\Lambda$CDM. We analyze present cosmological data from the Planck collaboration and baryon acoustic oscillations with priors given in terms of the axion-lepton couplings. For axions coupled to muons, the tension can be alleviated below the 3$\sigma$ level. Future experiments will measure $\Delta N_{\rm eff}$ with higher precision, providing an axion discovery channel and probing the role of hot axions in the $H_0$ tension.

We present a detailed overview of the cosmological surveys that will be carried out with Phase 1 of the Square Kilometre Array (SKA1), and the science that they will enable. We highlight three main surveys: a medium-deep continuum weak lensing and low-redshift spectroscopic HI galaxy survey over 5,000 sqdeg; a wide and deep continuum galaxy and HI intensity mapping survey over 20,000 sqdeg from z = 0.35 - 3; and a deep, high-redshift HI intensity mapping survey over 100 sqdeg from z = 3 - 6. Taken together, these surveys will achieve an array of important scientific goals: measuring the equation of state of dark energy out to z ~ 3 with percent-level precision measurements of the cosmic expansion rate; constraining possible deviations from General Relativity on cosmological scales by measuring the growth rate of structure through multiple independent methods; mapping the structure of the Universe on the largest accessible scales, thus constraining fundamental properties such as isotropy, homogeneity, and non-Gaussianity; and measuring the HI density and bias out to z = 6. These surveys will also provide highly complementary clustering and weak lensing measurements that have independent systematic uncertainties to those of optical surveys like LSST and Euclid, leading to a multitude of synergies that can improve constraints significantly beyond what optical or radio surveys can achieve on their own. This document, the 2018 Red Book, provides reference technical specifications, cosmological parameter forecasts, and an overview of relevant systematic effects for the three key surveys, and will be regularly updated by the Cosmology Science Working Group in the run up to start of operations and the Key Science Programme of SKA1.

We develop a practical methodology to remove modes from a galaxy survey power spectrum that are associated with systematic errors. We apply this to the BOSS CMASS sample, to see if it removes the excess power previously observed beyond the best-fit $\Lambda$CDM model on very large scales. We consider several possible sources of data contamination, and check whether they affect the number of targets that can be observed and the power spectrum measurements. We describe a general framework for how such knowledge can be transformed into template fields. Mode subtraction can then be used to remove these systematic contaminants at least as well as applying corrective weighting to the observed galaxies, but benefits from giving an unbiased power. Even after applying templates for all known systematics, we find a large-scale power excess, but this is reduced compared with that observed using the weights provided by the BOSS team. This excess is at much larger scales than the BAO scale and does not affect the main results of BOSS. However, it will be important for the measurement of a scale-dependent bias due to primordial non-Gaussianity. The excess is beyond that allowed by any simple model of non-Gaussianity matching Planck data, and is not matched in other surveys. We show that this power excess can further be reduced but is still present using "phenomenological" templates, designed to consider further potentially unknown sources of systematic contamination. As all discrepant angular modes can be removed using "phenomenological" templates, the potentially remaining contaminant acts radially.

When combining data sets to perform parameter inference, the results will be unreliable if there are unknown systematics in data or models. Here we introduce a flexible methodology, BACCUS: BAyesian Conservative Constraints and Unknown Systematics, which deals in a conservative way with the problem of data combination, for any degree of tension between experiments. We introduce hyperparameters that describe a bias in each model parameter for each class of experiments. A conservative posterior for the model parameters is then obtained by marginalization both over these unknown shifts and over the width of their prior. We contrast this approach with an existing hyperparameter method in which each individual likelihood is scaled, comparing the performance of each approach and their combination in application to some idealized models. Using only these rescaling hyperparameters is not a suitable approach for the current observational situation, in which internal null tests of the errors are passed, and yet different experiments prefer models that are in poor agreement. The possible existence of large shift systematics cannot be constrained with a small number of data sets, leading to extended tails on the conservative posterior distributions. We illustrate our method with the case of the $H_0$ tension between results from the cosmic distance ladder and physical measurements that rely on the standard cosmological model.

It is broadly accepted that Supermassive Black Holes (SMBHs) are located in the centers of most massive galaxies, although there is still no convincing scenario for the origin of their massive seeds. It has been suggested that primordial black holes (PBHs) of masses $\gtrsim 10^{2} M_\odot$ may provide such seeds, which would grow to become SMBHs. We suggest an observational test to constrain this hypothesis: gas accretion around PBHs during the cosmic dark ages powers the emission of high energy photons which would modify the spin temperature as measured by 21cm Intensity Mapping (IM) observations. We model and compute their contribution to the standard sky-averaged signal and power spectrum of 21cm IM, accounting for its substructure and angular dependence for the first time. If PBHs exist, the sky-averaged 21cm IM signal in absorption would be higher, while we expect an increase in the power spectrum for $\ell~\gtrsim 10^2-10^3$. We also forecast PBH detectability and measurement errors in the abundance and Eddington ratios for different fiducial parameter configurations for various future experiments, ranging from SKA to a futuristic radio array on the dark side of the Moon. While the SKA could provide a detection, only a more ambitious experiment would provide accurate measurements.

The model in which Primordial Black Holes (PBHs) constitute a non-negligible fraction of the dark matter has (re)gained popularity after the first detections of binary black hole mergers. Most of the observational constraints to date have been derived assuming a single mass for all the PBHs, although some more recent works tried to generalize constraints to the case of extended mass functions. Here we derive a general methodology to obtain constraints for any PBH Extended Mass Distribution (EMD) and any observables in the desired mass range. Starting from those obtained for a monochromatic distribution, we convert them into constraints for EMDs by using an equivalent, effective mass $M_{\rm eq}$ that depends on the specific observable. We highlight how limits of validity of the PBH modelling affect the EMD parameter space. Finally, we present converted constraints on the total abundance of PBH from microlensing, stellar distribution in ultra-faint dwarf galaxies and CMB accretion for Lognormal and Power Law mass distributions, finding that EMD constraints are generally stronger than monochromatic ones.

The possibility that a relevant fraction of the dark matter might be comprised of Primordial Black Holes (PBHs) has been seriously reconsidered after LIGO's detection of a $\sim 30 M_{\odot}$ binary black holes merger. Despite the strong interest in the model, there is a lack of studies on possible cosmological implications and effects on cosmological parameters inference. We investigate correlations with the other standard cosmological parameters using cosmic microwave background observations, finding significant degeneracies, especially with the tilt of the primordial power spectrum and the sound horizon at radiation drag. However, these degeneracies can be greatly reduced with the inclusion of small scale polarization data. We also explore if PBHs as dark matter in simple extensions of the standard $\Lambda$CDM cosmological model induces extra degeneracies, especially between the additional parameters and the PBH's ones. Finally, we present cosmic microwave background constraints on the fraction of dark matter in PBHs, not only for monochromatic PBH mass distributions but also for popular extended mass distributions. Our results show that extended mass distribution's constraints are tighter, but also that a considerable amount of constraining power comes from the high-$\ell$ polarization data. Moreover, we constrain the shape of such mass distributions in terms of the correspondent constraints on the PBH mass fraction.

DESI (Dark Energy Spectropic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. The DESI instrument is a robotically-actuated, fiber-fed spectrograph capable of taking up to 5,000 simultaneous spectra over a wavelength range from 360 nm to 980 nm. The fibers feed ten three-arm spectrographs with resolution $R= \lambda/\Delta\lambda$ between 2000 and 5500, depending on wavelength. The DESI instrument will be used to conduct a five-year survey designed to cover 14,000 deg$^2$. This powerful instrument will be installed at prime focus on the 4-m Mayall telescope in Kitt Peak, Arizona, along with a new optical corrector, which will provide a three-degree diameter field of view. The DESI collaboration will also deliver a spectroscopic pipeline and data management system to reduce and archive all data for eventual public use.

DESI (Dark Energy Spectroscopic Instrument) is a Stage IV ground-based dark energy experiment that will study baryon acoustic oscillations (BAO) and the growth of structure through redshift-space distortions with a wide-area galaxy and quasar redshift survey. To trace the underlying dark matter distribution, spectroscopic targets will be selected in four classes from imaging data. We will measure luminous red galaxies up to $z=1.0$. To probe the Universe out to even higher redshift, DESI will target bright [O II] emission line galaxies up to $z=1.7$. Quasars will be targeted both as direct tracers of the underlying dark matter distribution and, at higher redshifts ($ 2.1 < z < 3.5$), for the Ly-$\alpha$ forest absorption features in their spectra, which will be used to trace the distribution of neutral hydrogen. When moonlight prevents efficient observations of the faint targets of the baseline survey, DESI will conduct a magnitude-limited Bright Galaxy Survey comprising approximately 10 million galaxies with a median $z\approx 0.2$. In total, more than 30 million galaxy and quasar redshifts will be obtained to measure the BAO feature and determine the matter power spectrum, including redshift space distortions.

We perform a comprehensive cosmological study of the $H_0$ tension between the direct local measurement and the model-dependent value inferred from the Cosmic Microwave Background. With the recent measurement of $H_0$ this tension has raised to more than $3\sigma$. We consider changes in the early time physics without modifying the late time cosmology. We also reconstruct the late time expansion history in a model independent way with minimal assumptions using distances measures from Baryon Acoustic Oscillations and Type Ia Supernovae, finding that at $z<0.6$ the recovered shape of the expansion history is less than 5 % different than that of a standard LCDM model. These probes also provide a model insensitive constraint on the low-redshift standard ruler, measuring directly the combination $r_s h$ where $H_0=h \times 100$ km/s/Mpc and $r_s$ is the sound horizon at radiation drag (the standard ruler), traditionally constrained by CMB observations. Thus $r_s$ and $H_0$ provide absolute scales for distance measurements (anchors) at opposite ends of the observable Universe. We calibrate the cosmic distance ladder and obtain a model-independent determination of the standard ruler for acoustic scale, $r_s$. The tension in $H_0$ reflects a mismatch between our determination of $r_s$ and its standard, CMB-inferred value. Without including high-l Planck CMB polarization data (i.e., only considering the "recommended baseline" low-l polarisation and temperature and the high l temperature data), a modification of the early-time physics to include a component of dark radiation with an effective number of species around 0.4 would reconcile the CMB-inferred constraints, and the local $H_0$ and standard ruler determinations. The inclusion of the "preliminary" high-l Planck CMB polarisation data disfavours this solution.

Assuming the existence of standard rulers, standard candles and standard clocks, requiring only the cosmological principle, a metric theory of gravity, a smooth expansion history, and using state-of-the-art observations, we determine the length of the "low-redshift standard ruler". The data we use are a compilation of recent Baryon acoustic oscillation data (relying on the standard ruler), Type 1A supernov\ae\ (as standard candles), ages of early type galaxies (as standard clocks) and local determinations of the Hubble constant (as a local anchor of the cosmic distance scale). In a standard $\Lambda$CDM cosmology the "low-redshift standard ruler" coincides with the sound horizon at radiation drag, which can also be determined --in a model dependent way-- from CMB observations. However, in general, the two quantities need not coincide. We obtain constraints on the length of the low-redshift standard ruler: $r^h_{\rm s}=101.0 \pm 2.3 h^{-1}$ Mpc, when using only Type 1A supernov\ae\ and Baryon acoustic oscillations, and $r_{\rm s}=150.0\pm 4.7 $ Mpc when using clocks to set the Hubble normalisation, while $r_{\rm s}=141.0\pm 5.5 $ Mpc when using the local Hubble constant determination (using both yields $r_{\rm s}=143.9\pm 3.1 $ Mpc). The low-redshift determination of the standard ruler has an error which is competitive with the model-dependent determination from cosmic microwave background measurements made with the {\em Planck} satellite, which assumes it is the sound horizon at the end of baryon drag.

We perform an empirical consistency test of General Relativity/dark energy by disentangling expansion history and growth of structure constraints. We replace each late-universe parameter that describes the behavior of dark energy with two meta-parameters: one describing geometrical information in cosmological probes, and the other controlling the growth of structure. If the underlying model (a standard wCDM cosmology with General Relativity) is correct, that is under the null hypothesis, the two meta-parameters coincide. If they do not, it could indicate a failure of the model or systematics in the data. We present a global analysis using state-of-the-art cosmological data sets which points in the direction that cosmic structures prefer a weaker growth than that inferred by background probes. This result could signify inconsistencies of the model, the necessity of extensions to it or the presence of systematic errors in the data. We examine all these possibilities. The fact that the result is mostly driven by a specific sub-set of galaxy clusters abundance data, points to the need of a better understanding of this probe.