In the Halo Model, galaxies are hosted by dark matter halos, while the halos themselves are biased tracers of the underlying matter distribution. Measurements of galaxy correlation functions include contributions both from galaxies in different halos, and from galaxies in the same halo (the so-called 1-halo terms). We show that, for highly biased tracers, the 1-halo term of the power spectrum obeys a steep scaling relation in terms of bias. We also show that the 1-halo term of the trispectrum has a steep scaling with bias. The steepness of these scaling relations is such that the 1-halo terms can become key contributions to the $n$-point correlation functions, even at large scales. We interpret these results through analytical arguments and semi-analytical calculations in terms of the statistical properties of halos.

We investigate the non-linear evolution of Baryon Acoustic Oscillations (BAO) in the low-redshift matter power spectrum from the DEUS-FUR $\Lambda$CDM model simulation. This is the first cosmological N-body simulation encompassing the full observable cosmic volume, thus allowing cosmic variance limited predictions at BAO scales. We control the effect of numerical systematic errors using a series of large volume high-resolution simulations. The combined analysis allows us to measure the matter power spectrum between $z=0$ and $1$ to 1% over the entire BAO range, $0.03<k [\textrm{h Mpc}^{-1}]<0.3$, in bins of size $\Delta k/k\lesssim$ 1%. We define the BAO with respect to a non-linearly evolved wiggle-free spectrum and determine the characteristics of the BAO without recurring to extrapolation from global fitting functions. We quantify the effects of non-linearities on the position and amplitude of the BAO extrema, and the coupling to the broadband slope of the power spectrum. We use these estimates to test non-linear predictions from semi-analytical models. Quite remarkably from the analysis of the redshift evolution of BAO we find that the second dip and third peak remains unaltered by non-linear effects. Furthermore, we find that the square of the damping factor and the shift of the position of BAO extrema scale to good approximation as the square of the growth factor, in agreement with expectations from perturbation theory. This confirms the idea that, besides cosmic distances, an accurate measurement of BAO at different redshifts can directly probe the growth of cosmic structures.

We study the imprint of Dark Energy on the density profile of Dark Matter halos using a set of high-resolution large volume cosmological N-body simulations from the Dark Energy Universe Simulation Series (DEUSS). We first focus on the analysis of the goodness-of-fit of the Navarro-Frenk-White (NFW) profile which we find to vary with halo mass and redshift. We also find that the fraction of halos ill-fitted by NFW varies with cosmology, thus indicating that the mass assembly of halos with perturbed density profiles carries a characteristic signature of Dark Energy. To access this information independently of any parametric profile, we introduce a new observable quantity: the halo sparsity $s_\Delta$. This is defined as the mass ratio $M_{200}/M_\Delta$, i.e. the ratio of mass inside a sphere of radius $r_{200}$ to that contained within a radius $r_\Delta$, enclosing 200 and $\Delta$ times the mean matter density respectively. We find the average sparsity to be nearly independent of the total halo mass, while its value can be inferred to better than a few percent from the ratio of the integrated halo mass functions at overdensities $\Delta$ and 200 respectively. This provides a consistency relation that can validate observational measurements of the halo sparsity. Most importantly, the sparsity significantly varies with the underlying Dark Energy model, thus providing an alternative cosmological probe.

Bayesian model selection methods provide a self-consistent probabilistic framework to test the validity of competing scenarios given a set of data. We present a case study application to strong gravitational lens parametric models. Our goal is to select a homogeneous lens subsample suitable for cosmological parameter inference. To this end we apply a Bayes factor analysis to a synthetic catalog of 500 lenses with power-law potential and external shear. For simplicity we focus on double-image lenses (the largest fraction of lens in the simulated sample) and select a subsample for which astrometry and time-delays provide strong evidence for a simple power-law model description. Through a likelihood analysis we recover the input value of the Hubble constant to within 3\sigma statistical uncertainty. We apply this methodology to a sample of double image lensed quasars. In the case of B1600+434, SBS 1520+530 and SDSS J1650+4251 the Bayes' factor analysis favors a simple power-law model description with high statistical significance. Assuming a flat \LambdaCDM cosmology, the combined likelihood data analysis of such systems gives the Hubble constant H_0=76+15-5 km/s/Mpc having marginalized over the lens model parameters, the cosmic matter density and consistently propagated the observational errors on the angular position of the images. The next generation of cosmic structure surveys will provide larger lens datasets and the method described here can be particularly useful to select homogeneous lens subsamples adapted to perform unbiased cosmological parameter inference

Strong gravitational lenses are unique cosmological probes. These produce multiple images of a single source. Whether a single galaxy, a group or a cluster, extracting cosmologically relevant information requires an accurate modeling of the lens mass distribution. A variety of models are available to this purpose, nevertheless discrimination between them as primarely relied on the quality of fit without accounting for the size of the prior model parameter space. This is a problem of model selection that we address in the Bayesian statistics framework by evaluating Bayes' factors. Using simple test cases, we show that the assumption of more complicate lens models may not be justified given the level of accuracy of the available data.

Over the past decade advancements in the understanding of several astrophysical phenomena have allowed us to infer a concordance cosmological model that successfully accounts for most of the observations of our universe. This has opened up the way to studies that aim to better determine the constants of the model and confront its predictions with those of competing scenarios. Here, we use strong gravitational lenses as cosmological probes. Strong lensing, as opposed to weak lensing, produces multiple images of a single source. Extracting cosmologically relevant information requires accurate modeling of the lens mass distribution, the latter being a galaxy or a cluster. To this purpose a variety of models are available, but it is hard to distinguish between them, as the choice is mostly guided by the quality of the fit to the data without accounting for the number of additional parameters introduced. However, this is a model selection problem rather than one of parameter fitting that we address in the Bayesian framework. Using simple test cases, we show that the assumption of more complicate lens models may not be justified given the level of accuracy of the data.