Search

This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of recordMonthly Notices of the Royal Astronomical Society 511.2 (2022): 2897-2913 is available online at: https://academic.oup.com/mnras/article-abstract/511/2/2897/6517700?redirectedFrom=fulltext, Using the data set of the three hundred project, i.e. a suite of 324 hydrodynamical resimulations of cluster-sized haloes, we study galaxy cluster mergers and their effect on colour and luminosity changes of their brightest cluster galaxies (BCG). We track the main progenitor of each halo at z = 0 and search for merger situations based on its mass accretion history, defining mergers as very rapid increases in the halo mass. Based upon the evolution of the dynamical state of the cluster we define a pre-and post-merger phase. We create a list of all these events and statistically study their mass ratio and time-scales, with the former verifying that all instances are in fact major mergers. By comparing to a control sample of clusters without mergers, we study the effect mergers have on the stellar component of the BCG. Analysing the mass, age, and metallicity of the BCG stellar particles, we find that the stellar content of BCGs grows significantly during mergers and, even though the main growth mechanism is the accretion of older stars, there is even a burst in star formation induced by the merger. In our simulations, BCGs in mergers form in median around 70 per cent more stars than those normally growing, although this depends on the radius considered for defining the BCG. Regarding observable properties, we see an increase in SDSS-u luminosity of 20 per cent during mergers, accompanied by a slightly slower increase of the galaxy g-r colour as compared to the control sample

This is a pre-copyedited, author-produced PDF of an article accepted for publication in Monthly Notices of the Royal Astronomical Society following peer review. The version of record Monthly Notices of the Royal Astronomical Society 510.1 (2022): 292-308 is available online at: https://academic.oup.com/mnras/article-abstract/510/1/292/6442269?redirectedFrom=fulltext#no-access-message, We study the clustering of H i intensity maps (IM) produced from simulations with a focus on baryonic acoustic oscillations (BAOs) and the effects induced by telescope beam smoothing and foreground cleaning. We start by creating an H i catalogue at z = 1.321 based on the Semi-Analytic Galaxy Evolution (SAGE) model applied to the UNIT simulations. With this catalogue, we investigate the relation between model H i and the dark matter haloes and we also study the abundance of H i, Omega rm{H small {rm I}}, predicted by this model. We then create synthetic H i IM with a nearest-grid-point approach. In order to simulate the telescope beam effect, a Gaussian smoothing is applied on the plane perpendicular to the line of sight. The effect of foreground removal methods is simulated by exponentially damping the largest wavelength Fourier modes on the radial direction. We study the anisotropic 2-point correlation function (2PCF) ζ(rpdbl, rpdbl) and how it is affected by the aforementioned observational effects. In order to better isolate the BAO signal, we study several 2PCF μ-wedges (with a restricted range of orientations μ) tailored to address the systematics effects and we compare them with different definitions of radial 2PCFs. Finally, we discuss our findings in the context of an SKA-like survey, finding a clear BAO signal in most of the estimators here proposed

Next generation telescopes such as the James Webb Space Telescope (JWST) and the Nancy Grace Roman Space Telescope (NGRST) will enable us to study the first billion years of our Universe in unprecedented detail. In this work we use the ASTRAEUS (semi-numerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework, that couples galaxy formation and reionization (for a wide range of reionization feedback models), to estimate the cosmic variance expected in the UV Luminosity Function (UV LF) and the Stellar Mass Function (SMF) in JWST surveys. We find that different reionization scenarios play a minor role in the cosmic variance. Most of the cosmic variance is completely driven by the underlying density field and increases above $100\%$ for ${\rm M}_{\rm UV} \sim -17.5 (-20)$ at $z =12 (6)$ for the JADES-deep survey (the deep JWST Advanced Extragalactic Survey with an area of 46 arcmin$^2$); the cosmic variance decreases with an increasing survey area roughly independently of redshift. We find that the faint-end (${\rm M}_{\rm UV} > -17$) slope of the Lyman Break Galaxies (LBGs) UV LF becomes increasingly shallower with increasing reionization feedback and show how JWST observations will be able to distinguish between different models of reionization feedback at $z>9$, even accounting for cosmic variance. We also show the environments (in terms of density and ionization fields)of Lyman Break Galaxies during the EoR, finding that the underlying over-density and ionization fraction scale positively with the UV luminosity. Finally, we also provide a public software tool to allow interested readers to compute cosmic variance for different redshifts and survey areas., 16 pages, 13 figures. Accepted for publication in MNRAS

In this work, we use the {\sc astraeus} (seminumerical rAdiative tranSfer coupling of galaxy formaTion and Reionization in N-body dArk mattEr simUlationS) framework which couples galaxy formation and reionization in the first billion years. Exploring a number of models for reionization feedback and the escape fraction of ionizing radiation from the galactic environment ($f_\mathrm{esc}$), we quantify how the contribution of star-forming galaxies {(with halo masses $M_h>10^{8.2}$M$_\odot$)} to reionization depends on the radiative feedback model, $f_\mathrm{esc}$, and the environmental over-density. Our key findings are: (i) for constant $f_\mathrm{esc}$ models, intermediate-mass galaxies (with halo masses of $M_h\simeq10^{9-11}$M$_\odot$ and absolute UV magnitudes of $M_{UV} \sim -15$ to $-20$) in intermediate-density regions drive reionization; (ii) scenarios where $f_\mathrm{esc}$ increases with decreasing halo mass shift the galaxy population driving reionization to lower-mass galaxies ($M_h\lesssim10^{9.5}$M$_\odot$) with lower luminosities ($M_{UV} \gtrsim-16$) and over-densities; (iii) reionization imprints its topology on the ionizing emissivity of low-mass galaxies ($M_h\lesssim10^{9}$M$_\odot$) through radiative feedback. Low-mass galaxies experience a stronger suppression of star formation by radiative feedback and show lower ionizing emissivities in over-dense regions; (iv) a change in $f_\mathrm{esc}$ with galaxy properties has the largest impact on the sources of reionization and their detectability, with the radiative feedback strength and environmental over-density playing a sub-dominant role; (v) JWST-surveys (with a limiting magnitude of $M_{UV} = -16$) will be able to detect the galaxies providing $\sim 60-70\%$ ($\sim 10\%$) of reionization photons at $z=7$ for constant $f_\mathrm{esc}$ models (scenarios where $f_\mathrm{esc}$ increases with decreasing halo mass)., 14 pages, 13 figures, accepted for publication in MNRAS

Associations of dwarf galaxies are loose systems composed exclusively of dwarf galaxies. These systems were identified in the Local Volume for the first time more than 30 yr ago. We study these systems in the cosmological framework of the λ cold dark matter (λCDM) model.We consider the Small MultiDark Planck simulation and populate its dark matter haloes by applying the semi-analytic model of galaxy formation SAG. We identify galaxy systems using a friends-of-friends algorithm with a linking length equal to b = 0.4Mpc h-1to reproduce the size of dwarf galaxy associations detected in the Local Volume. Our samples of dwarf systems are built up removing those systems that have one or more galaxies with stellar mass larger than a maximum thresholdMmax.We analyse three different samples defined by log10(Mmax[M⊙ h-1]) = 8.5, 9.0, and 9.5. On average, our systems have typical sizes of ∼ 0.2Mpc h-1, velocity dispersion of ∼ 30km s-1, and estimated total mass of ∼ 1011M⊙ h-1. Such large typical sizes suggest that individual members of a given dwarf association reside in different dark matter haloes and are generally not substructures of any other halo. Indeed, in more than 90 per cent of our dwarf systems their individual members inhabit different dark matter haloes, while only in the remaining 10 per cent members do reside in the same halo. Our results indicate that the λCDM model can naturally reproduce the existence and properties of dwarf galaxies' associations without much difficulty. Fil: Yaryura, Claudia Yamila. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina Fil: Abadi, Mario Gabriel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina Fil: Gottlöber, Stefan. Leibniz Universitat Hannover; Alemania Fil: Libeskind, Noam I.. Leibniz Universitat Hannover; Alemania Fil: Cora, Sofia Alejandra. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; Argentina Fil: Ruiz, Andrés Nicolás. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Instituto de Astronomía Teórica y Experimental. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba. Instituto de Astronomía Teórica y Experimental; Argentina. Universidad Nacional de Córdoba. Observatorio Astronómico de Córdoba; Argentina Fil: Vega Martínez, Cristian Antonio. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - La Plata. Instituto de Astrofísica La Plata. Universidad Nacional de La Plata. Facultad de Ciencias Astronómicas y Geofísicas. Instituto de Astrofísica La Plata; Argentina. Universidad de La Serena; Chile Fil: Yepes, Gustavo. Universidad Autónoma de Madrid; España Fil: Behroozi, Peter. George Mason University. School Of Physics. Astronomy And Computational Sciences; Estados Unidos

We compare the statistics and morphology of giant arcs in galaxy clusters using N-body and non-radiative SPH simulations within the standard cold dark matter (CDM) model and simulations where dark matter (DM) has a non-negligible probability of interaction (parametrized by its cross-section), i.e self-interacting dark matter (SIDM). We use a ray-tracing technique to produce a statistically large number of arcs around six simulated galaxy clusters at different redshifts. Since DM is more likely to interact in colliding clusters than in relaxed clusters, and this probability of interaction is largest in denser regions, we focus our analysis on radial arcs (which trace the lensing potential in the central region better than tangential arcs) in galaxy clusters that underwent (or are undergoing) a major merger. We find that SIDM produces fewer radial arcs than standard CDM but they are on average more magnified. We also appreciate differences in the arc morphology that could be used to statistically favour one model versus the other., GY acknowledges financial support by the MINECO/FEDER under project grant AYA2015-63810-P and MICIU/FEDER under project grant PGC2018-094975-C21. WC acknowledges the supported by the European Research Council under grant number 670193. MM acknowledges support from PRIN-MIUR ‘Cosmology and Fundamental Physics: illuminating the dark universe with Euclid’, and from ASI through contract Euclid Phase D 1.05.04.37.01.

The James Webb Space Telescope (JWST) is expected to observe galaxies at $z>10$ that are presently inaccessible. Here, we use a self-consistent empirical model, the UniverseMachine, to generate mock galaxy catalogues and lightcones over the redshift range $z=0-15$. These data include realistic galaxy properties (stellar masses, star formation rates, and UV luminosities), galaxy-halo relationships, and galaxy-galaxy clustering. Mock observables are also provided for different model parameters spanning observational uncertainties at $z10^7 M_\odot$ and/or $M_{1500}12$ expand dramatically, so efforts to detect $z>12$ galaxies will provide the most valuable constraints on galaxy formation models. The faint-end slopes of the stellar mass/luminosity functions at a given mass/luminosity threshold steepen as redshift increases. This is because observable galaxies are hosted by haloes in the exponentially falling regime of the halo mass function at high redshifts. Hence, these faint-end slopes are robustly predicted to become shallower below current observable limits ($M_\ast < 10^7M_\odot$ or $M_\mathrm{1500}>-17$). For reionization models, extrapolating luminosity functions with a constant faint-end slope from $M_{1500}=-17$ down to $M_{1500}=-12$ gives the most reasonable upper limit for the total UV luminosity and cosmic star formation rate up to $z\sim 12$. We compare to three other empirical models and one semi-analytic model, showing that the range of predicted observables from our approach encompasses predictions from other techniques. Public catalogues and lightcones for common fields are available online., Comment: 17 pages, MNRAS submitted. Catalogs and lightcones available at https://www.peterbehroozi.com/data.html

Galaxy clusters can play a key role in modern cosmology provided their evolution is properly understood. However, observed clusters give us only a single timeframe of their dynamical state. Therefore, finding present observable data of clusters that are well correlated to their assembly history constitutes an inestimable tool for cosmology. Former studies correlating environmental descriptors of clusters to their formation history are dominated by halo mass - environment relations. This paper presents a mass-free correlation between the present neighbor distribution of cluster-size halos and the latter mass assembly history. From the Big Multidark simulation, we extract two large samples of random halos with masses ranging from Virgo to Coma cluster sizes. Additionally, to find the main environmental culprit for the formation history of the Virgo cluster, we compare the Virgo-size halos to 200 Virgo-like halos extracted from simulations that resemble the local Universe. The number of neighbors at different cluster-centric distances permits discriminating between clusters with different mass accretion histories. Similarly to Virgo-like halos, clusters with numerous neighbors within a distance of about 2 times their virial radius experience a transition at z~1 between an active period of mass accretion, relative to the mean, and a quiet history. On the contrary, clusters with few neighbors share an opposite trend: from passive to active assembly histories. Additionally, clusters with massive companions within about 4 times their virial radius tend to have recent active merging histories. Therefore, the radial distribution of cluster neighbors provides invaluable insights into the past history of these objects., Accepted for publication in MNRAS, 10 pages, 5 figures, 4 tables

Using the catalogues of galaxy clusters from The Three Hundred project, modelled with both hydrodynamic simulations, (Gadget-X and Gadget-MUSIC), and semi-analytic models (SAMs), we study the scatter and self-similarity of the profiles and distributions of the baryonic components of the clusters: the stellar and gas mass, metallicity, the stellar age, gas temperature, and the (specific) star formation rate. Through comparisons with observational results, we find that the shape and the scatter of the gas density profiles matches well the observed trends including the reduced scatter at large radii which is a signature of self-similarity suggested in previous studies. One of our simulated sets, Gadget-X, reproduces well the shape of the observed temperature profile, while Gadget-MUSIC has a higher and flatter profile in the cluster centre and a lower and steeper profile at large radii. The gas metallicity profiles from both simulation sets, despite following the observed trend, have a relatively lower normalisation. The cumulative stellar density profiles from SAMs are in better agreement with the observed result than both hydrodynamic simulations which show relatively higher profiles. The scatter in these physical profiles, especially in the cluster centre region, shows a dependence on the cluster dynamical state and on the cool-core/non-cool-core dichotomy. The stellar age, metallicity and (s)SFR show very large scatter, which are then presented in 2D maps. We also do not find any clear radial dependence of these properties. However, the brightest central galaxies have distinguishable features compared to the properties of the satellite galaxies., Comment: 20 pages, 18 figures

We perform a simulation with Galacticus, a semi-analytical galaxy formation model, to predict the number counts of H$\alpha$ and [OIII] emitting galaxies. With a state-of-the-art N-body simulation, UNIT, we first calibrate Galacticus with the current observation of H$\alpha$ luminosity function. The resulting model coupled with a dust attenuation model, can reproduce the current observations, including the H$\alpha$ luminosity function from HiZELS and number density from WISP. We extrapolate the model prediction to higher redshift and the result is found to be consistent with previous investigations. We then use the same galaxy formation model to predict the number counts for [OIII] emitting galaxies. The result provides further validation of our galaxy formation model and dust model. We present number counts of H$\alpha$ and [OIII] emission line galaxies for three different line flux limits: $5\times10^{-17}$erg/s/cm$^{2}$, $1\times10^{-16}$erg/s/cm$^{2}$ (6.5$\sigma$ nominal depth for WFIRST GRS), and $2\times10^{-16}$erg/s/cm$^{2}$ (3.5$\sigma$ depth of Euclid GRS). At redshift $2, Comment: 12 pages, 7+2 figures, 2 table; updated to match version accepted by MNRAS

We report the non-thermal pressure fraction (Pnt/Ptot) obtained from a three-dimensional triaxial analysis of 16 galaxy clusters in the CLASH sample using gravitational lensing (GL) data primarily from Subaru and HST, X-ray spectroscopic imaging from Chandra, and Sunyaev–Zel’dovich effect (SZE) data from Planck and Bolocam. Our results span the approximate radial range 0.015–0.4 R200m (∼35–1000 kpc). At cluster-centric radii smaller than 0.1 R200m the ensemble average Pnt/Ptot is consistent with zero with an upper limit of 9 per cent, indicating that heating from active galactic nuclei and other relevant processes does not produce significant deviations from hydrostatic equilibrium (HSE). The ensemble average Pnt/Ptot increases outside of this radius to approximately 20 per cent at 0.4 R200m, as expected from simulations, due to newly accreted material thermalizing via a series of shocks. Also in agreement with simulations, we find significant cluster-to-cluster variation in Pnt/Ptot and little difference in the ensemble average Pnt/Ptot based on dynamical state. We conclude that on average, even for diverse samples, HSE-derived masses in the very central regions of galaxy clusters require only modest corrections due to non-thermal motions.

We introduce the Lomonosov suite of high-resolution N-body cosmological simulations covering a full box of size 32 $h^{-1}$ Mpc with low-mass resolution particles ($2\times10^7$ $h^{-1}\,M_\odot$) and three zoom-in simulations of overdense, underdense and mean density regions at much higher particle resolution ($4\times10^4$ $h^{-1}\,M_\odot$). The main purpose of this simulation suite is to extend the concentration-mass relation of dark matter halos down to masses below those typically available in large cosmological simulations. The three different density regions available at higher resolution provide a better understanding of the effect of the local environment on halo concentration, known to be potentially important for small simulation boxes and small halo masses. Yet, we find the correction to be small in comparison with the scatter of halo concentrations. We conclude that zoom simulations, despite their limited representativity of the volume of the Universe, can be effectively used for the measurement of halo concentrations at least at the halo masses probed by our simulations. In any case, after a precise characterization of this effect, we develop a robust technique to extrapolate the concentration values found in zoom simulations to larger volumes with greater accuracy. All together, Lomonosov provides a measure of the concentration-mass relation in the halo mass range $10^7-10^{10}$ $h^{-1}\,M_\odot$ with superb halo statistics. This work represents a first important step to measure halo concentrations at intermediate, yet vastly unexplored halo mass scales, down to the smallest ones. All Lomonosov data and files are public for community's use., Comment: 11 pages, accepted by MNRAS

Reliable extraction of cosmological information from clustering measurements of galaxy surveys requires estimation of the error covariance matrices of observables. The accuracy of covariance matrices is limited by our ability to generate sufficiently large number of independent mock catalogs that can describe the physics of galaxy clustering across a wide range of scales. Furthermore, galaxy mock catalogs are required to study systematics in galaxy surveys and to test analysis tools. In this investigation, we present a fast and accurate approach for generation of mock catalogs for the upcoming galaxy surveys. Our method relies on low-resolution approximate gravity solvers to simulate the large scale dark matter field, which we then populate with halos according to a flexible nonlinear and stochastic bias model. In particular, we extend the \textsc{patchy} code with an efficient particle mesh algorithm to simulate the dark matter field (the \textsc{FastPM} code), and with a robust MCMC method relying on the \textsc{emcee} code for constraining the parameters of the bias model. Using the halos in the BigMultiDark high-resolution $N$-body simulation as a reference catalog, we demonstrate that our technique can model the bivariate probability distribution function (counts-in-cells), power spectrum, and bispectrum of halos in the reference catalog. Specifically, we show that the new ingredients permit us to reach percentage accuracy in the power spectrum up to $k\sim 0.4\; \,h\,{\rm Mpc}^{-1}$ (within 5\% up to $k\sim 0.6\; \,h\,{\rm Mpc}^{-1}$) with accurate bispectra improving previous results based on Lagrangian perturbation theory., 12 pages, 4 figures, published version

N-body cosmological simulations are an essential tool to understand the observed distribution of galaxies. We use the MultiDark simulation suite, run with the Planck cosmological parameters, to revisit the mass and velocity functions. At redshift z = 0, the simulations cover four orders of magnitude in halo mass from ~1011M¿ with 8783 874 distinct haloes and 532 533 subhaloes. The total volume used is ~515 Gpc3, more than eight times larger than in previous studies. We measure and model the halo mass function, its covariance matrix w.r.t halo mass and the large-scale halo bias.With the formalism of the excursion-set mass function, we explicit the tight interconnection between the covariance matrix, bias and halo mass function. We obtain a very accurate ( < 2 per cent level) model of the distinct halo mass function. We also model the subhalo mass function and its relation to the distinct halo mass function. The set of models obtained provides a complete and precise framework for the description of haloes in the concordance Planck cosmology. Finally, we provide precise analytical fits of the Vmax maximum velocity function up to redshift z < 2.3 to push for the development of halo occupation distribution using Vmax. The data and the analysis code are made publicly available in the Skies and Universes data base. © 2017 The Authors., JC and FP acknowledge support from the Spanish MICINNs Consolider-Ingenio 2010 Programme under grant MD CSD2009-00064, MINECO Centro de Excelencia Severo Ochoa Programme under the grants SEV-20120249, FPA2012-34694, and the projects AYA2014-60641-C2-1-P and AYA2012-31101. GY acknowledges financial support from MINECO/FEDER (Spain) under project number AYA2012-31101 and AYA2015-63810-P.

Following Cui et al. 2018 (hereafter Paper I) on the classification of large-scale environments (LSE) at z = 0, we push our analysis to higher redshifts and study the evolution of LSE and the baryon distributions in them. Our aim is to investigate how baryons affect the LSE as a function of redshift. In agreement with Paper I, the baryon models have negligible effect on the LSE over all investigated redshifts. We further validate the conclusion obtained in Paper I that the gas web is an unbiased tracer of total matter -- even better at high redshifts. By separating the gas mainly by temperature, we find that about 40 per cent of gas is in the so-called warm-hot intergalactic medium (WHIM). This fraction of gas mass in the WHIM decreases with redshift, especially from z = 1 (29 per cent) to z = 2.1 (10 per cent). By separating the whole WHIM gas mass into the four large-scale environments (i.e. voids, sheets, filaments, and knots), we find that about half of the WHIM gas is located in filaments. Although the total gas mass in WHIM decreases with redshift, the WHIM mass fractions in the different LSE seem unchanged.

Accurately predicting structural properties of dark matter halos is one of the fundamental goals of modern cosmology. We use the new suite of MultiDark cosmological simulations to study the evolution of dark matter halo density profiles, concentrations, and velocity anisotropies. The MultiDark simulations cover a large range of masses 1e10-1e15Msun and volumes upto 50Gpc**3. The total number of dark matter halos in all the simulations exceeds 60 billion. We find that in order to understand the structure of dark matter halos and to make ~1% accurate predictions for density profiles, one needs to realize that halo concentration is more complex than the traditional ratio of the virial radius to the core radius in the NFW profile. For massive halos the averge density profile is far from the NFW shape and the concentration is defined by both the core radius and the shape parameter alpha in the Einasto approximation. Combining results from different redshifts, masses and cosmologies, we show that halos progress through three stages of evolution. (1) They start as rare density peaks that experience very fast and nearly radial infall. This radial infall brings mass closer to the center producing a high concentrated halo. Here, the halo concentration increases with the increasing halo mass and the concentration is defined by the alpha parameter with nearly constant core radius. Later halos slide into (2) the plateau regime where the accretion becomes less radial, but frequent mergers still affect even the central region. Now the concentration does not depend on halo mass. (3) Once the rate of accretion slows down, halos move into the domain of declining concentration-mass relation because new accretion piles up mass close to the virial radius while the core radius is staying constant. We provide accurate analytical fits to the numerical results for halo density profiles and concentrations., 22 pages, 23 figures, accepted by MNRAS

Recent numerical studies of the dark matter density profiles of massive galaxy clusters ($M_{\rm halo} > 10^{15}$M$_{\odot}$) show that their median radial mass density profile remains unchanged up to $z > 1$, displaying a highly self-similar evolution. We verify this by using the data set of the THE THREE HUNDRED project, i.e. 324 cluster-sized haloes as found in full physics hydrodynamical simulations. We track the progenitors of the mass-complete sample of clusters at $z=0$, and find that their median shape is already in place by $z=2.5$. However, selecting a dynamically relaxed subsample ($\sim16$ per cent of the clusters), we observe a shift of the scale radius $r_s$ towards larger values at earlier times. Classifying the whole sample by formation time, this evolution is understood as a result of a two-phase halo mass accretion process. Early-forming clusters -- identified as relaxed today -- have already entered their slow accretion phase, hence their mass growth occurs mostly at the outskirts. Late-forming clusters -- which are still unrelaxed today -- are in their fast accretion phase, thus the central region of the clusters is still growing. We conclude that the density profile of galaxy clusters shows a profound self-similarity out to redshifts $z\sim2.5$. This result holds for both gas and total density profiles when including baryonic physics, as reported here for two rather distinct sub-grid models., accepted for publication in MNRAS; 15 pages, 10 figures, 3 tables

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2015 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, We study the relation between halo mass and its environment from a probabilistic perspective. We find that halo mass depends not only on local dark matter density, but also on non-local quantities such as the cosmic web environment and the halo- exclusion effect. Given these accurate relations, we have developed the HADRON-code (Halo mAss Distribution ReconstructiON), a technique which permits us to assign halo masses to a distribution of haloes in three-dimensional space. This can be ap- plied to the fast production of mock galaxy catalogues, by assigning halo masses, and reproducing accurately the bias for diferent mass cuts. The resulting clustering of the halo populations agree well with that drawn from the BigMultiDark N-body simulation: the power spectra are within 1-_ up to scales of k = 0:2 hMpc-1, when using augmented Lagrangian perturbation theory based mock catalogues. Only the most massive haloes show a larger deviation. For these, we finnd evidence of the halo- exclusion effect. A clear improvement is achieved when assigning the highest masses to haloes with a minimum distance separation. We also compute the 2- and 3-point correlation functions, and find an excellent agreement with N-body results. Our work represents a quantitative application of the cosmic web classification. It can have fur- ther interesting applications in the multi-tracer analysis of the large-scale structure for future galaxy surveys, CZ and CT acknowledge support from Tsinghua University, and 973 program No. 2013CB834906. FP were supported by the Spanish MICINNs Consolider-Ingenio 2010 Programme under grant MultiDark CSD2009-00064 and AYA2010-21231-C02-01 grant, the Comunidad de Madrid under grant HEPHACOS S2009/ESP-1473, and Spanish MINECOs Centro de Excelencia Severo Ochoa Programme under grant SEV-2012-0249. GY acknowledges support from the Spanish MINECO under research grants AYA2012-31101, FPA2012-34694 and Consolider Ingenio SyeC CSD2007-0050

This article has been accepted for publication in Monthly Notices of the Royal Astronomical Society ©: 2015 The Authors. Published by Oxford University Press on behalf of the Royal Astronomical Society. All rights reserved, We examine the integrated Sachs-Wolfe (ISW) imprint of voids and superclusters on the cosmic microwave background. We first study results from the Jubilee N-body simulation. From Jubilee, we obtain the full-sky ISW signal from structures out to redshift z=1.4 and a mock luminous red galaxy (LRG) catalogue. We confirm that the expected signal in the concordance \Lambda CDM model is very small and likely to always be much smaller than the anisotropies arising at the last scattering surface. Any current detections of such an imprint must, therefore, predominantly arise from something other than an ISW effect in a \Lambda CDM universe. Using the simulation as a guide, we then look for the signal using a catalogue of voids and superclusters from the Sloan Digital Sky Survey. We find a result that is consistent with the \Lambda CDM model, i.e. a signal consistent with zero., SH and ITI were supported by the Science and Technologies Facilities Council (grant numberST/I000976/1). The research leading to these results has received funding from the European Research Council under the European Unions Seventh Framework Programme (FP/20072013) / ERC Grant Agreement No. [308082]. SN acknowledges support from Academy of Finland grant 1263714. AK is supported by the Ministerio de Economía y Competividad (MINECO) in Spain through grant AYA2012-31101 as well as the Consolider-Ingenio 2010 Programme of the Spanish Ministerio de Ciencia e Innovación (MICINN) under grant MultiDark CSD2009-00064. He also acknowledges support from the Australian Research Council (ARC) grants DP130100117 and DP140100198. GY acknowledges support from MINECO (Spain) under research grants AYA2012-31101 and FPA2012-34694

We present the results of a series of adiabatic hydrodynamical simulations of several quintessence models (both with a free and an interacting scalar field) in comparison to a standard \LCDM\ cosmology. For each we use $2\times1024^3$ particles in a $250$\hMpc\ periodic box assuming WMAP7 cosmology. In this work we focus on the properties of haloes in the cosmic web at $z=0$. The web is classified into \emph{voids}, \emph{sheets}, \emph{filaments} and \emph{knots} depending on the eigenvalues of the velocity shear tensor, which are an excellent proxy for the underlying overdensity distribution. We find that the properties of objects classified according to their surrounding environment shows a substantial dependence on the underlying cosmology; for example, while $V_{\rm max}$ shows average deviations of $\approx5$ per cent across the different models when considering the full halo sample, comparing objects classified according to their environment, the size of the deviation can be as large as $20$ per cent. We also find that halo spin parameters are positively correlated to the coupling, whereas halo concentrations show the opposite behaviour. Furthermore, when studying the concentration-mass relation in different environments, we find that in all cosmologies underdense regions have a larger normalization and a shallower slope. While this behaviour is found to characterize all the models, differences in the best-fit relations are enhanced in (coupled) dark energy, thus providing a clearer prediction for this class of models., 16 pages, 10 figures. Accepted for publication in MNRAS

We present the first results from the Juropa Hubble Volume project, based a large N-body, dark matter-only cosmological simulation with a volume of V = (6 h-1 Gpc)3, containing 60003 particles, performed within the concordance cold dark matter cosmological model. The simulation volume is sufficient to probe extremely large length-scales in the universe, whilst at the same time the particle count is high enough so that dark matter haloes down to 1.5 × 1012 h-1M⊙ can be resolved. At z = 0, we identify over 400 million haloes. The cluster mass function is derived using three different halofinders and compared to fitting functions in the literature. The distribution of clusters of maximal mass across redshifts agrees well with predicted masses of extreme objects, and we explicitly confirm that the Poisson distribution is very good at describing the distribution of rare clusters. The Poisson distribution also matches well the level to which cosmic variance can be expected to affect number counts of high-mass clusters. We find that objects like the Bullet cluster exist in the far-tail of the distribution of mergers in terms of relative collisional speed. We also derive the number counts of voids in the simulation box for z = 0, 0.5 and 1. © 2013 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society., ITI was supported by the SEPNet and the Science and Technology Facilities Council grants ST/F002858/1 and ST/I000976/1. AK was supported by the Spanish Ministerio de Ciencia e Innovación (MICINN) in Spain through the Ramón y Cajal programme as well as the grants AYA 2009-13875-C03-02, AYA2009-12792-C03-03, CSD2009-00064, CAM S2009/ESP-1496 and the Ministerio de Economía y Competitividad (MINECO) through grantAYA2012-31101.GYacknowledges support fromMINECO(Spain) under research grantsAYA2009-13875-C03-02, FPA2009-08958, AYA2012-31101 and Consolider Ingenio SyeC CSD2007-0050. He also thanks Comunidad de Madrid for support under PRICIT project ASTROMADRID (S2009/ESP-146). JMD and EMG acknowledge support of the consolider project CSD2010-00064 funded by the Ministerio de Economía y Competitividad and AYA2012-39475-C02-01.

A new approach for the classification of the cosmic web is presented. In extension of the previous work of Hahn et al. (2007) and Forero-Romero et al. (2009) the new algorithm is based on the analysis of the velocity shear tensor rather than the gravitational tidal tensor. The procedure consists of the construction of the the shear tensor at each (grid) point in space and the evaluation of its three eigenvectors. A given point is classified to be either a void, sheet, filament or a knot according to the number of eigenvalues above a certain threshold, 0, 1, 2, or 3 respectively. The threshold is treated as a free parameter that defines the web. The algorithm has been applied to a dark matter only, high resolution simulation of a box of side-length 64$h^{-1}$Mpc and N = $1024^3$ particles with the framework of the WMAP5/LCDM model. The resulting velocity based cosmic web resolves structures down to

We present the results of a series of cosmological $N$-body simulations of a Vector Dark Energy (VDE) model, performed using a suitably modified version of the publicly available \texttt{GADGET}-2 code. The setups of our simulations were calibrated pursuing a twofold aim: 1) to analyze the large scale distribution of massive objects and 2) to determine the properties of halo structure in this different ramework.We observe that structure formation is enhanced in VDE, since the mass function at high redshift is boosted up to a factor of ten with respect to \LCDM, possibly alleviating tensions with the observations of massive clusters at high redshifts and early reionization epoch. Significant differences can also be found for the value of the growth factor, that in VDE shows a completely different behaviour, and in the distribution of voids, which in this cosmology are on average smaller and less abundant. We further studied the structure of dark matter haloes more massive than $5\times10^{13}$\hMsun, finding that no substantial difference emerges when comparing spin parameter, shape, triaxiality and profiles of structures evolved under different cosmological pictures. Nevertheless, minor differences can be found in the concentration-mass relation and the two point correlation function; both showing different amplitudes and steeper slopes.Using an additional series of simulations of a \LCDM\ scenario with the same $\Omega_M$ and $\sigma_8$ used in the VDE cosmology, we have been able to establish whether the modifications induced in the new cosmological picture were due to the particular nature of the dynamical dark energy or a straightforward consequence of the cosmological parameters.

We use the recently proposed scale-free mass estimators to determine the masses of the Milky Way (MW) and Andromeda (M31) galaxy in a dark matter only Constrained Local Universe Simulation. While these mass estimators work rather well for isolated spherical host systems, we examine here their applicability to a simulated binary system with a unique satellite population similar to the observed satellites of MW and M31. We confirm that the scale-free estimators work also very well in our simulated Local Group galaxies with the right number of satellites which follow the observed radial distribution. In the isotropic case and under the assumption that the satellites are tracking the total gravitating mass, the power-law index of the radial satellite distribution N(< r) ∝r3−γ is directly related to the host’s mass profile M(< r) ∝r1−α as α=γ− 2. The use of this relation for any given γ leads to highly accurate mass estimations which is a crucial point for observer, since they do not know a priori the mass profile of the MW and M31 haloes. We discuss possible bias in the mass estimators and conclude that the scale-free mass estimators can be satisfactorily applied to the real MW and M31 system.

Upcoming radio telescopes will allow to study the radio sky at low frequencies with unprecedent sensitivity and resolution. New surveys are expected to discover a large number of new radio sources. Here we investigate the abundance of radio relics, i.e. steep-spectrum diffuse radio emission coming from the periphery of galaxy clusters, which are believed to trace shock waves induced by cluster mergers. With the advent of comprehensive relic samples a framework is needed to analyze statistically the relic abundance. To this end, we introduce the probability to find a relic located in a galaxy cluster with given mass and redshift allowing us to relate the halo mass function of the Universe with the radio relic number counts. Up to date about 45 relics have been reported and we compile the resulting counts, N(>S_1.4). In principle, the parameters of the distribution could be determined using a sufficiently large relic sample. However, since the number of known relics is still small for that purpose we use the MareNostrum Universe simulation to determine the relic radio power scaling with cluster mass and redshift. Our model is able to reproduce the recently found tentative evidence for an increase in the fraction of clusters hosting relics, both with X-ray luminosity and redshift, using an X-ray flux limited cluster sample. Moreover, we find that a considerable fraction of faint relics (S_1.4 < ~10 mJy) reside in clusters with an X-ray flux below ~3e-12 erg/s/cm^2. Finally, we estimate the number of radio relics which await discovery by future low frequency surveys proposed for LOFAR and WSRT. We estimate that the WODAN survey proposed for WSRT may discover 900 relics and that the LOFAR-Tier 1-120 MHz survey may discover about 2500 relics. However, the actual number of newly discovered relics will crucially depend on the existence of sufficiently complete galaxy cluster catalogues.

We present theoretical results for the expected fraction of Lyman Break Galaxies (LBGs) to be detected as strong Lyman-α emitters (LAEs) in the redshift range 5 � z � 7. We base our analysis on the 2-billion particle SPH simulation MareNostrum High-z Universe. We approximate galaxies as static dusty slabs with an additional clumpy dust distribution affecting stellar populations younger than 25 Myr. The model for the Lyman-α escape fraction is based on the results of our Monte-Carlo radiative transfer code (CLARA) for a slab configuration. We also fix the transmission of Lyman-α photons through the intergalactic medium to a constant value of 50% at all redshifts. From the results of this model we calculate xLy�, the fraction of Lyman Break Galaxies with Lyα equivalent width (EW) larger than 50u We find a remarkable agreement with observational data at 4.5 < z < 6. For bright ( 22 < MUV < 20.5) and faint ( 20.5 < MUV < 18.5) galaxies our model predicts xLy� = 0.02 ± 0.01 and xLy� = 0.47±0.01 while observers report xLy� = 0.08±0.02 and xLy� = 0.47±0.16, respectively. Additional evolution of the extinction model at redshift z � 7, that decreases the intensity of transmitted Lyman-α radiation by a factor of fT = 0.4 as to match the LAE luminosity function at z � 6.5, naturally provides a good match for

The outer haloes of the Milky Way and Andromeda galaxies contain as much important information on their assembly and formation history as the properties of the discs resident in their centres. In this paper we have used the Constrained Local UniversE Simulation project to disentangle the stellar and DM component of three galaxies that resemble the MW, M31 and M33 using both DM only and DM + gas-dynamical simulations. Stars that are accreted in substructures and then stripped follow a completely different radial distribution than the stripped DM: the stellar halo is much more centrally concentrated than DM. In order to understand how the same physical process can lead to different z=0 radial profiles, we examined the potential at accretion of each stripped particle. We found that star particles sit at systematically higher potentials than DM, making them harder to strip. We then searched for a threshold in the potential of accreted particles \phi_[th], above which DM particles behave as star particles. We found such a threshold at >16 \phi_{subhalo}, where \phi_{subhalo} is the potential at a subhaloes edge at the time of accretion. Thus a rule as simple as selecting particles according to their potential at accretion is able to reproduce the effect that the complicated physics of star formation has on the stellar distribution. This niversal result reproduces the stellar halo to an accuracy of within ~2%. Studies which make use of DM particles as a proxy for stars will undoubtedly miscalculate their proper radial distribution and structure unless particles are selected according to their potential at accretion. Furthermore, we have examined the time it takes to strip a given star or DM particle after accretion. We find that, owing to their higher binding energies, stars take longer to be stripped than DM. The stripped DM halo is thus considerably older than the stripped stellar halo.

We make detailed theoretical predictions for the assembly properties of the Local Group (LG) in the standard Λ cold dark matter cosmological model. We use three cosmological N-body dark matter simulations from the Constrained Local Universe Simulations project, which are designed to reproduce the main dynamical features of the matter distribution down to the scale of a few Mpc around the LG. Additionally, we use the results of an unconstrained simulation with a 60 times larger volume to calibrate the influence of cosmic variance. We characterize the mass aggregation history (MAH) for each halo by three characteristic times: the formation, assembly and last major merger times. A major merger is defined by a minimal mass ratio of 10: 1. We find that the three LGs share a similar MAH with formation and last major merger epochs placed on average ≈10–12 Gyr ago. Between 12 and 17 per cent of the haloes in the mass range 5 × 1011 < Mh < 5 × 1012 h−1 M⊙ have a similar MAH. In a set of pairs of haloes within the same mass range, a fraction of 1–3 per cent share similar formation properties as both haloes in the simulated LG. An unsolved question posed by our results is the dynamical origin of the MAH of the LGs. The isolation criteria commonly used to define LG-like haloes in unconstrained simulations do not narrow down the halo population into a set with quiet MAHs, nor does a further constraint to reside in a low-density environment. The quiet MAH of the LGs provides a favourable environment for the formation of disc galaxies like the Milky Way and M31. The timing for the beginning of the last major merger in the Milky Way dark matter halo matches with the gas-rich merger origin for the thick component in the galactic disc. Our results support the view that the specific large- and mid-scale environments around the LG play a critical role in shaping its MAH and hence its baryonic structure at present.

Using CLARA (Code for Lyman Alpha Radiation Analysis) we constrain the escape fraction of Lyman-Alpha radiation in galaxies in the redshift range 5 ~5. The LAE Luminosity Function (LF) based on the extended clumpy model reproduces broadly the bright end of the LF derived from observations at z \sim 5 and z \sim 6. At z \sim 7 our model over-predicts the LF by roughly a factor of four, presumably because the effects of the neutral intergalactic medium are not taken into account. The remaining tension between the observed and simulated faint end of the LF, both in the UV-continuum and Lyman-Alpha at redshifts z \sim 5 and z \sim 6 points towards an overabundance of simulated LAEs hosted in haloes of masses 1.0x10^10h-1Msol < Mh < 4.0x10^10h-1Msol. Given the difficulties in explaining the observed overabundance by dust absorption, a probable origin of the mismatch are the high star formation rates in the simulated haloes around the quoted mass range. A more efficient supernova feedback should be able to regulate the star formation process in the shallow potential wells of these haloes.

We analyse the proximity zone of the intergalactic matter around high-redshift quasars in a cosmological environment. In a box of 64 h-1 Mpc base length we employ dark matter only simulations. For estimating the hydrogen temperature and density distribution we use the effective equation of state. Hydrogen is assumed to be in photoionisation equilibrium with a model background flux which is fit to recent observations of the mean optical depth and transmission flux statistics. At redshifts z = 3, 4, and 4.8, we select model quasar positions at the centre of the 20 most massive halos and 100 less massive halos identified in the simulation. From each assumed quasar position we cast 100 random lines of sight for two box length including the changes in the ionisation fractions by the QSO flux field and derive mock Ly{\alpha} spectra. The proximity effect describes the dependence of the mean normalised optical depth {\xi} = {\tau}eff, QSO/{\tau}eff, Ly{\alpha} as a function of the ratio of the ionisation rate by the QSO and the background field, {\omega} = {\Gamma}QSO/{\Gamma}UVB, i.e. the profile {\xi} = (1 + {\omega}/a)-0.5, where a strength parameter a is introduced. The strength parameter measures the deviation from the theoretical background model and is used to quantify any influence of the environmental density field. We reproduce an unbiased measurement of the proximity effect which is not affected by the host halo mass. The scatter between different lines of sight and different quasar host positions increases with decreasing redshift. Around the host halos, we find only a slight average overdensity in the proximity zone at comoving radii of 1 < rc < 10h-1 Mpc. However, a clear power-law correlation of the strength parameter with the average overdensity in rc is found, showing an overestimation of the ionising background in overdense regions and an underestimation in underdense regions.

[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).

We present the first detailed structure formation and radiative transfer simulations of the reionization history of our cosmic neighbourhood. To this end, we follow the formation of the Local Group of galaxies and nearby clusters by means of constrained simulations, which use the available observational constraints to construct a representation of those structures which reproduces their actual positions and properties at the present time. We find that the reionization history of the Local Group is strongly dependent on the assumed photon production efficiencies of the ionizing sources, which are still poorly constrained. If sources are relatively efficient, i.e. the process is 'photon-rich', the Local Group is primarily ionized externally by the nearby clusters. Alternatively, if the sources are inefficient, i.e. reionization is 'photon-poor' the Local Group evolves largely isolated and reionizes itself. The mode of reionization, external vs. internal, has important implications for the evolution of our neighbourhood, in terms of e.g. its satellite galaxy populations and primordial stellar populations. This therefore provides an important avenue for understanding the young universe by detailed studies of our nearby structures.

We study the differences and similarities in the luminosities of bound, infalling and the so-called backsplash (Gill et al. 2005) galaxies of the Milky Way and M31 using a hydrodynamical simulation performed within the Constrained Local UniversE Simulation (CLUES) project. The simulation models the formation of the Local Group within a self-consistent cosmological framework. We find that even though backsplash galaxies passed through the virial radius of their host halo and hence may have lost a (significant) fraction of their mass, their stellar populations are hardly affected. This leaves us with comparable luminosity functions for infalling and backsplash galaxies and hence little hope to decipher their past (and different) formation and evolutionary histories by luminosity measurements alone. Nevertheless, due to the tidal stripping of dark matter we find that the mass-to-light ratios have changed when comparing the various populations against each other: they are highest for the infalling galaxies and lowest for the bound satellites with the backsplash galaxies in-between.

Using a high resolution DM simulation of the Local Group (LG), conducted within the framework of the Constrained Local UniversE Simulation (CLUES) project, we investigate the nature of how satellites of the MW and M31 are accreted. Satellites of these galaxies are accreted anisotropically, entering the virial radius of their hosts from specific "spots" with respect to the large scale structure. Furthermore, material which is tidally stripped from these satellites is also, at z=0, distributed anisotropically and is characterized by an ellipsoidal sub-volume embedded in the halo. The angular pattern created by the locus of satellite infall points and the projected stripped DM is investigated within a coordinate system determined by the location of the LG companion and the simulated Virgo cluster across concentric shells ranging from 0.1 to 5 r_vir. Remarkably, the principal axis of the ellipsoidal sub-volume shows a coherent alignment extending from well within the halo to a few r_vir. A spherical harmonics transform applied to the angular distributions confirms the visual impression: namely, the angular distributions of both the satellites entry points and stripped DM is dominated by the l=2 quadrupole term, whose major principal axis is aligned across the shells considered. It follows that the structure of the outer (r>0.5 r_vir) regions of the main halos is closely related to the cosmic web within which it is embedded. Given the hypothesis that a large fraction of the stellar halo of the Milky Way has been accreted from satellites, the present results can be applied to the stellar halo of the MW and M31. We predict that the remnants of tidally stripped satellites should be embedded in streams of material composed of dark matter and stars. The present results can therefore shed light on the existence of satellites embedded within larger streams of matter, such as the Segue 2 satellite.

We use dark matter haloes identified in the MareNostrum Universe and galaxy groups identified in the Sloan Data Release 7 galaxy catalogue, to study the relation between halo shape and halo dynamics, parametrizing out the mass of the systems. A strong shape-dynamics, independent of mass, correlation is present in the simulation data, which we find it to be due to different halo formation times. Early formation time haloes are, at the present epoch, more spherical and have higher velocity dispersions than late forming-time haloes. The halo shape-dynamics correlation, albeit weaker, survives the projection in 2D (ie., among projected shape and 1-D velocity dispersion). A similar shape-dynamics correlation, independent of mass, is also found in the SDSS DR7 groups of galaxies and in order to investigate its cause we have tested and used, as a proxy of the group formation time, a concentration parameter. We have found, as in the case of the simulated haloes, that less concentrated groups, corresponding to late formation times, have lower velocity dispersions and higher elongations than groups with higher values of concentration, corresponding to early formation times.

We use a simulation performed within the Constrained Local UniversE Simulation (CLUES) project to study a realistic Local Group-like object. We employ this group as a numerical laboratory for studying the evolution of the population of its subhaloes from the point of view of the effects it may have on the origin of different types of dwarf galaxies. We focus on the processes of tidal stripping of the satellites, their interaction, merging and grouping before infall. The tidal stripping manifests itself in the transition between the phase of mass accretion and mass loss seen in most subhaloes, which occurs at the moment of infall on to the host halo, and the change of the shape of their mass function with redshift. Although the satellites often form groups, they are loosely bound within them and do not interact with each other. The infall of a large group could however explain the observed peculiar distribution of the Local Group satellites, but only if it occurred recently. Mergers between prospective subhaloes are significant only during an early stage of evolution, i.e. more than 7 Gyr ago, when they are still outside the host haloes. Such events could thus contribute to the formation of more distant early type Milky Way companions. Once the subhaloes enter the host halo the mergers become very rare., Comment: 14 pages, 16 figures, published in MNRAS, typos corrected to match published version

[Abridged] We study the thermal and kinetic Sunyaev-Zel'dovich (SZ) effect associated with superclusters of galaxies using the MareNostrum Universe SPH simulation. We consider superclusters similar to the Corona Borealis Supercluster (CrB-SC). This paper is motivated by the detection at 33GHz of a strong temperature decrement in the CMB towards the core of this supercluster. Multifrequency observations with VSA and MITO suggest the existence of a thermal SZ effect component in the spectrum of this cold spot, which would account for roughly 25% of the total observed decrement. We identify nine regions containing superclusters similar to CrB-SC, obtain the associated SZ maps and calculate the probability of finding such SZ signals arising from hot gas within the supercluster. Our results show that WHIM produces a thermal SZ effect much smaller than the observed value. Neither can summing the contribution of small clusters and galaxy groups in the region explain the amplitude of the SZ signal. When we take into account the actual posterior distribution from the observations, the probability that WHIM can cause a thermal SZ signal like the one observed is

Using high-resolution simulations within the Cold and Warm Dark Matter models we study the evolution of small scale structure in the Local Volume, a sphere of 8 Mpc radius around the Local Group. We compare the observed spectrum of mini-voids in the Local Volume with the spectrum of mini-voids determined from the simulations. We show that the \LWDM model can easily explain both the observed spectrum of mini-voids and the presence of low-mass galaxies observed in the Local Volume, provided that all haloes with circular velocities greater than 20 km/s host galaxies. On the contrary within the LCDM model the distribution of the simulated mini-voids reflects the observed one if haloes with maximal circular velocities larger than 35 km/s host galaxies. This assumption is in contradiction with observations of galaxies with circular velocities as low as 20 km/s in our Local Universe. A potential problem of the LWDM model could be the late formation of the haloes in which the gas can be efficiently photo-evaporated. Thus star formation is suppressed and low-mass haloes might not host any galaxy at all., 13 pages, 10 figures, version 2, subsection 3.1 added, accepted to MNRAS

Using a suite of N-body simulations in different Cold Dark Matter (CDM) scenarios, with cosmological constant (\LCDM) and without (OCDM, SCDM), we study the Hubble flow (\sigh) in Local Volumes (LV) around Local Group (LG) like objects found in these simulations, and compare the numerical results with the most recent observations. We show that \LCDM and OCDM models exhibit the same behavior of \sigh. Hence, we demonstrate that the observed coldness of the Hubble flow is not likely to be a manifestation of the dark energy, contrary to previous claims. The coldness does not constitute a problem by itself but it poses a problem to the standard \LCDM model only if the mean density within the Local Volume is greater than twice the mean matter cosmic density. The lack of blueshifted galaxies in the LV, outside of the LG can be considered as another manifestation of the coldness of the flow. Finally, we show that the main dynamical parameter that affects the coldness of the flow is the relative isolation of the LG, and the absence of nearby Milky Way like objects within a distance of about $3\mpc$., 13 pages, 8 figures, 5 tables. Submitted to MNRAS