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.

Dark matter may consist, at least partially, of primordial black holes formed during the radiation-dominated era. The radiation produced by accretion onto primordial black holes leaves characteristic signatures on the properties of the medium at high redshift, before and after Hydrogen recombination. Therefore, reliable modelling of accretion onto these objects is required to obtain robust constraints on their abundance. We investigate the effect of mechanical feedback, i.e. the impact of outflows (winds and/or jets) on the medium, on primordial black hole accretion, and thereby on the associated radiation. Using analytical and numerical calculations, we study for the first time whether outflows can reduce the accretion rate of primordial black holes with masses similar to those detected by the LIGO-Virgo collaboration. Despite the complexity of the accretion rate evolution, mechanical feedback is able to significantly reduce the primordial black hole accretion rate, at least by an order of magnitude, when outflows are aligned with the motion of the compact object. If the outflow is perpendicular to the direction of motion, the effect is less important but still non-negligible. Outflows from primordial black holes, even rather weak ones, can significantly decrease the accretion rate, effectively weakening abundance constraints on these objects. Our results motivate further numerical simulations with a more realistic setup, which would yield more precise quantitative predictions.

Are the stellar-mass merging binary black holes, recently detected by their gravitational wave signal, of stellar or primordial origin? Answering this question will have profound implications for our understanding of the Universe, including the nature of dark matter, the early Universe and stellar evolution. We build on the idea that the clustering properties of merging binary black holes can provide information about binary formation mechanisms and origin. The cross-correlation of galaxy with gravitational wave catalogues carries information about whether black hole mergers trace more closely the distribution of dark matter -- indicative of primordial origin -- or that of stars harboured in luminous and massive galaxies -- indicative of a stellar origin. We forecast the detectability of such signal for several forthcoming and future gravitational wave interferometers and galaxy surveys, including, for the first time in such analyses, an accurate modelling for the different merger rates, lensing magnification and other general relativistic effects. Our results show that forthcoming experiments could allow us to test most of the parameter space of the still viable models investigated, and shed more light on the issue of binary black hole origin and evolution.

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}$.

The detection and characterization of the Stochastic Gravitational Wave Background (SGWB) is one of the main goals of Gravitational Wave (GW) experiments. The observed SGWB will be the combination of GWs from cosmological (as predicted by many models describing the physics of the early Universe) and astrophysical origins, which will arise from the superposition of GWs from unresolved sources whose signal is too faint to be detected. Therefore, it is important to have a proper modeling of the astrophysical SGWB (ASGWB) in order to disentangle the two signals; moreover, this will provide additional information on astrophysical properties of compact objects. Applying the Cosmic Rulers formalism, we compute the observed ASGWB angular power spectrum, hence using gauge invariant quantities, accounting for all effects intervening between the source and the observer. These are the so-called projection effects, which include Kaiser, Doppler and gravitational potentials effect. Our results show that these projection effects are the most important at the largest scales, and they contribute to up to tens of percent of the angular power spectrum amplitude, with the Kaiser term being the largest at all scales. While the exact impact of these results will depend on instrumental and astrophysical details, a precise theoretical modeling of the ASGWB will necessarily need to include all these projection effects.

In the model where Primordial Black Holes (PBHs) form from large primordial curvature (C) perturbations, i.e., CPBHs, constraints on PBH abundance provide in principle constraints on the primordial curvature power spectrum. This connection however depends necessarily on the details of PBH formation mechanism. In this paper we provide, for the first time, constraints on the primordial curvature power spectrum from the latest limits on PBH abundance, taking into account all the steps from gravitational collapse in real space to PBH formation. In particular, we use results from numerical relativity simulations and peak theory to study the conditions for PBH formation for a range of perturbation shapes, including non-linearities, perturbation profile and a careful treatment of smoothing and filtering scales. We then obtain updated PBH formation conditions and translate that into primordial spectrum constraints for a wide range of shapes and abundances. These updated constraints cover a range of scales not probed by other cosmological observables. Our results show that the correct and accurate modelling of non-linearities, filtering and typical perturbation profile, is crucial for deriving meaningful cosmological implications.

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 determination of the inflationary energy scale represents one of the first step towards the understanding of the early Universe physics. The (very mild) non-Gaussian signals that arise from any inflation model carry information about the energy scale of inflation and may leave an imprint in some cosmological observables, for instance on the clustering of high-redshift, rare and massive collapsed structures. In particular, the graviton exchange contribution due to interactions between scalar and tensor fluctuations leaves a specific signature in the four-point function of curvature perturbations, thus on clustering properties of collapsed structures. We compute the contribution of graviton exchange on two- and three-point function of halos, showing that at large scales $k\sim 10^{-3}\ \mathrm{Mpc}^{-1}$ its magnitude is comparable or larger to that of other primordial non-Gaussian signals discussed in the literature. This provides a potential route to probe the existence of tensor fluctuations which is alternative and highly complementary to B-mode polarisation measurements of the cosmic microwave background radiation.

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.

Cosmological observables show a dependence with the neutrino mass, which is partially degenerate with parameters of extended models of gravity. We study and explore this degeneracy in Horndeski generalized scalar-tensor theories of gravity. Using forecasted cosmic microwave background and galaxy power spectrum datasets, we find that a single parameter in the linear regime of the effective theory dominates the correlation with the total neutrino mass. For any given mass, a particular value of this parameter approximately cancels the power suppression due to the neutrino mass at a given redshift. The extent of the cancellation of this degeneracy depends on the cosmological large-scale structure data used at different redshifts. We constrain the parameters and functions of the effective gravity theory and determine the influence of gravity on the determination of the neutrino mass from present and future surveys.