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.

Dark matter models involving a very light bosonic particle, generally known as Fuzzy Dark Matter (FDM), have been recently attracting great interest in the cosmology community, as their wave-like phenomenology would simultaneously explain the longstanding mis-detection of a dark matter particle and help easing the small-scale issues related to the standard Cold Dark Matter (CDM) scenario. With the present work, we initiate a series of papers aiming at investigating the evolution of FDM structures in a cosmological framework performed with our N-body code AX-GADGET, detailing for the first time in the literature how the actual scaling relations between solitonic cores and host haloes properties are significantly affected by the dynamical state, morphology and merger history of the individual systems. In particular, in this first paper we confirm the ability of AX-GADGET to correctly reproduce the typical FDM solitonic core and we employ it to study the non-linear evolution of eight FDM haloes in their cosmological context through the zoom-in simulation approach. We find that the scaling relations identified in previous works for isolated systems are generally modified for haloes evolving in a realistic cosmological environment, and appear to be valid only as a limit for the most relaxed and spherically symmetric systems.

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.

We use a new spherical accretion recipe tested on N-body simulations to measure the observed mass accretion rate (MAR) of 129 clusters in the CIRS and HeCS. The observed clusters cover the redshift range $0.01 < z < 0.30$ and the mass range $\sim 10^{14}-10^{15}h^{-1}\text{M}_\odot$. Based on three-dimensional mass profiles of simulated clusters reaching beyond the virial radius, our recipe returns MARs that agree with MARs based on merger trees. We adopt this recipe to estimate the MAR of real clusters based on measurement of the mass profile out to $\sim 3R_{200}$. We use the caustic method to measure the mass profiles to these large radii. We demonstrate the solidity of our estimates by applying the same approach to a set of mock redshift surveys of a sample of 2000 simulated clusters with median mass $M_{200}= 10^{14} {\,h^{-1}\rm{M_{\odot}}}$, and a sample of 50 simulated clusters with median mass $M_{200}= 10^{15} {\,h^{-1}\rm{M_{\odot}}}$: the median MARs based on the caustic mass profiles of the simulated clusters are unbiased and agree within $17\%$ with the median MARs based on the real mass profile of the clusters. The MAR of the CIRS and HeCS clusters increases with the mass and the redshift of the accreting cluster, in excellent agreement with the growth of clusters in the $\Lambda$CDM model.

The advent of a new generation of large-scale galaxy surveys is pushing cosmological numerical simulations in an uncharted territory. The simultaneous requirements of high resolution and very large volume pose serious technical challenges, due to their computational and data storage demand. In this paper, we present a novel approach dubbed Dynamic Zoom Simulations -- or DZS -- developed to tackle these issues. Our method is tailored to the production of lightcone outputs from N-body numerical simulations, which allow for a more efficient storage and post-processing compared to standard comoving snapshots, and more directly mimic the format of survey data. In DZS, the resolution of the simulation is dynamically decreased outside the lightcone surface, reducing the computational work load, while simultaneously preserving the accuracy inside the lightcone and the large-scale gravitational field. We show that our approach can achieve virtually identical results to traditional simulations at half of the computational cost for our largest box. We also forecast this speedup to increase up to a factor of 5 for larger and/or higher-resolution simulations. We assess the accuracy of the numerical integration by comparing pairs of identical simulations run with and without DZS. Deviations in the lightcone halo mass function, in the sky-projected lightcone, and in the 3D matter lightcone always remain below 0.1%. In summary, our results indicate that the DZS technique may provide a highly-valuable tool to address the technical challenges that will characterise the next generation of large-scale cosmological simulations.

Galaxy surveys are an important probe for superimposed oscillations on the primordial power spectrum of curvature perturbations, which are predicted in several theoretical models of inflation and its alternatives. In order to exploit the full cosmological information in galaxy surveys it is necessary to study the matter power spectrum to fully non-linear scales. We therefore study the non-linear clustering in models with superimposed linear and logarithmic oscillations to the primordial power spectrum by running high-resolution dark-matter-only N-body simulations. We fit a Gaussian envelope for the non-linear damping of superimposed oscillations in the matter power spectrum to the results of the N-body simulations for $k \lesssim 0.6\ h/$Mpc at $0 \leq z \leq 5$ with an accuracy below the percent. We finally use this fitting formula to forecast the capabilities of future galaxy surveys, such as Euclid and Subaru, to probe primordial oscillation down to non-linear scales alone and in combination with the information contained in CMB anisotropies.

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.

Accurate predictions for the non-linear matter power spectrum are needed to confront theory with observations in current and near future weak lensing and galaxy clustering surveys. We propose a computationally cheap method to create an emulator for modified gravity models by utilizing existing emulators for $\Lambda$CDM. Using a suite of $N$-body simulations we construct a fitting function for the enhancement of both the linear and non-linear matter power spectrum in the commonly studied Hu-Sawicki $f(R)$ gravity model valid for wave-numbers $k \lesssim 5-10\, h\text{Mpc}^{-1}$ and redshifts $z \lesssim 3$. We show that the cosmology dependence of this enhancement is relatively weak so that our fit, using simulations coming from only one cosmology, can be used to get accurate predictions for other cosmological parameters. We also show that the cosmology dependence can, if needed, be included by using linear theory, approximate $N$-body simulations (such as COLA) and semi-analytical tools like the halo model. Our final fit can easily be combined with any emulator or semi-analytical models for the non-linear $\Lambda$CDM power spectrum to accurately, and quickly, produce a non-linear power spectrum for this particular modified gravity model. The method we use can be applied to fairly cheaply construct an emulator for other modified gravity models. As an application of our fitting formula we use it to compute Fisher-forecasts for how well galaxy clustering and weak lensing in a Euclid-like survey will be at constraining modifications of gravity.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

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 and test a new numerical method to determine second-order Lagrangian displacement fields in the context of modified gravity (MG) theories. We start from the extension of Lagrangian Perturbation Theory to a class of MG models that can be described by a parametrized Poisson equation, with the introduction of a scale-dependent function. We exploit fast Fourier transforms to compute the full source term of the differential equation for the second-order Lagrangian displacement field. We compare its mean to the source term computed for specific configurations for which a k-dependent solution can be found numerically. We choose the configuration that best matches the full source term, thus obtaining an approximate factorization of the second-order displacement field as the space term valid for standard gravity times a k-dependent, second-order growth factor $D_2(k,t)$. This approximation is used to compute second order displacements for particles. The method is tested against N-body simulations run with standard and $f(R)$ gravity: we rely on the results of a friends-of-friends code run on the N-body snapshots to assign particles to halos, then compute the halo power spectrum. We find very consistent results for the two gravity theories: second-order LPT (2LPT) allows to recover the halo power spectrum of N-body simulations within $\sim 10\%$ precision to $k\sim 0.2-0.4\ h\ {\rm Mpc}^{-1}$, as well as halo positions, with an error that is a fraction of the inter-particle distance. We show that, when considering the same level of non-linearity in the density field, the performance of 2LPT with MG is the same (within $1\%$) as the one obtained for the standard $\Lambda$CDM model with General Relativity. When implemented in a computer code, this formulation of 2LPT can quickly generate dark matter distributions with $f(R)$ gravity, and can easily be extended to other MG theories.

Cosmic voids are large underdense regions that, together with galaxy clusters, filaments and walls, build up the large-scale structure of the Universe. The void size function provides a powerful probe to test the cosmological framework. However, to fully exploit this statistics, the void sample has to be properly cleaned from spurious objects. Furthermore, the bias of the mass tracers used to detect these regions has to be taken into account in the size function model. In our work we test a cleaning algorithm and a new void size function model on a set of simulated dark matter halo catalogues, with different mass and redshift selections, to investigate the statistics of voids identified in a biased mass density field. We then investigate how the density field tracers' bias affects the detected size of voids. The main result of this analysis is a new model of the size function, parameterised in terms of the linear effective bias of the tracers used, which is straightforwardly inferred from the large-scale two-point correlation function. This represents a crucial step to exploit the method on real data catalogues. The proposed size function model has been accurately calibrated on mock catalogues, and used to validate the possibility to provide forecasts on the cosmological constraints, namely on the matter density contrast, $\Omega_{\rm M}$, and on the normalisation of the linear matter power spectrum, $\sigma_8$, at different redshifts.

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.

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.

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.

Cosmic voids are the major volume component in the matter distribution of the Universe. They posses great potential for constraining dark energy as well as for testing theories of gravity. Nevertheless, in spite of their growing popularity as cosmological probes, a gap of knowledge between cosmic void observations and theory still persists. In particular, the void size function models proposed in literature have been proven unsuccessful in reproducing the results obtained from cosmological simulations in which cosmic voids are detected from biased tracers of the density field, undermining the possibility of using them as cosmological probes. The goal of this work is to cover this gap. In particular, we make use of the findings of a previous work in which we have improved the void selection procedure, presenting an algorithm that redefines the void ridges and, consequently, their radius. By applying this algorithm, we validate the volume conserving model of the void size function on a set of unbiased simulated density field tracers. We highlight the difference in the internal structure between voids selected in this way and those identified by the popular VIDE void finder. We also extend the validation of the model to the case of biased tracers. We find that a relation exists between the tracer used to sample the underlying dark matter density field and its unbiased counterpart. Moreover, we demonstrate that, as long as this relation is accounted for, the size function is a viable approach for studying cosmology with cosmic voids.

Searches for modified gravity in the large-scale structure try to detect the enhanced amplitude of density fluctuations caused by the fifth force present in many of these theories. Neutrinos, on the other hand, suppress structure growth below their free-streaming length. Both effects take place on comparable scales, and uncertainty in the neutrino mass leads to a degeneracy with modified gravity parameters for probes that are measuring the amplitude of the matter power spectrum. We explore the possibility to break the degeneracy between modified gravity and neutrino effects in the growth of structures by considering kinematic information related to either the growth rate on large scales or the virial velocities inside of collapsed structures. In order to study the degeneracy up to fully non-linear scales, we employ a suite of $N$-body simulations including both $f(R)$ modified gravity and massive neutrinos. Our results indicate that velocity information provides an excellent tool to distinguish massive neutrinos from modified gravity. Models with different values of neutrino masses and modified gravity parameters possessing a comparable matter power spectrum at a given time have different growth rates. This leaves imprints in the velocity divergence, which is therefore better suited than the amplitude of density fluctuations to tell the models apart. In such models with a power spectrum comparable to $\Lambda$CDM today, the growth rate is strictly enhanced. We also find the velocity dispersion of virialised clusters to be well suited to constrain deviations from general relativity without being affected by the uncertainty in the sum of neutrino masses.

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.

Fuzzy Dark Matter (FDM) represents an alternative and intriguing description of the standard Cold Dark Matter (CDM) fluid, able to explain the lack of direct detection of dark matter particles in the GeV sector and to alleviate small scales tensions in the cosmic large-scale structure formation. Cosmological simulations of FDM models in the literature were performed either with very expensive high-resolution grid-based simulations of individual haloes or through N-body simulations encompassing larger cosmic volumes but resorting on significant approximations in the FDM non-linear dynamics to reduce their computational cost. With the use of the new N-body cosmological hydrodynamical code AX-GADGET, we are now able not only to overcome such numerical problems, but also to combine a fully consistent treatment of FDM dynamics with the presence of gas particles and baryonic physical processes, in order to quantify the FDM impact on specific astrophysical observables. In particular, in this paper we perform and analyse several hydrodynamical simulations in order to constrain the FDM mass by quantifying the impact of FDM on Lyman-$\alpha$ forest observations, as obtained for the first time in the literature in a N-body setup without approximating the FDM dynamics. We also study the statistical properties of haloes, exploiting the large available sample, to extract information on how FDM affects the abundance, the shape, and density profiles of dark matter haloes.

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.

We present a new module of the parallel N-Body code P-GADGET3 for cosmological simulations of light bosonic non-thermal dark matter, often referred as Fuzzy Dark Matter (FDM). The dynamics of the FDM features a highly non-linear Quantum Potential (QP) that suppresses the growth of structures at small scales. Most of the previous attempts of FDM simulations either evolved suppressed initial conditions, completely neglecting the dynamical effects of QP throughout cosmic evolution, or resorted to numerically challenging full-wave solvers. The code provides an interesting alternative, following the FDM evolution without impairing the overall performance. This is done by computing the QP acceleration through the Smoothed Particle Hydrodynamics (SPH) routines, with improved schemes to ensure precise and stable derivatives. As an extension of the P-GADGET3 code, it inherits all the additional physics modules implemented up to date, opening a wide range of possibilities to constrain FDM models and explore its degeneracies with other physical phenomena. Simulations are compared with analytical predictions and results of other codes, validating the QP as a crucial player in structure formation at small scales.

We develop the Blooming Tree Algorithm, a new technique that uses spectroscopic redshift data alone to identify the substructures and the surrounding groups of galaxy clusters, along with their member galaxies. Based on the estimated binding energy of galaxy pairs, the algorithm builds a binary tree that hierarchically arranges all the galaxies in the field of view. The algorithm searches for buds, corresponding to gravitational potential minima on the binary tree branches; for each bud, the algorithm combines the number of galaxies, their velocity dispersion and their average pairwise distance into a parameter that discriminates between the buds that do not correspond to any substructure or group, and thus eventually die, and the buds that correspond to substructures and groups, and thus bloom into the identified structures. We test our new algorithm with a sample of 300 mock redshift surveys of clusters in different dynamical states; the clusters are extracted from a large cosmological $N$-body simulation of a $\Lambda$CDM model. We limit our analysis to substructures and surrounding groups identified in the simulation with mass larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within 6 $h^{-1}$~Mpc from the cluster center, the technique recovers $ 80$\% of the real substructures and $ 60$\% of the surrounding groups; in $57$\% of the identified structures, at least 60\% of the member galaxies of the substructures and groups belong to the same real structure. These results improve by roughly a factor of two the performance of the best substructure identification algorithm currently available, the $\sigma$ plateau algorithm, and suggest that our Blooming Tree Algorithm can be an invaluable tool for detecting substructures of galaxy clusters and investigating their complex dynamics.

We present a halo mass function accurate over the full relevant Hu-Sawicki $f(R)$ parameter space based on spherical collapse calculations and calibrated to a suite of modified gravity $N$-body simulations that include massive neutrinos. We investigate the ability of current and forthcoming galaxy cluster observations to detect deviations from general relativity while constraining the total neutrino mass and including systematic uncertainties. Our results indicate that the degeneracy between massive neutrino and modify gravity effects is a limiting factor for the current searches for new gravitational physics with clusters of galaxies, but future surveys will be able to break the degeneracy.

The key to understand the nature of dark energy relies in our ability to probe the distant Universe. In this framework, the recent detection of the kinematic Sunyaev-Zel'dovich (kSZ) effect signature in the cosmic microwave background obtained with the South Pole Telescope (SPT) is extremely useful since this observable is sensitive to the high-redshift diffuse plasma. We analyse a set of cosmological hydrodynamical simulation with 4 different realisations of a Hu & Sawicki $f(R)$ gravity model, parametrised by the values of $\overline{f}_{\rm R,0}=(0,-10^{-6},-10^{-5},-10^{-4})$, to compute the properties of the kSZ effect due to the ionized Universe and how they depend on $\overline{f}_{\rm R,0}$ and on the redshift of reionization, $z_{\rm re}$. In the standard General Relativity limit ($\overline{f}_{\rm R,0}$=0) we obtain an amplitude of the kSZ power spectrum of $\mathcal{D}^{\rm kSZ}_{3000}=4.1\,$$\mu$K$^2$ ($z_{\rm re}$=8.8), close to the $+1\sigma$ limit of the $\mathcal{D}^{\rm kSZ}_{3000}=(2.9\pm1.3)\,$$\mu$K$^2$ measurement by SPT. This corresponds to an upper limit on the kSZ contribute from patchy reionization of $\mathcal{D}^{\rm kSZ,patchy}_{3000}<0.9\,$$\mu$K$^2$ (95 per cent confidence level). Modified gravity boosts the kSZ signal by about 3, 12 and 50 per cent for $\overline{f}_{\rm R,0}=(-10^{-6},-10^{-5},-10^{-4})$, respectively, with almost no dependence on the angular scale. This means that with modified gravity the limits on patchy reionization shrink significantly: for $\overline{f}_{\rm R,0}=-10^{-5}$ we obtain $\mathcal{D}^{\rm kSZ,patchy}_{3000}<0.4\,$$\mu$K$^2$. Finally, we provide an analytical formula for the scaling of the kSZ power spectrum with $z_{\rm re}$ and $\overline{f}_{\rm R,0}$ at different multipoles: at $\ell=3000$ we obtain $\mathcal{D}^{\rm kSZ}_{3000}\propto z_{\rm re}^{0.24}\left(1+\sqrt{\left|\overline{f}_{\rm R,0}\right|}\right)^{41}$.

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 investigate the capabilities of perturbation theory in capturing non-linear effects of dark energy. We test constant and evolving $\omega$ models, as well as models involving momentum exchange between dark energy and dark matter. Specifically, we compare perturbative predictions at 1-loop level against N-body results for four non-standard equations of state as well as varying degrees of momentum exchange between dark energy and dark matter. The interaction is modelled phenomenologically using a time dependent drag term in the Euler equation. We make comparisons at the level of the matter power spectrum and the redshift space monopole and quadrupole. The multipoles are modelled using the Taruya, Nishimichi and Saito (TNS) redshift space spectrum. We find perturbation theory does very well in capturing non-linear effects coming from dark sector interaction. We isolate and quantify the 1-loop contribution coming from the interaction and from the non-standard equation of state. We find the interaction parameter $\xi$ amplifies scale dependent signatures in the range of scales considered. Non-standard equations of state also give scale dependent signatures within this same regime. In redshift space the match with N-body is improved at smaller scales by the addition of the TNS free parameter $\sigma_v$. To quantify the importance of modelling the interaction, we create mock data sets for varying values of $\xi$ using perturbation theory. This data is given errors typical of Stage IV surveys. We then perform a likelihood analysis using the first two multipoles on these sets and a $\xi=0$ modelling, ignoring the interaction. We find the fiducial growth parameter $f$ is generally recovered even for very large values of $\xi$ both at $z=0.5$ and $z=1$. The $\xi=0$ modelling is most biased in its estimation of $f$ for the phantom $\omega=-1.1$ case.

The large-scale structure of the universe can only be observed via luminous tracers of the dark matter. However, the clustering statistics of tracers are biased and depend on various properties, such as their host-halo mass and assembly history. On very large scales this tracer bias results in a constant offset in the clustering amplitude, known as linear bias. Towards smaller nonlinear scales, this is no longer the case and tracer bias becomes a complicated function of scale and time. We focus on tracer bias centred on cosmic voids, depressions of the density field that spatially dominate the universe. We consider three types of tracers: galaxies, galaxy clusters and AGN, extracted from the hydrodynamical simulation Magneticum Pathfinder. In contrast to common clustering statistics that focus on auto-correlations of tracers, we find that void-tracer cross-correlations are successfully described by a linear-bias relation. The tracer-density profile of voids can thus be related to their matter-density profile by a single number. We show that it coincides with the linear tracer bias extracted from the large-scale auto-correlation function and expectations from theory, if sufficiently large voids are considered. For smaller voids we observe a shift towards higher values. This has important consequences on cosmological parameter inference, as the problem of unknown tracer bias is alleviated up to a constant number. The smallest scales in existing datasets become accessible to simpler models, providing numerous modes of the density field that have been disregarded so far, but may help to further reduce statistical errors in constraining cosmology.

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.

The study of neutrinos in astrophysics requires the combination of different observational probes. The temperature anisotropies of the cosmic microwave background induced via the kinematic Sunyaev-Zel'dovich (kSZ) effect may provide interesting information since they are expected to receive significant contribution from high-redshift plasma. We present a set of cosmological hydrodynamical simulations that include a treatment of the neutrino component considering four different sum of neutrino masses: $\Sigma m_\nu=(0,0.15,0.3,0.6)$ eV. Using their outputs, we modelled the kSZ effect due to the large-scale structure after the reionization by producing mock maps, then computed the kSZ power spectrum and studied how it depends on $z_{\rm re}$ and $\Sigma m_\nu$. We also run a set of four simulations to study and correct possible systematics due to resolution, finite box size and astrophysics. With massless neutrinos we obtain $\mathcal{D}^{\rm kSZ}_{3000}=4.0$ $\mu {\rm K}^2$ ($z_{\rm re}$=8.8), enough to account for all of the kSZ signal of $\mathcal{D}^{\rm kSZ}_{3000}=(2.9\pm1.3)$ $\mu {\rm K}^2$ measured with the South Pole Telescope. This translates into an upper limit on the kSZ effect due to patchy reionization of $\mathcal{D}^{\rm kSZ,patchy}_{3000}<1.0$ $\mu {\rm K}^2$ (95 per cent confidence level). Massive neutrinos induce a damping of kSZ effect power of about 8, 12 and 40 per cent for $\Sigma m_\nu=(0.15,0.3,0.6)$ eV, respectively. We study the dependence of the kSZ signal with $z_{\rm re}$ and the neutrino mass fraction, $f_\nu$, and obtain $\mathcal{D}^{\rm kSZ}_{3000}\propto z_{\rm re}^{0.26}(1-f_\nu)^{14.3}$. Interestingly, the scaling with $f_\nu$ is significantly shallower with respect to the equivalent thermal SZ effect, and may be used to break the degeneracy with other cosmological parameters.

We consider cosmological models in which dark matter feels a fifth force mediated by the dark energy scalar field, also known as coupled dark energy. Our interest resides in estimating forecasts for future surveys like Euclid when we take into account non-linear effects, relying on new fitting functions that reproduce the non-linear matter power spectrum obtained from N-body simulations. We obtain fitting functions for models in which the dark matter-dark energy coupling is constant. Their validity is demonstrated for all available simulations in the redshift range $z=0-1.6$ and wave modes below $k=10 \text{h/Mpc}$. These fitting formulas can be used to test the predictions of the model in the non-linear regime without the need for additional computing-intensive N-body simulations. We then use these fitting functions to perform forecasts on the constraining power that future galaxy-redshift surveys like Euclid will have on the coupling parameter, using the Fisher matrix method for galaxy clustering (GC) and weak lensing (WL). We find that by using information in the non-linear power spectrum, and combining the GC and WL probes, we can constrain the dark matter-dark energy coupling constant squared, $\beta^{2}$, with precision smaller than 4\% and all other cosmological parameters better than 1\%, which is a considerable improvement of more than an order of magnitude compared to corresponding linear power spectrum forecasts with the same survey specifications.

We present for the first time the outcomes of a cosmological N-body simulation that simultaneously implements a Warm Dark Matter (WDM) particle candidate and a modified gravitational interaction in the form of $f(R)$ gravity, and compare its results with the individual effects of these two independent extensions of the standard $\Lambda $CDM scenario, and with the reference cosmology itself. We consider a rather extreme value of the WDM particle mass ($m_{\rm WDM}=0.4$ keV) and a single realisation of $f(R)$ gravity with $|\bar{f}_{R0}|=10^{-5}$, and we investigate the impact of these models and of their combination on a wide range of cosmological observables with the aim to identify possible observational degeneracies. In particular, we focus on the large-scale matter distribution, as well as on the statistical and structural properties of collapsed halos and cosmic voids. Differently from the case of combining $f(R)$ gravity with massive neutrinos -- previously investigated in Baldi et al. (2014) -- we find that most of the considered observables do not show any significant degeneracy due to the fact that WDM and $f(R)$ gravity are characterised by individual observational footprints with a very different functional dependence on cosmic scales and halo masses. In particular, this is the case for the nonlinear matter power spectrum in real space, for the halo and sub-halo mass functions, for the halo density profiles and for the concentration-mass relation. However, other observables -- like e.g. the halo bias -- do show some level of degeneracy between the two models, while a very strong degeneracy is observed for the nonlinear matter power spectrum in redshift space, for the density profiles of small cosmic voids -- with radius below $\approx 5$ Mpc$/h$ -- and for the voids abundance as a function of the void core density.

Euclid is a European Space Agency medium class mission selected for launch in 2020 within the Cosmic Vision 2015 2025 program. The main goal of Euclid is to understand the origin of the accelerated expansion of the universe. Euclid will explore the expansion history of the universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

In the current state of cosmology, where cosmological parameters are being measured to percent accuracy, it is essential to understand all sources of error to high precision. In this paper we present the results of a study of the local variations in the Hubble constant measured at the distance scale of the Coma Cluster, and test the validity of correcting for the peculiar velocities predicted by gravitational instability theory. The study is based on N-body simulations, and includes models featuring a coupling between dark energy and dark matter, as well as two $\Lambda$CDM simulations with different values of $\sigma_8$. It is found that the variance in the local flows is significantly larger in the coupled models, which increases the uncertainty in the local measurements of the Hubble constant in these scenarios. By comparing the results from the different simulations, it is found that most of the effect is caused by the higher value of $\sigma_8$ in the coupled cosmologies, though this cannot account for all of the additional variance. Given the discrepancy between different estimates of the Hubble constant in the universe today, cosmological models causing a greater cosmic variance is something that we should be aware of.

Persisting tensions between the cosmological constraints derived from low-redshift probes and the ones obtained from temperature and polarisation anisotropies of the Cosmic Microwave Background -- although not yet providing compelling evidence against the $\Lambda $CDM model -- seem to consistently indicate a slower growth of density perturbations as compared to the predictions of the standard cosmological scenario. Such behavior is not easily accommodated by the simplest extensions of General Relativity, such as f(R) models, which generically predict an enhanced growth rate. In the present work we present the outcomes of a suite of large N-body simulations carried out in the context of a cosmological model featuring a non-vanishing scattering cross section between the dark matter and the dark energy fields, for two different parameterisations of the dark energy equation of state. Our results indicate that these Dark Scattering models have very mild effects on many observables related to large-scale structures formation and evolution, while providing a significant suppression of the amplitude of linear density perturbations and the abundance of massive clusters. Our simulations therefore confirm that these models offer a promising route to alleviate existing tensions between low-redshift measurements and those of the CMB.

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.

We perform for the first time high-resolution zoom-in re-simulations of individual halos in the context of the Multi-coupled Dark Energy (McDE) scenario, which is characterised by the existence of two distinct dark matter particle species with opposite couplings to a Dark Energy scalar field. We compare the structural properties of the simulated halos to the standard Lambda-CDM results. The zoomed-in initial conditions are set up using a specifically designed code called ZInCo that we publicly release along with the present paper. Our numerical results allow to investigate in detail and with unprecedented resolution the halo segregation process that characterises McDE cosmologies from its very early stages. In particular, we find that in contrast to what could be inferred from previous numerical analysis at lower resolution, the segregation process is already in place at redshifts as high as z ~ 7. Most remarkably, we find that the subsequent evolution of the segregation leads to the formation of cored total matter density profiles with a core size that progressively increases in time. The shape of the cored profiles can be accurately predicted as the superposition of two NFW profiles with an increasing offset, thereby confirming the interpretation of the simulations results in terms of the segregation of the two dark matter components of the halo as a consequence of their different coupling to the Dark Energy field.

We test the imprint of f(R) modified gravity on the halo mass function, using N-body simulations and a theoretical model developed in (Kopp et al. 2013). We find a very good agreement between theory and simulations. We extend the theoretical model to the conditional mass function and apply it to the prediction of the linear halo bias in f(R) gravity. Using the halo model we obtain a prediction for the non-linear matter power spectrum accurate to ~10% at z=0 and up to k=2h/Mpc. We also study halo profiles for the f(R) models and find a deviation from the standard general relativity result up to 40%, depending on the halo masses and redshift. This has not been pointed out in previous analysis. Finally we study the number density and profiles of voids identified in these f(R) N-body simulations. We underline the effect of the bias and the sampling to identify voids. We find significant deviation from GR when measuring the f(R) void profiles with fR0<-10^{-6}.

Euclid is a European Space Agency medium class mission selected for launch in 2019 within the Cosmic Vision 2015-2025 programme. The main goal of Euclid is to understand the origin of the accelerated expansion of the Universe. Euclid will explore the expansion history of the Universe and the evolution of cosmic structures by measuring shapes and redshifts of galaxies as well as the distribution of clusters of galaxies over a large fraction of the sky. Although the main driver for Euclid is the nature of dark energy, Euclid science covers a vast range of topics, from cosmology to galaxy evolution to planetary research. In this review we focus on cosmology and fundamental physics, with a strong emphasis on science beyond the current standard models. We discuss five broad topics: dark energy and modified gravity, dark matter, initial conditions, basic assumptions and questions of methodology in the data analysis. This review has been planned and carried out within Euclid's Theory Working Group and is meant to provide a guide to the scientific themes that will underlie the activity of the group during the preparation of the Euclid mission.

Self-consistent ${\it N}$-body simulations of modified gravity models are a key ingredient to obtain rigorous constraints on deviations from General Relativity using large-scale structure observations. This paper provides the first detailed comparison of the results of different ${\it N}$-body codes for the $f(R)$, DGP, and Symmetron models, starting from the same initial conditions. We find that the fractional deviation of the matter power spectrum from $\Lambda$CDM agrees to better than $1\%$ up to $k \sim 5-10~h/{\rm Mpc}$ between the different codes. These codes are thus able to meet the stringent accuracy requirements of upcoming observational surveys. All codes are also in good agreement in their results for the velocity divergence power spectrum, halo abundances and halo profiles. We also test the quasi-static limit, which is employed in most modified gravity ${\it N}$-body codes, for the Symmetron model for which the most significant non-static effects among the models considered are expected. We conclude that this limit is a very good approximation for all of the observables considered here.

We investigate the power of the caustic technique for identifying substructures of galaxy clusters from optical redshift data alone. The caustic technique is designed to estimate the mass profile of galaxy clusters to radii well beyond the virial radius, where dynamical equilibrium does not hold. Two by-products of this technique are the identification of the cluster members and the identification of the cluster substructures. We test the caustic technique as a substructure detector on two samples of 150 mock redshift surveys of clusters; the clusters are extracted from a large cosmological $N$-body simulation of a $\Lambda$CDM model and have masses of $M_{200} \sim 10^{14} h^{-1} M_{\odot}$ and $M_{200} \sim 10^{15} h^{-1} M_{\odot}$ in the two samples. We limit our analysis to substructures identified in the simulation with masses larger than $10^{13} h^{-1} M_{\odot}$. With mock redshift surveys with 200 galaxies within $3R_{200}$, (1) the caustic technique recovers $\sim 30-50$\% of the real substructures, and (2) $\sim 15-20$\% of the substructures identified by the caustic technique correspond to real substructures of the central cluster, the remaining fraction being low-mass substructures, groups or substructures of clusters in the surrounding region, or chance alignments of unrelated galaxies. These encouraging results show that the caustic technique is a promising approach for investigating the complex dynamics of galaxy clusters.

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.

In this work we analyse the properties of cosmic voids in standard and coupled dark energy cosmologies. Using large numerical simulations, we investigate the effects produced by the dark energy coupling on three statistics: the filling factor, the size distribution and the stacked profiles of cosmic voids. We find that the bias of the tracers of the density field used to identify the voids strongly influences the properties of the void catalogues, and, consequently, the possibility of using the identified voids as a probe to distinguish coupled dark energy models from the standard $\Lambda $CDM cosmology. In fact, on one hand coupled dark energy models are characterised by an excess of large voids in the cold dark matter distribution as compared to the reference standard cosmology, due to their higher normalisation of linear perturbations at low redshifts. Specifically, these models present an excess of large voids with $R_{eff}>20, 15, 12$ Mpc h^{-1}, at $z=0, 0.55, 1$, respectively. On the other hand, we do not find any significant difference in the properties of the void detected in the distribution of collapsed dark matter halos. These results imply that the tracer bias has a significant impact on the possibility of using cosmic void catalogues to probe cosmology.

Among the possible alternatives to the standard cosmological model ($\Lambda$CDM), coupled Dark Energy models postulate that Dark Energy (DE), seen as a dynamical scalar field, may interact with Dark Matter (DM), giving rise to a "fifth-force", felt by DM particles only. In this paper, we study the impact of these cosmologies on the statistical properties of galaxy populations by combining high-resolution numerical simulations with semi-analytic models (SAM) of galaxy formation and evolution. New features have been implemented in the reference SAM in order to have it run self-consistently and calibrated on these cosmological simulations. They include an appropriate modification of the mass temperature relation and of the baryon fraction in DM haloes, due to the different virial scalings and to the gravitational bias, respectively. Our results show that the predictions of our coupled-DE SAM do not differ significantly from theoretical predictions obtained with standard SAMs applied to a reference $\Lambda$CDM simulation, implying that the statistical properties of galaxies provide only a weak probe for these alternative cosmological models. On the other hand, we show that both galaxy bias and the galaxy pairwise velocity distribution are sensitive to coupled DE models: this implies that these probes might be successfully applied to disentangle among quintessence, $f(R)$-Gravity and coupled DE models.

We present the first numerical simulations in coupled dark energy cosmologies with high enough resolution to investigate the effects of the coupling on galactic and sub-galactic scales. We choose two constant couplings and a time-varying coupling function and we run simulations of three Milky-Way-size halos ($\sim$10$^{12}$M$_{\odot}$), a lower mass halo (6$\times$10$^{11}$M$_{\odot}$) and a dwarf galaxy halo (5$\times$10$^{9}$M$_{\odot}$). We resolve each halo with several millions dark matter particles. On all scales the coupling causes lower halo concentrations and a reduced number of substructures with respect to LCDM. We show that the reduced concentrations are not due to different formation times, but they are related to the extra terms that appear in the equations describing the gravitational dynamics. On the scale of the Milky Way satellites, we show that the lower concentrations can help in reconciling observed and simulated rotation curves, but the coupling values necessary to have a significant difference from LCDM are outside the current observational constraints. On the other hand, if other modifications to the standard model allowing a higher coupling (e.g. massive neutrinos) are considered, coupled dark energy can become an interesting scenario to alleviate the small-scale issues of the LCDM model.

We study the gravitational lensing effect on the Cosmic Microwave Background (CMB) anisotropies performing a ray-tracing of the primordial CMB photons through intervening large-scale structures (LSS) distribution predicted by N-Body numerical simulations with a particular focus on the precise recovery of the lens-induced polarized counterpart of the source plane. We apply both a multiple plane ray-tracing and an effective deflection approach based on the Born approximation to deflect the CMB photons trajectories through the simulated lightcone. We discuss the results obtained with both these methods together with the impact of LSS non-linear evolution on the CMB temperature and polarization power spectra. We compare our results with semi-analytical approximations implemented in Boltzmann codes like, e.g., CAMB. We show that, with our current N-body setup, the predicted lensing power is recovered with good accuracy in a wide range of multipoles while excess power with respect to semi-analytic prescriptions is observed in the lensing potential on scales $\ell \gtrsim 3000$. We quantify the impact of the numerical effects connected to the resolution in the N-Body simulation together with the resolution and band-limit chosen to synthesise the CMB source plane. We found these quantities to be particularly important for the simulation of B-mode polarization power spectrum.

Cosmic voids are effective cosmological probes to discriminate among competing world models. Their identification is generally based on density or geometry criteria that, because of their very nature, are prone to shot noise. We propose two void finders that are based on dynamical criterion to select voids in Lagrangian coordinates and minimise the impact of sparse sampling. The first approach exploits the Zel'dovich approximation to trace back in time the orbits of galaxies located in voids and their surroundings, the second uses the observed galaxy-galaxy correlation function to relax the objects' spatial distribution to homogeneity and isotropy. In both cases voids are defined as regions of the negative velocity divergence, that can be regarded as sinks of the back-in-time streamlines of the mass tracers. To assess the performance of our methods we used a dark matter halo mock catalogue CoDECS, and compared the results with those obtained with the ZOBOV void finder. We find that the void divergence profiles are less scattered than the density ones and, therefore, their stacking constitutes a more accurate cosmological probe. The significance of the divergence signal in the central part of voids obtained from both our finders is 60% higher than for overdensity profiles in the ZOBOV case. The ellipticity of the stacked void measured in the divergence field is closer to unity, as expected, than what is found when using halo positions. Therefore our void finders are complementary to the existing methods, that should contribute to improve the accuracy of void-based cosmological tests.

Dark matter and dark energy are usually assumed to couple only gravitationally. An extension to this picture is to model dark energy as a scalar field coupled directly to cold dark matter. This coupling leads to new physical effects, such as a fifth-force and a time-dependent dark matter particle mass. In this work we examine the impact that coupling has on weak lensing statistics by constructing realistic simulated weak-lensing maps using raytracing techniques through N-body cosmological simulations. We construct maps for different lensing quantities, covering a range of scales from a few arcminutes to several degrees. The concordance $\Lambda$CDM model is compared to different coupled dark energy models, described either by an exponential scalar field potential (standard coupled dark energy scenario) or by a SUGRA potential (bouncing model). We analyse several statistical quantities and our results, with sources at low redshifts are largely consistent with previous work on CMB lensing by Carbone et al., 2013. The most significant differences from the $\Lambda$CDM model are due to the enhanced growth of the perturbations and to the effective friction term in non-linear dynamics. For the most extreme models, we see differences in the power spectra up to 40% compared to the $\Lambda$CDM model. The different time evolution of the linear matter overdensity can account for most of the differences, but when controlling for this using a $\Lambda$CDM model having the same normalization, the overall signal is smaller due to the effect of the friction term appearing in the equation of motion for dark matter particles.

Low energy interactions between particles are often characterised by elastic scattering. Just as electrons undergo Thomson scattering with photons, dark matter particles may experience an analogous form of momentum exchange with dark energy. We investigate the influence such an interaction has on the formation of linear and nonlinear cosmic structure, by running for the first time a suite of N-body simulations with different dark energy equations of state and scattering cross sections. In models where the linear matter power spectrum is suppressed by the scattering, we find that on nonlinear scales the power spectrum is strongly enhanced. This is due to the friction term increasing the efficiency of gravitational collapse, which also leads to a scale-independent amplification of the concentration and mass functions of halos. The opposite trend is found for models characterised by an increase of the linear matter power spectrum normalisation. More quantitatively, we find that power spectrum deviations at nonlinear scales ($k \approx 10\, h/$Mpc) are roughly ten times larger than their linear counterparts, exceeding $100%$ for the largest value of the scattering cross section considered in the present work. Similarly, the concentration-mass relation and the halo mass function show deviations up to $100%$ and $20%$, respectively, over a wide range of masses. Therefore, we conclude that nonlinear probes of structure formation might provide much tighter constraints on the scattering cross section between dark energy and dark matter as compared to the present bounds based on linear observables.

We investigate the possibility of constraining coupled dark energy (cDE) cosmologies using the three-point correlation function (3PCF). Making use of the CoDECS N-body simulations, we study the statistical properties of cold dark matter (CDM) haloes for a variety of models, including a fiducial $\Lambda$CDM scenario and five models in which dark energy (DE) and CDM mutually interact. We measure both the halo 3PCF, $\zeta(\theta)$, and the reduced 3PCF, $Q(\theta)$, at different scales ($2<r\,[$Mpc\h$]<40$) and redshifts ($0\leq z\leq2$). In all cDE models considered in this work, $Q(\theta)$ appears flat at small scales (for all redshifts) and at low redshifts (for all scales), while it builds up the characteristic V-shape anisotropy at increasing redshifts and scales. With respect to the $\Lambda $CDM predictions, cDE models show lower (higher) values of the halo 3PCF for perpendicular (elongated) configurations. The effect is also scale-dependent, with differences between $\Lambda$CDM and cDE models that increase at large scales. We made use of these measurements to estimate the halo bias, that results in fair agreement with the one computed from the two-point correlation function (2PCF). The main advantage of using both the 2PCF and 3PCF is to break the bias$-\sigma_{8}$ degeneracy. Moreover, we find that our bias estimates are approximately independent of the assumed strength of DE coupling. This study demonstrates the power of a higher-order clustering analysis in discriminating between alternative cosmological scenarios, for both present and forthcoming galaxy surveys, such as e.g. BOSS and Euclid.

The energy scale of Dark Energy, $\sim 2 \times 10^{-3}$ eV, is a long way off compared to all known fundamental scales - except for the neutrino masses. If Dark Energy is dynamical and couples to neutrinos, this is no longer a coincidence. The time at which Dark Energy starts to behave as an effective cosmological constant can be linked to the time at which the cosmic neutrinos become nonrelativistic. This naturally places the onset of the Universe's accelerated expansion in recent cosmic history, addressing the why-now problem of Dark Energy. We show that these mechanisms indeed work in the Growing Neutrino Quintessence model - even if the fully nonlinear structure formation and backreaction are taken into account, which were previously suspected of spoiling the cosmological evolution. The attractive force between neutrinos arising from their coupling to Dark Energy grows as large as $10^6$ times the gravitational strength. This induces very rapid dynamics of neutrino fluctuations which are nonlinear at redshift $z \approx 2$. Nevertheless, a nonlinear stabilization phenomenon ensures only mildly nonlinear oscillating neutrino overdensities with a large-scale gravitational potential substantially smaller than that of cold dark matter perturbations. Depending on model parameters, the signals of large-scale neutrino lumps may render the cosmic neutrino background observable.

In a recent work, Baldi et al. highlighted the issue of cosmic degeneracies, consisting in the fact that the standard statistics of the large-scale structure might not be sufficient to conclusively test cosmological models beyond $\Lambda $CDM when multiple extensions of the standard scenario coexist in nature. In particular, it was shown that the characteristic features of an $f(R)$ Modified Gravity theory and of massive neutrinos with an appreciable total mass $\Sigma _{i}m_{\nu _{i}}$ are suppressed in most of the basic large-scale structure observables for a specific combination of the main parameters of the two non-standard models. In the present work, we explore the possibility that the mean specific size of the supercluster spines -- which was recently proposed as a non-standard statistics by Shim and Lee to probe gravity at large scales -- can help to break this cosmic degeneracy. By analyzing the halo samples from N-body simulations featuring various combinations of $f(R)$ and $\Sigma _{i}m_{\nu _{i}}$ we find that -- at the present epoch -- the value of $\Sigma _{i}m_{\nu _{i}}$ required to maximally suppress the effects of $f(R)$ gravity on the specific sizes of the superclusters spines is different from that found for the other standard statistics. Furthermore, it is also shown that at higher redshifts ($z\ge 0.3$) the deviations of the mean specific sizes of the supercluster spines for all of the four considered combinations from its value for the standard $\Lambda$CDM case are statistically significant.

The recently proposed Multi-coupled Dark Energy (McDE) scenario - characterised by two distinct Cold Dark Matter (CDM) particle species with opposite couplings to a Dark Energy scalar field - introduces a number of novel features in the small-scale dynamics of cosmic structures, most noticeably the simultaneous existence of both attractive and repulsive fifth-forces. Such small-scale features are expected to imprint possibly observable footprints on nonlinear cosmic structures, that might provide a direct way to test the scenario. In order to unveil such footprints, we have performed the first suite of high-resolution N-body simulations of McDE cosmologies, covering the coupling range $|\beta |\leq 1$. We find that for coupling values corresponding to fifth-forces weaker than standard gravity, the impact on structure formation is very mild, thereby showing a new type of screening mechanism for long-range scalar interactions. On the contrary, for fifth-forces comparable to or stronger than standard gravity a number of effects appear in the statistical and structural properties of CDM halos. Collapsed structures start to fragment into pairs of smaller objects that move on different trajectories, providing a direct evidence of the violation of the weak equivalence principle. Consequently, the relative abundance of halos of different masses is significantly modified. For sufficiently large coupling values, the expected number of clusters is strongly suppressed, which might alleviate the present tension between CMB- and cluster-based cosmological constraints. Finally, the internal structure of halos is also modified, with a significant suppression of the inner overdensity, and a progressive segregation of the two CDM species.

We present the first suite of cosmological N-body simulations that simultaneously include the effects of two different and theoretically independent extensions of the standard $\Lambda$CDM cosmological scenario - namely an $f(R)$ theory of Modified Gravity (MG) and a cosmological background of massive neutrinos - with the aim to investigate their possible observational degeneracies. We focus on three basic statistics of the large-scale matter distribution, more specifically the nonlinear matter power spectrum, the halo mass function, and the halo bias, for which we determine the deviation with respect to the fiducial $\Lambda$CDM cosmology in the context of both separate and combined simulations of $f(R)$ MG and massive neutrinos scenarios. Our results show that while these two extended models separately determine very prominent and potentially detectable features in all the three statistics, when we allow them to be simultaneously at work these features are strongly suppressed, resulting in much weaker deviations from the standard model's predictions. In particular, when an $f(R)$ gravity model with $f_{R0}=-1\times 10^{-4}$ is combined with a total neutrino mass of $\Sigma _{i}m_{\nu _{i}}=0.4$ eV, the resulting matter power spectrum, halo mass function, and bias at z=0 are found to be consistent with the standard model's predictions at the 10%, 20%, and 5% accuracy levels, respectively. Therefore, our results imply an intrinsic theoretical limit to the effective discriminating power of present and future observational data sets with respect to these widely considered extensions of the standard cosmological scenario in the absence of independent measurements of the neutrino masses from laboratory experiments, even though the high-redshift evolution might still allow to partially break the degeneracy [Abridged].

The Multi-coupled Dark Energy (McDE) scenario has been recently proposed as a specific example of a cosmological model characterized by a non-standard physics of the dark sector of the universe that nevertheless gives an expansion history which does not significantly differ from the one of the standard $\Lambda $CDM model. In this work, we present the first constraints on the McDE scenario obtained by comparing the predicted evolution of linear density perturbations with a large compilation of recent data sets for the growth rate $f\sigma_{8}$, including 6dFGS, LRG, BOSS, WiggleZ and VIPERS. Confirming qualitative expectations, growth rate data provide much tighter bounds on the model parameters as compared to the extremely loose bounds that can be obtained when only the background expansion history is considered. In particular, the $95\%$ confidence level on the coupling strength $|\beta |$ is reduced from $|\beta |\leq 83$ (background constraints only) to $|\beta |\leq 0.88$ (background and linear perturbation constraints). We also investigate how these constraints further improve when using data from future wide-field surveys such as supernova data from LSST and growth rate data from Euclid-type missions. In this case the $95\%$ confidence level on the coupling further reduce to $|\beta |\leq 0.85$. Such constraints are in any case still consistent with a scalar fifth-force of gravitational strength, and we foresee that tighter bounds might be possibly obtained from the investigation of nonlinear structure formation in McDE cosmologies.[Abridged]

We produce lensing potential and deflection-angle maps in order to simulate CMB weak-lensing via ray-tracing through the COupled Dark Energy Cosmological Simulations (CoDECS). The constructed maps reflect the N-body cosmic structures on a range of scales going from the arcminute to the degree scale. We investigate the variation of the lensing pattern due to the DE dynamics, characterised by different background and perturbation behaviours as a consequence of the interaction between the DE field and CDM. We study the results from three models differing in the background and perturbations evolution, with the purpose to isolate their imprints in the lensing observables. The scenarios investigated include a reference LCDM cosmology, a standard coupled DE (cDE) scenario, and a "bouncing" cDE scenario. For the standard cDE scenario, we find that differences in the lensing potential result from two effects: the enhanced growth of linear CDM density fluctuations with respect to the LCDM case, and the modified nonlinear dynamics of collapsed structures induced by the DE-CDM interaction. As a consequence, CMB lensing highlights the DE impact in the cosmological expansion, even in the degenerate case where the amplitude of the linear matter density perturbations, parametrised through sigma_8, is the same in both the standard cDE and LCDM cosmologies. For the bouncing scenario, we find that the two opposite behaviours of the lens density contrast and of the matter abundance lead to a counter-intuitive effect, making the power of the lensing signal lower by 10% than in the LCDM scenario. Moreover, we compare the behaviour of CDM and baryons in CoDECS separately, in order to isolate effects coming from the coupling with the DE component. We find that, in the bouncing scenario, baryons show an opposite trend with respect to CDM, due to the coupling of the latter with the DE component. [abridged]

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.

The persisting consistency of ever more accurate observational data with the predictions of the standard LCDM cosmological model puts severe constraints on possible alternative scenarios, but still does not shed any light on the fundamental nature of the cosmic dark sector.As large deviations from a LCDM cosmology are ruled out by data, the path to detect possible features of alternative models goes necessarily through the definition of cosmological scenarios that leave almost unaffected the background and -- to a lesser extent -- the linear perturbations evolution of the universe. In this context,the Multi-coupled DE (McDE) model was proposed by Baldi 2012 as a particular realization of an interacting Dark Energy field characterized by an effective screening mechanism capable of suppressing the effects of the coupling at the background and linear perturbation level. In the present paper, for the first time, we challenge the McDE scenario through a direct comparison with real data, in particular with the luminosity distance of Type Ia supernovae. By studying the existence and stability conditions of the critical points of the associated background dynamical system, we select only the cosmologically consistent solutions, and confront their background expansion history with data. Confirming previous qualitative results, the McDE scenario appears to be fully consistent with the adopted sample of Type Ia supernovae, even for coupling values corresponding to an associated scalar fifth-force about four orders of magnitude stronger than standard gravity. Our analysis demonstrates the effectiveness of the McDE background screening, and shows some new non-trivial asymptotic solutions for the future evolution of the universe. Our results show how the background expansion history might be highly insensitive to the fundamental nature and to the internal complexity of the dark sector. [Abridged]

We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation solver on an adaptive mesh to efficiently solve for perturbations in the scalar degree of freedom of the modified gravity model. As this new algorithm is implemented as a module for the P-Gadget3 code, it can at the same time follow the baryonic physics included in P-Gadget3, such as hydrodynamics, radiative cooling and star formation. We demonstrate that the code works reliably by applying it to simple test problems that can be solved analytically, as well as by comparing cosmological simulations to results from the literature. Using the new code, we perform the first non-radiative and radiative cosmological hydrodynamical simulations of an f(R)-gravity model. We also discuss the impact of AGN feedback on the matter power spectrum, as well as degeneracies between the influence of baryonic processes and modifications of gravity.

In this talk, I have discussed the implications of a multi-component nature of cosmic Dark Matter for the observational bounds on possible long-range fifth-forces mediated by a Dark Energy scalar field. By assuming a simple internal symmetry of the Dark Matter component associated to opposite coupling "charges" of two different particle species, the effects of Dark Energy interactions on both the background and linear perturbations evolution are strongly suppressed during the whole matter dominated phase, thereby relaxing present bounds on the coupling strength. The associated attractive and repulsive fifth-forces, however, might still have a very significant impact on the nonlinear dynamics of collapsed structures. I have also described how some of these nonlinear effects are identified through dedicated cosmological N-body simulations as i) a possible fragmentation of bound Dark Matter halos into smaller objects, and ii) a consequent suppression of the nonlinear matter power at small scales. Both effects are potentially observable and might allow to further constrain the model.

The baryonic acoustic peak in the correlation function of galaxies and galaxy clusters provides a standard ruler to probe the space-time geometry of the Universe, jointly constraining the angular diameter distance and the Hubble expansion rate. Moreover, non-linear effects can systematically shift the peak position, giving us the opportunity to exploit this clustering feature also as a dynamical probe. We investigate the possibility of detecting interactions in the dark sector through an accurate determination of the baryonic acoustic scale. Making use of the public halo catalogues extracted from the CoDECS simulations -- the largest suite of N-body simulations of interacting dark energy models to date -- we determine the position of the baryonic scale fitting a band-filtered correlation function, specifically designed to amplify the signal at the sound horizon. We analyze the shifts due to non-linear dynamics, redshift-space distortions and Gaussian redshift errors, in the range 0 < z < 2. Since the coupling between dark energy and dark matter affects in a particular way the clustering properties of haloes and, specifically, the amplitude and location of the baryonic acoustic oscillations, the cosmic evolution of the baryonic peak position might provide a direct way to discriminate interacting dark energy models from the standard \Lambda CDM framework. To maximize the efficiency of the baryonic peak as a dynamic probe, the correlation function has to be measured in redshift-space, where the baryonic acoustic shift due to non-linearities is amplified. The typical redshift errors of spectroscopic galaxy surveys do not significantly impact these results.

(Abridged) The growing role played by numerical N-body simulations in cosmological studies as a fundamental connection between theoretical modeling and direct observations has led to impressive advancements also in the development and application of specific algorithms designed to probe a wide range of Dark Energy scenarios. Over the last decade, a large number of independent and complementary investigations have been carried out in the field of Dark Energy N-body simulations, starting from the simplest case of homogeneous Dark Energy models up to the recent development of highly sophisticated iterative solvers for a variety of Modified Gravity theories. In this Review - which is meant to be complementary to the general Review by Kuhlen et al. published in this Volume - I will discuss the range of scenarios for the cosmic acceleration that have been successfully investigated by means of dedicated N-body simulations, and I will provide a broad summary of the main results that have been obtained in this rather new research field. I will focus the discussion on a few selected studies that have led to particularly significant advancements in the field, and I will provide a comprehensive list of references for a larger number of related works. Due to the vastness of the topic, the discussion will not enter into the finest details of the different implementations and will mainly focus on the outcomes of the various simulations studies. Although quite recent, the field of Dark Energy simulations has witnessed huge developments in the last few years, and presently stands as a reliable approach to the investigation of the fundamental nature of Dark Energy.

We present non-linear weak lensing predictions for coupled dark energy models using the CoDECS simulations. We calculate the shear correlation function and error covariance expected for these models, for forthcoming ground-based (such as DES) and space-based (Euclid) weak lensing surveys. We obtain predictions for the discriminatory power of a ground-based survey similar to DES and a space-based survey such as Euclid in distinguishing between $\Lambda$CDM and coupled dark energy models; we show that using the non-linear lensing signal we could discriminate between $\Lambda$CDM and exponential constant coupling models with $\beta_0\geq0.1$ at $4\sigma$ confidence level with a DES-like survey, and $\beta_0\geq0.05$ at $5\sigma$ confidence level with Euclid. We also demonstrate that estimating the coupled dark energy models' non-linear power spectrum, using the $\Lambda$CDM Halofit fitting formula, results in biases in the shear correlation function that exceed the survey errors.

(Abridged) An interaction between Cold Dark Matter (CDM) and a classical scalar field playing the role of the cosmic dark energy (DE) might provide long-range dark interactions without conflicting with solar system bounds. Although presently available observations allow to constrain such interactions to a few percent of the gravitational strength, some recent studies have shown that if CDM is composed by two different particle species having opposite couplings to the DE field, such tight constraints can be considerably relaxed, allowing for long-range scalar forces of order gravity without significantly affecting observations both at the background and at the linear perturbations level. In the present work, we extend the investigation of such Multiple Dark Matter scenarios to the nonlinear regime of structure formation, by presenting the first N-body simulations ever performed for these cosmologies. Our results highlight some characteristic footprints of long-range scalar forces that arise only in the nonlinear regime for specific models that would be otherwise practically indistinguishable from the standard LCDM scenario both in the background and in the growth of linear density perturbations. Among these effects, the formation of "mirror" cosmic structures in the two CDM species, the suppression of the nonlinear matter power spectrum at k > 1 h/Mpc, and the fragmentation of collapsed halos, represent peculiar features that might provide a direct way to constrain this class of cosmological models.

We present a detailed investigation of the effects that a direct interaction between Dark Energy (DE) and Cold Dark Matter (CDM) particles imprints on the Halo Mass Function (HMF) of groups and clusters of galaxies. Making use of the public halo catalogs of the {\small CoDECS} simulations, we derive the HMF for several different types of coupled DE scenarios both based on the FoF algorithm and on the SO halo identification for different values of the overdensity threshold $\Delta_{c}$. We compare the computed HMFs for coupled DE cosmologies with $\Lambda $CDM as well as with the predictions of the standard analytic fitting functions. Our results show that the standard fitting functions still reproduce reasonably well both the FoF and the SO HMFs of interacting DE cosmologies at intermediate masses and at low redshifts, once rescaled to the characteristic amplitude of linear density perturbations of each specific model as given by $\sigma_{8}$. However, we also find that such apparent degeneracy with $\sigma_{8}$ is broken both by the high-mass tail and by the redshift evolution of our HMFs, with deviations beyond $\sim 10%$ for most of the models under investigation. Furthermore, the discrepancy with respect to the predictions of standard fitting functions rescaled with the characteristic value of $\sigma_{8}$ shows -- for some models -- a strong dependence on the spherical overdensity threshold $\Delta_{c}$ used for the halo identification. We find that such effect is due to a significant increase of halo concentration at low redshifts in these models, that is however absent in the majority of the cosmological scenarios considered in this work. We can therefore conclude that the universality of the HMF is violated by cosmological models that feature a direct interaction between DE and CDM.

(Abridged) Present cosmological constraints and the absence of a direct detection and identification of any dark matter particle candidate leave room to the possibility that the dark sector of the Universe be actually more complex than it is normally assumed. In particular, more than one new fundamental particle could be responsible for the observed dark matter density in the Universe, and possible new interactions between dark energy and dark matter might characterize the dark sector. In the present work, we investigate the possibility that two dark matter particles exist in nature, with identical physical properties except for the sign of their coupling constant to dark energy. Extending previous works on similar scenarios, we study the evolution of the background cosmology as well as the growth of linear density perturbations for a wide range of parameters of such model. Interestingly, our results show how the simple assumption that dark matter particles carry a "charge" with respect to their interaction with the dark energy field allows for new long-range scalar forces of gravitational strength in the dark sector without conflicting with present observations both at the background and linear levels. Our scenario does not introduce new parameters with respect to the case of a single dark matter species for which such strong dark interactions have been already ruled out.

We present the largest set of N-body and hydrodynamical simulations to date for cosmological models featuring a direct interaction between the Dark Energy (DE) scalar field, responsible of the observed cosmic acceleration, and the Cold Dark Matter (CDM) fluid. With respect to previous works, our simulations considerably extend the statistical significance of the simulated volume and cover a wider range of different realizations of the interacting DE scenario, including the recently proposed bouncing coupled DE model. Furthermore, all the simulations are normalized in order to be consistent with the present bounds on the amplitude of density perturbations at last scattering, thereby providing the first realistic determination of the effects of a DE coupling for cosmological growth histories fully compatible with the latest Cosmic Microwave Background data. As a first basic analysis, we have studied the impact of the coupling on the nonlinear matter power spectrum and on the bias between the CDM and baryon distributions, as a function of redshift and scale. For the former, we have addressed the issue of the degeneracy between the effects of the coupling and other standard cosmological parameters, as e.g sigma_8, showing how the redshift evolution of the linear amplitude or the scale dependence of the nonlinear power spectrum might provide a way to break the degeneracy. For the latter, instead, we have computed the redshift and scale dependence of the bias in all our different models showing how a growing coupling or a bouncing coupled DE scenario provide much stronger effects with respect to constant coupling models. We refer to this vast numerical initiative as the COupled Dark Energy Cosmological Simulations project, or CoDECS, and we hereby release all the CoDECS outputs for public use through a dedicated web database, providing information on how to access and interpret the data.

We investigate the spatial properties of the large scale structure (LSS) of the Universe in the framework of coupled dark energy (cDE) cosmologies. Using the public halo catalogues from the CoDECS simulations -- the largest set of N-body experiments to date for such cosmological scenarios -- we estimate the clustering and bias functions of cold dark matter (CDM) haloes, both in real- and redshift-space. Moreover, we investigate the effects of the dark energy (DE) coupling on the geometric and dynamic redshift-space distortions, quantifying the difference with respect to the concordance LambdaCDM model. At z~0, the spatial properties of CDM haloes in cDE models appear very similar to the LambdaCDM case, even if the cDE models are normalized at last scattering in order to be consistent with the latest Cosmic Microwave Background (CMB) data. At higher redshifts, we find that the DE coupling produces a significant scale-dependent suppression of the halo clustering and bias function. This effect, that strongly depends on the coupling strength, is not degenerate with sigma8 at scales r<5-10 Mpc/h. Moreover, we find that the coupled DE strongly affects both the linear distortion parameter, beta, and the pairwise peculiar velocity dispersion, sigma12. Although the models considered in this work are found to be all in agreement with presently available observational data, the next generation of galaxy surveys will be able to put strong constraints on the level of coupling between DE and CDM exploiting the shape of redshift-space clustering anisotropies.

[Abridged] High-resolution N-body simulations have recently shown that the structural properties of highly nonlinear cosmic structures, as e.g. their average concentration at a given mass, could be significantly modified in the presence of an interaction between Dark Energy and Dark Matter. While a constant interaction strength leads to less concentrated density profiles, a steep growth in time of the coupling function has been shown to determine a large increase of halo concentrations over a wide range of masses, including the typical halos hosting luminous spiral galaxies. This determines a substantial worsening of the "cusp-core" tension arising in the standard $\Lambda $CDM model and provides a direct way to constrain the form of the Dark Energy interaction. In the present paper we make use of the outcomes of some high-resolution N-body simulations of a specific class of interacting Dark Energy models to compare the predicted rotation curves of luminous spiral galaxies forming in these cosmologies against real observational data. Our results show how some specific interacting Dark Energy scenarios featuring a steep growth in time of the coupling function -- which are virtually indistinguishable from LCDM in the background -- cannot fit the observed rotation curves of luminous spiral galaxies and can therefore be ruled out only on the basis of dynamical properties of small-scale structures. Our study is a pilot investigation of the effects of a Dark Energy interaction at small scales, and demonstrates how the dynamical properties of visible galaxies can in some cases provide direct constraints on the nature of Dark Energy.

It has been recently shown that the observed morphological properties of the Bullet Cluster can be accurately reproduced in hydrodynamical simulations only when the infall pairwise velocity V_{c} of the system exceeds 3000km/s (or at least possibly 2500 km/s) at the pair separation of 2R_{vir}, where R_{vir} is the virial radius of the main cluster, and that the probability of finding such a bullet-like system is extremely low in the standard \Lambda CDM cosmology. We suggest here the fifth-force mediated by a coupled Dark Energy (cDE) as a possible velocity-enhancing mechanism and investigate its effect on the infall velocities of the bullet-like systems from the CoDECS (COupled Dark Energy Cosmological Simulations) public database. Five different cDE models are considered: three with constant coupling and exponential potential, one with exponential coupling and exponential potential, and one with constant coupling and supergravity potential. For each model, after identifying the bullet-like systems, we determine the probability density distribution of their infall velocities at the pair separations of (2-3)R_{vir}. Approximating each probability density distribution as a Gaussian, we calculate the cumulative probability of finding a bullet-like system with V_{c}>=3000 km/s or V_{c}>=2500 km/s. Our results show that in all of the five cDE models the cumulative probabilities increase compared to the \Lambda CDM case and that in the model with exponential coupling P(V_{c}>=2500 km/s) exceeds 10^{-4}. The physical interpretations and cosmological implications of our results are provided.

We investigate the non-linear evolution of the matter power spectrum by using a large set of high-resolution N-body/hydrodynamic simulations. The linear matter power in the initial conditions is consistently modified to accommodate warm dark matter particles which induce a small scale cut-off in the power as compared to standard cold dark matter scenarios. The impact of such thermal relics is addressed at small scales with k > 1 h/Mpc and at z < 5, which are particularly important for the next generation of Lyman-alpha forest, weak lensing and galaxy clustering surveys. We quantify the mass and redshift dependence of the warm dark matter non-linear matter power and we provide a fitting formula which is accurate at the ~2% level below z=3 and for masses m_wdm > 0.5 keV. The role of baryonic physics (cooling, star formation and feedback recipes) on the warm dark matter induced suppression is also quantified. Furthermore, we compare our findings with the halo model and show their impact on the cosmic shear power spectra.

The abundance of the most massive objects in the Universe at different epochs is a very sensitive probe of the cosmic background evolution and of the growth history of density perturbations, and could provide a powerful tool to distinguish between a cosmological constant and a dynamical dark energy field. In particular, the recent detection of very massive clusters of galaxies at high redshifts has attracted significant interest as a possible indication of a failure of the standard LCDM model. Several attempts have been made in order to explain such detections in the context of non-Gaussian scenarios or interacting dark energy models, showing that both these alternative cosmologies predict an enhanced number density of massive clusters at high redshifts, possibly alleviating the tension. However, all the models proposed so far also overpredict the abundance of massive clusters at the present epoch, and are therefore in contrast with observational bounds on the low-redshift halo mass function. In this paper we present for the first time a new class of interacting dark energy models that simultaneously account for an enhanced number density of massive clusters at high redshifts and for both the standard cluster abundance at the present time and the standard power spectrum normalization at CMB. The key feature of this new class of models is the "bounce" of the dark energy scalar field on the cosmological constant barrier at relatively recent epochs. We present the background and linear perturbations evolution of the model, showing that the standard amplitude of density perturbations is recovered both at CMB and at the present time, and we demonstrate by means of large N-body simulations that our scenario predicts an enhanced number of massive clusters at high redshifts without affecting the present halo abundance. (Abridged)

Growing Neutrino quintessence describes a form of dynamical dark energy that could explain why dark energy dominates the universe only in recent cosmological times. This scenario predicts the formation of large scale neutrino lumps which could allow for observational tests. We perform for the first time N-body simulations of the nonlinear growth of structures for cold dark matter and neutrino fluids in the context of Growing Neutrino cosmologies. Our analysis shows a pulsation - increase and subsequent decrease - of the neutrino density contrast. This could lead to interesting observational signatures, as an enhanced bulk flow in a situation where the dark matter density contrast only differs very mildly from the standard LCDM scenario. We also determine for the first time the statistical distribution of neutrino lumps as a function of mass at different redshifts. Such determination provides an essential ingredient for a realistic estimate of the observational signatures of Growing Neutrino cosmologies. Due to a breakdown of the non-relativistic Newtonian approximation our results are limited to redshifts z > 1.

We present the results of the first N-body simulations of the Growing Neutrino scenario, as recently discussed in Baldi et al. (2011). Our results have shown for the first time how neutrino lumps forming in the context of Growing Neutrino cosmologies are expected to pulsate as a consequence of the rapid oscillations of the dark energy scalar field. We have also computed for the first time a realistic statistical distribution of neutrino halos and determined their impact on the underlying Cold Dark Matter structures.

We investigate the effects of a coupled Dark Energy (cDE) scalar field on the alignment between satellites and matter distributions in galaxy clusters. Using high-resolution N-body simulations for LCDM and cDE cosmological models, we compute the probability density distribution for the alignment angle between the satellite galaxies and underlying matter distributions, finding a difference between the two scenarios. With respect to LCDM, in cDE cosmologies the satellite galaxies are less preferentially located along the major axis of the matter distribution, possibly reducing the tension with obersevational data. A physical explanation is that the coupling between dark matter and dark energy acts as an additional tidal force on the satellite galaxies diminishing the alignments between their distribution and the matter one. Through a Wald test based on the generalized chi-square statistics, the null hypothesis that the two probability distributions come from the same parent population is rejected at the 99 % confidence level. It is concluded that the galaxy-matter alignment in clusters may provide a unique probe of dark sector interactions as well as the nature of dark energy.

We present a complete numerical study of cosmological models with a time dependent coupling between the dark energy component driving the present accelerated expansion of the Universe and the Cold Dark Matter (CDM) fluid. Depending on the functional form of the coupling strength, these models show a range of possible intermediate behaviors between the standard LCDM background evolution and the widely studied case of interacting dark energy models with a constant coupling. These different background evolutions play a crucial role in the growth of cosmic structures, and determine strikingly different effects of the coupling on the internal dynamics of nonlinear objects. By means of a suitable modification of the cosmological N-body code GADGET-2 we have performed a series of high-resolution N-body simulations of structure formation in the context of interacting dark energy models with variable couplings. Depending on the type of background evolution, the halo density profiles are found to be either less or more concentrated with respect to LCDM, contrarily to what happens for constant coupling models where concentrations can only decrease. However, for some specific choice of the interaction function the reduction of halo concentrations can be larger than in constant coupling scenarios. In general, we find that time dependent interactions between dark energy and CDM can in some cases determine stronger effects on structure formation as compared to the constant coupling case, with a significantly weaker impact on the background evolution of the Universe, and might therefore provide a more viable possibility to alleviate the tensions between observations and the LCDM model on small scales than the constant coupling scenario. [Abridged]

The recent detection (Jee etal 2009) of the massive cluster XMMU J2235.3-2557 at a redshift z = 1.4, with an estimated mass M = 6.4 +- 1.2 X 10^14 M_sol, has been claimed to be a possible challenge to the standard LCDM cosmological model. More specifically, the probability to detect such a cluster has been estimated to be 0.005 if a LCDM model with gaussian initial conditions is assumed, resulting in a 3 sigma discrepancy from the standard cosmological model. In this paper we propose to use high redshift clusters as the one detected in Jee etal 2009 to compare the cosmological constant scenario with interacting dark energy models. We show that coupled dark energy models, where an interaction is present between dark energy and cold dark matter, can significantly enhance the probability to observe very massive clusters at high redshift.

We present an analysis of high-resolution hydrodynamical N-body simulations of coupled dark energy cosmologies which focusses on the statistical properties of the transmitted Lyman-alpha flux in the high-redshift intergalactic medium (IGM). In these models the growth of the diffuse cosmic web differs from the standard LCDM case: the density distribution is skewed towards underdense regions and the matter power spectra are typically larger (in a scale dependent way). These differences are also appreciable in the Lyman-alpha flux and are larger than 5% (10%) at z=2-4 in the flux probability distribution function (pdf) for high transmissivity regions and for values of the coupling parameter \beta = 0.08 (\beta = 0.2). The flux power spectrum is also affected at the ~2% (~ 5-10%) level for \beta = 0.08 (\beta = 0.2) in a redshift dependent way. We infer the behaviour of flux pdf and flux power for a reasonable range of couplings and present constraints using present high and low resolution data sets. We find an upper limit \beta < 0.15 (at 2 sigma confidence level), which is obtained using only IGM data and is competitive with those inferred from other large scale structure probes.

We present a detailed numerical study of the impact that cosmological models featuring a direct interaction between the Dark Energy component that drives the accelerated expansion of the Universe and Cold Dark Matter can have on the linear and nonlinear stages of structure formation. By means of a series of collisionless N-body simulations we study the influence that each of the different effects characterizing these cosmological models - which include among others a fifth force, a time variation of particle masses, and a velocity-dependent acceleration - separately have on the growth of density perturbations and on a series of observable quantities related to linear and nonlinear cosmic structures, as the matter power spectrum, the gravitational bias between baryons and Cold Dark Matter, the halo mass function and the halo density profiles. We perform our analysis applying and comparing different numerical approaches previously adopted in the literature, and we address the partial discrepancies recently claimed in a similar study by Li & Barrow (2010b) with respect to the first outcomes of Baldi et al. (2010), which are found to be related to the specific numerical approach adopted in the former work. Our results fully confirm the conclusions of Baldi et al. (2010) and show that when linear and nonlinear effects of the interaction between Dark Energy and Cold Dark Matter are properly disentangled, the velocity-dependent acceleration is the leading effect acting at nonlinear scales, and in particular is the most important mechanism in lowering the concentration of Cold Dark Matter halos.

If the accelerated expansion of the Universe at the present epoch is driven by a dark energy scalar field, there may well be a non-trivial coupling between the dark energy and the cold dark matter (CDM) fluid. Such interactions give rise to new features in cosmological structure growth, like an additional long-range attractive force between CDM particles, or variations of the dark matter particle mass with time. We have implemented these effects in the N-body code GADGET-2 and present results of a series of high-resolution N-body simulations where the dark energy component is directly interacting with the cold dark matter. As a consequence of the new physics, CDM and baryon distributions evolve differently both in the linear and in the nonlinear regime of structure formation. Already on large scales a linear bias develops between these two components, which is further enhanced by the nonlinear evolution. We also find, in contrast with previous work, that the density profiles of CDM halos are less concentrated in coupled dark energy cosmologies compared with LCDM, and that this feature does not depend on the initial conditions setup, but is a specific consequence of the extra physics induced by the coupling. Also, the baryon fraction in halos in the coupled models is significantly reduced below the universal baryon fraction. These features alleviate tensions between observations and the LCDM model on small scales. Our methodology is ideally suited to explore the predictions of coupled dark energy models in the fully non-linear regime, which can provide powerful constraints for the viable parameter space of such scenarios.

The evidence in favor of a dark energy component dominating the Universe, and driving its presently accelerated expansion, has progressively grown during the last decade of cosmological observations. If this dark energy is given by a dynamic scalar field, it may also have a direct interaction with other matter fields in the Universe, in particular with cold dark matter. Such interaction would imprint new features on the cosmological background evolution as well as on the growth of cosmic structure, like an additional long-range fifth-force between massive particles, or a variation in time of the dark matter particle mass. We review here the implementation of these new physical effects in the N-body code GADGET-2, and we discuss the outcomes of a series of high-resolution N-body simulations for a selected family of interacting dark energy models, as already presented in Baldi et al. [20]. We interestingly find, in contrast with previous claims, that the inner overdensity of dark matter halos decreases in these models with respect to LCDM, and consistently halo concentrations show a progressive reduction for increasing couplings. Furthermore, the coupling induces a bias in the overdensities of cold dark matter and baryons that determines a decrease of the halo baryon fraction below its cosmological value. These results go in the direction of alleviating tensions between astrophysical observations and the predictions of the LCDM model on small scales, thereby opening new room for coupled dark energy models as an alternative to the cosmological constant.

We investigate quintessence cosmologies with a matter component consisting of particles with an increasing mass. While negligible in early cosmology, the appearance of a growing matter component has stopped the evolution of the cosmon field at a redshift around six. In turn, this has triggered the accelerated expansion of the Universe. We propose to associate growing matter with neutrinos. Then the presently observed dark energy density and its equation of state are determined by the neutrino mass.

Inflation may have been driven by a component which violated the Null-Energy Condition, thereby leading to {\it super inflation}. We provide the formalism to study cosmological perturbations when such a component is described by a scalar field with arbitrary Lagrangian. Since the background curvature grows with time, gravitational waves always have a blue spectrum. Scalar perturbations may also have a blue spectrum, {\em albeit} in single field models. We apply our formalism to the case of phantom inflation with an exponential potential (whose pole-like inflationary stage is an attractor for inhomogeneous cosmological models for any value of the potential slope). We finally compare the predictions of super inflation with those of standard inflation stressing the role of gravitational waves.