We discuss the formation of dark compact objects in a dark matter environment in view of the possible mass dependence of pulsars on the distribution of dark matter in the Galaxy. Our results indicate that the pulsar masses should decrease going towards the center of the Milky Way due to dark matter capture, thus becoming a probe for the existence and nature of dark matter. We thus propose that the evolution of the pulsar mass in a dark matter rich environment can be used to put constraints, when combined with future experiments, on the characteristics of our Galaxy halo dark matter profile, on the dark matter particle mass and on the dark matter self-interaction strength.

We show that the pulsar mass depends on the environment, and that it decreases going towards the center of the Milky Way. This is due to two combined effects, the capture and accumulation of self-interacting, non-annihilating dark matter by pulsars, and the increase of the dark matter density going towards the galactic center. We show that mass decrease depends both on the density profile of dark matter, steeper profiles producing a faster and larger decrease of the pulsar mass, and on the strength of self-interaction. Once future observations will provide the pulsar mass in a dark matter rich environment, close to the galactic center, the present result will be able to put constraints on the characteristics of our Galaxy halo dark matter profile, on the nature of dark matter, namely on its annihilating or non-annihilating nature, on its strength of self-interaction, and on the particle mass.

Cameron {\it et al.}~2019 (hereafter C19) recommends a more cautious and rigorous approach to statistical analysis in astronomy. We welcome this particular side of their communication as it helps stimulating the effort towards the adoption of better statistical methods in galaxy rotation curves, an effort to which we contributed with Rodrigues {\it et al.}~2018 (hereafter R18). Indeed, R18 was the first work that, in order to conclude on the universality of the acceleration scale $a_0$, studied the posterior distributions on $a_0$ of a large set of galaxies. As C19 agrees, the credible intervals were found within the Bayesian framework, that is, the marginalized posteriors on $a_0$ were found using Bayes' theorem to update the priors in light of the observational data; this process was done without introducing any approximation. Considering R18, C19 also remarks that: i) better methods to select the nuisance parameters and the corresponding priors could be used; ii) a quality cut based on $\chi^2$ values should not be used, and iii) the compatibility of the posteriors should be assessed in a more robust way. In the following, after first clarifying the context of our work, we address these criticisms.

Recently, the discovery of the radial acceleration relation (RAR) in galaxies has been regarded as an indirect support of alternative theories of gravity such as Modified Newtonian Dynamics (MOND) and modified gravity. This relation indicates a tight correlation between dynamical mass and baryonic mass in galaxies with different sizes and morphology. However, if the RAR relation is scale-independent and could be explained by alternative theories of gravity, this relation should be universal and true for galaxy clusters as well. In this article, by using the x-ray data of a sample of galaxy clusters, we investigate if there exists any tight correlation between dynamical mass and baryonic mass in galaxy clusters, assuming hot gas mass distribution almost representing baryonic distribution and that the galaxy clusters are virialized. We show that the resulting RAR of 52 non-cool-core galaxy clusters scatters in a large parameter space, possibly due to our simplifying assumptions and unclear matter content in galaxy clusters. This might indicate that the RAR is unlikely to be universal and scale-independent.

The thermodynamic equilibrium condition for a static self-gravitating fluid in the Einstein theory is defined by the Tolman-Ehrenfest temperature law, $T{\sqrt {g_{00}(x^{i})}} = constant$, according to which the proper temperature depends explicitly on the position within the medium through the metric coefficient $g_{00}(x^{i})$. By assuming the validity of Tolman-Ehrenfest "pocket temperature", Klein also proved a similar relation for the chemical potential, namely, $\mu {\sqrt {g_{00}(x^{i})}} = constant$. In this letter we prove that a more general relation uniting both quantities holds regardless of the equation of state satisfied by the medium, and that the original Tolman-Ehrenfest law form is valid only if the chemical potential vanishes identically. In the general case of equilibrium, the temperature and the chemical potential are intertwined in such a way that only a definite (position dependent) relation uniting both quantities is obeyed. As an illustration of these results, the temperature expressions for an isothermal gas (finite spherical distribution) and a neutron star are also determined.

The evolution of the virial overdensity $\Delta_{\rm vir}$ for $\Lambda$CDM and seven dynamical dark-energy models is investigated in the extended spherical collapse model (SCM). Here the virialization process is naturally achieved by introducing shear and rotation instead of using the virial theorem. We generalise two approaches proposed in the literature and show that, regardless of the dark-energy model, the new virialization term can be calibrated on the peculiar velocity of the shell as measured from Einstein-de Sitter simulations. The two virialization recipes qualitatively reproduce the features of the ordinary SCM, i.e., a constant $\Delta_{\rm vir}$ for the EdS model and time-variation for dark-energy models, but without any mass dependence. Depending on the actual description of virialization and on the dark-energy model, the value of $\Delta_{\rm vir}$ varies between 10 and 40 percent. We use the new recipes to predict the surface-mass-density profile of dark matter haloes and the number of convergence density peaks for LSST- and Euclid-like weak lensing surveys.

We derive the mass-temperature relation using an improved top-hat model and a continuous formation model which takes into account the effects of the ordered angular momentum acquired through tidal-torque interaction between clusters, random angular momentum, dynamical friction, and modifications of the virial theorem to include an external pressure term usually neglected. We show that the mass-temperature relation differs from the classical self-similar behavior, $M \propto T^{3/2}$, and shows a break at $3--4$ keV, and a steepening with a decreasing cluster temperature. We then compare our mass-temperature relation with those obtained in the literature with $N$-body simulations for $f(R)$ and symmetron models. We find that the mass-temperature relation is not a good probe to test gravity theories beyond Einstein's general relativity, because the mass-temperature relation of the $\Lambda$CDM model is similar to that of the modified gravity theories.

We study the total and dark matter (DM) density profiles as well as their correlations for a sample of 15 high-mass galaxy clusters by extending our previous work on several clusters from Newman et al. Our analysis focuses on 15 CLASH X-ray-selected clusters that have high-quality weak- and strong-lensing measurements from combined Subaru and {\em Hubble Space Telescope} observations. The total density profiles derived from lensing are interpreted based on the two-phase scenario of cluster formation. In this context, the brightest cluster galaxy (BCG) forms in the first dissipative phase, followed by a dissipationless phase where baryonic physics flattens the inner DM distribution. This results in the formation of clusters with modified DM distribution and several correlations between characteristic quantities of the clusters. We find that the central DM density profiles of the clusters are strongly influenced by baryonic physics as found in our earlier work. The inner slope of the DM density for the CLASH clusters is found to be flatter than the Navarro--Frenk--White profile, ranging from $\alpha=0.30$ to $0.79$. We examine correlations of the DM density slope $\alpha$ with the effective radius $R_\mathrm{e}$ and stellar mass $M_\mathrm{e}$ of the BCG, finding that these quantities are anti-correlated with a Spearman correlation coefficient of $\sim -0.6$. We also study the correlation between $R_\mathrm{e}$ and the cluster halo mass $M_{500}$, and the correlation between the total masses inside 5\,kpc and 100\,kpc. We find that these quantities are correlated with Spearman coefficients of $0.68$ and $0.64$, respectively. These observed correlations are in support of the physical picture proposed by Newman et al.

We study the structure of compact objects that contain non-self annihilating, self-interacting dark matter admixed with ordinary matter made of neutron star and white dwarf materials. We extend the previous work Phys. Rev. D 92 123002 (2015) on these dark compact objects by analyzing the effect of weak and strongly interacting dark matter with particle masses in the range of 1-500 GeV, so as to set some constraints in the strength of the interaction and the mass of the dark matter particle. We find that the total mass of the compact objects increases with decreasing dark matter particle mass. In the strong interacting case and for dark matter particle masses in the range 1-10 GeV, the total mass of the compact objects largely exceeds the $2M_\odot$ constraint for neutron star masses and the nominal $1M_\odot$for white dwarfs, while for larger dark matter particle masses or in the weakly interacting case the compact objects show masses in agreement or smaller than these constraints, thus hinting at the exclusion of strongly self-interacting dark matter of masses 1-10 GeV in the interior of these compact objects. Moreover, we observe that the smaller the dark matter particle mass, the larger the quantity of dark matter captured is, putting constraints on the dark matter mass trapped in the compact objects so as to fullfill $\simeq 2M_\odot$ observations. Finally, the inhomogeneity of distribution of dark matter in the Galaxy implies a mass dependence of compact objects from the environment which can be used to put constraints on the characteristics of the Galaxy halo DM profile and on particle mass. In view of the these results, we discuss the formation of the dark compact objects in an homogeneous and non-homogeneous dark matter environment.

This correspondence answers the criticisms raised in McGaugh et al 2018 and Kroupa et al 2018.

Using a semianalytic code, we show how baryon physics in a $\Lambda$CDM cosmology could solve the discrepancy between numerical predictions of dark matter haloes and observations, ranging from dwarf galaxies to clusters, without the need of nonstandard dark matter models as advocated, for example, by [Kaplinghat et al., Phys. Rev. Lett. 116, 041302, (2016)]. Combining well established results, we show, for the first time, how accounting for baryon physics, in particular dynamical friction mechanisms, leads to flat galaxy-cluster profiles and correlations in several of their properties, solves the so-called `diversity problem' and reproduces very well the challenging, extremely low-rising rotation curve of IC2574. We therefore suggest treating baryonic physics properly before introducing new exotic features, albeit legitimate, in the standard cosmological model.

The Radial Acceleration Relation confirms that a nontrivial acceleration scale can be found in the average internal dynamics of galaxies. The existence of such a scale is not obvious as far as the standard cosmological model is concerned, and it has been interpreted as a possible sign of modified gravity. The implications could be profound: it could in principle explain galactic dynamics without large amounts of yet-undetected dark matter and address issues that the standard cosmological model faces at galactic scales. Here, we consider 193 disk galaxies from the SPARC and THINGS databases and, using Bayesian inference, we show that the probability of existence of a fundamental acceleration that is common to all the galaxies is essentially zero: the $p$-value is smaller than $10^{-20}$ or, equivalently, the null hypothesis is rejected at more than 10$\sigma$. We conclude that the acceleration scale unveiled by the Radial Acceleration Relation is of emergent nature, possibly caused by a complex interplay between baryons and dark matter. In particular, the MOND theory, or any other theory that behaves like it at galactic scales, is ruled out as a fundamental theory for galaxies at more than 10$\sigma$.

Aims. We aim to study the stochastic evolution of the smoothed overdensity $\delta$ at scale $S$ of the form $\delta(S) = \int_{0}^S K(S,u)\mathrm{d}W(u)$, where $K$ is a kernel and $\mathrm{d}W$ is the usual Wiener process. Methods. For a Gaussian density field, smoothed by the top-hat filter, in real space, we used a simple kernel that gives the correct correlation between scales. A Monte Carlo procedure was used to construct random walks and to calculate first crossing distributions and consequently mass functions for a constant barrier. Results. We show that the evolution considered here improves the agreement with the results of N-body simulations relative to analytical approximations which have been proposed from the same problem by other authors. In fact, we show that an evolution which is fully consistent with the ideas of the excursion set model, describes accurately the mass function of dark matter haloes for values of $\nu \leq 1$ and underestimates the number of larger haloes. Finally, we show that a constant threshold of collapse, lower than it is usually used, it is able to produce a mass function which approximates the results of N-body simulations for a variety of redshifts and for a wide range of masses. Conclusions. A mass function in good agreement with N-body simulations can be obtained analytically using a lower than usual constant collapse threshold.

Proper motions of collisionless pointlike objects in a spherically symmetric system, for example stars in a galaxy, can be used to test whether that system is in equilibrium, with no assumptions regarding isotropy. In particular, the fourth order spherical Jeans equation yields expressions for two observable quantities characterizing the departure from equilibrium, both of which can be expressed in terms of time derivatives of first and third moments of the velocities. As illustrations, we compute these quantities for tracer distributions drawn from an exact equilibrium configuration and also from near equilibrium configurations generated using the N-body code GALIC.

We develop and apply new techniques in order to uncover galaxy rotation curves (RC) systematics. Considering that an ideal dark matter (DM) profile should yield RCs that have no bias towards any particular radius, we find that the Burkert DM profile satisfies the test, while the Navarro-Frenk-While (NFW) profile has a tendency of better fitting the region between one and two disc scale lengths than the inner disc scale length region. Our sample indicates that this behaviour happens to more than 75% of the galaxies fitted with an NFW halo. Also, this tendency does not weaken by considering "large" galaxies, for instance those with $M_*\gtrsim 10^{10} M_\odot$. Besides the tests on the homogeneity of the fits, we also use a sample of 62 galaxies of diverse types to perform tests on the quality of the overall fit of each galaxy, and to search for correlations with stellar mass, gas mass and the disc scale length. In particular, we find that only 13 galaxies are better fitted by the NFW halo; and that even for the galaxies with $M_* \gtrsim 10^{10} M_\odot$ the Burkert profile either fits as good as, or better than, the NFW profile. This result is relevant since different baryonic effects important for the smaller galaxies, like supernova feedback and dynamical friction from baryonic clumps, indicate that at such large stellar masses the NFW profile should be preferred over the Burkert profile. Hence, our results either suggest a new baryonic effect or a change of the dark matter physics.

The $\Lambda$CDM model, or concordance cosmology, as it is often called, is a paradigm at its maturity. It is clearly able to describe the universe at large scale, even if some issues remain open, such as the cosmological constant problem , the small-scale problems in galaxy formation, or the unexplained anomalies in the CMB. $\Lambda$CDM clearly shows difficulty at small scales, which could be related to our scant understanding, from the nature of dark matter to that of gravity; or to the role of baryon physics, which is not well understood and implemented in simulation codes or in semi-analytic models. At this stage, it is of fundamental importance to understand whether the problems encountered by the $\Lambda$DCM model are a sign of its limits or a sign of our failures in getting the finer details right. In the present paper, we will review the small-scale problems of the $\Lambda$CDM model, and we will discuss the proposed solutions and to what extent they are able to give us a theory accurately describing the phenomena in the complete range of scale of the observed universe.

In this paper, extending past works of Del Popolo, we show how a high precision mass function (MF) can be obtained using the excursion set approach and an improved barrier taking implicitly into account a non-zero cosmological constant, the angular momentum acquired by tidal interaction of proto-structures and dynamical friction. In the case of the $\Lambda$CDM paradigm, we find that our MF is in agreement at the 3\% level to Klypin's Bolshoi simulation, in the mass range $M_{\rm vir} = 5 \times 10^9 h^{-1} M_{\odot} -- 5 \times 10^{14} h^{-1} M_{\odot}$ and redshift range $0 \lesssim z \lesssim 10$. For $z=0$ we also compared our MF to several fitting formulae, and found in particular agreement with Bhattacharya's within 3\% in the mass range $10^{12}-10^{16} h^{-1} M_{\odot}$. Moreover, we discuss our MF validity for different cosmologies.

We construct an integral equation for the first crossing distributions for fractional Brownian motion in the case of a constant barrier and we present an exact analytical solution. Additionally we present first crossing distributions derived by simulating paths from fractional Brownian motion. We compare the results of the analytical solutions with both those of simulations and those of some approximated solutions which have been used in the literature. Finally, we present multiplicity functions for dark matter structures resulting from our analytical approach and we compare with those resulting from N-body simulations. We show that the results of analytical solutions are in good agreement with those of path simulations but differ significantly from those derived from approximated solutions. Additionally, multiplicity functions derived from fractional Brownian motion are poor fits of the those which result from N-body simulations. We also present comparisons with other models which exist in the literature and we discuss different ways of improving the agreement between analytical results and N-body simulations.

A precise determination of the mass function is an important tool to verify cosmological predictions of the $\Lambda$CDM model and to infer more precisely the better model describing the evolution of the Universe. Galaxy clusters have been currently used to infer cosmological parameters, in particular the matter density parameter $\Omega_{\rm m}$, the matter power spectrum normalization $\sigma_8$ and the equation of state parameter $w_{\rm de}$ of the dark energy fluid. In this work, using data on massive galaxy clusters ($M>8\times 10^{14}~h^{-1}~M_{\odot}$) in the redshift range $0.05\lesssim z\lesssim 0.83$ we put constraints on the parameter $\alpha$ introduced within the formalism of the extended spherical collapse model to quantify deviations from sphericity due to shear and rotation. Since at the moment there is no physical model describing its functional shape, we assume it to be a logarithmic function of the cluster mass. By holding $\sigma_8$ fixed and restricting our analysis to a $\Lambda$CDM model, we find, at $1-\sigma$ confidence level, $\Omega_{\rm m}=0.284\pm0.0064$, $h=0.678\pm0.017$ and $\beta=0.0019^{+0.0008}_{-0.0015}$, where $\beta$ represents the slope of the parameter $\alpha$. This results translates into a $9\%$ decrement of the number of massive clusters with respect to a standard $\Lambda$CDM mass function, but better data are required to better constrain this quantity, since at the $2-\sigma$ and $3-\sigma$ confidence level we are only able to infer upper limits.

In the present paper, we extend the study of Del Popolo (2010) to determine the slope of the inner density profile of galaxy haloes with different morphologies. We study how galaxy morphology changes the relation between the inner slope of the galaxy halo density profile, $\alpha$, and the stellar mass, $M_{*}$, or rotation velocity $V_{\rm rot}$. For this, we use the model of Del Popolo (2009) in combination with observed data from the Romanowsky \& Fall (2012) sample of elliptical and spiral galaxies, the Local Group sample compiled by McConnachie (2012), and the simulation results by Cloet-Osselaer et al. (2014). We find that the slope $\alpha$ flattens monotonically, from $\alpha \simeq -1 $ at $V_{\rm rot} \simeq 250$ km/s, to $\alpha \simeq 0 $. After $V_{\rm rot}\simeq 25$ km/s the slope starts to steepen. The steepening happens in the mass range dominated by non-rotationally supported galaxies (e.g., dSphs) and depends on the level of offset in the angular momentum of rotationally and non-rotationally dominated galaxies. The steepening is a consequence of the decrease in baryons content, and angular momentum in spheroidal dwarf galaxies. We finally compare our result to the SPH simulations of Di Cintio. Our result is in qualitatively agreement with their simulations, with the main difference that the inner slope $\alpha$ at small stellar masses ($M_* \lesssim10^{8} M_{\odot}$) is flatter than that in their simulations. As a result, the claim that finding a core in dwarf galaxies with masses slightly smaller than $\simeq 10^6 M_{\odot}$, (as in the Di Cintio, or Governato, supernovae feedback mechanism) would be a problem for the $\Lambda$CDM model must be probably revised.

In the present paper, we compare the predictions of two well known mechanisms considered able to solve the cusp/core problem (a. supernova feedback; b. baryonic clumps-DM interaction) by comparing their theoretical predictions to recent observations of the inner slopes of galaxies with masses ranging from dSphs to normal spirals. We compare the $\alpha$-$V_{\rm rot}$ and the $\alpha$-$M_{\ast}$ relationships, predicted by the two models with high resolution data coming from \citep{Adams2014,Simon2005}, LITTLE THINGS \citep{Oh2015}, THINGS dwarves \citep{Oh2011a,Oh2011b}, THINGS spirals \citep{Oh2015}, Sculptor, Fornax and the Milky Way. The comparison of the theoretical predictions with the complete set of data shows that the two models perform similarly, while when we restrict the analysis to a smaller subsample of higher quality, we show that the method presented in this paper (baryonic clumps-DM interaction) performs better than the one based on supernova feedback. We also show that, contrarily to the first model prediction, dSphs of small mass could have cored profiles. This means that observations of cored inner profiles in dSphs having a stellar mass $<10^6 M_{\odot}$ not necessarily imply problems for the $\Lambda$CDM model.

The motion of a point like object of mass $M$ passing through the background potential of massive collisionless particles ($m << M$) suffers a steady deceleration named dynamical friction. In his classical work, Chandrasekhar assumed a Maxwellian velocity distribution in the halo and neglected the self gravity of the wake induced by the gravitational focusing of the mass $M$. In this paper, by relaxing the validity of the Maxwellian distribution due to the presence of long range forces, we derive an analytical formula for the dynamical friction in the context of the $q$-nonextensive kinetic theory. In the extensive limiting case ($q = 1$), the classical Gaussian Chandrasekhar result is recovered. As an application, the dynamical friction timescale for Globular Clusters spiraling to the galactic center is explicitly obtained. Our results suggest that the problem concerning the large timescale as derived by numerical $N$-body simulations or semi-analytical models can be understood as a departure from the standard extensive Maxwellian regime as measured by the Tsallis nonextensive $q$-parameter.

Extended dark matter (DM) substructures may play the role of microlenses in the Milky Way and in extragalactic gravitational lens systems (GLSs). We compare microlensing effects caused by point masses (Schwarzschild lenses) and extended clumps of matter using a simple model for the lens mapping. A superposition of the point mass and the extended clump is also considered. For special choices of the parameters, this model may represent a cusped clump of cold DM, a cored clump of self-interacting dark matter (SIDM) or an ultra compact minihalo of DM surrounding a massive point-like object. We built the resulting micro-amplification curves for various parameters of one clump moving with respect to the source in order to estimate differences between the light curves caused by clumps and by point lenses. The results show that it may be difficult to distinguish between these models. However, some region of the clump parameters can be restricted by considering the high amplification events at the present level of photometric accuracy. Then we estimate the statistical properties of the amplification curves in extragalactic GLSs. For this purpose, an ensemble of amplification curves is generated yielding the autocorrelation functions (ACFs) of the curves for different choices of the system parameters. We find that there can be a significant difference between these ACFs if the clump size is comparable with typical Einstein radii; as a rule, the contribution of clumps makes the ACFs less steep.

We continue the study of the impact of baryon physics on the small scale problems of the $\Lambda$CDM model, based on a semi-analytical model (Del Popolo, 2009). Withsuch model, we show how the cusp/core, missing satellite (MSP), Too Big to Fail (TBTF) problems and the angular momentum catastrophe can be reconciled with observations, adding parent-satellite interaction. Such interaction between darkmatter (DM) and baryons through dynamical friction (DF) can sufficiently flattenthe inner cusp of the density profiles to solve the cusp/core problem. Combining, in our model, a Zolotov et al. (2012)-like correction, similarly to Brooks et al. (2013), and effects of UV heating and tidal stripping, the number of massive, luminous satellites, as seen in the Via Lactea 2 (VL2) subhaloes,is in agreement with the numbers observed in the MW, thus resolving the MSP and TBTF problems. The model also produces a distribution of the angular spin parameter and angular momentum in agreement with observations of the dwarfs studied by van den Bosch, Burkert, \\& Swaters (2001).

In this article, we test the conclusion of Donato et al. (2009) concerning the universality of the DM halo surface density $\mu_{0D}=\rho_0r_0$. According to our study, the dispersion of values of $\mu_{0D}$ is twice higher than that found by Donato et al. (2009). We conclude, in contrast with Donato et al. (2009), that the DM surface density and its Newtonian acceleration are not constant but correlate with the luminosity, morphological type, $(B-V)_0$ colour index, and the content of neutral hydrogen. These DM parameters are higher for more luminous systems of early types with red colour and low gas content. We also found that the correlation of DM parameters with colour index appears to be the manifestation of a stronger relation between DM halo mass and the colour of a galaxy. This finding is in agreement with cosmological simulations (Guo et al, 2011). These results leave little room for the recently claimed universality of DM column density. We also found that isolated galaxies in our sample (contained in the Analysis of the interstellar Medium of Isolated GAlaxies (AMIGA) catalogue) do not differ significantly in their value of $\mu_{0D}$ from the entire sample. Thus, since the AMIGA catalogue gives a sample of galaxies that have not interacted with a significant mass neighbour in the past 3 Gyr, the difference between the systems with low and high values of $\mu_{0D}$ is not related to the merging events during this period of time.

In the framework of the spherical collapse model we study the influence of shear and rotation terms for dark matter fluid in clustering dark energy models. We evaluate, for different equations of state, the effects of these terms on the linear overdensity threshold parameter, $\delta_{\rm c}$, and on the virial overdensity, $\Delta_{\rm V}$. The evaluation of their effects on $\delta_{\rm c}$ allows us to infer the modifications occurring on the mass function. Due to ambiguities in the definition of the halo mass in the case of clustering dark energy, we consider two different situations: the first is the classical one where the mass is of the dark matter halo only, while the second one is given by the sum of the mass of dark matter and dark energy. As previously found, the spherical collapse model becomes mass dependant and the two additional terms oppose to the collapse of the perturbations, especially on galactic scales, with respect to the spherical non-rotating model, while on clusters scales the effects of shear and rotation become negligible. The values for $\delta_{\rm c}$ and $\Delta_{\rm V}$ are higher than the standard spherical model. Regarding the effects of the additional non-linear terms on the mass function, we evaluate the number density of halos. As expected, major differences appear at high masses and redshifts. In particular, quintessence (phantom) models predict more (less) objects with respect to the $\Lambda$CDM model and the mass correction due to the contribution of the dark energy component has negligible effects on the overall number of structures.

In the present paper, we studied by means of the SIM introduced in Del Popolo (2009), the total and DM density profiles, and the correlations among different quantities, observed by Newman et al. (2012a,b), in seven massive and relaxed clusters, namely MS2137, A963, A383, A611, A2537, A2667, A2390. Similarly to Newman et al. (2012a,b), the total density profile, in the radius range 0.003 - 0.03$r_{200}$, has a mean total density profile in agreement with dissipationless simulations. The slope of the DM profiles of all clusters is flatter than -1. The slope, $\alpha$, has a maximum value (including errors) of $\alpha=-0.88$ in the case of A2390, and minimum value $\alpha=-0.14$ for A2537. The baryonic component dominates the mass distribution at radii $< 5-10$ kpc, while the outer distribution is dark matter dominated. We found an anti-correlation among the slope $\alpha$, the effective radius, $R_e$, and the BCG mass, and a correlation among the core radius $r_{core}$, and $R_e$. Moreover, the mass in 100 kpc (mainly dark matter) is correlated with the mass inside 5 kpc (mainly baryons). The behavior of the total mass density profile, the DM density profile, and the quoted correlations can be understood in a double phase scenario. In the first dissipative phase the proto-BCG forms, and in the second dissipationless phase, dynamical friction between baryonic clumps (collapsing to the center) and the DM halo flattens the inner slope of the density profile. In simple terms, the large scatter in the inner slope from cluster to cluster, and the anti-correlation among the slope, $\alpha$ and $R_e$ is due to the fact that in order to have a total mass density profile which is NFW-like, clusters having more massive BCGs at their centers must contain less DM in their center. Consequently the inner profile has a flatter slope.

We study, by means of the model proposed in Del Popolo (2009), the effect of baryon physics on the small scale problems of the CDM model. We show that, using this model, the cusp/core problem, the missing satellite problem (MSP), the Too Big to Fail (TBTF) problem, and the angular momentum catastrophe can be reconciled with observations. Concerning the cusp/core problem, the interaction among dark matter (DM) and baryonic clumps of 1% the mass of the halo, through dynamical friction (DF), is able to flatten the inner cusp of the density profiles. We moreover assume that haloes form primarily through quiescent accretion, in agreement with the spherical collapse model (SCM)-secondary infall model (SIM) prescriptions. The results of this paper follow from the two assumptions above. Concerning the MSP and TBTF problem, applying to the Via Lactea II (VL2) subhaloes a series of corrections similar to those of Brooks et al. (2013), namely applying a Zolotov et al. (2012)-like correction obtained with our model, and further correcting for the UV heating and tidal stripping, we obtain that the number of massive, luminous satellites is in agreement with the number observed in the MW. The model also produces an angular momentum distribution in agreement with observations, that is with the distribution of the angular spin parameter and angular momentum of the dwarfs studied by van den Bosch, Burkert, & Swaters (2001). In conclusion, the small scale problems of the CDM model can all be solved by introducing baryon physics.

In this paper, we study how the presence of bulge formation in galaxies influence their inner density profile, by means of an extended version of the Del Popolo (2009) semi-analytical model. As in Del Popolo (2009), the model takes into account the effect of baryons adiabatic contraction, ordered and random angular momentum, dynamical friction, and adds to the previous the effect of gas cooling, star formation, supernova feedback, and reionization. Our model shows that dwarf galaxies are bulgeless, in agreement with observations showing that the large majority of them has no stellar bulges, and are characterized by a flat profile well described by a Burkert profile. {We then studied the effect of a bulge, added to the {cored} DM halo, on the density profile}. In the case of a galaxy having a mass $10^{11} M_{\odot}$ the inner density profile has a slope $\alpha \simeq 0.65$, for a bulge of $4.5 \times 10^{9} M_{\odot}$, while if bulge formation is not considered, the slope would be $\alpha \simeq 0.55$. If the bulge is larger, $6.5 \times 10^{9} M_{\odot}$ the slope is $\alpha \simeq 0.7$. In the case of a larger galaxy with $10^{12} M_{\odot}$ the slope is $\alpha \simeq 0.85$, while in absence of bulge it is $\alpha \simeq 0.75$. We finally study how the inner slope $\alpha$ changes with the bulge mass, and we find a correlation among the two quantities. The result shows that bulge formation has an important role in shaping the inner DM density profile in agreement with Inoue & Saitoh (2011) result. The result implies that going from Sc to SO Hubble type the slope is slightly steepening due to the bulge formation, and due to the fact that early type galaxies have larger bulges.

This paper is a broad-band review of the current status of non-baryonic dark matter research. I start with a historical overview of the evidences of dark matter existence, then I discuss how dark matter is distributed from small scale to large scale, and I then verge the attention to dark matter nature: dark matter candidates and their detection. I finally discuss some of the limits of the $\Lambda$CDM model, with particular emphasis on the small scale problems of the paradigm.

We study the effect of shear and rotation on results previously obtained dealing with the application of the spherical collapse model (SCM) to generalized Chaplygin gas (gCg) dominated universes. The system is composed of baryons and gCg and the collapse is studied for different values of the parameter $\alpha$ of the gCg. We show that the joint effect of shear and rotation is that of slowing down the collapse with respect to the simple SCM. This result is of utmost importance for the so-called unified dark matter models, since the described slow down in the growth of density perturbation can solve one of the main problems of the quoted models, namely the instability described in previous papers [e.g., H. B. Sandvik {\it et al.}, Phys. Rev. D {\bf 69}, 123524 (2004)] at the linear perturbation level.

In the paper the closed Friedmann-Robertson-Walker model with quantization in the presence of the positive cosmological constant and radiation is studied. For analysis of tunneling probability for birth of an asymptotically deSitter, inflationary Universe as a function of the radiation energy a new definition of a "free" wave propagating inside strong fields is proposed. On such a basis, tunneling boundary condition is corrected, penetrability and reflection concerning to the barrier are calculated in fully quantum stationary approach. For the first time non-zero interference between the incident and reflected waves has been taken into account which turns out to play important role inside cosmological potentials and could be explained by non-locality of barriers in quantum mechanics. Inside whole region of energy of radiation the tunneling probability for the birth of the inflationary Universe is found to be close to its value obtained in semiclassical approach. The reflection from the barrier is determined for the first time (which is essentially differs on 1 at the energy of radiation close to the barrier height). The proposed method could be easily generalized on the cosmological models with the barriers of arbitrary shape, that has been demonstrated for the FRW-model with included Chaplygin gas. Result is stable for variations of the studied barriers, accuracy are found to be 11--18 digits for all coefficients and energies below the barrier height.

In this paper, in the framework of the secondary infall model, the correlation between the central surface density and the halo core radius of galaxy, and cluster of galaxies, dark matter haloes was analyzed, this having recently been studied on a wide range of scales. We used Del Popolo (2009) secondary infall model taking into account ordered and random angular momentum, dynamical friction, and dark matter (DM) adiabatic contraction to calculate the density profile of haloes, and then these profiles are used to determine the surface density of DM haloes. The main result is that $r_\ast$ (the halo characteristic radius) is not an universal quantity as claimed by Donato et al. (2009) and Gentile et al. (2009). On the contrary, we find a correlation with the halo mass $M_{200}$ in agreement with Cardone & Tortora (2010), Boyarsky at al. (2009) and Napolitano et al. (2010), but with a significantly smaller scatter, namely $0.16 \pm 0.05$. We also consider the baryon column density finding this latter being indeed a constant for low mass systems such as dwarfs, but correlating with mass with a slope $\alpha= 0.18 \pm 0.05$. In the case of the surface density of dark matter for a system composed only of dark matter, as in dissipationless simulations, we get $\alpha=0.20 \pm 0.05$. These results leave little room for the recently claimed universality of (dark and stellar) column density.

We study, for the first time, how shear and angular momentum modify typical parameters of the spherical collapse model, in dark energy dominated universes. In particular, we study the linear density threshold for collapse $\delta_\mathrm{c}$ and the virial overdensity $\Delta_\mathrm{V}$, for several dark-energy models and its influence on the cumulative mass function. The equations of the spherical collapse are those obtained in Pace et al. (2010), who used the fully nonlinear differential equation for the evolution of the density contrast derived from Newtonian hydrodynamics, and assumed that dark energy is present only at the background level. With the introduction of the shear and rotation terms, the parameters of the spherical collapse model are now mass-dependant. The results of the paper show, as expected, that the new terms considered in the spherical collapse model oppose the collapse of perturbations on galactic scale giving rise to higher values of the linear overdensity parameter with respect to the non-rotating case. We find a similar effect also for the virial overdensity parameter. For what concerns the mass function, we find that its high mass tail is suppressed, while the low mass tail is slightly affected except in some cases, e.g. the Chaplygin gas case.

We present results based on Chandra observations of a large sample of 129 hot galaxy clusters. We measure the concentration parameter c_200, the dark mass M_200 and the baryonic mass content in all the objects of our sample, providing the largest dataset of mass parameters for galaxy clusters in the redshift range z = 0.01 - 1.4. We confirm a that a tight correlation between c_200 and M_200, c \propto M^a_vir /(1+z)^b with a = -0.56 +/- 0.15 and b =0.80 +/- 0.25 (68 per cent confidence limits), is present, in good agreement with the predictions from numerical simulations and previous observations. Fitting the mass profile with a generalized NFW model, we got the inner slope alpha, with alpha = 0.94 +/- 0.13. Finally, we show that the inner slope of the density profile, alpha correlates with the baryonic mass content, M_b : namely alpha is decreasing with increasing baryonic mass content.

We revisit the issue of the constancy of the dark matter (DM) and baryonic Newtonian acceleration scales within the DM scale radius by considering a large sample of late - type galaxies. We rely on a Markov Chain Monte Carlo (MCMC) method to estimate the parameters of the halo model and the stellar mass - to - light ratio and then propagate the uncertainties from the rotation curve data to the estimate of the acceleration scales. This procedure allows us to compile a catalog of 58 objects with estimated values of the $B$ band absolute magnitude $M_B$, the virial mass $M_{vir}$, the DM and baryonic Newtonian accelerations (denoted as $g_{DM}(r_0)$ and $g_{bar}(r_0)$, respectively) within the scale radius $r_0$ which we use to investigate whether it is possible to define a universal acceleration scale. We find a weak but statistically meaningful correlation with $M_{vir}$ thus making us argue against the universality of the acceleration scales. However, the results somewhat depend on the sample adopted so that a careful analysis of selection effects should be carried out before any definitive conclusion can be drawn.

The influence of the shear stress and angular momentum on the nonlinear spherical collapse model is discussed in the framework of the Einstein-de Sitter (EdS) and $\Lambda$CDM models. By assuming that the vacuum component is not clustering within the homogeneous nonspherical overdensities, we show how the local rotation and shear affects the linear density threshold for collapse of the non-relativistic component ($\delta_\mathrm{c}$) and its virial overdensity ($\Delta_\mathrm{V}$). It is also found that the net effect of shear and rotation in galactic scale is responsible for higher values of the linear overdensity parameter as compared with the standard spherical collapse model (no shear and rotation).

The present paper extends to clusters of galaxies the study of Del Popolo (2012), concerning how the baryon-dark matter (DM) interplay shapes the density profile of dwarf galaxies. Cluster density profiles are determined taking into account dynamical friction, random and ordered angular momentum and the response of dark matter halos to condensation of baryons. We find that halos containing only DM are characterized by Einasto's profiles, and that the profile flattens with increasing content of baryons, and increasing values of random angular momentum. The analytical results obtained in the first part of the paper were applied to well studied clusters whose inner profiles have slopes flatter than NFW predictions (A611, A383) or are characterized by profiles in agreement with the NFW model (MACS J1423.8+2404, RXJ1133). By using independently-measured baryonic fraction, a typical spin parameter value $\lambda \simeq 0.03$, and adjusting the random angular momentum, we re-obtain the mass and density profiles of the quoted clusters. Finally, we show that the baryonic mass inside $\simeq 10$ kpc, $M_{b,in}$ is correlated with the total mass of the clusters, %finding a correlation among the two quantities, as $M_{b,in} \propto M_{500}^{0.4}$.

In this paper, we give a possible solution to the cosmological constant problem. It is shown hat the traditional approach, based on volume weighting of probabilities, leads to an incoherent onclusion: the probability that a randomly chosen observer measures $\Lambda=0$ is exactly equal to 1. Using an alternative, volume averaging measure, instead of volume weighting can explain why the osmological constant is non-zero.

According to the now strongly supported concordance $\Lambda$CDM model, galaxies may be grossly described as a luminous component embedded in a dark matter halo. The density profile of these mass dominating haloes may be determined by N - body simulations which mimic the evolution of the tiny initial density perturbations during the process leading to the structures we observe today. Unfortunately, when the effect of baryons is taken into account, the situation gets much more complicated due to the difficulties in simulating their physics. As a consequence, a definitive prediction of how dark matter haloes should presently look like is still missing. We revisit here this issue from an observational point of view devoting our attention to dwarf galaxies. Being likely dark matter dominated, these systems are ideal candidates to investigate the present day halo density profiles and check whether dark matter related quantities correlate with the stellar ones or the environment. By fitting a large sample of well measured rotation curves, we infer constraints on both halo structural parameters (such as the logarithmic slope of the density profile and its concentration) and derived quantities (e.g., the mass fraction and the Newtonian acceleration) which could then be used to constrain galaxy formation scenarios. Moreover, we investigate whether the halo properties correlates with the environment the galaxy lives in thus offering a new tool to deepen our understanding of galaxy formation.

In the present paper, we study how the dark matter density profiles of dwarfs galaxies in the mass range $10^8-10^{10} M_{\odot}$ are modified by the interaction of the dwarf in study with the neighboring structures, and by changing baryon fraction in dwarfs. As already shown in Del Popolo (2009), the slope of density profile of inner halos flattens with decreasing halo mass and the profile is well approximated by a Burkert's profile. The analysis shows that dwarfs who suffered a smaller tidal torquing (consequently having smaller angular momentum) are characterized by steeper profiles with respect to dwarfs subject to higher torque, and similarly dwarfs having a smaller baryons fraction have also steeper profiles than those having a larger baryon fraction. In the case tidal torquing is shut down and baryons are not present, the density profile is very well approximated by an Einasto profile, similarly to dwarfs obtained in dissipationless N-body simulations. We then apply the result of the previous analysis to the dark matter halo rotation curves of three different dwarfs, namely NGC 2976, known to have a flat inner core, NGC 5949 having a profile intermediate between a cored and a cuspy one, and NGC 5963 having a cuspy profile. After calculating baryon fraction, which is $\simeq 0.1$ for the three galaxies, we fitted the rotation curves changing the value of angular momentum. NGC 2976, has an higher value of ordered angular momentum ($\lambda \simeq 0.04$) with respect to NGC 5949 ($\lambda \simeq 0.025$) and in the case of NGC 5963 the very steep profile can be obtained with a low value of $\lambda$ ($\lambda \simeq 0.02$) and also decreasing the value of the random angular momentum. In the case of NGC 2976 tidal interaction with M81 could have also influenced the inner part of the density profile.

We study the pseudo phase-space density, $\rho(r)/\sigma^3(r)$, of $\Lambda$CDM dark matter halos with and without baryons (baryons+DM, and pure DM), by using the model introduced in Del Popolo (2009), which takes into account the effect of dynamical friction, ordered and random angular momentum, baryons adiabatic contraction and dark matter baryons interplay. We examine the radial dependence of $\rho(r)/\sigma^3(r)$ over 9 orders of magnitude in radius for structures on galactic and cluster of galaxies scales. We find that $\rho(r)/\sigma^3(r)$ is approximately a power-law only in the range of halo radius resolved by current simulations (down to 0.1% of the virial radius) while it has a non-power law behavior below the quoted scale, with inner profiles changing with mass. The non-power-law behavior is more evident for halos constituted both of dark matter and baryons while halos constituted just of dark matter and with angular momentum chosen to reproduce a Navarro-Frenk-White (NFW) density profile, are characterized by an approximately power-law behavior. The results of the present paper lead to conclude that density profiles of the NFW type are compatible with a power-law behavior of $\rho(r)/\sigma^3(r)$, while those flattening to the halo center, like those found in Del Popolo (2009) or the Einasto profile, or the Burkert profile, cannot produce radial profile of the pseudo-phase-space density that are power-laws at all radii. The results argue against universality of the pseudo phase-space density and as a consequence argue against universality of density profiles constituted by dark matter and baryons as also discussed in Del Popolo (2009).

Scaling relations among dark matter (DM) and stellar quantities are a valuable tool to constrain formation scenarios and the evolution of galactic structures. However, most of the DM properties are actually not directly measured, but derived through model dependent mass mapping procedures. It is therefore crucial to adopt theoretically and observationally well founded models. We use here an updated version of the secondary infall model (SIM) to predict the halo density profile, taking into account the effects of angular momentum, dissipative friction and baryons collapse. The resulting family of halo profiles depends on one parameter only, the virial mass, and nicely fits the projected mass and aperture velocity dispersion of a sample of intermediate redshift lens galaxies. We derive DM related quantities (namely the column density and the Newtonian acceleration) and investigate their correlations with stellar mass, luminosity, effective radius and virial mass.

Considering galaxies as self - gravitating systems of many collisionless particles allows to use methods of statistical mechanics inferring the distribution function of these stellar systems. Actually, the long range nature of the gravitational force contrasts with the underlying assumptions of Boltzmann statistics where the interactions among particles are assumed to be short ranged. A particular generalization of the classical Boltzmann formalism is available within the nonextensive context of Tsallis q -statistics, subject to non -additivity of the entropies of sub - systems. Assuming stationarity and isotropy in the velocity space, it is possible solving the generalized collsionless Boltzmann equation to derive the galaxy distribution function and density profile. We present a particular set of nonextensive models and investigate their dynamical and observable properties. As a test of the viability of this generalized context, we fit the rotation curve of M33 showing that the proposed approach leads to dark matter haloes in excellent agreement with the observed data.

We use the secondary infall model described in Del Popolo (2009), which takes into account the effect of dynamical friction, ordered and random angular momentum, baryons adiabatic contraction and dark matter baryons interplay, to study how in- ner slopes of relaxed LCDM dark matter (DM) halos with and without baryons (baryons+DM, and pure DM) depend on redshift and on halo mass. We apply the quoted method to structures on galactic scales and clusters of galaxies scales. We find that the inner logarithmic density slope, of dark matter halos with baryons has a significant dependence on halo mass and redshift with slopes ranging from 0 for dwarf galaxies to 0.4 for objects of M = 10^13M_solar and 0.94 for M = 10^15M_solar clusters of galaxies. Structures slopes increase with increasing redshift and this trend reduces going from galaxies to clusters. In the case of density profiles constituted just of dark matter the mass and redshift dependence of slope is very slight. In this last case, we used the Merrit et al. (2006) analysis who compared N-body density profiles with various parametric models finding systematic variation in profile shape with halo mass. This last analysis suggests that the galaxy-sized halos obtained with our model have a different shape parameter, i.e. a different mass distribution, than the cluster-sized halos, obtained with the same model. The results of the present paper argue against universality of density profiles constituted by dark matter and baryons and confirm claims of a systematic variation in profile shape with halo mass, for dark matter halos.

We show how to improve constraints on \Omega_m, \sigma_8, and the dark-energy equation-of-state parameter, w, obtained by Mantz et al. (2008) from measurements of the X-ray luminosity function of galaxy clusters, namely MACS, the local BCS and the REFLEX galaxy cluster samples with luminosities L> 3 \times 10^{44} erg/s in the 0.1--2.4 keV band. To this aim, we use Tinker et al. (2008) mass function instead of Jenkins et al. (2001) and the M-L relationship obtained from Del Popolo (2002) and Del Popolo et al. (2005). Using the same methods and priors of Mantz et al. (2008), we find, for a \Lambda$CDM universe, \Omega_m=0.28^{+0.05}_{-0.04} and \sigma_8=0.78^{+0.04}_{-0.05}$ while the result of Mantz et al. (2008) gives less tight constraints $\Omega_m=0.28^{+0.11}_{-0.07}$ and \sigma_8=0.78^{+0.11}_{-0.13}. In the case of a wCDM model, we find \Omega_m=0.27^{+0.07}_{-0.06}, $\sigma_8=0.81^{+0.05}_{-0.06}$ and $w=-1.3^{+0.3}_{-0.4}$, while in Mantz et al. (2008) they are again less tight \Omega_m=0.24^{+0.15}_{-0.07}, \sigma_8=0.85^{+0.13}_{-0.20} and w=-1.4^{+0.4}_{-0.7}. Combining the XLF analysis with the f_{gas}+CMB+SNIa data set results in the constraint \Omega_m=0.269 \pm 0.012, \sigma_8=0.81 \pm 0.021 and w=-1.02 \pm 0.04, to be compared with Mantz et al. (2008), \Omega_m=0.269 \pm 0.016, \sigma_8=0.82 \pm 0.03 and w=-1.02 \pm 0.06. The tightness of the last constraints obtained by Mantz et al. (2008), are fundamentally due to the tightness of the $f_{gas}$+CMB+SNIa constraints and not to their XLF analysis. Our findings, consistent with w=-1, lend additional support to the cosmological-constant model.

Current theories on planetary formation establish that giant planet formation should be contextual to their quick migration towards the central star due to the protoplanets-disc interactions on a timescale of the order of $10^5$ years, for objects of nearly 10 terrestrial masses. Such a timescale should be smaller by an order of magnitude than that of gas accretion onto the protoplanet during the hierarchical growing-up of protoplanets by collisions with other minor objects. These arguments have recently been analysed using N-body and/or fluid-dynamics codes or a mixing of them. In this work, inviscid 2D simulations are performed, using the SPH method, to study the migration of one protoplanet, to evaluate the effectiveness of the accretion disc in the protoplanet dragging towards the central star, as a function of the mass of the planet itself, of disc tangential kinematics. To this purpose, the SPH scheme is considered suitable to study the roles of turbulence, kinematic and boundary conditions, due to its intrinsic advective turbulence, especially in 2D and in 3D codes. Simulations are performed both in disc sub-Keplerian and in Keplerian kinematic conditions as a parameter study of protoplanetary migration if moderate and consistent deviations from Keplerian Kinematics occur. Our results show migration times of a few orbital periods for Earth-like planets in sub-Keplerian conditions, while for Jupiter-like planets estimates give that about $10^4$ orbital periods are needed to half the orbital size. Timescales of planet migration are strongly dependent on the relative position of the planet with respect to the shock region near the centrifugal barrier of the disc flow.

We study the cusp/core problem using a secondary infall model (SIM) that takes into account the effect of ordered and random angular momentum, dynamical friction and baryons adiabatic contraction. The model is applied to structures on galactic scales (normal and dwarfs spiral galaxies) and on clusters of galaxies scales. Our analysis suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The slope of density profile of inner haloes flattens with decreasing halo mass and the profile is well approximated by a Burkert's profile. In order to obtain the NFW profile, starting from the profiles obtained from our model, the magnitude of angular momentum and dynamical friction must be reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from Gentile et al. (2004) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. The time evolution of the density profile of a galaxy of $10^8-10^9 M_{\odot}$ shows that after a transient steepening, due to the adiabatic contraction, the density profile flattens to $\alpha \simeq 0$. On cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $\alpha \simeq 0.6$ for a cluster of $\simeq 10^{14} M_{\odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.

We present constraints on the values of $\Omega_m$, $n$, $\sigma_8$, obtained from measurements of the X-ray luminosity function of galaxy clusters as compiled in EMSS, RDCS and BCS galaxy cluster samples. The values obtained $\Omega_m=0.25^{+0.15}_{-0.05}$, $n=-1^{+0.05}_{-0.05}$, and $\sigma_8=0.8^{+0.1}_{-0.1}$ are in agreement with WMAP data. We then put constraints on the quoted parameters and the dark-energy equation-of-state parameter, $w$, using Chandra observations of 37 clusters with $<z>$=0.55 derived from 400 deg$^2$ ROSAT serendipitous survey and 49 brightest $z \simeq 0.05$ clusters detected in the All-Sky Survey. In the case of $\Lambda$CDM model, we get $\Omega_m=0.25^{+0.1}_{-0.08}$ and $\sigma_8=0.75^{+0.08}_{-0.05}$, while in the case of the $w$CDM model, we get $\Omega_m=0.28^{+0.10}_{-0.10}$, $w=-1.25^{+0.30}_{-0.35}$, and $\sigma_8=0.8^{+0.09}_{-0.09}$. Our results are consistent with those from recent analyses of type Ia supernovae, cosmic microwave background anisotropies, the X-ray gas mass fraction of relaxed galaxy clusters, baryon acoustic oscillations and cosmic shear. The improvement in data quality from EMSS, RDCS, and BCS to Chandra observations leads to an improvement in the constraints even if not of the same entity of the improvement in data.

Early - type galaxies (ETGs) define a narrow strip in the size - mass plane because of the observed correlation between the effective radius $R_{eff}$ and the total stellar mass $M_{\star}$. When expressed in logarithmic units, a linear relation, $\log{R_{eff}} \propto \gamma \log{M_{\star}}$, is indeed observationally found, but the slope $\gamma$ deviates from the canonical $\gamma = 1/2$ value which can be naively predicted for a spherically symmetric isotropic system. We propose here that a transfer of angular momentum to the stellar component induces an anisotropy in the velocity space thus leading to a modified distribution function (DF). Assuming an Osipkov - Merritt like anisotropic DF, we derive an analytic relation between the slope $\gamma$ of the size - mass relation and the slope $\alpha$ of the angular momentum term in the DF. With this simple model, we are then able to recover the observed $\gamma$ value provided $\alpha$ is suitably set. It turns out that an anisotropy profile which is tangential inside $\sim 0.6 r_a$ and radial outside, with $r_a$ the anisotropy radius, is able to reproduce the observed size - mass relation observed for massive ($M_{\star} \ge 2 \times 10^{10} \ h^{-1} {\rm M}_{\odot}$) elliptical galaxies.

In the present paper, we improve the "Extended Secondary Infall Model" (ESIM) of Williams et al. (2004) to obtain further insights on the cusp/core problem. The model takes into account the effect of ordered and random angular momentum, dynamical friction and baryon adiabatic contraction in order to obtain a secondary infall model more close to the collapse reality. The model is applied to structures on galactic scales (normal and dwarf spiral galaxies) and on cluster of galaxies scales. The results obtained suggest that angular momentum and dynamical friction are able, on galactic scales, to overcome the competing effect of adiabatic contraction eliminating the cusp. The NFW profile can be reobtained, in our model only if the system is constituted just by dark matter and the magnitude of angular momentum and dynamical friction are reduced with respect to the values predicted by the model itself. The rotation curves of four LSB galaxies from de Blok & Bosma (2002) are compared to the rotation curves obtained by the model in the present paper obtaining a good fit to the observational data. On scales smaller than $\simeq 10^{11} h^{-1} M_{\odot}$ the slope $\alpha \simeq 0$ and on cluster scales we observe a similar evolution of the dark matter density profile but in this case the density profile slope flattens to $\alpha \simeq 0.6$ for a cluster of $\simeq 10^{14} h^{-1} M_{\odot}$. The total mass profile, differently from that of dark matter, shows a central cusp well fitted by a NFW model.

This paper is to explore a model of the ABS Algorithms dealing with the solution of a class of systems of linear stochastic equations $A\xi=\eta$ when $\eta$ is a $m$-dimensional normal distribution. It is shown that the stepsize $\alpha_i$ is distributed as $N(u_i,\sigma_i)$ (being $u_i$ the expected value of $\alpha_i$ and $\sigma_i$ its variance) and the approximation to the solutions $\xi_{i}$ is distributed as $N_n(U_i,\Sigma_i)$ (being $U_i$ the expected value of $\xi_i$ and $\Sigma_i$ its variance), for this algorithm model.

This paper provides a review of the variants of dark matter which are thought to be fundamental components of the universe and their role in origin and evolution of structures and some new original results concerning improvements to the spherical collapse model. In particular, I show how the spherical collapse model is modified when we take into account dynamical friction and tidal torques.

I study the joint effect of dynamical friction, tidal torques and cosmological constant on clusters of galaxies formation I show that within high-density environments, such as rich clusters of galaxies, both dynamical friction and tidal torques slows down the collapse of low-? peaks producing an observable variation in the time of collapse of the perturbation and, as a consequence, a reduction in the mass bound to the collapsed perturbation Moreover, the delay of the collapse produces a tendency for less dense regions to accrete less mass, with respect to a classical spherical model, inducing a biasing of over-dense regions toward higher mass I show how the threshold of collapse is modified if dynamical friction, tidal torques and a non-zero cosmological constant are taken into account and I use the Extended Press Schecter (EPS) approach to calculate the effects on the mass function Then, I compare the numerical mass function given in Reed et al (2003) with the theoretical mass function obtained in the present paper I show that the barrier obtained in the present paper gives rise to a better description of the mass function evolution with respect to other previous models (Sheth & Tormen 1999, MNRAS, 308, 119 (hereafter ST); Sheth & Tormen 2002, MNRAS, 329, 61 (hereafter ST1))

In this paper I show how the statistics of the gravitational field is changed when the system is characterized by a non-uniform distribution of particles. I show how the distribution functions W(dF/dt) giving the joint probability that a test particle is subject to a force F and an associated rate of change of F given by dF/dt, are modified by inhomogeneity. Then I calculate the first moment of dF/dt to study the effects of inhomogenity on dynamical friction. Finally I test, by N-Body simulations, that the theoretical W(F) and dF/dt describes correctly the experimental data and I find that the stochastic force distribution obtained for the evolved system is in good agreement with theory. Moreover, I find that in an inhomogeneous background the friction force is actually enhanced relative to the homogeneous case.

This paper provides, from one side, a review of the theory of the cosmological mass function from a theoretical point of view, starting from the seminal paper of Press & Shechter (1974) to the last developments (Del Popolo & Gambera (1998, 1999), Sheth & Tormen 1999 (ST), Sheth, Mo & Tormen 2001 (ST1), Jenkins et al. 2001 (J01), Shet & Tormen 2002 (ST2), Del Popolo 2002a, Yagi et al. 2004 (YNY)), and from another side some improvements on the multiplicity function models in literature. ...

We derive the luminosity-temperature relation for clusters of galaxies by means of a modification of the self-similar model to take account of angular momentum acquisition by protostructures and of an external pressure term in the virial theorem. The fundamental result of the model is that gravitational collapse, which takes account of angular momentum acquisition, can explain the non self-similarity in the L-T relation, in disagreement with the largely accepted assumption that heating/cooling processes and similar are fundamental in the originating the non-self similar behavior (shaping) of the L-T relation.

I compare the numerical multiplicity function given in Yahagi, Nagashima & Yoshii (2004) with the theoretical multiplicity function obtained by means of the excursion set model and an improved version of the barrier shape obtained in Del Popolo & Gambera (1998), which implicitly takes account of total angular momentum acquired by the proto-structure during evolution and of a non-zero cosmological constant. I show that the multiplicity function obtained in the present paper, is in better agreement with Yahagi, Nagashima & Yoshii (2004) simulations than other previous models (Sheth & Tormen 1999; Sheth, Mo & Tormen 2001; Sheth & Tormen 2002; Jenkins et al. 2001) and that differently from some previous multiplicity function models (Jenkins et al. 2001; Yahagi, Nagashima & Yoshii 2004) it was obtained from a sound theoretical background.

I calculate the mass function evolution in a $\Lambda$CDM model by means of the excursion set model and an improved version of the barrier shape obtained in Del Popolo & Gambera (1998), which implicitly takes account of the total angular momentum acquired by the proto-structure during evolution and of a non-zero cosmological constant. I compare the result with Reed et al. (2003), who used a high resolution $\Lambda$CDM numerical simulation to calculate the mass function of dark matter haloes down to the scale of dwarf galaxies, back to a redshift of fifteen. I show that the mass function obtained in the present paper, gives similar predictions to the Sheth & Tormen mass function but it does not show the overprediction of extremely rare objects shown by the Sheth and Tormen mass function. The results confirm previous findings that the simulated halo mass function can be described solely by the variance of the mass distribution, and thus has no explicit redshift dependence.

The constraints obtained by several authors (Eke et al. 1998; Henry 2000) on the estimated values of $\Omega_{\rm m}$, $n$ and $\sigma_8$ are revisited in the light of recent theoretical developments: 1) new theoretical mass functions; 2) a more accurate mass-temperature relation, also determined for arbitrary $\Omega_{\rm m}$ and $\Omega_{\rm \Lambda}$. I re-derive the X-ray Temperature Function (XTF), similarly to Henry (2000) and Eke et al. (1999), re-obtaining the constraints on $\Omega_{\rm m}$, $n$, $\sigma_8$. The result shows that changes in the mass function and M-T relation produces an increase in $\Omega_{\rm m}$ of $ \simeq 20%$ and similar results in $\sigma_8$ and $n$.

In this paper, I revisit the constraints obtained by several authors (Reichart et al. 1999; Eke et al. 1998; Henry 2000) on the estimated values of $\Omega_{\rm m}$, $n$ and $\sigma_8$ in the light of recent theoretical developments: 1) new theoretical mass functions (Sheth & Tormen 1999, Sheth, Mo & Tormen 2001, Del Popolo 2002b); 2) a more accurate mass-temperature relation, also determined for arbitrary $\Omega_{\rm m}$ and $\Omega_{\rm \Lambda}$ (Del Popolo 2002a). Firstly, using the quoted improvements, I re-derive an expression for the X-ray Luminosity Function (XLF), similarly to Reichart et al. (1999), and then I get some constraints to $\Omega_{\rm m}$ and $n$, by using the {\it ROSAT} BCS and {\rm EMSS} samples and maximum-likelihood analysis. Then I re-derive the X-ray Temperature Function (XTF), similarly to Henry (2000), re-obtaining the constraints on $\Omega_{\rm m}$, $n$, $\sigma_8$. Both in the case of the XLF and XTF, the changes in the mass function and M-T relation produces an increase in $\Omega_{\rm m}$ of $ \simeq 20%$ and similar results in $\sigma_8$ and $n$.

We use the model for the migration of planets introduced in Del Popolo, Yesilyurt & Ercan (2003) to calculate the observed mass and semimajor axis distribution of extra-solar planets. The assumption that the surface density in planetesimals is proportional to that of gas is relaxed, and in order to describe disc evolution we use a method which, using a series of simplifying assumptions, is able to simultaneously follow the evolution of gas and solid particles for up to $10^7 {\rm yr}$. The distribution of planetesimals obtained after $10^7 {\rm yr}$ is used to study the migration rate of a giant planet through the model of this paper. The disk and migration models are used to calculate the distribution of planets as function of mass and semimajor axis. The results show that the model can give a reasonable prediction of planets' semi-major axes and mass distribution. In particular there is a pile-up of planets at $a \simeq 0.05$ AU, a minimum near 0.3 AU, indicating a paucity of planets at that distance, and a rise for semi-major axes larger than 0.3 AU, out to 3 AU. The semi-major axis distribution shows that the more massive planets (typically, masses larger than $4 M_{\rm J}$) form preferentially in the outer regions and do not migrate much. Intermediate-mass objects migrate more easily whatever the distance they form, and that the lighter planets (masses from sub-Saturnian to Jovian) migrate easily.

A luminosity-temperature relation for clusters of galaxies is derived. The two models used, take into account the angular momentum acquisition by the proto-structures during their expansion and collapse. The first one is a modification of the self-similar model (SSM) while the second one is a modification of the Punctuated Equilibria Model (Cavaliere et al. 1999). In both models the mass-temperature relation (M-T) used is based on the calculations of Del Popolo (2002b). We show that the above models lead, in X-rays, to a luminosity-temperature relation that scales as L \propto T^5, at scale of groups, flattening to L \propto T^3 for rich clusters and converging to L \propto T^2 at higher temperatures. However a fundamental result of our paper is that the non-similarity in the L-T relation, can be explained by a simple model that takes into account the amount of the angular momentum of a proto-structure. This result is in disagreement with the widely accepted idea that the above non-similarity is due to non-gravitating processes as those of heating/cooling.

In this paper, we discuss improvements of the Suto et al. (2000) model, in the light of recent theoretical developments (new theoretical mass functions, a more accurate mass-temperature relation and an improved bias model) to predict the clustering properties of galaxy clusters and to obtain constraints on cosmological parameters. We re-derive the two-point correlation function of clusters of galaxies for OCDM and LambdaCDM cosmological models, and we compare these results with the observed spatial correlation function for clusters in RASS1 (ROSAT All-Sky Survey 1), and in XBACs (X-RAY Brighest Abell-Type) samples. The comparison shows that the best agreement is obtained for the LambdaCDM model with Omega=0.3. The values of the correlation length obtained, (r_\simeq 28.2 \pm 5.2 \rm h^{-1}} Mpc for LambdaCDM), are larger than those found in the literature and comparable with the results found in Borgani, Plionis & Kolokotronis (1999). (REST IN THE PAPER ABSTRACT)

We have used merger trees realizations to study the formation of dark matter haloes. The construction of merger-trees is based on three different pictures about the formation of structures in the Universe. These pictures include: the spherical collapse (SC), the ellipsoidal collapse (EC) and the non-radial collapse (NR). The reliability of merger-trees has been examined comparing their predictions related to the distribution of the number of progenitors, as well as the distribution of formation times, with the predictions of analytical relations. The comparison yields a very satisfactory agreement. Subsequently, >.........

In a planetary or satellite system, idealized as n small bodies in initially coplanar, concentric orbits around a large central body, obeying Newtonian point-particle mechanics, resonant perturbations will cause dynamical evolution of the orbital radii except under highly specific mutual relationships, here derived analytically.

In this paper, I revisit the constraints obtained by several authors (Reichart et al. 1999; Eke et al. 1998; Henry 2000) on the estimated values of Omegam, n and sigma_8 in the light of recent theoretical developments: 1) new theoretical mass functions (Sheth & Tormen 1999, Sheth, Mo & Tormen 2001, Del Popolo 2002b); 2) a more accurate mass-temperature relation, also determined for arbitrary Omega_m and Omega_Lambda (Del Popolo 2002a).

We consider the non-radial motions, originating in the outskirts of clusters of galaxies and we show how it may reduce the discrepancy between Cold Dark Matter (CDM) predicted X-ray temperature distribution function of clusters of galaxies and that really observed, so resolving the problem of the X-ray clusters abundance over-production predicted by the CDM model. We construct the X-ray temperature distribution function using Press-Schechter's (1974) theory and Evrard's (1990) prescriptions for the mass-temperature relation, taking also account of the non-radial motions originating from the gravitational interaction of the quadrupole moment of the protocluster. Then, we compare our X-ray temperature distribution function with the tidal field of the matter of the neighboring protostructures. We find that the model produces a reasonable cluster temperature distribution.

We calculate the density profiles of virialized halos in the case of structure evolving hierarchically from a scale-free Gaussian delta-field having a power spectrum $P(k)\propto k^n$ in a Omega=1 Universe; we suppose that the initial density contrast profile around local maxima is given by the mean peak profile introduced by Bardeen et al. (1986 hereafter BBKS). We show both that the density profiles are not power-laws but have a logarithmic slope that increases from the inner halo to its outer parts and for n \geq -1 are well approximated by Navarro et al. (1995, 1996, 1997) profile and the radius a, at which the slope alpha=-2, is a function of the mass of the halo and of the spectral index n.

In this paper, I revisit the constraints obtained by several authors (Reichart et al. 1999; Eke et al. 1998; Henry 2000) on the estimated values of Omega_m, n and sigma_8 in the light of recent theoretical developments: 1) new theoretical mass functions (Sheth & Tormen 1999, Sheth, Mo & Tormen 1999, Del Popolo 2002b); 2) a more accurate mass-temperature relation, also determined for arbitrary Omega_m and Omega_{\Lambda} (Voit 2000, Pierpaoli et al. 2001, Del Popolo 2002a). Firstly, using the quoted improvements, I re-derive an expression for the X-ray Luminosity Function (XLF), similarly to Reichart et al. (1999), and then I get some constraints to \Omega_m and n, by using the ROSAT BCS and EMSS samples and maximum-likelihood analysis. Then I re-derive the X-ray Temperature Function (XTF), similarly to Henry (2000) and Eke et al. (1999), re-obtaining the constraints on Omega_m, n, sigma_8. Both in the case of the XLF and XTF, the changes in the mass function and M-T relation produces an increase in Omega_m of \simeq 20% and similar results in sigma_8 and n.

In the current paper, we further improved the model for the migration of planets introduced in Del Popolo et al. (2001) and extended to time-dependent planetesimal accretion disks in Del Popolo and Eksi (2002). In the current study, the assumption of Del Popolo and Eksi (2002), that the surface density in planetesimals is proportional to that of gas, is released. In order to obtain the evolution of planetesimal density, we use a method developed in Stepinski and Valageas (1997) which is able to simultaneously follow the evolution of gas and solid particles for up to 10^7 yrs. Then, the disk model is coupled to migration model introduced in Del Popolo et al. (2001) in order to obtain the migration rate of the planet in the planetesimal. We find that the properties of solids known to exist in protoplanetary systems, together with reasonable density profiles for the disk, lead to a characteristic radius in the range 0.03-0.2 AU for the final semi-major axis of the giant planet.

In this paper, I test by numerical simulations the results of Del Popolo & Gambera (1998),dealing with the extension of Chandrasekhar and von Neumann's analysis of the statistics of the gravitational field to systems in which particles (e.g., stars, galaxies) are inhomogeneously distributed. The paper is an extension of that of Ahmad & Cohen (1974), in which the authors tested some results of the stochastic theory of dynamical friction developed by Chandrasekhar & von Neumann (1943) in the case of homogeneous gravitational systems. It is also a continuation of the work developed in Del Popolo (1996a,b), which extended the results of Ahmad & Cohen (1973), (dealing with the study of the probability distribution of the stochastic force in homogeneous gravitational systems) to inhomogeneous gravitational systems. Similarly to what was done by Ahmad & Cohen (1974) in the case of homogeneous systems, I test, by means of the evolution of an inhomogeneous system of particles, that the theoretical rate of force fluctuation d F/dt describes correctly the experimental one, I find that the stochastic force distribution obtained for the evolved system is in good agreement with the Del Popolo & Gambera (1998) theory. Moreover, in an inhomogeneous background the friction force is actually enhanced relative to the homogeneous case.

In the current paper, we further develop the model for the migration of planets introduced in Del Popolo et al. (2001) and extended to time-dependent accretion discs in Del Popolo and Eksi (2002). We use a method developed by Stepinski and Valageas (1996, 1997), that is able to simultaneously follow the evolution of gas and solid particles for up to $10^7 {\rm yr}$. The disc model is coupled to the migration model introduced in Del Popolo et al. (2001) in order to obtain the migration rate of the planet in the planetesimal disc. We find that in the case of discs having total mass of $10^{-3}-0.1 M_{\odot}$, and $0.1<\alpha<0.0001$, planets can migrate inward a large distance while if $M<10^{-3} M_{\odot}$ the planets remain almost in their initial position for $0.1<\alpha<0.01$ and only in the case $\alpha<0.001$ the planets move to a minimum value of orbital radius of $\simeq 2 {\rm AU}$. The model gives a good description of the observed distribution of planets in the period range 0-20 days.

In this paper, we further develop the model for the migration of planets introduced in Del Popolo, Gambera and Ercan, and extended to time-dependent planetesimal accretion disks in Del Popolo and Eksi. More precisely, the assumption of Del Popolo and Eksi that the surface density in planetesimals is proportional to that of gas is released. Indeed, the evolution of the radial distribution of solids is governed by many processes: gas-solid coupling, coagulation, sedimentation, evaporation/condensation, so that the distribution of planetesimals emerging from a turbulent disk does not necessarily reflect that of gas. In order to describe this evolution we use a method developed in Stepinski and Valageas which, using a series of simplifying assumptions, is able to simultaneously follow the evolution of gas and solid particles for up to $10^7 {\rm yr}$. Then, the distribution of planetesimals obtained after $10^7 {\rm yr}$ is used to study the migration rate of a giant planet through the migration model introduced in \cite{DP1}. This allows us to investigate the dependence of the migration rate on the disk mass, on its time evolution and on the value of the dimensionless viscosity parameter alpha. We find that in the case of disks having a total mass of $10^{-3}-10^{-1} M_{\odot}$, and $10^{-4}<\alpha<10^{-1}$, planets can migrate inward over a large distance while if $M_{\rm d}<10^{-3} M_{\odot}$ the planets remain almost at their initial position for $\alpha>10^{-3}$ and only in the case $\alpha<10^{-3}$ the planets move to a minimum value of orbital radius of $\simeq 2 {\rm AU}$. Moreover, the observed distribution of planets in the period range 0-20 days can be easily obtained from our model.

ABS methods are a large class of methods, based upon the Egervary rank reducing algebraic process, first introduced in 1984 by Abaffy, Broyden and Spedicato for solving linear algebraic systems, and later extended to nonlinear algebraic equations, to optimization problems and other fields; software based upon ABS methods is now under development. Current ABS literature consists of about 400 papers. ABS methods provide a unification of several classes of classical algorithms and more efficient new solvers for a number of problems. In this paper we review ABS methods for linear systems and optimization, from both the point of view of theory and the numerical performance of ABSPACK.

This paper provides a review of the variants of dark matter (CDM, HDM) which are thought to be fundamental components of the universe and their role in origin and evolution of structures.

We present some non-standard approaches to the N-body problem in an attempt to overcome its epistemological limits. We provide, in a preliminary way, not-ordinary insights and breakthroughs.

We study the effect of the interaction between small scale clumps having velocities larger than the average escape velocity v_{\rm esc} and bound clumps (having v<v_{\rm esc}), on clusters of galaxies evolution. We find that the collapse of the bound clumps is accelerated with respect to the homogeneous case (Gunn & Gott's model - Gunn & Gott 1972). A comparison between the acceleration produced by this effect to the slowing down effect studied in Del Popolo & Gambera (1998) shows that the magnitude of the slowing down effect is larger than the acceleration produced by the effect studyed in this paper, only in the outskirts of the cluster.

In this paper, we study the role of shear fields on the evolution of density perturbations embedded in a Friedmann flat background universe, by studying the evolution of a homogeneous ellipsoid model. In this context, we show that while the effect of the shear is that of increasing the growth rate of the density contrast of a mass element, the angular momentum acquired by the ellipsoid has the right magnitude to counterbalance the shear. Finally, our result show that initial asphericities and tidal interaction induce a slowing down of the collapse after the system has broken away from the general expansion, in perfect agreement with the previrialization conjecture (Peebles & Groth 1976; Davis & Peebles 1977).

In this paper, I use the extension of the excursion set model of Sheth & Tormen (2002) and the barrier shape obtained in Del Popolo & Gambera (1998) to calculate the unconditional halo mass function, and the conditional mass function in several cosmological models. I show that the barrier obtained in Del Popolo & Gambera (1998), which takes account of tidal interaction between proto-haloes, is a better description of the mass functions than the spherical collapse and is in good agreement with numerical simulations (Tozzi & Governato 1998, and Governato et al. 1999). The results are also in good agreement with those obtained by Sheth & Tormen (2002), only slight differences are observed expecially at the low mass end. I moreover calculate, and compare with simulations, the temperature function obtained by means of the mass functions previously calculated and also using an improved version of the M-T relation, which accounts for the fact that massive clusters accrete matter quasi-continuously, and finally taking account of the tidal interaction with neighboring clusters. Even in this case the discrepancy between the Press-Schecter predictions and simulations is considerably reduced.

I derive the mass-temperature relation and its time evolution for clusters of galaxies in different cosmologies by means of two different models. The first one is a modification and improvement of a model by Del Popolo & Gambera(1999), namely based upon a modification of the top-hat model in order to take account of angular momentum acquisition by protostructures and of an external pressure term in the virial theorem. The second one is based on the merging-halo formalism of Lacey & Cole (1993), accounting for the fact that massive clusters accrete matter quasi-continuously, and is an improvement of a model proposed by Voit (2000). The final result is that, in both models, the M-T relation shows a break at T \sim 3-4 keV. The behavior of the M-T relation is as usual, M \propto T^{3/2}, at the high mass end, and M \propto T^{\gamma}, with a value of \gamma>3/2 depending on the chosen cosmology. The evolution of the M-T relation, for a given M_{\rm vir}, is more modest both in flat and open universes in comparison to previous estimate found in literature, even more modest than what found by V2000. Moreover the time evolution is more rapid in models with L=0 than in models in which the angular momentum acquisition by protostructures is taken into account (L \neq 0). The effect of a non-zero cosmological constant is that of slightly increasing the evolution of the M-T relation with respect to open models with L \neq 0. The evolution is more rapid for larger values (in absolute value) of the spectral index, n. The comparison of the mass-temperature relation with the data by Finoguenov, Reiprich & Bohringer (2001), shows that the FRB data is able to rule out very low Omega_0 models (<0.3), particularly in the open case, and that better fit are obtained by LambdaCDM models and by CDM models with Omega_0>0.3.

In the Cold Dark Matter (hereafter CDM) scenario even isolated density peaks contain a high fraction of small scale clumps having velocities larger than the average escape velocity from the structure. These clumps populate protoclusters, especially in the peripheral regions, r > R_f (where R_f is the filtering scale). During the cluster collapse and the subsequent secondary infall, collapsing or infalling clumps (having v<v_esc) interact with the quoted unbound clumps (or high speed clumps, as we also call them) having v > v_esc. We study the interaction between these two kind of clumps by means of the impulse approximation and we find that the collapse of bound clumps is accelerated with respect to the homogeneous case (Gunn & Gott's model). The acceleration of the collapse increases with decreasing height of the peak, \nu. We finally compare the acceleration produced by this effect to the slowing down effect produced by the gravitational interaction of the quadrupole moment of the system with the tidal field of the matter of the neighboring proto-clusters studied in Del Popolo & Gambera. We find that the magnitude of the slowing down effect is larger than the acceleration produced by the effect studied in this paper, only in the outskirts of the cluster. We want to stress that the one which we study in this paper is also present in an isolated protocluster, being produced by the interaction of the collapsing clumps with the unbound substructure internal to the collapsing clumps itself while that studied in Del Popolo & Gambera (1998) is produced by substructure external to the density peak.

I study the role of shear fields by using an analytical approximate solution for the equations of motion of homogeneous ellipsoids embedded in a homogeneous background. The equations of motion of a homogeneous ellipsoid (Icke 1973; White & Silk 1979 (hereafter WS)) are modified in order to take account of the tidal field as done in Watanabe 1993 and then are integrated analytically, similarly to what done in WS. The comparison of the analytical solution with numerical simulations shows that it is a good approximation to the numerical one. This solution is used to study the evolution of the configuration ofthe ellipsoids, to calculate the evolution of the density contrast and that of the axial peculiar velocity of the ellipsoids for several values of the amplitude of the external tidal field, and compared again with numerical simulations. In order to calculate the evolution of the density contrast at turn-around and collapse velocity at the epoch of collapse, as a function of the ratio of the initial value of the semi-axes, I use the previously obtained approximate solution to modify the analytical model proposed by Barrow & Silk (1981) for the ellipsoids evolution in the non-linear regime. The density contrast at turn-around and the collapse velocity are found to be reduced with respect to that found by means of the spherical model. The reduction increases with increasing strength of the external tidal field and with increasing initial asymmetry of the ellipsoids. These last calculations are also compared with numerical solutions and they are again in good agreement with the numerical ones.

In this paper, we further develop the model for the migration of planets introduced in Del Popolo et al. (2001). We first model the protoplanetary nebula as a time-dependent accretion disc and find self-similar solutions to the equations of the accretion disc that give to us explicit formulas for the spatial structure and the temporal evolution of the nebula. These equations are then used to obtain the migration rate of the planet in the planetesimal disc and to study how the migration rate depends on the disc mass, on its time evolution and on some values of the dimensionless viscosity parameter alpha. We find that planets that are embedded in planetesimal discs, having total mass of 10^{-4}-0.1 M_{\odot}, can migrate inward a large distance for low values of alpha (e.g., alpha \simeq 10^{-3}-10^{-2}) and/or large disc mass and can survive only if the inner disc is truncated or because of tidal interaction with the star. Orbits with larger $a$ are obtained for smaller value of the disc mass and/or for larger values of alpha. This model may explain several orbital features of the recently discovered giant planets orbiting nearby stars.

In this paper we extend Chandrasekhar and von Neumann's analysis of the statistics of the gravitational field to systems in which particles (e.g. stars, galaxies) are not homogeneously distributed. We derive a distribution function W(F,d F/dt) giving the joint probability that a test particle is subject to a force F and an associated rate of change of F given by d F/dt. We calculate the first moment of d F/dt to study the effects of inhomogenity on dynamical friction.

We study the non-linear evolution of a dust ellipsoid,embedded in a Friedmann flat background universe, in order to determine the evolution of the density of the ellipsoid as the perturbation to it related detaches from general expansion and begins to collapse. We show that while the growth rate of the density contrast of a mass element is enhanced by the shear in agreement with Hoffman (1986a), the angular momentum acquired by the ellipsoid has the right magnitude to counterbalance the effect of the shear. The result confirms the previrialization conjecture (Peebles & Groth 1976; Davis & Peebles 1977; Peebles 1990) by showing that initial asphericities and tidal interactions begin to slow down the collapse after the system has broken away from the general expansion.

Numerical experiments are performed in order to study the performance of ABS codes in solving non-linear systems of equations.

We associated, in a classical formulation of "strong gravity", hadron constituents with suitable stationary, axisymmetric solutions of some new Einstein-type equations supposed to describe the strong field inside hadrons. These new equations can be obtained by the Einstein equations with cosmological term Lambda. As a consequence, Lambda and the masses M result in our theory to be scaled up, and transformed into a "hadronic constant" and into "strong masses", respectively. Due to the unusual range of Lambda and M values considered, we met a series of solutions of the Kerr-Newman-de Sitter (hereafter KNdS) type with rather interesting properties. The requirement that those solutions be stable, i.e., that their temperature (or surface gravity) be vanishingly small, implies the coincidence of at least two of their (in general, three) horizons. Imposing the stability condition of a certain horizon does yield (once chosen the values of J, q and Lambda) mass and radius of the associated black-hole (hereafter BH). In the case of ordinary Einstein equations and for stable BHs of the KNdS type, we get in particular Regge-like (hereafter RL) relations among mass M, angular momentum J, charge q and cosmological constant Lambda; which did not receive enough attention in the previous literature. Besides, we show some particular and interesting cases of these relations. Another interesting point is that, with few exceptions, all such relations (among M, J, q, Lambda) lead to solutions that can be regarded as (stable) cosmological models.

We study the joint effect of tidal torques and dynamical friction on the collapse of density peaks solving numerically the equations of motion of a shell of barionic matter falling into the central regions of a cluster of galaxies. We calculate the evolution of the expansion parameter, a(t), of the perturbation using a coefficient of dynamical friction eta_{cl} obtained from a clustered system and taking into account the gravitational interaction of the quadrupole moment of the system with the tidal field of the matter of the neighboring proto-galaxies. We show that within high-density environments, such as rich clusters of galaxies, tidal torques and dynamical friction slow down the collapse of low-nu peaks producing an observable variation of the parameter of expansion of the shell. As a consequence a bias of dynamical nature arises because high-density peaks preferentially collapse to form halos within which visible objects eventually will condense. For a standard Cold Dark Matter model this dynamical bias can account for a substantial part of the total bias required by observations on cluster scales.

We study the effect of non-radial motions on the mass function, the VDF and on the shape of clusters of galaxies using the model introduced in Del Popolo & Gambera (1998a,b; 1999). The mass function of clusters, obtained using the quoted model, is compared with the statistical data by Bahcall & Cen (1992) and Girardi et al. (1998), while the VDF is compared with the CfA data by Zabludoff et al. (1993) for local clusters. In both cases the model predictions are in good agreement with the observational data showing once more how non-radial motions can reduce many of the discrepancies between CDM model predictions and observational data. Besides we study the effect of non-radial motions on the intrinsic shape of clusters of galaxies showing that non-radial motions produce clusters less elongated with respect to CDM model.

In this paper I extend the results of Ahmad & Cohen (1973), regarding the study of the probability distribution of the stochastic force in homogeneous gravitational systems, to inhomogeneous gravitational ones. To this aim, I study the stochastic force distribution using N-body realizations of Plummer's spherically symmetric models. I find that the stochastic force distribution obtained for the evolved system is in good agreement with Kandrup's (1980) theory of stochastic force in inhomogeneous systems. Correlation effects that arise during the evolution of the system of particles are well described by Antonuccio-Delogu & Atrio-Barandela's (1992) theory.

This report is devoted to some C++ codes implementing the implicit LU class algorithms for solving linear determined, and undetermined systems with $n$ variables and $m$ equations. A main program used in part of the numerical test is given in the last section.

This paper is to explore a model of the ABS Algorithms for dealing with a class of systems of linear stochastic equations A xi=eta satisfying eta sim N_m(v, I_{m}). It is shown that the iteration step alpha_{i} is N(V,\pi) and approximation solutions is xi_{i} \sim N_n(U,\Sigma) for this algorithm model. And some properties of (V,\pi)$ and $(U,\Sigma) are given.

In this paper I study the probability distribution of the gravitational force in gravitational systems through numerical experiments. I show that Kandrup's (1980) and Antonuccio-Delogu & Atrio-Barandela's (1992) theories describe correctly the stochastic force probability distribution respectively in inhomogeneous and clustered systems. I find equations for the probability distribution of stochastic forces in finite systems, both homogeneous and clustered, which I use to compare the theoretical predictions with Montecarlo simulations of spherically symmetric systems. The agreement between theoretical predictions and simulations proves to be quite satisfactory.

We present a review and bibliography of the main results obtained during a research on ABS (Abaffy, Broyden, Spedicato) methods.

The results of computational experiments with ABS algorithms for KKT linear systems are reported.

The results of computational experiments with ABS algorithms for overdetermined linear systems are reported.

We present the main results obtained during a research on ABS methods in the framework of the project Analisi Numerica e Matematica Computazionale.

I discuss the effect of non-radial motions on the small-scale pairwise peculiar velocity dispersions of galaxies (PVD) in a CDM model. I calculate the PVD for the SCDM model by means of the refined cosmic virial theorem (CVT) (Suto & Jing 1997) and taking account of non-radial motions by means of Del Popolo & Gambera (1998) model. I compare the results of the present model with the data from Davis & Peebles (1983), the IRAS value at 1 h{-1} Mpc of Fisher et al. (1993) and Marzke et al. (1995). I show that while the SCDM model disagrees with the observed values, as pointed out by several authors (Peebles 1976, 1980; Davis & Peebles 1983; Mo et. al 1993; Jing et al. 1998), taking account of non-radial motions produce smaller values for the PVD. At r <=1 h^{-1} Mpc the result is in agreement with Bartlett & Blanchard (1996) (hereafter BB96). At the light of this last paper, the result may be also read as a strong dependence of the CVT prediction on the model chosen to describe the mass distribution around galaxies, suggesting that the CVT cannot be taken as a direct evidence of a low density universe. Similarly to what shown in Del Popolo & Gambera (1999, 2000) (hereafter DG99, DG00), Del Popolo et al. (1999), the agreement of our model to the observational data is due to a scale dependent bias induced by the presence of non-radial motions. Since the assumptions on which CVT is based have been questioned by several authors (BB96; SJ97), I also calculated the PVD using the redshift distortion in the redshift-space correlation function, \xi_z(r_p,\pi), and I compared it with the PVD measured from the Las Campanas Redshift Survey by J98. The result confirms that non-radial motions influence the PVD making them better agree with observed data.

I study the evolution of haloes density profiles as a function of time in the SCDM and LCDM cosmologies. Following Del Popolo et al. (1999), I calculate the concentration parameter c=r_v/a and study its time evolution. For a given halo mass, I find that c(z) ~1/(1+z) both in the LCDM and SCDM cosmology, in agreement with Bullock et al. (1999) analytic model and N-body simulations. In both models, a(z) is roughly constant. The present model predicts a stronger evolution of c(z) with respect to Navarro et al. (1997) model. Finally I show some consequences of the result on galaxy modelling.

Planets orbiting a planetesimal circumstellar disc can migrate inward from their initial positions because of dynamical friction between planets and planetesimals. The migration rate depends on the disc mass and on its time evolution. Planets that are embedded in long-lived planetesimal discs, having total mass of $10^{-4}-0.01 M_{\odot}$, can migrate inward a large distance and can survive only if the inner disc is truncated or because of tidal interaction with the star. In this case the semi-major axis, a, of the planetary orbit is less than 0.1 AU. Orbits with larger $a$ are obtained for smaller value of the disc mass or for a rapid evolution (depletion) of the disc. This model may explain several of the orbital features of the giant planets that were discovered in last years orbiting nearby stars as well as the metallicity enhancement found in several stars associated with short-period planets.

We study the effect of non-radial motions on the mass function, the velocity dispersion function (hereafter VDF) and on the shape of clusters of galaxies using the model introduced in Del Popolo & Gambera (1998a,b; 1999). The mass function of clusters, obtained using the quoted model, is compared with the statistical data by Bahcall & Cen (1992a,b) and Girardi et al. (1998), while the VDF is compared with the Center for Astrophysics (hereafter CfA) data by Zabludoff et al. (1993) for local clusters and those of Mazure et al. (1996) and Fadda et al. (1996). In both cases the model predictions are in good agreement with the observational data showing once more how non-radial motions can reduce many of the discrepancies between Cold Dark Matter (hereafter CDM) model predictions and observational data. Finally we study the effect of non-radial motions on the intrinsic shape of clusters of galaxies showing that non-radial motions produce clusters less elongated with respect to CDM model in agreement with de Theije et al. (1995, 1997) results.

In this paper we calculate the density profiles of virialized halos both in the case of structure evolving hierarchically from a scale-free Gaussian delta-field having a power spectrum P(k)=A k^n in a Omega=1 Universe and in the case of the CDM model, by using a modified version of Hoffman & Shaham's (1985) (hereafter HS) and Hoffman's (1988) model. We suppose that the initial density contrast profile around local maxima is given by the mean peak profile introduced by Bardeen et al. (1986) (hereafter BBKS), and is not just proportional to the two-point correlation function, as assumed by HS. We show that the density profiles, both for scale-free Universes and the CDM model, are not power-laws but have a logarithmic slope that increases from the inner halo to its outer parts. Both scale-free, for n >=-1, and CDM density profiles are well approximated by Navarro et al. (1995, 1996, 1997) profile. The radius a, at which the slope alpha=-2, is a function of the mass of the halo and in the scale-free models also of the spectral index n.

We study the effect of non-radial motions, originating from the gravitational interaction of the quadrupole moment of a protogalaxy with the tidal field of the matter of the neighboring protostructures, on the angular correlation function of galaxies. We calculate the angular correlation function using a Standard Cold Dark Matter (hereafter SCDM) model (Omega=1, h=0.5, n=1) and we compare it with the angular correlation function of the APM galaxy survey (Maddox et al. 1990; Maddox et al. 1996). We find that taking account of non-radial motions in the calculation of the angular correlation function gives a better agreement of the theoretical prediction of the SCDM model to the observed estimates of large-scale power in the galaxy distribution.

In this paper we study the role of dynamical friction on the evolution of a population of large objects ($m>10^{22}$ g) at heliocentric distances $>70$ AU in the Kuiper Belt. We show that the already flat distribution of these objects must flatten further due to non-spherically symmetric distribution of matter in the Kuiper Belt. Moreover the dynamical drag, produced by dynamical friction, causes objects of masses $\geq 10^{24} g$ to lose angular momentum and to fall through more central regions in a timescale $\approx 10^9 yr$. This mechanism is able to transport inwards objects of the size of Pluto, supposing it was created beyond 50 AU, according to a Stern & Colwell's (1997b) suggestion.

In this paper we show how non-radial motions, originating from the tidal interaction of the irregular mass distribution within and around protoclusters, can solve some of the problems of the CDM model. Firstly the discrepancy between the CDM predicted two-points correlation function of clusters and the observed one. We compare the two-points correlation function, that we obtain taking account of non-radial motions, with that obtained by Sutherland & Efstathiou (1991) from the analysis of Geller & Hucra's (1988) deep redshift survey and with the data points for the APM clusters obtained by Efstathiou et al. (1992). Secondly the problem of the X-ray clusters abundance over-production predicted by the CDM model. In this case we compare the X-ray temperature distribution function, calculated using Press-Schechter theory and Evrard's (1990) prescriptions for the mass-temperature relation, taking also account of the non-radial motions, with Henry & Arnaud (1991) and Edge et al. (1990) X-ray temperature distributions for local clusters. We find that in both cases the model is in good agreement with experimental data. Finally we calculate the bias coefficient using a selection function that takes into account the effects of non-radial motions, and we show that the bias so obtained can account for a substantial part of the total bias required by observations on cluster scales.

In this paper we extend Chandrasekhar and von Neumann's analysis of the statistics of the gravitational field to systems in which particles (e.g. stars, galaxies) are not homogeneously distributed. We derive a distribution function W({F},dF/dt) giving the joint probability that a test particle is subject to a force F and an associated rate of change of F given by dF/dt. We calculate the first moment of dF/dt to study the effects of inhomogenity on dynamical friction.

We show how non-radial motions, originating in the outskirts of clusters of galaxies, may reduce the discrepancy between the Cold Dark Matter (CDM) predicted X-ray temperature distribution function of clusters of galaxies and the observed one and also the discrepancy between the CDM predicted two-point correlation function of clusters of galaxies and that observed. We compare Edge et al. (1990) and Henry & Arnaud (1991) data with the distribution function of X-ray temperature, calculated using Press- Schechter's (1974 - hereafter PS) theory and Evrard's (1990) prescriptions for the mass-temperature relation and taking account of the non-radial motions originating from the gravitational interaction of the quadrupole moment of the protocluster with the tidal field of the matter of the neighboring protostructures. We find that the model produces a reasonable clusters temperature distribution. We compare the two-point cluster correlation function which takes account of the non-radial motions both with that obtained by Sutherland & Efstathiou (1991), from the analysis of Huchra's et al. (1990) deep redshift survey, and with the data points for the Automatic Plate Measuring (APM) clusters, computed by Efstathiou et al. (1992a), showing that non-radial motions reduce the discrepancy between the theoretical and the obsservational correlation function.

We study the effect of tidal torques on the collapse of density peaks through the equations of motion of a shell of barionic matter falling into the central regions of a cluster of galaxies. We calculate the time of collapse of the perturbation taking into account the gravitational interaction of the quadrupole moment of the system with the tidal field of the matter of the neighbouring proto-clusters. We show that within high-density environments, such as rich clusters of galaxies, tidal torques slow down the collapse of low \nu peaks producing an observable variation in the time of collapse of the shell and, as a consequence, a reduction in the mass bound to the collapsed perturbation. Moreover, the delay of the collapse produces a tendency for less dense regions to accrete less mass, with respect to a classical spherical model, inducing a biasing of over-dense regions toward higher mass. Finally we calculate the bias coefficient using a selection function properly defined showing that for a Standard Cold Dark Matter (SCDM) model this bias can account for a substantial part of the total bias required by observations on cluster scales.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter falling into the central regions of a cluster of galaxies, taking into account of the presence of the substructure inducing dynamical friction. The evolution of the expansion parameter a(t) of the perturbation is calculated in spherical systems. The effect of dynamical friction is to reduce the binding radius and the total mass accreted by the central regions. Using a peak density profile given by Bardeen et al. (1986) we show how the binding radius of the perturbation is modified by dinamical friction. We show how dynamical friction modifies the collapse parameter of the perturbation slowing down the collapse.

We solve numerically the equations of motion for the collapse of a shell of baryonic matter, made of galaxies and substructure of $ 10^{6} M_{\odot}- 10^{9} M_{\odot}$, falling into the central regions of a cluster of galaxies taking account of dynamical friction. We calculate the evolution of the expansion parameter, a(t), of the perturbation using a coefficient of dynamical friction, $ \eta_{0}$, calculated for a perturbation in which clustering is absent and a coefficient $ \eta_{cl}$ obtained from a clustered one. The effect of the dynamical friction is to slow down the collapse (V. Antonuccio-Delogu \& S. Colafrancesco 1994, hereafter AC) producing an observable variation of the parameter of expansion of the shell. The effect increases with increasing $ \eta $ and with the increasing of clustering. Finally, we show how the collapse time depends on $ \eta_{0}$ and $ \eta_{cl}$. keywords: {cosmology: theory- galaxies: formation}

We show how the results given by several authors relatively to the mass of a density peak are changed when small scale substructure induced by dynamical friction are taken into account. The peak mass obtained is compared to the result of Peacock \& Heavens (1990) and to the peak mass when dynamical friction is absent to show how these effects conspire to reduce the mass accreted by the peak.

We study the effect of the dynamical friction induced by the presence of substructure on the statistics of the collapse of density peaks. Applying the results of a former paper we show that within high density environments, like rich clusters of galaxies, the collapse of smaller peaks is strongly delayed until very late epochs. A bias of dynamical nature thus naturally arises because high density peaks preferentially collapse For a standard CDM model we find that this dynamical bias can account for a substantial part of the total bias required by observations on cluster scales.