Pith. Sign in

REVIEW

The Inner Structure of LambdaCDM Halos III: Universality and Asymptotic Slopes

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Julio F. Navarro (1), Eric Hayashi (1), Chris Power (2), Adrian Jenkins (2), Carlos S. Frenk (2), Simon D.M. White (3), Volker Springel (3), Joachim Stadel (4)
show 5 more authors
Thomas R. Quinn (5) ((1) University of Victoria, (2) University fo Durham, (3) Max-Planck Institute for Astrophysics, (4) University of Zurich, (5) University of Washington)

Astrophysics · Submitted 11 Nov 2003 · 37 references

Pith's verdict: The mass-budget bound at the convergence radius is the real contribution here — it kills the r^-1.5 cusp without leaning on any fitting formula, and that argument survives scrutiny even though the paper oversells it by a notch.

T0 · AI review

How this review was made: An AI referee (selfhosted/deepseek-v4.1-flash) read the paper on October 29, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 10 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv IDastro-ph/0311231 (version 2)
Pith IDS3MNWOPQ, Pith's permanent record of this paper and its review
Citation(2) University fo Durham, (3) Max-Planck Institute for Astrophysics, (4) University of Zurich, (5) University of Washington), Adrian Jenkins (2), Carlos S. Frenk (2), Chris Power (2), Eric Hayashi (1), Joachim Stadel (4), Julio F. Navarro (1), Simon D.M. White (3), Thomas R. Quinn (5) ((1) University of Victoria, Volker Springel (3). "The Inner Structure of LambdaCDM Halos III: Universality and Asymptotic Slopes." arXiv:astro-ph/0311231 (2003). https://pith.science/paper/S3MNWOPQ
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Nineteen simulated dark matter halos, from dwarf galaxies to galaxy clusters, show inner density profiles that keep flattening toward the center and never settle onto a fixed power law, which rules out the steepest proposed cusp and points,

The reading

The paper asks a narrow but consequential question: as you walk toward the center of a cold dark matter halo, does the density profile approach a definite power law $r^{-1}$, $r^{-1.5}$, or something else? Earlier simulations disagreed, and the disagreement mattered because a steeper cusp puts more dark matter where gamma-ray annihilation signals would be brightest and where rotation curves are hardest to reconcile with observed cores. Using nineteen resimulated halos spanning five decades in mass, each with roughly a million particles inside the virial radius and resolved down to about half a percent of it, the authors argue that the density simply keeps getting shallower inward, with no sign of settling. Their sharpest claim is a mass-budget argument: the local density and the enclosed mean density at the innermost resolved radius together cap how steep the central slope can possibly be, and that cap sits below $1.5$ for almost every halo in the sample. That eliminates the $r^{-1.5}$ profile as a description of these halos while leaving the original $r^{-1}$ Navarro-Frenk-White slope still allowed. The authors then offer a replacement fitting formula whose slope is a power of radius, which fits individual halos to better than ten percent with no leftover radial trend.

What carries the argument

Two devices carry the argument. The first is the upper bound on the asymptotic inner slope, written as $d$, evaluated at the innermost converged radius. It uses only two measured quantities, the local density and the enclosed mean density, and it turns the question of the central cusp into a statement about how much mass sits inside the smallest resolved sphere. The second is a replacement fitting formula in which the logarithmic slope is a power law in radius, $d$, which integrates to a density profile with a finite central value but a very slow approach to it. With a single free parameter $d$ in the range $0.1$ to $0.2$ (mean $0.172$, scatter $0.032$), this form reproduces individual halo profiles to better than about ten percent across the resolved range and leaves no systematic radial residual, unlike the Navarro-Frenk-White and Moore fits whose worst deviations sit at the innermost point. The same formula, at larger $d$, also reproduces the parabolic circular-velocity shapes used to describe substructure halos.

Extended reading notes

Core claim

The central claim is that the inner slope of the ΛCDM halo density profile does not converge to a well-defined asymptotic value anywhere the simulations resolve. The authors measure $d$, the logarithmic slope $d$, as a function of radius for all nineteen halos and find it decreasing monotonically inward, reaching only about $1.1$ for clusters, $1.2$ for galaxy-mass halos, and $1.35$ for dwarfs at the innermost converged radius. That is steeper than the Navarro-Frenk-White formula predicts at the same scaled radius but shallower than the $r^{-1.5}$ form proposed by Moore and collaborators. The load-bearing step is a bound derived from the density budget: at any radius, $d$, where the bound is the ratio of local density to mean enclosed density. Under the assumption that the slope varies monotonically with radius, that number is an upper limit on the central slope. Applied at the innermost resolved point, it excludes $r^{-1.5}$ in essentially all cases while still permitting the $r^{-1}$ Navarro-Frenk-White asymptote. The authors read the widely reported inner density excess relative to Navarro-Frenk-White fits not as evidence of a steeper cusp but as a symptom of a fitting formula whose slope changes too abruptly with radius.

Load-bearing premise

The exclusion of a steep central cusp rests entirely on the assumption that the density slope changes monotonically with radius; if the slope were to turn back up inside the smallest resolved radius, the mass measured there would say nothing about the true center.

Editorial extensions

If this is right

  • The $r^{-1.5}$ Moore profile should no longer be used as the default description of ΛCDM halo centers when estimating annihilation signals or comparing to rotation-curve data.
  • Extrapolating Navarro-Frenk-White fits inward of the resolved radius will systematically underestimate the dark matter density, and the error grows as one goes deeper.
  • The new slope-power-law formula gives a safer inward extrapolation, because its residuals show no radial trend over the resolved range.
  • Because the profiles of dwarfs and clusters coincide once scaled to the radius where the slope equals two, a single shape plus a mass-dependent concentration is enough to describe halos across five decades in mass.
  • The concentration-mass relation predicted by the Eke-Navarro-Steinmetz and Bullock et al. formalisms reproduces the measured characteristic radii and densities reasonably well, so either can supply the scaling parameters the new formula needs.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • If the slope really keeps flattening without limit, then any finite-resolution simulation will always report a slope that depends on its resolution, which would explain why published values drifted from $1$ to $1.5$ and back as codes improved; the disagreement would be a resolution artifact rather than a physical disagreement.
  • The bound on the central slope is a statement about enclosed mass, so it can be sharpened cheaply: doubling the resolved particle count inside the convergence radius tightens the cap directly, without needing a better fit to the profile shape.
  • A finite central density, as the new formula implies, would soften the strongest small-scale objections to cold dark matter from observed constant-density cores, since the profile would be cuspy over all astrophysically relevant radii yet not divergent.
  • The same slope-power-law form could be tested against the density profiles of halos in warm or self-interacting dark matter simulations, where a steeper or shallower inner behavior would show up as a different fitted exponent rather than a different functional form.
Original abstract

We investigate the mass profile of LambdaCDM halos using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a power law and depends little on mass. The logarithmic slope of the spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r), decreases monotonically towards the center and becomes shallower than isothermal (beta<2) inside a characteristic radius, r_{-2}. Although the fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual halos, systematic deviations from the best NFW fits are also noticeable. Inside r_{-2}, the profile of simulated halos gets shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. This discrepancy has been interpreted as indicating a steeply divergent cusp, but our results suggest a different interpretation. We use the density and enclosed mass at our innermost resolved radii to place strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are inconsistent with most of our simulations, although beta_{0}=1 is still consistent with our data. Our density profiles show no sign of converging to a well-defined asymptotic inner power law. We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-09-23: 1 request, 7,202 input tokens.

What the abstract itself states, as probabilities: claims first 0.03 · states limitation 0.89 · human participants 0.01 · controlled comparison 0.15 · randomized assignment 0.01 · states quantitative result 0.96 · compares with prior results 0.97 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v2. Results and the main claim are underlined; .

We investigate the mass profile of LambdaCDM halos using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a power law and depends little on mass. The logarithmic slope of the spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r), decreases monotonically towards the center and becomes shallower than isothermal (beta<2) inside a characteristic radius, r_{-2}. Although the fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual halos, systematic deviations from the best NFW fits are also noticeable. Inside r_{-2}, the profile of simulated halos gets shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. This discrepancy has been interpreted as indicating a steeply divergent cusp, but our results suggest a different interpretation. We use the density and enclosed mass at our innermost resolved radii to place strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are inconsistent with most of our simulations, although beta_{0}=1 is still consistent with our data. Our density profiles show no sign of converging to a well-defined asymptotic inner power law. We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 We investigate the mass profile of LambdaCDM halos using a suite of numerical simulations spanning five decades in halo mass, from dwarf galaxies to rich galaxy clusters. method 0.57 0.96 0.03 asserted 0.90
s1 Our analysis confirms the proposal of Navarro, Frenk & White (NFW) that the shape of LambdaCDM halo mass profiles differs strongly from a power law and depends little on mass. result 0.97 0.75 0.27 asserted 1.00
s2 The logarithmic slope of the spherically-averaged density profile, as measured by beta=-dln(rho)/dln(r), decreases monotonically towards the center and becomes shallower than isothermal (beta<2) inside a characteristic radius, r_{-2}. result 1.00 0.89 0.19 asserted 1.00
s3 Although the fitting formula proposed by NFW provides a reasonably good approximation to the density and circular velocity profiles of individual halos, systematic deviations from the best NFW fits are also noticeable. result 0.91 0.82 0.06 asserted 0.95
s4 Inside r_{-2}, the profile of simulated halos gets shallower with radius more gradually than predicted and, as a result, NFW fits tend to underestimate the dark matter density in these regions. result 0.99 0.88 0.10 asserted 1.00
s5 This discrepancy has been interpreted as indicating a steeply divergent cusp, but our results suggest a different interpretation. interpretation 0.90 0.22 0.09 hedged possibility 0.69
s6 We use the density and enclosed mass at our innermost resolved radii to place strong constraints on beta_{0}: density cusps as steep as r^{-1.5} are inconsistent with most of our simulations, although beta_{0}=1 is still consistent with our data. result 0.92 0.96 0.21 hedged possibility 0.68
s7 Our density profiles show no sign of converging to a well-defined asymptotic inner power law. result 0.99 0.95 0.43 asserted 1.00
s8 We propose a simple formula that reproduces the radial dependence of the slope better than the NFW profile, and so may minimize errors when extrapolating our results inward to radii not yet reliably probed by numerical simulations. contribution 0.96 0.92 0.17 hedged possibility 0.94

What kind of paper the classifier reports this is: simulation study 1.00 · evidence from simulation or computation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v2
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Spherically-averaged density profiles of all our simulated halos. Densities are computed in radial bins of equal logarithmic width and are shown from the innermost converged radius (rconv) out to about the virial radius of each halo (r200). Our simulations target halos in three distinct mass groups: “dwarf”, “galaxy”, and “cluster” halos. These groups span more than five decades in mass. Thick solid lines in the top panels illustrate the expected halo profile for each mass range according to the fitting formula proposed by NFW (top-left) or M99 (top-right). Bottom panels indicate…
Figure 2. Spherically-averaged circular velocity (Vc(r) = p GM(r)/r) profiles of all our simulated halos. As in
Figure 3. Logarithmic slope of the density profile of all halos in our sample, plotted versus radius. Thick solid and dotted curves illustrate the radial dependence of the slope expected from the NFW profile (eq. 1) and the modification proposed by M99 (eq. 2), respectively. Note that although both fitting formulae have well-defined asymptotic inner slopes (−1 and −1.5, respectively) there is no sign of convergence to a well-defined value of the central slope in the simulated halos. At the innermost converged radius, the simulated halo profiles are shallower than −1.5, in disagreement with…
Figure 4. Maximum asymptotic inner slope compatible with the mean density interior to radius r, ¯ρ(r), and with the local density at that radius, ρ(r). This provides a robust limit to the central slope, β0 < βmax(r) = −3(1 − ρ(r)/ρ¯(r)), under the plausible assumption that β is monotonic with radius. Note that there is not enough mass within the innermost converged radius in our simulations to support density cusps as steep as r −1.5 . The asymptotic slope of the NFW profile, β0 = 1, is still compatible with the simulated halos, although there is no convincing evidence for convergence to a…
Figure 5. (a-left) Density profiles of all halos in our series, scaled to the radius, r−2, where the local logarithmic slope of the density profile takes the isothermal value of β = −d log ρ/d log r = 2. Densities are scaled to ρ−2 = ρ(r−2). This figure shows that, with proper scaling, there is little difference in the shape of the density profile of halos of different mass, confirming the “universal” nature of the mass profile of ΛCDM halos. The NFW profile (eq. 1) is a fixed curve in these scaled units, and is shown with a thick solid line. The M99 formula (eq. 2) is shown with a dashed…
Figure 7. Comparison between the density (top) and circular velocity (bottom) profiles corresponding to four different fitting formulae: NFW (solid curves, eq. 1), M99 (dotted curves, eq. 2), SWTS (dashed curves, eq. 6), and ρα (dot-dashed curves, eq. 5). Circular velocity profiles are scaled to the maximum, Vmax, and to the radius where that is reached, rmax. Density profiles are scaled as in
Figure 8. The radius, r−2, where the logarithmic slope of the density profile takes the “isothermal” value, β(r−2) = 2, plotted versus the local density at that radius, ρ−2 = ρ(r−2), for all simulated halos in our series. This figure illustrates the mass dependence of the central concentration of dark matter halos: low mass halos are systematically denser than their more massive counterparts. Solid and dotted lines indicate the scale radius-characteristic density correlation predicted by the formalisms presented by Eke, Navarro & Steinmetz (2001) and Bullock et al (2001). These parameters may…

Showing 7 of the paper's 8 figures. See all figures in the PDF ↗ Figures reproduced from arXiv:astro-ph/0311231 by the authors.

References
How later papers used it

195 later papers cite it

Evidence: substantial. 178 of the citing papers share no author with it, and 24 of those carry a confident label from Pith's citation classifier.

Pith read 15 of the 248 citing sentences, choosing only those whose citation markers point at this paper alone, and wrote the paraphrases from them. The summary, the counts, the evidence level, the links and the quotes come from code; the summary counts the papers listed below. A later paper is listed under a label when Pith's citation classifier gives one of its sentences that label at probability 0.8 or more (0.95 for support), Pith's reading finds that the citing paper itself did it, and a second check confirms the paraphrase against the sentence. This section lists uses only. It does not list criticism or disagreement, so their absence here says nothing about how the paper has been received. Papers that share an author with this one are listed apart, under "By the same authors". A quote is the citing authors' own sentence, shown when their paper's licence is CC BY, CC BY-SA or CC0. JSON

Uses its method or tool 3 of 15 papers with this label · 1 by the same authors

By the same authors

  • Adopted an approach originally introduced in this paper to constrain the inner asymptotic slopes of dark matter density profiles.

    Micro galaxies in LCDM (2023) (same authors) · 1.00 · paraphrased: the licence does not allow quoting

Uses its data or values 3 of 9 papers with this label
Compares with it 1 of 3 papers with this label
How the literature uses it

86 arXiv papers cite this in 113 sentences

Each citing sentence is taken from the citing paper's LaTeX source and typed by Pith's classifier. A label appears at probability 0.8 or more (0.95 for support), a contest at 0.7 or more; the quoted words are the citing authors'. How each label was checked · JSON

Papers citing it · 87 sentences without a confident label

Contests it 1 paper
Uses its method or tool 3 papers
Uses its data or values 7 papers
Compares with it 1 paper
Cites it as background 14 papers

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.