We explore the possibility that Dark Matter (DM) may be explained by a non-uniform background of approximately stellar-mass clusters of Primordial Black Holes (PBHs), by simulating the evolution them from recombination to the present with over 5000 realisations using a Newtonian $ N $-body code. We compute the cluster rate of evaporation, and extract the binary and merged sub-populations along with their parent and merger tree histories, lifetimes and formation rates; the dynamical and orbital parameter profiles, the degree of mass segregation and dynamical friction, and power spectrum of close encounters. Overall, we find that PBHs can constitute a viable DM candidate, and that their clustering presents a rich phenomenology throughout the history of the Universe. We show that binary systems constitute about 9.5% of all PBHs at present, with mass ratios of $ \bar{q}_{\mathrm{B}} = 0.154 $, and total masses of $ \bar{m}_{\mathrm{T},\mathrm{B}} = 303 \mathrm{M_\odot}$. Merged PBHs are rare, about 0.0023% of all PBHs at present, with mass ratios of $ \bar{q}_{\mathrm{M}} = 0.965 $ with total and chirp masses of $ \bar{m}_{\mathrm{T},\mathrm{M}} = 1670 \mathrm{M_\odot} $ and $ \bar{m}_{c,\mathrm{M}} = 642 \mathrm{M_\odot} $ respectively. We find that cluster puffing up and evaporation leads to bubbles of these PBHs of order 1 kpc containing at present times about 36% of objects and mass, with hundred pc sized cores. We also find that these PBH sub-haloes are distributed in wider PBH haloes of order hundreds of kpc, containing about 63% of objects and mass, coinciding with the sizes of galactic halos. We find at last high rates of close encounters of massive Black Holes ($ M \sim 1000 \mathrm{M_\odot} $), with $ \Gamma^{\mathrm{S}} = (1.2 \substack {+5.9 \\ -0.9}) \times 10^{7} \mathrm{yr^{-1} Gpc^{-3}}$ and mergers with $\Gamma^{\mathrm{M}} = 1337 \pm 41 \mathrm{yr^{-1} Gpc^{-3}} $.

The Universe on comoving scales larger than 100 Mpc/h is assumed to be statistically homogeneous, with the transition scale where we have $1\%$ deviation from homogeneity being known as the homogeneity scale $R_H$. The latter was recently proposed in Ntelis et al. (2018) [arXiv:1810.09362] to be a standard ruler. In this paper we perform a comparison between the baryon acoustic oscillations and the homogeneity scales $R_H$, assuming a spatially flat universe, linear cosmological perturbation theory and the cosmological constant cold dark matter $\Lambda$CDM model. By direct theoretical calculations, we demonstrate that the answer to the question of whether the homogeneity scale can be used as a standard ruler, is clearly negative due to the non-monotonicity of the homogeneity scale $R_H$ with the matter density parameter $\Omega_{m0}$ of the $\Lambda$CDM model. Furthermore, we also consider the effect of redshift space distortions and we find that they do not break the degeneracies but instead only change the value of $R_H$ by a few percent. Finally, we also confirm our findings with the help of an N-Body simulation.

We test the viability of the Higgs-dilaton model (HDM) compared to the evolving dark energy ($w_0 w_a$CDM) model, in which the cosmological constant model $\Lambda$CDM is also nested, by using the latest cosmological data that includes the cosmic microwave background temperature, polarization and lensing data from the \textit{Planck} satellite (2015 data release), the BICEP and Keck Array experiments, the Type Ia supernovae from the JLA catalog, the baryon acoustic oscillations from CMASS, LOWZ and 6dF, the weak lensing data from the CFHTLenS survey and the matter power Spectrum measurements from the SDSS (data release 7). We find that the values of all cosmological parameters allowed by the Higgs-dilaton inflation model are well within the \textit{Planck} satellite (2015 data release) constraints. In particular, we have that $w_0 = -1.0001^{+0.0072}_{-0.0074}$, $w_a = 0.00^{+0.15}_{-0.16}$, $n_s = 0.9693^{+0.0083}_{-0.0082}$, $\alpha_s = -0.001^{+0.013}_{-0.014}$ and $r_{0.05} = 0.0025^{+0.0017}_{-0.0016}$ (95.5\%C.L.). We also place new stringent constraints on the couplings of the Higgs-dilaton model and we find that $\xi_{\chi} < 0.00328$ and $\xi_h / \sqrt{\lambda} = 59200^{+30000}_{-20000}$ (95.5\%C.L.). Furthermore, we report that the HDM is at a slightly better footing than the $w_0 w_a$CDM model, as they both have practically the same chi-square, i.e. $\Delta \chi^2 = \chi^2_{w_0 w_a\mathrm{CDM}}-\chi^2_{\mathrm{HDM}}=0.18$, with the HDM model having two fewer parameters. Finally Bayesian evidence favors equally the two models, with the HDM being preferred by the AIC and DIC information criteria.

Gravitational wave (GW) source counts have been recently shown to be able to test how gravitational radiation propagates with the distance from the source. Here, we extend this formalism to cosmological scales, i.e. the high redshift regime, and we discuss the complications of applying this methodology to high redshift sources. We also allow for models with compactified extra dimensions like in the Kaluza-Klein model. Furthermore, we also consider the case of intermediate redshifts, i.e. $0<z\lesssim1$, where we show it is possible to find an analytical approximation for the source counts $\frac{dN}{d(S/N)}$. This can be done in terms of cosmological parameters, such as the matter density $\Omega_{m,0}$ of the cosmological constant model or the cosmographic parameters for a general dark energy model. Our analysis is as general as possible, but it depends on two important factors: a source model for the black hole binary mergers and the GW source to galaxy bias. This methodology also allows us to obtain the higher order corrections of the source counts in terms of the signal-to-noise $S/N$. We then forecast the sensitivity of future observations in constraining GW physics but also the underlying cosmology by simulating sources distributed over a finite range of signal-to-noise with a number of sources ranging from 10 to 500 sources as expected from future detectors. We find that with 500 events it will be possible to provide constraints on the matter density parameter at present $\Omega_{m,0}$ on the order of a few percent and with the precision growing fast with the number of events. In the case of extra dimensions we find that depending on the degeneracies of the model, with 500 events it may be possible to provide stringent limits on the existence of the extra dimensions if the aforementioned degeneracies can be broken.