The observed galaxy distribution via galaxy redshift surveys appears distorted due to redshift-space distortions (RSD). While one dominant contribution to RSD comes from the Doppler effect induced by the peculiar velocity of galaxies, the relativistic effects, including the gravitational redshift effect, are recently recognized to give small but important contributions. Such contributions lead to an asymmetric galaxy clustering along the line of sight, and produce non-vanishing odd multipoles when cross-correlating between different biased objects. However, non-zero odd multipoles are also generated by the Doppler effect beyond the distant-observer approximation, known as the wide-angle effect, and at quasi-linear scales, the interplay between wide-angle and relativistic effects becomes significant. In this paper, based on the formalism developed by Taruya et al., we present a quasi-linear model of the cross-correlation function taking a proper account of both the wide-angle and gravitational redshift effects, as one of the major relativistic effects. Our quasi-linear predictions of the dipole agree well with simulations even at the scales below $20\,h^{-1}\,$Mpc, where non-perturbative contributions from the halo potential play an important role, flipping the sign of the dipole amplitude. When increasing the bias difference and redshift, the scale where the sign flip happens is shifted to a larger scale. We derive a simple approximate formula to quantitatively account for the behaviors of the sign flip.

Using the motion of accreting particles onto halos in cosmological $N$-body simulations, we study the radial phase-space structures of cold dark matter (CDM) halos. In CDM cosmology, formation of virialized halos generically produces radial caustics, followed by multi-stream flows of accreted dark matter inside the halos. In particular, the radius of the outermost caustic called the splashback radius exhibits a sharp drop in the slope of the density profile. Here, we focus on the multi-stream structure of CDM halos inside the splashback radius. To analyze this, we use and extend the SPARTA algorithm developed by Diemer. By tracking the particle trajectories accreting onto the halos, we count their number of apocenter passages, which is then used to reveal the multi-stream flows of the dark matter particles. The resultant multi-stream structure in radial phase space is compared with the prediction of the self-similar solution by Fillmore & Goldreich for each halo. We find that $\sim30\%$ of the simulated halos satisfy our criteria to be regarded as being well fitted to the self-similar solution. The fitting parameters in the self-similar solution characterizes physical properties of the halos, including the mass accretion rate and the size of the outermost caustic (i.e., the splashback radius). We discuss in detail the correlation of these fitting parameters and other measures directly extracted from the $N$-body simulation.

The late integrated Sachs-Wolfe (ISW) effect correlates the Cosmic Microwave Background (CMB) temperature anisotropies with foreground cosmic large-scale structures. As the correlation depends crucially on the growth history in the era of dark energy, it is a key observational probe for constraining the cosmological model. Here we present a detailed study based on full-sky and deep light cones generated from very large volume numerical N-body simulations, which allow us to avoid the use of standard replica techniques, while capturing the entirety of the late ISW effect on the large scales. We post-process the light cones using an accurate ray-tracing method and construct full-sky maps of the ISW temperature anisotropy for three different dark energy models. We quantify in detail the extent to which the ISW effect can be used to discriminate between different dark energy scenarios when cross-correlated with the matter distribution or the CMB lensing potential. We also investigate the onset of non-linearities, the so-called Rees-Sciama effect which provides a complementary probe of the dark sector. We find the signal of the lensing-lensing and ISW-lensing correlation of the three dark energy models to be consistent with measurements from the Planck satellite. Future surveys of the large-scale structures may provide cross-correlation measurements that are sufficiently precise to distinguish the signal of these models. Our methodology is very general and can be applied to any dark energy or modified gravity scenario as long as the metric seen by photons can still be characterized by a Weyl potential.

Redshift-space distortions (RSD) in galaxy redshift surveys generally break both the isotropy and homogeneity of galaxy distribution. While the former aspect is particularly highlighted as a probe of growth of structure induced by gravity, the latter aspect, often quoted as wide-angle RSD but ignored in most of the cases, will become important and critical to account for as increasing the statistical precision in next-generation surveys. However, the impact of wide-angle RSD has been mostly studied using linear perturbation theory. In this paper, employing the Zel'dovich approximation, i.e., first-order Lagrangian perturbation theory for gravitational evolution of matter fluctuations, we present a quasi-linear treatment of wide-angle RSD, and compute the cross-correlation function. The present formalism consistently reproduces linear theory results, and can be easily extended to incorporate relativistic corrections (e.g., gravitational redshift).

We simulate the flux emitted from galaxy halos in order to quantify the brightness of the circumgalactic medium (CGM). We use dedicated zoom-in cosmological simulations with the hydrodynamical Adaptive Mesh Refinement code RAMSES, which are evolved down to z=0 and reach a maximum spatial resolution of 380 $h^{-1}$pc and a gas mass resolution up to 1.8$\times 10^{5} h^{-1} \rm{M}_{\odot}$ in the densest regions. We compute the expected emission from the gas in the CGM using CLOUDY emissivity models for different lines (e.g. Ly$\alpha$, CIV, OVI, CVI, OVIII) considering UV background fluorescence, gravitational cooling and continuum emission. In the case of Ly$\alpha$ we additionally consider the scattering of continuum photons. We compare our predictions to current observations and find them to be in good agreement at any redshift after adjusting the Ly$\alpha$ escape fraction. We combine our mock observations with instrument models for FIREBall-2 (UV balloon spectrograph) and HARMONI (visible and NIR IFU on the ELT) to predict CGM observations with either instrument and optimise target selections and observing strategies. Our results show that Ly$\alpha$ emission from the CGM at a redshift of 0.7 will be observable with FIREBall-2 for bright galaxies (NUV$\sim$18 mag), while metal lines like OVI and CIV will remain challenging to detect. HARMONI is found to be well suited to study the CGM at different redshifts with various tracers.

The dark matter halo sparsity provides a direct observational proxy of the halo mass profile, characterizing halos in terms of the ratio of masses within radii which enclose two different overdensities. Previous numerical simulation analyses have shown that at a given redshift the halo sparsity carries cosmological information encoded in the halo mass profile. Moreover, its ensemble averaged value can be inferred from prior knowledge of the halo mass function at the overdensities of interest. Here, we present a detailed study of the ensemble average properties of the halo sparsity. In particular, using halo catalogs from high-resolution N-body simulations, we show that its ensemble average value can be estimated from the ratio of the averages of the inverse halo masses as well as the ratio of the averages of the halo masses at the overdensity of interests. This can be relevant for galaxy clusters data analyses. As an example, we have estimated the average sparsity properties of galaxy clusters from the LoCuSS and HIFLUGCS datasets respectively. The results suggest that the expected consistency of the different average sparsity estimates can provide a test of the robustness of mass measurements in galaxy cluster samples.

The apparent distribution of large-scale structures in the universe is sensitive to the velocity/potential of the sources as well as the potential along the line-of-sight through the mapping from real space to redshift space (redshift-space distortions, RSD). Since odd multipoles of the halo cross-correlation function vanish when considering standard Doppler RSD, the dipole is a sensitive probe of relativistic and wide-angle effects. We build a catalogue of ten million haloes (Milky-Way size to galaxy-cluster size) from the full-sky light-cone of a new "RayGalGroupSims" N-body simulation which covers a volume of ($2.625~h^{-1}$Gpc)$^3$ with $4096^3$ particles. Using ray-tracing techniques, we find the null geodesics connecting all the sources to the observer. We then self-consistently derive all the relativistic contributions (in the weak-field approximation) to RSD: Doppler, transverse Doppler, gravitational, lensing and integrated Sachs-Wolfe. It allows us, for the first time, to disentangle all contributions to the dipole from linear to non-linear scales. At large scale, we recover the linear predictions dominated by a contribution from the divergence of neighbouring line-of-sights. While the linear theory remains a reasonable approximation of the velocity contribution to the dipole at non-linear scales it fails to reproduce the potential contribution below $30-60~h^{-1}$Mpc (depending on the halo mass). At scales smaller than $\sim 10~h^{-1}$Mpc, the dipole is dominated by the asymmetry caused by the gravitational redshift. The transition between the two regimes is mass dependent as well. We also identify a new non-trivial contribution from the non-linear coupling between potential and velocity terms.

The environmental dependence of galaxy clustering encodes information about the physical processes governing the growth of cosmic structure. We analyze the baryon acoustic peak as a function of environment in the galaxy correlation function of the Baryon Oscillation Spectroscopic Survey CMASS sample. Dividing the sample into three subsets by smoothed local overdensity, we detect acoustic peaks in the six separate auto-correlation and cross-correlation functions of the sub-samples. Fitting models to these correlation functions, calibrated by mock galaxy and dark matter catalogues, we find that the inferred distance scale is independent of environment, and consistent with the result of analyzing the combined sample. The shape of the baryon acoustic feature, and the accuracy of density-field reconstruction in the Zeldovich approximation, varies with environment. By up-weighting underdense regions and down-weighting overdense regions in their contribution to the full-sample correlation function, by up to 50%, we achieve a fractional improvement of a few per cent in the precision of baryon acoustic oscillation fits to the CMASS data and mock catalogues: the scatter in the preferred-scale fits to the ensemble of mocks improves from 1.45% to 1.34% (pre-reconstruction) and 1.03% to 1.00% (post-reconstruction). These results are consistent with the notion that the acoustic peak is sharper in underdense environments.

We present a new cosmological probe for galaxy clusters, the halo sparsity. This characterises halos in terms of the ratio of halo masses measured at two different radii and carries cosmological information encoded in the halo mass profile. Building upon the work of Balmes et al. (2014) we test the properties of the sparsity using halo catalogs from a numerical N-body simulation of ($2.6$ Gpc/h)$^3$ volume with $4096^3$ particles. We show that at a given redshift the average sparsity can be predicted from prior knowledge of the halo mass function. This provides a quantitative framework to infer cosmological parameter constraints using measurements of the sparsity of galaxy clusters. We show this point by performing a likelihood analysis of synthetic datasets with no systematics, from which we recover the input fiducial cosmology. We also perform a preliminary analysis of potential systematic errors and provide an estimate of the impact of baryonic effects on sparsity measurements. We evaluate the sparsity for a sample of 104 clusters with hydrostatic masses from X-ray observations and derive constraints on the cosmic matter density $\Omega_m$ and the normalisation amplitude of density fluctuations at the $8$ Mpc h$^{-1}$ scale, $\sigma_8$. Assuming no systematics, we find $\Omega_m=0.42\pm 0.17$ and $\sigma_8=0.80\pm 0.31$ at $1\sigma$, corresponding to $S_8\equiv \sigma_8\sqrt{\Omega_m}=0.48\pm 0.11$. Future cluster surveys may provide opportunities for precise measurements of the sparsity. A sample of a few hundreds clusters with mass estimate errors at a few percent level can provide competitive cosmological parameter constraints complementary to those inferred from other cosmic probes.

The large-scale matter distribution in the late-time Universe exhibits gravity-induced non-Gaussianity, and the bispectrum, three-point cumulant is expected to contain significant cosmological information. In particular, the measurement of the bispectrum helps to tighten the constraints on dark energy and modified gravity through the redshift-space distortions (RSD). In this paper, extending the work by Taruya, Nishimichi & Saito (2010, Phys.Rev.D 82, 063522), we present a perturbation theory (PT) based model of redshift-space matter bispectrum that can keep the non-perturbative damping effect under control. Characterizing this non-perturbative damping by a univariate function with single free parameter, the PT model of the redshift-space bispectrum is tested against a large set of cosmological $N$-body simulations, finding that the predicted monopole and quadrupole moments are in a good agreement with simulations at the scales of baryon acoustic oscillations (well beyond the range of agreement of standard PT). The validity of the univariate ansatz of the damping effect is also examined, and with the PT calculation at next-to-leading order, the fitted values of the free parameter is shown to consistently match those obtained from the PT model of power spectrum by Taruya, Nishimichi & Saito (2010).

We perform a study of the nonlinear clustering of matter in the late-forming dark matter (LFDM) scenario in which dark matter results from the transition of a nonminimally coupled scalar field from radiation to collisionless matter. A distinct feature of this model is the presence of a damped oscillatory cutoff in the linear matter power spectrum at small scales. We use a suite of high-resolution N-body simulations to study the imprints of LFDM on the nonlinear matter power spectrum, the halo mass and velocity functions and the halo density profiles. The model largely satisfies high-redshift matter power spectrum constraints from Lyman-$\alpha$ forest measurements, while it predicts suppressed abundance of low-mass halos ($\sim 10^{9}-10^{10}$ h$^{-1}$ M$_\odot$) at all redshifts compared to a vanilla $\Lambda$CDM model. The analysis of the LFDM halo velocity function shows a better agreement than the $\Lambda$CDM prediction with the observed abundance of low-velocity galaxies in the local volume. Halos with mass $M\gtrsim 10^{11}$ h$^{-1}$ M$_\odot$ show minor departures of the density profiles from $\Lambda$CDM expectations, while smaller-mass halos are less dense, consistent with the fact that they form later than their $\Lambda$CDM counterparts.

In this article we compare the halo mass function predicted by the excursion set theory with a drifting diffusive barrier against the results of N-body simulations for several cosmological models. This includes the standard LCDM case for a large range of halo masses, models with different types of primordial non-Gaussianity, and the Ratra-Peebles quintessence model of Dark Energy. We show that in all those cosmological scenarios, the abundance of dark matter halos can be described by a drifting diffusive barrier, where the two parameters describing the barrier have physical content. In the case of the Gaussian LCDM, the statistics are precise enough to actually predict those parameters at different redshifts from the initial conditions. Furthermore, we found that the stochasticity in the barrier is nonnegligible making the simple deterministic spherical collapse model a bad approximation even at very high halo masses. We also show that using the standard excursion set approach with a barrier inspired by peak patches leads to inconsistent predictions of the halo mass function.

The excursion set approach provides a framework for predicting how the abundance of dark matter halos depends on the initial conditions. A key ingredient of this formalism comes from the physics of halo formation: the specification of a critical overdensity threshold (barrier) which protohalos must exceed if they are to form bound virialized halos at a later time. Another ingredient is statistical, as it requires the specification of the appropriate statistical ensemble over which to average when making predictions. The excursion set approach explicitly averages over all initial positions, thus implicitly assuming that the appropriate ensemble is that associated with randomly chosen positions in space, rather than special positions such as peaks of the initial density field. Since halos are known to collapse around special positions, it is not clear that the physical and statistical assumptions which underlie the excursion set approach are self-consistent. We argue that they are at least for low mass halos, and illustrate by comparing our excursion set predictions with numerical data from the DEUS simulations.

The upcoming generation of galaxy surveys will probe the distribution of matter in the universe with unprecedented accuracy. Measurements of the matter power spectrum at different scales and redshifts will provide stringent constraints on the cosmological parameters. However, on non-linear scales this will require an accurate evaluation of the covariance matrix. Here, we compute the covariance matrix of the 3D matter density power spectrum for the concordance $\Lambda$CDM cosmology from an ensemble of N-body simulations of the Dark Energy Universe Simulation - Parallel Universe Runs (DEUS-PUR). This consists of 12288 realisations of a $(656\,h^{-1}\,\textrm{Mpc})^3$ simulation box with $256^3$ particles. We combine this set with an auxiliary sample of 96 simulations of the same volume with $1024^3$ particles. We find N-body mass resolution effect to be an important source of systematic errors on the covariance at high redshift and small intermediate scales. We correct for this effect by introducing an empirical statistical method which provide an accurate determination of the covariance matrix over a wide range of scales including the Baryon Oscillations interval. Contrary to previous studies that used smaller N-body ensembles, we find the power spectrum distribution to significantly deviate from expectations of a Gaussian random density field at $k\gtrsim 0.25\,h\,\textrm{Mpc}^{-1}$ and $z<0.5$. This suggests that in the case of finite volume surveys an unbiased estimate of the ensemble averaged band power at these scales and redshifts may require a careful assessment of non-Gaussian errors more than previously considered.

Observations of colliding galaxy clusters with high relative velocity probe the tail of the halo pairwise velocity distribution with the potential of providing a powerful test of cosmology. As an example it has been argued that the discovery of the Bullet Cluster challenges standard $\Lambda$CDM model predictions. Halo catalogs from N-body simulations have been used to estimate the probability of Bullet-like clusters. However, due to simulation volume effects previous studies had to rely on a Gaussian extrapolation of the pairwise velocity distribution to high velocities. Here, we perform a detail analysis using the halo catalogs from the Dark Energy Universe Simulation Full Universe Runs (DEUS-FUR), which enables us to resolve the high-velocity tail of the distribution and study its dependence on the halo mass definition, redshift and cosmology. Building upon these results we estimate the probability of Bullet-like systems in the framework of Extreme Value Statistics. We show that the tail of extreme pairwise velocities significantly deviates from that of a Gaussian, moreover it carries an imprint of the underlying cosmology. We find the Bullet Cluster probability to be two orders of magnitude larger than previous estimates, thus easing the tension with the $\Lambda$CDM model. Finally, the comparison of the inferred probabilities for the different DEUS-FUR cosmologies suggests that observations of extreme interacting clusters can provide constraints on dark energy models complementary to standard cosmological tests.

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

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

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

We have performed the first-ever numerical N- body simulation of the full observable universe (DEUS "Dark Energy Universe Simulation" FUR "Full Universe Run"). This has evolved 550 billion particles on an Adaptive Mesh Refinement grid with more than two trillion computing points along the entire evolutionary history of the universe and across 6 order of magnitudes length scales, from the size of the Milky Way to that of the whole observable universe. To date, this is the largest and most advanced cosmological simulation ever run. It provides unique information on the formation and evolution of the largest structure in the universe and an exceptional support to future observational programs dedicated to mapping the distribution of matter and galaxies in the universe. The simulation has run on 4752 (of 5040) thin nodes of BULL supercomputer CURIE, using more than 300 TB of memory for 10 million hours of computing time. About 50 PBytes of data were generated throughout the run. Using an advanced and innovative reduction workflow the amount of useful stored data has been reduced to 500 TBytes.

We present for the first time predictions for UV line emission of intergalactic and circumgalactic gas from Adaptive Mesh Resolution (AMR) Large Scale Structure (LSS) simulations at redshifts 0.3<z<1.2, with specific emphasis on observability with current and near-future UV instrumentation. In three transitions of interest (Lya, OVI and CIV) there is a clear bimodality in the type of objects : the overwhelming majority of flux stems from discrete, compact sources, while a much larger volume fraction is filled by more tenuous gas. We characterise both object types with regard to number density, physical size and shape, brightness, luminosity, velocity structure, mass, temperature, ionisation state, and metal content. Degrading AMR grids to characteristic resolutions of available (such as FIREBall) or foreseeable instrumentation, allows to assess which inferences can be drawn from currently possible observations, and set foundations to prepare observing strategies for future missions. In general, the faint emission of the IGM and filamentary structure remains beyond capabilities of instruments with only short duration exposure potential (stratospheric balloons), even for optimistic assumptions for Lya, while the yet fainter metal line transitions for these structures will remain challenging for long duration exposures (space-based telescopes), mostly due to low metallicity pushing them more than three orders of magnitudes in brightness below Lya radiation. For the circum-galactic medium (CGM) the situation is more promising, and it is foreseeable that in the near future we will not only just dectect such sources, but the combination of all three lines in addition to velocity information will yield valuable insight into the physical processes at hand, illuminating important mechanisms during the formation of galaxies and their backreaction onto the IGM from whence they formed. (abrigded)

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

The universality of the halo mass function is investigated in the context of dark energy cosmologies. This widely used approximation assumes that the mass function can be expressed as a function of the matter density omega_m and the rms linear density fluctuation sigma only, with no explicit dependence on the properties of dark energy or redshift. In order to test this hypothesis we run a series of 15 high-resolution N-body simulations for different cosmological models. These consists of three LCDM cosmologies best fitting WMAP-1, 3 and 5 years data, and three toy-models characterized by a Ratra-Peebles quintessence potential with different slopes and amounts of dark energy density. These toy models have very different evolutionary histories at the background and linear level, but share the same sigma8 value. For each of these models we measure the mass function from catalogues of halos identified in the simulations using the Friend-of-Friend (FoF) algorithm. We find redshift dependent deviations from a universal behaviour, well above numerical uncertainties and of non-stochastic origin, which are correlated with the linear growth factor of the investigated cosmologies. Using the spherical collapse as guidance, we show that such deviations are caused by the cosmology dependence of the non-linear collapse and virialization process. For practical applications, we provide a fitting formula of the mass function accurate to 5 percents over the all range of investigated cosmologies. We also derive an empirical relation between the FoF linking parameter and the virial overdensity which can account for most of the deviations from an exact universal behavior. Overall these results suggest that the halo mass function contains unique cosmological information since it carries a fossil record of the past cosmic evolution.

In this "Invisible Universe" proceedings, we introduce the Dark Energy Universe Simulation Series (DEUSS) which aim at investigating the imprints of realistic dark energy models on cosmic structure formation. It represents the largest dynamical dark energy simulation suite to date in term of spatial dynamics. We first present the 3 realistic dark energy models (calibrated on latest SNIa and CMB data): LambdaCDM, quintessence with Ratra-Peebles potential, and quintessence with Sugra potential. We then isolate various contributions for non-linear matter power spectra from a series of pre-DEUSS high-resolution simulations (130 million particles). Finally, we introduce DEUSS which consist in 9 Grand Challenge runs with 1 billion particles each thus probing scales from 4 Gpc down to 3 kpc at z=0. Our goal is to make these simulations available to the community through the "Dark Energy Universe Virtual Observatory" (DEUVO), and the "Dark Energy Universe Simulations" (DEUS) consortium.

The largest galaxies acquire their mass early on, when the Universe is still youthful. Cold streams violently feed these young galaxies a vast amount of fresh gas, resulting in very efficient star formation. Using a well resolved hydrodynamical simulation of galaxy formation, we demonstrate that these mammoth galaxies are already in place a couple of billion years after the Big Bang. Contrary to local starforming galaxies, where dust re-emits a large part of the stellar ultraviolet (UV) light at infrared and sub-millimetre wavelengths, our self-consistent modelling of dust extinction predicts that a substantial fraction of UV photons should escape from primordial galaxies. Such a model allows us to compute reliably the number of high redshift objects as a function of luminosity, and yields galaxies whose UV luminosities closely match those measured in the deepest observational surveys available. This agreement is remarkably good considering our admittedly still simple modelling of the interstellar medium (ISM) physics. The luminosity functions (LF) of virtual UV luminous galaxies coincide with the existing data over the whole redshift range from 4 to 7, provided cosmological parameters are set to their currently favoured values. Despite their considerable emission at short wavelengths, we anticipate that the counterparts of the brightest UV galaxies will be detected by future sub-millimetre facilities like ALMA

Dark energy as a quintessence component causes a typical modification of the background cosmic expansion, which in addition to its clustering properties, can leave a potentially distinctive signature on large scale structures. Many previous studies have investigated this topic, particularly in relation to the non-linear regime of structure formation. However, no careful pre-selection of viable quintessence models with high precision cosmological data was performed. Here we show that this has led to a misinterpretation (and underestimation) of the imprint of quintessence on the distribution of large scale structures. To this purpose we perform a likelihood analysis of the combined Supernova Ia UNION dataset and WMAP5-years data to identify realistic quintessence models. Differences from the vanilla LambdaCDM are especially manifest in the predicted amplitude and shape of the linear matter power spectrum, though these remain within the uncertainties of the SDSS data. We use these models as benchmark for studying the clustering properties of dark matter halos by performing a series of high resolution N-body simulations. We find that realistic quintessence models allow for relevant differences of the dark matter distribution with the respect to the LambdaCDM scenario well into the non-linear regime, with deviations up to 40% in the non-linear power spectrum. Such differences are shown to depend on the nature of DE, as well as the scale and epoch considered. At small scales (k~1-5 h Mpc^{-1}, depending on the redshift) the structure formation process is about 20% more efficient than in LambdaCDM. We show that these imprints are a specific record of the cosmic structure formation history in DE cosmologies and therefore cannot be accounted in standard fitting functions of the non-linear matter power spectrum.

In clusters of galaxies, the specific entropy of intracluster plasma increases outwards. Nevertheless, a number of recent studies have shown that the intracluster medium is subject to buoyancy instabilities due to the effects of cosmic rays and anisotropic thermal conduction. In this paper, we present a new numerical algorithm for simulating such instabilities. This numerical method treats the cosmic rays as a fluid, accounts for the diffusion of heat and cosmic rays along magnetic field lines, and enforces the condition that the temperature and cosmic-ray pressure remain positive. We carry out several tests to ensure the accuracy of the code, including the detailed matching of analytic results for the eigenfunctions and growth rates of linear buoyancy instabilities. This numerical scheme will be useful for simulating convection driven by cosmic-ray buoyancy in galaxy cluster plasmas and may also be useful for other applications, including fusion plasmas, the interstellar medium, and supernovae remnants.

Central peaks in the iron abundance of intracluster plasma are a common feature of cooling-core galaxy clusters. These abundance peaks have a much broader profile than the stars of the central brightest cluster galaxy (BCG), which produce the excess iron, indicating that metal-enriched plasma is transported out of the BCG by some process such as turbulent diffusion. At the same time, cooling-core clusters are likely heated by central active galactic nuclei (AGNs) by means of jets, cosmic-ray bubbles, and convection. The recent AGN-driven convection model of Chandran&Rasera predicts the turbulent velocity profile in a steady-state cluster in which radiative cooling is balanced by heating from a combination of AGN-driven convection and thermal conduction. We use the velocity profiles as input into an advection/diffusion model for the transport of metals in the intracluster medium, taking the iron to be injected by the BCG. We compare the results of our model to XMM and Chandra observations of eight clusters. Assuming a constant turbulence level over a cluster's lifetime, the turbulent velocities in the model can explain the observed abundance profiles in only five of the eight clusters. However, we go on to develop an analytic fit of the turbulent velocity profile as a function of the AGN power. We then deduce for each cluster the average AGN power (during the past 10 Gyr) required to match the abundance profiles. The required average values are between 10^43 and 2.10^44 erg/s, while the present AGN powers span a much larger range from 6.10^41 (Virgo) to 2.10^44 erg/s (Hydra A). Our results suggest that AGN-driven convection can account for the observed abundance profiles if the AGN power varies over a cluster's lifetime between Virgo-like and Hydra-A-like values, with average values in the above-quoted range.

A number of studies have shown that the convective stability criterion for the intracluster medium (ICM) is very different from the Schwarzchild criterion due to the effects of anisotropic thermal conduction and cosmic rays. Building on these studies, we develop a model of the ICM in which a central active galactic nucleus (AGN) accretes hot intracluster plasma at the Bondi rate and produces cosmic rays that cause the ICM to become convectively unstable. The resulting convection heats the intracluster plasma and regulates its temperature and density profiles. By adjusting a single parameter in the model (the size of the cosmic-ray acceleration region), we are able to achieve a good match to the observed density and temperature profiles in a sample of eight clusters. Our results suggest that convection is an important process in cluster cores. An interesting feature of our solutions is that the cooling rate is more sharply peaked about the cluster center than is the convective heating rate. As a result, in several of the clusters in our sample, a compact cooling flow arises in the central region with a size R that is typically a few kpc. The cooling flow matches onto a Bondi flow at smaller radii. The mass accretion rate in the Bondi flow is equal to, and controlled by, the rate at which mass flows in through the cooling flow. Our solutions suggest that the AGN regulates the mass accretion rate in these clusters by controlling R: if the AGN power rises above the equilibrium level, R decreases, the mass accretion rate drops, and the AGN power drops back down to the equilibrium level.

The bulk of the extragalactic background between 10 keV and 10 GeV is likely to be explained by the emission of Seyfert galaxies, type Ia supernovae, and blazars. However, as revealed by the INTEGRAL satellite, the bulge of our galaxy is an intense source of a 511 keV gamma-ray line, indicating the production of a large number of positrons that annihilate. The origin of the latter is debated, and they could be produced, in particular, by the (S- or P-wave) annihilation of light Dark Matter particles into e+e-. In any case, the cumulated effect of similar sources at all redshifts could lead to a new background of hard X-ray and soft gamma-ray photons. On the basis of the hierarchical model of galaxy formation, we compute analytically the SNIa contribution to the background, and add it to Seyfert and blazars emission models. We find that any extra contribution to this unresolved background at 511 keV should be lower than about 4 keV/cm^2/s/sr. We also estimate analytically the extragalactic background due to Dark Matter annihilation, increasing the accuracy of the earlier computations. Indeed, we take into account the large positron escape fraction from low mass dark matter halos, unable to confine a dense and magnetized interstellar medium. Our new background estimate turns out to be one order of magnitude lower, so that the hypothesis of a light Dark Matter candidate remains compatible with the observed extragalactic background for a wider range of particle masses and cross-sections.

Using a series of high-resolution N-body hydrodynamical numerical simulations, we investigate several scenarios for the evolution of the baryon budget in galactic halos. We derive individual halo star formation history (SFH), as well as the global star formation rate in the universe. We develop a simple analytical model that allows us to compute surprisingly accurate predictions, when compared to our simulations, but also to other simulations presented in Springel & Hernquist (2003). The model depends on two main parameters: the star formation time scale t* and the wind efficiency eta_w. We also compute, for halos of a given mass, the baryon fraction in each of the following phases: cold disc gas, hot halo gas and stars. Here again, our analytical model predictions are in good agreement with simulation results, if one correctly takes into account finite resolution effect. We compare predictions of our analytical model to several observational constraints, and conclude that a very narrow range of the model parameters is allowed. The important role played by galactic winds is outlined, as well as a possible `superwind' scenario in groups and clusters. The `anti-hierarchical' behavior of observed SFH is well reproduced by our best model with t*=3Gyr and eta_w=1.5. We obtain in this case a present-day cosmic baryon budget of Omega*= 0.004, Omega_cold=0.0004, Omega_hot=0.01 and Omega_back=0.02 (diffuse background).