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

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

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

We use high-resolution N-body simulations to follow the formation and evolution of tidal streams associated to dwarf spheroidal galaxies (dSphs). The dSph models are embedded in dark matter (DM) haloes with either a centrally-divergent 'cusp', or an homogeneous-density 'core'. In agreement with previous studies, we find that as tides strip the galaxy the evolution of the half-light radius and the averaged velocity dispersion follows well-defined tracks that are mainly controlled by the amount of mass lost. Crucially, the evolutionary tracks behave differently depending on the shape of the DM profile: at a fixed remnant mass, dSphs embedded in cored haloes have larger sizes and higher velocity dispersions than their cuspy counterparts. The divergent evolution is particularly pronounced in galaxies whose stellar component is strongly segregated within their DM halo and becomes more disparate as the remnant mass decreases. Our analysis indicates that the DM profile plays an important role in defining the internal dynamics of tidal streams. We find that stellar streams associated to cored DM models have velocity dispersions that lie systematically above their cuspy counterparts. Our results suggest that the dynamics of streams with known dSph progenitors may provide strong constraints on the distribution of DM on the smallest galactic scales.

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

In this paper, we study the formation and dynamical evolution of black hole-black hole (BH-BH) binaries in young star clusters (YSCs), by means of N-body simulations. The simulations include metallicity-dependent recipes for stellar evolution and stellar winds, and have been run for three different metallicities (Z = 0.01, 0.1 and 1 Zsun). Following recent theoretical models of wind mass-loss and core-collapse supernovae, we assume that the mass of the stellar remnants depends on the metallicity of the progenitor stars. We find that BH-BH binaries form efficiently because of dynamical exchanges: in our simulations, we find about 10 times more BH-BH binaries than double neutron star binaries. The simulated BH-BH binaries form earlier in metal-poor YSCs, which host more massive black holes (BHs) than in metal-rich YSCs. The simulated BH-BH binaries have very large chirp masses (up to 80 Msun), because the BH mass is assumed to depend on metallicity, and because BHs can grow in mass due to the merger with stars. The simulated BH-BH binaries span a wide range of orbital periods (10^-3-10^7 yr), and only a small fraction of them (0.3 per cent) is expected to merge within a Hubble time. We discuss the estimated merger rate from our simulations and the implications for Advanced VIRGO and LIGO.

We describe an algorithm for identifying ellipsoidal haloes in numerical simulations, and quantify how the resulting estimates of halo mass and shape differ with respect to spherical halo finders. Haloes become more prolate when fit with ellipsoids, the difference being most pronounced for the more aspherical objects. Although the ellipsoidal mass is systematically larger, this is less than 10% for most of the haloes. However, even this small difference in mass corresponds to a significant difference in shape. We quantify these effects also on the initial mass and deformation tensors, on which most models of triaxial collapse are based. By studying the properties of protohaloes in the initial conditions, we find that models in which protohaloes are identified in Lagrangian space by three positive eigenvalues of the deformation tensor are tenable only at the masses well-above $M_*$. The overdensity $\delta$ within almost any protohalo is larger than the critical value associated with spherical collapse (increasing as mass decreases); this is in good qualitative agreement with models which identify haloes requiring that collapse have occured along all three principal axes, each axis having turned around from the universal expansion at a different time. The distributions of initial values are in agreement with the simplest predictions associated with ellipsoidal collapse, assuming initially spherical protohaloes, collapsed around random positions which were sufficiently overdense. However, most protohaloes are not spherical and departures from sphericity increase as protohalo mass decreases. [Abridged]

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

We combine the physics of the ellipsoidal collapse model with the excursion set theory to study the shapes of dark matter halos. In particular, we develop an analytic approximation to the nonlinear evolution that is more accurate than the Zeldovich approximation; we introduce a planar representation of halo axis ratios, which allows a concise and intuitive description of the dynamics of collapsing regions and allows one to relate the final shape of a halo to its initial shape; we provide simple physical explanations for some empirical fitting formulae obtained from numerical studies. Comparison with simulations is challenging, as there is no agreement about how to define a non-spherical gravitationally bound object. Nevertheless, we find that our model matches the conditional minor-to-intermediate axis ratio distribution rather well, although it disagrees with the numerical results in reproducing the minor-to-major axis ratio distribution. In particular, the mass dependence of the minor-to-major axis distribution appears to be the opposite to what is found in many previous numerical studies, where low-mass halos are preferentially more spherical than high-mass halos. In our model, the high-mass halos are predicted to be more spherical, consistent with results based on a more recent and elaborate halo finding algorithm, and with observations of the mass dependence of the shapes of early-type galaxies. We suggest that some of the disagreement with some previous numerical studies may be alleviated if we consider only isolated halos.

We compare the apparent axial ratio distributions of Brightest Cluster Galaxies (BCGs) and normal ellipticals (Es) in our sample of 75 galaxy clusters from the WINGS survey. Most BCGs in our clusters (69%) are classified as cD galaxies. The sample of cDs has been completed by 14 additional cDs (non-BCGs) we found in our clusters. We find that: (i) Es have triaxial shape, the triaxiality sharing almost evenly the intrinsic axial ratios parameter space, with a weak preference for prolateness; (ii) the BCGs have triaxial shape as well. However, their tendence towards prolateness is much stronger than in the case of Es. Such a strong prolateness appears entirely due to the sizeable (dominant) component of cDs inside the WINGS sample of BCGs. In fact, while the 'normal'(non-cD) BCGs do not differ from Es, as far as the shape distribution is concerned, the axial ratio distribution of BCG_cD galaxies is found to support quite prolate shapes; (iii) our result turns out to be strongly at variance with the only similar, previous analysis by Ryden et al.(1993)(RLP93), where BCGs and Es were found to share the same axial ratio distribution; (iv) our data suggest that the above discrepancy is mainly caused by the different criteria that RLP93 and ourselves use to select the cluster samples, coupled with a preference of cDs to reside in powerful X-ray emitting clusters; (v) the GIF2 N-body results suggest that the prolateness of the BCGs (in particular the cDs) could reflect the shape of the associated dark matter halos.

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

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

We reassess the estimate of the cross-correlation of the spatial distribution of the NRAO VLA Sky Survey (NVSS) radio sources with that of Cosmic Microwave Background (CMB) anisotropies from the Wilkinson Microwave Anisotropy Probe (WMAP). This re-analysis is motivated by the fact that most previous studies adopted a redshift distribution of NVSS sources inconsistent with recent data. We find that the constraints on the bias-weighted redshift distribution, b(z)xN(z), of NVSS sources, set by the observed angular correlation function, w(theta), strongly mitigate the effect of the choice of N(z). If such constraints are met, even highly discrepant redshift distributions yield NVSS-WMAP cross-correlation functions consistent with each other within statistical errors. The models favoured by recent data imply a bias factor, b(z), decreasing with increasing z, rather than constant, as assumed by most previous analyses. As a consequence, the function b(z)xN(z) has more weight at z<1, i.e. in the redshift range yielding the maximum contribution to the ISW in a standard LambdaCDM cosmology. On the whole, the NVSS turns out to be better suited for ISW studies than generally believed, even in the absence of an observational determination of the redshift distribution. The NVSS-WMAP cross-correlation function is found to be fully consistent with the prediction of the standard LambdaCDM cosmology.

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

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

We present a precise estimate of the bulk virial scaling relation of halos formed via hierarchical clustering in an ensemble of simulated cold dark matter cosmologies. The result is insensitive to cosmological parameters, the presence of a trace, dissipationless gas component, and numerical resolution down to a limit of ~1000 particles. The dark matter velocity dispersion scales with total mass as log(sigma_{DM}(M,z)) = log(1082.9 +- 4.0 \kms) + (0.3361 +- 0.0026) log(h(z)M_{200}/10^{15} Msun), with h(z) the dimensionless Hubble parameter. At fixed mass, the velocity dispersion likelihood is nearly log-normal, with scatter sigma_{ln sigma} = 0.0426 +- 0.015, except for a tail to higher dispersions containing 10% of the population that are merger transients. We combine this relation with the halo mass function in LCDM models, and show that a low normalization condition, S_8 = sigma_8(Omega_m/0.3)^{0.35} \= 0.69, favored by recent WMAP and SDSS analysis requires that galaxy and gas specific energies in rich clusters be 50% larger than that of the underlying dark matter. Such large energetic biases are in conflict with the current generation of direct simulations of cluster formation. A higher normalization, S_8 = 0.80, alleviates this tension and implies that the hot gas fraction within r_{500} is (0.71 +- 0.09) h_{70}^{-3/2} Omega_b/\Omega_m, a value consistent with recent Sunyaev-Zel'dovich observations.

The polarization sensitivity of the upcoming millimetric observatories will open new possibilities for studying the properties of galaxy clusters and for using them as powerful cosmological probes. For this reason it is necessary to investigate in detail the characteristics of the polarization signals produced by their highly ionized intra-cluster medium (ICM). This work is focussed on the polarization effect induced by the ICM bulk motions, the so-called kpSZ signal, which has an amplitude proportional to the optical depth and to the square of the tangential velocity. In particular we study how this polarization signal is affected by the internal dynamics of galaxy clusters and what is its dependence on the physical modelling adopted to describe the baryonic component. This is done by producing realistic kpSZ maps starting from the outputs of two different sets of high-resolution hydrodynamical N-body simulations. The first set (17 objects) follows only non-radiative hydrodynamics, while for each of 9 objects of the second set we implement four different kinds of physical processes. Our results shows that the kpSZ signal turns out to be a very sensitive probe of the dynamical status of galaxy clusters. We find that major merger events can amplify the signal up to one order of magnitude with respect to relaxed clusters, reaching amplitude up to about 100 nuK. This result implies that the internal ICM dynamics must be taken into account when evaluating this signal because simplicistic models, based on spherical rigid bodies, may provide wrong estimates. Finally we find that the dependence on the physical modelling of the baryonic component is relevant only in the very inner regions of clusters.

We use high-resolution hydrodynamic re-simulations to investigate the properties of the thermal Sunyaev-Zel'dovich (SZ) effect from galaxy clusters. We compare results obtained using different physical models for the intracluster medium (ICM), and show how they modify the SZ emission in terms of cluster profiles and scaling relations. We also produce realistic mock observations to verify whether the results from hydrodynamic simulations can be confirmed. We find that SZ profiles depend marginally on the modelled physical processes, while they exhibit a strong dependence on cluster mass. The central and total SZ emission strongly correlate with the cluster X-ray luminosity and temperature. The logarithmic slopes of these scaling relations differ from the self-similar predictions by less than 0.2; the normalization of the relations is lower for simulations including radiative cooling. The observational test suggests that SZ cluster profiles are unlikely to be able to probe the ICM physics. The total SZ decrement appears to be an observable much more robust than the central intensity, and we suggest using the former to investigate scaling relations.

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

We present a study of the turbulent velocity fields in the Intra Cluster Medium of a sample of 21 galaxy clusters simulated by the SPH--code Gadget2, using a new numerical scheme where the artificial viscosity is suppressed outside shocks. The turbulent motions in the ICM of our simulated clusters are detected with a novel method devised to better disentangle laminar bulk motions from chaotic ones. We focus on the scaling law between the turbulent energy content of the gas particles and the total mass, and find that the energy in the form of turbulence scales approximatively with the thermal energy of clusters. We follow the evolution with time of the scaling laws and discuss the physical origin of the observed trends. The simulated data are in agreement with independent semi--analytical calculations, and the combination between the two methods allows to constrain the scaling law over more than two decades in cluster mass.

We examine the systematics affecting the X-ray mass estimators applied to a set of five simulated galaxy clusters. They have been processed through the X-ray Map Simulator, X-MAS, to provide Chandra-like long exposures that are analyzed to reconstruct the gas temperature, density, and mass profiles used as input. We find that at R_2500 the mass profile obtained via a direct application of the hydrostatic equilibrium equation is consistent within 1 sigma with the actual mass; although we notice this estimator shows high statistical errors due to high level of Chandra background. Instead, the poorness of the beta-model in describing the gas density profile makes the evaluated masses to be underestimated by \sim 40 per cent with respect to the true mass, both with an isothermal and a polytropic temperature profile. We also test ways to recover the mass by adopting an analytic mass model, such as those proposed by Navarro et al. (1997) and Rasia et al. (2004), and fitting the temperature profile expected from the hydrostatic equilibrium equation to the observed one. These methods and the one of the hydrostatic equilibrium equation provide a more robust mass estimation than the ones based on the beta-model. In the present work the main limitation for a precise mass reconstruction is to ascribe to the relatively high level of the background chosen to reproduce the Chandra one. After artificially reducing it by a factor of 100, we find that the estimated mass significantly underestimates the true mass profiles. This is manly due (i) to the neglected contribution of the gas bulk motions to the total energy budget and (ii) to the bias towards lower values of the X-ray temperature measurements because of the complex thermal structure of the emitting plasma.

We discuss an extended set of Tree+SPH simulations of galaxy clusters, with the goal of investigating the interplay between numerical resolution effects and star-formation/feedback processes. The simulated clusters span the mass range (0.1-2.3) 10^{15}Msun/h, with mass resolution varying by several decades. At the highest achieved resolution, we resolve the virial region of a Virgo-like cluster with more than 2 million gas particles and with at least as many dark-matter (DM) particles. Our resolution study confirms that, in the absence of an efficient feedback mechanism, runaway cooling leads to about 35 per cent of baryons in clusters to be locked up in long lived stars at our highest resolution, with no evidence of convergence. However, including feedback causes the fraction of cooled baryons to converge at about 15 per cent already at modest resolution. Feedback also stabilizes other gas-related quantities, such as radial profiles of entropy, gas density and temperature, against variations due to changes in resolution. We also investigate the influence of the gravitational force softening length, and that of numerical heating of the gas induced by two-body encounters between DM and lighter gas particles. We show that simulations where more DM than gas particles are used, show a significantly enhanced efficiency of star formation at z>3. Our results are important for establishing and delineating the regime of numerical reliability of the present generation of hydrodynamical simulations of galaxy clusters.

Smoothed particle hydrodynamics (SPH) employs an artificial viscosity to properly capture hydrodynamical shock waves. In its original formulation, the resulting numerical viscosity is large enough to suppress structure in the velocity field on scales well above the nominal resolution limit, and to damp the generation of turbulence by fluid instabilities. This could artificially suppress random gas motions in the intracluster medium (ICM), which are driven by infalling structures during the hierarchical structure formation process. We show that this is indeed the case by analysing results obtained with an SPH formulation where an individual, time-variable viscosity is used for each particle (Monaghan 1997). Using test calculations involving strong shocks, we demonstrate that this scheme captures shocks as well as the original formulation of SPH, but, in regions away from shocks, the numerical viscosity is much smaller. In a set of nine high-resolution simulations of cosmological galaxy cluster formation, we find that this low--viscosity formulation of SPH produces substantially higher levels of turbulent gas motions in the ICM, reaching a kinetic energy content in random gas motions (measured within a 1Mpc cube) of up to 5%-30% of the thermal energy content, depending on cluster mass. This has also significant effects on radial gas profile. We find a central flattening of the entropy profile and a reduction of the central gas density in the low--viscosity scheme. Interestingly, this tends to reduce the differences seen in SPH and adaptive mesh refinement simulations of cluster formation. Finally, invoking a model for particle acceleration by MHD waves driven by turbulence, we find efficient electron acceleration to power diffuse radio emission.

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

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. Here we show that, if the intra-cluster medium is thermally complex, T_ew fails at reproducing T_spec. We propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we also find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and beta-model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X--ray temperature function could bias low sigma_8 by 10-20 per cent. This may alleviate the tension between the value of sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

The fluctuations in the Cosmic Microwave Background (CMB) intensity due to the Sunyaev-Zeldovich (SZ) effect are the sum of a thermal and a kinetic contribution. Separating the two components to measure the peculiar velocity of galaxy clusters requires radio and microwave observations at three or more frequencies, and knowledge of the temperature T_e of the intracluster medium weighted by the electron number density. To quantify the systematics of this procedure, we extract a sample of 117 massive clusters at redshift z=0 from an N-body hydrodynamical simulation, with 2x480^3 particles, of a cosmological volume 192 Mpc/h on a side of a flat Cold Dark Matter model with Omega_0=0.3 and Lambda=0.7. Our simulation includes radiative cooling, star formation and the effect of feedback and galactic winds from supernovae. We find that (1) our simulated clusters reproduce the observed scaling relations between X-ray and SZ properties; (2) bulk flows internal to the intracluster medium affect the velocity estimate by less than 200 km/s in 93 per cent of the cases; (3) using the X-ray emission weighted temperature, as an estimate of T_e, can overestimate the peculiar velocity by 20-50 per cent, if the microwave observations do not spatially resolve the cluster. For spatially resolved clusters, the assumptions on the spatial distribution of the ICM, required to separate the two SZ components, still produce a velocity overestimate of 10-20 per cent, even with an unbiased measure of T_e. Thanks to the large size of our cluster samples, these results set a robust lower limit of 200 km/s to the systematic errors that will affect upcoming measures of cluster peculiar velocities with the SZ effect.

The detection of excess of soft X-ray or Extreme Ultraviolet (EUV) radiation, above the thermal contribution from the hot intracluster medium (ICM), has been a controversial subject ever since the initial discovery of this phenomenon. We use a large--scale hydrodynamical simulation of a concordance $\Lambda$CDM model, to investigate the possible thermal origin for such an excess in a set of 20 simulated clusters having temperatures in the range 1--7 keV. Simulated clusters are analysed by mimicking the observational procedure applied to ROSAT--PSPC data, which for the first time showed evidences for the soft X-ray excess. For cluster--centric distances $0.4< R/R_{\rm vir}< 0.7$ we detect a significant excess in most of the simulated clusters, whose relative amount changes from cluster to cluster and, for the same cluster, by changing the projection direction. In about 30 per cent of the cases, the soft X-ray flux is measured to be at least 50 per cent larger than predicted by the one--temperature plasma model. We find that this excess is generated in most cases within the cluster virialized regions. It is mainly contributed by low--entropy and high--density gas associated with merging sub--halos, rather than to diffuse warm gas. Only in a few cases the excess arises from fore/background groups observed in projection, while no evidence is found for a significant contribution from gas lying within large--scale filaments. We compute the distribution of the relative soft excess, as a function of the cluster--centric distance, and compare it with the observational result by Bonamente et al. (2003) for the Coma cluster. Similar to observations, we find that the relative excess increases with the distance from the cluster center, with no significant excess detected for $R<0.4R_{\rm vir}$. (abridged)

Theoretical studies of the physical processes in clusters of galaxies are mainly based on the results of numerical simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We show that the projected spectroscopic temperature of clusters obtained from X-ray observations is always lower than the emission-weighed temperature. This bias is related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. In real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps. We propose a new formula, the spectroscopic-like temperature function that better approximates the spectroscopic temperature, making simulations more directly comparable to observations

The thermal properties of hydrodynamical simulations of galaxy clusters are usually compared to observations by relying on the emission-weighted temperature T_ew, instead of on the spectroscopic X-ray temperature T_spec, which is obtained by actual observational data. In a recent paper Mazzotta et al. show that, if the cluster is thermally complex, T_ew fails at reproducing T_spec, and propose a new formula, the spectroscopic-like temperature, T_sl, which approximates T_spec better than a few per cent. By analyzing a set of hydrodynamical simulations of galaxy clusters, we find that T_sl is lower than T_ew by 20-30 per cent. As a consequence, the normalization of the M-T_sl relation from the simulations is larger than the observed one by about 50 per cent. If masses in simulated clusters are estimated by following the same assumptions of hydrostatic equilibrium and \beta--model gas density profile, as often done for observed clusters, then the M-T relation decreases by about 40 per cent, and significantly reduces its scatter. Based on this result, we conclude that using the observed M-T relation to infer the amplitude of the power spectrum from the X-ray temperature function could bias low \sigma_8 by 10-20 per cent. This may alleviate the tension between the value of \sigma_8 inferred from the cluster number density and those from cosmic microwave background and large scale structure.

(ABRIDGED) We analyze the X-ray properties of a sample of local and high redshift galaxy clusters extracted from a large cosmological hydrodynamical simulation. This simulation has been realized using the Tree+SPH code GADGET-2 for a LambdaCDM model. In our analysis, we consider only objects with T_ew >2 keV and adopt an approach that mimics observations, associating with each measurement an error comparable with recent observations and providing best-fit results via robust techniques. Within the clusters, baryons are distributed among (i) a cold neutral phase, with a relative contribution that increases from less than 1 to 3 per cent at higher redshift, (ii) stars which contribute with about 20 per cent and (iii) the X-ray emitting plasma that contributes by 80 (76) per cent at z=0 (1) to the total baryonic budget. A depletion of the cosmic baryon fraction of ~7 (at z=0) and 5 (at z=1) per cent is measured at the virial radius, R_vir, in good agreement with adiabatic hydrodynamical simulations. We confirm that, also at redshift >0.5, power-law relations hold between gas temperature, T, bolometric luminosity, L, central entropy, S, gas mass, M_gas, and total gravitating mass, M_tot and that these relations are steeper than predicted by simple gravitational collapse. A significant, negative evolution in the L-T and L-M_tot relations and positive evolution in the S-T relation are detected at 0.5 < z < 1 in this set of simulated galaxy clusters. This is partially consistent with recent analyses of the observed properties of z>0.5 X-ray galaxy clusters. By fixing the slope to the values predicted by simple gravitational collapse, we measure at high redshift normalizations lower by 10-40 per cent in the L-T, M_tot-T, M_gas-T, f_gas-T and L-M_tot relations than the observed estimates.

In this paper we present the analysis of Beppo-SAX observations of the cluster pairs A2061-A2067 and A2122-A2124, located in the Corona Borealis supercluster, which have been selected as candidate merging clusters. The aim of this work is to study the physics of the intracluster medium and to look for the possible presence of merging signatures. We derived the global temperatures and abundances and the temperature profiles and maps for these clusters. We do not find any significant evidence of interaction between the clusters forming the pairs, but we detect a candidate shock inside A2061. On the basis of the X-ray and optical properties of this cluster we propose a scenario in which a group is falling inside A2061. This interaction is in the phase in which the cores have not encountered yet and in which the formation of a shock is expected.

We study the properties of the diffuse light in galaxy clusters forming in a large hydrodynamical cosmological simulation of the Lambda-CDM cosmology. The simulation includes a model for radiative cooling, star formation in dense cold gas, and feedback by SN-II explosions. We select clusters having mass M>10^(14) h^(-1) Msun and study the spatial distribution of their star particles. While most stellar light is concentrated in gravitationally bound galaxies orbiting in the cluster potential, we find evidence for a substantial diffuse component, which may account for the extended halos of light observed around central cD galaxies. We find that more massive simulated clusters have a larger fraction of stars in the diffuse light than the less massive ones. The intracluster light is more centrally concentrated than the galaxy light, and the stars in the diffuse component are on average older than the stars in cluster galaxies, supporting the view that the diffuse light is not a random sampling of the stellar population in the cluster galaxies. We thus expect that at least ~10% of the stars in a cluster may be distributed as intracluster light, largely hidden thus far due to its very low surface brightness.

We study the efficiency of numerically simulated galaxy clusters for producing strong gravitational lensing events in different dark-energy cosmologies with constant and time-variable equation of state and we compare it with ``standard'' $\Lambda$CDM and OCDM models. Our main results are that: (1) the expected abundance of gravitational arcs with large length-to-width ratio depends on the equation of state of dark energy at the epoch of formation of the halo; (2) as found in previous studies, the strong lensing cross section of galaxy clusters is very sensitive to dynamical processes like mergers, whose impact however is different for different cosmologies, being stronger for models in which halos are less concentrated. As expected, the largest differences in the lensing optical depth occur at moderate and high redshift.

High-resolution, numerical simulations of 17 cluster-sized dark-matter haloes in eight different cosmologies with and without dynamical dark energy confirm the picture that core halo densities are imprinted early during their formation by the mean cosmological density. Quite independent of cosmology, halo concentrations have a log-normal distribution with a scatter of ~0.2 about the mean. We propose a simple scaling relation for halo concentrations in dark-energy cosmologies.

Theoretical studies of the physical processes guiding the formation and evolution of galaxies and galaxy clusters in the X-ray are mainly based on the results of numerical hydrodynamical N-body simulations, which in turn are often directly compared to X-ray observations. Although trivial in principle, these comparisons are not always simple. We demonstrate that the projected spectroscopic temperature of thermally complex clusters obtained from X-ray observations is always lower than the emission-weighed temperature, which is widely used in the analysis of numerical simulations. We show that this temperature bias is mainly related to the fact that the emission-weighted temperature does not reflect the actual spectral properties of the observed source. This has important implications for the study of thermal structures in clusters, especially when strong temperature gradients, like shock fronts, are present. Because of this bias, in real observations shock fronts appear much weaker than what is predicted by emission-weighted temperature maps, and may even not be detected. This may explain why, although numerical simulations predict that shock fronts are a quite common feature in clusters of galaxies, to date there are very few observations of objects in which they are clearly seen. To fix this problem we propose a new formula, the spectroscopic-like temperature function, and show that, for temperature larger than 3 keV, it approximates the spectroscopic temperature better than few per cent, making simulations more directly comparable to observations.

We analyze the dynamical and thermal evolution of dark matter and ICM in hydrodynamical Tree-SPH simulations of galaxy clusters. Starting from a sample of 17 high-resolution objects, with virial mass ranging from 3 10^{14} to 1.7 10^{15} M_sun, we follow the build-up of the systems in dark matter and hot gas through the repeated merging of satellites along their merging history trees. We measure the self-bound mass fraction of subhaloes as a function of time after the merging, estimate the satellite mean orbital properties as a function of the mass ratio with the main cluster at merging time, and study the evolution of their internal velocity dispersion, gas temperature and entropy as the substructure is disrupted by various dynamical processes, reaching eventually thermo-dynamic equilibrium in the gravitational potential of the main cluster. We model some relevant properties of subhalo orbits, as the time of the first pericentric and apocentric passages, and the typical distances and velocities at the corresponding times. This survival study can be used to interpret the dynamics of observed merging clusters; as an example we apply our results to the system 1E0657-56. We show that, in the light of our results, the most likely interpretation of the data for this cluster points to the merger of a small group with mass M\approx 1 10^{13} M_sun with a massive cluster with M\approx 1.3 10^{15} M_sun.

Using high-resolution non-radiative hydrodynamic simulations of galaxy clusters we obtain simple analytic formulae for DM and gas distribution in the spherical approximation. We derive fits for the DM density, velocity dispersion and velocity anisotropy. We use these models to test the dynamical equilibrium hypothesis through the Jeans equation: we find that this is satisfied to good accuracy. This result show that our fits constitute a self-consistent dynamical model. We then extend our analysis to the gas, studying its density, temperature and velocity structure, with no further hypothesis on the dynamical status or state equation. Gas and DM show self-similar density profiles down to 0.06 R_vir, while at smaller radii the gas produces a central core. Gas temperatures are almost isothermal out to 0.2 R_vir, then steeply decrease reaching at R_vir a value a factor of 2 lower. We find that the gas is not at rest inside R_vir: velocity dispersions are increasing functions of the radius, motions are isotropic to slightly tangential, and contribute to the total pressure support. We test our model using a generalization of the hydrostatic equilibrium equation, where the gas motion is properly taken into account. We find that our fits provide an accurate description of the system: the gas is in equilibrium and is a good tracer the overall cluster potential if all terms are taken into account, while simpler assumptions, as the beta-model, cause systematic mass underestimates. We also find that, if gas velocities are neglected, then a simple isothermal model fares better at large radii than a non-isothermal one. (abridged)

A generic prediction of hierarchical gravitational clustering models is that the distribution of halo formation times should depend relatively strongly on halo mass, massive haloes forming more recently, and depend only weakly, if at all, on the large scale environment of the haloes. We present a novel test of this assumption which uses a statistic that proves to be particularly well-suited to detecting and quantifying weak correlations with environment. We find that close pairs of haloes form at slightly higher redshifts than do more widely separated halo pairs, suggesting that haloes in dense regions form at slightly earlier times than do haloes of the same mass in less dense regions. The environmental trends we find are useful for models which relate the properties of galaxies to the formation histories of the haloes which surround them.

The most commonly used definition of halo formation is the time when a halo's most massive progenitor first contains at least half the final mass of its parent. Reasonably accurate formulae for the distribution of formation times of haloes of fixed mass have been available for some time. We use numerical simulations of hierarchical gravitational clustering to test the accuracy of formulae for the mass at formation. We also derive and test a formula for the joint distribution of formation masses and times. The structure of a halo is expected to be related to its accretion history. Our tests show that our formulae for formation masses and times are reasonably accurate, so we expect that they will aid future analytic studies of halo structure.

If the dark matter in the Universe is made of weakly self-interacting particles, they may self-annihilate and emit gamma-rays. We use high resolution numerical simulations to estimate directly the annihilation flux from the central regions of the Milky Way and from dark matter substructures in its halo. Although such estimates remain uncertain because of their strong dependence on the structure of the densest regions, our numerical experiments suggest that less direct calculations have overestimated the emission both from the centre and from halo substructure. We estimate a maximal enhancement of at most a factor of a few with respect to a smooth spherical halo of standard Navarro-Frenk-White (NFW) structure. We discuss detection strategies for the next generation of gamma-ray detectors and find that the annihilation flux may be detectable, regardless of uncertainties about the densest regions, for the annihilation cross-sections predicted by currently popular elementary particle models for the dark matter.

By means of high-resolution N-body simulations in a LambdaCDM cosmology, we verify that scaling relations similar to those observed for nearby galaxy clusters are also defined by their dark matter hosts; the slopes, however, are not the same. We then show that the scaling relations of galaxy clusters can be explained as the result of the cosmological collapse of density fluctuations at the appropriate scales, plus a systematic trend of the M/L ratio with cluster mass. The empirical fact that the exponent of the Faber-Jackson relation of elliptical galaxies is significantly different (higher) than that of clusters, force us to conclude that the galaxy scaling laws might derive from the cosmological collapse of density fluctuations at the epoch when galactic scales became non-linear, plus modifications afterward due to early-time dissipative merging.

We analyse the properties of substructures within dark matter halos (subhalos) using a set of high-resolution numerical simulations of the formation of structure in a Lambda-CDM Universe. Our simulation set includes 11 high-resolution simulations of massive clusters as well as a region of mean density, allowing us to study the spatial and mass distribution of substructures down to a mass resolution limit of 10^9 h^(-1)Mo. We also investigate how the properties of substructures vary as a function of the mass of the `parent' halo in which they are located. We find that the substructure mass function depends at most weakly on the mass of the parent halo and is well described by a power-law. The radial number density profiles of substructures are steeper in low mass halos than in high mass halos. More massive substructures tend to avoid the centres of halos and are preferentially located in the external regions of their parent halos. We also study the mass accretion and merging histories of substructures, which we find to be largely independent of environment. We find that a significant fraction of the substructures residing in clusters at the present day were accreted at redshifts z < 1. This implies that a significant fraction of present-day `passive' cluster galaxies should have been still outside the cluster progenitor and more active at z~1.

We study the impact of merger events on the strong lensing properties of galaxy clusters. Previous lensing simulations were not able to resolve dynamical time scales of cluster lenses, which arise on time scales which are of order a Gyr. In this case study, we first describe qualitatively with an analytic model how some of the lensing properties of clusters are expected to change during merging events. We then analyse a numerically simulated lens model for the variation in its efficiency for producing both tangential and radial arcs while a massive substructure falls onto the main cluster body. We find that: (1) during the merger, the shape of the critical lines and caustics changes substantially; (2) the lensing cross sections for long and thin arcs can grow by one order of magnitude and reach their maxima when the extent of the critical curves is largest; (3) the cross section for radial arcs also grows, but the cluster can efficiently produce this kind of arcs only while the merging substructure crosses the main cluster centre; (4) while the arc cross sections pass through their maxima as the merger proceeds, the cluster's X-ray emission increases by a factor of $\sim5$. Thus, we conclude that accounting for these dynamical processes is very important for arc statistics studies. In particular, they may provide a possible explanation for the arc statistics problem.

Numerical hydro-N-body simulations are very important tools for making theoretical predictions for the formation of galaxy clusters. They show that the atmospheres of clusters of galaxies have quite complex angular and thermal structures. The full understanding of the physical processes behind these features can be only achieved by direct comparison of observations to hydro-N-body simulations. Although simple in principle, these comparisons are not always trivial. In fact, real data are convolved with the instrument response which may substantially influence the apparent properties of the studied features. To overcome this problem we build the software package X-MAS devoted to simulate X-ray observations of galaxy clusters obtained from hydro-N-body simulations. In this paper we present how this software package works and discuss its application to the simulation of Chandra ACIS-S3 observations. We compare some of the main physical properties of the input data to the ones derived from simulated observations after performing a standard imaging and spectral analysis. We show that, if the thermal structure of the cluster along a particular line of sight is quite complex, the projected spectroscopic temperature obtained from the observation is significantly lower than the emission-weighed value inferred directly from hydrodynamical simulation. This implies that much attention must be paid in the theoretical interpretation of observational temperatures.

We present results on the X-ray properties of clusters and groups of galaxies, extracted from a large hydrodynamical simulation. We used the GADGET code to simulate a LambdaCDM model within a box of 192 Mpc/h on a side, with 480^3 dark matter particles and as many gas particles. The simulation includes radiative cooling, star formation and supernova feedback. The simulated M-T relation is consistent with observations once we mimic the procedure for mass estimates applied to real clusters. Also, with the adopted choices of Omega_m=0.3 and sigma_8=0.8 the resulting XTF agrees with observational determinations. The L-T relation also agrees with observations for clusters with T>2 keV, with no change of slope at the scale of groups. The entropy in central cluster regions is higher than predicted by gravitational heating alone, the excess being almost the same for clusters and groups. The simulated clusters appear to have suffered some overcooling, with f*~0.2, thus about twice as large as the value observed. Interestingly, temperature profiles are found to steadily increase toward cluster centers. They decrease in the outer regions, much like observational data do at r>0.2r_vir, while not showing an isothermal regime followed by a smooth temperature decline in the innermost regions.

We study the concentration parameters, their mass dependence and redshift evolution, of dark-matter halos in different dark-energy cosmologies with constant and time-variable equation of state, and compare them with "standard'' Lambda-CDM and OCDM models. We find that previously proposed algorithms for predicting halo concentrations can be well adapted to dark-energy models. When centred on the analytically expected values, halo concentrations show a log-normal distribution with a uniform standard deviation of ~0.2. The dependence of averaged halo concentrations on mass and redshift permits a simple fit of the form (1+z) c=c0 (M/M0)^a, with a~-0.1 throughout. We find that the cluster concentration depends on the dark energy equation of state at the cluster formation redshift z_{coll} through the linear growth factor D_+(z_{coll}). As a simple correction accounting for dark-energy cosmologies, we propose scaling c0 from Lambda-CDM with the ratio of linear growth factors, c0 -> c0 D_+(z_{coll})/D_{+,Lambda-CDM}(z_{coll}).

Like early-type galaxies, also nearby galaxy clusters define a Fundamental Plane, a luminosity-radius, and a luminosity-velocity dispersion relations, whose physical origin is still unclear. By means of high resolution N--body simulations of massive dark matter halos in a Lambda-CDM cosmology, we find that scaling relations similar to those observed for galaxy clusters are already defined by their dark matter hosts. The slopes however are not the same, and among the various possibilities in principle able to bring the simulated and the observed scaling relations in mutual agreement, we show that the preferred solution is a luminosity dependent mass-to-light ratio (M/L ~ L^0.3), that well corresponds to what inferred observationally. We then show that at galactic scales there is a conflict between the cosmological predictions of structure formation, the observed trend of the mass-to-light ratio in ellipticals, and the slope of their luminosity-velocity dispersion relation (that significantly differs from the analogous one followed by clusters). The conclusion is that the scaling laws of elliptical galaxies might be the combined result of the cosmological collapse of density fluctuations at the epoch when galactic scales became non-linear, plus important modifications afterward due to early-time dissipative merging. Finally, we briefly discuss the possible evolution of the cluster scaling relations with redshift.

Arc statistics is known to be a powerful cosmological tool. Numerical lensing simulations show that orders of magnitude differences in the number of {\em giant} arcs on the whole sky are expected in different cosmological models. In this paper, we discuss the analytic and numerical methods in arc statistics and show that analytic models fail to reproduce the efficiency for strong lensing of more realistic numerical cluster models. Then, we discuss two recent extensions of the lensing simulations, i.e. the effects of cD galaxies in the lensing clusters and the impact of cluster mergers on arc statistics. We show that cD galaxies can increase the lensing cross sections for long and thin arcs by perhaps up to $\sim50%$, while major mergers can change the cluster efficiency for producing such arcs by up to one order of magnitude.

We compare the structure and kinematics of the 11 known satellites of the Milky Way with high resolution simulations of the formation of its dark halo in a LambdaCDM universe. In contrast to earlier work, we find excellent agreement. The observed kinematics are exactly those predicted for stellar populations with the observed spatial structure orbiting within the most massive ``satellite'' substructures in our simulations. Less massive substructures have weaker potential wells than those hosting the observed satellites. If there is a halo substructure ``problem'', it consists in understanding why halo substructures have been so inefficient in making stars. Suggested modifications of dark matter properties (for example, self-interacting or warm dark matter) may well spoil the good agreement found for standard Cold Dark Matter.

The excursion set approach allows one to estimate the abundance and spatial distribution of virialized dark matter haloes efficiently and accurately. The predictions of this approach depend on how the nonlinear processes of collapse and virialization are modelled. We present simple analytic approximations which allow us to compare the excursion set predictions associated with spherical and ellipsoidal collapse. In particular, we present formulae for the universal unconditional mass function of bound objects and the conditional mass function which describes the mass function of the progenitors of haloes in a given mass range today. We show that the ellipsoidal collapse based moving barrier model provides a better description of what we measure in the numerical simulations than the spherical collapse based constant barrier model, although the agreement between model and simulations is better at large lookback times. Our results for the conditional mass function can be used to compute accurate approximations to the local-density mass function which quantifies the tendency for massive haloes to populate denser regions than less massive haloes. This happens because low density regions can be thought of as being collapsed haloes viewed at large lookback times, whereas high density regions are collapsed haloes viewed at small lookback times. Although we have only applied our analytic formulae to two simple barrier shapes, we show that they are, in fact, accurate for a wide variety of moving barriers. We suggest how they can be used to study the case in which the initial dark matter distribution is not completely cold.

We carry out ray tracing through five high resolution simulations of a galaxy cluster to study how its ability to produce giant gravitationally lensed arcs is influenced by the collision cross-section of its dark matter. In three cases typical dark matter particles in the cluster core undergo between 1 and 100 collisions per Hubble time; two more explore the long (``collisionless'') and short (``fluid'') mean free path limits. We study the size and shape distributions of arcs and compute the cross-section for producing ``extreme'' arcs of various sizes. Even a few collisions per particle modify the core structure enough to destroy the cluster's ability to produce long, thin arcs. For larger collision frequencies the cluster must be scaled up to unrealistically large masses before it regains the ability to produce giant arcs. None of our models with self-interacting dark matter (except the ``fluid'' limit) is able to produce radial arcs; even the case with the smallest scattering cross-section must be scaled to the upper limit of observed cluster masses before it produces radial arcs. Apparently the elastic collision cross-section of dark matter in clusters must be very small, below 0.1 cm^2 g^-1, to be compatible with the observed ability of clusters to produce both radial arcs and giant arcs.

We simulate the assembly of a massive rich cluster and the formation of its constituent galaxies in a flat, low-density universe. Our most accurate model follows the collapse, the star-formation history and the orbital motion of all galaxies more luminous than the Fornax dwarf spheroidal, while dark halo structure is tracked consistently throughout the cluster for all galaxies more luminous than the SMC. Within its virial radius this model contains about 2.0e7 dark matter particles and almost 5000 distinct dynamically resolved galaxies. Simulations of this same cluster at a variety of resolutions allow us to check explicitly for numerical convergence both of the dark matter structures produced by our new parallel N-body and substructure identification codes, and of the galaxy populations produced by the phenomenological models we use to follow cooling, star formation, feedback and stellar aging. This baryonic modelling is tuned so that our simulations reproduce the observed properties of isolated spirals outside clusters. Without further parameter adjustment our simulations then produce a luminosity function, a mass-to-light ratio, luminosity, number and velocity dispersion profiles, and a morphology-radius relation which are similar to those observed in real clusters. In particular, since our simulations follow galaxy merging explicitly, we can demonstrate that it accounts quantitatively for the observed cluster population of bulges and elliptical galaxies.

We present the results of a set of numerical simulations aiming at evaluating the effects of cluster galaxies on the arc statistics. At this goal we use nine different galaxy clusters obtained from N-body simulations. We mimic the presence of a population of galaxies inside each cluster, trying to reproduce the observed luminosity function and the spatial distribution. We compare the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of the properties of long arcs.

We study the formation of dark halos in a $\Lambda$CDM universe under the assumption that Cold Dark Matter particles have a finite cross-section for elastic collisions. We compare evolution when CDM mean free paths are comparable to halo sizes with the collisionless and fluid limits. We show that a few collisions per particle per Hubble time at halo centre can substantially affect the central density profile. Cross-sections an order of magnitude larger produce sufficient relaxation for rich clusters to develop core radii in the range 100-200 $h^{-1}$kpc. The structural evolution of halos is a competition between collisional relaxation caused by individual particle interactions and violent relaxation resulting from the infall and merging processes by which clusters grow. Although our simulations concentrate on systems of cluster size, we can scale our results to address the halo structure expected for dwarf galaxies. We find that collision cross-sections sufficiently large to significantly modify the cores of such galaxies produce cluster cores which are too large and/or too round to be consistent with observation. Thus the simplest model for self-interacting dark matter is unable to improve fits to published dwarf galaxy rotation curves without violating other observational constraints.

We study how the internal structure of dark halos is affected if Cold Dark Matter particles are assumed to have a large cross-section for elastic collisions. We identify a cluster halo in a large cosmological N-body simulation and resimulate its formation with progressively increasing resolution. We compare the structure found in the two cases where dark matter is treated as collisionless or as a fluid. For the collisionless case our results agree with those of other workers. Collisional dark matter results in a cluster which is more nearly spherical and has a more singular central density profile. Substructure within the cluster is only weakly suppressed relative to the collisionless case. The observed structure of dwarf galaxies argues against self-interacting dark matter if, as seems likely, intermediate cross-sections produce structure lying between the extremes we have simulated.

We present the results of a set of numerical simulations evaluating the effect of cluster galaxies on arc statistics. We perform a first set of gravitational lensing simulations using three independent projections for each of nine different galaxy clusters obtained from N-body simulations. The simulated clusters consist of dark matter only. We add a population of galaxies to each cluster, mimicking the observed luminosity function and the spatial galaxy distribution, and repeat the lensing simulations including the effects of cluster galaxies, which themselves act as individual lenses. Each galaxy is represented by a spherical Navarro, Frenk & White (1997) density profile. We consider the statistical distributions of the properties of the gravitational arcs produced by our clusters with and without galaxies. We find that the cluster galaxies do not introduce perturbations strong enough to significantly change the number of arcs and the distributions of lengths, widths, curvature radii and length-to-width ratios of long arcs. We find some changes to the distribution of short-arc properties in presence of cluster galaxies. The differences appear in the distribution of curvature radii for arc lengths smaller than 12'', while the distributions of lengths, widths and length-to-width ratios are significantly changed only for arcs shorter than 4''.

The Press--Schechter, excursion set approach allows one to make predictions about the shape and evolution of the mass function of bound objects. It combines the assumption that objects collapse spherically with the assumption that the initial density fluctuations were Gaussian and small. While the predicted mass function is reasonably accurate at the high mass end, it has more low mass objects than are seen in simulations of hierarchical clustering. This discrepancy can be reduced substantially if bound structures are assumed to form from an ellipsoidal, rather than a spherical collapse. In the spherical model, a region collapses if the initial density within it exceeds a threshold value, delta_sc. This value is independent of the initial size of the region, and since the mass of the collapsed object is related to its initial size, delta_sc is independent of final mass. In the ellipsoidal model, the collapse of a region depends on the surrounding shear field, as well as on its initial overdensity. Therefore, there is a relation between the density threshold value required for collapse, and the mass of the final object. We provide a fitting function to this delta_ec(m) relation for initially Gaussian fields which simplifies the inclusion of ellipsoidal dynamics in the excursion set approach. We discuss the relation between the excursion set predictions and the halo distribution in high resolution N-body simulations, and show that our simple parametrization of the ellipsoidal collapse model represents a significant improvement on the spherical model on an object-by-object basis. Finally, we show that the associated statistical predictions, the mass function and the large scale halo-to-mass bias relation, are also more accurate than the standard predictions.

Dark matter haloes are biased tracers of the underlying dark matter distribution. We use a simple model to provide a relation between the abundance of dark matter haloes and their spatial distribution on large scales. Our model shows that knowledge of the unconditional mass function alone is sufficient to provide an accurate estimate of the large scale bias factor. Then we use the mass function measured in numerical simulations of SCDM, OCDM and LCDM to compute this bias. Comparison with these simulations shows that this simple way of estimating the bias relation and its evolution is accurate for less massive haloes as well as massive ones. In particular, we show that haloes which are less/more massive than typical M* haloes at the time they form are more/less strongly clustered than formulae based on the standard Press-Schechter mass function predict.

We use N-body simulations to study the infall of dark matter haloes onto rich clusters of galaxies. After identification of all cluster progenitors in the simulations, we select those haloes which accrete directly onto the main cluster progenitor. We construct the mass function of these merging satellites, and calculate the main orbital parameters for the accreted lumps. The average circularity of the orbits is epsilon = 0.5, while either radial or almost circular orbits are equally avoided. More massive satellites move along slightly more eccentric orbits, with lower specific angular momentum and a smaller pericentre. We find that the infall of satellites onto the main cluster progenitor has a very anisotropic distribution. This anisotropy is to a large extent responsible for the shape and orientation of the final cluster and of its velocity ellipsoid. At the end of the simulations, the major axis of the cluster is aligned both with that of its velocity ellipsoid, and with the major axis of the ellipsoid defined by the satellite infall pattern, to 30 degrees on average. We also find that, in lower mass clusters, a higher fraction of the final virial mass is provided by small, dense satellites. These sink to the centre of the parent cluster and so enhance its central density. This mechanism is found to be partially responsible for the correlation between halo masses and characteristic overdensities, recently highlighted by Navarro, Frenk & White (1996).

We study the merging history of dark matter haloes that end up in rich clusters, using N-body simulations of a scale-free universe. We compare the predictions of the extended Press & Schechter (PS) formalism (Bond et al. 1991; Bower 1991; Lacey & Cole 1993) with several conditional statistics of the proto-cluster matter: the mass distribution and relative abundance of progenitor haloes at different redshifts, the infall rate of progenitors within the proto-cluster, the formation redshift of the most massive cluster progenitor, and the accretion rates of other haloes onto it. The high quality of our simulations allows an unprecedented resolution in the mass range of the studied distributions. We also present the global mass function for the same cosmological model. We find that the PS formalism and its extensions cannot simultaneously describe the global evolution of clustering and its evolution in a proto-cluster environment. The best-fit PS model for the global mass function is a poor fit to the statistics of cluster progenitors. This discrepancy is in the sense of underpredicting the number of high-mass progenitors at high redshift. Although the PS formalism can provide a good qualitative description of the global evolution of hierarchical clustering, particular attention is needed when applying the theory to the mass distribution of progenitor objects at high redshift.

Using high resolution cosmological N-body simulations, we investigate the survival of dark matter satellites falling into larger haloes. Satellites preserve their identity for some time after merging. We compute their loss of mass, energy and angular momentum as dynamical friction, tidal forces and collisions with other satellites dissolve them. We also analyse the evolution of their internal structure. Satellites with less than a few per cent the mass of the main halo may survive for several billion years, whereas larger satellites rapidly sink into the center of the main halo potential well and lose their identity. Penetrating encounters between satellites are frequent and may lead to significant mass loss and disruption. Only a minor fraction of cluster mass (10 per cent on average) is bound to substructure at most redshifts of interest. We discuss the application of these results to the survival and extent of dark matter haloes associated with cluster galaxies, and to interactions between galaxies in clusters. We find that 35-40 per cent of galaxy dark matter haloes are disrupted by the present time. The fraction of satellites undergoing close encounters is similar to the fraction of interacting or merging galaxies in clusters at moderate redshift.

(Shortened) We use N-body simulations to investigate the structure and dynamical evolution of dark matter halos in galaxy clusters. Our sample consists of nine massive halos from an EdS universe with scale free power spectrum and n = -1. Halos are resolved by ~20000 particles each, with a dynamical resolution of 20-25 kpc. Large scale tidal fields are included up to L=150 Mpc using background particles. The halo formation process can be characterized by the alternation of two dynamical configurations: a merging phase and a relaxation phase, defined by their signature on the evolution of the total mass and rms velocity. Halos spend on average one 1/3 of their evolution in the merging phase and 2/3 in the relaxation phase. Using this definition, we study the density profiles and their change during the halo history. The average density profiles are fitted by the NFW analytical model with an rms residual of 17% between the virial radius Rv and 0.01 Rv. The Hernquist (1990) profiles fits the same halos with an rms residual of 26%. The trend with mass of the scale radius of these fits is marginally consistent with that found by Cole & Lacey (1996): in comparison our halos are more centrally concentrated, and the relation between scale radius and halo mass is slightly steeper. We find a moderately large scatter in this relation, due both to dynamical evolution within halos and to fluctuations in the halo population. We analyze the dynamical equilibrium of our halos using the Jeans' equation, and find that on average they are approximately in equilibrium within their virial radius. Finally, we find that the projected mass profiles of our simulated halos are in very good agreement with the profiles of three rich galaxy clusters derived from strong and weak gravitational lensing observations.

We present a new method to add long wavelength power to an evolved $N$-body simulation, making use of the Zel'dovich (1970) approximation to change positions and velocities of particles. We describe the theoretical framework of our technique and apply it to a P$^3$M cosmological simulation performed on a cube of $100$ Mpc on a side, obtaining a new ``simulation'' of $800$ Mpc on a side. We study the effect of the power added by long waves by mean of several statistics of the density and velocity field, and suggest possible applications of our method to the study of the large-scale structure of the universe.

This review is a short introduction to numerical hydrodynamics in a cosmological context, intended for the non specialist. The main processes relevant to galaxy formation are first presented. The fluid equations are then introduced, and their implementation in numerical codes by Eulerian grid based methods and by Smooth Particle Hydrodynamics is sketched. As an application, I finally show some results from an SPH simulation of a galaxy cluster.

I use $N$-body simulations to investigate the morphology and the dynamical evolution and properties of dark matter halos in clusters of galaxies. My sample consists of nine massive halos, coming from an Einstein-de Sitter universe with scale free power spectrum and spectral index $n = -1$. Halos are resolved by 20,000 particles each, on average, and have a dynamical resolution of 20-25 kpc, as shown by extensive tests. I find that the {\em average} density profiles of the halos are fitted by the Navarro, Frenk and White (1995) analytical fit with a root mean square (rms) residual of 20\% within the Virial radius. The Hernquist (1990) analytical density profiles fits the same halos with an rms residual of 30\%. The projected mass profiles of the simulated halos are in very good agreement with the profiles of rich galaxy clusters derived from strong and weak gravitational lensing observations.

We analyse the large--scale velocity field obtained by N--body simulations of cold dark matter (CDM) models with non--Gaussian primordial density fluctuations, considering models with both positive and negative primordial skewness in the density fluctuation distribution. We study the velocity probability distribution and calculate the dependence of the bulk flow, one--point velocity dispersion and Cosmic Mach Number on the filtering size. We find that the sign of the primordial skewness of the density field provides poor discriminatory power on the evolved velocity field. All non--Gaussian models here considered tend to have lower velocity dispersion and bulk flow than the standard Gaussian CDM model, while the Cosmic Mach Number turns out to be a poor statistic in characterizing the models. Next, we compare the large--scale velocity field of a composite sample of optically selected galaxies as described by the Local Group properties, bulk flow, velocity correlation function and Cosmic Mach Number, with the velocity field of mock catalogues extracted from the N--body simulations. The comparison does not clearly permit to single out a best model: the standard Gaussian model is however marginally preferred by the maximum likelihood analysis.

The H magnitude aperture data published by the Aaronson et al. collaboration over a 10-year period is collected into a homogeneous data set of 1731 observations of 665 galaxies. 96\% of these galaxies have isophotal diameters and axial ratios determined by the Third Reference Catalogue of Bright Galaxies (de Vaucouleurs et al. 1991), the most self-consistent set of optical data currently available. The precepts governing the optical data in the RC3 are systematically different from those of the Second Reference Catalogue (de Vaucouleurs, de Vaucouleurs and Corwin 1976), which were used by Aaronson et al. for their original analyses of galaxy peculiar motions. This in turn leads to systematic differences in growth curves and fiducial H magnitudes, prompting the present recalibration of the near-infrared Tully Fisher relationship. New optically-normalized H magnitude growth curves are defined for galaxies of types S0 to Im, from which new values of fiducial H magnitudes, H$\rm _{-0.5}^g$, are measured for the 665 galaxies. A series of internal tests show that these four standard growth curves are defined to an accuracy of 0.05 mag over the interval $\rm -1.5 \leq \log (A/D_g) \leq -0.2$. Comparisons with the Aaronson et al. values of diameters, axial ratios and fiducial H magnitudes show the expected differences, given the different definitions of these parameters. The values of H$\rm _{-0.5}^g$ are assigned quality indices; a value of 1 indicates an accuracy of $<$0.2 magnitude; Quality 2 indicates an accuracy of 0.2 to 0.35 mag; and Quality 3 indicates $>$0.35 mag. Revised values of corrected HI velocity widths are also given, based on the new set of axial ratios defined by the RC3.

We analyze the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We perform a statistical analysis, by calculating the bulk flow, Cosmic Mach Number and velocity correlation function for this sample and for mock catalogs extracted from a set of N--body simulations. We run four cold dark matter (CDM) simulations: two tilted models (with spectral index $n=0.6$ and $n=0.8$), the standard model ($n=1$) and a ``blue" one ($n=1.2$), with different values of the linear bias parameter $b$. By means of a Maximum Likelihood analysis we estimate the ability of our models to fit the observations, as measured by the above statistics, and to reproduce the Local Group properties. On the basis of this analysis we conclude that the best model is the unbiased standard model $(n,~b)=(1,~1)$, even though the overall flatness of the joint likelihood function implies that one cannot strongly discriminate models in the range $0.8 \le n \le 1$, and $1 \le b \le 1.5$. Models with $b \geq 2.5$ are rejected at the $95\%$ confidence level. For $n=0.8$ the values of $b$ preferred by the present analysis, together with the {\em COBE} data, require a negligible contribution to $\Delta T/T$ by gravitational waves. Finally, the blue model, normalized to {\em COBE}, does not provide a good fit to the velocity data.

We study the cosmic peculiar velocity field as traced by a sample of 1184 spiral, elliptical and S0 galaxies, grouped in 704 objects. We carry out a statistical analysis, by calculating bulk flows and velocity correlation functions for this sample and for mock catalogs which we extract from N--body simulations. For the simulations we consider tilted (i.e. with spectral index $n\leq 1$) CDM models with different values of the linear bias parameter $b$. By mean of a maximum likelihood analysis we estimate the ability of the models in fitting the observations as measured by the above statistics and in reproducing the Local Group properties.