Cosmological Simulations of Stellar Halos with Gaia Sausage-Enceladus Analogues: Two Sausages, One Bun?
Pith reviewed 2026-05-23 21:55 UTC · model grok-4.3
The pith
Roughly one third of Gaia Sausage-Enceladus analogues in simulations arise from two separate mergers.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Using the TNG50 simulation, the authors identify GSE-like stellar debris in 32 Milky Way analogue galaxies out of 98 examined. In one third of these cases the debris originates from two distinct merger events rather than a single accretion. The two-merger cases involve galaxies that fell in earlier, around 10.7 Gyr ago compared to 5.9 Gyr ago for single mergers, producing observable differences in abundances and ages.
What carries the argument
Classification of accretion events as single-merger or two-merger GSE analogues based on the number of progenitor galaxies contributing to the radially biased stellar debris in simulated halos.
If this is right
- Kinematic spaces alone are insufficient to distinguish single-merger GSEs from two-merger GSEs.
- Chemical abundances and star formation histories reveal more evident differences between the two scenarios.
- Single-merger GSEs are typically accreted more recently than the galaxies in two-merger GSEs.
- Systematic shifts in abundances and ages make modeling these aspects of the stellar halo important for understanding its assembly.
Where Pith is reading between the lines
- If the real Milky Way GSE matches the two-merger population, its formation would involve an earlier pair of accretions whose debris later overlapped.
- Better age and abundance measurements of halo stars could test whether multiple mergers are common enough to explain a substantial fraction of radially biased debris.
- The same identification approach could be applied to other simulation volumes to check whether the one-third fraction holds across different model assumptions.
Load-bearing premise
The specific criteria used to tag GSE-like debris in the simulated galaxies produce populations whose properties map reliably onto the real Milky Way GSE.
What would settle it
An observational measurement of the median infall time or abundance patterns of the actual GSE stars that lies outside the reported ranges and distributions for both the single-merger and two-merger simulation populations.
Figures
read the original abstract
Observations of the Milky Way's stellar halo find that it is predominantly comprised of a radially biased population of stars, dubbed the Gaia Sausage--Enceladus, or GSE. These stars are thought to be debris from dwarf galaxy accretion early in the Milky Way's history. Though typically considered to be from a single merger, it is possible that the GSE debris has multiple sources. To investigate this possibility, we use the TNG50 simulation to identify stellar accretion histories in 98 Milky Way analogues -- the largest sample for which such an identification has been performed -- and find GSE-like debris in 32, with two-merger GSEs accounting for a third of these cases. Distinguishing single-merger GSEs from two-merger GSEs is difficult in common kinematic spaces, but differences are more evident through chemical abundances and star formation histories. This is because single-merger GSEs are typically accreted more recently than the galaxies in two-merger GSEs: the median infall times (with 16th and 84th percentiles) are $5.9^{+3.3}_{-2.0}$ and $10.7^{+1.2}_{-3.7}$ Gyr ago for single- and two-merger scenarios, respectively. The systematic shifts in abundances and ages that occur as a result suggest that efforts in modeling these aspects of the stellar halo prove ever-important in understanding its assembly.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript uses the TNG50 simulation to identify GSE-like stellar debris in 98 Milky Way analogue galaxies, the largest such sample to date. GSE-like debris is found in 32 analogues, with two-merger histories accounting for approximately one third of these cases. Single-merger and two-merger GSEs are difficult to distinguish using common kinematic spaces but differ in chemical abundances and star formation histories; the reported median infall times are 5.9^{+3.3}_{-2.0} Gyr ago for single-merger and 10.7^{+1.2}_{-3.7} Gyr ago for two-merger scenarios.
Significance. The use of the largest sample of 98 MW analogues constitutes a clear strength, enabling a statistical view of GSE-like accretion that goes beyond earlier smaller samples. If the GSE identification procedure is shown to be faithful to the observed Milky Way GSE, the result that two-merger cases are common (~1/3) and produce systematic shifts in abundances and ages would be significant for halo assembly models and for guiding observational efforts to distinguish merger histories.
major comments (2)
- [GSE identification method] Section on GSE identification method: the criteria used to tag GSE-like debris are not demonstrated to reproduce the full observed GSE phase-space density, metallicity distribution function, and age spread. Without this validation, the central quantitative claim (32/98 analogues, two-merger fraction ~1/3) risks contamination by unrelated accretion events, rendering the reported infall-time offset a possible selection artifact.
- [Results on infall times] Results on infall times and fractions: the median values (5.9^{+3.3}_{-2.0} and 10.7^{+1.2}_{-3.7} Gyr) and the one-third two-merger fraction are presented without quantitative tests of robustness to variations in the tagging thresholds, simulation resolution, or baryonic physics choices.
minor comments (1)
- The abstract would be improved by a one-sentence summary of the precise kinematic/chemical cuts used to define GSE-like debris.
Simulated Author's Rebuttal
We thank the referee for their constructive comments, which highlight important aspects of our GSE identification and robustness analysis. We address each major comment below and indicate planned revisions.
read point-by-point responses
-
Referee: Section on GSE identification method: the criteria used to tag GSE-like debris are not demonstrated to reproduce the full observed GSE phase-space density, metallicity distribution function, and age spread. Without this validation, the central quantitative claim (32/98 analogues, two-merger fraction ~1/3) risks contamination by unrelated accretion events, rendering the reported infall-time offset a possible selection artifact.
Authors: We agree that the original manuscript did not include a direct, quantitative comparison of the tagged GSE-like debris to the full observed GSE phase-space density, MDF, and age spread. Our identification criteria (radial anisotropy β > 0.5 combined with a metallicity window motivated by Gaia observations) were chosen to select analogues rather than exact replicas. To address the concern, we will add a dedicated validation subsection that overlays the simulated debris distributions against published observational constraints from Gaia DR3 and spectroscopic surveys. This addition will quantify the degree of fidelity and explicitly discuss any residual risk of contamination from unrelated events. revision: yes
-
Referee: Results on infall times and fractions: the median values (5.9^{+3.3}_{-2.0} and 10.7^{+1.2}_{-3.7} Gyr) and the one-third two-merger fraction are presented without quantitative tests of robustness to variations in the tagging thresholds, simulation resolution, or baryonic physics choices.
Authors: We acknowledge that explicit robustness tests were absent from the submitted version. We will add a new analysis varying the anisotropy and metallicity thresholds by ±20 % and report the resulting changes to the two-merger fraction and median infall times; these results will be presented in a revised methods/results section. However, tests involving different simulation resolutions or alternative baryonic physics implementations are not possible within the current study, as they would require entirely separate simulation suites beyond TNG50. revision: partial
- Quantitative robustness tests to variations in simulation resolution and baryonic physics choices cannot be performed, as they require access to or computation of simulation data outside the TNG50 suite.
Circularity Check
No circularity: results extracted directly from TNG50 simulation outputs
full rationale
The paper reports counts and statistics (32/98 analogues with GSE-like debris, ~1/3 two-merger cases, median infall times 5.9 vs 10.7 Gyr) obtained by applying an identification procedure to the external TNG50 simulation volume. No equations, fitted functional forms, or ansatzes are invoked that reduce the reported fractions or time offsets to the inputs by construction. No load-bearing self-citations appear in the abstract or described chain; the simulation data and tagging criteria are independent of the present authors' prior results. This is the expected non-circular outcome for a direct simulation analysis paper.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption TNG50 baryonic physics and resolution are adequate to produce realistic GSE analogues
Forward citations
Cited by 1 Pith paper
-
Galactic Amnesia: The Information Washout of the Milky Way Merger History
Mutual information analysis of TNG50 simulations shows gravitational potential and total energy retain merger mass and infall time information longest, while radial velocity loses it within ~5 Gyr, with washout depend...
Reference graph
Works this paper leans on
-
[1]
Abdurro'uf , Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35, 10.3847/1538-4365/ac4414
-
[2]
Akins, H. B., Christensen, C. R., Brooks, A. M., et al. 2021, ApJ, 909, 139, 10.3847/1538-4357/abe2ab
-
[3]
Amarante, J. A. S., Beraldo e Silva , L., Debattista, V. P., & Smith, M. C. 2020, ApJ, 891, L30, 10.3847/2041-8213/ab78a4
-
[4]
Amarante, J. A. S., Debattista, V. P., Beraldo e Silva , L., Laporte, C. F. P., & Deg, N. 2022, ApJ, 937, 12, 10.3847/1538-4357/ac8b0d
-
[5]
Astropy Collaboration , Price-Whelan , A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, 10.3847/1538-4357/ac7c74
work page internal anchor Pith review doi:10.3847/1538-4357/ac7c74 2022
-
[6]
Belokurov, V., Erkal, D., Evans, N. W., Koposov, S. E., & Deason, A. J. 2018, MNRAS, 478, 611, 10.1093/mnras/sty982
-
[7]
Belokurov, V., Sanders, J. L., Fattahi, A., et al. 2020, MNRAS, 494, 3880, 10.1093/mnras/staa876
-
[8]
Belokurov, V., Vasiliev, E., Deason, A. J., et al. 2023, MNRAS, 518, 6200, 10.1093/mnras/stac3436
-
[9]
Benisty, D., Vasiliev, E., Evans, N. W., et al. 2022, ApJ, 928, L5, 10.3847/2041-8213/ac5c42
-
[10]
Bignone, L. A., Helmi, A., & Tissera, P. B. 2019, ApJ, 883, L5, 10.3847/2041-8213/ab3e0e
-
[11]
2016, ARA&A, 54, 529, 10.1146/annurev-astro-081915-023441
Bland-Hawthorn , J., & Gerhard, O. 2016, ARA&A, 54, 529, 10.1146/annurev-astro-081915-023441
-
[12]
Bonaca, A., Conroy, C., Cargile, P. A., et al. 2020, ApJ, 897, L18, 10.3847/2041-8213/ab9caa
-
[13]
Bose, S., Deason, A. J., Belokurov, V., & Frenk, C. S. 2020, MNRAS, 495, 743, 10.1093/mnras/staa1199
-
[14]
Boylan-Kolchin , M., Springel, V., White, S. D. M., & Jenkins, A. 2010, MNRAS, 406, 896, 10.1111/j.1365-2966.2010.16774.x
-
[15]
Buch, D., Nadler, E. O., Wechsler, R. H., & Mao, Y.-Y. 2024, preprint, 10.48550/arXiv.2404.08043
-
[16]
Buder, S., Sharma, S., Kos, J., et al. 2021, MNRAS, 506, 150, 10.1093/mnras/stab1242
-
[17]
Bullock, J. S., & Johnston, K. V. 2005, ApJ, 635, 931, 10.1086/497422
-
[18]
Busha, M. T., Wechsler, R. H., Behroozi, P. S., et al. 2011, ApJ, 743, 117, 10.1088/0004-637X/743/2/117
-
[19]
Cautun, M., Ben \'i tez-Llambay , A., Deason, A. J., et al. 2020, MNRAS, 494, 4291, 10.1093/mnras/staa1017
-
[20]
J., Basu, S., Huber, D., et al
Chaplin, W. J., Basu, S., Huber, D., et al. 2014, ApJS, 210, 1, 2020022814070166
work page 2014
- [21]
-
[22]
2019, ApJ, 883, 107, doi: 10.3847/1538-4357/ab38b8
Conroy, C., Bonaca, A., Cargile, P., et al. 2019, ApJ, 883, 107, 10.3847/1538-4357/ab38b8
-
[23]
Cooper, A. P., Koposov, S. E., Allende Prieto, C., et al. 2023, ApJ, 947, 37, 10.3847/1538-4357/acb3c0
-
[24]
Research in Astronomy and Astrophysics , year = 2012, month = sep, volume =
Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197, 10.1088/1674-4527/12/9/003
-
[25]
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. M. 1985, ApJ, 292, 371, 10.1086/163168
-
[26]
Deason, A. J., & Belokurov, V. 2024, preprint, 10.48550/arXiv.2402.12443
-
[27]
Deason, A. J., Belokurov, V., Koposov, S. E., & Lancaster, L. 2018, ApJ, 862, L1, 10.3847/2041-8213/aad0ee
-
[28]
Deason, A. J., Belokurov, V., & Weisz, D. R. 2015, MNRAS, 448, L77, 10.1093/mnrasl/slv001
-
[29]
Deason, A. J., Mao, Y.-Y., & Wechsler, R. H. 2016, ApJ, 821, 5, 10.3847/0004-637X/821/1/5
-
[30]
Di Cintio, A., Mostoghiu, R., Knebe, A., & Navarro, J. F. 2021, MNRAS, 506, 531, 10.1093/mnras/stab1682
-
[31]
Dillamore, A. M., Belokurov, V., Font, A. S., & McCarthy, I. G. 2022, MNRAS, 513, 1867, 10.1093/mnras/stac1038
-
[32]
Dolag, K., Borgani, S., Murante, G., & Springel, V. 2009, MNRAS, 399, 497, 10.1111/j.1365-2966.2009.15034.x
-
[33]
Donlon, T., Newberg, H. J., Sanderson, R., et al. 2024, MNRAS, 531, 1422, 10.1093/mnras/stae1264
-
[34]
Donlon, II, T., Newberg, H. J., Kim, B., & L \'e pine, S. 2022, ApJ, 932, L16, 10.3847/2041-8213/ac7531
-
[35]
J., Sanderson, R., & Widrow, L
Donlon, II, T., Newberg, H. J., Sanderson, R., & Widrow, L. M. 2020, ApJ, 902, 119, 10.3847/1538-4357/abb5f6
-
[36]
D'Souza, R., & Bell, E. F. 2018, MNRAS, 474, 5300, 10.1093/mnras/stx3081
-
[37]
Dancing in the dark: galactic properties trace spin swings along the cosmic web
Dubois, Y., Pichon, C., Welker, C., et al. 2014, MNRAS, 444, 1453, 10.1093/mnras/stu1227
work page internal anchor Pith review doi:10.1093/mnras/stu1227 2014
-
[38]
Elias, L. M., Sales, L. V., Helmi, A., & Hernquist, L. 2020, MNRAS, 495, 29, 10.1093/mnras/staa1090
-
[39]
Engler, C., Pillepich, A., Joshi, G. D., et al. 2023, MNRAS, 522, 5946, 10.1093/mnras/stad1357
-
[40]
Ernandes, H., Feuillet, D., Feltzing, S., & Sk \'u lad \'o ttir, \'A . 2024, preprint
work page 2024
-
[41]
Evans, T. A., Fattahi, A., Deason, A. J., & Frenk, C. S. 2020, MNRAS, 497, 4311, 10.1093/mnras/staa2202
-
[42]
Fattahi, A., Belokurov, V., Deason, A. J., et al. 2019, MNRAS, 484, 4471, 10.1093/mnras/stz159
-
[43]
Fattahi, A., Deason, A. J., Frenk, C. S., et al. 2020, MNRAS, 497, 4459, 10.1093/mnras/staa2221
-
[44]
2023, MNRAS, 522, 3831, 10.1093/mnras/stad1205
Feldmann, R., Quataert, E., Faucher-Gigu \`e re , C.-A., et al. 2023, MNRAS, 522, 3831, 10.1093/mnras/stad1205
-
[45]
Feuillet, D. K., Sahlholdt, C. L., Feltzing, S., & Casagrande, L. 2021, MNRAS, 508, 1489, 10.1093/mnras/stab2614
-
[46]
Flynn, C., Holmberg, J., Portinari, L., Fuchs, B., & Jahrei , H. 2006, MNRAS, 372, 1149, 10.1111/j.1365-2966.2006.10911.x
-
[47]
Font, A. S., Johnston, K. V., Bullock, J. S., & Robertson, B. E. 2006, ApJ, 646, 886, 10.1086/505131
-
[48]
Font, A. S., McCarthy, I. G., Crain, R. A., et al. 2011, MNRAS, 416, 2802, 10.1111/j.1365-2966.2011.19227.x
-
[49]
Gaia Collaboration , Vallenari, A., Brown, A. G. A., et al. 2022, Gaia Data Release 3: Summary of the Content and Survey Properties, 10.48550/arXiv.2208.00211
-
[50]
Gallart, C., Bernard, E. J., Brook, C. B., et al. 2019, Nature Astronomy, 3, 932, 10.1038/s41550-019-0829-5
-
[51]
Garavito-Camargo , N., Besla, G., Laporte, C. F. P., et al. 2021, ApJ, 919, 109, 10.3847/1538-4357/ac0b44
-
[52]
Garc \'i a-Bethencourt , G., Brook, C. B., Grand, R. J. J., & Kawata, D. 2023, MNRAS, 526, 1190, 10.1093/mnras/stad2832
-
[53]
G \'o mez, F. A., & Helmi, A. 2010, MNRAS, 401, 2285, 10.1111/j.1365-2966.2009.15841.x
-
[54]
2016, Pydata/Bottleneck, Python for Data
Goodman, K. 2016, Pydata/Bottleneck, Python for Data. https://github.com/pydata/bottleneck
work page 2016
-
[55]
Grand, R. J. J., Kawata, D., Belokurov, V., et al. 2020, MNRAS, 497, 1603, 10.1093/mnras/staa2057
-
[56]
Grunblatt, S. K., Zinn, J. C., Price-Whelan , A. M., et al. 2021, ApJ, 916, 88, 10.3847/1538-4357/ac0532
-
[57]
J., Conroy, C., & Hernquist, L
Han, J. J., Conroy, C., & Hernquist, L. 2023 a , Nature Astronomy, 7, 1481, 10.1038/s41550-023-02076-9
-
[58]
J., Semenov, V., Conroy, C., & Hernquist, L
Han, J. J., Semenov, V., Conroy, C., & Hernquist, L. 2023 b , ApJ, 957, L24, 10.3847/2041-8213/ad0641
-
[59]
Han, J. J., Naidu, R. P., Conroy, C., et al. 2022, ApJ, 934, 14, 10.3847/1538-4357/ac795f
-
[60]
Harris, C. R., Millman, K. J., van der Walt , S. J., et al. 2020, Nature, 585, 357, 10.1038/s41586-020-2649-2
-
[61]
Hasselquist, S., Hayes, C. R., Lian, J., et al. 2021, ApJ, 923, 172, 10.3847/1538-4357/ac25f9
-
[62]
Haywood, M., Di Matteo, P., Lehnert, M. D., et al. 2018, ApJ, 863, 113, 10.3847/1538-4357/aad235
-
[63]
2020, ARA&A, 58, 205, doi: 10.1146/annurev-astro-032620-021917
Helmi, A. 2020, ARA&A, 58, 205, 10.1146/annurev-astro-032620-021917
-
[64]
Helmi, A., Babusiaux, C., Koppelman, H. H., et al. 2018, Nature, 563, 85, 10.1038/s41586-018-0625-x
-
[65]
2019, The Messenger, 175, 23, 10.18727/0722-6691/5120
Helmi, A., Irwin, M., Deason, A., et al. 2019, The Messenger, 175, 23, 10.18727/0722-6691/5120
-
[66]
2024, preprint, 10.48550/arXiv.2403.09777
Horta, D., Yuxi , Lu , et al. 2024, preprint, 10.48550/arXiv.2403.09777
-
[67]
Horta, D., Schiavon, R. P., Mackereth, J. T., et al. 2021, MNRAS, 500, 1385, 10.1093/mnras/staa2987
-
[68]
---. 2023, MNRAS, 520, 5671, 10.1093/mnras/stac3179
-
[69]
Hunter, J. D. 2007, Comput. Sci. Eng., 9, 90, 10.1109/MCSE.2007.55
-
[70]
2024, ApJ, 967, 89, 10.3847/1538-4357/ad382d
Ibata, R., Malhan, K., Tenachi, W., et al. 2024, ApJ, 967, 89, 10.3847/1538-4357/ad382d
-
[71]
2021, MNRAS, 502, 5686, 10.1093/mnras/stab005
Iorio, G., & Belokurov, V. 2021, MNRAS, 502, 5686, 10.1093/mnras/stab005
-
[72]
LSST: From Science Drivers to Reference Design and Anticipated Data Products , journal =
Ivezi \'c , Z ., Kahn, S. M., Tyson, J. A., et al. 2019, ApJ, 873, 111, 10.3847/1538-4357/ab042c
-
[73]
2017, A&A, 604, A106, 10.1051/0004-6361/201629691
Jean-Baptiste , I., Di Matteo, P., Haywood, M., et al. 2017, A&A, 604, A106, 10.1051/0004-6361/201629691
-
[74]
Jin, S., Trager, S. C., Dalton, G. B., et al. 2024, MNRAS, 530, 2688, 10.1093/mnras/stad557
-
[75]
Johnson, J. W., Conroy, C., Johnson, B. D., et al. 2023, MNRAS, 526, 5084, 10.1093/mnras/stad2985
-
[76]
Johnston, K. V. 1998, ApJ, 495, 297, 10.1086/305273
-
[77]
Johnston, K. V., Hernquist, L., & Bolte, M. 1996, ApJ, 465, 278, 10.1086/177418
-
[78]
Kafle, P. R., Sharma, S., Lewis, G. F., Robotham, A. S. G., & Driver, S. P. 2018, MNRAS, 475, 4043, 10.1093/mnras/sty082
-
[79]
Karachentsev, I. D. 2005, AJ, 129, 178, 10.1086/426368
-
[80]
2016, in Positioning and Power in Academic Publishing: Players , Agents and Agendas, ed
Kluyver, T., Ragan-Kelley , B., P \'e rez, F., et al. 2016, in Positioning and Power in Academic Publishing: Players , Agents and Agendas, ed. F. Loizides & B. Scmidt (IOS Press), 87--90. https://ui.adsabs.harvard.edu/abs/2016ppap.book...87K
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.