The distribution of galaxy morphological types is a key test for models of galaxy formation and evolution, providing strong constraints on the relative contribution of different physical processes responsible for the growth of the spheroidal components. In this paper, we make use of a suite of semi-analytic models to study the efficiency of galaxy mergers in disrupting galaxy discs and building galaxy bulges. In particular, we compare standard prescriptions usually adopted in semi-analytic models, with new prescriptions proposed by Kannan et al., based on results from high-resolution hydrodynamical simulations, and we show that these new implementations reduce the efficiency of bulge formation through mergers. In addition, we compare our model results with a variety of observational measurements of the fraction of spheroid-dominated galaxies as a function of stellar and halo mass, showing that the present uncertainties in the data represent an important limitation to our understanding of spheroid formation. Our results indicate that the main tension between theoretical models and observations does not stem from the survival of purely disc structures (i.e. bulgeless galaxies), rather from the distribution of galaxies of different morphological types, as a function of their stellar mass.

We study the effect of mergers on the morphology of galaxies by means of the simulated merger tree approach first proposed by Moster et al. This method combines N-body cosmological simulations and semi-analytic techniques to extract realistic initial conditions for galaxy mergers. These are then evolved using high resolution hydrodynamical simulations, which include dark matter, stars, cold gas in the disc and hot gas in the halo. We show that the satellite mass accretion is not as effective as previously thought, as there is substantial stellar stripping before the final merger. The fraction of stellar disc mass transferred to the bulge is quite low, even in the case of a major merger, mainly due to the dispersion of part of the stellar disc mass into the halo. We confirm the findings of Hopkins et al., that a gas rich disc is able to survive major mergers more efficiently. The enhanced star formation associated with the merger is not localised to the bulge of galaxy, but a substantial fraction takes place in the disc too. The inclusion of the hot gas reservoir in the galaxy model contributes to reducing the efficiency of bulge formation. Overall, our findings suggest that mergers are not as efficient as previously thought in transforming discs into bulges. This possibly alleviates some of the tensions between observations of bulgeless galaxies and the hierarchical scenario for structure formation.

Models including an energy transfer from CDM to DE are widely considered in the literature, namely to allow DE a significant high-z density. Strongly Coupled cosmologies assume a much larger coupling between DE and CDM, together with the presence of an uncoupled warm DM component, as the role of CDM is mostly restricted to radiative eras. This allows us to preserve small scale fluctuations even if the warm particle, possibly a sterile neutrino, is quite light, O(100 eV). Linear theory and numerical simulations show that these cosmologies agree with LCDM on supergalactic scales; e.g., CMB spectra are substantially identical. Simultaneously, simulations show that they significantly ease problems related to the properties of MW satellites and cores in dwarfs. SC cosmologies also open new perspectives on early black hole formation, and possibly lead towards unificating DE and inflationary scalar fields.

The annihilation of dark matter particles in the halo of galaxies may end up into gamma rays, which travel almost unperturbed till to their detection at Earth. This annihilation signal can exhibit an anisotropic behavior quantified by the angular power spectrum, whose properties strongly depend on the dark matter distribution and its clumpiness. We use high resolution pure dark matter N-body simulations to quantify the contribution of different components (main halo and satellites) to the global signal as a function of the analytical profile adopted to describe the numerical results. We find that the smooth main halo dominates the angular power spectrum of the gamma-ray signal up to quite large multipoles, where the sub-haloes anisotropy signal starts to emerge, but the transition multipole strongly depends on the assumed radial profile. The extrapolation down to radii not resolved by current numerical simulations can affect both the normalization and the shape of the gamma-ray angular power spectrum. For the sub-haloes described by an asymptotically cored dark matter distribution, the angular power spectrum shows an overall smaller normalization and a flattening at high multipoles. Our results show the criticality of the dark matter density profile shape in gamma-ray anisotropy searches, and evaluate quantitatively the intrinsic errors occurring when extrapolating the dark matter radial profiles down to spatial scales not yet explored by numerical simulations.

We use the relations between aperture stellar velocity dispersion (\sigma_ap), stellar mass (M_sps), and galaxy size (R_e) for a sample of \sim 150,000 early-type galaxies from SDSS/DR7 to place constraints on the stellar initial mass function (IMF) and dark halo response to galaxy formation. We build LCDM based mass models that reproduce, by construction, the relations between galaxy size, light concentration and stellar mass, and use the spherical Jeans equations to predict \sigma_ap. Given our model assumptions (including those in the stellar population synthesis models), we find that reproducing the median \sigma_ap vs M_sps relation is not possible with {\it both} a universal IMF and a universal dark halo response. Significant departures from a universal IMF and/or dark halo response are required, but there is a degeneracy between these two solutions. We show that this degeneracy can be broken using the strength of the correlation between residuals of the velocity-mass (\Delta log \sigma_ap) and size-mass (\Delta log R_e) relations. The slope of this correlation, d_vr \equiv \Delta log \sigma_ap/\Delta log R_e, varies systematically with galaxy mass from d_vr \simeq -0.45 at M_sps \sim 10^{10}M_sun, to d_vr \simeq -0.15 at M_sps \sim 10^{11.6} M_sun. The virial fundamental plane (FP) has d_vr=-1/2, and thus we find the tilt of the observed FP is mass dependent. Reproducing this tilt requires {\it both} a non-universal IMF and a non-universal halo response. Our best model has mass-follows-light at low masses (Msps < 10^{11.2}M_sun) and unmodified NFW haloes at M_sps \sim 10^{11.5} M_sun. The stellar masses imply a mass dependent IMF which is "lighter" than Salpeter at low masses and "heavier" than Salpeter at high masses.

In this paper we study the tidal stripping process for satellite galaxies orbiting around a massive host galaxy, and focus on its dependence on the morphology of both satellite and host galaxy. For this purpose, we use three different morphologies for the satellites: pure disc, pure bulge and a mixture bulge+disc. Two morphologies are used for the host galaxies: bulge+disc and pure bulge. We find that while the spheroidal stellar component experiences a constant power-law like mass removal, the disc is exposed to an exponential mass loss when the tidal radius of the satellite is of the same order of the disc scale length. This dramatic mass loss is able to completely remove the stellar component on time scale of 100 Myears. As a consequence two satellites with the same stellar and dark matter masses, on the same orbit could either retain considerable fraction of their stellar mass after 10 Gyrs or being completely destroyed, depending on their initial stellar morphology. We find that there are two characteristic time scales describing the beginning and the end of the disc removal, whose values are related to the size of the disc. This result can be easily incorporated in semi-analytical models. We also find that the host morphology and the orbital parameters also have an effect on the determining the mass removal, but they are of secondary importance with respect to satellite morphology. We conclude that satellite morphology has a very strong effect on the efficiency of stellar stripping and should be taken into account in modeling galaxy formation and evolution.

We analyze the properties of dark matter halos in the cold-plus-warm dark matter cosmologies (CWDM). We study their dependence on the fraction and velocity dispersion of the warm particle, keeping the free-streaming scale fixed. To this end we consider three models with the same free-streaming: (1) a mixture of 90% of CDM and 10% of WDM with the mass 1 keV; (2) a mixture of 50% of CDM and 50% of WDM with the mass 5 keV; and (3) pure WDM with the mass 10 keV. "Warm" particles have rescaled Fermi-Dirac spectrum of primordial velocities (as non-resonantly produced sterile neutrinos would have). We compare the properties of halos among these models and with a LCDM with the same cosmological parameters. We demonstrate, that although these models have the same free-streaming length and the suppression of matter spectra are similar at scales probed by the Lyman-alpha forest comoving wave-numbers k<3-5 h/Mpc), the resulting properties of halos with masses below 1e11 Msun are different due to the different behaviour of matter power spectra at smaller scales. In particular, we find that while the number of galaxies remains the same as in LCDM case, their density profiles become much less concentrated, and hence in better agreement with current observational constraints. Our results imply that a single parameter (e.g. free streaming length) description of these models is not enough to fully capture their effects on the structure formation process.

The free streaming of warm dark matter particles dampens the fluctuation spectrum, flattens the mass function of haloes and imprints a fine grained phase density limit for dark matter structures. The phase space density limit is expected to imprint a constant density core at the halo center on the contrary to what happens for cold dark matter. We explore these effects using high resolution simulations of structure formation in different warm dark matter scenarios. We find that the size of the core we obtain in simulated haloes is in good agreement with theoretical expectations based on Liouville's theorem. However, our simulations show that in order to create a significant core, (r_c~1 kpc), in a dwarf galaxy (M~1e10 Msun), a thermal candidate with a mass as low as 0.1 keV is required. This would fully prevent the formation of the dwarf galaxy in the first place. For candidates satisfying large scale structure constrains (m_wdm larger than 1-2 keV) the expected size of the core is of the order of 10 (20) pc for a dark matter halo with a mass of 1e10 (1e8) Msun. We conclude that "standard" warm dark matter is not viable solution for explaining the presence of cored density profiles in low mass galaxies.

A clear prediction of the Cold Dark Matter model is the existence of cuspy dark matter halo density profiles on all mass scales. This is not in agreement with the observed rotation curves of spiral galaxies, challenging on small scales the otherwise successful CDM paradigm. In this work we employ high resolution cosmological hydro-dynamical simulations to study the effects of dissipative processes on the inner distribution of dark matter in Milky-Way like objects (M~1e12 Msun). Our simulations include supernova feedback, and the effects of the radiation pressure of massive stars before they explode as supernovae. The increased stellar feedback results in the expansion of the dark matter halo instead of contraction with respect to N-body simulations. Baryons are able to erase the dark matter cuspy distribution creating a flat, cored, dark matter density profile in the central several kpc of a massive Milky-Way like halo. The profile is well fit by a Burkert profile, with fitting parameters consistent with the observations. In addition, we obtain flat rotation curves as well as extended, exponential stellar disk profiles. While the stellar disk we obtain is still partially too thick to resemble the MW thin disk, this pilot study shows that there is enough energy available in the baryonic component to alter the dark matter distribution even in massive disc galaxies, providing a possible solution to the long standing problem of cusps vs. cores.

We employ hydrodynamical simulations to study the effects of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. For the first time we present a suite of simulations that includes a diffuse, rotating, cooling, hot gaseous halo, as predicted by cosmological hydrodynamical simulations as well as models of galaxy formation. We study the effect of this new gaseous component on the vertical structure of a Milky Way-like stellar disc during 1:10 and 1:5 mergers. For 1:10 mergers we find no increased final thin disc scale height compared to the isolated simulation, leading to the conclusion that thin discs can be present even after a 1:10 merger if a reasonable amount of hot gas is present. The reason for this is the accretion of new cold gas, leading to the formation of a massive new thin stellar disc that dominates the surface brightness profile. In a previous study, in which we included only cold gas in the disk, we showed that the presence of cold gas decreased the thickening by a minor merger relative to the no-gas case. Here, we show that the evolution of the scale height in the presences of a cooling hot halo is dominated by the formation of the new stellar disc. In this scenario, the thick disc is the old stellar disc that has been thickened in a minor merger at z>1, while the thin disc is the new stellar disc that reforms after this merger. In addition, we study the evolution of the scale height during a 1:5 merger and find that a thin disc can be present even after this merger, provided enough hot gas is available. The final scale height in our simulations depends on the mass of the hot gaseous halo, the efficiency of the winds and the merger mass ratio. We find post-merger values in the range 0.5<z0<1.0 kpc in good agreement with observational constraints by local galaxies.

We study the correlations between the structural parameters of dark matter haloes using Principal Component Analysis (PCA). We consider a set of eight parameters, six of which are commonly used to characterize dark matter halo properties: mass, concentration, spin, shape, overdensity, and the angle ($\Phi_L$) between the major axis and the angular momentum vector. Two additional parameters ($\x_{off}$ and $\rho_{rms}$) are used to describe the degree of `relaxedness' of the halo. We find that we can account for much of the variance of these properties with halo mass and concentration, on the one hand, and halo relaxedness on the other. Nonetheless, three principle components are usually required to account for most of the variance. We argue that halo mass is not as dominant as expected, which is a challenge for halo occupation models and semi-analytic models that assume that mass determines other halo (and galaxy) properties. In addition, we find that the angle $\Phi_L$ is not significantly correlated with other halo parameters, which may present a difficulty for models in which galaxy disks are oriented in haloes in a particular way. Finally, at fixed mass, we find that a halo's environment (quantified by the large-scale overdensity) is relatively unimportant.

Cosmological hydrodynamical simulations as well as observations indicate that spiral galaxies are comprised of five different components: dark matter halo, stellar disc, stellar bulge, gaseous disc and gaseous halo. While the first four components have been extensively considered in numerical simulations of binary galaxy mergers, the effect of a hot gaseous halo has usually been neglected even though it can contain up to 80% of the total gas within the galaxy virial radius. We present a series of hydrodynamic simulations of major mergers of disc galaxies, that for the first time include a diffuse, rotating, hot gaseous halo. Through cooling and accretion, the hot halo can dissipate and refuel the cold gas disc before and after a merger. This cold gas can subsequently form stars, thus impacting the morphology and kinematics of the remnant. Simulations of isolated systems with total mass M~10^12Msun show a nearly constant star formation rate of ~5Msun/yr if the hot gaseous halo is included, while the star formation rate declines exponentially if it is neglected. We conduct a detailed study of the star formation efficiency during mergers and find that the presence of a hot gaseous halo reduces the starburst efficiency (e=0.5) compared to simulations without a hot halo (e=0.68). Moreover we find cases where the stellar mass of the merger remnant is lower than the sum of the stellar mass of the two progenitor galaxies when evolved in isolation. This suggests a revision to semi-analytic galaxy formation models which assume that a merger always leads to enhanced star formation. We show that adding the hot gas component has a significant effect on the kinematics and internal structure of the merger remnants, like an increased abundance of fast rotators and an r^(1/4) surface brightness profile at small scales.

In this work, we examine the different properties of galactic satellites in hydrodynamical and pure dark matter simulations. We use three pairs of simulations (collisional and collision-less) starting from identical initial conditions. We concentrate our analysis on pairs of satellites in the hydro and Nbody runs that form from the same Lagrangian region. We look at the radial positions, mass loss as a function of time and orbital parameters of these "twin" satellites. We confirm an overall higher radial density of satellites in the hydrodynamical runs, but find that trends in the mass loss and radial position of these satellites in the inner and outer region of the parent halo differ from the pure dark matter case. In the outskirts of the halo (~70% of the virial radius) satellites experience a stronger mass loss and higher dynamical friction in pure dark matter runs. The situation is reversed in the central region of the halo, where hydrodynamical satellites have smaller apocenter distances and suffer higher mass stripping. We partially ascribe this bimodal behaviour to the delayed infall time for hydro satellites, which on average cross the virial radius of the parent halo 0.7 Gyrs after their dark matter twins. Finally, we briefly discuss the implications of the different set of satellite orbital parameters and mass loss rates in hydrodynamical simulations within the context of thin discs heating and destruction.

The next generation mass probes will obtain information on non--linear power spectra P(k,z) and their evolution, allowing us to investigate the nature of Dark Energy. To exploit such data we need high precision simulations, extending at least up to scales of k 10 h/Mpc, where the effects of baryons can no longer be neglected. In this paper, we present a series of large scale hydrodynamical simulations for LCDM and dynamical Dark Energy (dDE) models, in which the equation of state parameter is z-dependent. The simulations include gas cooling, star formation and Supernovae feedback. They closely approximate the observed star formation rate and the observationally derived star/Dark Matter mass ratio in collapsed systems. Baryon dynamics cause spectral shifts exceeding 1% at k > 2-3 h/Mpc compared to pure n-body simulations in the LCDM simulations. This agrees with previous studies, although we find a smaller effect (~50%) on the power spectrum amplitude at higher k's. dDE exhibits similar behavior, even though the dDE simulations produce ~20% less stars than the analogous LCDM cosmologies. Finally, we show that the technique introduced in Casarini et al. to obtain spectra for any $w(z)$ cosmology from constant-w models at any redshift still holds when gas physics is taken into account. While this relieves the need to explore the entire functional space of dark energy state equations, we illustrate a severe risk that future data analysis could lead to misinterpretation of the DE state equation.

The next generation mass probes will investigate DE nature by measuring non-linear power spectra at various z, and comparing them with high precision simulations. Producing a complete set of them, taking into account baryon physics and for any DE state equation w(z), would really be numerically expensive. Regularities reducing such duty are essential. This paper presents further n-body tests of a relation we found, linking models with DE state parameter w(z) to const.-w models, and also tests the relation in hydro simulations.

(Abridged) We apply a very general statistical theorem introduced by Cramer (1936) to study the origin of the deviations of the halo spin PDF from the reference lognormal shape. We find that these deviations originate from correlations between two quantities entering the definition of spin, namely the ratio $J/M^{5/2}$ (which depends only on mass) and the total gravitational binding energy $E$. To reach this conclusion, we have made usage of the results deduced from two high spatial- and mass resolution simulations. Our simulations cover a relatively small volume and produce a sample of more than 16.000 gravitationally bound halos, each traced by at least 300 particles. We verify that our results are stable to different systematics, by comparing our results with those derived by the GIF2 and by a more recent simulation performed by Maccio' et al. We find that the spin probability distribution function shows systematic deviations from a lognormal, at all redshifts z <= 1. These deviations depend on mass and redshift: at small masses they change little with redshift, and also the best lognormal fits are more stable. The J-M relationship is well described by a power law of exponent $\alpha$ very near to the linear theory prediction (alpha=5/3), but systematically lower than this at z<= 0.3. We argue that the fact that deviations from a lognormal PDF are present only for high-spin halos could point to a role of large-scale tidal fields in the evolution of the spin PDF.

We test the luminosity function of Milky Way satellites as a constraint for the nature of Dark Matter particles. We perform dissipationless high-resolution N-body simulations of the evolution of Galaxy-sized halo in the standard Cold Dark Matter (CDM) model and in four Warm Dark Matter (WDM) scenarios, with a different choice for the WDM particle mass (m_w). We then combine the results of the numerical simulations with semi-analytic models for galaxy formation, to infer the properties of the satellite population. Quite surprisingly we find that even WDM models with relatively low m_w values (2-5 keV) are able to reproduce the observed abundance of ultra faint (Mv<-9) dwarf galaxies, as well as the observed relation between Luminosity and mass within 300 pc. Our results suggest a lower limit of 1 keV for thermal warm dark matter, in broad agreement with previous results from other astrophysical observations like Lyman-alpha forest and gravitational lensing.

We lay out the framework to numerically study nonlinear structure formation in the context of scalar-field-coupled cold dark matter models (phiCDM models) where the scalar field phi serves as dynamical dark energy. Adopting parameters for the scalar field which yield a realistic CMB spectrum, we generate the initial conditions for our Nbody simulations, which follow the spatial distributions of dark matter and the scalar field, by solving their equations of motion using the multilevel adaptive grid technique. We show that the spatial configuration of the scalar field tracks both the voids and clusters of dark matter. The phiCDM model differs from standard LCDM at small scales with observable modifications of, e.g., the mass function of halos as well as the matter power spectrum. Nevertheless, the predictions of both models for the Hubble expansion and the CMB spectrum are virtually indistinguishable. Hence, galaxy cluster counts and weak lensing observations, which probe structure formation at small scales, are needed to falsify this class of models.

We have obtained deep photometry in two 1x1 degree fields covering the close pair of dwarf spheroidal galaxies (dSph) Leo IV and Leo V and part of the area in between. We find that both systems are significantly larger than indicated by previous measurements based on shallower data and also significantly elongated. With half-light radii of r_h=4'.6 +- 0'.8 (206 +- 36 pc) and r_h=2'.6 +- 0'.6 (133 +- 31 pc), respectively, they are now well within the physical size bracket of typical Milky Way dSph satellites. Their ellipticities of epsilon ~0.5 are shared by many faint (M_V>-8) Milky Way dSphs. The large spatial extent of our survey allows us to search for extra-tidal features with unprecedented sensitivity. The spatial distribution of candidate red giant branch and horizontal branch stars is found to be non-uniform at the ~3 sigma level. This substructure is aligned along the direction connecting the two systems, indicative of a possible `bridge' of extra-tidal material. Fitting the stellar distribution with a linear Gaussian model yields a significance of 4 sigma for this overdensity, a most likely FWHM of ~16 arcmin and a central surface brightness of ~32 mag arcsec^{-2}. We investigate different scenarios to explain the close proximity of Leo IV and Leo V and the possible tidal bridge between them. Orbit calculations demonstrate that they are unlikely to be remnants of a single disrupted progenitor, while a comparison with cosmological simulations shows that a chance collision between unrelated subhalos is negligibly small. Leo IV and Leo V could, however, be a bound `tumbling pair' if their combined mass exceeds 8 +- 4 x 10^9 M_sun. The scenario of an internally interacting pair appears to be the most viable explanation for this close celestial companionship. (abridged)

We use a statistical approach to determine the relationship between the stellar masses of galaxies and the masses of the dark matter halos in which they reside. We obtain a parameterized stellar-to-halo mass (SHM) relation by populating halos and subhalos in an N-body simulation with galaxies and requiring that the observed stellar mass function be reproduced. We find good agreement with constraints from galaxy-galaxy lensing and predictions of semi-analytic models. Using this mapping, and the positions of the halos and subhalos obtained from the simulation, we find that our model predictions for the galaxy two-point correlation function (CF) as a function of stellar mass are in excellent agreement with the observed clustering properties in the SDSS at z=0. We show that the clustering data do not provide additional strong constraints on the SHM function and conclude that our model can therefore predict clustering as a function of stellar mass. We compute the conditional mass function, which yields the average number of galaxies with stellar masses in the range [m, m+dm] that reside in a halo of mass M. We study the redshift dependence of the SHM relation and show that, for low mass halos, the SHM ratio is lower at higher redshift. The derived SHM relation is used to predict the stellar mass dependent galaxy CF and bias at high redshift. Our model predicts that not only are massive galaxies more biased than low mass ones at all redshifts, but the bias increases more rapidly with increasing redshift for massive galaxies than for low mass ones. We present convenient fitting functions for the SHM relation as a function of redshift, the conditional mass function, and the bias as a function of stellar mass and redshift.

We study the effect of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. We produce a suite of minor merger simulations for Milky Way-like galaxies. This suite consists of collisionless simulations as well as hydrodynamical runs including a gaseous component in the galactic disc. We find that in dissipationless simulations minor mergers cause the scale height of the disc to increase by up to a factor of ~2. When the presence of gas in the disc is taken into account this thickening is reduced by 25% (50%) for an initial disc gas fraction of 20% (40%), leading to a final scale height z0 between 0.6 and 0.7 kpc, regardless of the initial scale height. We argue that the presence of gas reduces disc heating via two mechanisms: absorption of kinetic impact energy by the gas and/or formation of a new thin stellar disc that can cause heated stars to recontract towards the disc plane. We show that in our simulations most of the gas is consumed during the merger and thus the regrowth of a new thin disc has a negligible impact on the z0 of the post merger galaxy. Final disc scale heights found in our simulations are in good agreement with studies of the vertical structure of spiral galaxies where the majority of the systems are found to have scale heights of 0.4 kpc < z0 < 0.8 kpc. We also found no tension between recent measurements of the scale height of the Milky Way thin disc and results coming from our hydrodynamical simulations. We conclude that the existence of a thin disc in the Milky Way and in external galaxies is not in obvious conflict with the predictions of the CDM model.

We study merger histories of dark-matter haloes in a suite of N-body simulations that span different cosmological models. The simulated cases include the up-to-date WMAP5 cosmology and other test cases based on the Einstein-deSitter cosmology with different power spectra. We provide a robust fitting function for the conditional mass function (CMF) of progenitor haloes of a given halo. This fit is valid for the different cosmological models and for different halo masses and redshifts, and it is a significant improvement over earlier estimates. Based on this fit, we develop a simple and accurate technique for transforming the merger history of a given simulated halo into haloes of different mass, redshift and cosmology. Other statistics such as main-progenitor history and merger rates are accurately transformed as well. This method can serve as a useful tool for studying galaxy formation. It is less sensitive to the low accuracy of the fit at small time-steps, and it can thus replace the more elaborate task of construction Monte-Carlo realizations. As an alternative approach, we confirm the earlier finding by Neistein & Dekel that the main-progenitor follows a log-normal distribution. This property of merger trees allows us to better capture their behaviour as a function of time and descendant mass, but a broader suite of simulations is required for evaluating the dependence of the log-normal parameters on the cosmological model.

We study the luminosity function and the radial distribution of satellite galaxies within Milky Way sized haloes as predicted in Cold Dark Matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic model (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the semi-analytic models. We present a detailed comparison of our predictions with the observational data recently obtained on the Milky Way satellite luminosity function (LF). We find that semi-analytic models with rather standard astrophysical ingredients are able to reproduce the observed luminosity function over six orders of magnitude in luminosity, down to magnitudes as faint as M_V=-2. We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photo-ionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of Milky Way satellites can be naturally accounted for within the (Lambda) Cold Dark Matter paradigm, and this should no longer be considered a problem.

We present a new universal relation, satisfied by matter distributions at all observed scales, and show its amazingly good and detailed agreement with the predictions of the most up-to-date pure dark matter simulations of structure formation in the Universe. This work extends the previous analysis [0904.4054; 0909.5203] to a larger range of masses, demonstrates a different scaling law, and compares it with numerical simulations. This behaviour seems to be insensitive to the complicated feedback of baryons on dark matter. Therefore, it potentially allows to compare theoretical predictions directly with observations, thus providing a new tool to constrain the properties of dark matter. Such a universal property, observed in structures of all sizes (from dwarf spheroidal galaxies to galaxy clusters), is difficult to explain without dark matter, thus providing new evidence for its existence.

Cosmological simulations of dark matter structures have shown that the equilibrated dark matter structures have a fairly small angular momentum. It appears from these N-body simulations that the radial profile of the angular momentum has an almost universal behavior, even if the different dark matter structures have experienced very different formation and merger histories. We suggest a perturbed Jeans equation, which includes a rotational term. This is done under a reasonable assumed form of the change in the distribution function. By conjecturing that the (new) subdominant rotation term must be proportional to the (old) dominant mass term, we find a clear connection, which is in rather good agreement with the results of recent high resolution simulations. We also present a new connection between the radial profiles of the angular momentum and the velocity anisotropy, which is also in fair agreement with numerical findings. Finally we show how the spin parameter $\lambda$ increases as a function of radius.

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.

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.

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.

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.

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.

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.

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.

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)

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<T_{IGM}<1e6 K) at a given redshift (z_r > 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.

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

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.

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.

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.

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.

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.

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 <delta_rho/rho < 0.6 on scales of 7 Mpc, in agreement with existing observational estimates on the local matter density. These results provide new, independent evidence for the presence of dark energy on scales of few Mpc, corroborating the evidence gathered from observations of distant objects and the early Universe.

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.

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.

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.

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.

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.

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.

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.