Pith. sign in

REVIEW 3 major objections 4 minor 40 references

The Epoch of the GSE Merger: Insights from the Splash and Thick Disk Age Distributions

T0 review · 3 major / 4 minor · reviewed 2026-07-12 · grok-4.5

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

Pith's one-line read Age uncertainties make Splash stars look older than the thick disk; that shift dates the GSE merger to 10.1 Gyr ago.

desk verdict Clean Eddington-bias explanation of the Splash age paradox that yields a useful ~10 Gyr GSE clock; the advertised ±0.2 Gyr is model-dependent and selection already moves the answer by ~0.7 Gyr. read the letter →

arxiv 2607.02940 v1 pith:6YLLKOVD submitted 2026-07-03 astro-ph.GA

Sihan Lu , Yang Huang , Qikang Feng , Huawei Zhang This is my paper
classification astro-ph.GA
keywords Gaia-Sausage-EnceladusSplashthickdiskstellaragesMilkyWaymergerEddingtonbiasstarformationhistory
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

The last major merger the Milky Way experienced, Gaia-Sausage-Enceladus, heated part of the early thick disk and threw stars onto low-angular-momentum orbits that we call the Splash. Those Splash stars appear systematically older than the thick-disk population they are thought to come from, which has been a long-standing puzzle. This paper shows the offset is an observational bias: the merger cuts off Splash formation at a sharp epoch, and age measurement errors then shift the truncated distribution toward older ages (an Eddington-like bias). Because the size of that shift depends on when the merger happened relative to the peak of thick-disk star formation, the observed age offset becomes a clock. Combining the measured peak offset with a model of the thick-disk age distribution yields a merger epoch of 10.1 Gyr ago with a half-Gyr-wide uncertainty, one of the tightest timings available. A refined chemical-plus-kinematic selection and a constant Splash-to-disk number ratio at the oldest ages independently support the same window.

What carries the argument

Eddington-like peak shift of a truncated age distribution: the merger imposes a hard upper cutoff on Splash ages; Gaussian age errors then move the observed peak older by an amount that is a monotonic function of merger epoch relative to the thick-disk star-formation peak.

What would settle it

A large sample of Splash and thick-disk ages with substantially smaller, well-calibrated uncertainties (for example from asteroseismology) whose measured peak offset no longer matches the predicted Delta-t versus t_merger curve under the same selection.

Watch

Extended reading notes

Core claim

The apparent older age of Splash stars relative to the thick disk is produced by truncation at the GSE merger plus age uncertainties; the magnitude of the resulting peak offset maps directly onto merger time, giving t_GSE = 10.1^{+0.2}_{-0.2} Gyr ago.

Load-bearing premise

That stars of every age in the pre-merger thick disk were equally likely to be heated into the Splash, and that the thick-disk star-formation history can be treated as roughly Gaussian when converting the observed peak offset into a merger time.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. The paper argues that the observed older age distribution of the Splash relative to the high-α thick disk is an Eddington-like bias: GSE truncates the Splash age distribution at the merger epoch, and convolution with age errors then shifts the truncated peak older. Using the XR22 subgiant catalog, chemically selected high-α stars, and age-dependent kinematic cuts in Vφ–[Fe/H], the authors measure a peak offset Δt = 0.30^{+0.06}_{-0.06} Gyr relative to a thick-disk SFH peak at 11.58 Gyr. Forward-modeling a Gaussian thick-disk SFH truncated at t_merger and convolved with the catalog age errors maps this offset to t_GSE = 10.1^{+0.2}_{-0.2} Gyr. An independent Splash-to-thick-disk number-ratio diagnostic is roughly constant at the oldest ages under the refined selection, supporting a merger near 10–11 Gyr.

Significance. If the mapping from observed peak offset to merger epoch is robust, the result supplies one of the tightest empirical anchors for the last major merger of the Milky Way and cleanly resolves a long-standing age paradox between Splash and thick-disk stars. The work is carefully executed on a high-quality public catalog, combines two complementary diagnostics (number ratio and peak offset), and makes the bias mechanism transparent with a simple forward model. Even if the quoted 0.2 Gyr precision is optimistic, the qualitative explanation of the age offset and the ~10–11 Gyr window remain useful contributions to Galactic archaeology.

major comments (3)
  1. §3.3 and Fig. 5: the quoted posterior width ±0.2 Gyr is generated under an explicit Gaussian SFH assumption (Fig. 4 and the quadrature-width construction). The manuscript itself notes that a more realistic SFH may refine the epoch, yet no alternative shapes (e.g., skewed or exponentially declining) are propagated into the Δt–t_merger map. Because the functional form of f(t_merger | t_peak) depends on the SFH shape, the statistical error alone understates the systematic uncertainty on t_GSE.
  2. Appendix B vs. main text: switching from the refined to the classical Vφ–[Fe/H] selection changes Δt from 0.30 to 0.70 Gyr and shifts the inferred merger epoch from 10.1 to 10.83 Gyr—already ~3.5 times the quoted statistical error. The refined boundaries were chosen so that the Splash-to-disk ratio becomes flat at old ages (§3.1, right panel of Fig. 2), which is precisely the age-independent heating assumption used in the forward model. Selection choice and the heating assumption are therefore not independent; the final uncertainty budget should include this systematic or justify why the refined selection is uniquely preferred beyond the ratio plateau.
  3. §3.1: the age-independent heating probability is stated as an assumption and is only indirectly supported by the flat number ratio under the refined cuts. If heating efficiency varied with stellar age or orbital properties (e.g., older, hotter stars more easily scattered), the intrinsic Splash distribution would already differ from the thick-disk SFH before truncation, altering the predicted peak shift. A brief test or literature-based bound on age-dependent heating efficiency would strengthen the central claim.
minor comments (4)
  1. Fig. 3 and §3.2: the KDE peaks are reported to high precision; a short statement of the KDE bandwidth choice and its effect on the measured Δt would improve reproducibility.
  2. §2.2: the age-dependent selection boundaries (green lines in Fig. 2) are described qualitatively; an explicit functional form or tabulated thresholds per age bin would allow others to reproduce the refined sample.
  3. Appendix C: the GSE dynamical box is given, but the resulting sample size and any residual contamination estimate would help the reader judge the comparison of peak ages.
  4. Throughout: a few typographical inconsistencies (e.g., spacing around ± and superscripts in the abstract and §3.3) should be cleaned for the final version.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: t_GSE is obtained by inverting a forward model of truncation plus age errors against the observed peak offset, not by definitional identity or fitted-input renaming.

full rationale

The central result (t_GSE = 10.1^{+0.2}_{-0.2} Gyr) is produced by measuring an empirical peak offset Δt between the Splash and thick-disk age distributions, then finding the merger epoch whose truncated-and-convolved model reproduces that offset (Eq. 1 and Fig. 5). The mapping Δt ~ f(t_merger | t_peak_thick) is generated from an explicit physical model (Gaussian SFH truncated at the merger, then convolved with the catalog age uncertainties); it is not an algebraic identity that forces t_GSE to equal any input quantity. The refined kinematic selection is motivated by the same age-independent-heating assumption used in the model, and the classical selection yields a shifted but still consistent value (Appendix B); this is a robustness issue, not a circular reduction. No self-citation supplies a uniqueness theorem or load-bearing ansatz that closes the loop, and the number-ratio diagnostic is independent of the peak-offset fit. The Gaussian SFH and constant-heating assumptions are under-tested and affect the quoted precision, but those are model-dependence concerns outside the circularity criteria. The derivation is therefore self-contained parameter estimation against external data (XR22 ages), warranting only a minimal score for the mild interdependence of selection and modeling assumption.

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

The result rests on standard Galactic-archaeology assumptions plus two modeling choices (Gaussian SFH, age-independent heating) and data-driven free parameters (SFH peak/width, selection boundaries, the fitted t_GSE itself). No new physical entities are invented; the ledger is dominated by domain assumptions and fitted descriptors of the observed age distributions.

free parameters (5)
  • thick-disk SFH Gaussian peak = 11.58^{+0.02}_{-0.02} Gyr
    Measured from KDE/Monte-Carlo of the high-α sample (11.58 Gyr); used as the fixed prior when mapping Δt to t_merger.
  • thick-disk SFH intrinsic width = data-derived (observed width minus median σ_age)
    Estimated by subtracting median age uncertainty in quadrature from the observed distribution; controls how much truncation shifts the peak.
  • observed peak offset Δt = 0.30^{+0.06}_{-0.06} Gyr
    Monte-Carlo peak difference between Splash and thick-disk KDEs; the observable that is matched by the model.
  • age-dependent Splash selection boundaries = age-binned Vφ–[Fe/H] rectangles (visual)
    Green lines in Fig. 2 chosen so that Splash/thick number ratio is flat at old ages; changes sample purity and therefore measured Δt.
  • t_GSE (merger epoch) = 10.1^{+0.2}_{-0.2} Gyr
    Target parameter obtained by intersecting the model Δt(t_merger) curve with the observed Δt band (Fig. 5).
assumptions (4)
  • domain assumption Pre-merger thick-disk stars of different ages were heated into the Splash with approximately equal probability, so their intrinsic age distributions were identical before truncation.
    Stated in §3.1; required for both the flat number-ratio expectation and the shared pre-merger SFH used in the peak-offset model.
  • ad hoc to paper Thick-disk star-formation history is adequately described by a single Gaussian for the purpose of predicting the peak shift.
    Adopted in §3.3 and Fig. 4 for simplicity; authors explicitly leave more realistic SFHs to future work.
  • domain assumption The GSE merger can be treated as an instantaneous event that fully truncates further Splash formation.
    §3.2–3.3; produces the sharp cutoff whose convolution generates the Eddington bias.
  • domain assumption Stellar age uncertainties in the XR22 subgiant catalog are well-characterized (approximately Gaussian) and correctly reported.
    Used for all convolutions and Monte-Carlo draws; supported by the Shariat et al. (2025) wide-binary check cited in §4.

how reviews work

0 comments
Cite this review

Pith. "Pith review of The Epoch of the GSE Merger: Insights from the Splash and Thick Disk Age Distributions." pith.science (2026). https://pith.science/paper/6YLLKOVD

@misc{pith2026260702940,
  author       = {Pith},
  title        = {Pith review of: The Epoch of the GSE Merger: Insights from the Splash and Thick Disk Age Distributions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/6YLLKOVD}},
  note         = {Machine review of arXiv:2607.02940}
}
abstract

The epoch of the Gaia-Sausage-Enceladus (GSE) merger, the last major merger experienced by the Milky Way, is crucial for reconstructing the Galaxy's evolutionary history. This event dynamically heated the disk, scattering some stars into the halo and producing the so-called Splash. Yet the observed age distribution of the Splash is systematically older than that of the thick disk from which it is thought to originate, posing a puzzling inconsistency. In this work, we show that this apparent discrepancy can be naturally explained once stellar age uncertainties are taken into account. The GSE truncated the intrinsic age distribution of the Splash at the merger epoch, while measurement errors introduce an Eddington-like bias that systematically shifts the truncated distribution toward older ages. Importantly, the magnitude of this shift depends on both the thick-disk star formation history and the merger epoch, thereby offering a new way to constrain the timing of the merger. By combining the observed Splash age distribution with the peak of thick-disk star formation, we infer a GSE merger epoch of $10.1^{+0.2}_{-0.2}$ Gyr ago, providing one of the tightest constraints to date. This result offers new insight into the Milky Way's accretion history and opens a path toward more robust reconstructions of its early evolution.

Figures

Figures reproduced from arXiv: 2607.02940 by the authors.

Figure 1
Figure 1. [α/Fe]–[Fe/H] diagram for the adopted subgiant sample, showing a significant bimodality, with the low-α thin-disk and high-α thick-disk sequences. High-α sequence stars were then selected using [α/Fe] > −0.12×[Fe/H] + 0.1, [α/Fe] > 0, and [Fe/H] > −1.5 (red lines). Stars below the red lines are dominated by the traditional thin-disk popula￾tion (Imig et al. 2023). from the effective temperature–absolute magnitude di… view at source ↗
Figure 2
Figure 2. Left: Row-normalized rotation-velocity–[Fe/H] diagrams comparing the classical selection (red) with the refined criteria (green). The first panel shows the population boundaries from Belokurov et al. (2020b), and the second panel presents the distribution of all high-α subsample. The remaining panels show the distributions in different stellar age bins. In the oldest bin (bottom-right), the refined criteria adopt a … view at source ↗
Figure 3
Figure 3. Age distribution histograms and their corre￾sponding KDEs for the thick disk (orange) and Splash (blue) populations. The majority of stars in both components have ages older than 10 Gyr, indicating that the present-day thick disk formed at early times. Nevertheless, the Splash popu￾lation shows a relatively higher fraction at the oldest ages compared to the thick disk. previous estimates of the merger epoch and supp… view at source ↗
Figures from the paper (1 more)
Figure 4
Figure 4. Figure 4: Illustration of why the Splash appears older than the thick disk. We assume that the thick-disk star-formation history (orange) follows a Gaussian distribution. Prior to the GSE merger, the intrinsic age distribution of Splash stars (blue) is identical to that of the t…

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

40 extracted references · 3 linked inside Pith

  1. [1]

    2022, ApJS, 259, 35

    Abdurro’uf, Accetta, K., Aerts, C., et al. 2022, ApJS, 259, 35

    2022

  2. [2]

    Amarante, J. A. S., Beraldo e Silva, L., Debattista, V. P., & Smith, M. C. 2020, ApJL, 891, L30

    2020

  3. [3]

    Amarante, J. A. S., Laporte, C. F. P., Debattista, V. P., et al. 2025, arXiv e-prints, arXiv:2512.01293

    arXiv 2025

  4. [4]

    S., Hirai, Y., & Beers, T

    An, D., Lee, Y. S., Hirai, Y., & Beers, T. C. 2025, A&A, 700, A273

    2025

  5. [5]

    Deason, A. J. 2018, MNRAS, 478, 611

    2018

  6. [6]

    2012, MNRAS, 427, 127

    Bressan, A., Marigo, P., Girardi, L., et al. 2012, MNRAS, 427, 127

    2012

  7. [7]

    2026, MNRAS, 548, stag729

    Buder, S., Buck, T., Sk´ ulad´ ottir,´A., et al. 2026, MNRAS, 548, stag729

    2026

  8. [8]

    2018, Monthly Notices of the Royal Astronomical Society, 478, 4513

    Buder, S., Asplund, M., Duong, L., et al. 2018, Monthly Notices of the Royal Astronomical Society, 478, 4513

    2018

Show all 40 references
  1. [9]

    K., et al

    Buder, S., Lind, K., Ness, M. K., et al. 2022, MNRAS, 510, 2407

    2022

  2. [10]

    J., Fattahi, A., Callingham, T

    Carrillo, A., Deason, A. J., Fattahi, A., Callingham, T. M., & Grand, R. J. J. 2024, MNRAS, 527, 2165

    2024

  3. [11]

    J., Serenelli, A

    Chaplin, W. J., Serenelli, A. M., Miglio, A., et al. 2020, Nature Astronomy, 4, 382

    2020

  4. [12]

    H., Naidu, R

    Conroy, C., Weinberg, D. H., Naidu, R. P., et al. 2022, arXiv e-prints, arXiv:2204.02989

    2022 arXiv

  5. [13]

    2012, Research in Astronomy and Astrophysics, 12, 1197

    Cui, X.-Q., Zhao, Y.-H., Chu, Y.-Q., et al. 2012, Research in Astronomy and Astrophysics, 12, 1197

    2012

  6. [14]

    2020, MNRAS, 493, 5195

    Das, P., Hawkins, K., & Jofr´ e, P. 2020, MNRAS, 493, 5195

    2020

  7. [15]

    J., & Belokurov, V

    Deason, A. J., & Belokurov, V. 2024, NewAR, 99, 101706

    2024

  8. [16]

    J., Belokurov, V., Evans, N

    Deason, A. J., Belokurov, V., Evans, N. W., & Johnston, K. V. 2013, The Astrophysical Journal, 763, 113

    2013

  9. [17]

    Demarque, P., Woo, J.-H., Kim, Y.-C., & Yi, S. K. 2004, ApJS, 155, 667 Di Matteo, P., Haywood, M., Lehnert, M. D., et al. 2019, A&A, 632, A4

    2004

  10. [18]

    M., Belokurov, V., Evans, N

    Dillamore, A. M., Belokurov, V., Evans, N. W., & Davies, E. Y. 2023, MNRAS, 524, 3596

    2023

  11. [19]

    Eddington, A. S. 1913, MNRAS, 73, 359

    1913

  12. [20]

    2021, Monthly Notices of the Royal Astronomical Society, 508, 1489 8 Gaia Collaboration, Prusti, T., de Bruijne, J

    Casagrande, L. 2021, Monthly Notices of the Royal Astronomical Society, 508, 1489 8 Gaia Collaboration, Prusti, T., de Bruijne, J. H. J., et al. 2016, A&A, 595, A1 Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, A&A, 616, A1 —. 2021, A&A, 649, A1 Gaia Collabor...

    2021

  13. [21]

    J., Brook, C

    Gallart, C., Bernard, E. J., Brook, C. B., et al. 2019, Nature Astronomy, 3, 932

    2019

  14. [22]

    2015, MNRAS, 453, 758

    Hawkins, K., Jofr´ e, P., Masseron, T., & Gilmore, G. 2015, MNRAS, 453, 758

    2015

  15. [23]

    D., et al

    Haywood, M., Di Matteo, P., Lehnert, M. D., et al. 2018, ApJ, 863, 113

    2018

  16. [24]

    H., et al

    Helmi, A., Babusiaux, C., Koppelman, H. H., et al. 2018, Nature, 563, 85

    2018

  17. [25]

    Helmi, A., & White, S. D. M. 1999, Monthly Notices of the Royal Astronomical Society, 307, 495

    1999

  18. [26]

    2020, ApJS, 249, 29

    Huang, Y., Sch¨ onrich, R., Zhang, H., et al. 2020, ApJS, 249, 29

    2020

  19. [27]

    A., et al

    Imig, J., Price, C., Holtzman, J. A., et al. 2023, The Astrophysical Journal, 954, 124

    2023

  20. [28]

    S., Beers, T

    Lee, Y. S., Beers, T. C., An, D., et al. 2011, The Astrophysical Journal, 738, 187

    2011

  21. [29]

    L., Zhao, Y.-H., Zhao, G., et al

    Luo, A. L., Zhao, Y.-H., Zhao, G., et al. 2015, Research in Astronomy and Astrophysics, 15, 1095

    2015

  22. [30]

    R., Schiavon, R

    Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, The Astronomical Journal, 154, 94

    2017

  23. [31]

    A., et al

    Malhan, K., Yuan, Z., Ibata, R. A., et al. 2021, The Astrophysical Journal, 920, 51

    2021

  24. [32]

    S., Pilkington, K., Gibson, B

    Miranda, M. S., Pilkington, K., Gibson, B. K., et al. 2016, A&A, 587, A10 Montalb´ an, J., Mackereth, J. T., Miglio, A., et al. 2021, Nature Astronomy, 5, 640

    2016

  25. [33]

    2019, MNRAS, 488, 1235

    Belokurov, V. 2019, MNRAS, 488, 1235

    2019

  26. [34]

    P., Conroy, C., Bonaca, A., et al

    Naidu, R. P., Conroy, C., Bonaca, A., et al. 2020, The Astrophysical Journal, 901, 48 —. 2021, The Astrophysical Journal, 923, 92

    2020

  27. [35]

    2025, arXiv e-prints, arXiv:2510.08675

    Shariat, C., El-Badry, K., & Bhattacharjee, S. 2025, arXiv e-prints, arXiv:2510.08675

    2025 arXiv

  28. [36]

    2023, A&A, 675, A26

    Wang, C., Huang, Y., Zhou, Y., & Zhang, H. 2023, A&A, 675, A26

    2023

  29. [37]

    2022, Nature, 603, 599

    Xiang, M., & Rix, H.-W. 2022, Nature, 603, 599

    2022

  30. [38]

    2019, ApJS, 245, 34

    Xiang, M., Ting, Y.-S., Rix, H.-W., et al. 2019, ApJS, 245, 34

    2019

  31. [39]

    2001, The Astrophysical Journal Supplement Series, 136, 417

    Yi, S., Demarque, P., Kim, Y.-C., et al. 2001, The Astrophysical Journal Supplement Series, 136, 417

    2001

  32. [40]

    Zhao, G., Zhao, Y.-H., Chu, Y.-Q., Jing, Y.-P., & Deng, L.-C. 2012, Research in Astronomy and Astrophysics, 12, 723 9 APPENDIX A.IMPACT OF ERROR ON NUMBER RATIO Figure A1.Age dependence of the Splash-to-thick-disk number ratio under our merger-ejection toy model. The green das...

    2012

Pith tools

Reviewed July 12, 2026 · model on record in the stance chip above.