- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
Advances in Astronomy
Volume 2011 (2011), Article ID 903429, 8 pages
Dark Matter Halos: The Dynamical Basis of Effective Empirical Models
1Department of Physics, University of Rome ‘Tor Vergata’, Via della Ricerca Scientifica 1, 00133 Rome, Italy
2Astrophysics Sector, SISSA/ISAS, Via Bonomea 265, 34136 Trieste, Italy
Received 27 May 2010; Accepted 31 August 2010
Academic Editor: Alberto J. Castro-Tirado
Copyright © 2011 Andrea Lapi and Alfonso Cavaliere. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
We investigate the dynamical basis of the classic empirical models (specifically, Sérsic-Einasto and generalized NFW) that are widely used to describe the distributions of collisionless matter in galaxies. We submit that such a basis is provided by our -profiles, shown to constitute solutions of the Jeans dynamical equilibrium with physical boundary conditions. We show how to set the parameters of the empirical in terms of the dynamical models; we find the empirical models, and specifically Sérsic-Einasto, to constitute a simple and close approximation to the dynamical models. Finally, we discuss how these provide a useful baseline for assessing the impact of the small-scale dynamics that may modulate the density slope in the central galaxy regions.
The classic Sérsic  models met a wide and lasting success as empirical representations of the projected (-dimensional) light distributions in spheroidal galaxies (for a review, see Kormendy et al. ). Einasto  developed and used a similar shape to describe in simple terms -dimensional stellar mass profiles.
On the other hand, recent extensive -body simulations (e.g., [4–8]) indicate that the Sérsic and Einasto functional forms also provide good patterns to represent the spherically-averaged mass distributions in dark matter (DM) halos ranging from galaxies to galaxy clusters. These hold at levels comparable to, or even better than the popular NFW formula .
Still, no agreed understanding is available to explain the value in both the real and the virtual worlds of the Sérsic and Einasto representations (see [2, 10]). Can we identify the underlying astrophysical basis?
2. Empirical Models
2.1. Density Runs
The density runs of the SE family may be represented in the form Here, quantities are normalized to their value at , the reference radius where the logarithmic slope takes on the value ; typically, in nearby elliptical galaxies, corresponds to sizes of order kpc, a few times the half-light radius .
The parameters and describe the inner slope and the middle curvature of the density run, respectively. The original Einasto profile belongs to this family and is obtained when . Note, however, that by deprojecting from the plane of the sky, a Sérsic -dimensional run with index (suited for normal ellipticals, see ) produces a cuspy inner run as in (1) with significantly different from and less than , as shown by Prugniel and Simien .
On the DM side, recent simulations (see [4, 6, 8]) only provide an upper bound for the inner slope. When the original Einasto profile (with ) is adopted, the best-fit to simulated DM halos obtains for ; we will come back to this value later on.
In turn, the density runs of the gNFW family may be written in the following form(In the literature these runs are sometimes referred to as -models, and equivalently defined via the parameters , , and .): the parameters , , and describe the central slope, the middle curvature, and the outer decline of the density run, respectively. Note that familiar empirical profiles are recovered for specific values of the triple (); for example, Plummer's  corresponds to , Jaffe's  to , Herquist's  to , and NFW to .
2.2. Toward a Single Family
The main apparent difference between SE and gNFW is constituted by the former's exponential decline versus the latter's powerlaw falloff for large .
On the other hand, (2) is to be considered for large values of anyway, since a steep density run in the halo outskirts is indicated by observations of light distribution in spheroidal galaxies (other than cDs, see ) and of DM distributions from weak lensing in galaxies and galaxy clusters (e.g., [18–20]).
The circumstance is easily translated into the formal statement that the gNFW family converges to the SE for large . This is seen on recasting from (2) in exponential form, to read for approximating the middle and last terms, we have used the circumstance that implies and so applies. Thus the two families in (1) and (2) actually become one in this limit.
Thus in the following we focus mainly on the SE family, and proceed to discuss its dynamical basis in terms of the Jeans equation.
3. The Dynamical Model
The dynamical model of DM halos hinges upon the radial Jeans equation that expresses the self-gravitating, equilibrium of collisionless matter (see ). The Jeans equation reads in terms of the density , the related cumulative mass , and the radial velocity dispersion . The last term on the r.h.s. describes the effects of anisotropic random velocities via the standard Binney  parameter .
Note that the Jeans equation is designed to describe a (quasi-)static equilibrium, away from extreme major merger events like in the case of the Bullet Cluster (see ). But even in relaxed conditions, solving Jeans requires an “equation of state,” that is, a functional relation expressing the DM pressure in terms of density (and possibly radius) only.
3.1. Equation of State
In seeking for such a relation, one can make contact with the classic theory of the nonlinear collapse for DM halos in an expanding Universe; here self-similar arguments play the role of a pivotal pattern (see [24–26]). This still applies to modern views of the halo development (e.g., [27–31]), that comprise two stages: an early collapse that builds up the halo main body via a few major merger events and sets its phase-space structure by dynamical relaxation of DM particle orbits; this tails off into a secular development of the outskirts by smooth accretion and minor mergers.
The essence of the macroscopic equilibrium is conveyed by the self-similar scaling adding to the geometric relation . The macroscopic import of the halo phase-space structure is conveyed by combining these two quantities into the “phase-space density” or equivalently into the functional often referred to as DM “entropy” (see [24, 26]). For the latter quantity, one easily derives the scaling implying whence one expects a slope slightly exceeding .
To focus the values of , in Lapi and Cavaliere [27, 28] we have developed a full semianalytic treatment of the halo growth in the standard accelerating Universe (see ). We found constant values of , that fall within the narrow range 1.25–1.3; on average, such values grow weakly with the mass of the halo body, from galaxies to rich clusters.
The halo development process has been probed, and the two-stage view confirmed by intensive -body simulations (e.g., [33–39]). These also confirm that (i) a (quasi-)static macroscopic equilibrium is attained at the end of the fast collapse and is retained during the subsequent stage of secular, smooth mass addition; (ii) a persistent feature of such an equilibrium is constituted by powerlaw correlations holding in the form , although it is still widely debated whether the radial or the total velocity dispersion best applies (see also the discussions by [6, 36]).
In building up our dynamical models, we focus on the quantity that involves the radial dispersion (see also ). Operationally, this provides a direct expression for the radial pressure term in the Jeans (4); anisotropies are accounted for by the last term on the r.h.s., as discussed in Section 3.3 below.
3.2. The DM -Profiles
In terms of , the Jeans equation may be recast into the compact form with representing the logarithmic density slope and the circular velocity. Remarkably, by double differentiation, this integrodifferential equation for reduces to a handy second order differential equation for (see [40, 41]).
Tackling first the isotropic case , we recall that the solution space of (6) spans the range : the specific solution for the upper bound, and the behaviors of the others have been analytically investigated by Taylor and Navarro , Austin et al. , and Dehnen and McLaughlin . In Lapi and Cavaliere , we explicitly derived the solutions in the full range , that are marked by a monotonically decreasing run and satisfy physical boundary conditions: a finite central pressure or energy density (equivalent to a round minimum of the gravitational potential) and a steep outer run implying a finite and rapidly converging (hence a definite) overall mass. We dubbed -profiles these physical solutions.
We shall use the following basic features of the latter. In the halo body at the point the -profile is tangent to the pure powerlaw solution of the Jeans equation; there applies, to imply from (6) This is consistent for with the self-similar slope. Note that the point lies in the neighborhood of the radius (see Section 2.1); specifically holds for .
On the other hand, a monotonic density run implies the term to vanish at the center; this results in an inner powerlaw with the slope This differs from zero as long as the entropy run grows from the center following with .
Finally, a finite mass implies to hold in the outskirts, so as to yield a typical outer decline with slope This exceeds the value and so constitutes the hallmark of a rapidly saturating mass; the circumstance makes less compelling here the role of a virial boundary.
Thus, compared to NFW the inner slope of the dynamical model is considerably flatter and the outer slope steeper; compared to the original Einasto profile, the main difference occurs in the inner regions where the dynamical model is (moderately) steeper.
It is clear from (6) that anisotropies will steepen the density run for positive , and flatten it for negative . The latter condition is expected to prevail in the inner region, where tangential components develop from the angular momentum barrier [30, 42]. Moving outwards, radial motions are expected to prevail, so raising up to values around at ; outwards of this, is expected to saturate or even decrease, as one enters a region increasingly populated by DM particles on eccentric orbits with vanishing radial dispersions at their apocenters (see ).
This view is supported by numerical simulations (see [6, 40, 41, 43]), which in detail suggest the average anisotropy-density slope relation to hold with parameters , , and the constraint . Note that at all radii the inequality is satisfied; this has been conjectured to constitute a necessary condition for a self-consistent spherical model with positive distribution function (see ).
In Lapi and Cavaliere , we extended the dynamical model to such anisotropic conditions in the full range , inspired by the analysis of Dehnen and McLaughlin  for the upper bound of . We note that the latter is now slightly modified to ; likewise, the point where applies moves slightly inwards, so that now holds.
The main outcome, however, is that the density profile is somewhat flattened at the center relative to the isotropic case; the inner slope now reads In particular, even a limited central anisotropy (corresponding to values ) causes an appreciable flattening down to for .
On the other hand, we stress that such small phenomenological anisotropies near the center imply the radial and the total dispersions to be very close.
4. From Dynamical to Empirical Models
Here we discuss how the parameters of the empirical profiles (see Section 2) can be set based on our dynamical model (see Section 3); in such conditions, it will turn out that such profiles constitute close approximations to the model over a wide radial range.
4.1. Parameters from Dynamics
First we consider the original Einasto profile ( in (1)), since this has been widely used in the context of DM halo simulations. Here, is fixed to , and the only free parameter is the curvature . This we set by requiring the logarithmic density slope to equal at the point . So we find the expression that takes on values ; see Tables 1 and 2; remarkably, these turn out to agree with those derived from fitting the outcomes of state-of-the-art -body simulations in terms of the same Einasto density run, as performed by Navarro et al. .
On the other hand, the flat central slope of the Einasto profile is at variance with the value given by our dynamical models based on Jeans; to wit, consistency between pure Einasto and Jeans would require at the center a flat entropy distribution and a vanishing pressure.
The latter features two parameters, the inner slope and the middle curvature . These we set by requiring the logarithmic density slope to equal for , and at ; so we find Thus we predict the central slope to take on values and the corresponding curvature parameter to take on values ; see Tables 1 and 2. It will be worth fitting the outcomes of -body simulations based on these extended SE profiles with .
Finally, we report the corresponding results for the empirical gNFW family. This features three parameters: inner slope , middle curvature , and strength of the outer decline ; these we set by requiring the logarithmic density slope to equal for , at , and for . So we find
4.2. Results and Comments
With the parameters focused as discussed in the previous subsection, Figure 1 illustrates how the empirical compare with our dynamical models. We plot the density run of the latter (specifically for the -profile with suitable to galactic halos), compared to those of the Einasto, SE, and gNFW models. The left and right panels refer to isotropic and anisotropic conditions, respectively; the popular NFW profile is also shown for reference. To make comparisons easier, we plot in the lower panels the corresponding logarithmic density slopes.
It turns out that the closest approximation to the dynamical model is provided by SE, which shares with it not only the central slope by construction, but also the body and the outer behaviors. The original Einasto profile provides an acceptable approximation in the middle and outer ranges, but not at the center, because of its flatness. On the other hand, the gNFW family provides an acceptable approximation in the inner and middle ranges, but not in the outskirts where its slope is too flat. Finally, the NFW profile provides an acceptable approximation in the middle range.
Similar conclusions concern the profiles of circular velocities , that are analytically dealt with in the Appendix and illustrated in Figure 2.
We stress that the handy SE representation is convenient in analyzing data in several contexts, including the DM particle annihilation signal expected from the Galactic Center (see ), rotation curves of dwarf and normal spiral galaxies (see ); individual and statistical properties of elliptical and spiral galaxies (see ), strong and weak gravitational lensing (see )—currently observed in clusters (e.g., ) and soon in massive elliptical galaxies (see discussion by ) and X-ray emission from the intracluster plasma (see [50–52]).
We first stress that the dynamical model (as well as its approximations in terms of empirical models) is in keeping with the basic features of standard DM, that is, its cold and collisionless nature. In fact, it implies for large , a behavior expected in the outskirts for cold matter dominating the potential well.
At the inner end, with decreasing we expect to increase toward a maximum, corresponding to effective conversion of inflow kinetic into random energy. In fact, toward the center Jeans requires to hold as the gravitational force vanishes there, to the effect that .
Concerning the collisionless nature of the DM, the boundary conditions at the center imply a finite, non-zero pressure (and energy density), while a long collisional mean free path allows the pressure gradient to diverge. Conversely, with a short mean free path the pressure gradient cannot diverge on scales , where a finite and a flatter apply. In fact, weakly collisional conditions have been proposed to explain the cored light profiles observed in many spheroidal galaxies (see ).
On approaching the center of a galactic halo, one expects the basic dynamical model from large-scale Jeans equilibrium to be altered to an increasing degree by small-scale dynamics and/or energetics related to baryons. These processes are specifically related to following issues: transfer of energy/angular momentum from baryons to DM during galaxy formation, scouring baryons by the energy feedback from central active galactic nuclei and any “adiabatic” contraction of the baryons. Such issues will be briefly discussed in turn, with a warning that they enter increasingly debated grounds.
5.1. Energy/Angular Momentum Transfers
In detail, upon transfer of tangential random motions from the baryons to an initially isotropic DM structure, the density in the inner region is expected to behave as  Thus for the profile is flattened relative to the original , down to the point of developing a core for .
However, a reliable assessment of the amount of angular momentum transferred from the baryons to the DM is still wanting, and would require aimed numerical simulations of better resolution than presently achieved.
5.2. Other Processes on Inner Scales
Less agreed processes may affect galactic scales pc, for example, at the formation of a spheroid, central starbursts and a supermassive black hole may easily discharge enough energy ( erg for a black hole mass ) with sufficient coupling () to blow most of the gaseous baryonic mass out of the gravitational potential well. This will cause an expansion of the DM and of the stellar distributions (see ), that flattens the central slope.
In addition, binary black hole dynamics following a substantial merger may eject on longer timescales formed stars from radii pc containing an overall mass of a few times the black hole's, and so may cause a light deficit in some galaxy cores (see [2, 58, 59]). A full discussion of the issue concerning cored versus cusped ellipticals is beyond the scope of the present paper.
5.3. Adiabatic Contraction?
On the other hand, some steepening of the inner density profile may be induced by any “adiabatic” contraction of the diffuse star-forming baryons into a disc-like structure, as proposed by Blumenthal et al.  and Mo et al.  but currently under scrutiny, see Abadi et al. .
On the basis of the standard treatments, it is easily shown that in the inner region an initial powerlaw is modified into this yields typical slopes around , steeper than the original but still significantly flatter than .
However, recent numerical simulations (see discussion by ) suggest that the treatment of adiabatic contraction leading to (19) is likely to be extreme; actually, in the inner region the contraction is ineffective and the density slope hardly modified. Again, highly resolved -body experiments are needed to clarify the issue.
We have discussed the dynamical basis of the Sérsic-Einasto empirical models, in terms of well-behaved solutions of the Jeans equation with physical boundary conditions comprising: a finite central energy density, a closely self-similar body, and a finite (definite) overall mass.
We find the SE profile to be particularly suitable to accommodate the general run of the dynamical solution. Specifically, we have discussed how to tune the parameters of SE in terms of the dynamical model; in such conditions, we find the former to constitute a simple and close approximation to the latter.
The resulting SE profile shares with the dynamical model the following features: an outer steep decline, hence a definite overall mass, a closely self-similar body with slope , and an inner slope around , hence flatter than . The latter slope provides an useful baseline for discussing alterations of the inner behavior caused by additional baryonic processes.
In conclusion, we submit that the dynamical models discussed here, namely, the -profiles, provide the astrophysical basis for understanding the empirical success of the SE profiles in fitting the real and the virtual observables, from galaxies to galaxy clusters.
The paper is supported by ASI and INAF. It is a pleasure to acknowledge inspiring correspondence with A. Graham, and stimulating discussions with L. Danese, P. Salucci, and Y. Rephaeli. A. Lapi thanks SISSA and INAF-OATS for warm hospitality.
- J.-L. Sérsic, “Influence of the atmospheric and instrumental dispersion on the brightness distribution in a galaxy,” Boletin de la Asociacion Argentina de Astronomia, vol. 6, p. 41, 1963.
- J. Kormendy, D. B. Fisher, M. E. Cornell, and R. Bender, “Structure and formation of elliptical and spheroidal galaxies,” Astrophysical Journal, Supplement Series, vol. 182, no. 1, pp. 216–309, 2009.
- J. Einasto, “Influence of the atmospheric and instrumental dispersion on the brightness distribution in a galaxy,” Trudy Inst. Astroz. Alma-Ata, vol. 5, p. 87, 1965.
- L. Gao, J. F. Navarro, S. Cole et al., “The redshift dependence of the structure of massive Λ cold dark matter haloes,” Monthly Notices of the Royal Astronomical Society, vol. 387, no. 2, pp. 536–544, 2008.
- D. Merritt, J. F. Navarro, A. Ludlow, and A. Jenkins, “A universal density profile for dark and luminous matter?” Astrophysical Journal, vol. 624, no. 2, pp. L85–L88, 2005.
- J. F. Navarro, A. Ludlow, V. Springel et al., “The diversity and similarity of simulated cold dark matter haloes,” Monthly Notices of the Royal Astronomical Society, vol. 402, no. 1, pp. 21–34, 2010.
- J. F. Navarro, E. Hayashi, C. Power et al., “The inner structure of ΛCDM haloes—III: universality and asymptotic slopes,” Monthly Notices of the Royal Astronomical Society, vol. 349, no. 3, pp. 1039–1051, 2004.
- J. Stadel, et al., “Quantifying the heart of darkness with GHALO—a multibillion particle simulation of a galactic halo,” Monthly Notices of the Royal Astronomical Society, vol. 398, pp. L21–L25, 2009.
- J. F. Navarro, C. S. Frenk, and S. D. M. White, “A universal density profile from hierarchical clustering,” Astrophysical Journal, vol. 490, no. 2, pp. 493–508, 1997.
- A. W. Graham, D. Merritt, B. Moore, J. Diemand, and B. Terzić, “Empirical models for dark matter halos. II. Inner profile slopes, dynamical profiles, and ρ/σ3,” Astronomical Journal, vol. 132, no. 6, pp. 2701–2710, 2006.
- L. Hernquist, “An analytical model for spherical galaxies and bulges,” Astrophysical Journal, vol. 356, no. 2, pp. 359–364, 1990.
- L. M. Widrow, “Distribution functions for cuspy dark matter density profiles,” Astrophysical Journal, Supplement Series, vol. 131, no. 1, pp. 39–46, 2000.
- H. Zhao, “Analytical models for galactic nuclei,” Monthly Notices of the Royal Astronomical Society, vol. 278, no. 2, pp. 488–496, 1996.
- D. Merritt, A. W. Graham, B. Moore, J. Diemand, and B. Terzić, “Empirical models for dark matter halos. I. Nonparametric construction of density profiles and comparison with parametric models,” Astronomical Journal, vol. 132, no. 6, pp. 2685–2700, 2006.
- PH. Prugniel and F. Simien, “The fundamental plane of early-type galaxies: non-homology of the spatial structure,” Astronomy and Astrophysics, vol. 321, no. 1, pp. 111–122, 1997.
- H. C. Plummer, “On the problem of distribution in globular star clusters,” Monthly Notices of the Royal Astronomical Society, vol. 71, pp. 460–470, 1911.
- W. Jaffe, “A simple model for the distribution of light in spherical galaxies,” Monthly Notices of the Royal Astronomical Society, vol. 202, pp. 995–999, 1983.
- T. Broadhurst, K. Umetsu, E. Medezinski, M. Oguri, and Y. Rephaeli, “Comparison of cluster lensing profiles with ACDM predictions,” Astrophysical Journal, vol. 685, no. 1, pp. L9–L12, 2008.
- A. B. Newman, T. Treu, R. S. Ellis et al., “The distribution of dark matter over three decades in radius in the lensing cluster abell 611,” Astrophysical Journal, vol. 706, no. 2, pp. 1078–1094, 2009.
- M. Oguri, J. F. Hennawi, M. D. Gladders et al., “Subaru weak lensing measurements of four strong lensing clusters: are lensing clusters overconcentrated?” Astrophysical Journal, vol. 699, no. 2, pp. 1038–1052, 2009.
- J. Binney and S. Tremaine, Galactic Dynamics, Princeton University Press, Princeton, NJ, USA, 2nd edition, 2008.
- J. Binney, “Twisted and warped disks as consequences of heavy halos,” Monthly Notices of the Royal Astronomical Society, vol. 183, p. 779, 1978.
- D. Clowe, M. Bradač, A. H. Gonzalez et al., “A direct empirical proof of the existence of dark matter,” Astrophysical Journal, vol. 648, no. 2, pp. L109–L113, 2006.
- E. Bertschinger, “Self-similar secondary infall and accretion in an Einstein-de Sitter universe,” The Astrophysical Journal Supplement Series, vol. 58, p. 39, 1985.
- J. A. Fillmore and P. Goldreich, “Self-similar gravitational collapse in an expanding universe,” Astrophysical Journal, vol. 281, pp. 1–8, 1984.
- J. E. Taylor and J. F. Navarro, “The phase-space density profiles of cold dark matter halos,” Astrophysical Journal, vol. 563, no. 2, pp. 483–488, 2001.
- A. Lapi and A. Cavaliere, “Structure and history of dark matter halos probed with gravitational lensing,” Astrophysical Journal, vol. 695, no. 2, pp. L125–L129, 2009.
- A. Lapi and A. Cavaliere, “Dark matter equilibria in galaxies and galaxy systems,” Astrophysical Journal, vol. 692, pp. 174–186, 2009.
- Y. Li, H. J. Mo, F. C. van den Bosch, and W. P. Lin, “On the assembly history of dark matter haloes,” Monthly Notices of the Royal Astronomical Society, vol. 379, no. 2, pp. 689–701, 2007.
- Y. Lu, H. J. Mo, N. Katz, and M. D. Weinberg, “On the origin of cold dark matter halo density profiles,” Monthly Notices of the Royal Astronomical Society, vol. 368, no. 4, pp. 1931–1940, 2006.
- H. J. Mo and S. Mao, “Galaxy formation in pre-processed dark haloes,” Monthly Notices of the Royal Astronomical Society, vol. 353, no. 3, pp. 829–840, 2004.
- E. Komatsu, et al., “Seven-year wilkinson microwave anisotropy probe (WMAP) observations: cosmological interpretation,” submitted to The Astrophysical Journal Supplement Series, http://arxiv.org/abs/1001.4538.
- Y. Ascasibar and S. Gottlöber, “The dynamical structure of dark matter haloes,” Monthly Notices of the Royal Astronomical Society, vol. 386, no. 4, pp. 2022–2030, 2008.
- J. Diemand, M. Kuhlen, and P. Madau, “Formation and evolution of galaxy dark matter halos and their substructure,” Astrophysical Journal, vol. 667, no. 2, pp. 859–877, 2007.
- Y. Hoffman, E. Romano-Díaz, I. Shlosman, and C. Heller, “Evolution of the phase-space density in dark matter halos,” Astrophysical Journal, vol. 671, no. 2, pp. 1108–1114, 2007.
- K. B. Schmidt, S. H. Hansen, and A. V. Macciò, “Alas, the dark matter structures were not that trivial,” Astrophysical Journal, vol. 689, no. 1, pp. L33–L36, 2008.
- I. M. Vass, M. Valluri, A. V. Kravtsov, and S. Kazantzidis, “Evolution of the dark matter phase-space density distributions of ΛcDM haloes,” Monthly Notices of the Royal Astronomical Society, vol. 395, no. 3, pp. 1225–1236, 2009.
- R. H. Wechsler, A. R. Zentner, J. S. Bullock, A. V. Kravtsov, and B. Allgood, “The dependence of halo clustering on halo formation history, concentration, and occupation,” Astrophysical Journal, vol. 652, no. 1 I, pp. 71–84, 2006.
- D. H. Zhao, H. J. Mo, Y. P. Jing, and G. Börner, “The growth and structure of dark matter haloes,” Monthly Notices of the Royal Astronomical Society, vol. 339, no. 1, pp. 12–24, 2003.
- W. Dehnen and D. E. McLaughlin, “Dynamical insight into dark matter haloes,” Monthly Notices of the Royal Astronomical Society, vol. 363, no. 4, pp. 1057–1068, 2005.
- C. G. Austin, L. L. R. Williams, E. I. Barnes, A. Babul, and J. J. Dalcanton, “Semianalytical dark matter halos and the Jeans equation,” Astrophysical Journal, vol. 634, no. 2, pp. 756–774, 2005.
- A. Nusser, “Self-similar spherical collapse with non-radial motions,” Monthly Notices of the Royal Astronomical Society, vol. 325, no. 4, pp. 1397–1401, 2001.
- S. H. Hansen and B. Moore, “A universal density slope—velocity anisotropy relation for relaxed structures,” New Astronomy, vol. 11, no. 5, pp. 333–338, 2006.
- L. Ciotti and L. Morganti, “Consistency criteria for generalized Cuddeford systems,” Monthly Notices of the Royal Astronomical Society, vol. 401, no. 2, pp. 1091–1098, 2010.
- A. Lapi, et al., “Gamma rays from annihilations at the galactic center in a physical dark matter distribution,” Astronomy & Astrophysics, vol. 510, p. 90, 2010.
- P. Salucci, A. Lapi, C. Tonini, G. Gentile, I. Yegorova, and U. Klein, “The universal rotation curve of spiral galaxies—II: the dark matter distribution out to the virial radius,” Monthly Notices of the Royal Astronomical Society, vol. 378, no. 1, pp. 41–47, 2007.
- M. Cook, A. Lapi, and G. L. Granato, “Two-phase galaxy formation,” Monthly Notices of the Royal Astronomical Society, vol. 397, no. 1, pp. 534–547, 2009.
- A. Zitrin, et al., “Strong-lensing analysis of a complete sample of 12 MACS clusters at : mass models and Einstein radii,” submitted to Monthly Notices of the Royal Astronomical Society, http://arxiv.org/abs/1002.0521.
- M. Bradač, T. Treu, D. Applegate et al., “Focusing cosmic telescopes: exploring redshift z 5-6 galaxies with the bullet cluster 1E0657 56,” Astrophysical Journal, vol. 706, no. 2, pp. 1201–1212, 2009.
- A. Cavaliere, A. Lapi, and R. Fusco-Femiano, “Galaxy clusters: a novel look at diffuse baryons withstanding dark matter gravity,” Astrophysical Journal, vol. 698, no. 1, pp. 580–593, 2009.
- R. Fusco-Femiano, A. Cavaliere, and A. Lapi, “Supermodel analysis of galaxy clusters,” Astrophysical Journal, vol. 705, no. 1, pp. 1019–1030, 2009.
- A. Lapi, et al., “Probing the astrophysics of cluster outskirts,” Astronomy & Astrophysics, vol. 516, p. 34, 2010.
- J. P. Ostriker, “Collisional dark matter and the origin of massive black holes,” Physical Review Letters, vol. 84, no. 23, pp. 5258–5260, 2000.
- A. El-Zant, I. Shlosman, and Y. Hoffman, “Dark halos: the flattening of the density cusp by dynamical friction,” Astrophysical Journal, vol. 560, no. 2, pp. 636–643, 2001.
- E. Romano-Díaz, I. Shlosman, Y. Hoffman, and C. Heller, “Erasing dark matter cusps in cosmological galactic halos with baryons,” Astrophysical Journal, vol. 685, no. 2, pp. L105–L108, 2008.
- C. Tonini, A. Lapi, and P. Salucci, “Angular momentum transfer in dark matter halos: erasing the CUSP,” Astrophysical Journal, vol. 649, no. 2, pp. 591–598, 2006.
- L. Fan, A. Lapi, G. De Zotti, and L. Danese, “The dramatic size evolution of elliptical galaxies and the quasar feedback,” Astrophysical Journal, vol. 689, no. 2, pp. L101–L104, 2008.
- T. R. Lauer, K. Gebhardt, S. M. Faber et al., “The centers of early-type galaxies with hubble space telescope. VI. Bimodal central surface brightness profiles,” Astrophysical Journal, vol. 664, no. 1, pp. 226–256, 2007.
- D. Merritt, “Single and binary black holes and their influence on nuclear structure,” in Carnegie Observatories Astrophysics Series Vol. 1: Coevolution of Black Holes and Galaxies, L. C. Ho, Ed., p. 263, Cambridge Univ. Press, Cambridge, UK, 2004.
- G. R. Blumenthal, S. M. Faber, R. Flores, and J. R. Primack, “Contraction of dark matter galactic halos due to baryonic infall,” The Astrophysical Journal, vol. 301, pp. 27–34, 1986.
- H. J. Mo, S. Mao, and S. D. M. White, “The formation of galactic discs,” Monthly Notices of the Royal Astronomical Society, vol. 295, no. 2, pp. 319–336, 1998.
- M. G. Abadi, et al., “Galaxy-induced transformation of dark matter halos,” Monthly Notices of the Royal Astronomical Society, vol. 407, no. 1, pp. 435–446, 2010.