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.

The TOV equation appears as the relativistic counterpart of the classical condition for hydrostatic equilibrium. In the present work we aim at showing that a generalised TOV equation also characterises the equilibrium of models endowed with other symmetries besides spherical. We apply the dual null formalism to spacetimes with two dimensional spherical, planar and hyperbolic symmetries with a perfect fluid as the source. We also assume a Killing vector field orthogonal to the surfaces of symmetry, which gives us static solutions, in the timelike Killing field case, and homogeneous dynamical solutions in the case the Killing field is spacelike. In order to treat equally all the aforementioned cases, we discuss the definition of a quasi-local energy for the spacetimes with planar and hyperbolic foliations, since the Hawking-Hayward definition only applies to compact foliations. After this procedure, we are able to translate our geometrical formalism to the fluid dynamics language in a unified way, to find the generalized TOV equation, for the three cases when the solution is static, and to obtain the evolution equation, for the homogeneous spacetime cases. Remarkably, we show that the static solutions which are not spherically symmetric violate the weak energy condition (WEC). We have also shown that the counterpart of the TOV equation for the spatially homogeneous models is just the familiar equation \r{ho} + P = 0, defining a cosmological constant-type behaviour, both in the hyperbolic and spherical cases. This implies a violation of the strong energy condition in both cases, added to the above mentioned violation of the weak energy condition in the hyperbolic case. We illustrate our unified treatment obtaining analogs of Schwarzschild interior solution, for an incompressible fluid $\rho = \rho_0$ constant.

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.

As the dark sector remains unknown in composition and interaction between dark energy and dark matter stand out as natural, observations of galaxy clusters out of equilibrium abound, opening a promising window on these questions. We continue here the exploration of dark sector interaction detection via clusters virial equilibrium state for all clusters configurations. The dynamics of clusters is evaluated with the Layzer-Irvine equation, a simple model of an interacting dark sector and some simplifying assumptions to obtain the time-dependent part of the virial dynamics. The clusters' data are concentrated in optical weak lensing and X-ray observations that evaluate, respectively, the clusters' mass profiles and temperatures. The global inconsistency of available X-ray data led us to constitute `gold' cluster samples. Through a Bayesian analysis, they are processed to obtain consistent interaction detected up to 3$\sigma$, in compounded interaction strength for 11 clusters at $-0.027 \pm 0.009$ that translate in compounded universal equilibrium virial ratio of $-0.61^{+0.04}_{-0.03}$. The level of detection and inconsistency of X-ray data call for caution, although future instruments promise a clearer detection soon.

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 Birkhoff theorem is a well-known result in general relativity and it is used in many applications. However, its most general version, due to Bona, is almost unknown and presented in a form less accessible to the relativist and cosmologist community. Moreover, many wield it mistakenly as a simple transposition of Newton's iron sphere theorem. In the present work, we propose a modern, dual null, presentation --- useful in many explorations, including black holes --- of the theorem that renders accessible most of the results of Bona's version. In addition, we discuss the fluid contents admissible for the application of the theorem, beyond a vacuum, and we demonstrate how the formalism greatly simplifies solving the dynamical equations and allows one to express the solution as a power expansion in $r$. We present a family of solutions that share the properties predicted by the Birkhoff theorem and discuss the existence of trapped and antitrapped regions. The formalism manifestly shows how the type of region --- trapped or untrapped --- determines the character of the Killing vector.

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.

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.

The observation of galaxy and gas distributions, as well as cosmological simulations in a $\Lambda$CDM Universe, suggests that clusters of galaxies are still accreting mass and are not expected to be in equilibrium. In this work, we investigate the possibility to evaluate the departure from virial equilibrium in order to detect, in that balance, effects from a Dark Matter--Dark Energy interaction. We continue, from previous works, using a simple model of interacting dark sector, the Layzer--Irvine equation for dynamical virial evolution, and employ optical observations in order to obtain the mass profiles through weak lensing and X-ray observations giving the intracluster gas temperatures. Through a Monte Carlo method, we generate, for a set of clusters, measurements of observed virial ratios, interaction strength, rest virial ratio and departure from equilibrium factors. We found a compounded interaction strength of $-1.99^{+2.56}_{-16.00}$, compatible with no interaction, but also a compounded rest virial ratio of $-0.79 \pm 0.13$, which would entail a $2\sigma$ detection. We confirm quantitatively that clusters of galaxies are out of equilibrium but further investigation is needed to constrain a possible interaction in the dark sector.

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 paper we consider spherically symmetric general fluids with heat flux, motivated by causal thermodynamics, and give the appropriate set of conditions that define separating shells defining the divide between expansion and collapse. To do so we add the new requirement that heat flux and its evolution vanish at the separating surface. We extend previous works with a fully nonlinear analysis in the 1+3 splitting, and present gauge-invariant results. The definition of the separating surface is inspired by the conservation of the Misner-Sharp mass, and is obtained by generalizing the Tolman-Oppenheimer-Volkoff equilibrium and turnaround conditions. We emphasize the nonlocal character of these conditions as found in previous works and discuss connections to the phenomena of spacetime cracking and thermal peeling.

We investigate spherically symmetric spacetimes with an anisotropic fluid and discuss the existence and stability of a dividing shell separating expanding and collapsing regions. We resort to a 3+1 splitting and obtain gauge invariant conditions relating intrinsic spacetimes quantities to properties of the matter source. We find that the dividing shell is defined by a generalization of the Tolman-Oppenheimer-Volkoff equilibrium condition. The latter establishes a balance between the pressure gradients, both isotropic and anisotropic, and the strength of the fields induced by the Misner-Sharp mass inside the separating shell and by the pressure fluxes. This defines a local equilibrium condition, but conveys also a non-local character given the definition of the Misner-Sharp mass. By the same token, it is also a generalized thermodynamical equation of state as usually interpreted for the perfect fluid case, which now has the novel feature of involving both the isotropic and the anisotropic stress. We have cast the governing equations in terms of local, gauge invariant quantities which are revealing of the role played by the anisotropic pressures and inhomogeneous electric part of the Weyl tensor. We analyse a particular solution with dust and radiation that provides an illustration of our conditions. In addition, our gauge invariant formalism not only encompasses the cracking process from Herrera and coworkers but also reveals transparently the interplay and importance of the shear and of the anisotropic stresses.

We use the spherical collapse model of structure formation to investigate the separation in the collapse of uncoupled matter (essentially baryons) and coupled dark matter in an interacting dark energy scenario. Following the usual assumption of a single radius of collapse for all species, we show that we only need to evolve the uncoupled matter sector to obtain the evolution for all matter components. This gives us more information on the collapse with a simplified set of evolution equations compared with the usual approaches. We then apply these results to four quintessence potentials and show how we can discriminate between different quintessence models.

Interaction between dark energy and dark matter is probed through deviation from the virial equilibrium for two relaxed clusters: A586 and A1689. The evaluation of the virial equilibrium is performed using realistic density profiles. The virial ratios found for the more realistic density profiles are consistent with the absence of interaction.

We analyse the dynamics of trapped matter shells in spherically symmetric inhomogeneous \Lambda-CDM models. The investigation uses a Generalised Lema\^itre-Tolman-Bondi description with initial conditions subject to the constraints of having spatially asymptotic cosmological expansion, initial Hubble-type flow and a regular initial density distribution. We discuss the effects of shell crossing and use a qualitative description of the local trapped matter shells to explore global properties of the models. Once shell crossing occurs, we find a splitting of the global shells separating expansion from collapse into, at most, two global shells: an inner and an outer limit trapped matter shell. In the case of expanding models, the outer limit trapped matter shell necessarily exists. We also study the role of shear in this process, compare our analysis with the Newtonian framework and give concrete examples using density profile models of structure formation in cosmology.

In this paper, we study the distribution functions that arise naturally during self-similar radial infall of collisionless matter. Such matter may be thought of either as stars or as dark matter particles. If a rigorous steady state is assumed, then the system is infinite and is described by a universal distribution function given the self-similar index. The steady logarithmic potential case is exceptional and yields the familiar Gaussian for an infinite system with an inverse-square density profile. We show subsequently that for time-dependent radial self-similar infall, the logarithmic case is accurately described by the Fridmann and Polyachenko distribution function. The system in this case is finite but growing. We are able to embed a central mass in the universal steady distribution only by iteration, except in the case of massless particles. The iteration yields logarithmic corrections to the massless particle case and requires a `renormalization' of the central mass. A central spherical mass may be accurately embedded in the Fridmann and Polyachenko growing distribution however. Some speculation is given concerning the importance of radial collisionless infall in actual galaxy formation.

Aims. In this paper we continue our study of density cusps that may contain central black holes. Methods. We recall our attempts to use distribution functions with a memory of self-similar relaxation, but mostly they apply only in restricted regions of the global system. We are forced to consider related distribution functions that are steady but not self-similar. Results. One remarkably simple distribution function that has a filled loss cone describes a bulge that transits from a near black hole domain to an outer 'zero flux' regime where$\rho\propto r^{-7/4}$. The transition passes from an initial inverse square profile through a region having a 1/r density profile. The structure is likely to be developed at an early stage in the growth of a galaxy. A central black hole is shown to grow exponentially in this background with an e-folding time of a few million years. Conclusions. We derive our results from first principles, using only the angular momentum integral in spherical symmetry. The initial relaxation probably requires bar instabilities and clump-clump interactions.

Aims: In this paper we study density cusps that may contain central black holes. The actual co-eval self-similar growth would not distinguish between the central object and the surroundings. Methods: To study the environment of a growing black hole we seek descriptions of steady `cusps' that may contain a black hole and that retain at least a memory of self-similarity. We refer to the environment in brief as the `bulge' and on smaller scales, the `halo'. Results: We find simple descriptions of the simulations of collisionless matter by comparing predicted densities, velocity dispersions and distribution functions with the simulations. In some cases central point masses may be included by iteration. We emphasize that the co-eval self-similar growth allows an explanation of the black hole bulge mass correlation between approximately similar collisionless systems. Conclusions: We have derived our results from first principles assuming adiabatic self-similarity and either self-similar virialisation or normal steady virialisation. We conclude that distribution functions that retain a memory of self-similar evolution provide an understanding of collisionless systems. The implied energy relaxation of the collisionless matter is due to the time dependence. Phase mixing relaxation may be enhanced by clump-clump interactions.

We investigate spherically symmetric perfect-fluid spacetimes and discuss the existence and stability of a dividing shell separating expanding and collapsing regions. We perform a 3+1 splitting and obtain gauge invariant conditions relating the intrinsic spatial curvature of the shells to the Misner-Sharp mass and to a function of the pressure that we introduce and that generalizes the Tolman-Oppenheimer-Volkoff equilibrium condition. We find that surfaces fulfilling those two conditions fit, locally, the requirements of a dividing shell and we argue that cosmological initial conditions should allow its global validity. We analyze the particular cases of the Lema\^itre-Tolman-Bondi dust models with a cosmological constant as an example of a cold dark matter model with a cosmological constant (\Lambda-CDM) and its generalization to contain a central perfect-fluid core. These models provide simple, but physically interesting illustrations of our results.

We discuss the existence of a dividing shell separating expanding and collapsing regions in spherically symmetric solutions with pressure. We obtain gauge invariant conditions relating not only the intrinsic spatial curvature of the shells to the ADM mass, but also a function of the pressure which we introduce that generalises the Tolman-Oppenheimer-Volkoff equilibrium condition, in the framework of a 3+1 spacetime splitting. We consider the particular case of a Lema\^itre-Tolman-Bondi dust models with a cosmological constant (a $\Lambda$-CDM model) as an example of our results.

We discuss the current bounds on the Equivalence Principle, in particular from structure formation and, reexamine in this context, the recent claim on the evidence of the interaction between dark matter and dark energy in the Abell Cluster A586 and the ensued violation of the Equivalence Principle.

We investigate spherically symmetric solutions with pressure and discuss the existence of a dividing shell separating expanding and collapsing regions. We perform a 3+1 splitting and obtain gauge invariant conditions relating not only the intrinsic spatial curvature of the shells to the ADM mass, but also a function of the pressure which we introduce that generalises the Tolman-Oppenheimer-Volkoff equilibrium condition. We consider the particular case of a Lema\^itre-Tolman dust models with a cosmological constant (a $\Lambda$-CDM model) as an example of our results.

The Self-Similar Secondary Infall Model (SSIM) is modified to simulate a merger event. The model encompass spherical versions of tidal stripping and dynamical friction that agrees with the Syer & White merger paradigm's behaviour. The SSIM shows robustness in absorbing even comparable mass perturbations and returning to its original state. It suggests the approach to be invertible and allows to consider accretion as smooth mass inflow merging and mergers as intermittent mass inflow accretion.

We find that deviation from the virial equilibrium of the Abell Cluster A586 yields evidence of the interaction between dark matter and dark energy. We argue that this interaction might imply a violation of the Equivalence Principle. Our analysis show that evidence is found in the context of two different models of dark energy-dark matter interaction.

We find that the Abell Cluster A586 exhibits evidence of the interaction between dark matter and dark energy and argue that this interaction suggests evidence of violation of the Equivalence Principle. This violation is found in the context of two different models of dark energy-dark matter interaction.

We show that the Abell Cluster A586 exhibits evidence of the interaction between dark matter and dark energy and argue that this interaction implies a violation of the Equivalence Principle. This violation is found in the context of two different models of dark energy-dark matter interaction. We also argue, based on the spherical symmetry of the Abell Cluster A586 that skewness is not the most general quantity to test the Equivalence Principle.

Measurement of accelerated expansion in the Universe led to propose a new cosmic fluid as its cause: dark energy. Its various incarnations offer a wealth of models whose relevance it is important to discriminate via contacts with observations. I will present my investigations on the influence of dynamical dark energy models on the formation of non-linear dark matter structures. In particular, I will focus on structures traced by the mass function of dark matter haloes.

We present an investigations on the influence of dynamical quintessence field models on the formation of non-linear structures. In particular, we focus on the structure traced by the mass function. Our contribution builds on previous studies by considering two potentials not treated previousy and investigating the impact of the free parameters of two other models. Our approach emphasises the physical insight into the key role of the evolution of the equation of state. We use the variations in the spherical non-linear collapse caused by quintessence as a function of the model's potential in a Press Schechter scheme to obtain the differences in halo mass functions. The comparison is also done with the more usual scenario $\Lambda$CDM. We conclude the method should be promissing for discrimination using cluster mass determination and that the key role lies within the evolution of the equation of state.

We present detailed predictions for the properties of Ly-alpha-emitting galaxies in the framework of the Lambda-CDM cosmology, calculated using the semi-analytical galaxy formation model GALFORM. We explore a model which assumes a top-heavy IMF in starbursts, and which has previously been shown to explain the sub-mm number counts and the luminosity function of Lyman-break galaxies at high redshift. We show that this model, with the simple assumption that a fixed fraction of Ly-alpha photons escape from each galaxy, is remarkably successful at explaining the observed luminosity function of Ly-alpha emitters over the redshift range 3<z<6.6. We also examine the distribution of Ly-alpha equivalent widths and the broad-band continuum magnitudes of emitters, which are in good agreement with the available observations. We look more deeply into the nature of Ly-alpha emitters, presenting predictions for fundamental properties such as the stellar mass and radius of the emitting galaxy and the mass of the host dark matter halo. The model predicts that the clustering of Ly-alpha emitters at high redshifts should be strongly biased relative to the dark matter, in agreement with observational estimates. We also present predictions for the luminosity function of Ly-alpha emitters at z>7, a redshift range which is starting to be be probed by near-IR surveys and using new instruments such as DAzLE.

We present predictions for the abundance of Ly-$\alpha$ emitters in hierarchical structure formation models. We use the {\tt GALFORM} semi-analytical model to explore the impact on the predicted counts of varying assumptions about the escape fraction of Ly-$\alpha$ photons, the redshift at which the universe reionised and the cosmological density parameter. A model with a fixed escape fraction gives a remarkably good match to the observed counts over a wide redshift interval. We present predictions for the expected counts in a typical observation with the Multi Unit Spectroscopic Explorer instrument proposed for the Very Large Telescope.

The self-similar infall model (SSIM) is normally discussed in the context of radial orbits in spherical symmetry. However it is possible to retain the spherical symmetry while permitting the particles to move in Keplerian ellipses, each having the squared angular momentum peculiar to their 'shell'. The spherical 'shell', defined for example by the particles turning at a given radius, then moves according to the radial equation of motion of a 'shell' particle. The 'shell' itself has no physical existence except as an ensemble of particles, but it is convenient to sometimes refer to the shells since it is they that are followed by a shell code. In this note we find the distribution of squared angular momentum as a function of radius that yields the NFW density profile for the final dark matter halo. It transpires that this distribution is amply motivated dimensionally. An effective 'lambda' spin parameter is roughly constant over the shells. We also study the effects of angular momentum on the relaxation of a dark matter system using a three dimensional representation of the relaxed phase space.

The density profiles established in the self-similar dynamical phase of dark matter haloes recollapse are all close to isothermal. This is steeper than the predictions of some n-body simulations for the central regions of the halo, which are in turn steeper than the profiles observed in the centre of some galaxies; particularly dwarfs and low surface brightness galaxies. The outer regions of galaxies both as observed and as simulated have density profiles steeper than the self-similar profile. Nevertheless there appears to be an intermediate region in most galaxies in which the inverse square behaviour is a good description. The outer deviations can be explained plausibly in terms of the transition from a self-gravitating extended halo to a Keplerian flow onto a dominant central mass (the isothermal distribution can not be complete), but the inner deviations are more problematic. Rather than attack this question directly, we use in this paper a novel coarse-graining technique combined with a shell code to establish both the distribution function associated with the self-similar density profile and the nature of the possible deviations in the central regions. In spherical symmetry we find that both in the case of purely radial orbits and in the case of orbits with non-zero angular momentum the self-similar density profile should flatten progressively near the centre of the system. The NFW limit of -1 seems possible. In a section aimed at demonstrating our technique for a spherically symmetric steady state, we argue that a Gaussian distribution function is the best approximation near the centre of the system.