pith. sign in

arxiv: 2512.14897 · v2 · submitted 2025-12-16 · 🌌 astro-ph.GA

Portrait of a Galaxy on FIRE: Is the α-bimodality a natural consequence of inside-out disc growth in a hierarchical formation scenario?

Mar\'ia Benito , Annaliina Aavik , Giuseppina Battaglia , Salvador Cardona-Barrero , Ele-Liis Evestus , Emma Fern\'andez-Alvar , Sven P\~oder , Heleri Ramler
show 2 more authors
Boris Deshev Elmo Tempel
This is my paper

Pith reviewed 2026-05-16 21:28 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords alpha bimodalitydisc galaxy formationchemical evolutioninside-out growthminor mergersstellar abundancesgas accretionhot gas corona
0
0 comments X p. Extension

The pith

A Milky Way-mass galaxy develops distinct high- and low-alpha stellar sequences through inside-out disc growth and minor gas accretion alone.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper tests whether the chemical split between high-alpha and low-alpha stars, long observed in the Milky Way disc, can arise in a galaxy that builds its disc gradually from inside out without major mergers or stars drifting radially across large distances. In the simulation a thick high-alpha disc assembles early during a chaotic phase of star formation; a later minor merger dilutes the surrounding hot gas reservoir, pauses the rise in iron abundance, and creates a temporary dip in star formation. Fresh gas brought in by subsequent minor mergers then fuels a low-alpha population that begins in the outer disc and flows inward, while incomplete mixing leaves chemical spreads even among stars born at the same radius. A sympathetic reader would care because this sequence suggests the observed bimodality is a routine outcome of ordinary hierarchical gas accretion rather than a signature of rare violent events.

Core claim

The simulated galaxy exhibits high and low-α sequences without having experienced major mergers nor significant radial migration. A high-α thick disc forms during the early chaotic clustering phase. Afterwards, as the star formation rate declines, a dip in the stellar number density appears, coinciding with the dilution of the galactic corona by a minor merger, which subsequently halts the rise of [Fe/H] in the disc. Later, accreted gas onto the disc from minor mergers mildly enhances the star formation rate and generates the low-α sequence in the outer disc, with radial inward flows of this material feeding the low-α inner disc. Even at fixed radii, newly formed stars retain a sizablespread

What carries the argument

The interplay of early chaotic star formation, minor-merger gas accretion, corona dilution, and radial gas flows that together set the timing and radial distribution of high-alpha and low-alpha star formation.

If this is right

  • The alpha bimodality can form without major mergers or significant radial migration.
  • A temporary pause in iron enrichment caused by corona dilution creates the gap between the two sequences.
  • Low-alpha stars first appear in the outer disc and reach the inner disc via inward gas flows.
  • Chemical abundance spreads persist at fixed radii because in-situ and infalling gas retain distinct compositions.
  • Vertical and radial gas flows redistribute metals and must be included in models of disc chemical evolution.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same inside-out accretion sequence could produce observable alpha bimodality in other disc galaxies of similar mass.
  • External galaxies showing clear alpha sequences but little evidence of major mergers would support this pathway over merger-driven explanations.
  • Simulations that vary the timing or mass of minor mergers could test how robust the bimodality is to changes in accretion history.

Load-bearing premise

The detailed gas accretion, dilution, and mixing behavior seen in this single simulated galaxy occurs in the same way in typical real Milky Way-mass galaxies.

What would settle it

A Milky Way-mass galaxy observed to lack any alpha bimodality despite showing no major mergers and a minor-accretion history similar to the simulation would falsify the claim.

Figures

Figures reproduced from arXiv: 2512.14897 by Annaliina Aavik, Boris Deshev, Ele-Liis Evestus, Elmo Tempel, Emma Fern\'andez-Alvar, Giuseppina Battaglia, Heleri Ramler, Mar\'ia Benito, Salvador Cardona-Barrero, Sven P\~oder.

Figure 1
Figure 1. Figure 1: Top: Stellar surface density profiles of the Romeo and Vintergatan (Agertz et al. 2021) simulated galaxies, together with the Milky Way profile taken from Lian et al. (2024) and the best-fitting morphology of McMillan (2011). The former Milky Way profile has been normalised to enclose a stellar mass of 4 × 1010 M⊙. The dot-dashed orange line depicts a falling expo￾nential profile with scale-length of 2.6 k… view at source ↗
Figure 2
Figure 2. Figure 2: One-dimensional age distribution of all disc star particles [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: Two-dimensional distribution of the azimuthal or rota [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: Top panel: two-dimensional number density of a sample of giants in Gaia DR3 crossmatched with APOGEE DR17 that spanned a Galactocentric range of 4-16 kpc (see text for details). The red dashed lines delimit the thick disc/bridge/low-α regions used to derive the kinematic properties in [PITH_FULL_IMAGE:figures/full_fig_p005_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: Time evolution of radial stellar surface density (top) and vertical density (bottom) profiles from [PITH_FULL_IMAGE:figures/full_fig_p006_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Top: two-dimensional number density of disc star particles within concentric rings with |z| < 1 kpc. The number density dis￾tribution is normalised within each ring. Bottom: Mean age distribution of disc stars. The red lines depict the chemical evolutionary track of stars formed from in-situ gas in each of the rings. In the case of the outer rings, the sequence is discontinuous, as there are periods when n… view at source ↗
Figure 9
Figure 9. Figure 9: Temperature versus number density distribution of the gas [PITH_FULL_IMAGE:figures/full_fig_p008_9.png] view at source ↗
Figure 7
Figure 7. Figure 7: Top panel: gas chemical tracks in the inner, central, and outer disc regions. Each track is constructed by, at each snapshot, selecting the gas within the corresponding region and computing its (mass weighted) average chemical properties. Grey dashed lines indicate schematically the three general trends exhibited by the chemical tracks. Bottom panels: dynamically-evolved stellar chemical tracks for the inn… view at source ↗
Figure 8
Figure 8. Figure 8: Chemical evolutionary tracks of stars, separated accord [PITH_FULL_IMAGE:figures/full_fig_p008_8.png] view at source ↗
Figure 10
Figure 10. Figure 10: Vertical (top panels) and radial (bottom panels) net (inflow minus outflow) gas flow rates through an inner (0 [PITH_FULL_IMAGE:figures/full_fig_p009_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Gas temperature at different lookback times. These are number density-weighted profiles along the perpendicular direc￾tion, enhancing the visibility of underdense regions by reducing the dominance of high-density areas. The reference frames are aligned with the principal axes of Romeo at the last simulated snapshot. The yellow circle marks the virial radius R200,c at each time. 5. Discussion & conclusions… view at source ↗
read the original abstract

The chemical dichotomy in the [$\alpha$/Fe]-[Fe/H] plane is a consequence of the complex processes underlying the formation and evolution of disc galaxies such as observed in the stellar Milky Way disc. We determine what can drive an $\alpha$-bimodality of the disc in a zoom-in hydrodynamical simulated galaxy which has had no major mergers and negligible radial migration. Using a Milky Way-mass galaxy from the FIRE-2 suite of simulations, we analyse gas flows in the disc together with its star formation and merger history, as well as the chemical evolution of the hot corona, to investigate their connection to transitions in the chemo-dynamical structure of the stellar disc and its radial distribution. The simulated galaxy exhibits high and low-$\alpha$ sequences without having experienced major mergers nor significant radial migration. A high-$\alpha$ thick disc forms during the early chaotic clustering phase. Afterwards, as the star formation rate declines, a dip in the stellar number density appears, coinciding with the dilution of the galactic corona by a minor merger, which subsequently halts the rise of [Fe/H] in the disc. Later, accreted gas onto the disc from minor mergers, mildly enhances the star formation rate and generates the low-$\alpha$ sequence in the outer disc, with radial inward flows of this material feeding the low-$\alpha$ inner disc. Furthermore, we find that even at fixed radii, newly formed stars retain a sizable spread in their chemical abundances, reflecting chemical differences between the in-situ and the infalling gas from which they formed, further indicating that instantaneous gas mixing is invalid. Understanding the chemical evolution of stellar discs requires accounting for their accretion merger history and interaction with the surrounding hot corona, as well as the vertical and radial gas flows that redistribute metals within the disc.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript analyzes a single Milky Way-mass galaxy (m12) from the FIRE-2 zoom-in hydrodynamical simulations that experiences no major mergers and negligible radial migration. It traces the high-α sequence to an early chaotic clustering phase of star formation, a subsequent dip in stellar density coinciding with dilution of the hot corona by a minor merger that stalls the rise in [Fe/H], and the low-α sequence to later gas accretion from minor mergers that mildly boosts star formation (primarily in the outer disc) with inward radial flows supplying the inner disc. The authors further report sizable [α/Fe] and [Fe/H] spreads at fixed radii arising from differences between in-situ and infalling gas, arguing that instantaneous mixing is invalid. The central claim is that α-bimodality is a natural outcome of inside-out disc growth plus hierarchical accretion and corona interactions.

Significance. If the sequence of events holds, the work supplies a concrete, merger-minor pathway for the observed high- and low-α sequences that does not invoke major mergers or substantial radial migration. It directly links gas flows, corona metallicity evolution, and radial abundance gradients to the chemo-dynamical structure of the disc, reinforcing the value of tracking sub-grid metal diffusion and accretion history in cosmological simulations. The demonstration that abundance spreads persist even at fixed radii provides a falsifiable signature against instantaneous mixing assumptions.

major comments (2)
  1. [§3 and §4] §3 (merger and star-formation history) and §4 (corona dilution): the pivotal minor-merger dilution event that produces the density dip and stalls [Fe/H] is shown for only this one accretion history. Because the bimodality is tied to the timing and mass ratio of that specific merger, the claim that the sequences form as a 'natural consequence' of inside-out growth requires explicit discussion of how sensitive the outcome is to variations in merger timing or to the sub-grid metal-diffusion coefficient; otherwise the result risks being realization-dependent.
  2. [§5] §5 (radial abundance distributions and gas flows): the reported spread in abundances at fixed radii is generated by the FIRE-2 sub-grid metal-mixing implementation. While the paper correctly notes that this spread contradicts instantaneous mixing, no test is presented of how the bimodality changes if the diffusion coefficient is varied within its plausible range; this leaves open whether the low-α sequence is robust or an artifact of the chosen sub-grid prescription.
minor comments (2)
  1. [Figure 3] Figure 3 (or equivalent radial profiles): the distinction between in-situ and accreted gas contributions to the low-α sequence should be shown with separate curves or shaded regions to make the inward-flow argument visually quantitative.
  2. [§4] Notation: the definition of 'corona dilution' and the precise metallicity threshold used to identify the minor-merger event should be stated explicitly in the text rather than only in a figure caption.

Simulated Author's Rebuttal

2 responses · 1 unresolved

We thank the referee for their constructive and detailed report. We address each major comment below, providing clarifications and noting revisions made to the manuscript.

read point-by-point responses

  1. Referee: [§3 and §4] §3 (merger and star-formation history) and §4 (corona dilution): the pivotal minor-merger dilution event that produces the density dip and stalls [Fe/H] is shown for only this one accretion history. Because the bimodality is tied to the timing and mass ratio of that specific merger, the claim that the sequences form as a 'natural consequence' of inside-out growth requires explicit discussion of how sensitive the outcome is to variations in merger timing or to the sub-grid metal-diffusion coefficient; otherwise the result risks being realization-dependent.

    Authors: We agree that the analysis is based on a single galaxy realization and that the timing of the minor merger is important for the exact sequence of events. The manuscript demonstrates that α-bimodality can emerge from inside-out growth and hierarchical accretion without major mergers or radial migration, rather than claiming this is the universal mechanism. In the revised manuscript we have added discussion in §6 noting that minor-merger-driven corona dilution is a recurring feature in FIRE-2 galaxies and that the overall pathway is expected to be robust across similar accretion histories, while acknowledging that the precise timing of the density dip may shift with different merger schedules. We have also noted that the fiducial metal-diffusion coefficient is fixed by the FIRE-2 calibration; varying it would require new simulations outside the present scope. revision: partial

  2. Referee: [§5] §5 (radial abundance distributions and gas flows): the reported spread in abundances at fixed radii is generated by the FIRE-2 sub-grid metal-mixing implementation. While the paper correctly notes that this spread contradicts instantaneous mixing, no test is presented of how the bimodality changes if the diffusion coefficient is varied within its plausible range; this leaves open whether the low-α sequence is robust or an artifact of the chosen sub-grid prescription.

    Authors: We thank the referee for this point. The abundance spreads at fixed radii are a direct outcome of the explicit sub-grid metal mixing in FIRE-2, which permits chemical differences between in-situ and accreted gas and thereby falsifies instantaneous mixing. The low-α sequence itself is driven by the timing of gas accretion from minor mergers and the associated star-formation enhancement, not by the precise diffusion strength. In the revised §5 we have added text clarifying this distinction and included a caveat that the quantitative width of the sequences could change with different diffusion coefficients, while the qualitative bimodality tied to the accretion history is expected to persist. A systematic variation study is not performed here as it would require a new simulation suite. revision: partial

standing simulated objections not resolved
  • Quantitative sensitivity of the α-bimodality to variations in minor-merger timing or the sub-grid metal-diffusion coefficient, which would require additional simulations beyond the single m12 realization analyzed in the manuscript.

Circularity Check

0 steps flagged

No circularity; results are direct outputs of hydrodynamical simulation post-processing

full rationale

The paper analyzes gas flows, star formation history, merger events, and chemical evolution directly from a single FIRE-2 zoom-in run. The α-bimodality, high-α thick disc, density dip, corona dilution, and low-α sequence are presented as emergent features of the simulated accretion and inside-out growth, without any fitted parameters tuned to match observations or any equation that reduces the outcome to a self-defined input. No self-citations are load-bearing for the central claim, and the derivation chain contains no self-definitional, fitted-prediction, or ansatz-smuggling steps.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on the fidelity of the FIRE-2 sub-grid physics for star formation, feedback, and metal transport; no new free parameters are introduced or fitted to produce the bimodality.

axioms (1)
  • domain assumption The FIRE-2 hydrodynamical prescriptions and metal-mixing model correctly capture gas accretion, corona dilution, and radial flows on galactic scales.
    Invoked throughout the analysis of gas flows and chemical evolution in sections describing the simulation outputs.

pith-pipeline@v0.9.0 · 5696 in / 1395 out tokens · 52433 ms · 2026-05-16T21:28:41.939486+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

3 extracted references · 3 canonical work pages

  1. [1]

    arXiv:2503.19536

    Agertz, O., Renaud, F., Feltzing, S., et al. 2021, MNRAS, 503, 5826 Andrae, R., Fouesneau, M., Sordo, R., et al. 2023, A&A, 674, A27 Ansar, S., Pearson, S., Sanderson, R. E., et al. 2025, ApJ, 978, 37 Barbani, F., Pascale, R., Marinacci, F., et al. 2023, MNRAS, 524, 4091 Belokurov, V . & Kravtsov, A. 2022, MNRAS, 514, 689 Belokurov, V . & Kravtsov, A. 202...

    work page arXiv 2021

  2. [2]

    In addition, we select only stars with 3500<T eff[K]<6500 and 1<logg[dex]<3.5

    We also excluded those objects with theapogee1_target1,apogee1_target2, apogee2_target1, andapogee2_target2classifying them as clus- ters, streams, dwarf galaxies, sky, telluric, binaries, radial veloc- ity variables, the bar, extended objects, and the TriAnd, GASS, and A13 disc structures. In addition, we select only stars with 3500<T eff[K]<6500 and 1<l...

    work page 2019

  3. [3]

    Cells withv z z>0 (v z z<0) were counted as outflows (inflows)

    In partic- ular, for each annular ring defined byR min <R<R max and |z|<2 kpc, vertical fluxes were measured through the top and bottom surfaces by selecting cells withinR min <R<R max and z∈[±2−∆/2,±2+ ∆/2]. Cells withv z z>0 (v z z<0) were counted as outflows (inflows). Additionally, radial fluxes were measured through the outer cylindrical surface by s...

    work page 2023