pith. sign in

arxiv: 2511.10582 · v1 · submitted 2025-11-13 · 🌌 astro-ph.GA

Co-evolution of baryons and dark matter halos of LYRA dwarf galaxies

Joaquin Sureda , Shaun T. Brown , Azadeh Fattahi , Thales Gutcke , Sownak Bose , Jessica E. Doppel , R\"udiger Pakmor This is my paper

Pith reviewed 2026-05-17 22:14 UTC · model grok-4.3

classification 🌌 astro-ph.GA
keywords dwarf galaxiesreionizationstar formation historydark matter halosmetallicity distributionsgalaxy shapescosmological simulations
0
0 comments X

The pith

Dwarf galaxies that restart star formation after reionization develop rounder shapes and distinct metal distributions compared to those that remain quenched.

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

The paper examines six high-resolution simulations of dwarf galaxies to track how their stars formed before and after reionization. Some galaxies, called rejuvenators, restart star formation later, while others, relics, stay dark. It shows that the amount of post-reionization star formation links to the galaxies' metal content and their overall shape at the present day. Rejuvenators have higher median iron abundance and rounder distributions, and their dark matter halos also become rounder. This suggests that the history of star formation leaves observable marks on both the stars and the dark matter in these small galaxies.

Core claim

In the LYRA simulations, the fraction of stellar mass formed after reionization, denoted f_post-reio^star, correlates with the metallicity distribution such that reionization relics have lower median [Fe/H] and a more prominent low-metallicity tail. The shape of the stellar component at z=0 is more spherical for rejuvenators than for relics, with the difference emerging after reionization. Inner dark matter halos are also rounder in rejuvenators. All simulated halos become rounder than their dark-matter-only counterparts, with stronger evolution in rejuvenators. No correlation is found between star formation activity and the formation of shallow dark matter density cores.

What carries the argument

The fraction of stars formed post-reionization (f_post-reio^star), which distinguishes rejuvenators from reionization relics and drives correlations with metallicity distributions and galaxy shapes.

If this is right

  • Reionization relics exhibit lower median iron abundances and more prominent low-metallicity tails than rejuvenators.
  • Rejuvenators display more spherical stellar distributions at z=0 due to post-reionization star formation activity.
  • Inner dark matter halos are rounder in rejuvenators than in relics, with all halos rounder than dark-matter-only counterparts.
  • The scatter in the low-mass mass-metallicity relation correlates with star formation histories and galaxy shapes.
  • These trends offer predictions testable with upcoming observational data on dwarf galaxies.

Where Pith is reading between the lines

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

  • Observers could use measurements of galaxy sphericity and metallicity distributions in dwarfs to estimate the fraction of stars formed after reionization.
  • Baryonic processes appear to influence dark matter halo shapes on small scales beyond core formation.
  • The co-evolution pattern may help interpret shape variations in observed samples of low-mass galaxies.

Load-bearing premise

The six simulated galaxies and the subgrid physics prescriptions accurately represent the real population of dwarf galaxies and the effects of reionization, feedback, and mergers without major systematic biases.

What would settle it

Finding no difference in median [Fe/H], low-metallicity tail prominence, or shape between observed dwarf galaxies with and without post-reionization star formation in large samples.

Figures

Figures reproduced from arXiv: 2511.10582 by Azadeh Fattahi, Jessica E. Doppel, Joaquin Sureda, R\"udiger Pakmor, Shaun T. Brown, Sownak Bose, Thales Gutcke.

Figure 1
Figure 1. Figure 1: presents the stellar mass-halo mass1 (SMHM) relation for the six galaxies in the sample, indicated by the diamonds. The LYRA halos, which span a halo mass range of 7 × 108 − 5 × 109 M⊙, host galaxies with stellar masses in the range 4 × 105 − 1 × 107 M⊙ at 𝑧 = 0. In the same figure, there are samples of other simulated dwarf galaxies from other projects, namely, FIRE-2 (Hopkins et al. 2018), NIHAO (Wang et… view at source ↗
Figure 2
Figure 2. Figure 2: Stellar (top) and cold gas (mid) mass growth vs z at different radius, where cold gas refers to gas cells with temperature 𝑇gas < 103 K. The bottom panel shows the SFH for each halo in 50 Myr bins. The left panels correspond to rejuvenated galaxies whereas the right ones are reionization relics. The grey vertical band illustrates reionization while the red translucent vertical lines indicate the correspond… view at source ↗
Figure 4
Figure 4. Figure 4: Ellipsoidal axis ratios for the stellar component of the LYRA dwarf galaxies at 𝑧 = 0, measured enclosed to the stellar half mass radius, 𝑟1/2 (see Sec. 2.3 for the details). The colour of the symbols represents their SFH category, rejuvenator (blue) or relic (red), and the size of each symbol is proportional to the fraction of stellar mass formed after reionization. This figure illustrates that SFH activi… view at source ↗
Figure 6
Figure 6. Figure 6: Evolution of the gas ellipsoidal axis ratios for the rejuvenated halos, A, F, and D. Shapes are measured in the same way as in [PITH_FULL_IMAGE:figures/full_fig_p008_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: Axes ratios for the dark matter halo measured at a fixed radius of 𝑟1/2,★. The coloured points represent the rejuvenated and relic dwarfs in the hydro runs whereas the cross symbols show the corresponding DM-only counterparts. The symbol sizes are proportional to the fraction of stellar mass formed after reionization. are displayed as circles, and are connected to their DM-only coun￾terparts, shown as cros… view at source ↗
Figure 9
Figure 9. Figure 9: Left: Dark matter density profiles for all halos. The panel at the bottom shows the ratio between the profile in the Hydro run over the one in the DM-only run. Lines are drawn fainter below the Power et al. (2003) convergence criterion. The inset shows the region corresponding to the inner 100 pc to highlight the contraction feature of Halo A. Right: Inner slope of the density profile, measured between 1 a… view at source ↗
Figure 10
Figure 10. Figure 10: Evolution of total gas mass to DM mass ratio within 1 kpc for all simulated dwarfs. This shows that gas never dominates the mass budget in any of the galaxies, hence explaining the lack of a dark matter core. baryonic evolution in [PITH_FULL_IMAGE:figures/full_fig_p011_10.png] view at source ↗
Figure 12
Figure 12. Figure 12: From top to bottom the panels show the circular velocity profile (left) and surface density of stars and gas (right) of Halo A on different snapshots during the starburst at 𝑧 ∼ 5. This illustrates the contraction of the dark matter halo being attributed to the condensation of the gas in the centre followed by a strong star formation event. scale, mostly on the Milky Way mass regime. The difference in the… view at source ↗
read the original abstract

We use the extremely high-resolution ($m_{\rm bary}=4\rm{M}_\odot$) LYRA cosmological galaxy formation simulations of six dwarf galaxies with $M_{\rm 200c}\sim10^9\rm{M}_\odot$ at $z=0$ to investigate their stellar assembly histories. Based on the age of stars in these galaxies at $z=0$, $40-100\%$ of their stellar mass was formed by the time of reionization, when star formation (SF) abruptly shuts down. Depending on their halo mass evolution, some of the dwarfs reignite SF post-reionization (rejuvenators), while others remain quenched for the rest of cosmic time (reionization relics). However, the stellar mass of relics can still grow by more than $50\%$ through mergers post-reionisation. We find clear correlations between metallicity distributions of the galaxies and the fraction of stars formed post-reionization ($f_{\rm post-reio}^\star$) such that relics have lower median $\rm [Fe/H]$ with a more prominent low metallicity tail. Moreover, the shape of the galaxies at $z=0$ correlates with their $f_{\rm post-reio}^\star$, with rejuvenators showing more spherical stellar distribution than relics. This difference arises only post-reionization when rejuvenators become rounder with more SF activity. Similarly, the shape of dark matter (DM) halos in the inner regions display more spherical distributions in rejuvenators than in relics. The shape evolution shows that DM haloes in all galaxy formation simulations become rounder in comparison to their collision-less, DM-only counterparts. However, DM haloes of rejuvenators evolve more significantly. We do not find any correlation between SF activity and formation of shallow DM density cores in these galaxies. These predictions can be tested using upcoming observational data. In particular, our results indicate that the scatter in the mass-metallicity relation in the low mass regime is correlated with SF histories and the shape of galaxies.

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

1 major / 2 minor

Summary. The manuscript analyzes six high-resolution (m_bary = 4 M_⊙) LYRA zoom-in simulations of dwarf galaxies with M_200c ≈ 10^9 M_⊙ at z=0. Galaxies are classified as reionization relics (quenched after reionization) or rejuvenators (resuming star formation post-reionization) based on halo mass evolution and stellar assembly. The central results are empirical correlations: relics exhibit lower median [Fe/H] and a more prominent low-metallicity tail as a function of f_post-reio^star; rejuvenators display more spherical stellar and inner-DM distributions at z=0; DM halos in all hydro runs become rounder than their DM-only counterparts, with stronger evolution in rejuvenators. No correlation is reported between star-formation activity and shallow DM density cores. These trends are offered as testable predictions for the scatter in the low-mass mass-metallicity relation.

Significance. If the reported trends prove robust, the work would link post-reionization baryonic processes directly to observable chemo-morphological properties and to the shape evolution of inner dark-matter halos in the dwarf regime. The extreme mass resolution permits detailed tracking of individual stellar populations and merger-driven growth, which is a clear technical strength. The explicit, falsifiable predictions for metallicity-shape correlations provide a concrete bridge to upcoming resolved-star observations.

major comments (1)
  1. [§4 and §5] §4 (metallicity results) and §5 (shape results): All reported correlations (median [Fe/H], low-metallicity tail, stellar and inner-DM sphericity) are derived from trends across exactly six galaxies. The manuscript presents these as clear correlations without quantified statistical measures (e.g., Spearman coefficients with bootstrap uncertainties) or tests against modest variations in subgrid feedback efficiency or different zoom-in initial conditions. Because the relic/rejuvenator classification itself depends on the fixed reionization and feedback prescriptions, the load-bearing claim that these differences constitute general predictions for the dwarf population requires either larger statistics or explicit robustness checks.
minor comments (2)
  1. [§5] The definition of 'inner regions' for DM shape measurements and the precise algorithm used to compute axis ratios (inertia tensor, radial weighting, etc.) should be stated explicitly so that the reported rounding relative to DM-only runs can be reproduced or compared to other suites.
  2. [Abstract and §6] The abstract and conclusion state that the results 'can be tested using upcoming observational data' but do not name specific observables (e.g., metallicity distribution functions from resolved stars in Local Group dwarfs or shape measurements from IFU surveys) that would directly confront the predicted f_post-reio^star correlations.

Simulated Author's Rebuttal

1 responses · 1 unresolved

We thank the referee for their constructive comments and for recommending minor revision. We address the single major comment below, incorporating quantitative measures where feasible while being transparent about the inherent limitations of the current simulation suite.

read point-by-point responses

  1. Referee: [§4 and §5] §4 (metallicity results) and §5 (shape results): All reported correlations (median [Fe/H], low-metallicity tail, stellar and inner-DM sphericity) are derived from trends across exactly six galaxies. The manuscript presents these as clear correlations without quantified statistical measures (e.g., Spearman coefficients with bootstrap uncertainties) or tests against modest variations in subgrid feedback efficiency or different zoom-in initial conditions. Because the relic/rejuvenator classification itself depends on the fixed reionization and feedback prescriptions, the load-bearing claim that these differences constitute general predictions for the dwarf population requires either larger statistics or explicit robustness checks.

    Authors: We agree that the sample of six galaxies limits the statistical robustness of the reported trends and that the relic/rejuvenator classification is tied to the specific reionization and feedback implementation used. The small sample size is a direct consequence of the extreme baryonic mass resolution (4 M_⊙) required to track individual stellar populations and merger histories in cosmological zoom-ins. In the revised manuscript we have added Spearman rank correlation coefficients together with bootstrap uncertainties for the relations between f_post-reio^star and both the metallicity diagnostics (median [Fe/H] and low-metallicity tail) and the shape parameters (stellar and inner-DM sphericity). These quantitative measures are now reported in §§4 and 5. We have also clarified in the text that the trends are presented as testable predictions from this particular high-resolution suite rather than as statistically general results for the full dwarf population. Explicit robustness checks against changes in subgrid feedback efficiency or different zoom-in initial conditions would require a new set of simulations that lies outside the scope of the present study. revision: partial

standing simulated objections not resolved
  • We cannot expand the sample size or perform robustness tests against variations in feedback prescriptions and initial conditions within the current work owing to the prohibitive computational cost of the 4 M_⊙ resolution.

Circularity Check

0 steps flagged

No significant circularity; results are direct simulation outputs

full rationale

The paper reports empirical correlations and trends extracted from the outputs of the LYRA cosmological simulations run on six specific dwarf galaxies. The relic/rejuvenator classification follows directly from the simulated star-formation histories and halo-mass evolution; the reported differences in metallicity distributions, stellar shapes, and inner DM halo shapes are measured quantities from those runs rather than quantities derived via equations that reduce to the inputs by construction. No fitted parameters are renamed as predictions, no self-citation chains supply the central claims, and no ansatz or uniqueness theorem is invoked to force the results. The analysis is therefore self-contained against the simulation data themselves.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

The central claims rest on the accuracy of the LYRA cosmological hydrodynamical simulation framework, which incorporates standard cosmological assumptions plus numerous subgrid parameters for star formation and feedback that are typically calibrated to observations.

free parameters (1)
  • subgrid star formation and feedback parameters
    Standard in galaxy formation simulations; values chosen to reproduce observed galaxy properties at low masses.
axioms (1)
  • standard math Lambda-CDM cosmology and a specific reionization model
    Invoked to set the timing of reionization and the background evolution for the six dwarf halos.

pith-pipeline@v0.9.0 · 5701 in / 1322 out tokens · 37067 ms · 2026-05-17T22:14:39.926405+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

1 extracted references · 1 canonical work pages

  1. [1]

    V., 2015, ApJ, 804, 18 Agertz O., et al., 2020, MNRAS, 491, 1656 Allgood B., Flores R

    Agertz O., Kravtsov A. V., 2015, ApJ, 804, 18 Agertz O., et al., 2020, MNRAS, 491, 1656 Allgood B., Flores R. A., Primack J. R., Kravtsov A. V., Wechsler R. H., Faltenbacher A., Bullock J. S., 2006, MNRAS, 367, 1781 Andersson E. P., Rey M. P., Pontzen A., Cadiou C., Agertz O., Read J. I., Martin N. F., 2025, ApJ, 978, 129 Behroozi P. S., Wechsler R. H., C...

    work page arXiv 2015