Stellar halos of bright central galaxies II: Scaling relations, colors and metallicity evolution with redshift
Pith reviewed 2026-05-18 02:07 UTC · model grok-4.3
The pith
Stellar halos act as transition regions between bright central galaxies and intracluster light.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Stellar halos emerge as transition regions between bright central galaxies and the intracluster light. They are dynamically and chemically coupled to both, with their properties depending on halo concentration, intracluster light formation efficiency, and the progenitor mass spectrum. The stellar halo mass correlates strongly with both the bright central galaxy and intracluster light masses, with tighter scatter in the stellar halo-intracluster light relation. Stellar halos and intracluster light exhibit nearly identical color distributions that redden toward the present day, while the metallicity gap between stellar halos and bright central galaxies narrows over cosmic time.
What carries the argument
The transition radius, defined as the stellar component boundary linked to dark matter halo concentration, which separates stellar halos from intracluster light and tracks their coupled mass, color, and metallicity evolution in the updated FEGA25 model.
If this is right
- Stellar halo mass scales more tightly with intracluster light mass than with bright central galaxy mass.
- Stellar halos and intracluster light maintain nearly identical colors across all redshifts and both redden toward the present.
- The metallicity difference between stellar halos and bright central galaxies decreases from about 0.4 dex at redshift 2 to 0.1 dex at redshift 0.
- Transition radii typically peak near 30 to 40 kpc but can extend to hundreds of kiloparsecs in the most massive halos after redshift 0.5.
Where Pith is reading between the lines
- Stellar halo observations could indirectly constrain dark matter halo concentrations in clusters through the transition radius relation.
- The metallicity mismatch implies that future models need stronger contributions from low-mass disrupted galaxies to match data.
- Kinematic mapping in large samples could test whether the dynamical coupling between stellar halos, central galaxies, and intracluster light holds in real clusters.
Load-bearing premise
Stellar halos are correctly identified as the stellar material inside a transition radius set by halo concentration, and the model's updated treatment of intracluster light formation captures the dominant physical processes.
What would settle it
Deep imaging and spectroscopy of many galaxy clusters that show stellar halo colors or metallicities diverging strongly from those of the intracluster light, or that find no link between transition radius and halo concentration.
Figures
read the original abstract
We study the formation and evolution of stellar halos (SHs) around bright central galaxies (BCGs), focusing on their scaling relations, colors, and metallicities across cosmic time, and compare model predictions with ultra--deep imaging data. We use the semianalytic model \textsc{FEGA25}, applied to merger trees from high--resolution dark matter simulations, including an updated treatment of intracluster light (ICL) formation. SHs are defined as the stellar component within the transition radius, linked to halo concentration. Predictions are compared with observations from the VST Early-type GAlaxy Survey (VEGAS) and Fornax Deep Survey (FDS). The SH mass correlates strongly with both BCG and ICL masses, with tighter scatter in the SH--ICL relation. The transition radius peaks at 30--40 kpc nearly independent of redshift, but can reach $\sim400$ kpc in the most massive halos, after z=0.5. SHs and ICL show nearly identical color distributions at all epochs, both reddening toward $z=0$. At $z=2$, SHs and the ICL are $\sim0.4$ dex more metal--poor than BCGs, but the gap shrinks to $\sim0.1$ dex by the present time. Observed colors are consistent with model predictions, while observed metallicities are lower, suggesting a larger contribution from disrupted dwarfs. SHs emerge as transition regions between BCGs and the ICL, dynamically and chemically coupled to both. Their properties depend on halo concentration, ICL formation efficiency, and the progenitor mass spectrum. Upcoming wide--field photometric and spectroscopic surveys (e.g. LSST, WEAVE, 4MOST) will provide crucial tests by mapping structure, metallicity, and kinematics in large galaxy samples.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper uses the semianalytic model FEGA25 applied to high-resolution dark matter merger trees, with an updated ICL formation treatment, to study stellar halos (SHs) of bright central galaxies (BCGs). SHs are defined as the stellar component inside a transition radius tied to halo concentration. The work reports strong SH mass correlations with BCG and ICL masses (tighter for SH-ICL), a transition radius peaking at 30-40 kpc independent of redshift (extending to ~400 kpc in massive halos post-z=0.5), nearly identical SH and ICL color distributions that redden toward z=0, and a metallicity gap shrinking from ~0.4 dex at z=2 to ~0.1 dex at z=0. Model colors match VEGAS and FDS observations while metallicities are overpredicted, implying greater dwarf disruption; the authors conclude SHs are dynamically and chemically coupled transition regions whose properties depend on halo concentration, ICL efficiency, and progenitor mass spectrum, with forecasts for LSST, WEAVE, and 4MOST.
Significance. If the central trends hold, the results advance understanding of BCG, SH, and ICL co-evolution by linking observable scaling relations and chemical properties to halo concentration and merger history across cosmic time. The direct comparison to ultra-deep VEGAS and FDS photometry provides a useful observational anchor, especially for colors, while the noted metallicity offset highlights the role of low-mass progenitors. The redshift-independent transition radius and survey predictions constitute falsifiable outputs that can be tested with upcoming wide-field data.
major comments (2)
- [§2] §2 (model description): The SH definition as the stellar component within the transition radius, which is itself tied to halo concentration from the SAM, renders the claim of dynamical and chemical coupling between SHs, BCGs, and ICL dependent on the internal partitioning of FEGA25; without explicit cross-validation against hydrodynamical simulations that independently resolve stripping and merger debris, the reported color similarity and shrinking metallicity gap could partly reflect shared model assumptions rather than physical coupling.
- [Abstract and results on metallicity evolution] Abstract and results on metallicity evolution: The model predicts SHs and ICL ~0.4 dex more metal-poor than BCGs at z=2 shrinking to 0.1 dex at z=0, yet observed metallicities are lower than predicted; this offset is load-bearing for the stated dependence on the progenitor mass spectrum, but the manuscript does not quantify how varying the ICL formation efficiency parameters would alter the gap or the SH-ICL correlation tightness.
minor comments (2)
- Figure captions should explicitly state the redshift bins and sample selection criteria used for the color and metallicity histograms to allow direct reproduction of the reported agreement with VEGAS/FDS.
- The manuscript would benefit from a short table summarizing the free parameters in the updated ICL formation prescription and their adopted values.
Simulated Author's Rebuttal
We thank the referee for the constructive review and the recommendation for minor revision. The comments raise valid points about model dependence and parameter sensitivity that we address below. We will incorporate targeted revisions to strengthen the manuscript without altering its core conclusions.
read point-by-point responses
-
Referee: §2 (model description): The SH definition as the stellar component within the transition radius, which is itself tied to halo concentration from the SAM, renders the claim of dynamical and chemical coupling between SHs, BCGs, and ICL dependent on the internal partitioning of FEGA25; without explicit cross-validation against hydrodynamical simulations that independently resolve stripping and merger debris, the reported color similarity and shrinking metallicity gap could partly reflect shared model assumptions rather than physical coupling.
Authors: We agree that the SH definition relies on the transition radius computed from halo concentration within FEGA25, making the reported coupling inherently tied to the model's partitioning scheme. This choice is physically motivated by the expected transition between the central galaxy and ICL regimes and is consistent with observational surface-brightness-based definitions. While direct cross-validation against hydrodynamical simulations is not included here, the close match between predicted and observed colors from VEGAS and FDS provides empirical support for the physical relevance of the results. The metallicity evolution follows from the shared merger history and progenitor mass spectrum in the underlying dark-matter trees. In the revised manuscript we will expand the discussion in §2 to explicitly note these model assumptions, their potential influence on the reported similarities, and the value of future hydrodynamical comparisons. revision: partial
-
Referee: Abstract and results on metallicity evolution: The model predicts SHs and ICL ~0.4 dex more metal-poor than BCGs at z=2 shrinking to 0.1 dex at z=0, yet observed metallicities are lower than predicted; this offset is load-bearing for the stated dependence on the progenitor mass spectrum, but the manuscript does not quantify how varying the ICL formation efficiency parameters would alter the gap or the SH-ICL correlation tightness.
Authors: We concur that a quantitative exploration of how ICL formation efficiency parameters affect the metallicity gap and the tightness of the SH–ICL relation would better substantiate the link to the progenitor mass spectrum. The current work presents results for the fiducial model calibrated to observations. We will add a short sensitivity analysis (new subsection or appendix) that varies the ICL efficiency within the range explored during model calibration and shows the resulting changes to the gap and correlation scatter. This will demonstrate that the main evolutionary trends remain robust. revision: yes
Circularity Check
No significant circularity: predictions derived from external merger trees and SAM
full rationale
The paper applies the FEGA25 semianalytic model to merger trees extracted from high-resolution dark matter simulations, defines SHs as the stellar component inside a transition radius tied to halo concentration, and generates predictions for scaling relations, colors, and metallicities that are compared against independent observational datasets (VEGAS, FDS). No quoted step shows a prediction reducing by construction to a fitted input, a self-citation chain, or a definitional tautology; the central claims about dynamical/chemical coupling emerge as model outputs rather than inputs. This is the normal self-contained case for simulation-based work.
Axiom & Free-Parameter Ledger
free parameters (1)
- ICL formation efficiency parameters
axioms (2)
- domain assumption Merger trees from high-resolution dark matter simulations accurately represent the assembly history of halos.
- domain assumption Stellar halos are the stellar component inside the transition radius tied to halo concentration.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
SHs are defined as the stellar component within the transition radius, linked to halo concentration... R_trans = R_200 / c_ICL (Eq. 2)
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
ICL arises through stellar stripping, mergers, and pre-processing... SH forms directly from stars originally belonging to the ICL
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Amorisco, N. C. 2017, MNRAS, 464, 2882 Bahé, Y . M., Barnes, D. J., Dalla Vecchia, C., et al. 2017, MNRAS, 470, 4186
work page 2017
- [2]
-
[3]
C., Carollo, D., Ivezi´c, Ž., et al
Beers, T. C., Carollo, D., Ivezi´c, Ž., et al. 2012, ApJ, 746, 34
work page 2012
-
[4]
Belokurov, V ., Erkal, D., Evans, N. W., Koposov, S. E., & Deason, A. J. 2018, MNRAS, 478, 611
work page 2018
-
[5]
Belokurov, V ., Sanders, J. L., Fattahi, A., et al. 2020, MNRAS, 494, 3880
work page 2020
-
[6]
Beltrand, C., Monachesi, A., D’Souza, R., et al. 2024, A&A, 690, A115
work page 2024
-
[7]
Bullock, J. S. & Johnston, K. V . 2005, ApJ, 635, 931
work page 2005
-
[8]
Capaccioli, M., Spavone, M., Grado, A., et al. 2015, A&A, 581, A10
work page 2015
-
[9]
Carollo, D., Beers, T. C., Chiba, M., et al. 2010, ApJ, 712, 692
work page 2010
- [10]
- [11]
-
[12]
L., Habib, S., Heitmann, K., et al
Child, H. L., Habib, S., Heitmann, K., et al. 2018, ApJ, 859, 55
work page 2018
-
[13]
Contini, E., De Lucia, G., Villalobos, Á., & Borgani, S. 2014, MNRAS, 437, 3787
work page 2014
-
[14]
Contini, E. & Gu, Q. 2020, ApJ, 901, 128
work page 2020
-
[15]
Contini, E., Jeon, S., Rhee, J., Han, S., & Yi, S. K. 2023, ApJ, 958, 72
work page 2023
-
[16]
Contini, E., Yi, S. K., & Jeon, S. 2024c, arXiv e-prints, arXiv:2404.01560
- [17]
- [18]
-
[19]
Contreras-Santos, A., Knebe, A., Cui, W., et al. 2024, A&A, 683, A59
work page 2024
-
[20]
Cooper, A. P., Cole, S., Frenk, C. S., et al. 2010, MNRAS, 406, 744
work page 2010
-
[21]
P., Martínez-Delgado, D., Helly, J., et al
Cooper, A. P., Martínez-Delgado, D., Helly, J., et al. 2011, ApJ, 743, L21
work page 2011
-
[22]
Cooper, A. P., Parry, O. H., Lowing, B., Cole, S., & Frenk, C. 2015, MNRAS, 454, 3185
work page 2015
-
[23]
Correa, C. A., Wyithe, J. S. B., Schaye, J., & Duffy, A. R. 2015, MNRAS, 452, 1217
work page 2015
-
[24]
Crain, R. A., Schaye, J., Bower, R. G., et al. 2015, MNRAS, 450, 1937
work page 2015
-
[25]
2018, MNRAS, 480, 2898 de Jong, R
Cui, W., Knebe, A., Yepes, G., et al. 2018, MNRAS, 480, 2898 de Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, The Messenger, 175, 3 de Vaucouleurs, G. 1959, Handbuch der Physik, 53, 275
work page 2018
-
[26]
J., Belokurov, V ., & Evans, N
Deason, A. J., Belokurov, V ., & Evans, N. W. 2013, ApJ, 763, 113
work page 2013
-
[27]
J., Belokurov, V ., Koposov, S
Deason, A. J., Belokurov, V ., Koposov, S. E., & Rockosi, C. M. 2014, ApJ, 787, 30
work page 2014
-
[28]
J., Mao, Y .-Y ., & Wechsler, R
Deason, A. J., Mao, Y .-Y ., & Wechsler, R. H. 2016, ApJ, 821, 5
work page 2016
-
[29]
Drinkwater, M. J., Gregg, M. D., & Colless, M. 2001, ApJ, 548, L139 D’Souza, R. & Bell, E. F. 2018, MNRAS, 474, 5300
work page 2001
-
[30]
Duc, P.-A., Cuillandre, J.-C., Karabal, E., et al. 2015, MNRAS, 446, 120
work page 2015
- [31]
-
[32]
Elias, L. M., Sales, L. V ., Helmi, A., & Hernquist, L. 2020, MNRAS, 495, 29
work page 2020
-
[33]
Font, A. S., McCarthy, I. G., Crain, R. A., et al. 2011, MNRAS, 416, 2802
work page 2011
-
[34]
Gallart, C., Bernard, E. J., Brook, C. B., et al. 2019, Nature Astronomy, 3, 932
work page 2019
- [35]
-
[36]
Gilbert, K. M., Kalirai, J. S., Guhathakurta, P., et al. 2014, ApJ, 796, 76
work page 2014
- [37]
-
[38]
Grand, R. J. J., Gómez, F. A., Marinacci, F., et al. 2017, MNRAS, 467, 179
work page 2017
- [39]
- [40]
- [41]
- [42]
-
[43]
F., Kereš, D., Oñorbe, J., et al
Hopkins, P. F., Kereš, D., Oñorbe, J., et al. 2014, MNRAS, 445, 581
work page 2014
-
[44]
Ibata, R. A., Lewis, G. F., McConnachie, A. W., et al. 2014, ApJ, 780, 128
work page 2014
- [45]
- [46]
-
[47]
Iodice, E., Spavone, M., Capaccioli, M., et al. 2019, A&A, 623, A1
work page 2019
-
[48]
2021, The Messenger, 183, 25 Ivezi´c, Ž., Kahn, S
Iodice, E., Spavone, M., Capaccioli, M., et al. 2021, The Messenger, 183, 25 Ivezi´c, Ž., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 111
work page 2021
- [49]
- [50]
-
[51]
Kalirai, J. S., Gilbert, K. M., Guhathakurta, P., et al. 2006, ApJ, 648, 389 Krajnovi´c, D., Emsellem, E., den Brok, M., et al. 2018, MNRAS, 477, 5327
work page 2006
-
[52]
Lane, J. M. M., Bovy, J., & Mackereth, J. T. 2022, MNRAS, 510, 5119
work page 2022
-
[53]
Lee, J. & Yi, S. K. 2017, ApJ, 836, 161
work page 2017
-
[54]
Merritt, A., van Dokkum, P., Abraham, R., & Zhang, J. 2016, ApJ, 830, 62
work page 2016
-
[55]
Mihos, J. C., Harding, P., Feldmeier, J. J., et al. 2017, ApJ, 834, 16
work page 2017
-
[56]
Monachesi, A., Gómez, F. A., Grand, R. J. J., et al. 2016, MNRAS, 459, L46
work page 2016
-
[57]
Monachesi, A., Gómez, F. A., Grand, R. J. J., et al. 2019, MNRAS, 485, 2589
work page 2019
-
[58]
Montenegro-Taborda, D., Rodriguez-Gomez, V ., Pillepich, A., et al. 2023, MN- RAS, 521, 800
work page 2023
- [59]
-
[60]
Mouhcine, M., Ferguson, H. C., Rich, R. M., Brown, T. M., & Smith, T. E. 2005, ApJ, 633, 810
work page 2005
-
[61]
P., Conroy, C., Bonaca, A., et al
Naidu, R. P., Conroy, C., Bonaca, A., et al. 2020, ApJ, 901, 48
work page 2020
- [62]
-
[63]
Nelson, D., Pillepich, A., Ayromlou, M., et al. 2024, A&A, 686, A157
work page 2024
-
[64]
Nelson, D., Pillepich, A., Springel, V ., et al. 2018, MNRAS, 475, 624
work page 2018
-
[65]
2019, Computational Astrophysics and Cosmology, 6, 2
Nelson, D., Springel, V ., Pillepich, A., et al. 2019, Computational Astrophysics and Cosmology, 6, 2
work page 2019
-
[66]
Oliva-Altamirano, P., Brough, S., Lidman, C., et al. 2014, MNRAS, 440, 762
work page 2014
-
[67]
Pillepich, A., Nelson, D., Hernquist, L., et al. 2018, MNRAS, 475, 648
work page 2018
-
[68]
Pinna, F., Walo-Martín, D., Grand, R. J. J., et al. 2024, A&A, 683, A236 Planck Collaboration, Aghanim, N., Akrami, Y ., et al. 2020, A&A, 641, A6
work page 2024
-
[69]
Prada, F., Klypin, A. A., Cuesta, A. J., Betancort-Rijo, J. E., & Primack, J. 2012, MNRAS, 423, 3018
work page 2012
-
[70]
Proctor, K. L., Lagos, C. d. P., Ludlow, A. D., & Robotham, A. S. G. 2024, MNRAS, 527, 2624
work page 2024
- [71]
-
[72]
2022, Frontiers in Astronomy and Space Sciences, 9, 852810
Ragusa, R., Mirabile, M., Spavone, M., et al. 2022, Frontiers in Astronomy and Space Sciences, 9, 852810
work page 2022
- [73]
- [74]
- [75]
- [76]
-
[77]
Spavone, M., Iodice, E., Capaccioli, M., et al. 2018, ApJ, 864, 149
work page 2018
- [78]
- [79]
-
[80]
Spavone, M., Iodice, E., van de Ven, G., et al. 2020, A&A, 639, A14
work page 2020
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.