REVIEW 2 major objections 1 minor 64 references
Mild radial gas flows of 1.5 km/s let a chemical evolution model match the Milky Way disc's [O/Fe]-[Fe/H] distribution, stellar density profile, and abundance gradients at once.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.3
2026-07-01 05:05 UTC pith:2U7UJXKH
load-bearing objection The Lagrangian method of characteristics is a clean technical step for radial flows, but the simultaneous fits rest on the two-infall accretion history. the 2 major comments →
Chemical evolution of the Milky Way disc with radial gas flows: a Lagrangian approach
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
When radial gas flows are present, the chemical abundances of the gas at a given radius result from its whole inward journey in the disc, reflecting the star formation and accretion experienced at every radius it crossed; models with mild flows of v = 1.5 km/s reproduce simultaneously the observed [O/Fe]-[Fe/H] distribution across the disc, the present-day stellar surface-density profile, and the [Fe/H] and [O/H] gradients.
What carries the argument
The Lagrangian method of characteristics that reduces the gas surface density and abundance equations to one-dimensional integrals while accounting for radial flows and enrichment from core-collapse and Type Ia supernovae.
Load-bearing premise
The two-infall scenario for gas accretion history is assumed to hold for the Milky Way disc.
What would settle it
A direct measurement of the radial dependence of the stellar mass formed per Type Ia supernova, or of the [O/Fe] offset at fixed [Fe/H] that contradicts the 50 percent departure predicted along trajectories at 1.5 km/s, would falsify the claim.
If this is right
- The integrated stellar mass along each characteristic drops by up to an order of magnitude at v = 1.5 km/s compared with the local value.
- The [O/Fe] ratio departs from its in-situ value by up to 50 percent, making alpha-enhancement the strongest signature of the flows.
- Even mild inflows require that the full gas trajectory be tracked to recover the correct enrichment history.
- The same flow speed that fits abundances also improves the match to the observed age-abundance relations.
Where Pith is reading between the lines
- The same characteristic-tracking method could be applied to external disc galaxies that show flat or inverted abundance gradients.
- Stronger flows would demand compensatory changes in the accretion timeline or star-formation efficiency to keep the same observables in agreement.
- Comparing the model's predicted radial variation in [O/Fe] at fixed [Fe/H] against large spectroscopic surveys would test whether the 50 percent offset is observed.
- The public code allows direct substitution of different supernova yield tables to check how sensitive the flow signature remains.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript develops a semi-analytic Lagrangian model for Milky Way disc chemical evolution that incorporates radial gas flows via the method of characteristics, reducing the PDEs for gas surface density and abundances of alpha-elements and iron to one-dimensional integrals. Assuming a two-infall gas accretion history, the authors report that a constant radial flow speed of v = 1.5 km/s simultaneously reproduces the observed [O/Fe]-[Fe/H] distribution across the disc, the present-day stellar surface-density profile, the [Fe/H] and [O/H] gradients, and improves agreement with age-abundance relations. The integrated stellar mass along each characteristic is reduced relative to the local value, and the stellar mass formed per Type Ia supernova (which sets [O/Fe]) departs from its in-situ value by up to ~50%. The code is made publicly available.
Significance. If the two-infall assumption holds and the simultaneous reproduction is quantitatively verified, the work demonstrates that even mild radial flows leave a strong signature on alpha-enhancement by integrating enrichment along gas trajectories rather than at fixed radius. The reduction to one-dimensional integrals and the public code are clear strengths that support reproducibility. This approach could be adopted in other galactic chemical evolution studies to efficiently include radial mixing effects.
major comments (2)
- [Abstract] Abstract: The headline claim that v = 1.5 km/s simultaneously reproduces the [O/Fe]-[Fe/H] distribution, stellar surface-density profile, [Fe/H] and [O/H] gradients, and age-abundance relations is load-bearing on the two-infall accretion history, which is adopted without any alternative histories or sensitivity tests. Because abundances are obtained from integrals along characteristics whose source terms are set by the accretion rate at every radius crossed, a qualitatively different infall timeline would change those integrals and could invalidate the simultaneous reproduction.
- [Abstract] Abstract: No quantitative fit metrics, residual statistics, or comparison tables are presented for any of the claimed reproductions (e.g., no R² values, chi-squared, or explicit data-model overlays for the gradients or [O/Fe]-[Fe/H] plane). Without these, the strength of the evidence for the central claim cannot be assessed beyond the qualitative statement in the abstract.
minor comments (1)
- The manuscript would benefit from a dedicated section or table listing all model parameters, their adopted values, and which are fixed versus tuned.
Simulated Author's Rebuttal
We thank the referee for the constructive comments. We address each major comment below.
read point-by-point responses
-
Referee: [Abstract] Abstract: The headline claim that v = 1.5 km/s simultaneously reproduces the [O/Fe]-[Fe/H] distribution, stellar surface-density profile, [Fe/H] and [O/H] gradients, and age-abundance relations is load-bearing on the two-infall accretion history, which is adopted without any alternative histories or sensitivity tests. Because abundances are obtained from integrals along characteristics whose source terms are set by the accretion rate at every radius crossed, a qualitatively different infall timeline would change those integrals and could invalidate the simultaneous reproduction.
Authors: We agree that the results depend on the adopted two-infall accretion history. While this framework is standard in Milky Way chemical evolution studies, we will add a sensitivity analysis to variations in the infall parameters in the revised manuscript to illustrate how different timelines affect the integrals along characteristics. revision: yes
-
Referee: [Abstract] Abstract: No quantitative fit metrics, residual statistics, or comparison tables are presented for any of the claimed reproductions (e.g., no R² values, chi-squared, or explicit data-model overlays for the gradients or [O/Fe]-[Fe/H] plane). Without these, the strength of the evidence for the central claim cannot be assessed beyond the qualitative statement in the abstract.
Authors: We agree that quantitative metrics would strengthen the evidence presented. In the revised manuscript we will add chi-squared values, residual statistics, and comparison tables quantifying the agreement with data for the [O/Fe]-[Fe/H] distribution, stellar surface-density profile, and abundance gradients. revision: yes
Circularity Check
No circularity: Lagrangian integrals are independent math; two-infall is explicit assumption and v is tuned to external data
full rationale
The derivation reduces abundances to explicit one-dimensional integrals along characteristics whose source terms are the two-infall accretion history (an external premise stated in the abstract and methods). The value v = 1.5 km/s is introduced as a free parameter whose effect is then compared to observed [O/Fe]-[Fe/H] distributions, gradients, and surface-density profiles; the match is therefore a fit to independent benchmarks rather than a quantity defined from those same benchmarks. No self-citation is invoked as a uniqueness theorem, no ansatz is smuggled, and no fitted quantity is relabeled a first-principles prediction. The model is therefore self-contained against external data.
Axiom & Free-Parameter Ledger
free parameters (1)
- radial gas flow velocity v =
1.5 km/s
axioms (1)
- domain assumption Two-infall scenario accurately describes the gas accretion history of the Milky Way disc
read the original abstract
Chemical abundance patterns result from the interplay between gas accretion, star formation, and radial mixing of gas and stars. Disentangling these processes is crucial to recover the mechanisms shaping the formation and evolution of galaxies. We model the chemical evolution of the Galactic disc in the presence of radial gas flows, to assess their impact on the [O/Fe]-[Fe/H] abundance patterns and on the radial gradients of [Fe/H] and [O/H]. We develop fast, semi-analytic solutions for the gas surface mass density and the abundances of alpha-elements and iron, accounting for radial gas flows and chemical enrichment from core-collapse and Type Ia supernovae. The model follows a Lagrangian approach, using the method of characteristics, reducing the solutions to one-dimensional integrals. We apply our model to the Milky Way disc assuming a two-infall scenario. When radial gas flows are present, the chemical abundances of the gas at a given radius result from its whole inward journey in the disc, reflecting the star formation and accretion experienced at every radius it crossed. The integrated stellar mass along the characteristic is lower than the local value by up to an order of magnitude at v = 1.5 km/s. Models with mild flows of v = 1.5 km/s reproduce simultaneously the observed [O/Fe]-[Fe/H] distribution across the disc, the present-day stellar surface-density profile, and the [Fe/H] and [O/H] gradients, improving also the agreement with the observed age-abundance relations. The stellar mass formed per Type Ia supernova sets the [O/Fe] ratio and departs from its in-situ value by up to ~50 per cent, making the alpha-enhancement the quantity on which radial flows leave their strongest signature. Following the gas along its trajectory is essential to recover the correct enrichment history even for models with mild radial gas inflows. The code is made publicly available.
Figures
Reference graph
Works this paper leans on
- [1]
-
[2]
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
work page 2009
-
[3]
Bell, E. F., McIntosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289
work page 2003
-
[4]
Bensby, T., Bergemann, M., Rybizki, J., et al. 2019, The Messenger, 175, 35
work page 2019
- [5]
- [6]
- [7]
-
[8]
Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, 154, 28
work page 2017
-
[9]
Buder, S., Asplund, M., Duong, L., et al. 2018, MNRAS, 478, 4513
work page 2018
- [10]
- [11]
-
[12]
Chiappini, C., Minchev, I., Starkenburg, E., et al. 2019, The Messenger, 175, 30
work page 2019
-
[13]
2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, V ol
Cirasuolo, M., Afonso, J., Carollo, M., et al. 2014, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, V ol. 9147, Ground- based and Airborne Instrumentation for Astronomy V , ed. S. K. Ramsay, I. S. McLean, & H. Takami, 91470N de Jong, R. S., Agertz, O., Berbel, A. A., et al. 2019, The Messenger, 175, 3
work page 2014
-
[14]
Deason, A. J. & Belokurov, V . 2024, New A Rev., 99, 101706
work page 2024
-
[15]
Dubay, L. O., Johnson, J. A., Johnson, J. W., & Roberts, J. D. 2026, ApJ, 1000, 244
work page 2026
-
[16]
Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306
work page 2013
- [17]
-
[18]
1998, A&A, 338, 161 Gaia Collaboration, Vallenari, A., Brown, A
Fuhrmann, K. 1998, A&A, 338, 161 Gaia Collaboration, Vallenari, A., Brown, A. G. A., et al. 2023, A&A, 674, A1
work page 1998
-
[19]
A., Mucciarelli, A., Origlia, L., et al
Gonzalez, O. A., Mucciarelli, A., Origlia, L., et al. 2020, The Messenger, 180, 18
work page 2020
-
[20]
Goodman, J. & Weare, J. 2010, Communications in Applied Mathematics and Computational Science, 5, 65
work page 2010
-
[21]
R., Bland-Hawthorn, J., Sharma, S., et al
Hayden, M. R., Bland-Hawthorn, J., Sharma, S., et al. 2020, MNRAS, 493, 2952
work page 2020
-
[22]
Hayden, M. R., Bovy, J., Holtzman, J. A., et al. 2015, ApJ, 808, 132
work page 2015
-
[23]
R., Recio-Blanco, A., de Laverny, P., et al
Hayden, M. R., Recio-Blanco, A., de Laverny, P., et al. 2018, A&A, 609, A79
work page 2018
-
[24]
Iwamoto, K., Brachwitz, F., Nomoto, K., et al. 1999, ApJS, 125, 439
work page 1999
- [25]
- [26]
-
[27]
Johnson, J. W., Weinberg, D. H., Vincenzo, F., et al. 2021, MNRAS, 508, 4484
work page 2021
-
[28]
Kennicutt, Jr., R. C. 1998, ApJ, 498, 541
work page 1998
- [29]
-
[30]
Kubryk, M., Prantzos, N., & Athanassoula, E. 2015, A&A, 580, A126
work page 2015
-
[31]
Lacey, C. G. & Fall, S. M. 1985, ApJ, 290, 154
work page 1985
-
[32]
Leung, H. W. & Bovy, J. 2019, MNRAS, 483, 3255
work page 2019
-
[33]
Licquia, T. C. & Newman, J. A. 2015, ApJ, 806, 96
work page 2015
-
[34]
2015, Research in Astronomy and As- trophysics, 15, 1095
Luo, A.-L., Zhao, Y .-H., Zhao, G., et al. 2015, Research in Astronomy and As- trophysics, 15, 1095
work page 2015
-
[35]
Mackereth, J. T., Bovy, J., Leung, H. W., et al. 2019, MNRAS, 489, 176
work page 2019
-
[36]
Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94
work page 2017
- [37]
- [38]
- [39]
-
[40]
Matteucci, F. & Greggio, L. 1986, A&A, 154, 279 Mollá, M., Díaz, Á. I., Gibson, B. K., Cavichia, O., & López-Sánchez, Á.-R. 2016, MNRAS, 462, 1329
work page 1986
- [41]
- [42]
-
[43]
Palla, M., Matteucci, F., Spitoni, E., Vincenzo, F., & Grisoni, V . 2020, MNRAS, 498, 1710
work page 2020
- [44]
- [45]
-
[46]
2022, A&A, 666, A121 Schönrich, R
Randich, S., Gilmore, G., Magrini, L., et al. 2022, A&A, 666, A121 Schönrich, R. & Binney, J. 2009, MNRAS, 396, 203
work page 2022
-
[47]
Sharma, S., Hayden, M. R., & Bland-Hawthorn, J. 2021, MNRAS, 507, 5882
work page 2021
- [48]
- [49]
-
[50]
Spitoni, E., Matteucci, F., & Marcon-Uchida, M. M. 2013, A&A, 551, A123
work page 2013
-
[51]
Spitoni, E., Romano, D., Matteucci, F., & Ciotti, L. 2015, ApJ, 802, 129
work page 2015
-
[52]
Spitoni, E., Silva Aguirre, V ., Matteucci, F., Calura, F., & Grisoni, V . 2019, A&A, 623, A60
work page 2019
-
[53]
Spitoni, E., Verma, K., Silva Aguirre, V ., & Calura, F. 2020, A&A, 635, A58
work page 2020
-
[54]
Spitoni, E., Verma, K., Silva Aguirre, V ., et al. 2021, A&A, 647, A73
work page 2021
- [55]
-
[56]
Tinsley, B. M. 1980, Fund. Cosmic Phys., 5, 287
work page 1980
-
[57]
Vincenzo, F., Matteucci, F., Belfiore, F., & Maiolino, R. 2016, MNRAS, 455, 4183
work page 2016
-
[58]
Vincenzo, F., Matteucci, F., & Spitoni, E. 2017, MNRAS, 466, 2939
work page 2017
-
[59]
Vincenzo, F., Thompson, T. A., Weinberg, D. H., et al. 2021, MNRAS, 508, 3499
work page 2021
- [60]
-
[61]
Weinberg, D. H., Griffith, E. J., Johnson, J. W., & Thompson, T. A. 2024, ApJ, 973, 122
work page 2024
-
[62]
Weinberg, D. H., Holtzman, J. A., Hasselquist, S., et al. 2019, ApJ, 874, 102
work page 2019
-
[63]
Weinberg, D. H., Holtzman, J. A., Johnson, J. A., et al. 2022, ApJS, 260, 32
work page 2022
-
[64]
2019, ApJS, 245, 34 Article number, page 15
Xiang, M., Ting, Y .-S., Rix, H.-W., et al. 2019, ApJS, 245, 34 Article number, page 15
work page 2019
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.