The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input $N$-body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra and peak counts. We find that the convergence power spectra of the lensing codes agree to $\lesssim 2\%$ out to scales $\ell \approx 4000$. For lensing peak counts, the agreement is better than $5\%$ for peaks with signal-to-noise $\lesssim 6$. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they nonetheless display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test.

The generation of simulated convergence maps is of key importance in fully exploiting weak lensing by Large Scale Structure (LSS) from which cosmological parameters can be derived. In this paper we present an extension of the PINOCCHIO code which produces catalogues of dark matter haloes so that it is capable of simulating weak lensing by LSS. Like WL-MOKA, the method starts with a random realisation of cosmological initial conditions, creates a halo catalogue and projects it onto the past-light-cone, and paints in haloes assuming parametric models for the mass density distribution within them. Large scale modes that are not accounted for by the haloes are constructed using linear theory. We discuss the systematic errors affecting the convergence power spectra when Lagrangian Perturbation Theory at increasing order is used to displace the haloes within PINOCCHIO, and how they depend on the grid resolution. Our approximate method is shown to be very fast when compared to full ray-tracing simulations from an N-Body run and able to recover the weak lensing signal, at different redshifts, with a few percent accuracy. It also allows for quickly constructing weak lensing covariance matrices, complementing PINOCCHIO's ability of generating the cluster mass function and galaxy clustering covariances and thus paving the way for calculating cross covariances between the different probes. This work advances these approximate methods as tools for simulating and analysing surveys data for cosmological purposes.

Modified Gravity (MG) scenarios have been advocated to account for the dark energy phenomenon in the universe. These models predict departures from General Relativity on large cosmic scales that can be tested through a variety of probes such as observations of galaxy clusters among others. Here, we investigate the imprint of MG models on the internal mass distribution of cluster-like halos as probed by the dark matter halo sparsity. To this purpose we perform a comparative analysis of the properties of the halo sparsity using N-body simulation halo catalogs of a standard flat $\Lambda$CDM model and MG scenarios from the DUSTGRAIN-pathfinder simulation suite. We find that the onset of the screening mechanism leaves a distinct signature in the redshift evolution of the ensemble average halos sparsity. Measurements of the sparsity of galaxy clusters from currently available mass estimates are unable to test MG models due to the large uncertainties on the cluster masses. We show that this should be possible in the future provided large cluster samples with cluster masses determined to better than $30\%$ accuracy level.

The degree of turbulent pressure support by residual gas motions in galaxy clusters is not well known. Mass modelling of combined X-ray and Sunyaev Zel'dovich observations provides an estimate of turbulent pressure support in the outer regions of several galaxy clusters. Here, we test two different filtering techniques to disentangle bulk from turbulent motions in non-radiative high-resolution cosmological simulations of galaxy clusters using the cosmological hydro code ENZO. We find that the radial behavior of the ratio of non-thermal pressure to total gas pressure as a function of cluster-centric distance can be described by a simple polynomial function. The typical non-thermal pressure support in the centre of clusters is $\sim$5%, increasing to $\sim$15% in the outskirts, in line with the pressure excess found in recent X-ray observations. While the complex dynamics of the ICM makes it impossible to reconstruct a simple correlation between turbulent motions and hydrostatic bias, we find that a relation between them can be established using the median properties of a sample of objects. Moreover, we estimate the contribution of radial accelerations to the non-thermal pressure support and conclude that it decreases moving outwards from 40% (in the core) to 15% (in the cluster's outskirts). Adding this contribution to one provided by turbulence, we show that it might account for the entire observed hydrostatic bias in the innermost regions of the clusters, and for less than 80% of it at $r > 0.8 r_{200, m}$.

Understanding how galaxy properties are linked to the dark matter halos they reside in, and how they co-evolve is a powerful tool to constrain the processes related to galaxy formation. The stellar-to-halo mass relation (SHMR) and its evolution over the history of the Universe provides insights on galaxy formation models and allows to assign galaxy masses to halos in N-body dark matter simulations. We use a statistical approach to link the observed galaxy stellar mass functions on the COSMOS field to dark matter halo mass functions from the DUSTGRAIN simulation and from a theoretical parametrization from z=0 to z=4. We also propose an empirical model to describe the evolution of the stellar-to-halo mass relation as a function of redshift. We calculate the star-formation efficiency (SFE) of galaxies and compare results with previous works and semi-analytical models.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

We present an improved determination of the total mass distribution of three CLASH/HFF massive clusters, MACS J1206.2-0847 (z=0.44), MACS J0416.1-2403 (z=0.40), Abell S1063 (z=0.35). We specifically reconstruct the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We use VLT/MUSE integral-field spectroscopy in the cluster cores to measure the stellar velocity dispersion, $\sigma$, of 40-60 member galaxies per cluster, covering 4-5 magnitudes to $m_{F160W}\simeq 21.5$. We verify the robustness and quantify the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determine the normalization and slope of the galaxy $L\mbox{-}\sigma$ Faber-Jackson relation in each cluster and use these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, $\sigma_0$, and truncation radius, $r_{cut}$, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the $\sigma_0\mbox{-}r_{cut}$ scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of $\sigma_0$, or halo masses. We derive galaxy central velocity dispersion functions of the three clusters and found them in agreement with each other. Sub-halo mass functions determined with this method can be compared with those obtained from cosmological simulations.

Modified gravity and massive neutrino cosmologies are two of the most interesting scenarios that have been recently explored to account for possible observational deviations from the concordance $\Lambda$-cold dark matter ($\Lambda$CDM) model. In this context, we investigated the large-scale structure of the Universe by exploiting the \dustp simulations that implement, simultaneously, the effects of $f(R)$ gravity and massive neutrinos. To study the possibility of breaking the degeneracy between these two effects, we analysed the redshift-space distortions in the clustering of dark matter haloes at different redshifts. Specifically, we focused on the monopole and quadrupole of the two-point correlation function, both in real and redshift space. The deviations with respect to $\Lambda$CDM model have been quantified in terms of the linear growth rate parameter. We found that redshift-space distortions provide a powerful probe to discriminate between $\Lambda$CDM and modified gravity models, especially at high redshifts ($z \gtrsim 1$), even in the presence of massive neutrinos.

Galaxy cluster counts in bins of mass and redshift have been shown to be a competitive probe to test cosmological models. This method requires an efficient blind detection of clusters from surveys with a well-known selection function and robust mass estimates. The Euclid wide survey will cover 15000 deg$^2$ of the sky in the optical and near-infrared bands, down to magnitude 24 in the $H$-band. The resulting data will make it possible to detect a large number of galaxy clusters spanning a wide-range of masses up to redshift $\sim 2$. This paper presents the final results of the Euclid Cluster Finder Challenge (CFC). The objective of these challenges was to select the cluster detection algorithms that best meet the requirements of the Euclid mission. The final CFC included six independent detection algorithms, based on different techniques, such as photometric redshift tomography, optimal filtering, hierarchical approach, wavelet and friend-of-friends algorithms. These algorithms were blindly applied to a mock galaxy catalog with representative Euclid-like properties. The relative performance of the algorithms was assessed by matching the resulting detections to known clusters in the simulations. Several matching procedures were tested, thus making it possible to estimate the associated systematic effects on completeness to $<3$%. All the tested algorithms are very competitive in terms of performance, with three of them reaching $>80$% completeness for a mean purity of 80% down to masses of $10^{14}$ M$_{\odot}$ and up to redshift $z=2$. Based on these results, two algorithms were selected to be implemented in the Euclid pipeline, the AMICO code, based on matched filtering, and the PZWav code, based on an adaptive wavelet approach. [abridged]

We present a convolutional neural network to classify distinct cosmological scenarios based on the statistically similar weak-lensing maps they generate. Modified gravity (MG) models that include massive neutrinos can mimic the standard concordance model ($\Lambda$CDM) in terms of Gaussian weak-lensing observables. An inability to distinguish viable models that are based on different physics potentially limits a deeper understanding of the fundamental nature of cosmic acceleration. For a fixed redshift of sources, we demonstrate that a machine learning network trained on simulated convergence maps can discriminate between such models better than conventional higher-order statistics. Results improve further when multiple source redshifts are combined. To accelerate training, we implement a novel data compression strategy that incorporates our prior knowledge of the morphology of typical convergence map features. Our method fully distinguishes $\Lambda$CDM from its most similar MG model on noise-free data, and it correctly identifies among the MG models with at least 80% accuracy when using the full redshift information. Adding noise lowers the correct classification rate of all models, but the neural network still significantly outperforms the peak statistics used in a previous analysis.

Satellite galaxies in galaxy clusters represent a significant fraction of the global galaxy population. Because of the unusual dense environment of clusters, their evolution is driven by different mechanisms than the ones affecting field or central galaxies. Understanding the different interactions they are subject to, and how they are influenced by them, is therefore an important step towards explaining the global picture of galaxy evolution. In this paper, we use the publicly-available high resolution hydrodynamical simulation Illustris-1 to study satellite galaxies in the three most massive host haloes (with masses $M_{200} > 10^{14}\,h^{-1}\rm{M}_{\odot}$) at $z=0$. We measure the Stellar-to-Halo Mass Relation (hereafter SHMR) of the galaxies, and find that for satellites it is shifted towards lower halo masses compared to the SHMR of central galaxies. We provide simple fitting functions for both the central and satellite SHMR. To explain the shift between the two, we follow the satellite galaxies since their time of accretion into the clusters, and quantify the impact of dark matter stripping and star formation. We find that subhaloes start losing their dark matter as soon as they get closer than $\sim 1.5\times R_{\rm{vir}}$ to the centre of their host, and that up to 80\% of their dark matter content gets stripped during infall. On the other hand, star formation quenching appears to be delayed, and galaxies continue to form stars for a few Gyr after accretion. The combination of these two effects impacts the ratio of stellar to dark matter mass which varies drastically during infall, from 0.03 to 0.3.

We present the UNIT $N$-body cosmological simulations project, designed to provide precise predictions for nonlinear statistics of the galaxy distribution. We focus on characterizing statistics relevant to emission line and luminous red galaxies in the current and upcoming generation of galaxy surveys. We use a suite of precise particle mesh simulations (FastPM) as well as with full $N$-body calculations with a mass resolution of $\sim 1.2\times10^9\,h^{-1}$M$_{\odot}$ to investigate the recently suggested technique of Angulo & Pontzen 2016 to suppress the variance of cosmological simulations We study redshift space distortions, cosmic voids, higher order statistics from $z=2$ down to $z=0$. We find that both two- and three-point statistics are unbiased. Over the scales of interest for baryon acoustic oscillations and redshift-space distortions, we find that the variance is greatly reduced in the two-point statistics and in the cross correlation between halos and cosmic voids, but is not reduced significantly for the three-point statistics. We demonstrate that the accuracy of the two-point correlation function for a galaxy survey with effective volume of 20 ($h^{-1}$Gpc)$^3$ is improved by about a factor of 40, indicating that two pairs of simulations with a volume of 1 ($h^{-1}$Gpc)$^3$ lead to the equivalent variance of $\sim$150 such simulations. The $N$-body simulations presented here thus provide an effective survey volume of about seven times the effective survey volume of DESI or Euclid. The data from this project, including dark matter fields, halo catalogues, and their clustering statistics, are publicly available at http://www.unitsims.org.

Based on the DUSTGRAIN-pathfinder suite of simulations, we investigate observational degeneracies between nine models of modified gravity and massive neutrinos. Three types of machine learning techniques are tested for their ability to discriminate lensing convergence maps by extracting dimensional reduced representations of the data. Classical map descriptors such as the power spectrum, peak counts and Minkowski functionals are combined into a joint feature vector and compared to the descriptors and statistics that are common to the field of digital image processing. To learn new features directly from the data we use a Convolutional Neural Network (CNN). For the mapping between feature vectors and the predictions of their underlying model, we implement two different classifiers; one based on a nearest-neighbour search and one that is based on a fully connected neural network. We find that the neural network provides a much more robust classification than the nearest-neighbour approach and that the CNN provides the most discriminating representation of the data. It achieves the cleanest separation between the different models and the highest classification success rate of 59% for a single source redshift. Once we perform a tomographic CNN analysis, the total classification accuracy increases significantly to 76% with no observational degeneracies remaining. Visualising the filter responses of the CNN at different network depths provides us with the unique opportunity to learn from very complex models and to understand better why they perform so well.

The combination of Galaxy-Galaxy Lensing (GGL) and Redshift Space Distortion of galaxy clustering (RSD) is a privileged technique to test General Relativity predictions, and break degeneracies between the growth rate of structure parameter $f$ and the amplitude of the linear power-spectrum $\sigma_8$. We perform a joint GGL and RSD analysis on 250 sq. degrees using shape catalogues from CFHTLenS and CFHT-Stripe 82, and spectroscopic redshifts from the BOSS CMASS sample. We adjust a model that includes non-linear biasing, RSD and Alcock-Paczynski effects. We find $f(z=0.57) =0.95\pm0.23$, $\sigma_8(z=0.57)=0.55\pm0.07$ and $\Omega_{\rm m} = 0.31\pm0.08$, in agreement with Planck cosmological results 2018. We also estimate the probe of gravity $E_{\rm G} = 0.43\pm0.10$ in agreement with $\Lambda$CDM-GR predictions of $E_{\rm G} = 0.40$. This analysis reveals that RSD efficiently decreases the GGL uncertainty on $\Omega_{\rm m}$ by a factor of 4, and by 30\% on $\sigma_8$. We use an N-body simulation supplemented by an abundance matching prescription for CMASS to build a set of overlapping lensing and clustering mocks. Together with additional spectroscopic data, this helps us to quantify and correct several systematic errors, such as photometric redshifts. We make our mock catalogues available on the Skies and Universe database.

General relativity (GR) has been well tested up to solar system scales, but it is much less certain that standard gravity remains an accurate description on the largest, that is, cosmological, scales. Many extensions to GR have been studied that are not yet ruled out by the data, including by that of the recent direct gravitational wave detections. Degeneracies among the standard model ($\Lambda$CDM) and modified gravity (MG) models, as well as among different MG parameters, must be addressed in order to best exploit information from current and future surveys and to unveil the nature of dark energy. We propose various higher-order statistics in the weak-lensing signal as a new set of observables able to break degeneracies between massive neutrinos and MG parameters. We have tested our methodology on so-called $f(R)$ models, which constitute a class of viable models that can explain the accelerated universal expansion by a modification of the fundamental gravitational interaction. We have explored a range of these models that still fit current observations at the background and linear level, and we show using numerical simulations that certain models which include massive neutrinos are able to mimic $\Lambda$CDM in terms of the 3D power spectrum of matter density fluctuations. We find that depending on the redshift and angular scale of observation, non-Gaussian information accessed by higher-order weak-lensing statistics can be used to break the degeneracy between $f(R)$ models and $\Lambda$CDM. In particular, peak counts computed in aperture mass maps outperform third- and fourth-order moments.

We perform for the first time N-body simulations of Interacting Dark Energy assuming non-Gaussian initial conditions, with the aim of investigating possible degeneracies of these two theoretically independent phenomena in different observational probes. We focus on the large-scale matter distribution, as well as on the statistical and structural properties of collapsed halos and cosmic voids. On very large scales, we show that it is possible to choose the Interaction and non-Gaussian parameters such that their effects on the halo power spectrum cancel, and the power spectrum is indistinguishable from a $\mathrm{\Lambda CDM}$ model. On small scales, measurements of the non-linear matter power spectrum, halo-matter bias, halo and subhalo mass function and cosmic void number function validate the degeneracy determined on large scales. However, the internal structural properties of halos and cosmic voids, namely halo concentration-mass relation and void density profile, are very different from those measured in the $\mathrm{\Lambda CDM}$ model, thereby breaking the degeneracy. In practice, the values of $\mathrm{f_{NL}}$ required to cancel the effect of interaction are already ruled by observations. Our results show in principle that the combination of large- and small-scale probes is needed to constrain Interacting Dark Energy and Primordial non-Gaussianity separately.

We present a novel suite of cosmological N-body simulations called the DUSTGRAIN-pathfinder, implementing simultaneously the effects of an extension to General Relativity in the form of $f(R)$ gravity and of a non-negligible fraction of massive neutrinos. We describe the generation of simulated weak lensing and cluster counts observables within a past light-cone extracted from these simulations. The simulations have been performed by means of a combination of the MG-GADGET code and a particle-based implementation of massive neutrinos, while the light-cones have been generated using the MapSim pipeline allowing us to compute weak lensing maps through a ray-tracing algorithm for different values of the source plane redshift. The mock observables extracted from our simulations will be employed for a series of papers focussed on understanding and possibly breaking the well-known observational degeneracy between $f(R)$ gravity and massive neutrinos, i.e. the fact that some specific combinations of the characteristic parameters for these two phenomena (the $f_{R0}$ scalar amplitude and the total neutrino mass $\Sigma m_{\nu}$) may result indistinguishable from the standard $\mathrm{\Lambda CDM}$ cosmology through several standard observational probes. In particular, in the present work we show how a tomographic approach to weak lensing statistics could allow - especially for the next generation of wide-field surveys - to disentangle some of the models that appear statistically indistinguishable through standard single-redshift weak lensing probe.

Galaxy clusters form at the highest density nodes of the cosmic web. The clustering of massive halos is enhanced relative to the general mass distribution and matter beyond the virial region is strongly correlated to the halo mass (halo bias). Clustering can be further enhanced depending on halo properties other than mass (secondary bias). The questions of how much and why the regions surrounding rich clusters are over-dense are still unanswered. Here, we report the analysis of the environment bias in a sample of very massive clusters, selected through the Sunyaev-Zel'dovich effect by the Planck mission. We present the first detection of the correlated dark matter associated to a single cluster, PSZ2 G099.86+58.45. The system is extremely rare in the current paradigm of structure formation. The gravitational lensing signal was traced up to 30 megaparsecs with high signal-to-noise ratio ~3.4. The measured shear is very large and points at environment matter density in notable excess of the cosmological mean. The boosting of the correlated dark matter density around high mass halos can be very effective. Together with ensemble studies of the large scale structure, lensing surveys can picture the surroundings of single haloes.

Observational cosmology is passing through a unique moment of grandeur with the amount of quality data growing fast. However, in order to better take advantage of this moment, data analysis tools have to keep up the pace. Understanding the effect of baryonic matter on the large-scale structure is one of the challenges to be faced in cosmology. In this work, we have thoroughly studied the effect of baryonic physics on different lensing statistics. Making use of the Magneticum Pathfinder suite of simulations we show that the influence of luminous matter on the 1-point lensing statistics of point sources is significant, enhancing the probability of magnified objects with $\mu>3$ by a factor of $2$ and the occurrence of multiple-images by a factor $5-500$ depending on the source redshift and size. We also discuss the dependence of the lensing statistics on the angular resolution of sources. Our results and methodology were carefully tested in order to guarantee that our uncertainties are much smaller than the effects here presented.

In this paper, we study the statistical properties of weak lensing peaks in light-cones generated from cosmological simulations. In order to assess the prospects of such observable as a cosmological probe, we consider simulations that include interacting Dark Energy (hereafter DE) models with coupling term between DE and Dark Matter. Cosmological models that produce a larger population of massive clusters have more numerous high signal-to-noise peaks; among models with comparable numbers of clusters those with more concentrated haloes produce more peaks. The most extreme model under investigation shows a difference in peak counts of about $20\%$ with respect to the reference $\mathrm{\Lambda}$CDM model. We find that peak statistics can be used to distinguish a coupling DE model from a reference one with the same power spectrum normalisation. The differences in the expansion history and the growth rate of structure formation are reflected in their halo counts, non-linear scale features and, through them, in the properties of the lensing peaks. For a source redshift distribution consistent with the expectations of future space-based wide field surveys, we find that typically seventy percent of the cluster population contributes to weak-lensing peaks with signal-to-noise ratios larger than two, and that the fraction of clusters in peaks approaches one-hundred percent for haloes with redshift z$\leq$0.5. Our analysis demonstrates that peak statistics are an important tool for disentangling DE models by accurately tracing the structure formation processes as a function of the cosmic time.

We study the statistical properties of mergers between central and satellite galaxies in galaxy clusters in the redshift range $0<z<1$, using a sample of dark-matter only cosmological N-body simulations from Le SBARBINE dataset. Using a spherical overdensity algorithm to identify dark-matter haloes, we construct halo merger trees for different values of the over-density $\Delta_c$. While the virial overdensity definition allows us to probe the accretion of satellites at the cluster virial radius $r_{vir}$, higher overdensities probe satellite mergers in the central region of the cluster, down to $\approx 0.06 r_{vir}$, which can be considered a proxy for the accretion of satellite galaxies onto central galaxies. We find that the characteristic merger mass ratio increases for increasing values of $\Delta_c$: more than $60\%$ of the mass accreted by central galaxies since $z\approx 1$ comes from major mergers. The orbits of satellites accreting onto central galaxies tend to be more tangential and more bound than orbits of haloes accreting at the virial radius. The obtained distributions of merger mass ratios and orbital parameters are useful to model the evolution of the high-mass end of the galaxy scaling relations without resorting to hydrodynamic cosmological simulations.

We investigate how Einstein rings and magnified arcs are affected by small-mass dark-matter haloes placed along the line-of-sight to gravitational lens systems. By comparing the gravitational signature of line-of-sight haloes with that of substructures within the lensing galaxy, we derive a mass-redshift relation that allows us to rescale the detection threshold (i.e. lowest detectable mass) for substructures to a detection threshold for line-of-sight haloes at any redshift. We then quantify the line-of-sight contribution to the total number density of low-mass objects that can be detected through strong gravitational lensing. Finally, we assess the degeneracy between substructures and line-of-sight haloes of different mass and redshift to provide a statistical interpretation of current and future detections, with the aim of distinguishing between CDM and WDM. We find that line-of-sight haloes statistically dominate with respect to substructures, by an amount that strongly depends on the source and lens redshifts, and on the chosen dark matter model. Substructures represent about 30 percent of the total number of perturbers for low lens and source redshifts (as for the SLACS lenses), but less than 10 per cent for high redshift systems. We also find that for data with high enough signal-to-noise ratio and angular resolution, the non-linear effects arising from a double-lens-plane configuration are such that one is able to observationally recover the line-of-sight halo redshift with an absolute error precision of 0.15 at the 68 per cent confidence level.

In the hierarchical formation model, galaxy clusters grow by accretion of smaller groups or isolated galaxies. During the infall into the centre of a cluster, the properties of accreted galaxies change. In particular, both observations and numerical simulations suggest that its dark matter halo is stripped by the tidal forces of the host. We use galaxy-galaxy weak lensing to measure the average mass of dark matter haloes of satellite galaxies as a function of projected distance to the centre of the host, for different stellar mass bins. Assuming that the stellar component of the galaxy is less disrupted by tidal stripping, stellar mass can be used as a proxy of the infall mass. We study the stellar to halo mass relation of satellites as a function of the cluster-centric distance to measure tidal stripping. We use the shear catalogues of the DES science verification archive, the CFHTLenS and the CFHT Stripe 82 (CS82) surveys, and we select satellites from the redMaPPer catalogue of clusters. For galaxies located in the outskirts of clusters, we find a stellar to halo mass relation in good agreement with the theoretical expectations from \citet{moster2013} for central galaxies. In the centre of the cluster, we find that this relation is shifted to smaller halo mass for a given stellar mass. We interpret this finding as further evidence for tidal stripping of dark matter haloes in high density environments.

We carry out a joint analysis of redshift-space distortions and galaxy-galaxy lensing, with the aim of measuring the growth rate of structure; this is a key quantity for understanding the nature of gravity on cosmological scales and late-time cosmic acceleration. We make use of the final VIPERS redshift survey dataset, which maps a portion of the Universe at a redshift of $z \simeq 0.8$, and the lensing data from the CFHTLenS survey over the same area of the sky. We build a consistent theoretical model that combines non-linear galaxy biasing and redshift-space distortion models, and confront it with observations. The two probes are combined in a Bayesian maximum likelihood analysis to determine the growth rate of structure at two redshifts $z=0.6$ and $z=0.86$. We obtain measurements of $f\sigma_8(0.6) = 0.48 \pm 0.12$ and $f\sigma_8(0.86) = 0.48 \pm 0.10$. The additional galaxy-galaxylensing constraint alleviates galaxy bias and $\sigma_8$ degeneracies, providing direct measurements of $[f(0.6),\sigma_8(0.6)] = [0.93 \pm 0.22, 0.52 \pm 0.06]$ and $f(0.86),\sigma_8(0.86)] = [0.99 \pm 0.19, 0.48 \pm 0.04]$. These measurements are statistically consistent with a Universe where the gravitational interactions can be described by General Relativity, although they are not yet accurate enough to rule out some commonly considered alternatives. Finally, as a complementary test we measure the gravitational slip parameter, $E_G$ , for the first time at $z>0.6$. We find values of $\smash{\overline{E}_G}(0.6) = 0.16 \pm 0.09$ and $\smash{\overline{E}_G}(0.86) = 0.09 \pm 0.07$, when $E_G$ is averaged over scales above $3 h^{-1} \rm{Mpc}$. We find that our $E_G$ measurements exhibit slightly lower values than expected for standard relativistic gravity in a {\Lambda}CDM background, although the results are consistent within $1-2\sigma$.

Strong gravitational lensing by galaxy clusters is a fundamental tool to study dark matter and constrain the geometry of the Universe. Recently, the Hubble Space Telescope Frontier Fields programme has allowed a significant improvement of mass and magnification measurements but lensing models still have a residual root mean square between 0.2 arcsec and few arcsec- onds, not yet completely understood. Systematic errors have to be better understood and treated in order to use strong lensing clusters as reliable cosmological probes. We have analysed two simulated Hubble-Frontier-Fields-like clusters from the Hubble Frontier Fields Comparison Challenge, Ares and Hera. We use several estimators (relative bias on magnification, den- sity profiles, ellipticity and orientation) to quantify the goodness of our reconstructions by comparing our multiple models, optimized with the parametric software LENSTOOL , with the input models. We have quantified the impact of systematic errors arising, first, from the choice of different density profiles and configurations and, secondly, from the availability of con- straints (spectroscopic or photometric redshifts, redshift ranges of the background sources) in the parametric modelling of strong lensing galaxy clusters and therefore on the retrieval of cosmological parameters. We find that substructures in the outskirts have a significant im- pact on the position of the multiple images, yielding tighter cosmological contours. The need for wide-field imaging around massive clusters is thus reinforced. We show that competitive cosmological constraints can be obtained also with complex multimodal clusters and that photometric redshifts improve the constraints on cosmological parameters when considering a narrow range of (spectroscopic) redshifts for the sources.

Full ray-tracing maps of gravitational lensing, constructed from N-Body simulations, represent a fundamental tool to interpret present and future weak lensing data. However the limitation of computational resources and storage capabilities severely restrict the number of realizations that can be performed in order to accurately sample both the cosmic shear models and covariance matrices. In this paper we present a halo model formalism for weak gravitational lensing that alleviates these issues by producing weak-lensing mocks at a reduced computational cost. Our model takes as input the halo population within a desired light-cone and the linear power spectrum of the underlined cosmological model. We examine the contribution given by the presence of substructures within haloes to the cosmic shear power spectrum and quantify it to the percent level. Our method allows us to reconstruct high-resolution convergence maps, for any desired source redshifts, of light-cones that realistically trace the matter density distribution in the universe, account for masked area and sample selections. We compare our analysis on the same large scale structures constructed using ray-tracing techniques and find very good agreements both in the linear and non-linear regimes up to few percent levels. The accuracy and speed of our method demonstrate the potential of our halo model for weak lensing statistics and the possibility to generate a large sample of convergence maps for different cosmological models as needed for the analysis of large galaxy redshift surveys.

In this paper we study the triaxial properties of dark matter haloes of a wide range of masses extracted from a set of cosmological N-body simulations. We measure the shape at different distances from the halo centre (characterised by different overdensity thresholds), both in three and in two dimensions. We discuss how halo triaxiality increases with mass, redshift and distance from the halo centre. We also examine how the orientation of the different ellipsoids are aligned with each other and what is the gradient in internal shapes for halos with different virial configurations. Our findings highlight that the internal part of the halo retains memory of the violent formation process keeping the major axis oriented toward the preferential direction of the in-falling material while the outer part becomes rounder due to continuous isotropic merging events. This effect is clearly evident in high mass haloes - which formed more recently - while it is more blurred in low mass haloes. We present simple distributions that may be used as priors for various mass reconstruction algorithms, operating in different wavelengths, in order to recover a more complex and realistic dark matter distribution of isolated and relaxed systems.

In this paper we investigate the strong lensing statistics in galaxy clusters. We extract dark matter haloes from the Millennium-XXL simulation, compute their Einstein radius distribution, and find a very good agreement with Monte Carlo predictions produced with the MOKA code. The distribution of the Einstein radii is well described by a log-normal distribution, with a considerable fraction of the largest systems boosted by different projection effects. We discuss the importance of substructures and triaxiality in shaping the size of the critical lines for cluster size haloes. We then model and interpret the different deviations, accounting for the presence of a Bright Central Galaxy (BCG) and two different stellar mass density profiles. We present scaling relations between weak lensing quantities and the size of the Einstein radii. Finally we discuss how sensible is the distribution of the Einstein radii on the cosmological parameters {\Omega}_M-{\sigma}_8 finding that cosmologies with higher {\Omega}_M and {\sigma}_8 possess a large sample of strong lensing clusters. The Einstein radius distribution may help distinguish Planck13 and WMAP7 cosmology at 3{\sigma}.

Radio and X-ray observations of galaxy clusters probe a direct link between cluster mergers and giant radio halos (RH), suggesting that these sources can be used as probes of the cluster merging rate with cosmic time. In this paper we carry out an explorative study that combines the observed fractions of merging clusters (fm) and RH (fRH) with the merging rate predicted by cosmological simulations and attempt to infer constraints on merger properties of clusters that appear disturbed in X-rays and of clusters with RH. We use morphological parameters to identify merging systems and analyze the currently largest sample of clusters with radio and X-ray data (M500>6d14 Msun, and 0.2<z<0.33, from the Planck SZ cluster catalogue). We found that in this sample fm~62-67% while fRH~44-51%. The comparison of the theoretical f_m with the observed one allows to constrain the combination (xi_m,tau_m), where xi_m and tau_m are the minimum merger mass ratio and the timescale of merger-induced disturbance. Assuming tau_m~ 2-3 Gyr, as constrained by simulations, we find that the observed f_m matches the theoretical one for xi_m~0.1-0.18. This is consistent with optical and near-IR observations of clusters in the sample (xi_m~0.14-0.16). The fact that RH are found only in a fraction of merging clusters may suggest that merger events generating RH are characterized by larger mass ratio; this seems supported by optical/near-IR observations of RH clusters in the sample (xi_min~0.2-0.25). Alternatively, RH may be generated in all mergers but their lifetime is shorter than \tau_m (by ~ fRH/fm). This is an explorative study, however it suggests that follow up studies using the forthcoming radio surveys and adequate numerical simulations have the potential to derive quantitative constraints on the link between cluster merging rate and RH at different cosmic epochs and for different cluster masses.

Gravitational lensing by clusters of galaxies offers a powerful probe of their structure and mass distribution. Deriving a lens magnification map for a galaxy cluster is a classic inversion problem and many methods have been developed over the past two decades to solve it. Several research groups have developed techniques independently to map the predominantly dark matter distribution in cluster lenses. While these methods have all provided remarkably high precision mass maps, particularly with exquisite imaging data from the Hubble Space Telescope (HST), the reconstructions themselves have never been directly compared. In this paper, we report the results of comparing various independent lens modeling techniques employed by individual research groups in the community. Here we present for the first time a detailed and robust comparison of methodologies for fidelity, accuracy and precision. For this collaborative exercise, the lens modeling community was provided simulated cluster images -- of two clusters Ares and Hera -- that mimic the depth and resolution of the ongoing HST Frontier Fields. The results of the submitted reconstructions with the un-blinded true mass profile of these two clusters are presented here. Parametric, free-form and hybrid techniques have been deployed by the participating groups and we detail the strengths and trade-offs in accuracy and systematics that arise for each methodology. We note in conclusion that lensing reconstruction methods produce reliable mass distributions that enable the use of clusters as extremely valuable astrophysical laboratories and cosmological probes.

We explore the possibility of measuring the mass accretion rate (MAR) of galaxy clusters from their mass profiles beyond the virial radius $R_{200}$. We derive the accretion rate from the mass of a spherical shell whose inner radius is $2R_{200}$, whose thickness changes with redshift, and whose infall velocity is assumed to be equal to the mean infall velocity of the spherical shells of dark matter halos extracted from $N$-body simulations. This approximation is rather crude in hierarchical clustering scenarios where both smooth accretion and aggregation of smaller dark matter halos contribute to the mass accretion of clusters.Nevertheless, in the redshift range $z=[0,2]$, our prescription returns an average MAR within $20-40 \%$ of the average rate derived from the merger trees of dark matter halos extracted from $N$-body simulations. The MAR of galaxy clusters has been the topic of numerous detailed numerical and theoretical investigations, but so far it has remained inaccessible to measurements in the real universe. Since the measurement of the mass profile of clusters beyond their virial radius can be performed with the caustic technique applied to dense redshift surveys of the cluster outer regions, our result suggests that measuring the mean MAR of a sample of galaxy clusters is actually feasible. We thus provide a new potential observational test of the cosmological and structure formation models.

In this paper, we present a new method to select the faint, background galaxies used to derive the mass of galaxy clusters by weak lensing. The method is based on the simultaneous analysis of the shear signal, that should be consistent with zero for the foreground, unlensed galaxies, and of the colors of the galaxies: photometric data from the COSMic evOlution Survey are used to train the color selection. In order to validate this methodology, we test it against a set of state-of-the-art image simulations of mock galaxy clusters in different redshift [$0.23-0.45$] and mass [$0.5-1.55\times10^{15}M_\odot$] ranges, mimicking medium-deep multicolor imaging observations (e.g. SUBARU, LBT). The performance of our method in terms of contamination by unlensed sources is comparable to a selection based on photometric redshifts, which however requires a good spectral coverage and is thus much more observationally demanding. The application of our method to simulations gives an average ratio between estimated and true masses of $\sim 0.98 \pm 0.09$. As a further test, we finally apply our method to real data, and compare our results with other weak lensing mass estimates in the literature: for this purpose we choose the cluster Abell 2219 ($z=0.228$), for which multi-band (BVRi) data are publicly available.

In this paper we present a large database of weak lensing light cones constructed using different snapshots from the Big MultiDark simulation (BigMDPL). The ray-tracing through different multiple planes has been performed with the GLAMER code accounting both for single source redshifts and for sources distributed along the cosmic time. This first paper presents weak lensing forecasts and results according to the geometry of the VIPERS-W1 and VIPERS-W4 field of view. Additional fields will be available on our database and new ones can be run upon request. Our database also contains some tools for lensing analysis. In this paper we present results for convergence power spectra, one point and high order weak lensing statistics useful for forecasts and for cosmological studies. Covariance matrices have also been computed for the different realisations of the W1 and W4 fields. In addition we compute also galaxy-shear and projected density contrasts for different halo masses at two lens redshifts according to the CFHTLS source redshift distribution both using stacking and cross-correlation techniques, finding very good agreement.

[abridged] We present a strong-lensing analysis of MACSJ0717.5+3745, based on the full depth of the Hubble Frontier Field (HFF) observations, which brings the number of multiply imaged systems to 61, ten of which are spectroscopically confirmed. The total number of images comprised in these systems rises to 165. Our analysis uses a parametric mass reconstruction technique, as implemented in the Lenstool software, to constrain a mass distribution composed of four large-scale mass components + galaxy-scale perturbers. We find a superposition of cored isothermal mass components to provide a good fit to the observational constraints, resulting in a very shallow mass distribution for the smooth (large-scale) component. Given the implications of such a flat mass profile, we investigate whether a model composed of "peaky" non-cored mass components can also reproduce the observational constraints. We find that such a non-cored mass model reproduces the observational constraints equally well. Although the total mass distributions of both models are consistent, as well as the integrated two dimensional mass profiles, we find that the smooth and the galaxy-scale components are very different. We conclude that, even in the HFF era, the generic degeneracy between smooth and galaxy-scale components is not broken, in particular in such a complex galaxy cluster. Consequently, insights into the mass distribution of MACS J0717 remain limited, underlining the need for additional probes beyond strong lensing. Our findings also have implications for estimates of the lensing magnification: we show that the amplification difference between the two models is larger than the error associated with either model. This uncertainty decreases the area of the image plane where we can reliably study the high-redshift Universe by 50 to 70%.

The next generation of large sky photometric surveys will finally be able to use arc statistics as a cosmological probe. Here we present the first of a series of papers on this topic. In particular, we study how arc counts are sensitive to the variation of two cosmological parameters: the (total) matter density parameter, $\Omega_m$, and the normalisation of the primordial power spectrum, expressed in terms of $\sigma_8$. Both these parameters influence the abundances of collapsed structures and their internal structure. We compute the expected number of gravitational arcs with various length-to-width ratios in mock light cones, by varying these cosmological parameters in the ranges $0.1\leq\Omega_m\leq0.5$ and $0.6\leq\sigma_8\leq 1$. We find that the arc counts dependence on $\Omega_m$ and $\sigma_8$ is similar, but not identical, to that of the halo counts. We investigate how the precision of the constraints on the cosmological parameters based on arc counts depends on the survey area. We find that the constraining power of arc statistics degrades critically only for surveys covering an area smaller than $10\%$ of the whole sky. Finally, we consider the case in which the search for arcs is done only in frames where galaxy clusters have been previously identified. Adopting the selection function for galaxy clusters expected to be detected from photometric data in future wide surveys, we find that less than $10\%$ of the arcs will be missed, with only a small degradation of the corresponding cosmological constraints.

The abundance of galaxy clusters can constrain both the geometry and growth of structure in our Universe. However, this probe could be significantly complicated by recent claims of nonuniversality -- non-trivial dependences with respect to the cosmological model and redshift. In this work we analyse the dependance of the mass function on the way haloes are identified and establish if this can cause departures from universality. In order to explore this dependance, we use a set of different N-body cosmological simulations (Le SBARBINE simulations), with the latest cosmological parameters from the Planck collaboration; this first suite of simulations is followed by a lower resolution set, carried out with different cosmological parameters. We identify dark matter haloes using a Spherical Overdensity algorithm with varying overdensity thresholds (virial, 2000rho_c, 1000rho_c, 500rho_c, 200rho_c and 200rho_b) at all redshifts. We notice that, when expressed in term of the rescaled variable nu, the mass functionfor virial haloes is a nearly universal as a function of redshift and cosmology, while this is clearly not the case for the other overdensities we considered. We provide fitting functions for the halo mass function parameters as a function of overdensity, that allow to predict, to within a few percent accuracy, the halo mass function for a wide range of halo definitions, redshifts and cosmological models. We then show how the departures from universality associated with other halo definitions can be derived by combining the universality of the virial definition with the expected shape of the density profile of halos.

A number of alternatives to general relativity exhibit gravitational screening in the non-linear regime of structure formation. We describe a set of algorithms that can produce weak lensing maps of large scale structure in such theories and can be used to generate mock surveys for cosmological analysis. By analysing a few basic statistics we indicate how these alternatives can be distinguished from general relativity with future weak lensing surveys.

We propose a new approach to the missing baryons problem. Building on the common assumption that the missing baryons are in the form of the Warm Hot Intergalactic Medium (WHIM), we further assumed here that the galaxy luminosity density can be used as a tracer of the WHIM. The latter assumption is supported by our finding of a significant correlation between the WHIM density and the galaxy luminosity density in the hydrodynamical simulations of Cui et al. (2012). We further found that the fraction of the gas mass in the WHIM phase is substantially (by a factor of $\sim$1.6) higher within the large scale galactic filaments, i.e. $\sim$70\%, compared to the average in the full simulation volume of $\sim$0.1\,Gpc$^3$. The relation between the WHIM overdensity and the galaxy luminosity overdensity within the galactic filaments is consistent with linear: $\delta_{\rm whim}\,=\,0.7\,\pm\,0.1\,\times\,\delta_\mathrm{LD}^{0.9 \pm 0.2}$. We applied our procedure to the line of sight to the blazar H2356-309 and found evidence for the WHIM in correspondence of the Sculptor Wall (z $\sim$0.03 and $\log{N_H}$ = $19.9^{+0.1}_{-0.3}$) and Pisces-Cetus superclusters (z $\sim$0.06 and $\log{N_H}$ = $19.7^{+0.2}_{-0.3}$), in agreement with the redshifts and column densities of the X-ray absorbers identified and studied by Fang et al. (2010) and Zappacosta et al. (2010). This agreement indicates that the galaxy luminosity density and galactic filaments are reliable signposts for the WHIM and that our method is robust in estimating the WHIM density. The signal that we detected cannot originate from the halos of the nearby galaxies since they cannot account for the large WHIM column densities that our method and X-ray analysis consistently find in the Sculptor Wall and Pisces-Cetus superclusters.

We perform multi-plane ray-tracing using the GLAMER gravitational lensing code within high-resolution light-cones extracted from the CoDECS simulations: a suite of cosmological runs featuring a coupling between Dark Energy and Cold Dark Matter. We show that the presence of the coupling is evident not only in the redshift evolution of the normalisation of the convergence power spectrum, but also in differences in non-linear structure formation with respect to {\Lambda}CDM. Using a tomographic approach under the assumption of a {\Lambda}CDM cosmology, we demonstrate that weak lensing measurements would result in a {\sigma}8 value that changes with the source redshift if the true underlying cosmology is a coupled Dark Energy one. This provides a generic null test for these types of models. We also find that different models of coupled Dark Energy can show either an enhanced or a suppressed correlation between convergence maps with differing source redshifts as compared to {\Lambda}CDM. This would provide a direct way to discriminate between different possible realisations of the coupled Dark Energy scenario. Finally, we discuss the impact of the coupling on several lensing observables for different source redshifts and angular scales with realistic source redshift distributions for current ground-based and future space-based lensing surveys.

We study the characteristics of the galaxy cluster samples expected from the European Space Agency's Euclid satellite and forecast constraints on cosmological parameters describing a variety of cosmological models. The method used in this paper, based on the Fisher Matrix approach, is the same one used to provide the constraints presented in the Euclid Red Book (Laureijs et al.2011). We describe the analytical approach to compute the selection function of the photometric and spectroscopic cluster surveys. Based on the photometric selection function, we forecast the constraints on a number of cosmological parameter sets corresponding to different extensions of the standard LambdaCDM model. The dynamical evolution of dark energy will be constrained to Delta w_0=0.03 and Delta w_a=0.2 with free curvature Omega_k, resulting in a (w_0,w_a) Figure of Merit (FoM) of 291. Including the Planck CMB covariance matrix improves the constraints to Delta w_0=0.02, Delta w_a=0.07 and a FoM=802. The amplitude of primordial non-Gaussianity, parametrised by f_NL, will be constrained to \Delta f_NL ~ 6.6 for the local shape scenario, from Euclid clusters alone. Using only Euclid clusters, the growth factor parameter \gamma, which signals deviations from GR, will be constrained to Delta \gamma=0.02, and the neutrino density parameter to Delta Omega_\nu=0.0013 (or Delta \sum m_\nu=0.01). We emphasise that knowledge of the observable--mass scaling relation will be crucial to constrain cosmological parameters from a cluster catalogue. The Euclid mission will have a clear advantage in this respect, thanks to its imaging and spectroscopic capabilities that will enable internal mass calibration from weak lensing and the dynamics of cluster galaxies. This information will be further complemented by wide-area multi-wavelength external cluster surveys that will already be available when Euclid flies. [Abridged]

The relation between mass and concentration of galaxy clusters traces their formation and evolution. Massive lensing clusters were observed to be over-concentrated and following a steep scaling in tension with predictions from the concordance $\Lambda$CDM paradigm. We critically revise the relation in the CLASH, the SGAS, the LOCUSS, and the high-redshift samples of weak lensing clusters. Measurements of mass and concentration are anti-correlated, which can bias the observed relation towards steeper values. We corrected for this bias and compared the measured relation to theoretical predictions accounting for halo triaxiality, adiabatic contraction of the halo, presence of a dominant BCG and, mostly, selection effects in the observed sample. The normalisation, the slope and the scatter of the expected relation are strongly sample-dependent. For the considered samples, the predicted slope is much steeper than that of the underlying relation characterising dark-matter only clusters. We found that the correction for statistical and selection biases in observed relations mostly solve the tension with the $\Lambda$CDM model.

Giant Radio Halos (RH) are diffuse, Mpc-sized, synchrotron radio sources observed in a fraction of merging galaxy clusters. The current scenario for the origin of RHs assumes that turbulence generated during cluster mergers re-accelerates pre-existing fossil and/or secondary electrons in the intra-cluster-medium (ICM) to the energies necessary to produce the observed radio emission. Moreover, more relaxed clusters could host diffuse "off state" halos produced by secondary electrons. In this Chapter we use Monte Carlo simulations, that combine turbulent-acceleration physics and the generation of secondaries in the ICM, to calculate the occurrence of RHs in the Universe, their spectral properties and connection with properties of the hosting clusters. Predictions for SKA1 surveys are presented at low (100-300 MHz) and mid (1-2 GHz) frequencies assuming the expected sensitivities and spatial resolutions of SKA1. SKA1 will step into an unexplored territory allowing us to study the formation and evolution of RHs in a totally new range of cluster masses and redshift, allowing firm tests of the current theoretical hypothesis. In particular, the combination of SKA1-LOW and SUR will allow the discovery of ~1000 ultrasteep- spectrum halos and to detect for the very first time "off state" RHs. We expect that at least ~2500 giant RHs will be discovered by SKA1-LOW surveys up to z~0.6. Remarkably these surveys will be sensitive to RHs in a cluster mass range (down to ~10^14 solar masses) and redshifts (up to ~1) that are unexplored by current observations. SKA1 surveys will be highly competitive with present and future SZ-surveys in the detection of high-redshift massive objects.

We have developed a new method to approach the missing baryons problem. We assume that the missing baryons reside in a form of Warm Hot Intergalactic Medium, i.e. the WHIM. Our method consists of (a) detecting the coherent large scale structure in the spatial distribution of galaxies that traces the Cosmic Web and that in hydrodynamical simulations is associated to the WHIM, (b) map its luminosity into a galaxy luminosity density field, (c) use numerical simulations to relate the luminosity density to the density of the WHIM, (d) apply this relation to real data to trace the WHIM using the observed galaxy luminosities in the Sloan Digital Sky Survey and 2dF redshift surveys. In our application we find evidence for the WHIM along the line of sight to the Sculptor Wall, at redshifts consistent with the recently reported X-ray absorption line detections. Our indirect WHIM detection technique complements the standard method based on the detection of characteristic X-ray absorption lines, showing that the galaxy luminosity density is a reliable signpost for the WHIM. For this reason, our method could be applied to current galaxy surveys to optimise the observational strategies for detecting and studying the WHIM and its properties. Our estimates of the WHIM hydrogen column density in Sculptor agree with those obtained via the X-ray analysis. Due to the additional column density estimate, our method has potential for improving the constrains of the physical parameters of the WHIM as derived with X-ray absorption, and thus for improving the understanding of the missing baryons problem.

For the last 30 years many observational and theoretical evidences have shown that galaxy clusters are not spherical objects, and that their shape is much better described by a triaxial geometry. With the advent of multi-wavelength data of increasing quality, triaxial investigations of galaxy clusters is gathering a growing interest from the community, especially in the time of "precision cosmology". In this work, we aim to provide the first statistically significant predictions in the unexplored mass range above 3x10^14 Mo/h, using haloes from two redshifts (z=0 and z=1) of the Millennium XXL simulation. The size of this cosmological dark matter only simulation (4.1 Gpc) allows the formation of a statistically significant number of massive cluster scale haloes (about 500 with M>2x10^15 Mo/h and 780000 with M>10^14 Mo/h). Besides, we aim to extend this investigation to lower masses in order to look for universal predictions across nearly six orders of magnitude in mass, from 10^10 to almost 10^16 Mo/h. For this purpose we use the SBARBINE simulations, allowing to model haloes of masses starting from 10^10 Mo/h. We use an elliptical overdensity method to select haloes and compute the shapes of the unimodal ones (approximately 50%), while we discard the unrelaxed. The minor to major and intermediate to major axis ratio are found to be well described by simple functional forms. For a given mass we can fully characterize the shape of a halo and give predictions about the distribution of axis ratios for a given cosmology and redshift. Moreover, these results are in some disagreement with the findings of Jing & Suto (2002) which are widely used in the community even though they have to be extrapolated far beyond their original mass range. This "recipe" is made available to the community in this paper and in a dedicated web page.

We present a detailed analysis of dark matter halo shapes, studying how the distributions of ellipticity, prolateness and axial ratios evolve as a function of time and mass. With this purpose in mind, we analysed the results of three cosmological simulations, running an ellipsoidal halo finder to measure triaxial halo shapes. The simulations have different scales, mass limits and cosmological parameters, which allows us to ensure a good resolution and statistics in a wide mass range, and to investigate the dependence of halo properties on the cosmological model. We confirm the tendency of haloes to be prolate at all times, even if they become more triaxial going to higher redshifts. Regarding the dependence on mass, more massive haloes are also less spherical at all redshifts, since they are the most recent forming systems and so still retain memory of their original shape at the moment of collapse. We then propose a rescaling of the shape-mass relations, using the variable $\nu = \delta_{c}/\sigma$ to represent the mass, which absorbs the dependence on both cosmology and time, allowing to find universal relations between halo masses and shape parameters (ellipticity, prolateness and the axial ratios) which hold at any redshift. This may be very useful to determine prior distributions of halo shapes for observational studies.

We present the results of a numerical study based on the analysis of the MUSIC-2 simulations, aimed at estimating the expected concentration-mass relation for the CLASH cluster sample. We study nearly 1400 halos simulated at high spatial and mass resolution, which were projected along many lines-of-sight each. We study the shape of both their density and surface-density profiles and fit them with a variety of radial functions, including the Navarro-Frenk-White, the generalised Navarro-Frenk-White, and the Einasto density profiles. We derive concentrations and masses from these fits and investigate their distributions as a function of redshift and halo relaxation. We use the X-ray image simulator X-MAS to produce simulated Chandra observations of the halos and we use them to identify objects resembling the X-ray morphologies and masses of the clusters in the CLASH X-ray selected sample. We also derive a concentration-mass relation for strong-lensing clusters. We find that the sample of simulated halos which resemble the X-ray morphology of the CLASH clusters is composed mainly by relaxed halos, but it also contains a significant fraction of un-relaxed systems. For such a sample we measure an average 2D concentration which is ~11% higher than found for the full sample of simulated halos. After accounting for projection and selection effects, the average NFW concentrations of CLASH clusters are expected to be intermediate between those predicted in 3D for relaxed and super-relaxed halos. Matching the simulations to the individual CLASH clusters on the basis of the X-ray morphology, we expect that the NFW concentrations recovered from the lensing analysis of the CLASH clusters are in the range [3-6], with an average value of 3.87 and a standard deviation of 0.61. Simulated halos with X-ray morphologies similar to those of the CLASH clusters are affected by a modest orientation bias.

We present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

We study how well halo properties of galaxy clusters, like mass and concentration, are recovered using lensing data. In order to generate a large sample of systems at different redshifts we use the code MOKA. We measure halo mass and concentration using weak lensing data alone (WL), fitting to an NFW profile the reduced tangential shear profile, or by combining weak and strong lensing data, by adding information about the size of the Einstein radius (WL+SL). For different redshifts, we measure the mass and the concentration biases and find that these are mainly caused by the random orientation of the halo ellipsoid with respect to the line-of-sight. Since our simulations account for the presence of a bright central galaxy, we perform mass and concentration measurements using a generalized NFW profile which allows for a free inner slope. This reduces both the mass and the concentration biases. We discuss how the mass function and the concentration mass relation change when using WL and WL+SL estimates. We investigate how selection effects impact the measured concentration-mass relation showing that strong lens clusters may have a concentration 20-30% higher than the average, at fixed mass, considering also the particular case of strong lensing selected samples of relaxed clusters. Finally, we notice that selecting a sample of relaxed galaxy clusters, as is done in some cluster surveys, explain the concentration-mass relation biases.

We present a study of the substructure finder dependence of subhalo clustering in the Aquarius Simulation. We run 11 different subhalo finders on the haloes of the Aquarius Simulation and we study their differences in the density profile, mass fraction and 2-point correlation function of subhaloes in haloes. We also study the mass and vmax dependence of subhalo clustering. As the Aquarius Simulation has been run at different resolutions, we study the convergence with higher resolutions. We find that the agreement between finders is at around the 10% level inside R200 and at intermediate resolutions when a mass threshold is applied, and better than 5% when vmax is restricted instead of mass. However, some discrepancies appear in the highest resolution, underlined by an observed resolution dependence of subhalo clustering. This dependence is stronger for the smallest subhaloes, which are more clustered in the highest resolution, due to the detection of subhaloes within subhaloes (the sub-subhalo term). This effect modifies the mass dependence of clustering in the highest resolutions. We discuss implications of our results for models of subhalo clustering and their relation with galaxy clustering.

We present an extension to multiple planes of the gravitational lensing code {\small GLAMER}. The method entails projecting the mass in the observed light-cone onto a discrete number of lens planes and inverse ray-shooting from the image to the source plane. The mass on each plane can be represented as halos, simulation particles, a projected mass map extracted form a numerical simulation or any combination of these. The image finding is done in a source oriented fashion, where only regions of interest are iteratively refined on an initially coarse image plane grid. The calculations are performed in parallel on shared memory machines. The code is able to handle different types of analytic halos (NFW, NSIE, power-law, etc.), haloes extracted from numerical simulations and clusters constructed from semi-analytic models ({\small MOKA}). Likewise, there are several different options for modeling the source(s) which can be distributed throughout the light-cone. The distribution of matter in the light-cone can be either taken from a pre-existing N-body numerical simulations, from halo catalogs, or are generated from an analytic mass function. We present several tests of the code and demonstrate some of its applications such as generating mock images of galaxy and galaxy cluster lenses.

We study the halo mass accretion history (MAH) and its correlation with the internal structural properties in coupled dark energy cosmologies (cDE). To accurately predict all the non-linear effects caused by dark interactions, we use the COupled Dark Energy Cosmological Simulations (CoDECS). We measure the halo concentration at z=0 and the number of substructures above a mass resolution threshold for each halo. Tracing the halo merging history trees back in time, following the mass of the main halo, we develope a MAH model that accurately reproduces the halo growth in term of M_{200} in the {\Lambda}CDM Universe; we then compare the MAH in different cosmological scenarios. For cDE models with a weak constant coupling, our MAH model can reproduce the simulation results, within 10% of accuracy, by suitably rescaling the normalization of the linear matter power spectrum at z=0, {\sigma}_8. However, this is not the case for more complex scenarios, like the "bouncing" cDE model, for which the numerical analysis shows a rapid growth of haloes at high redshifts, that cannot be reproduced by simply rescaling the value of {\sigma}_8. Moreover, at fixed value of {\sigma}_8, cold dark matter (CDM) haloes in these cDE scenarios tend to be more concentrated and have a larger amount of substructures with respect to {\Lambda}CDM predictions. Finally, we present an accurate model that relates the halo concentration to the time at which it assembles half or 4% of its mass. Combining this with our MAH model, we show how halo concentrations change while varying only {\sigma}_8 in a {\Lambda}CDM Universe, at fixed halo mass.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

In this paper we estimate the number of gravitational arcs detectable in a wide-field survey such as that which will be operated by the Euclid space mission, assuming a {\Lambda}CDM cosmology. We use the publicly available code MOKA to obtain realistic deflection angle maps of mock gravitational lenses. The maps are processed by a ray-tracing code to estimate the strong lensing cross sections of each lens. Our procedure involves 1) the generation of a light-cone which is populated with lenses drawn from a theoretical mass-function; 2) the modeling of each single lens using a triaxial halo with a NFW (Navarro-Frenk-White) density profile and theoretical concentration-mass relation, including substructures, 3) the determination of the lensing cross section as a function of redshift for each lens in the light-cone, 4) the simulation of mock observations to characterize the redshift distribution of sources that will be detectable in the Euclid images. We focus on the so-called giant arcs, i.e. gravitational arcs characterized by large length-to-width ratios (l/w > 5, 7.5 and 10). We quantify the arc detectability at different significances above the level of the background. Performing 128 different realizations of a 15,000 sq. degree survey, we find that the number of arcs detectable at 1{\sigma} above the local background will be 8912,2914, and 1275 for l/w>5, 7.5 and 10, respectively. The expected arc numbers decrease to 2409, 790, and 346 for a detection limit at 3{\sigma} above the background level. From our analysis, we find that most of the lenses which contribute to the lensing optical depth are located at redshifts 0.4<zl<0.7 and that the 50% of the arcs are images of sources at zs > 3. This is the first step towards the full characterization of the population of strong lenses that will be observed by Euclid. [abridged]

The aim of this work is to perform a systematic study of the measures of the mass and concentration estimated by fitting the convergence profile of a large sample of mock galaxy cluster size lenses, created with the publicly available code MOKA. We found that the main contribution to the bias in mass and in concentration is due to the halo triaxiality and second to the presence of substructures within the host halo virial radius. We show that knowing the cluster elongation along the line of sight helps in correcting the mass bias, but still keeps a small negative bias for the concentration. If these mass and concentration biases will characterize the galaxy cluster sample of a wide field survey it will be difficult to well recover within one sigma the cosmological parameters that mainly influence the c - M relation, using as reference a 3D c - M relation measured in cosmological N-body simulation. In this work we propose how to correct the c - M relation for projection effects and for adiabatic contraction and suggest to use these as reference for real observed data. Correcting mass and concentration estimates, as we propose, gives a measurement of the cosmological parameter within 1 - {\sigma} confidence contours.

We develop a simple model for estimating the mass growth histories of dark matter halos. The model is based on a fit to the formation time distribution, where formation is defined as the earliest time that the main branch of the merger tree contains a fraction f of the final mass M . Our analysis exploits the fact that the median formation time as a function of f is the same as the median of the main progenitor mass distribution as a function of time. When coupled with previous work showing that the concentration c of the final halo is related to the formation time tf associated with f \sim 0.04, our approach provides a simple algorithm for estimating how the distribution of halo concentrations may be expected to depend on mass, redshift and the expansion history of the background cosmology. We also show that one can predict log10 c with a precision of about 0.13 and 0.09 dex if only its mass, or both mass and tf are known. And, conversely, one can predict log10 tf from mass or c with a precision of 0.12 and 0.09 dex, approximately independent of f. Adding the mass to the c-based estimate does not result in further improvement. These latter results may be useful for studies which seek to compare the age of the stars in the central galaxy in a halo with the time the core was first assembled.

Strong gravitational lensing is a powerful tool for probing the matter distribution in the cores of massive dark matter haloes. Recent and ongoing analyses of galaxy cluster surveys (MACS, CFHTLS, SDSS, SGAS, CLASH, LoCuSS) have adressed the question of the nature of the dark matter distribution in clusters. N-body simulations of cold dark-matter haloes consistently find that haloes should be characterized by a concentration-mass relation that decreases monotonically with halo mass, and populated by a large amount of substructures, representing the cores of accreted progenitor halos. It is important for our understanding of dark matter to test these predictions. We present MOKA, a new algorithm for simulating the gravitational lensing signal from cluster-sized haloes. It implements the most recent results from numerical simulations to create realistic cluster-scale lenses with properties independent of numerical resolution. We perform systematic studies of the strong lensing cross section as a function of halo structures. We find that the strong lensing cross sections depend most strongly on the concentration and on the inner slope of the density profile of a halo, followed in order of importance by halo triaxiality and the presence of a bright central galaxy.

[Abridged] We present a comparison between weak-lensing (WL) and X-ray mass estimates of a sample of numerically simulated clusters. The sample consists on the 20 most massive objects at redshift z=0.25 and Mvir > 5 x 10^{14} Msun h^{-1}. They were found in a cosmological simulation of volume 1 h^{-3} Gpc^3, evolved in the framework of a WMAP-7 normalized cosmology. Each cluster has been resimulated at higher resolution and with more complex gas physics. We processed it thought Skylens and X-MAS to generate optical and X-ray mock observations along three orthogonal projections. The optical simulations include lensing effects on background sources. Standard observational tools and methods of analysis are used to recover the mass profiles of each cluster projection from the mock catalogues. Given the size of our sample, we could also investigate the dependence of the results on cluster morphology, environment, temperature inhomogeneity, and mass. We confirm previous results showing that WL masses obtained from the fit of the cluster tangential shear profiles with NFW functionals are biased low by ~ 5-10% with a large scatter (~10-25%). We show that scatter could be reduced by optimally selecting clusters either having regular morphology or living in substructure-poor environment. The X-ray masses are biased low by a large amount (~25-35%), evidencing the presence of both non-thermal sources of pressure in the ICM and temperature inhomogeneity, but they show a significantly lower scatter than weak-lensing-derived masses. The X-ray mass bias grows from the inner to the outer regions of the clusters. We find that both biases are weakly correlated with the third-order power ratio, while a stronger correlation exists with the centroid shift. Finally, the X-ray bias is strongly connected with temperature inhomogeneities.

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

(Abridged) We apply a very general statistical theorem introduced by Cramer (1936) to study the origin of the deviations of the halo spin PDF from the reference lognormal shape. We find that these deviations originate from correlations between two quantities entering the definition of spin, namely the ratio $J/M^{5/2}$ (which depends only on mass) and the total gravitational binding energy $E$. To reach this conclusion, we have made usage of the results deduced from two high spatial- and mass resolution simulations. Our simulations cover a relatively small volume and produce a sample of more than 16.000 gravitationally bound halos, each traced by at least 300 particles. We verify that our results are stable to different systematics, by comparing our results with those derived by the GIF2 and by a more recent simulation performed by Maccio' et al. We find that the spin probability distribution function shows systematic deviations from a lognormal, at all redshifts z <= 1. These deviations depend on mass and redshift: at small masses they change little with redshift, and also the best lognormal fits are more stable. The J-M relationship is well described by a power law of exponent $\alpha$ very near to the linear theory prediction (alpha=5/3), but systematically lower than this at z<= 0.3. We argue that the fact that deviations from a lognormal PDF are present only for high-spin halos could point to a role of large-scale tidal fields in the evolution of the spin PDF.

Accurate knowledge of the non-linear dark-matter power spectrum is important for understanding the large-scale structure of the Universe, the statistics of dark-matter haloes and their evolution, and cosmological gravitational lensing. We analytically model the dark-matter power spectrum and its cross-power spectrum with dark-matter haloes. Our model extends the halo-model formalism, including realistic substructure population within individual dark-matter haloes and the scatter of the concentration parameter at fixed halo mass. We consider three prescriptions for the mass-concentration relation and two for the substructure distribution in dark-matter haloes. We show that this extension of the halo model mainly increases the predicted power on the small scales, and is crucial for proper modeling the cosmological weak-lensing signal due to low-mass haloes. Our extended formalism shows how the halo model approach can be improved in accuracy as one increases the number of ingredients that are calibrated from n-body simulations.

We present a new algorithm for identifying the substructure within simulated dark matter haloes. The method is an extension of that proposed by Tormen et al. (2004) and Giocoli et al. (2008a), which identifies a subhalo as a group of self-bound particles that prior to being accreted by the main progenitor of the host halo belonged to one and the same progenitor halo (hereafter satellite). However, this definition does not account for the fact that these satellite haloes themselves may also have substructure, which thus gives rise to sub-subhaloes, etc. Our new algorithm identifies substructures at all levels of this hierarchy, and we use it to determine the mass function of all substructure (counting sub-haloes, sub-subhaloes, etc.). On average, haloes which formed more recently tend to have a larger mass fraction in substructure and to be less concentrated than average haloes of the same mass. We provide quantitative fits to these correlations. Even though our algorithm is very different from that of Gao et al. (2004), we too find that the subhalo mass function per unit mass at redshift z = 0 is universal. This universality extends to any redshift only if one accounts for the fact that host haloes of a given mass are less concentrated at higher redshifts, and concentration and substructure abundance are anti-correlated. This universality allows a simple parametrization of the subhalo mass function integrated over all host halo masses, at any given time. We provide analytic fits to this function which should be useful in halo model analyses which equate galaxies with halo substructure when interpreting clustering in large sky surveys. Finally, we discuss systematic differences in the subhalo mass function that arise from different definitions of (host) halo mass.

In the standard model of structure formation galaxies reside in virialized dark matter haloes which extend much beyond the observational radius of the central system. The dark matter halo formation process is hierarchical, small systems collapse at high redshift and then merge together forming larger ones. In this work we study the mass assembly history of host haloes at different observation redshifts and the mass function of accreted satellites (haloes that merge directly on the main halo progenitor). We show that the satellite mass function is universal, both independent on the host halo mass and observation redshift. The satellite mass function also turn out to be universal once only satellites before or after the host halo formation redshift (time at which the main halo progenitor assembles half of its final mass) are considered. We show that the normalizations of these distributions are directly related to the main halo progenitor mass distributions before and after its formation, while their slope and the exponential high mass cut-off remain unchanged.

We present measurements of the small-scale (0.1<~ r <~ 1 Mpc/h) quasar two-point correlation function at z>2.9, for a flux-limited (i<21) sample of 15 binary quasars compiled by Hennawi et al. (2009). The amplitude of the small-scale clustering increases from z ~ 3 to z ~ 4. The small-scale clustering amplitude is comparable to or lower than power-law extrapolations (with slope gamma=2) from the large-scale correlation function of the i<20.2 quasar sample from the Sloan Digital Sky Survey. Using simple prescriptions relating quasars to dark matter halos, we model the observed small-scale clustering with halo occupation models. Reproducing the large-scale clustering amplitude requires that the active fraction of the black holes in the central galaxies of halos is near unity, but the level of small-scale clustering favors an active fraction of black holes in satellite galaxies 0.1 <~ f_s <~ 0.5 at z >~ 3.

In hierarchical models, the time derivative of the halo mass function may be thought of as the difference of two terms - a creation term, which describes the increase in the number of haloes of mass m from mergers of less massive objects, and a destruction term, which describes the decrease in the number of m-haloes as these merge with other haloes, creating more massive haloes as a result. In models where haloes form from a spherical collapse, the distribution of halo creation times can be estimated from the same formalism which is used to estimate halo abundances: the constant-barrier excursion-set approach. In the excursion set approach, moving, rather than constant-barriers, are necessary for estimating halo abundances when the collapse is triaxial. We show how the corresponding estimates are modified by providing analytic estimates of the creation rate for a wide range of halo masses. We then show that these moving-barrier based predictions are in better agreement with measurements in numerical simulations than are the corresponding predictions of the spherical collapse model. These results should be useful for studies of merger-driven star-formation rates and AGN activity. We also present a similar study of creation of haloes conditioned on belonging to an object of a certain mass today, and reach similar conclusions - the moving barrier based estimates are in substantially better agreement with the simulations. This part of the study may be useful for understanding the tendency for the oldest stars to exist in the most massive objects.

The hierarchical growth of dark matter haloes, in which galaxies are hosted, has been studied and modeled using various approaches. In this paper we use a modified version the Sheth & Lemson algorithm for a $\mathrm{\Lambda}$ cold dark matter power spectrum, and model the growth of a Milky-Way sized halo with microsolar mass resolution, corresponding to the typical Jeans mass for a dark matter Weakly Interacting Massive Particle with mass of 100 GeV. We then compute the \emph{unevolved} subhalo mass function and build-up a Milky-Way halo placing and evolving its satellites. This subhalo population is used to study the $\gamma$-ray emission from dark matter annihilation. In this case, the subhaloes which populate the host halo have been computed considering only progenitor haloes accreted by the main branch of the tree, so as to correctly treat the embedding of sub-subhaloes inside subhaloes. Each subhalo will indeed host atthe present-time sub-subhaloes accreted when it was an isolated system. In order to compute the sub-subhalo population of a Milky-Way dwarf galaxy, like Draco, and to study its $\gamma$-ray emission, we first estimate the Draco virial mass at merging redshift $z_m$ and then we run the merger tree from $z_m$ following the halo down to the dark matter Jeans mass. We then study the effect on the Fermi-LAT (GLAST) detectability for both subhaloes in the Milky-Way and in Draco, and we show how subhaloes cannot be responsible for the boost factor needed for detection.

We present an algorithm for generating merger histories of dark matter haloes. The algorithm is based on the excursion set approach with moving barriers whose shape is motivated by the ellipsoidal collapse model of halo formation. In contrast to most other merger-tree algorithms, ours takes discrete steps in mass rather than time. This allows us to quantify effects which arise from the fact that outputs from numerical simulations are usually in discrete time bins. In addition, it suggests a natural set of scaling variables for describing the abundance of halo progenitors; this scaling is not as general as that associated with a spherical collapse. We test our algorithm by comparing its predictions with measurements in numerical simulations. The progenitor mass fractions and mass functions are in good agreement, as is the predicted scaling law. We also test the formation-redshift distribution, the mass distribution at formation, and the redshift distribution of the most recent major merger; all are in reasonable agreement with N-body simulation data, over a broad range of masses and redshifts. Finally, we study the effects of sampling in discrete time snapshots. In all cases, the improvement over algorithms based on the spherical collapse assumption is significant.

Using a cosmological N-Body simulation and a sample of re-simulated cluster-like haloes, we study the mass loss rates of dark matter subhaloes, and interpret the mass function of subhaloes at redshift zero in terms of the evolution of the mass function of systems accreted by the main halo progenitor. When expressed in terms of the ratio between the mass of the subhalo at the time of accretion and the present day host mass the unevolved subhalo mass function is found to be universal. However, the subhalo mass function at redshift zero clearly depends on $M_0$, in that more massive host haloes host more subhaloes. To relate the unevolved and evolved subhalo mass functions, we measure the subhalo mass loss rate as a function of host mass and redshift. We find that the average, specific mass loss rate of dark matter subhaloes depends mainly on redshift. These results suggest a pleasingly simple picture for the evolution and mass dependence of the evolved subhalo mass function. Less massive host haloes accrete their subhaloes earlier, which are thus subjected to mass loss for a longer time. In addition, their subhaloes are typically accreted by denser hosts, which causes an additional boost of the mass loss rate. To test the self-consistency of this picture, we use a merger trees constructed using the extended Press-Schechter formalism, and evolve the subhalo populations using the average mass loss rates obtained from our simulations, finding the subhalo mass functions to be in good agreement with the simulations. [abridged]

In the standard model of cosmic structure formation, dark matter haloes form by gravitational instability. The process is hierarchical: smaller systems collapse earlier, and later merge to form larger haloes. The galaxy clusters, hosted by the largest dark matter haloes, are at the top of this hierarchy representing the largest as well as the last structures formed in the universe, while the smaller and first haloes are those Earth-sized dark subhaloes which have been both predicted by theoretical considerations and found in numerical simulations, though it does not exist any observational hints of their existence. The probability that a halo of mass $m$ at redshift $z$ will be part of a larger halo of mass $M$ at the present time can be described in the frame of the extended Press & Schecter theory making use of the progenitor (conditional) mass function. Using the progenitor mass function we calculate analytically, at redshift zero, the distribution of subhaloes in mass, formation epoch and rarity of the peak of the density field at the formation epoch. That is done for a Milky Way-size system, assuming both a spherical and an ellipsoidal collapse model. Our calculation assumes that small progenitors do not lose mass due to dynamical processes after entering the parent halo, and that they do not interact with other subhaloes. For a $\mathrm{\Lambda}$CDM power spectrum we obtain a subhalo mass function $\mathrm{d}n/\mathrm{d}m$ proportional to $m^{- \alpha}$ with a model-independent $\alpha \sim 2$. Assuming the dark matter is a weakly interacting massive particle, the inferred distributions is used to test the feasibility of an indirect detection in the $\gamma$-rays energy band of such a population of subhaloes with a GLAST-like satellite.

A dark matter halo is said to have formed when at least half its mass hass been assembled into a single progenitor. With this definition, it is possible to derive a simple but useful analytic estimate of the distribution of halo formation times. The standard estimate of this distribution depends on the shape of the conditional mass function--the distribution of progenitor masses of a halo as a function of time. If the spherical collapse model is used to estimate the progenitor mass function, then the formation times one infers systematically underestimate those seen in numerical simulations of hierarchical gravitational clustering. We provide estimates of halo formation which may be related to an ellipsoidal collapse model. These estimates provide a substantially better description of the simulations. We also provide an alternative derivation of the formation time distribution which is based on the assumption that haloes increase their mass through binary mergers only. Our results are useful for models which relate halo structure to halo formation.

We present a simple, semi-analytical model to compute the mass functions of dark matter subhaloes. The masses of subhaloes at their time of accretion are obtained from a standard merger tree. During the subsequent evolution, the subhaloes experience mass loss due to the combined effect of dynamical friction, tidal stripping, and tidal heating. Rather than integrating these effects along individual subhalo orbits, we consider the average mass loss rate, where the average is taken over all possible orbital configurations. This allows us to write the average mass loss rate as a simple function that depends only on redshift and on the instantaneous mass ratio of subhalo and parent halo. After calibrating the model by matching the subhalo mass function (SHMF) of cluster-sized dark matter haloes obtained from numerical simulations, we investigate the predicted mass and redshift dependence of the SHMF.We find that, contrary to previous claims, the subhalo mass function is not universal. Instead, both the slope and the normalization depend on the ratio of the parent halo mass, M, and the characteristic non-linear mass M*. This simply reflects a halo formation time dependence; more massive parent haloes form later, thus allowing less time for mass loss to operate. We analyze the halo-to-halo scatter, and show that the subhalo mass fraction of individual haloes depends most strongly on their accretion history in the last Gyr. Finally we provide a simple fitting function for the average SHMF of a parent halo of any mass at any redshift and for any cosmology, and briefly discuss several implications of our findings.