REVIEW 4 major objections 5 minor 28 references
The Density Profile of Dynamical Halos
T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash
Pith's one-line read This paper claims that individual dynamical halos of fixed mass have orbiting density profiles controlled by a single variable, the halo radius $r_{\rm h}$, with halo radii scattered by about 16% at fixed mass and by about 11% after…
desk verdict Plausible but unproven: the claim that individual dynamical halo profiles have one degree of freedom needs independent validation and covariance handling. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the halo radius $r_{\rm h}$, the length scale at which the exponential truncation $\exp(-r^2/2r_{\rm h}^2)$ turns on in the adopted profile model [9]. With the inner scaling $\epsilon$ fixed at 0.037, the slope $\alpha_\infty$ is replaced by the deterministic relation $\Delta\alpha_\infty = 0.05 + \ln R$, where $R \equiv r_{\rm h}/r_{\rm h,st}(M_{\rm orb})$, turning the two-parameter fit into a one-parameter fit. Then $r_{\rm h}$ is linked to the relative formation time $a_{\rm RF} = a_{60}/{\rm median}(a_{60}|M_{\rm orb})$ by $r_{\rm h,mod} = (1.91-a_{\rm RF})\,r_{\rm h,st}(M_{\rm orb})$, which tightens the predicted radius. The slope–radius relation is what makes the profile a single-degree-of-freedom family; the formation-time relation is what makes the radius partly predictable from accretion history.
What would settle it
Fit every halo in an independent simulation with both $\alpha_\infty$ and $r_{\rm h}$ free, then refit with $\alpha_\infty$ forced by Equation (12); if the forced fits degrade significantly, or if the free fits scatter about $\Delta\alpha_\infty = 0.05 + \ln R$ much more widely than the quoted errors, the single-degree-of-freedom claim is not supported.
Extended reading notes
Core claim
The central discovery is that the orbiting profiles of dynamical halos form an effectively one-parameter family at fixed mass. Fitting each halo with the truncated power-law model $\rho_{\rm orb}(x) = A(x/\epsilon)^{-\alpha(x)}\exp(-x^2/2)$ yields two shape parameters, the asymptotic slope $\alpha_\infty$ and the halo radius $r_{\rm h}$; the paper shows these are tightly correlated within a mass bin, $\Delta\alpha_\infty = 0.05 + \ln(r_{\rm h}/r_{\rm h,st}(M_{\rm orb}))$. This correlation lets every halo's profile be described by $r_{\rm h}$ alone, with the slope imposed by Equation (12). The halo radius at fixed orbiting mass is lognormally distributed with intrinsic scatter $\sigma_{\ln r_{\rm h}|M_{\rm orb}} = 0.16$; late-forming halos are more compact, and using the relative formation time through Equation (18) reduces the scatter to $\sigma_{\ln r_{\rm h}|M_{\rm orb},a_{\rm RF}} = 0.11$. The paper takes these two scalings as evidence that the orbiting profile of a dynamical halo is a single-degree-of-freedom object whose extent encodes mass and formation history.
Load-bearing premise
The central assumption is that the tight correlation between fitted inner slope and halo radius, $\Delta\alpha_\infty = 0.05 + \ln R$, reflects a physical single degree of freedom and not a degeneracy in the two-parameter fitting procedure, because the same fits that define the relation are then refit under it.
Editorial extensions
If this is right
- A halo's orbiting density profile at fixed mass is fully characterized by one number, $r_{\rm h}$, so the profile shape carries no independent information beyond its spatial extent.
- From a measured profile slope and radius, Equations (11) and (18) can be inverted to estimate both the orbiting mass and the relative formation time of the halo.
- The intrinsic scatter in $r_{\rm h}$ at fixed mass is about 16%; adding formation time reduces it to about 11%, setting the accuracy with which mass, and mass plus formation time, predict the profile extent.
- Late-forming halos are more compact and have shallower inner slopes, extending the known concentration–formation-time connection to the orbiting profiles of dynamical halos.
- The same trends appear with an alternative four-parameter profile parameterization, suggesting the single-degree-of-freedom behavior is not an artifact of the adopted fitting function alone.
Reading between the lines
- If the relation $\Delta\alpha_\infty = 0.05 + \ln R$ holds across redshifts and simulation resolutions, the orbiting profile could serve as a low-cost estimator of halo mass and formation time that bypasses full merger-tree construction.
- Because the residual scatter after including formation time is still 11%, a natural next test is whether combining $a_{\rm RF}$ with the recent accretion rate removes more of the scatter than either variable alone; the paper compares but does not combine the two.
- If profile shape determines mass and formation time, then lensing or satellite-kinematics measurements of the inner slope may be able to infer halo radius, turning the single-degree-of-freedom claim into a route from observed profiles to halo assembly history.
- The 16% scatter in halo radius at fixed mass implies a corresponding spread in the truncation scale of the orbiting profile, which should be propagated into stacked halo-model predictions of large-scale structure.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper analyzes the dark matter density profiles of dynamical halos in a single N-body simulation. It fits individual halo profiles with the two-parameter model of Salazar et al., identifies a tight correlation between the fitted slope offset Δα∞ and the scaled halo radius at fixed mass, adopts this correlation as a one-parameter model, measures the scatter in halo radius, and studies its dependence on formation time and accretion rate. The central claims are that individual profiles at fixed mass have one effective degree of freedom, the halo radius, and that the profile shape can be used to read off both mass and formation time.
Significance. If the single-degree-of-freedom claim is established, this would be a useful empirical simplification of halo structure: the orbiting density profile of dynamical halos would be a one-parameter family at fixed mass, with a modest additional dependence on formation time. The paper is transparent about its methods, makes the analysis code publicly available, and includes a useful comparison with the independent Diemer parameterization. The main limitation is that the central correlation is calibrated and then imposed on the same data, so the paper currently falls short of demonstrating that the correlation is physical rather than a fitting degeneracy.
major comments (4)
- [III.B (Eq. 11), Fig. 2] The central claim that individual profiles have a single effective degree of freedom rests entirely on the tight Δα∞–ln R relation, but this relation is measured from the same two-parameter fits that define both variables, and the subsequent one-parameter refits in Section III.C impose Eq. (11) and therefore cannot validate it. The fixed value ε=0.037 places the slope transition at x∼0.037, below the fitted radial range r>0.09 h−1 Mpc, so the data constrain only the approach of α(x) to α∞ over the fitted interval; this is exactly the regime in which a degeneracy between rh and α∞ could masquerade as a physical correlation. I ask the authors to demonstrate, for example by fitting mock halos with independent rh and α∞, by cross-validating Eq. (11) on an independent subsample or simulation, or by reporting the individual-fit covariance and showing that its noise direction is not aligned with Eq. (11), that the correlation is not an artifact of the fitting procedure.
- [III.C (Eq. 14)] The intrinsic scatter σln rh|Morb = 0.16 is derived by fitting Var(ln rh/rh,st|N) = σ0^2 + k/N to the one-parameter refits, but the fit is not shown, no uncertainty on σ0 is quoted, and the variance model is not tested. Since these refits already assume Eq. (11), any systematic error in the assumed relation enters the quoted scatter. Please report the fitted k, the uncertainty on σ0, and validate the particle-noise decomposition with mocks or by checking stability across radial binning and the choice of δ in Eq. (5).
- [IV (Eqs. 17-18)] The formation-time relation is calibrated on the same high-mass sample used to define aRF via median(a60|Morb), and the result is summarized by the residual scatter σ = 0.11 without an uncertainty or residual diagnostics. The paper should show that the residuals are consistent with zero mean and constant variance as a function of mass and aRF, and quote an uncertainty on the 0.11 value. The abstract's wording that only a small fraction of the scatter is due to formation time should also be quantified relative to the variance: (0.16)^2 versus (0.11)^2 leaves roughly half the variance unexplained, which is not obviously a small fraction.
- [IV (Eq. 22)] The internal consistency check between Eq. (20) and Eq. (21) yields a 0.029 difference between the right-hand side (−0.094) and the directly measured value (−0.064). The attribution of this offset to non-Gaussian tails is plausible, but it should be verified, for instance by showing that the tail correction has the required sign and size. As it stands, the model has a small unmodeled offset in a relation that is used to support the formation-time dependence of the halo radius.
minor comments (5)
- [VI, bullet list] The second bullet says that 'at fixed mass, the slope of the profile (α∞) is tightly correlated with halo mass'; this appears to be a typo, as the text and Fig. 2 show a correlation with halo radius (or scaled radius R), not with halo mass.
- [VI, first comparison paragraph] The phrase 'characterized the orbiting orbiting profiles of spherical halos' contains a duplicated word and should be corrected.
- [Throughout] The manuscript mixes the spellings 'halos' and 'haloes'; please choose one and use it consistently.
- [III.C (Eq. 14)] The best-fit value of k in Eq. (14) is not reported, which makes the variance model difficult to reproduce; please add the fitted value and its uncertainty.
- [III.B (Eq. 10)] The uncertainties on the fitted line come from jackknifing the box, but the number of jackknife realizations and the scatter around the line are not stated; reporting these would help readers assess the tightness of the correlation in Fig. 2.
Circularity Check
No significant circularity: the central claim is an empirical correlation between independently fitted profile parameters, not a definitional identity or a fitted input renamed as a prediction.
full rationale
This paper is an empirical characterization rather than a derivation, and no step reduces by construction to its own inputs. The key relation Δα∞ = 0.05 + ln R (Eq. 11) is fit to the same individual two-parameter fits used to define Δα∞ and R, but those two parameters are independently fitted degrees of freedom in Eq. (1); neither is defined in terms of the other. The one-parameter refits in Sec. III.C impose Eq. (12), so they do not independently validate the relation, but the paper does not present those refits as the validation; the inference of a single effective degree of freedom rests on the small scatter about the fitted correlation in Fig. 2. The quoted scatters σln rh|Morb = 0.16 and σln rh|Morb,aRF = 0.11 are measured residual widths, not quantities forced to be small by the calibration, and the formation-time relation Eq. (18) is fit before the residual scatter is evaluated. The Salazar et al. [9] parameterization is adopted from prior work with overlapping authors, but it is an input model, not the target result, and the qualitative α∞–rh correlation is cross-checked against Diemer's [17] independent analysis. The concern that the correlation might reflect a fitting degeneracy rather than an intrinsic degree of freedom is a robustness/correctness question requiring covariance or mock tests; it is not a circularity of the paper's logic.
Assumptions & free parameters
free parameters (10)
- epsilon =
0.037
- Stacked mass-scaling parameters (rh,p, rh,s, aInf,p, aInf,s) =
840.3 h^-1 kpc, 0.226, 2.018, -0.050
- Delta-aInf intercept =
0.05
- Delta-aInf slope =
1.0
- Mean log-radius offset =
-0.064 (ratio 0.937)
- Intrinsic scatter at fixed mass =
0.16
- Formation-time relation intercept a0 =
1.91
- Formation-time relation slope sa =
-1.0
- Residual scatter at fixed mass and formation time =
0.11
- Cost-function percent error delta =
0.05
assumptions (6)
- domain assumption The Quijote simulation faithfully represents CDM halo structure at the resolved scales used.
- domain assumption The Garcia et al. orbiting/infalling classification correctly separates particles and produces accurate accretion times.
- domain assumption The Salazar et al. profile model with fixed epsilon = 0.037 is a valid description of individual halo orbiting profiles.
- ad hoc to paper The CDF level a60 = 0.60 and its median normalization define a useful formation-time proxy.
- ad hoc to paper The fitted aInf-ln R relation (Eq. 11) is exact enough to impose in one-parameter refits.
- domain assumption Rockstar/ConsistentTrees halo masses and accretion rates are accurate for the comparison in Figure 5.
Cite this review
Pith. "Pith review of The Density Profile of Dynamical Halos." pith.science (2026). https://pith.science/paper/56HNOAO3
@misc{pith2026250700410,
author = {Pith},
title = {Pith review of: The Density Profile of Dynamical Halos},
year = {2026},
howpublished = {\url{https://pith.science/paper/56HNOAO3}},
note = {Machine review of arXiv:2507.00410}
}
abstract
Among the most fundamental properties of a dark matter halo is its density profile. Motivated by the recent proposal by Garc\'ia et al. [R. Garc\'ia et. al., MNRAS 521, 2464 (2023)] to define a dynamical halo as the collection of orbiting particles in a gravitationally bound structure, we characterize the mean and scatter of the orbiting profile of dynamical halos as a function of their orbiting mass. We demonstrate that the orbiting profile of individual halos at fixed mass depends on a single dynamical variable -- the halo radius $r_{\rm h}$ -- which characterizes the spatial extent of the profile. The scatter in halo radius at fixed orbiting mass is $\approx 16\%$. Only a small fraction of this scatter arises due to differences in halo formation time, with late-forming halos being more compact (smaller halo radii). Accounting for this additional correlation results in an $\approx 11\%$ scatter in halo radius at fixed mass and halo formation time.
Figures
Figures from the paper (3 more)
Reference graph
Works this paper leans on
-
[1]
Diemer focuses on the properties of ρorb for spheri- cal overdensity halos, whereas we focus specifically on dynamical halos. Even though the difference in algorithms has a minimal impact on the orbit- ing profiles of individual halos [18], the choice of definition matters with regards to the scatter and correlation of halo properties with mass
-
[2]
in setting ϵ = 0.037, leaving only two free parameters with which to fit the orbiting profile of dynamical halos. Salazar et al. [9] fit the stacked density profile of dynam- ical halos, finding both α∞ and rh vary with mass as Table I. Maximum likelihood parameter values from Salazar et al. [9] from Equations (1)-(4). Parameter Description Value rh,p Hal...
work page 2025
-
[3]
J. Einasto, On the Construction of a Composite Model for the Galaxy and on the Determination of the System of Galactic Parameters, Trudy Astrofizicheskogo Instituta Alma-Ata 5, 87 (1965)
work page 1965
-
[4]
Diemer relies on the parametric model of Diemer [11], which takes the form ρorb(r) = ρs exp(S(r)) (23) where S(r) ≡ −2 α r rs α − 1 − 1 β " r rt β − rs rt β# .(24) We note that while the Diemer model has four degrees of freedom ( rt, rs, α, and β), when fitting individual halo profiles Diemer holds the parameters α and β fixed to their best-fit values in ...
-
[5]
S. More, B. Diemer, and A. V. Kravtsov, The Splashback- jacob.d.santiago2.mil@us.navy.mil Radius as a Physical Halo Boundary and the Growth of Halo Mass, Astrophys. J. 810, 36 (2015), arXiv:1504.05591 [astro-ph.CO]
arXiv 2015
-
[6]
R. Garc ´ ıa, E. Salazar, E. Rozo, S. Adhikari, H. Aung, B. Diemer, D. Nagai, and B. Wolfe, A better way to define dark matter haloes, MNRAS 521, 2464 (2023), arXiv:2207.11827 [astro-ph.CO]
arXiv 2023
-
[7]
J. F. Navarro, C. S. Frenk, and S. D. M. White, A Univer- sal Density Profile from Hierarchical Clustering, Astro- phys. J. 490, 493 (1997), arXiv:astro-ph/9611107 [astro- ph]
arXiv 1997
-
[8]
M. Fong and J. Han, A natural boundary of dark matter haloes revealed around the minimum bias and maximum infall locations, MNRAS 503, 4250 (2021), arXiv:2008.03477 [astro-ph.CO]
work page Pith review arXiv 2021
Show all 28 references
-
[9]
orbiting,
as our fiducial parameterization because it has fewer free parameters than that of Diemer [11]. We comment on how this choice impacts our results in Section V. Our paper is organized as follows: in Section II, we de- scribe our simulation data set, and how the particles near h...
2025 arXiv
-
[10]
Diemer and A
B. Diemer and A. V. Kravtsov, Dependence of the Outer Density Profiles of Halos on Their Mass Accretion Rate, Astrophys. J. 789, 1 (2014), arXiv:1401.1216 [astro- ph.CO]
2014 arXiv
-
[11]
Mansfield, A
P. Mansfield, A. V. Kravtsov, and B. Diemer, Splashback Shells of Cold Dark Matter Halos, Astrophys. J. 841, 34 (2017), arXiv:1612.01531 [astro-ph.CO]
2017 arXiv
-
[12]
Garc ´ ıa, E
R. Garc ´ ıa, E. Rozo, M. R. Becker, and S. More, A redefi- nition of the halo boundary leads to a simple yet accurate halo model of large-scale structure, MNRAS 505, 1195 (2021), arXiv:2006.12751 [astro-ph.CO]
2021 arXiv
-
[13]
Springel, The cosmological simulation code GADGET-2, MNRAS 364, 1105 (2005), arXiv:astro- ph/0505010 [astro-ph]
V. Springel, The cosmological simulation code GADGET-2, MNRAS 364, 1105 (2005), arXiv:astro- ph/0505010 [astro-ph]
2005
-
[14]
E. M. Salazar, E. Rozo, R. Garc ´ ıa, N. Kokron, S. Ad- hikari, B. Diemer, and C. Osinga, Dynamics-based halo model for large scale structure, Phys. Rev. D111, 043527 (2025), arXiv:2406.04054 [astro-ph.CO]
2025 arXiv
-
[15]
Diemer, A dynamics-based density profile for dark haloes - I
B. Diemer, A dynamics-based density profile for dark haloes - I. Algorithm and basic results, MNRAS 513, 573 (2022), arXiv:2112.03921 [astro-ph.CO]
2022 arXiv
-
[16]
Diemer, A dynamics-based density profile for dark haloes - II
B. Diemer, A dynamics-based density profile for dark haloes - II. Fitting function, MNRAS 519, 3292 (2023), arXiv:2205.03420 [astro-ph.CO]
2023 arXiv
-
[17]
read off
for details). Similar conclusions hold for the splash- back radius in simulations Shin and Diemer [22]. The right panel in Figure 5 shows that the variation in halo radius at fixed halo mass (i.e. R ≡ rh/rh,st) is corre- lated with accretion rate per unit mass as measured usin...
-
[18]
Villaescusa-Navarro, C
F. Villaescusa-Navarro, C. Hahn, E. Massara, A. Baner- jee, A. M. Delgado, D. K. Ramanah, T. Charnock, E. Giusarma, Y. Li, E. Allys, A. Brochard, C. Uhlemann, C.-T. Chiang, S. He, A. Pisani, A. Obuljen, Y. Feng, E. Castorina, G. Contardo, C. D. Kreisch, A. Nicola, J. Alsing, R...
2020 arXiv
-
[19]
Springel, N-GenIC: Cosmological structure initial conditions, Astrophysics Source Code Library, record ascl:1502.003 (2015)
V. Springel, N-GenIC: Cosmological structure initial conditions, Astrophysics Source Code Library, record ascl:1502.003 (2015)
2015
-
[20]
P. S. Behroozi, R. H. Wechsler, and H.-Y. Wu, The ROCKSTAR Phase-space Temporal Halo Finder and the Velocity Offsets of Cluster Cores, Astrophys. J. 762, 109 (2013), arXiv:1110.4372 [astro-ph.CO]. 9
2013 arXiv
-
[21]
R. H. Wechsler, J. S. Bullock, J. R. Primack, A. V. Kravtsov, and A. Dekel, Concentrations of Dark Halos from Their Assembly Histories, Astrophys. J. 568, 52 (2002), arXiv:astro-ph/0108151 [astro-ph]
2002 arXiv
-
[22]
Diemer, A dynamics-based density profile for dark haloes - III
B. Diemer, A dynamics-based density profile for dark haloes - III. Parameter space, MNRAS 536, 1718 (2025), arXiv:2410.17324 [astro-ph.CO]
2025 arXiv
-
[23]
Vladimir, C
Z. Vladimir, C. Osinga, B. Diemer, E. M. Salazar, and E. Rozo, Distinguishing Orbiting and Infalling Dark Matter Particles with Machine Learning, arXiv e-prints , arXiv:2506.09146 (2025), arXiv:2506.09146 [astro-ph.CO]
2025
-
[24]
Bertschinger, Self-similar secondary infall and ac- cretion in an Einstein-de Sitter universe, ApJS 58, 39 (1985)
E. Bertschinger, Self-similar secondary infall and ac- cretion in an Einstein-de Sitter universe, ApJS 58, 39 (1985)
1985
-
[25]
J. A. Fillmore and P. Goldreich, Self-similar gravitational collapse in an expanding universe, Astrophys. J. 281, 1 (1984)
1984
-
[26]
Shi, The outer profile of dark matter haloes: an analytical approach, MNRAS 459, 3711 (2016), arXiv:1603.01742 [astro-ph.CO]
X. Shi, The outer profile of dark matter haloes: an analytical approach, MNRAS 459, 3711 (2016), arXiv:1603.01742 [astro-ph.CO]
2016 arXiv
-
[27]
Shin and B
T.-h. Shin and B. Diemer, What sets the splashback ra- dius of dark matter haloes: accretion history or other properties?, MNRAS 521, 5570 (2023), arXiv:2210.14262 [astro-ph.CO]
2023 arXiv
-
[28]
P. S. Behroozi, R. H. Wechsler, H.-Y. Wu, M. T. Busha, A. A. Klypin, and J. R. Primack, Gravitationally Consis- tent Halo Catalogs and Merger Trees for Precision Cos- mology, Astrophys. J. 763, 18 (2013), arXiv:1110.4370 [astro-ph.CO]
2013 arXiv
Reviewed August 6, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.