Search

We measure the mass function of dark matter halos in a large set of collisionless cosmological simulations of flat LCDM cosmology and investigate its evolution at z<~2. Halos are identified as isolated density peaks, and their masses are measured within a series of radii enclosing specific overdensities. We argue that these spherical overdensity masses are more directly linked to cluster observables than masses measured using the friends-of-friends algorithm (FOF), and are therefore preferable for accurate forecasts of halo abundances. Our simulation set allows us to calibrate the mass function at z=0 for virial masses in the range 10^{11} Msol/h < M < 10^{15} Msol/h, to <~ 5%. We derive fitting functions for the halo mass function in this mass range for a wide range of overdensities, both at z=0 and earlier epochs. In addition to these formulae, which improve on previous approximations by 10-20%, our main finding is that the mass function cannot be represented by a universal fitting function at this level of accuracy. The amplitude of the "universal" function decreases monotonically by ~20-50%, depending on the mass definition, from z=0 to 2.5. We also find evidence for redshift evolution in the overall shape of the mass function., Comment: 19 emulateapj pages, 15 figures, submitted to Apj

We explore the dependence of the radial alignment of subhalos on the mass of the host halo they orbit in. As the effect is seen on a broad range of scales including massive clusters as well as galactic systems it only appears natural to explore this phenomenon by means of cosmological simulations covering the same range in masses. We have 25 well resolved host dark matter halos at our disposal ranging from 10^15 Msun/h down to 10^12 Msun/h each consisting of order of a couple of million particles within the virial radius. We observe that subhalos tend to be more spherical than isolated objects. Both the distributions of sphericity and triaxiality of subhalos are Gaussian distributed with peak values of s approx. 0.80 and T approx. 0.56, irrespective of host mass. Interestingly we note that the radial alignment is independent of host halo mass and the distribution of \cos\theta (i.e. the angle between the major axis E_a of each subhalo and the radius vector of the subhalo in the reference frame of the host) is well fitted by a simple power law P(\cos\theta) proportional to \cos^4\theta with the same fitting parameters for all host halos., Comment: 5 pages, 2 figures. MNRAS Letter, in press

To study the local Hubble flow, we have run constrained dark matter (DM) simulations of the Local Group (LG) in the concordance LCDM and OCDM cosmologies, with identical cosmological parameters apart from the Lambda term. The simulations were performed within a computational box of 64 h^{-1}Mpc centred on the LG. The initial conditions were constrained by the observed peculiar velocities of galaxies and positions of X-ray nearby clusters of galaxies. The simulations faithfully reproduce the nearby large scale structure, and in particular the Local Supercluster and the Virgo cluster. LG-like objects have been selected from the DM halos so as to closely resemble the dynamical properties of the LG. Both the LCDM and OCDM simulations show very similar local Hubble flow around the LG-like objects. It follows that, contrary to recent statements, the dark energy (DE) does not manifest itself in the local dynamics., Comment: 10 pages, 7 figures, submitted to MNRAS

The cosmic ultraviolet background (UVB) heats the intergalactic medium (IGM), as a result the gas in dark matter halos below a certain mass is too hot to cool within a Hubble time. The UVB effectively suppresses the formation of dwarf galaxies. Using high resolution cosmological hydrodynamical simulations we show that photo heating leads to small baryon fractions in halos below ~ 6x10^9 h^{-1}M_sun, independent of the cosmic environment. The simulations are carried out assuming a homogeneous UVB with flux densities as given by Haardt & Madau (1996). A halo may stop to condense gas significantly after the universe is reionised, namely when its mass falls below the characteristic mass scale set by the photo heating. Assuming a spherical halo model we derive this characteristic mass analytically and identify the main mechanisms that prevent the gas from cooling in small halos. The theoretically derived characteristic mass is smaller than the one obtained from observations. Increasing the energy per ionising photon by a factor between four and eight would be sufficient to reconcile both. This is equivalent to an average temperature of the IGM of ~ 10^4 K. In this sense the faint end of the luminosity function may serve as a calorimeter for the IGM., Comment: To appear in Proceedings of IAU Symp #244, "Dark Galaxies and Lost Baryons", June 2007, 5 pages including 3 figures

We present the mass and X-ray temperature functions derived from a sample of more than 15,000 galaxy clusters of the MareNostrum Universe cosmological SPH simulations. In these simulations, we follow structure formation in a cubic volume of 500/h Mpc on a side assuming cosmological parameters consistent with either the first or third year WMAP data and gaussian initial conditions. We compare our numerical predictions with the most recent observational estimates of the cluster X-ray temperature functions and find that the low normalization cosmological model inferred from the 3 year WMAP data results is barely compatible with the present epoch X-ray cluster abundances. We can only reconcile the simulations with the observational data if we assume a normalization of the Mass-Temperature relation which is a factor of 2.5--3 smaller than our non-radiative simulations predict. This deviation seems to be too large to be accounted by the effects of star formation or cooling in the ICM, not taken into account in these simulations., Comment: 4 pages, 3 figures. Accepted for publication in The Astrophysical Journal Letters

We study the distinct effects of Dark Matter and Dark Energy on the future evolution of nearby large scale structures using constrained N-body simulations. We contrast a model of Cold Dark Matter and a Cosmological Constant (LCDM) with an Open CDM (OCDM) model with the same matter density Omega_m =0.3 and the same Hubble constant h=0.7. Already by the time the scale factor increased by a factor of 6 (29 Gyr from now in LCDM; 78 Gyr from now in OCDM) the comoving position of the Local Group is frozen. Well before that epoch the two most massive members of the Local Group, the Milky Way and Andromeda, will merge. However, as the expansion rates of the scale factor in the two models are different, the Local Group will be receding in physical coordinates from Virgo exponentially in a LCDM model and at a roughly constant velocity in an OCDM model. More generally, in comoving coordinates the future large scale structure will look like a sharpened image of the present structure: the skeleton of the cosmic web will remain the same, but clusters will be more `isolated' and the filaments will become thinner. This implies that the long-term fate of large scale structure as seen in comoving coordinates is determined primarily by the matter density. We conclude that although the LCDM model is accelerating at present due to its Dark Energy component while the OCDM model is non accelerating, their large scale structure in the future will look very similar in comoving coordinates., Comment: 14 pages and 4 figures. New version, accepted for publication in JCAP

It has been recently claimed that there is no correlation between the distribution of galaxies and their peculiar velocities within the Local Volume (LV), namely a sphere of R=7/h Mpc around the Local Group (LG). It has been then stated that this implies that either locally dark matter is not distributed in the same way as luminous matter, or peculiar velocities are not due to fluctuations in mass. To test that statement a set of constrained N-body cosmological simulations, designed to reproduce the main observed large scale structure, have been analyzed. The simulations were performed within the flat-Lambda, open and flat matter only CDM cosmogonies. Two unconstrained simulations of the flat-Lambda and open CDM models were performed for comparison. LG-like objects have been selected so as to mimic the real LG environment. The local gravitational field due to all halos found within each LV is compared with the exact gravitational field induced by all matter in the simulation. We conclude that there is no correlation between the exact and the local gravitational field obtained by pairwise newtonian forces between halos. Moreover, the local gravitational field is uncorrelated with the peculiar velocities of halos. The exact gravitational field has a linear correlation with peculiar velocities but the proportionality constant relating the velocity with gravitational field falls below the prediction of the linear theory. Upon considering all matter inside the LVs, the exact and local gravitational accelerations show a much better correlation, but with a considerable scatter independent on the cosmological models. The main conclusion is that the lack of correlation between the local gravitation and the peculiar velocity fields around LG-like objects is naturally expected in the CDM cosmologies., Comment: 10 pages, 19 figures. Accepted for publication in MNRAS

The MareNostrum Universe is one of the largest cosmological SPH simulation done so far. It consists of $1024^3$ dark and $1024^3$ gas particles in a box of 500 $h^{-1}$ Mpc on a side. Here we study the shapes and spins of the dark matter and gas components of the 10,000 most massive objects extracted from the simulation as well as the gas fraction in those objects. We find that the shapes of objects tend to be prolate both in the dark matter and gas. There is a clear dependence of shape on halo mass, the more massive ones being less spherical than the less massive objects. The gas distribution is nevertheless much more spherical than the dark matter, although the triaxiality parameters of gas and dark matter differ only by a few percent and it increases with cluster mass. The spin parameters of gas and dark matter can be well fitted by a lognormal distribution function. On average, the spin of gas is 1.4 larger than the spin of dark matter. We find a similar behavior for the spins at higher redshifts, with a slightly decrease of the spin ratios to 1.16 at $z=1.$ The cosmic normalized baryon fraction in the entire cluster sample ranges from $Y_b = 0.94$, at $z=1$ to $Y_b = 0.92$ at $z=0$. At both redshifts we find a slightly, but statistically significant decrease of $Y_b$ with cluster mass., Comment: 7 pages, 6 figures. Accepted for publication in The Astrophysical Journal

A large scale SPH+N-body simulation (GADGET) of the concordance LCDM universe is used to investigate orientation and angular momentum of galaxy clusters at z=0 in connection with their recent accretion histories. The basic cluster sample comprises the 3000 most massive friends-of-friends halos found in the 500 Mpc/h simulation box. Two disjoint sub-samples are constructed, using the mass ratio of the two most massive progenitors at z=0.5 m_2 / m_1 (m_1 < m_2), namely a recent major merger sample and a steady accretion mode sample. The mass of clusters in the merger sample is on average ~43% larger than the mass of the two progenitors (m_1 + m_2), whereas in the steady accretion mode sample a smaller increase of ~25% is found. The separation vector connecting the two most massive progenitor halos at z=0.5 is strongly correlated with the orientation of the cluster at z=0. The angular momentum of the clusters in the recent major merger sample tends to be parallel to orbital angular momentum of the two progenitors, whereas the angular momentum of the steady accretion mode sample is mainly determined by the angular momentum of the most massive progenitor. The long range correlations for the major and the minor principal axes of cluster pairs extend to distances of ~100 Mpc/h. Weak angular momentum correlations are found for distances < 20 Mpc/h. Within these ranges the major axes tend to be aligned with the connecting line of the cluster pairs whereas minor axes and angular momenta tend to be perpendicular to this line. A separate analysis of the two sub-samples reveals that the long range correlations are independent of the mass accretion mode. Thus orientation and angular momentum of galaxy clusters is mainly determined by the accretion along the filaments independently of the particular accretion mode., Comment: 11 pages, 8 figures, replaced to match version accepted for publication in MNRAS

We study structure formation in cosmological void regions using high-resolution hydrodynamical simulations. Despite being significantly underdense, voids are populated abundantly with small dark matter halos which should appear as dwarf galaxies if their star formation is not suppressed significantly. We here investigate to which extent the cosmological UV-background photo-evaporates baryons out of halos of dwarf galaxies, and thereby limits their cooling and star formation rates. Assuming a Haardt & Madau UV-background with reionisation at redshift z=6, our samples of simulated galaxies show that halos with masses below a characteristic mass of M_c(z=0) = 6.5 x 10^9 h^{-1} M_sun are baryon-poor, but in general not completely empty, because baryons that are in the condensed cold phase or are already locked up in stars resist evaporation. In halos with mass M < M_c, we find that photo-heating suppresses further cooling of gas. The redshift and UV-background dependent characteristic mass M_c(z) can be understood from the equilibrium temperature between heating and cooling at a characteristic overdensity of \delta ~ 1000. If a halo is massive enough to compress gas to this density despite the presence of the UV background, gas is free to `enter' the condensed phase and cooling continues in the halo, otherwise it stalls. By analysing the mass accretion histories of dwarf galaxies in voids, we show that they can build up a significant amount of condensed mass at early times before the epoch of reionisation. Later on, the amount of mass in this phase remains roughly constant, but the masses of the dark matter halos continue to increase. (abridged), Comment: revised version as accepted by MNRAS, 15 pages, 15 figures, new simulation results and a significantly extended discussion have been included

Motivated by the discovery of a number of radio relics we investigate the fate of fossil radio plasma during a merger of clusters of galaxies using cosmological smoothed-particle hydrodynamics simulations. Radio relics are extended, steep-spectrum radio sources that do not seem to be associated with a host galaxy. One proposed scenario whereby these relics form is through the compression of fossil radio plasma during a merger between clusters. The ensuing compression of the plasma can lead to a substantial increase in synchrotron luminosity and this appears as a radio relic. Our simulations show that relics are most likely to be found at the periphery of the cluster at the positions of the outgoing merger shock waves. Relics are expected to be very rare in the centre of the cluster where the life time of relativistic electrons is short and shock waves are weaker than in the cooler, peripheral regions of the cluster. These predictions can soon be tested with upcoming low-frequency radio telescopes., Comment: 6 pages, 3 figures (high-resolution figures are available at http://www.faculty.iu-bremen.de/~mbrueggen/), submitted to MNRAS

We present estimates of the energy input from supernovae (SNe) into the intergalactic medium using (i) recent measurements of Si and Fe abundances in the intracluster medium (ICM) and (ii) self-consistent gasdynamical galaxy formation simulations that include processes of cooling, star formation, SNe feedback, and a multi-phase model of the interstellar medium. We estimate the energy input from observed abundances using two different assumptions: (i) spatial uniformity of metal abundances in the ICM and (ii) radial abundance gradients. We show that these two cases lead to energy input estimates which are different by an order of magnitude, highlighting a need for observational data on large-scale abundance gradients in clusters. Analysis of galaxy formation results and estimates from observed Fe and Si abundances indicates that the SNe energy input can be important for heating of the entire ICM (providing energy of ~1 keV per particle) only if the ICM abundances are uniform and the efficiency of gas heating by SN explosions is close to 100% (implying that all of the initial kinetic energy of the explosion goes into heating of the ICM). We conclude that unless these most favorable conditions are met, SNe alone are unlikely to provide sufficient energy input to heat all of the cluster ICM and may need to be supplemented or even substituted by some other heating process(es). (Abridged), Comment: submitted to MNRAS, LaTeX, 15 pages, 5 figures, 2 tables

Formation of luminous matter in the Universe is a complicated process, which includes many processes and components. It is the vastly different scales involved in the process (from star formation on few parsec scales to galaxy clusters and superclusters on megaparsecs scales) and numerous ill-understood processes, which make the whole field a maze of unsolved, but exciting problems. We present new approximations for numerical treatment of multiphase ISM forming stars. The approximations were tested and calibrated using N-body+fluid numerical simulations. We specifically target issues related with effects of unresolved lumpinesses of the gas. We show that the degree of freedom is much smaller than naively expected because of self-regulating nature of the process of global star formation. One of the problems of numerical simulations is related with the uncertainties of approximation of the supernovae (SN) feedback. It is often assumed that the feedback is mainly due to momentum transferred by SN in to the ISM. We argue that this may not be true. We present a realistic example of gas actively forming stars with short cooling time for which the SF feedback is important, but the kinetic energy of the gas motion due to SN is negligibly small as compared with the thermal energy of gas., Comment: 5 pages, 2 PS figures, uses epsf.sty. Talk presented at the 12-th Postdam Cosmology Workshop, ``Large Scale Structure: Tracks and Traces'', World Scientific 1998. Available at http://astro.nmsu.edu/~aklypin/index.html