Search

The radial density profile of dark matter structures has been observed to have an almost universal behaviour in numerical simulations, however, the physical reason for this behaviour remains unclear. It has previously been shown that if the pseudo phase-space density, rho/sigma_d^epsilon, is a beautifully simple power-law in radius, with the "golden values" epsilon=3 and d=r (i.e., the phase-space density is only dependent on the radial component of the velocity dispersion), then one can analytically derive the radial variation of the mass profile, dispersion profile etc. That would imply, if correct, that we just have to explain why rho/sigma^3_r ~r^{-alpha}, and then we would understand everything about equilibrated DM structures. Here we use a set of simulated galaxies and clusters of galaxies to demonstrate that there are no such golden values, but that each structure instead has its own set of values. Considering the same structure at different redshifts shows no evolution of the phase-space parameters towards fixed points. There is also no clear connection between the halo virialized mass and these parameters. This implies that we still do not understand the origin of the profiles of dark matter structures., Comment: 4 pages, 3 figures, accepted for publication in ApJL

It has been pointed out that the Galactic satellites all have a common mass around 1e7 Msun within 300 pc (M03), while they span almost four order of magnitudes in luminosity (Mateo et al. 1993, Strigari et al. 2008). It is argued that this may reflect a specific scale for galaxy formation or a scale for dark matter clustering. Here we use numerical simulations coupled with a semi-analytic model for galaxy formation, to predict the central mass and luminosity of galactic satellites in the LCDM model. We show that this common mass scale can be explained within the Cold Dark Matter scenario when the physics of galaxy formation is taken into account. The narrow range of M03 comes from the narrow distribution of circular velocities at time of accretion (peaking around 20 km/s) for satellites able to form stars and the not tight correlation between halo concentration and circular velocity. The wide range of satellite luminosities is due to a combination of the mass at time of accretion and the broad distribution of accretion redshifts for a given mass. This causes the satellites baryonic content to be suppressed by photo-ionization to very different extents. Our results favor the argument that the common mass M03 reflects a specific scale (circular velocity ~ 20 km/s) for star formation., Comment: 5 pages, 3 figures. References added, discussion enlarged, new version of Figure 3. Minor changes to match the version accepted for publication on ApJ Letters

Accurate predictions on non--linear power spectra, at various redshift z, will be a basic tool to interpret cosmological data from next generation mass probes, so obtaining key information on Dark Energy nature. This calls for high precision simulations, covering the whole functional space of w(z) state equations and taking also into account the admitted ranges of other cosmological parameters; surely a difficult task. A procedure was however suggested, able to match the spectra at z=0, up to k~3, hMpc^{-1}, in cosmologies with an (almost) arbitrary w(z), by making recourse to the results of N-body simulations with w = const. In this paper we extend such procedure to high redshift and test our approach through a series of N-body gravitational simulations of various models, including a model closely fitting WMAP5 and complementary data. Our approach detects w= const. models, whose spectra meet the requirement within 1% at z=0 and perform even better at higher redshift, where they are close to a permil precision. Available Halofit expressions, extended to (constant) w \neq -1 are unfortunately unsuitable to fit the spectra of the physical models considered here. Their extension to cover the desired range should be however feasible, and this will enable us to match spectra from any DE state equation., Comment: method definitely improved in semplicity and efficacy,accepted for publication on JCAP

We investigate the effects of changes in the cosmological parameters between the WMAP 1st, 3rd, and 5th year results on the structure of dark matter haloes. We use a set of simulations that cover 5 decades in halo mass ranging from the scales of dwarf galaxies (V_c ~30 km/s) to clusters of galaxies (V_c ~ 1000 km/s). We find that the concentration mass relation is a power law in all three cosmologies. However the slope is shallower and the zero point is lower moving from WMAP1 to WMAP5 to WMAP3. For haloes of mass log(M_200/Msun) = 10, 12, and 14 the differences in the concentration parameter between WMAP1 and WMAP3 are a factor of 1.55, 1.41, and 1.29, respectively. As we show, this brings the central densities of dark matter haloes in good agreement with the central densities of dwarf and low surface brightness galaxies inferred from their rotation curves, for both the WMAP3 and WMAP5 cosmologies. We also show that none of the existing toy models for the concentration-mass relation can reproduce our simulation results over the entire range of masses probed. In particular, the model of Bullock et al (B01) fails at the higher mass end (M > 1e13 Msun), while the NFW model of Navarro, Frenk & White (1997) fails dramatically at the low mass end (M < 1e12 Msun). We present a new model, based on a simple modification of that of B01, which reproduces the concentration-mass relations in our simulations over the entire range of masses probed (1e10 Msun < M < 1e15 Msun). Haloes in the WMAP3 cosmology (at a fixed mass) are more flatted compared to the WMAP1 cosmology, with a medium to long axis ration reduced by ~10 %. Finally, we show that the distribution of halo spin parameters is the same for all three cosmologies., Comment: 16 pages, 16 figures, references updated, minor changes. Accepted for publication on MNRAS. WMAP5 simulations available upon request

Ultra compact dwarf galaxies (UCDs) have similar properties as massive globular clusters or the nuclei of nucleated galaxies. Recent observations suggesting a high dark matter content and a steep spatial distribution within groups and clusters provide new clues as to their origins. We perform high-resolution N-body / smoothed particle hydrodynamics simulations designed to elucidate two possible formation mechanisms for these systems: the merging of globular clusters in the centre of a dark matter halo, or the massively stripped remnant of a nucleated galaxy. Both models produce density profiles as well as the half light radii that can fit the observational constraints. However, we show that the first scenario results to UCDs that are underluminous and contain no dark matter. This is because the sinking process ejects most of the dark matter particles from the halo centre. Stripped nuclei give a more promising explanation, especially if the nuclei form via the sinking of gas, funneled down inner galactic bars, since this process enhances the central dark matter content. Even when the entire disk is tidally stripped away, the nucleus stays intact and can remain dark matter dominated even after severe stripping. Total galaxy disruption beyond the nuclei only occurs on certain orbits and depends on the amount of dissipation during nuclei formation. By comparing the total disruption of CDM subhaloes in a cluster potential we demonstrate that this model also leads to the observed spatial distribution of UCDs which can be tested in more detail with larger data sets., Comment: 8 pages, 8 figures, final version accepted for publication in MNRAS

Warm Dark Matter (WDM) has been invoked to resolve apparent conflicts of Cold Dark Matter (CDM) models with observations on subgalactic scales. In this work we provide a new and independent lower limit for the WDM particle mass (e.g. sterile neutrino) through the analysis of image fluxes in gravitationally lensed QSOs. Starting from a theoretical unperturbed cusp configuration we analyze the effects of intergalactic haloes in modifying the fluxes of QSO multiple images, giving rise to the so-called anomalous flux ratio. We found that the global effect of such haloes strongly depends on their mass/abundance ratio and it is maximized for haloes in the mass range $10^6-10^8 \Msun$. This result opens up a new possibility to constrain CDM predictions on small scales and test different warm candidates, since free streaming of warm dark matter particles can considerably dampen the matter power spectrum in this mass range. As a consequence, while a ($\Lambda$)CDM model is able to produce flux anomalies at a level similar to those observed, a WDM model, with an insufficiently massive particle, fails to reproduce the observational evidences. Our analysis suggests a lower limit of a few keV ($m_{\nu} \sim 10$) for the mass of warm dark matter candidates in the form of a sterile neutrino. This result makes sterile neutrino Warm Dark Matter less attractive as an alternative to Cold Dark Matter, in good agreement with previous findings from Lyman-$\alpha$ forest and Cosmic Microwave Background analysis., Comment: One equation added, typo in eq: 5 corrected, minor changes to match the accepted version by MNRAS

The reionization of the Universe has profound effects on the way galaxies form and on their observed properties at later times. Of particular importance is the relative timing of the reionization history of a region and its halo assembly history, which can affect the nature of the first stars formed in that region, the properties and radial distribution of its stellar halo, globular cluster population and its satellite galaxies. We distinguish two basic cases for the reionization of a halo - internal reionization, whereby the stars forming in situ reionize their host galaxy, and external reionization, whereby the progenitor of a galaxy is reionized by external radiation before its own stars are able to form in sufficient numbers. We use a set of large-scale radiative transfer and structure formation simulations, based on cosmologies derived from both WMAP 1-year and WMAP 3-year data, to evaluate the mean reionization redshifts and the probability of internal/external reionization for Local Group-like systems, galaxies in the field and central cD galaxies in clusters. We find that these probabilities are strongly dependent on the underlying cosmology and the efficiency of photon production, but also on the halo mass. There is a rapid transition between predominantly external and predominantly internal reionization at a mass scale of 1.0e12 Msun (corresponding roughly to L*galaxies), with haloes less massive than this being reionized preferentially from distant sources. We provide a fit for the reionization redshift as a function of halo mass, which could be helpful to parameterize reionization in semi-analytical models of galaxy formation on cosmological scales. We find no statistical correlation between the reionization history of field galaxies and their environment., Comment: 11 pages, 7 figures. Accepted for publication in MNRAS. Methodology section rewritten for clarity, results and conclusions unchanged. A preprint with high-resolution figures is available at http://www-theorie.physik.unizh.ch/~andrea/Reion/

We present non-parametric radial mass profiles for ten QSO strong lensing galaxies. Five of the galaxies have profiles close to $\rho(r)\propto r^{-2}$, while the rest are closer to r^{-1}, consistent with an NFW profile. The former are all relatively isolated early-types and dominated by their stellar light. The latter --though the modeling code did not know this-- are either in clusters, or have very high mass-to-light, suggesting dark-matter dominant lenses (one is a actually pair of merging galaxies). The same models give $H_0^{-1} = 15.2_{-1.7}^{+2.5}\Gyr$ ($H_0 = 64_{-9}^{+8} \legacy$), consistent with a previous determination. When tested on simulated lenses taken from a cosmological hydrodynamical simulation, our modeling pipeline recovers both H_0 and $\rho(r)$ within estimated uncertainties. Our result is contrary to some recent claims that lensing time delays imply either a low H_0 or galaxy profiles much steeper than r^{-2}. We diagnose these claims as resulting from an invalid modeling approximation: that small deviations from a power-law profile have a small effect on lensing time-delays. In fact, as we show using using both perturbation theory and numerical computation from a galaxy-formation simulation, a first-order perturbation of an isothermal lens can produce a zeroth-order change in the time delays., Comment: Replaced with final version accepted for publication in ApJ; very minor changes to text; high resolution figures may be obtained at justinread.net

It has been shown that the dissipative gas infall during galaxy formation has the capability to modify the shape of dark matter halos. In this paper we perform the first detailed analysis of particle orbits in a cosmological dark matter halos to understand {\it how} and {\it why} baryons alter its shape. We perform a series of numerical experiments where we grow a baryonic core inside a live dark matter halo extracted from a cosmological simulation. We follow the orbits of more than 300 particles with more that 50000 timesteps. Our results clearly show that the dissipational component is responsible for repeatedly deflecting orbits which visit often the center of the system. Moreover the gravitational potential time dependence associated with the growth of the baryonic mass, shifts the frequencies of the orbits, making them extremely chaotic. This randomization makes the orbits explore a large phase space. When this effect takes place for a significant number of orbits it will be manifested in the density distribution as an approach to a rounder configuration. As a consequence, the influence of the central mass on the shape of the phase space decreases with increasing distance from the center. We discuss the importance of future analysis of controlled experiments (i.e. using analytic potentials instead of live DM halos) to better decipher the dynamics of this phenomenon., Comment: 14 pages, 8 figures, MNRAS submitted, a preprint with high-resolution figures is available at http://www-theorie.physik.unizh.ch/~andrea/Orbits/

We use large cosmological N-body simulations to study the subhalo population in galaxy group sized halos. In particular, we look for fossil group candidates with typical masses ~10-25% of Virgo cluster but with an order of magnitude less substructure. We examine recent claims that the earliest systems to form are deficient enough in substructure to explain the luminosity function found in fossil groups. Although our simulations show a correlation between the halo formation time and the number of subhalos, the maximum suppression of subhalos is a factor of 2-2.5, whereas a factor of 6 is required to match fossil groups and galaxies. While the number of subhalos depends weakly on the formation time, the slope of the halo substructure velocity function does not. The satellite population within Cold Dark Matter (CDM) halos is self-similar at scales between galaxies and galaxy clusters regardless of mass, whereas current observations show a break in self-similarity at a mass scale corresponding to group of galaxies., Comment: Submitted for publication in MNRAS, 8 pages, 8 figures

We use the conditional luminosity function (CLF) and data from the 2dFGRS to constrain the average relation between light and mass in a LCDM cosmology with Omega_m=0.23 and sigma_8=0.74 (hereafter WMAP3 cosmology). Reproducing the observed luminosity dependence of the galaxy two-point correlation function results in average mass-to-light ratios that are about 35 percent lower than in a LCDM cosmology with Omega_m=0.3 and sigma_8=0.9 (hereafter WMAP1 cosmology). This removes an important problem with previous halo occupation models which had a tendency to predict cluster mass-to-light ratios that were too high. For the WMAP3 cosmology our model yields average mass-to-light ratios, central galaxy luminosities, halo occupation numbers, satellite fractions, and luminosity-gap statistics, that are all in excellent agreement with those obtained from a 2dFGRS group catalogue and from other independent studies. We also use our CLF model to compute the probability distribution P(M|L_cen), that a central galaxy of luminosity L_cen resides in a halo of mass M. We find this distribution to be much broader than what is typically assumed in HOD models, which has important implications for the interpretation of galaxy-galaxy lensing data. Finally, reproducing the luminosity dependence of the pairwise velocity dispersions in the 2dFGRS requires relatively low mass-to-light ratios for clusters and a satellite fraction that decreases strongly with increasing luminosity. This is only marginally consistent with our CLF constraints. We argue that a cosmology with parameters between those of the WMAP1 and WMAP3 cosmologies is likely to yield results with a higher level of consistency., Comment: 23 pages, 11 figures. Submitted for publication in MNRAS

A possible explanation of the present discrepancy between the abundance of galactic subhaloes predicted by N-Body simulations with those observed in the Local Group is presented. We study the impact of an early reionization on the baryonic component of the Universe using SPH simulations on group and galactic scales. We use a simplified model for reionization described as an instantaneous increment of the IGM temperature (with 1e5 8). We find that a key role is played by compton cooling (interaction between hot electrons and the CMB photons); at high redshift (z_r>10) this cooling is very efficient and it is able to counteract any heating of the gas within few Myrs. This means that a late reionization is needed (z_r<9) to sufficiently reduce the number of luminous dwarf satellites around our Galaxy. For a reionization redshift z_r=8 and a reionization temperature of T_{IGM}~1e5 K we are able to reproduce the observed number of Local Group dwarf galaxies in our simulations., Comment: 6 pages, 8 figures. To appear in the proceedings of the XLIst Rencontres de Moriond, XXVIth Astrophysics Moriond Meeting: "From dark halos to light", Eds. L.Tresse, S. Maurogordato and J. Tran Thanh Van (Editions Frontieres)

We use a series of cosmological N-body simulations for a flat LCDM cosmology to investigate the properties of dark matter haloes in the mass range 3.0e9-3.0e13 Msun. These properties include the concentration parameter (c), the spin parameter (lambda) and the mean axis ratio (q). For the concentration-mass relation we find c~M^(-0.11) in agreement with the model proposed by Bullock et al. even if we find a lower normalization (15%). The results for lambda and q are in good agreement with previous studies, while c and lambda are anti-correlated. In an attempt to remove unrelaxed haloes, we use the offset parameter (xoff), defined as the distance between the most bound particle and the center of mass. Removing haloes with large xoff increases the c by ~10%, lowers the lambda by ~15%, and removes the most prolate haloes. In addition, it largely removes the anti-correlation between c and lambda though not entirely. We also investigate the effects of the large-scale environment. We find that more concentrated haloes live in denser environments. Note, however, that the trend is weak compared to the scatter. For the spin parameters we find no environment dependence, while there is a weak indication that the most spherical haloes reside in denser region. Finally, using a simple model for disk galaxy formation we show that haloes that host low surface brightness galaxies are expected to be hosted by a biased sub-set of haloes. Not only do these haloes have spin parameters that are larger than average, they also have c that are 15% percent lower than the average at a given halo mass. We discuss the implications of all these findings for the claimed disagreement between halo concentrations inferred from LSB rotation curves, and those expected for a LCDM cosmology. (abridged), Comment: 23 pages, 14 figure. Resolution and environment effects discussed in more details. Conclusions unchanged. References added. Accepted for publication by MNRAS. A preprint with high-resolution figures is available at http://www-theorie.physik.unizh.ch/~andrea/ConcParam/

We present a simultaneous analysis of 10 galaxy lenses having time-delay measurements. For each lens we derive a detailed free-form mass map, with uncertainties, and with the additional requirement of a shared value of the Hubble parameter across all the lenses. We test the prior involved in the lens reconstruction against a galaxy-formation simulation. Assuming a concordance cosmology, we obtain 1/H_0 = 13.5 (+2.5/-1.3) Gyr, Comment: To appear in ApJ Lett

It has been argued that the flux anomalies detected in gravitationally lensed QSOs are evidence for substructures in the foreground lensing haloes. In this paper we investigate this issue in greater detail focusing on the Cusp relation which corresponds to images of a source located to the cusp of the inner caustic curve. We use numerical simulations combined with a Monte Carlo approach to study the effects of the expected power law distribution of substructures within LCDM haloes on the multiple images. Generally, the high number of anomalous flux ratios in the cusp configurations is unlikely explained by 'simple' perturbers (subhaloes) inside the lensing galaxy, either modeled by point masses or extended NFW subhaloes. We considered in our analysis a mass range of 10^5-10^7 Msun for the subhaloes. We also demonstrate that including the effects of the surrounding mass distribution, such as other galaxies close to the primary lens, does not change the results. We conclude that triple images of lensed QSOs do not show any direct evidence for dark dwarf galaxies such as cold dark matter substructure., Comment: 10 pages, 19 figures, Effects of different subhalos concentrations discussed, analysis improved, accepted by MNRAS

Polar ring galaxies are flattened stellar systems with an extended ring of gas and stars rotating in a plane almost perpendicular to the central galaxy. We show that their formation can occur naturally in a hierarchical universe where most low mass galaxies are assembled through the accretion of cold gas infalling along megaparsec scale filamentary structures. Within a large cosmological hydrodynamical simulation we find a system that closely resembles the classic polar ring galaxy NGC 4650A. How galaxies acquire their gas is a major uncertainty in models of galaxy formation and recent theoretical work has argued that cold accretion plays a major role. This idea is supported by our numerical simulations and the fact that polar ring galaxies are typically low mass systems., Comment: 4 pages, 5 figures, stability of the ring discussed, minor changes to match the accepted version by ApJL. A preprint with high-resolution figures is available at http://krone.physik.unizh.ch/~andrea/PolarRing/PolarRing.ps

We analyse the number density and radial distribution of substructures and satellite galaxies using cosmological simulations that follow the gas dynamics of a baryonic component, including shock heating, radiative cooling and star formation within the hierarchical concordance LCDM model. We find that the dissipation of the baryons greatly enhances the survival of subhaloes, expecially in the galaxy core, resulting in a radial distribution of satellite galaxies that closely follows the overall mass distribution in the inner part of the halo. Hydrodynamical simulations are necessary to resolve the adiabatic contraction and dense cores of galaxies, resulting in a total number of satellites a factor of two larger than found in pure dark matter simulation, in good agreement with the observed spatial distribution of satellite galaxies within galaxies and clusters. Convergence tests show that the cored distribution found by previous authors in pure N-body simulations was due to physical overmerging of dark matter only structures. We proceed to use a ray-shooting technique in order to study the impact of these additional substructures on the number of violations of the cusp caustic magnification relation. We develop a new approach to try to disentangle the effect of substructures from the intrinsic discreteness of N-Body simulations. Even with the increased number of substructures in the centres of galaxies, we are not able to reproduce the observed high numbers of discrepancies observed in the flux ratios of multiply lensed quasars., Comment: 11 pages, 15 figures, comparison with previous works updated, one more plot added, minor changes to match the accepted version by MNRAS

Dynamical Dark Energy (DE) is a viable alternative to the cosmological constant. Yet, constructing tests to discriminate between Lambda and dynamical DE models is difficult because the differences are not large. In this paper we explore tests based on the galaxy mass function, the void probability function (VPF), and the number of galaxy clusters. At high z the number density of clusters shows large differences between DE models, but geometrical factors reduce the differences substantially. We find that detecting a model dependence in the cluster redshift distribution is a hard challenge. We show that the galaxy redshift distribution is potentially a more sensitive characteristics. We do so by populating dark matter halos in Nbody simulations with galaxies using well-tested Halo Occupation Distribution (HOD). We also estimate the Void Probability Function and find that, in samples with the same angular surface density of galaxies in different models, the VPF is almost model independent and cannot be used as a test for DE. Once again, geometry and cosmic evolution compensate each other. By comparing VPF's for samples with fixed galaxy mass limits, we find measurable differences., Comment: 12 pages, 11 figures, dependence on mass-luminosity relation discussed, minor changes to match the accepted version by MNRAS

Using N-body simulations of flat, dark energy dominated cosmologies, we show that galaxies around simulated binary systems resembling the Local Group (LG) have low peculiar velocities, in good agreement with observational data. We have compared results for LG-like systems selected from large, high resolution simulations of three cosmologies: a LCDM model, a LWDM model with a 2 keV warm dark matter candidate and a quintessence model (QCDM) with an equation of state parameter w=-0.6. The Hubble flow is significant colder around LGs selected in a flat, Lambda dominated cosmology than around LGs in open or critical models, showing that a dark energy component manifests itself on the scales of nearby galaxies, cooling galaxy peculiar motions. Flows in the LWDM and QCDM models are marginally colder than in the LCDM one. The results of our simulations have been compared to existing data and a new data set of 28 nearby galaxies with robust distance measures (Cepheids and Surface Brightness Fluctuations). The measured line-of-sight velocity dispersion is sigma = 88 +- 20 km/sec x (R/7 Mpc). The best agreement with observations is found for LGs selected in the $\Lambda$CDM cosmology in environments with -0.1

We study the strong gravitational lensing properties of galaxy clusters obtained from N-body simulations with different kind of Dark Energy (DE). We consider both dynamical DE, due to a scalar field self--interacting through Ratra-Peebles (RP) or SUGRA potentials, and DE with constant negative w=p/rho= -1 (LCDM). We have 12 high resolution lensing systems for each cosmological model with a mass greater than 5x10^{14} solar masses. Using a Ray Shooting technique we make a detailed analysis of the lensing properties of these clusters with particular attention to the number of arcs and their properties (magnification, length and width). We found that the number of giant arcs produced by galaxy clusters changes in a considerable way from LCDM models to Dynamical Dark Energy models with a RP or SUGRA potentials. These differences originate from the different epochs of cluster formation and from the non-linearity of the strong lensing effect. We suggest the Strong lensing is one of the best tool to discriminate among different kind of Dark Energy., Comment: 8 pages, 11 figures. Lensing map resolution improved and effects of resolution discussed. One more RP model analysed. Accepted for publication by MNRAS

We present the results of a series of N-body simulations in cosmologies where dark matter (DM) is coupled to dark energy (DE), so easing the cosmic coincidence problem. The dark-dark coupling introduces two novel effects in N-body dynamics: (i) DM particle masses vary with time; (ii) gravity between DM particles is ruled by a constant $G^{*}$, greater than Newton's constant $G$, holding in other 2-body interactions. As a consequence, baryons and DM particle distributions develop a large scale bias. Here we investigate DE models with Ratra-Peebles (RP) potentials; the dark-dark coupling is set in a parametric range compatible with observations, for as concern background and linear perturbation properties. We study the halo mass function, the halo density profile and the behavior of the non-linear bias. We find that non-linear dynamics puts additional constraints to the coupling parameter. They mostly arise from density profiles, that we find to yield higher concentrations, in coupled RP models, with respect to (uncoupled) dynamical DE cosmologies. Such enhancement, although being a strong effect in some coupling parameter range, is just a minor change for smaller but significant values of the coupling parameter. With these further restrictions, coupled DE models with RP potential are consistent with non-linear observables., Comment: 12 pages, 10 figures, v2, added some references, submitted to Phys.Rev.D

We perform N-body simulations for models with a DE component. Besides of DE with constant negative state parameter w, we consider DE due to scalar fields, self-interacting through RP or SUGRA potentials. According to our post-linear analysis, at z=0, DM power spectra and halo mass functions do not depend on DE nature. This is welcome, as LCDM fits observations. Halo profiles, instead, are denser than LCDM. For example, the density at 15 kpc of a DE halo with M=10^13 exceeds LCDM by ~45%. Differences, therefore, are small but, however, DE does not ease the problem with cuspy DM profiles. We study also subhalos and find that, at $z=0$, the number of satellites coincides in all DE models. At higher z, DE models show increasing differences from LCDM and among themselves (i.e. in the mass function evolution); this is the obvious pattern to distinguish between different DE state equations., Comment: 4 Pages, 7 Figures, to appear in the proceedings of the congress "Multiwavelenght Cosmology", Mykonos June 17-20

Cosmological models with different types of Dark Energy are becoming viable alternatives for standard models with the cosmological constant. Yet, such models are more difficult to analyze and to simulate. We present analytical approximations and discuss ways of making simulations for two families of models, which cover a wide range of possibilities and include models with both slow and fast changing ratio w=p\rho. More specifically, we give analytical expressions for the evolution of the matter density parameter Omega_m(z) and the virial density contrast Delta_c at any redshift z. The latter is used to identify halos and to find their virial masses. We also provide an approximation for the linear growth factor of linear fluctuations between redshift z=40 and z=0. This is needed to set the normalization of the spectrum of fluctuations. Finally, we discuss the expected behavior of the halo mass function and its time evolution., Comment: 10 pages, 10 figures ApJ submitted

We study properties of dark matter halos in a variety of models which include Dark Energy (DE). We consider both DE due to a scalar field self--interacting through Ratra-Peebles or SUGRA potentials, and DE with constant negative w=p\rho >-1. We find that at redshift zero the nonlinear power spectrum of the dark matter, and the mass function of halos, practically do not depend on DE state equation and are almost indistinguishable from predictions of the LCDM model. This is consistent with the nonlinear analysis presented in the accompanying paper. It is also a welcome feature because LCDM models fit a large variety of data. On the other hand, at high redshifts DE models show substantial differences from LCDM and substantial differences among themselves. Halo profiles differ even at z=0. DE halos are denser than LCDM in their central parts because the DE halos collapse earlier. Nevertheless, differences between the models are not so large. For example, the density at 10 kpc of a DE \~10^{13}Msun halo deviates from LCDM by not more than 50%. This, however, means that DE is not a way to ease the problem with cuspy dark matter profiles. Addressing another cosmological problem - abundance of subhalos -- we find that the number of satellites of halos in various DE models does not change relative to the LCDM, when normalized to the same circular velocity of the parent halo. To summarize, the best way to find which DE model fits the observed Universe is to look for evolution of halo properties. For example, the abundance of galaxy groups with mass larger than 10^{13}\Msun at z> 2 can be used to discriminate between the models, and, thus, to constrain the nature of DE., Comment: 10 pages, 8 figures ApJ Submitted

We show that DM halos, in n--body simulations, have a boundary layer (BL), separating bound from unbound mass. Let T(r) and W(r) be the kinetic and potential energies in shells of halos. We find that, in almost all halos: (i) The virial ratio R=-2T/W has at least one persistent (i.e. resolution independent) minimum r_c, such that, close to it, (ii) the function w=-d log W/d log r has a maximum and (iii) the relation R(r_c)=w(r_c) is almost exactly fulfilled. The radius r_c is the position of BL's, in halos found in simulations of TCDM and LCDM models, run using ART and GADGET codes at various resolutions. We find that 97% of the ~300 halos (per model) with M>4.2 10^14 M_s h^-1 owns a BL. Those with no BL are undergoing major mergings. The mass M_c enclosed in r_c almost coincides with the mass evaluated from velocities, according to virial theorem. Particles at r>r_c are not in virial equilibrium. Using r_c we have a density contrast \Delta_c for each halo. For each mass scale, \Delta_v=178 Omega_m^0.45 is within the range of \Delta_c's, but the spread in \Delta_c is wide and the average \Delta_c is ~25% smaller than the corresponding \Delta_v. The matching of properties derived under the assumption of spherical symmetry is a consequence of violent relaxation destroying features related to previous ellipsoidal non-linear growth. In turn, the spread of final \Delta_c's is an imprint of the different 3-D geometries and of the variable environment during collapses, as suggested by a comparison with Sheth & Tormen analysis., Comment: 19 Pages, 15 figures, ApJ in press

We study the cluster mass function and its evolution in different models with Dark Energy arising from a self--interacting scalar field, with Ratra-Peebles and SUGRA potentials. Computations are based on a Press & Schechter approximation. The mass functions we obtain are compared with results holding for open models or models with Dark Energy due to a cosmological constant. Evolution results, in the Dark Energy cases, closely approach open models., Comment: 13 pages, 3 new figures included, references added. Accepted for pubblication in New Astronomy

We show that dark matter haloes, in n--body simulations, have a boundary layer (BL) with precise features. In particular, it encloses all dynamically stable mass while, outside it, dynamical stability is lost soon. Particles can pass through such BL, which however acts as a confinement barrier for dynamical properties. BL is set by evaluating kinetic and potential energies (T(r) and W(r)) and calculating R=-2T/W. Then, on BL, R has a minimum which closely approaches a maximum of w= -dlog W/dlog r. Such $Rw$ ``requirement'' is consistent with virial equilibrium, but implies further regularities. We test the presence of a BL around haloes in spatially flat CDM simulations, with or without cosmological constant. We find that the mass M_c, enclosed within the radius r_c, where the $Rw$ requirement is fulfilled, closely approaches the mass M_{dyn}, evaluated from the velocities of all particles within r_c, according to the virial theorem. Using r_c we can then determine an individual density contrast Delta_c for each virialized halo, which can be compared with the "virial" density contrast $\Delta_v ~178 \Omega_m^{0.45}$ (Omega_m: matter density parameter) obtained assuming a spherically symmetric and unperturbed fluctuation growth. The spread in Delta_c is wide, and cannot be neglected when global physical quantities related to the clusters are calculated, while the average Delta_c is ~25 % smaller than the corresponding Delta_v; moreover if $M_{dyn}$ is defined from the radius linked to Delta_v, we have a much worse fit with particle mass then starting from {\it Rw} requirement., Comment: 4 pages, 5 figures, contribution to the XXXVIIth Rencontres de Moriond, The Cosmological Model, Les Arc March 16-23 2002, to appear in the proceedings

We examine how the statistics of the quadrupoles of (projected) cluster masses can discriminate between flat cold dark matter (CDM) universes with or without a cosmological constant term. Even in the era of high precision cosmology that cosmic microwave background experiments should open soon, it is important to devise self consistency tests of cosmogonic theories tuned at the matter radiation decoupling epoch using data from the non--linear evolved universe. We build cluster catalogs from two large volume simulations of a ``tilted'' CDM model and a $\Lambda$CDM model with cosmic density parameter $\Omega_m=0.35$ and cosmological constant contribution $\Omega_\Lambda=0.65$. From the projected mass distribution of the clusters we work out the quadrupoles Q and examine their dependence on cluster mass and the cosmological model. We find that TCDM clusters have systematically larger quadrupoles than their \lcdm counterpart. The effect is mass dependent: massive clusters ($M\gsim 10^{15}\Msunh$) have quadrupoles differing by more than 30% in the two models, while for $M\lsim 4\times10^{14}\Msunh$ the difference rapidly drops to $\sim 1$%. Performing a K-S test of the Q distributions, we estimate that using just the 15 most massive clusters in the simulation volume ($360\Mpch$ a side) we can discriminate between TCDM and \lcdm at a confidence level better than 99.9%. In the volume probed by exhisting observations, there are potentially several hundred clusters with masses above the threshold for which the differences in the quadrupoles become relevant. Should weak lensing data become available for this whole set, a quadrupole analysis may be expected to discriminate among different values of $\Lambda$., Comment: 8 pages, 4 figures. Accepted for pubblication in ApJ

The passage from angular to spatial correlations, in the case of spatial clustering length depending on the average distance between nearby objects is studied. We show that, in a number of cases, the scaling law of angular correlation amplitudes, which holds for constant spatial clustering length, is still true also for a luminosity dependent spatial correlation. If the Limber equation is then naively used to obtain `the' spatial clustering length from the angular function amplitude, a quantity close to the average object separation is obtained. The case of cluster clustering is explicitly considered., Comment: 10 pages, LaTeX, no figures, Accepted for publication in New Astronomy