REVIEW
Baryonic Imprints on DM Halos: the concentration-mass relation and its dependence on 28 model parameters in the CAMELS suite
Cosmology and Nongalactic Astrophysics · Submitted 2 Oct 2026 · 299 references
Pith's verdict: Solid, reproducible extension of S23's KLLR program from 6 to 28 parameters, with a genuinely useful three-component compression; the slope catalog is trustworthy, but the compression is validated only in-sample and the full weight matrix is missing, so the model deliverable needs completing before it is used as advertised.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The paper claims that the impact of 23 astrophysical parameters on the halo concentration–mass relation reduces to three principal components—cubic polynomials in halo mass—that capture roughly 95% of the variation in the measured slopes.
The reading
What carries the argument
Two tools carry the argument. The first is kernel-localized linear regression (KLLR), which fits the multivariate relation $\log_{10}c_{\rm vir} = \pi(M_{\rm vir},z) + \alpha_M(M_{\rm vir},z)\log_{10}M_{\rm vir} + \sum_X \alpha_X(M_{\rm vir},z)\,f_X(X)$, letting every slope vary with halo mass and redshift by regressing on kernel-weighted halo samples. The second is a principal-component decomposition of the 23 astrophysical slope functions, which separates the mass-dependent behaviors into three dominant modes; each mode is then approximated by a cubic polynomial in log mass (Equations 6–8). The deliverable is this three-component phenomenological model, which compresses roughly 95% of the baryonic imprint on $c_{\rm vir}$ into three one-dimensional curves.
Extended reading notes
Core claim
The paper's core claim has two parts. First, the sensitivity of the concentration–mass relation to each astrophysical parameter is itself strongly mass-dependent: SNe wind parameters dominate near Milky Way mass scales, AGN parameters matter mainly at higher masses, and several slopes show a zero crossing that the authors attribute to nonlinear SNe–AGN coupling. Second, and more important for practical use, the 23 astrophysical slope functions live in a low-dimensional space: three principal components account for roughly 95% of their variation, and each mode's mass dependence is well described by a cubic polynomial in $\mu = \log_{10}(M_{\rm vir}/M_{\rm piv})$ with pivot mass $M_{\rm piv}=5.6\times 10^{12}\,M_\odot/h$. The first mode is dominated by supernova wind parameters, the second and third by AGN parameters, and the same three cubics fit the slopes across $0<z<3$ even though they are constructed from $z=0$ data. The authors present these cubics as an adequate representation of the effective impact of astrophysics on halo concentration as determined by the TNG model.
Load-bearing premise
Equation (2) models $c_{\rm vir}$ as locally linear in each of the 28 parameters with no interaction terms, yet the paper itself argues that supernova and AGN feedback couple nonlinearly; if those couplings vary strongly inside the sampled parameter ranges, the recovered slopes are biased, and the three-component compression built on them inherits the bias.
Editorial extensions
If this is right
- Halo-model analyses can marginalize over unknown baryonic feedback by varying three fitted curves instead of 23 astrophysical parameters.
- Because the same three cubics fit the slopes across $0<z<3$, redshift-dependent predictions need only carry the three components, not the full parameter set.
- The mass at which each component acts—SNe near $10^{12}\,M_\odot/h$, AGN at group scales—gives a physical separation that isolates which feedback channel shapes the profile at each scale.
- Cosmological parameters, known to roughly 1%, contribute under 1% uncertainty to concentration, so the three-component model targets the dominant, astrophysics-driven uncertainty.
Reading between the lines
- The low dimensionality suggests baryonic feedback acts on the outer profile through a small number of physical channels (wind energetics, AGN thermal-versus-kinetic balance); a testable prediction is that the same three mode shapes reappear in other galaxy-formation models, with coefficients changing but shapes persisting.
- The no-interaction regression of Equation (2) is the load-bearing hinge: if SNe–AGN coupling is strong inside the sampled ranges, the recovered slopes—and the cubics built on them—are biased, and the compression would need interaction terms.
- A direct application is to recalibrate the cubics at higher resolution, where star formation is more efficient than at the CAMELS resolution, and check whether the mode shapes shift.
Original abstract
The evolution of baryons in halos, as well as their gravitational impact on the total-matter phase space, has a rich phenomenology and is vital in answering many questions in both cosmology and astrophysics. We use the CAMELS simulation suite to quantify the impact of twenty-three different galaxy formation parameters/processes (as well as five cosmological parameters) on the concentration--mass relation, $c_{\rm vir}$ - $M_{\rm vir}$, which characterizes the radial density distribution in halos. We construct a simulation-informed nonlinear model for concentration as a function of halo mass, redshift, and twenty-eight cosmological/astrophysical parameters. This is done using 2048 simulations of the IllustrisTNG galaxy formation model from the CAMELS suite. We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir} \in [10^{11}, 10^{14.5}]$ $M_{\odot}/h$ and in redshift $z \in [0, 3]$, finding clear mass- and redshift-dependent features in each of the different feedback mechanisms. From our measurements, we construct a three-component phenomenological model that can include astrophysical imprints in the concentration--mass relation in a fast and simple manner.
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.
The abstract, sentence by sentence
The evolution of baryons in halos, as well as their gravitational impact on the total-matter phase space, has a rich phenomenology and is vital in answering many questions in both cosmology and astrophysics. We use the CAMELS simulation suite to quantify the impact of twenty-three different galaxy formation parameters/processes (as well as five cosmological parameters) on the concentration--mass relation, $c_{\rm vir}$ - $M_{\rm vir}$, which characterizes the radial density distribution in halos. We construct a simulation-informed nonlinear model for concentration as a function of halo mass, redshift, and twenty-eight cosmological/astrophysical parameters. This is done using 2048 simulations of the IllustrisTNG galaxy formation model from the CAMELS suite. We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir} \in [10^{11}, 10^{14.5}]$ $M_{\odot}/h$ and in redshift $z \in [0, 3]$, finding clear mass- and redshift-dependent features in each of the different feedback mechanisms. From our measurements, we construct a three-component phenomenological model that can include astrophysical imprints in the concentration--mass relation in a fast and simple manner.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | The evolution of baryons in halos, as well as their gravitational impact on the total-matter phase space, has a rich phenomenology and is vital in answering many questions in both cosmology and astrophysics. | background 1.00 | 0.05 | 0.03 | asserted 1.00 |
| s1 | We use the CAMELS simulation suite to quantify the impact of twenty-three different galaxy formation parameters/processes (as well as five cosmological parameters) on the concentration--mass relation, $c_{\rm vir}$ - $M_{\rm vir}$, which characterizes the radial density distribution in halos. | contribution 0.65 | 0.96 | 0.09 | asserted 0.99 |
| s2 | We construct a simulation-informed nonlinear model for concentration as a function of halo mass, redshift, and twenty-eight cosmological/astrophysical parameters. | method 0.50 | 0.96 | 0.22 | asserted 1.00 |
| s3 | This is done using 2048 simulations of the IllustrisTNG galaxy formation model from the CAMELS suite. | method 1.00 | 0.67 | 0.03 | asserted 0.99 |
| s4 | We extract the imprints of galaxy formation across a wide range in mass $M_{\rm vir} \in [10^{11}, 10^{14.5}]$ $M_{\odot}/h$ and in redshift $z \in [0, 3]$, finding clear mass- and redshift-dependent features in each of the different feedback mechanisms. | result 0.97 | 0.97 | 0.18 | asserted 1.00 |
| s5 | From our measurements, we construct a three-component phenomenological model that can include astrophysical imprints in the concentration--mass relation in a fast and simple manner. | contribution 0.76 | 0.95 | 0.32 | asserted 1.00 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2610.03544 by the authors.
References
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.