Pith. sign in

REVIEW 2 major objections 1 minor 64 references

Chemical evolution of the Milky Way disc with radial gas flows: a Lagrangian approach

T0 review · 2 major / 1 minor · reviewed 2026-07-01 · grok-4.3

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

Pith's one-line read 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.

desk verdict The Lagrangian method of characteristics is a clean technical step for radial flows, but the simultaneous fits rest on the two-infall accretion history. read the letter →

arxiv 2606.31161 v1 pith:2U7UJXKH submitted 2026-06-30 astro-ph.GA

Fiorenzo Vincenzo This is my paper · ORCID
classification astro-ph.GA
keywords chemicalevolutionMilkyWaydiscradialgasflowsLagrangianapproachabundancegradientstwo-infallscenarioalpha-elementsTypeIasupernovae
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 paper builds a semi-analytic chemical evolution model that follows gas parcels inward under radial flows using the method of characteristics. Gas at any radius carries the enrichment history from every radius it crossed, so the stellar mass formed along each trajectory is lower than the local in-situ value. With mild inflows at 1.5 km/s the same model simultaneously fits the observed abundance plane, the present-day stellar surface density, the [Fe/H] and [O/H] gradients, and improves the age-abundance relations. The [O/Fe] ratio is the quantity most sensitive to the flows because the stellar mass per Type Ia supernova changes by up to 50 percent along the path.

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.

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.

Watch

Extended reading notes

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.

Load-bearing premise

The two-infall scenario for gas accretion history is assumed to hold for the Milky Way disc.

Editorial extensions

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.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • 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.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 1 minor

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)
  1. [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.
  2. [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)
  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

2 responses · 0 unresolved

We thank the referee for the constructive comments. We address each major comment below.

read point-by-point responses

  1. 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

  2. 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

0 steps flagged · score 0.0 of 10

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.

Assumptions & free parameters 1 free parameters · 1 assumptions · 0 invented entities

The model rests on the two-infall accretion history and a chosen radial flow speed; no new particles or forces are introduced.

free parameters (1)
  • radial gas flow velocity v = 1.5 km/s
    Chosen at 1.5 km/s to reproduce multiple Milky Way observables simultaneously
assumptions (1)
  • domain assumption Two-infall scenario accurately describes the gas accretion history of the Milky Way disc
    Invoked to set the boundary conditions for the Lagrangian characteristics

how reviews work

0 comments
Cite this review

Pith. "Pith review of Chemical evolution of the Milky Way disc with radial gas flows: a Lagrangian approach." pith.science (2026). https://pith.science/paper/2U7UJXKH

@misc{pith2026260631161,
  author       = {Pith},
  title        = {Pith review of: Chemical evolution of the Milky Way disc with radial gas flows: a Lagrangian approach},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/2U7UJXKH}},
  note         = {Machine review of arXiv:2606.31161}
}
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

Figures reproduced from arXiv: 2606.31161 by the authors.

Figure 1
Figure 1. The [O/Fe]–[Fe/H] abundance plane in eight radial bins spanning 3.5 ≤ R/kpc < 11.5. The background colour map shows the stellar number density of APOGEE DR17 stars estimated via a 2D Gaussian kernel density estimator, with regions below 2% of the peak density shown in white. Stars are selected within |z| ≤ 0.5 kpc of the Galactic plane. The coloured tracks show the predictions of the chemical evolution model for con… view at source ↗
Figure 2
Figure 2. Same as Fig [PITH_FULL_IMAGE:figures/full_fig_p009_2.png] view at source ↗
Figure 3
Figure 3. Chemical evolution tracks in the [O/Fe]–[Fe/H] plane of the gas residing at fixed Galactocentric radii R = 6, 8, 10 kpc (rows) for constant inflow velocities v = 0.5, 1.0, 1.5 km s−1 (columns). Each track is the Eulerian abundance history of the gas located at R, coloured by the radius ξ from which that gas originated, i.e. the Lagrangian label of the characteristic passing through R at each time (ξ = R + v τ for a … view at source ↗
Figures from the paper (6 more)
Figure 4
Figure 4. Figure 4: Left: present-day stellar surface-density profile [PITH_FULL_IMAGE:figures/full_fig_p011_4.png]
Figure 5
Figure 5. Figure 5: Age–abundance relations across the disc. Each column corresponds to a Galactocentric ring (labelled on top), and the three [PITH_FULL_IMAGE:figures/full_fig_p011_5.png]
Figure 6
Figure 6. Figure 6: Radial gradients of [Fe/H] (top row) and [O/H] (bottom row) for the OCCAM sample of open clusters (Myers et al. 2022), divided into the same four age bins used in their analysis (black stars with error bars; for [O/H] the uncertainty is propagated from [Fe/H] and [O/Fe…
Figure 7
Figure 7. Figure 7: Top row: the total stellar mass formed Mchar ⋆ , the total number of Type Ia SNe N char Ia , and the total accreted gas mass Mchar inf , integrated along the characteristic of the gas that reaches each radius R today (equations 26), as a function of R, for inflow veloc…
Figure 8
Figure 8. Figure 8: Ratio between the quantities integrated along the flow characteristic, [PITH_FULL_IMAGE:figures/full_fig_p013_8.png]
Figure 9
Figure 9. Figure 9: Effect of radial gas flows on the ratios that drive the abundances, shown as the ratio between each driver evaluated along the flow characteristic and the same driver that would be measured locally at the same radius, as a function of Galactocentric radius R, for inflo…

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

64 extracted references · 64 canonical work pages

  1. [1]

    2022, ApJS, 259, 35

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

    work page 2022

  2. [2]

    J., & Scott, P

    Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481

    work page 2009

  3. [3]

    F., McIntosh, D

    Bell, E. F., McIntosh, D. H., Katz, N., & Weinberg, M. D. 2003, ApJS, 149, 289

    work page 2003

  4. [4]

    2019, The Messenger, 175, 35

    Bensby, T., Bergemann, M., Rybizki, J., et al. 2019, The Messenger, 175, 35

    work page 2019

  5. [5]

    2003, A&A, 410, 527

    Bensby, T., Feltzing, S., & Lundström, I. 2003, A&A, 410, 527

    work page 2003

  6. [6]

    & Schönrich, R

    Bilitewski, T. & Schönrich, R. 2012, MNRAS, 426, 2266

    work page 2012

  7. [7]

    & Gerhard, O

    Bland-Hawthorn, J. & Gerhard, O. 2016, ARA&A, 54, 529

    work page 2016

  8. [8]

    R., Bershady, M

    Blanton, M. R., Bershady, M. A., Abolfathi, B., et al. 2017, AJ, 154, 28

    work page 2017

Show all 64 references
  1. [9]

    2018, MNRAS, 478, 4513

    Buder, S., Asplund, M., Duong, L., et al. 2018, MNRAS, 478, 4513

    2018

  2. [10]

    E., et al

    Buder, S., Kos, J., Wang, X. E., et al. 2025, PASA, 42, e051

    2025

  3. [11]

    1997, ApJ, 477, 765

    Chiappini, C., Matteucci, F., & Gratton, R. 1997, ApJ, 477, 765

    1997

  4. [12]

    2019, The Messenger, 175, 30

    Chiappini, C., Minchev, I., Starkenburg, E., et al. 2019, The Messenger, 175, 30

    2019

  5. [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., Age...

    2014

  6. [14]

    Deason, A. J. & Belokurov, V . 2024, New A Rev., 99, 101706

    2024

  7. [15]

    O., Johnson, J

    Dubay, L. O., Johnson, J. A., Johnson, J. W., & Roberts, J. D. 2026, ApJ, 1000, 244

    2026

  8. [16]

    W., Lang, D., & Goodman, J

    Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306

    2013

  9. [17]

    & Bland-Hawthorn, J

    Freeman, K. & Bland-Hawthorn, J. 2002, ARA&A, 40, 487

    2002

  10. [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

    1998

  11. [19]

    A., Mucciarelli, A., Origlia, L., et al

    Gonzalez, O. A., Mucciarelli, A., Origlia, L., et al. 2020, The Messenger, 180, 18

    2020

  12. [20]

    & Weare, J

    Goodman, J. & Weare, J. 2010, Communications in Applied Mathematics and Computational Science, 5, 65

    2010

  13. [21]

    R., Bland-Hawthorn, J., Sharma, S., et al

    Hayden, M. R., Bland-Hawthorn, J., Sharma, S., et al. 2020, MNRAS, 493, 2952

    2020

  14. [22]

    R., Bovy, J., Holtzman, J

    Hayden, M. R., Bovy, J., Holtzman, J. A., et al. 2015, ApJ, 808, 132

    2015

  15. [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

    2018

  16. [24]

    1999, ApJS, 125, 439

    Iwamoto, K., Brachwitz, F., Nomoto, K., et al. 1999, ApJS, 125, 439

    1999

  17. [25]

    C., Dalton, G

    Jin, S., Trager, S. C., Dalton, G. B., et al. 2024, MNRAS, 530, 2688

    2024

  18. [26]

    Johnson, J. W. 2025, arXiv e-prints, arXiv:2510.05223

    2025

  19. [27]

    W., Weinberg, D

    Johnson, J. W., Weinberg, D. H., Vincenzo, F., et al. 2021, MNRAS, 508, 4484

    2021

  20. [28]

    Kennicutt, Jr., R. C. 1998, ApJ, 498, 541

    1998

  21. [29]

    A., & Gilmore, G

    Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545

    1993

  22. [30]

    2015, A&A, 580, A126

    Kubryk, M., Prantzos, N., & Athanassoula, E. 2015, A&A, 580, A126

    2015

  23. [31]

    Lacey, C. G. & Fall, S. M. 1985, ApJ, 290, 154

    1985

  24. [32]

    Leung, H. W. & Bovy, J. 2019, MNRAS, 483, 3255

    2019

  25. [33]

    Licquia, T. C. & Newman, J. A. 2015, ApJ, 806, 96

    2015

  26. [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

    2015

  27. [35]

    T., Bovy, J., Leung, H

    Mackereth, J. T., Bovy, J., Leung, H. W., et al. 2019, MNRAS, 489, 176

    2019

  28. [36]

    R., Schiavon, R

    Majewski, S. R., Schiavon, R. P., Frinchaboy, P. M., et al. 2017, AJ, 154, 94

    2017

  29. [37]

    & Graur, O

    Maoz, D. & Graur, O. 2017, ApJ, 848, 25

    2017

  30. [38]

    2014, ARA&A, 52, 107

    Maoz, D., Mannucci, F., & Nelemans, G. 2014, ARA&A, 52, 107

    2014

  31. [39]

    & Francois, P

    Matteucci, F. & Francois, P. 1989, MNRAS, 239, 885

    1989

  32. [40]

    & Greggio, L

    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

    1986

  33. [41]

    2013, MNRAS, 435, 2918

    Mott, A., Spitoni, E., & Matteucci, F. 2013, MNRAS, 435, 2918

    2013

  34. [42]

    2022, AJ, 164, 85

    Myers, N., Donor, J., Spoo, T., et al. 2022, AJ, 164, 85

    2022

  35. [43]

    2020, MNRAS, 498, 1710

    Palla, M., Matteucci, F., Spitoni, E., Vincenzo, F., & Grisoni, V . 2020, MNRAS, 498, 1710

    2020

  36. [44]

    & Fraternali, F

    Pezzulli, G. & Fraternali, F. 2016, MNRAS, 455, 2308

    2016

  37. [45]

    & Chiosi, C

    Portinari, L. & Chiosi, C. 2000, A&A, 355, 929

    2000

  38. [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

    2022

  39. [47]

    R., & Bland-Hawthorn, J

    Sharma, S., Hayden, M. R., & Bland-Hawthorn, J. 2021, MNRAS, 507, 5882

    2021

  40. [48]

    2017, A&A, 605, A38

    Spitoni, E., Gioannini, L., & Matteucci, F. 2017, A&A, 605, A38

    2017

  41. [49]

    & Matteucci, F

    Spitoni, E. & Matteucci, F. 2011, A&A, 531, A72

    2011

  42. [50]

    Spitoni, E., Matteucci, F., & Marcon-Uchida, M. M. 2013, A&A, 551, A123

    2013

  43. [51]

    2015, ApJ, 802, 129

    Spitoni, E., Romano, D., Matteucci, F., & Ciotti, L. 2015, ApJ, 802, 129

    2015

  44. [52]

    2019, A&A, 623, A60

    Spitoni, E., Silva Aguirre, V ., Matteucci, F., Calura, F., & Grisoni, V . 2019, A&A, 623, A60

    2019

  45. [53]

    2020, A&A, 635, A58

    Spitoni, E., Verma, K., Silva Aguirre, V ., & Calura, F. 2020, A&A, 635, A58

    2020

  46. [54]

    2021, A&A, 647, A73

    Spitoni, E., Verma, K., Silva Aguirre, V ., et al. 2021, A&A, 647, A73

    2021

  47. [55]

    & Price, K

    Storn, R. & Price, K. 1997, Journal of Global Optimization, 11, 341

    1997

  48. [56]

    Tinsley, B. M. 1980, Fund. Cosmic Phys., 5, 287

    1980

  49. [57]

    2016, MNRAS, 455, 4183

    Vincenzo, F., Matteucci, F., Belfiore, F., & Maiolino, R. 2016, MNRAS, 455, 4183

    2016

  50. [58]

    2017, MNRAS, 466, 2939

    Vincenzo, F., Matteucci, F., & Spitoni, E. 2017, MNRAS, 466, 2939

    2017

  51. [59]

    A., Weinberg, D

    Vincenzo, F., Thompson, T. A., Weinberg, D. H., et al. 2021, MNRAS, 508, 3499

    2021

  52. [60]

    E., et al

    Virtanen, P., Gommers, R., Oliphant, T. E., et al. 2020, Nature Methods, 17, 261

    2020

  53. [61]

    H., Griffith, E

    Weinberg, D. H., Griffith, E. J., Johnson, J. W., & Thompson, T. A. 2024, ApJ, 973, 122

    2024

  54. [62]

    H., Holtzman, J

    Weinberg, D. H., Holtzman, J. A., Hasselquist, S., et al. 2019, ApJ, 874, 102

    2019

  55. [63]

    H., Holtzman, J

    Weinberg, D. H., Holtzman, J. A., Johnson, J. A., et al. 2022, ApJS, 260, 32

    2022

  56. [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

    2019

Pith tools

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