Pith. sign in

REVIEW 3 major objections 6 minor 2 cited by

The Impact of Star Formation Histories on the Inner Dark Matter Density Slopes of Galaxies

T0 review · 3 major / 6 minor · reviewed 2026-08-03 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Star formation histories, not just galaxy mass, set whether dark matter centers are cored or cuspy.

desk verdict Plausible qualitative result, but the quantitative improvement claim needs out-of-sample validation and resolution checks. read the letter →

arxiv 2601.17113 v2 pith:MA7B77QG submitted 2026-01-23 astro-ph.GA

Jorge Sarrato-Alós , James Bullock , Arianna Di Cintio , Christopher Brook , Fernando Valenciano , Andrea V. Macciò This is my paper · ORCID
classification astro-ph.GA
keywords starformationhistoriesdarkmattercusp/coreinnerdensityslopegalaxysimulationsstellarfeedbackburstinessreionizationstellar-to-halomassratio
topics Dark Matter
open problems Dark Matter
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The reading

This paper argues that the shape of a galaxy's inner dark matter profile—whether it has a steep cusp or a flat core—depends on more than how much stellar mass it has relative to its halo. By comparing two large sets of simulated galaxies analysed with identical methods, the authors show that galaxies with burstier and more temporally extended star formation histories develop shallower dark matter slopes at fixed stellar-to-halo mass ratio. They distill this into a new two-variable fitting formula that predicts the inner slope more accurately than the standard one-parameter relation. If true, this resolves apparent discrepancies between simulation suites and gives observers a direct prediction: galaxies with old, early-peaked stellar populations should retain cusps, while those with extended, late star formation should host cores.

What carries the argument

The central objects are the burstiness metric—the fraction of stellar mass formed in starburst phases, defined by comparing star formation rates averaged over 50 and 500 Myr—and the post-to-pre reionization stellar mass ratio M_post/M_pre, which measures how temporally concentrated star formation is relative to the reionization epoch at z~6.5. Armed with these, the paper defines 'burst deviation' and 'SFH deviation' (residuals from the mass trends) and shows they order the scatter in the inner slope relation. The final fitting formula (Eq. 6) combines the stellar-to-halo mass ratio x with a power of y=M_post/M_pre to predict the inner slope.

What would settle it

A concrete test would be to recompute the same inner slopes from the same simulated galaxies at two or more resolution levels (or with different binning and centering choices) and check whether the burstiness and SFH-deviation trends survive; alternatively, an observational survey of dwarf galaxies with well-measured star formation histories and dark matter cores could check whether early-concentrated star formation systematically correlates with cusps, as the paper predicts. If the correlation vanishes or reverses, the claim is falsified.

Watch

Extended reading notes

Core claim

The central discovery is that the inner logarithmic slope of the dark matter density profile, measured between 1% and 2% of the virial radius, correlates with the star formation history in a way that is independent of the stellar-to-halo mass ratio. Galaxies whose star formation rate is burstier than average for their mass—quantified by a bursty mass fraction—and galaxies that formed a larger fraction of their stars after reionization (high M_post/M_pre) systematically have shallower (more cored) profiles. The paper packages this in a modified fitting formula (their Eq. 6) that adds a term proportional to a power of M_post/M_pre to the existing mass-ratio fit, cutting the mean squared predic

Load-bearing premise

The paper's central claim rests on the assumption that the inner slope measured by a linear fit to the density profile between 1% and 2% of the virial radius is a converged, resolution-independent measure of dark matter cusp/core structure; no convergence test is given, and Appendix B notes that for some FIRE-2 galaxies the analogous 150 pc radius lies below the resolution limit, requiring extrapolation.

Editorial extensions

If this is right

  • At fixed stellar-to-halo mass ratio, galaxies with burstier star formation develop shallower dark matter cores.
  • Galaxies with more extended SFHs (higher post-to-pre reionization ratio) are more efficient at creating cores; earlier-concentrated star formation preserves cusps.
  • The new two-variable fitting formula (Eq. 6) reduces the mean squared error of inner-slope prediction by roughly one-quarter in both simulation suites (FIRE-2: 0.073 to 0.055; NIHAO: 0.063 to 0.048).
  • The apparent disagreement between previous NIHAO and FIRE-2 inner-slope trends largely disappears when analysis is homogenised (same halo definition, binning), implying that much of the scatter is physical rather than numerical.
  • Observational prediction: among galaxies near the peak core-formation mass (M_star~10^8 Msun), those with early-peaked stellar age distributions should retain cusps, while extended, late star formation should correlate with large cores; finding large cores in early-forming dwarfs would challenge feedback-driven core formation.

Reading between the lines

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

  • The burstiness metric relies on a specific threshold (1.5 times the ratio of SFR averaged over 50 Myr vs 500 Myr); the paper does not test sensitivity to this choice, so a natural extension is to verify with other burst definitions or timescales.
  • If the trend holds, the scatter in observed core/cusp properties of dwarf galaxies could be used as a probe of their star formation histories, even where direct SFH measurement is difficult.
  • The correlation between burstiness and M_post/M_pre may mean the two metrics are not independent; disentangling their relative importance (e.g., with galaxies that are bursty but early-forming) could sharpen the causal claim.
  • The paper's exclusion of a small subsample lacking pre-reionization stars and one outlier in the fitting leaves open how universal the formula is for extremely low-mass galaxies.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

3 major / 6 minor

Summary. The paper analyzes 93 NIHAO and 109 FIRE-2 zoom-in galaxies, measuring inner dark-matter logarithmic slopes between 1% and 2% of R_vir and characterizing star formation histories via a bursty mass fraction (defined by SFR(50 Myr) > 1.5 SFR(500 Myr)) and the post-to-pre reionisation stellar mass ratio M_post/M_pre. It argues that, after homogenizing profile construction and overdensity definition, the two simulation suites agree much better than previously reported, and that residual burstiness and SFH concentration at fixed stellar mass correlate with cored versus cuspy inner profiles. It presents Eq. (6), which adds an M_post/M_pre correction to the standard inner-slope fitting formula, and reports lower mean squared errors in both suites. It also reproduces Muni et al.'s 150 pc density test for FIRE-2 and offers a falsifiable prediction for M* ~ 10^8 M_sun galaxies.

Significance. If the central claims hold, this is a valuable contribution: it demonstrates that previous NIHAO/FIRE-2 differences were partly due to differing overdensity definitions and binning, it identifies a plausible second parameter (burstiness / SFH timing) that reduces scatter in the inner-slope relation, and it gives a concrete observational prediction tied to resolved stellar populations. The homogeneous re-analysis and the explicit reproduction of Muni et al.'s EDGE test are strengths. However, the quantitative claims are not yet established: the inner-slope measurement lacks resolution validation, and the reported improvement of Eq. (6) is in-sample and, for NIHAO, effectively a constant offset rather than an SFH-dependent term. The qualitative correlation between bursty/extended SFHs and cores is plausible, but the paper currently overstates the predictive gain.

major comments (3)
  1. [§3, Figs. 1/3/5, Appendix B] The inner slope is a linear fit to the logarithmic density profile between 0.01 and 0.02 R_vir, using 35 bins down to 0.005 R_vir. The paper never establishes that this window is resolved for every galaxy: there is no requirement such as r_inner > r_200 or > 3× softening, no convergence test, and no per-galaxy slope uncertainties. The sample spans M* ~ 10^5–10^11 M_sun, so 0.01 R_vir ranges from tens of pc in dwarfs to kpc in massive halos. Appendix B admits that for some FIRE-2 galaxies even 150 pc lies below the r_200 resolution limit and requires extrapolation; the analogous risk for the 1–2% R_vir window is not assessed. If unresolved or marginally resolved systems are biased toward shallower slopes, and if the resolved fraction correlates with M*/Mvir or with SFH (low-mass dwarfs being the burstiest), then the split-sample trends in Figs. 3 and 5 and the MSE gains in §3.3 could be n
  2. [§3.3, Eq. (6), Table 4] The claimed improvement in prediction is in-sample: Eq. (6) is fitted and evaluated on the same galaxies, so a lower MSE is expected merely from adding parameters (n2, α) and removing an outlier. No cross-validation or held-out sample is presented. More seriously, the NIHAO fit gives α = 2.06×10^-3; over any plausible M_post/M_pre range, y^α ≈ 1, so the correction term is nearly a constant offset. The NIHAO MSE reduction therefore cannot be attributed to an SFH dependence of the slope. The FIRE-2 value α = 0.14 is also weak, and neither fit reports parameter uncertainties or covariance. Please add out-of-sample validation, report the actual dynamic range of y^α, and show the sensitivity of the results to the removed outlier.
  3. [§3.1–3.2, Figs. 3 and 5] The 'burst deviation' and 'SFH deviation' are residuals from second-order polynomial fits to the same galaxies that are then split into positive/negative deviation groups for slope fits. This construction guarantees that the two groups differ by construction; it can create apparent secondary correlations when the underlying relation already has mass-dependent scatter. To support the claim that burstiness and SFH concentration explain scatter at fixed M*/Mvir, the authors should either use independent variables (e.g., direct SFR metrics not de-trended against M*) or demonstrate with cross-validation that the split-sample fits predict slopes of galaxies not used in the polynomial fit. As written, the explanatory claim is partly in-sample.
minor comments (6)
  1. [§4] Duplicate word in 'the temporal temporal concentration of SFHs'.
  2. [Appendix B] 'median relative deviation over 0.15%' conflicts with the earlier 'values above ~0.15'; if 0.15 is a fraction, it should be written as 0.15 or 15%.
  3. [Fig. 7] No indication of the distribution of y=M_post/M_pre is given. Please show the range and percentiles for each suite, otherwise the plotted curves may suggest extrapolation beyond the data.
  4. [Table 2] The FIRE-2 (<0) fit has β=348, clearly an unconstrained/degenerate parameter. Report parameter uncertainties or use a more robust fitting procedure for the split samples.
  5. [Eq. (6)] The equation display does not define x and y; x=M*/Mvir and y=M_post/M_pre are given in the text but should be repeated in the equation caption or below the display.
  6. [Appendix A] The sentence about the NIHAO peak value rounding to the same value at one significant digit is confusing; report more significant digits.

Circularity Check

1 steps flagged · score 6.0 of 10

Eq. 6's advertised 'improved predictions' are in-sample fits: adding the fitted term n2·y^α to the same galaxies guarantees the reported MSE drop, so the prediction claim is statistically forced.

  1. fitted input called prediction [Section 3.3, Eq. 6 and Table 4]
    "We report a two variable fitting formula for the inner slope described by Eq. 6, where x=M⋆/Mvir and y=M⋆,post/M⋆,pre. The inclusion of the additional correction lowers the mean squared error in the prediction of the inner slope from 0.073 to 0.055 for FIRE-2 galaxies, and from 0.063 to 0.048 for NIHAO galaxies."

    The MSE values are computed on the same NIHAO and FIRE-2 galaxies used to fit the parameters of Eq. 6 (Table 4: n, n1, n2, x1, x0, β, γ, α). Since Eq. 6 contains Eq. 2 as the n2=0 special case and all parameters are jointly fitted to these data, the in-sample mean squared error cannot increase when the extra fitted term n2·y^α is added. The reported reduction is therefore a guaranteed property of adding a fitted parameter, not evidence of out-of-sample predictive improvement. The paper calls this 'prediction of the inner slope' and 'higher accuracy,' but no holdout, cross-validation, or AIC/BIC-type penalty is described, so the improvement is a training-set gain rather than a validated prediction.

full rationale

The core SFH–slope correlations in Figs. 3 and 5 are not definitionally circular: the inner slope is measured from density profiles, while burstiness and M⋆,post/M⋆,pre are constructed from SFR histories, and the target variable is not used to define these metrics. The baseline Eq. 2 fits are refitted to the data, not assumed, so the self-citations to Di Cintio et al. (2014b), Tollet et al. (2016), and Lazar et al. (2020) are not load-bearing circularity. The main circular step is in Section 3.3, where Eq. 6 is fitted to the same galaxies whose slopes it is then said to 'predict'; with the additional fitted term, the reported MSE reduction is an in-sample artifact. The paper also defines 'burst deviation' and 'SFH deviation' as residuals from polynomial fits to the same data, making the subsequent 'secondary predictors' language partly a two-stage descriptive fit, though the deviations are not derived from the inner slope itself. The unresolved-slope and resolution concerns in Appendix B are robustness/correctness issues, not circularity. Overall, the qualitative SFH dependence is not forced by construction, but the headline quantitative claim of improved predictive accuracy is statistically forced, giving a partial circularity score of 6.

Assumptions & free parameters 5 free parameters · 4 assumptions · 0 invented entities

The paper introduces several hand-chosen thresholds (burstiness ratio, timescales) and ad hoc functional forms for defining deviations and the final correction term. No new physical entities are postulated. The assumptions are typical for simulation-based studies, but the lack of convergence or error-bar analysis makes the slope measurement assumptions load-bearing.

free parameters (5)
  • burstiness threshold ratio = 1.5
    Eq. 3 defines a burst when SFR(50 Myr) > 1.5 * SFR(500 Myr); the 1.5 factor is chosen without derivation, adapting Sparre et al. (2017).
  • SFR averaging timescales = 50 Myr and 500 Myr
    Eq. 3 uses 50 and 500 Myr windows, longer than Sparre et al.'s 10/200 Myr, with no tuning justification.
  • polynomial degree for f_M,burst vs M* = 2
    Section 3.1 fits a second-degree polynomial to define burst deviation; no model comparison.
  • polynomial degree for M_post/M_pre vs M* = 2
    Section 3.2 fits a second-degree polynomial to define SFH deviation.
  • functional form of Eq. 6 correction = n2 * y^alpha
    The additive power-law correction is chosen ad hoc; in the NIHAO fit alpha ≈ 2e-3, making the term nearly constant, so the improvement is not from SFH variation.
assumptions (4)
  • domain assumption NIHAO and FIRE-2 subgrid feedback models faithfully reproduce the physical core-forming processes
    The entire comparison rests on these simulations being valid; the paper cites their prior validation but does not test convergence.
  • domain assumption The inner slope at 1-2% R_vir is a robust measure of cusp/core
    Used throughout; no convergence tests or error bars are given for this statistic.
  • domain assumption A single reionization epoch at z=6.5 applies to all galaxies
    M_post/M_pre is defined with treion=0.84 Gyr for every galaxy, ignoring patchy reionization.
  • domain assumption Burstiness and M_post/M_pre deviations from polynomial fits are independent secondary parameters
    The deviations are residuals from fits that include correlated scatter; no test ensures they are not proxies for other variables.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Impact of Star Formation Histories on the Inner Dark Matter Density Slopes of Galaxies." pith.science (2026). https://pith.science/paper/MA7B77QG

@misc{pith2026260117113,
  author       = {Pith},
  title        = {Pith review of: The Impact of Star Formation Histories on the Inner Dark Matter Density Slopes of Galaxies},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MA7B77QG}},
  note         = {Machine review of arXiv:2601.17113}
}
abstract

Aims. We aim to investigate the connection between star formation histories (SFHs) and the inner dark matter density profiles of simulated galaxies. In particular, we test whether the burstiness and temporal distribution of star formation influence the formation of cored versus cuspy dark matter profiles. Methods. We homogeneously analysed simulated galaxies from the NIHAO and FIRE-2 projects. For each galaxy, we derived dark matter density profiles and measured the logarithmic slope in the inner region of the dark matter halo (1-2% of R$_{\rm vir}$). To characterise star formation burstiness, we introduced a criterion based on comparing the star formation rate (SFR) averaged over two distinct timescales. We further quantified the temporal concentration of SFHs by computing $M_{\star, \rm post}$ / $M_{\star, \rm pre}$, the ratio of stellar mass formed after versus before the epoch of reionisation at redshift z $\sim$ 6.5. Results. Homogeneous analysis reveals that inner slope versus stellar-to-halo mass ratio trends for NIHAO and FIRE-2 galaxies are in much better agreement than reported in previous works. The burstiness and post-to-pre reionisation stellar mass ratio of the SFH explain the scatter in the inner slope versus stellar-to-halo mass ratio relation, revealing that galaxies with above average burstiness and more extended SFHs are more efficient at developing cored dark matter profiles. In contrast, galaxies with smoother SFHs and earlier stellar mass assembly tend to maintain cuspier dark matter profiles. We present an analytic expression that improves predictions for the inner slope using the parameter $M_{\star \rm,post}$ / $M_{\star \rm,pre}$, which reduces the mean squared error in both simulation suites relative to previous formulations based solely on the stellar-to-halo mass ratio.

Figures

Figures reproduced from arXiv: 2601.17113 by the authors.

Figure 1
Figure 1. Inner slope of the dark matter density profile, measured between 1% and 2% of the virial radius, as a function of the stellar-to-halo mass ratio. Results are shown for galaxies from the NIHAO (blue) and FIRE￾2 (orange) simulations. Solid lines represent fits following Eq. 2 with parameters from [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. Inner slope of the dark matter density profile, measured between 1% and 2% of the virial radius, as a function of the stellar-to-halo mass ratio. The top panel shows results for NIHAO galaxies and the bot￾tom panel shows results for FIRE-2 galaxies. Circles for each galaxy are colour coded by their burst deviation (the difference between actual and expected bursty mass fraction at a given stellar mass). Separate fit… view at source ↗
Figure 5
Figure 5. , with fit parameters for Eq. 2 listed in [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figures from the paper (1 more)
Figure 6
Figure 6. Figure 6: Relation between M⋆,post / M⋆,pre and the bursty mass fraction for NIHAO and FIRE-2 galaxies. The Pearson correlation coefficient is indicated in the legend. One outlier NIHAO galaxy, marked with an asterisk, is removed from the calculation. We report a two variable fi…

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Dark matter density profiles of the Milky Way satellite population: reconciling simulations and observations

    astro-ph.GA 2026-07 conditional novelty 6.0 of 10

    Fixed-physical-radius slope and full-profile comparisons of NIHAO/FIRE-2 satellites to 16 MW dwarfs show mass-dependent core formation and agreement within ~1σ for most systems, outperforming NFW.

  2. It's Not Just Star Formation: A trend of low dark matter densities in the Andromeda dwarf galaxy system

    astro-ph.GA 2026-05 unverdicted novelty 4.0 of 10

    Five of seven modeled M31 dwarf spheroidals show anomalously low central DM densities at 150 pc, with star formation heating disfavored as the sole cause.

Reference graph

Works this paper leans on

4 extracted references · cited by 2 Pith papers

  1. [1]

    I., et al

    Agertz, O., Pontzen, A., Read, J. I., et al. 2019, MNRAS, 491, 1656 Brook, C. B. & Di Cintio, A. 2015, MNRAS, 453, 2133 Brook, C. B., Stinson, G., Gibson, B. K., Wadsley, J., & Quinn, T. 2012, MN- RAS, 424, 1275 Bryan, G. L. & Norman, M. L. 1998, ApJ, 495, 80 Chabrier, G. 2003, PASP, 115, 763 Chan, T. K., Kereš, D., Wetzel, A., et al. 2018, MNRAS, 478, 90...

    2019

  2. [4]

    Nevertheless, we confirm via visual inspection that the inner slope of the fitted profile gener- ally matches the data

    We note that the fixed core radius does not provide satisfactory fits for some galaxies in our sample, which span a broader range of stellar and halo proper- ties than the EDGE simulations. Nevertheless, we confirm via visual inspection that the inner slope of the fitted profile gener- ally matches the data. To quantify the fit quality, we compute the med...

    2024

  3. [2016]

    To facilitate a direct comparison to these works, we repeat our analysis using∆ =200c

    use a fixed overdensity of∆ =200c to investigate the relation between the inner slope of the dark matter density profile and the stellar-to-halo mass ra- tio. To facilitate a direct comparison to these works, we repeat our analysis using∆ =200c. For this purpose, we construct new density profiles analogous to those presented in Fig. 2 of Tollet et al. (20...

    2016

  4. [2025]

    They reported a tight, decreasing linear relation between the inner dark matter density and the ratioM ⋆,post/M⋆,pre

    by measuring the dark matter density at 150 pc from galaxy cen- ters. They reported a tight, decreasing linear relation between the inner dark matter density and the ratioM ⋆,post/M⋆,pre. We repli- cate this analysis using FIRE-2 galaxies to calculate their density within a spherical shell ranging from 125 to 175 pc, with results shown in Fig. B.1. NIHAO ...

    2003

Pith tools

Reviewed August 3, 2026 · model on record in the stance chip above.