Pith. sign in

REVIEW 4 major objections 6 minor 2 cited by

One set of equations links star formation history shape to feedback

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 →

A unified empirical equation, fitted separately to three simulation suites, links galaxy star formation history shape to halo mass, baryon fraction, black hole mass, and feedback strength.

T0 review reviewed 2026-08-05 challenge →

load-bearing objection A careful and useful survey of how feedback shapes SFHs across three CAMELS models, but the headline unified equations are in-sample fits to an emulator and should not be treated as validated predictions. the 4 major comments →

arxiv 2508.21152 v1 pith:MXCGDHLN submitted 2025-08-28 astro-ph.GA astro-ph.CO

How does feedback affect the star formation histories of galaxies?

classification astro-ph.GA astro-ph.CO
keywords galaxy star formation historiesstellar feedbackAGN feedbackcosmological simulationsbaryon cyclingdouble power-law SFHCAMELSnormalizing flows
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

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

The reading

This paper tries to show that average galaxy star formation histories can be summarized by a double power-law curve whose five shape parameters depend on a small set of physical quantities: halo mass, baryon fraction, black hole mass, and feedback strength. Using the three CAMELS simulation models, a single set of equations with model-specific coefficients reproduces the average SFHs across all three. In the equations, cosmology sets the early rise, halo mass sets the overall scale, and feedback plus baryon cycling set the late-time decline and width. This matters because it offers a path from observations of galaxy samples back to feedback physics.

Core claim

The paper's central discovery is a system of empirical equations, Eqn 18, describing the average SFH shape parameters in terms of Omega_m, sigma_8, halo mass, baryon fraction, relative black hole mass, and the CAMELS feedback parameters. The same functional form applies to IllustrisTNG, SIMBA, and ASTRID, with only coefficients changing (Table 2). The rising slope beta depends mainly on Omega_m; the falling slope alpha on halo mass with AGN and baryon terms; the peak/width tau on cosmology, halo mass, baryon fraction, and black hole mass; the normalization phi on halo mass times feedback corrections; and the start time eta on Omega_m. The paper further finds that stellar feedback is the domi

What carries the argument

The load-bearing machinery is a double power-law SFH parameterization plus a normalizing-flow emulator. The double power law, SFH(t) = phi * ([(t-eta)/tau]^alpha + [(t-eta)/tau]^(-beta))^(-1), compresses each average SFH into five interpretable numbers. The normalizing flow, trained on the CAMELS LH simulations, generates average SFHs anywhere in parameter space, enabling regression of those five numbers against halo mass, baryon fraction, black hole mass, and feedback scalings. Symbolic regression term frequencies and generalized additive model losses guide which physical variables enter each equation, with a bias toward galaxy state variables over direct feedback parameters.

Load-bearing premise

The entire analysis uses SFHs sampled from a machine-learning emulator rather than directly from the simulations, and the emulator is validated only qualitatively against a small single-parameter set; if it is biased, or if the double power-law fails for a substantial fraction of average SFHs, the derived equations and conclusions inherit that distortion.

What would settle it

Recompute average SFHs directly from the CAMELS LH simulation particle data, fit double power laws, and compare the resulting alpha, beta, tau, phi, eta to Eqn 18 on a grid of parameter values; systematic deviations beyond sampling noise would disprove universality. Specifically examine ASTRID galaxies with late-time star-formation tails: if those SFHs are not well described by a double power law, the falling-slope equation is biased.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • Observed average SFHs could constrain cosmology and stellar feedback strength using Eqn 18 without rerunning simulations.
  • The same equation form should apply to other simulation codes, turning cross-model calibration into a coefficient-fitting exercise.
  • Because stellar feedback changes black hole growth, SFH-based constraints on stellar and AGN feedback will remain partially degenerate.
  • Reparameterizing winds by mass loading and energy per unit SFR makes the three models' SFH responses qualitatively consistent, clarifying interpretations.

Where Pith is reading between the lines

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

  • Eqn 18 is fitted on emulated SFHs; applying it to other simulation suites like SWIFT-EAGLE or CAMELS-SAM would test whether the shared functional form reflects physical regularity or a property of these three models.
  • AGN feedback is almost unconstrained from SFH shape alone, so combining these equations with baryon fraction or black hole mass observations should sharpen late-time feedback constraints.
  • The double power-law's poor fit to ASTRID's sustained late-time tails suggests a modified form with an added plateau could alter the predicted falling slopes and the conclusion that ASTRID galaxies quench fastest.
Share X Bluesky LinkedIn Reddit HN

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

4 major / 6 minor

Summary. The paper studies how variations in stellar feedback, AGN feedback, and cosmology affect the average star formation histories (SFHs) of galaxies at z~0 in the CAMELS IllustrisTNG, SIMBA, and ASTRID simulations. The authors train normalizing-flow emulators on the CAMELS Latin Hypercube (LH) sets, use the emulators to sample average SFHs, and fit these with a double power-law (Eq. 13) to obtain shape parameters {α, β, τ, φ, η}. They then use random forests, generalized additive models, and symbolic regression to construct a common set of equations (Eq. 18) with model-specific coefficients (Table 2) that relate SFH shape parameters to halo mass, baryon fraction, black hole-to-halo mass ratio, and CAMELS feedback/cosmology parameters. The paper also presents interaction analyses, a reparameterization of supernova feedback in terms of mass/energy loading, an SBI-based proof-of-concept inference of CAMELS parameters from SFHs, and several supplementary observational diagnostics.

Significance. If the central claim holds—that one functional form with model-dependent coefficients describes SFH shape across three different hydrodynamical codes—this would be a useful empirical framework for connecting SFH observations to feedback physics and for interpreting CAMELS parameter-space studies. The paper's strengths include its use of the public CAMELS suite, explicit treatment of galaxy selection and resolution limits, the combination of multiple ML methods, and the extensive diagnostic figures. However, the equations in Eq. 18 are fitted to emulator-sampled SFHs, the emulator validation is only qualitative, and the coefficients carry no reported uncertainties or out-of-sample tests. These gaps are load-bearing for the 'single set of equations' claim, so the result is currently a promising but not fully supported empirical summary.

major comments (4)
  1. [§3.1, Appendix A (Fig. 19)] All SFHs used to build Eqn. 18 are generated by sampling the trained normalizing flows ('Unless otherwise mentioned, all the SFHs in the following sections are generated by sampling the trained normalizing flows', §3.1). The only comparison to actual simulation outputs is the qualitative 1P validation in Fig. 19, which states 'qualitative agreement' but reports no quantitative coverage, calibration, or held-out error statistic. Since Eqn. 16 minimizes loss on emulator-sampled SFHs, any emulator bias in under-sampled regions of the 6D parameter space propagates directly into the DPL parameters and hence into every coefficient in Table 2. I request a quantitative validation of the emulator against held-out LH boxes and/or the 1P runs, with per-parameter residuals and coverage diagnostics, and a demonstration that the derived Eqn. 18 coefficients are stable when the emulator is retrained or
  2. [§6.1.1, Eq. (18), Table 2] The central equations are presented without any measure of predictive accuracy. The coefficients in Table 2 have no uncertainties, and the loss defined in Eq. (16) is minimized on the same emulator-generated data used to select the terms. No residual plots, R² values, or held-out predictions are shown for the DPL parameters. The paper even uses the GAM loss as a 'proxy of the Bayes risk' (§3.3.2), but does not compare Eq. 18 against that benchmark. Without out-of-sample validation, the claim that 'a single set of equations ... can describe the SFHs across all three CAMELS models' is supported only by in-sample agreement. Please provide bootstrap/subsample coefficient uncertainties and a held-out evaluation (e.g., fitting coefficients on a training subset and evaluating on a withheld LH subset, or predicting one model's coefficients from another).
  3. [§3.2, §4.2, §6.1.1] The paper acknowledges that the double power-law 'does not describe a subset of ASTRID SFHs with sustained late-time tails' (§3.2). These ASTRID SFHs are nevertheless included in the fits that determine the ASTRID coefficients in Table 2, where α and τ are the least constrained parameters. Because Eqn. 18 is claimed to hold across all three models, the fraction of average SFHs that are poorly represented by the DPL form must be quantified, and the sensitivity of the derived coefficients to excluding or reparameterizing these cases should be shown. If the DPL failure is non-negligible in ASTRID, the corresponding rows of Table 2 may encode an artifact of the fitting form rather than the feedback response.
  4. [§6.1.1, Eq. (16), Eq. (18)] The claim of a 'single set of equations' is weakened by the fact that the coefficients are free to vary per model in Eq. (16). If the functional form is the same but every coefficient differs, the statement reduces to 'each model can be fit by a member of a parametric family'. The paper needs to demonstrate what is shared beyond the functional form—for example, that the same terms remain important across models, that coefficients can be predicted from model properties, or that the equations generalize to held-out models/parameters. As written, the evidence for universality is largely the symbolic-regression term frequencies in Fig. 12, which are qualitative and in-sample.
minor comments (6)
  1. [§3.1, Appendix A] The validation text in Appendix A repeats 'qualitative agreement' twice. Please state explicitly which quantitative metrics were computed (e.g., coverage, calibration, chi-square) and whether any failed.
  2. [§3.2, Eq. (13)] The text defines α as the falling slope and β as the rising slope, but Figure 1 and Eq. (13) can be misread because the two power-law terms are symmetric. Consider adding a sentence explicitly defining the relation between α, β and the t<τ versus t>τ behavior.
  3. [§4.3, Eq. (17)] The χ² metric in Eq. (17) has unusual units (SFR² over SFR² integrated over time). Please clarify whether the integrand is intended to be a dimensionless ratio or whether the normalization is meant to produce a time-averaged statistic.
  4. [§3.3.1] The random forest feature importances in Fig. 13 are shown without error bars or sensitivity checks. Given that they are used as a 'sanity check' for Eq. 18, a bootstrap estimate would strengthen the comparison.
  5. [§6.1.1, Eq. (18)] The equation for η is written separately with no coefficients in Table 2. For completeness, state explicitly that η = 12.5Ωm − 3 is used for all three models, and whether a coefficient uncertainty was estimated.
  6. [General] There are several typographical and formatting issues, including 'early rimes' in §6.1.1, inconsistent use of 'Mhalo' (log Mhalo vs Mhalo) in Eq. (18) and Table 2, and missing figure cross-references in the text. A careful proofread is recommended.

Circularity Check

0 steps flagged

No circularity: Eqn 18 is an explicitly empirical fit, not a first-principles derivation; the real weaknesses are missing out-of-sample validation and qualitative emulator checks, not a definitional loop.

full rationale

The claimed 'derivation' of Eqn 18 is openly empirical: Section 3.3.3 leaves coefficients free and minimizes L = sum_model (Theta_SFH - sum_i c_i,model f(Theta_sfh,i))^2 (Eqn 16) on Theta_SFH parameters obtained by MCMC-fitting the same normalizing-flow-sampled SFHs, and Section 6.1.1 calls the result 'empirical equations.' Fitting a response surface to its own training data is in-sample regression, not circularity; the paper does not present Eqn 18 as a held-out prediction or as a consequence of a self-citation or uniqueness theorem. The normalizing flow itself is anchored to the LH simulations and has an admittedly qualitative 1P validation in Appendix A (Fig 19), and the acknowledged DPL failure for some ASTRID SFHs with late-time tails (Section 3.2) is a stated modeling limitation. The feedback-trend conclusions in Section 4 and the mass/energy-loading reparametrization in Section 6.2 are separate empirical findings. The genuine gap is the absence of a quantitative held-out test of Eqn 18 against direct simulation SFHs; that is a validation/correctness weakness, not a circular step. No load-bearing self-citation or ansatz-smuggling was found.

Axiom & Free-Parameter Ledger

1 free parameters · 5 axioms · 0 invented entities

The central equations rest on freshly fitted coefficients (Table 2), the emulator, the double power-law functional form, and the sufficiency of present-day state variables. No new physical entities are introduced. The free parameters are the equation coefficients; the axioms include the averaging procedure, the functional-form choice, emulator reliability, and the state-variable reduction.

free parameters (1)
  • SFH equation coefficients (Table 2) = Multiple values, e.g., tau_TNG: c1=4.571, c2=0.891, c3=-0.173, c4=0.930, c5=3.236
    Coefficients c1-c5 for each SFH shape parameter (tau, alpha, beta, phi, eta) are fitted per model by minimizing the joint loss across CAMELS models (Eqn 16). These are the principal fitted quantities in the paper.
axioms (5)
  • domain assumption The average SFH of 100 randomly selected galaxies in a box is a stable proxy for the mean SFH in that box and mass range.
    Introduced in Section 2.4; the entire analysis uses this averaged quantity rather than individual galaxies.
  • ad hoc to paper The double power-law form (Eqn 13) adequately captures the shape of average SFHs.
    Adopted in Section 3.2; the paper notes it does not describe a small fraction of ASTRID SFHs with sustained late-time tails.
  • domain assumption Simulation-based inference with normalizing flows trained on the LH datasets provides unbiased SFH samples across the parameter space.
    Used throughout Section 4 and in deriving Eqn 18; validation against 1P runs (Appendix A) is only qualitative.
  • ad hoc to paper The z=0 values of baryon fraction and black hole-to-halo mass ratio are sufficient state variables to characterize the integrated feedback history relevant to SFH shape.
    Used in Section 5 and in Eqn 18; these state variables are themselves outputs of the same simulations, so using them as inputs assumes they summarize the relevant physics without needing full time evolution.
  • domain assumption The CAMELS parameter variations (ASN1, ASN2, AAGN1, AAGN2, Omega_m, sigma_8) span the physically interesting feedback and cosmology space.
    Implicit in using the CAMELS LH sets; if the parameter priors miss relevant physics, the derived equations will not generalize.

reviewed 2026-08-05 · how reviews work

0 comments
Cite this review

Pith. "Pith review of How does feedback affect the star formation histories of galaxies?." pith.science (2026). https://pith.science/paper/MXCGDHLN

@misc{pith2026250821152,
  author       = {Pith},
  title        = {Pith review of: How does feedback affect the star formation histories of galaxies?},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/MXCGDHLN}},
  note         = {Machine review of arXiv:2508.21152}
}
Share X Bluesky LinkedIn Reddit HN
abstract

Star formation in galaxies is regulated by the interplay of a range of processes that shape the multiphase gas in the interstellar and circumgalactic media. Using the CAMELS suite of cosmological simulations, we study the effects of varying feedback and cosmology on the average star formation histories (SFHs) of galaxies at $z\sim0$ across the IllustrisTNG, SIMBA and ASTRID galaxy formation models. We find that galaxy SFHs in all three models are sensitive to changes in stellar feedback, which affects the efficiency of baryon cycling and the rates at which central black holes grow, while effects of varying AGN feedback depend on model-dependent implementations of black hole seeding, accretion and feedback. We also find strong interaction terms that couple stellar and AGN feedback, usually by regulating the amount of gas available for the central black hole to accrete. Using a double power-law to describe the average SFHs, we derive a general set of equations relating the shape of the SFHs to physical quantities like baryon fraction and black hole mass across all three models. We find that a single set of equations (albeit with different coefficients) can describe the SFHs across all three CAMELS models, with cosmology dominating the SFH at early times, followed by halo accretion, and feedback and baryon cycling at late times. Galaxy SFHs provide a novel, complementary probe to constrain cosmology and feedback, and can connect the observational constraints from current and upcoming galaxy surveys with the physical mechanisms responsible for regulating galaxy growth and quenching.

Figures

Figures reproduced from arXiv: 2508.21152 by Austen Gabrielpillai, Boon Kiat Oh, Camilla Pacifici, Christian Kragh Jespersen, Christopher C. Lovell, Daniel Angl\'es-Alc\'azar, Francisco Villaescusa-Navarro, Greg L. Bryan, John F. Wu, Joshua S. Speagle, Juan Pablo Alfonzo, Kartheik G. Iyer, Lars Hernquist, Laura Sommovigo, Lucia A. Perez, Megan T. Tillman, Rachel S. Somerville, Romeel Dav\'e, Sandro Tacchella, Shy Genel, Suchetha Cooray, Sultan Hassan, Tjitske K. Starkenburg.

Figure 1
Figure 1. Figure 1: Motivation: The star formation histories (SFHs) of galaxies act as tracers of the overall gas regulation in the galaxy, which depend sensitively on the adopted prescription and strength of stellar and AGN feedback. By analyzing how the shape of the average SFH at a given mass (described using normalization (ϕ), peak/width (τ ), and rising and falling slopes (α, β)) depends on the strength and nature of the… view at source ↗
Figure 2
Figure 2. Figure 2: Median SFHs in different stellar mass bins from the three CAMELS models we use in this work, compiled using the CV datasets. The plot shows the SFHs of galaxies starting at 109.25M⊙ in increments of 0.2 dex up to 1010.85M⊙, with darker colors indicating higher mass bins. Feedback type TNG SIMBA ASTRID ASN1 Energy per unit SFR Mass loading factor Energy per unit SFR (Galactic Winds; range: (0.25,4)) (Galact… view at source ↗
Figure 3
Figure 3. Figure 3: Changes in the average SFHs of galaxies with variations in the double power-law parameters τ (peak/width), α (falling slope), β (rising slope), and ϕ (overall normalization), as described in Eqn. 13. plain them needs to account for both the effects of in￾dividual parameter variations as well as interactions. In this section, we briefly describe the three methods we use to examine the relation between the a… view at source ↗
Figure 4
Figure 4. Figure 4: Variation in the average SFHs of 100 randomly selected galaxies with stellar masses 109.5 <M∗ < 1011.5M⊙ caused by varying the cosmology and feedback strengths across the three CAMELS models (ASTRID, TNG, SIMBA), predicted using the trained normalizing flows. Solid lines and shaded regions show the median SFH and corresponding 1σ uncertainties from the bootstrap sampling used to train the normalizing flow.… view at source ↗
Figure 5
Figure 5. Figure 5: Similar to [PITH_FULL_IMAGE:figures/full_fig_p015_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: Similar to [PITH_FULL_IMAGE:figures/full_fig_p016_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: Similar to [PITH_FULL_IMAGE:figures/full_fig_p017_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: Interactions between different feedback and cosmology channels (thicker lines) can sometimes affect the SFHs beyond the linear additive effects of variations in the individual parameters (thin dashed lines). The three examples illustrate the effects of the interactions between the galactic wind parameters (ASN1 and ASN2) in TNG (top), between cosmology and feedback (Ωm and ASN1) in SIMBA (middle) and betwe… view at source ↗
Figure 9
Figure 9. Figure 9: Baryon content as a function of halo mass across the three CAMELS models [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: Similar to [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Similar to [PITH_FULL_IMAGE:figures/full_fig_p024_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: The frequency with which different terms occur in the symbolic regression solutions to the SFH parameters across the three CAMELS models. While different from traditional metrics of feature importance (e.g. using saliency maps or random forests), the occurrence of certain recurring terms (e.g. Mhalo being an important feature for ϕ across all three models) in the set of equations generated by symbolic reg… view at source ↗
Figure 13
Figure 13. Figure 13: (Top) Summary of the various factors impacting average galaxy SFHs across the three CAMELS models. (Bottom) Feature importance of the cosmology, feedback scaling factors and galaxy properties as defined in Eqn. (18) on the star formation rate as a function of lookback time. Cosmology dominates the average SFR of galaxies at early times, giving way to halo mass as they assemble. At late times, the SFR is a… view at source ↗
Figure 14
Figure 14. Figure 14: Changes in the SFHs for all three CAMELS models while varying ASN1 and ASN2 such that the energy per unit SFR is varied while keeping the mass loading factor fixed, at low (top), intermediate (middle) and high (bottom row) halo masses (using the normalizing flow emulator developed in Section 3.1). Note that (except for the high mass halos where BH feedback dominates), increasing the energy per units SFR (… view at source ↗
Figure 15
Figure 15. Figure 15: Same as [PITH_FULL_IMAGE:figures/full_fig_p032_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Variations in the baryon fraction (top panels) and black hole mass (bottom panels) for all three simulations (going from left to right: CAMELS/TNG, CAMELS/SIMBA, and CAMELS/ASTRID) as a function of halo mass. Within each set of six panels, the top column shows the result of varying the mass loading while keeping the energy per SFR constant while the bottom column is the opposite. time star formation. This… view at source ↗
Figure 17
Figure 17. Figure 17: Varying the energy per unit SFR (ASN1) and wind speed (ASN2) in CAMELS/TNG have distinct effects on the baryon fraction and stellar-to-gas mass ratio (SGMR). This allows us to disentangle the effects of factors that change the gas content of halos from factors that change the efficiency of forming stars from that gas. The black lines show lines of constant baryon-fraction to stellar-to-gas mass ratio. Whi… view at source ↗
Figure 18
Figure 18. Figure 18: Inferring the CAMELS box parameters (both cosmology and feedback strengths) from a sample of 100 SFHs using SBI across the three models. Top: Corner plots showing the inferred posteriors for a single sample of SFHs corresponding to a realization of cosmology and feedback strength. Bottom: Average results across 1000 realizations across the full CAMELS parameter space. The cosmology is well constrained, as… view at source ↗
Figure 19
Figure 19. Figure 19: Validation plots comparing median (solid lines) and 16-84th percentiles (shaded regions) of the SFHs sampled from the trained normalizing flows (labeled NN) compared against SFHs from the CAMELS 1P dataset that vary a single parameter in the CAMELS box while keeping the others fixed at the fiducial value. The different panels show the 1P dataset for all six CAMELS box parameters. While the 1P datasets sho… view at source ↗
Figure 20
Figure 20. Figure 20: The distribution of χ 2 SFH values across the CAMELS parameter space when comparing the SFHs varying 2 parameters to the linear combination of the corresponding 1 parameter variations for the three CAMELS models: TNG (blue), SIMBA (red) and ASTRID (green). Values in each cell show the χ 2 SFH value, which also correspond the intensity of the colormaps [PITH_FULL_IMAGE:figures/full_fig_p042_20.png] view at source ↗
Figure 21
Figure 21. Figure 21: Similar to [PITH_FULL_IMAGE:figures/full_fig_p042_21.png] view at source ↗
Figure 22
Figure 22. Figure 22: The impact of varying cosmology and feedback strengths on the gas phase metallicity of galaxies as a function of halo mass at z ∼ 0 in the three models. Astropy Collaboration, Price-Whelan, A. M., Lim, P. L., et al. 2022, ApJ, 935, 167, doi: 10.3847/1538-4357/ac7c74 Ayromlou, M., Nelson, D., & Pillepich, A. 2023, MNRAS, 524, 5391, doi: 10.1093/mnras/stad2046 Bagley, M. B., Pirzkal, N., Finkelstein, S. L.,… view at source ↗
Figure 23
Figure 23. Figure 23: Similar to [PITH_FULL_IMAGE:figures/full_fig_p044_23.png] view at source ↗
Figure 24
Figure 24. Figure 24: Similar to [PITH_FULL_IMAGE:figures/full_fig_p044_24.png] view at source ↗

discussion (0)

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

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

  1. POLAR-II: modeling star formation history of galaxies on the 21-cm signal from Epoch of Reionization

    astro-ph.GA 2026-03 unverdicted novelty 4.0

    Varying the star formation history of galaxies while keeping total ionizing photons fixed changes the size and warmth of ionized regions, altering the topology of reionization and the 21-cm global signal and power spectrum.

  2. Quenching of Star Formation in Massive Galaxies

    astro-ph.GA 2026-06 unverdicted novelty 2.0

    Massive galaxies quench via either rapid supermassive black hole outflow-driven shutdown or gradual gas exhaustion, virial heating, or preventative feedback, each with distinct signatures.

Reference graph

Works this paper leans on

169 extracted references · 14 canonical work pages · cited by 2 Pith papers · 6 internal anchors

  1. [1]

    , " * write output.state after.block = add.period write newline

    ENTRY address archivePrefix author booktitle chapter doi edition editor eprint howpublished institution journal key month number organization pages publisher school series title misctitle type volume year version url label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts ...

  2. [2]

    write newline

    " write newline "" before.all 'output.state := FUNCTION format.url url empty "" new.block "" url * "" * if FUNCTION format.eprint eprint empty "" archivePrefix empty "" archivePrefix "arXiv" = new.block " " eprint * " " * new.block " " eprint * " " * if if if FUNCTION format.doi doi empty "" " " doi * " " * if FUNCTION format.pid doi empty eprint empty ur...

  3. [3]

    " nf| N3 t1 Q1 Q1 Q1 T|k_׿I ǸBDD p 7Tm߾ ; [h:x [._Y2X C| G ;C!

    thebibliography [1] 20pt to REFERENCES 6pt =0pt -12pt 10pt plus 3pt =0pt =0pt =1pt plus 1pt =0pt =0pt -12pt =13pt plus 1pt =20pt =13pt plus 1pt \@M =10000 =-1.0em =0pt =0pt 0pt =0pt =1.0em @enumiv\@empty 10000 10000 `\.\@m \@noitemerr \@latex@warning Empty `thebibliography' environment \@ifnextchar \@reference \@latexerr Missing key on reference command E...

  4. [4]

    E., Gladders , M

    Abramson , L. E., Gladders , M. D., Dressler , A., et al. 2016, , 832, 7, 10.3847/0004-637X/832/1/7

  5. [5]

    Angl \'e s-Alc \'a zar , D., Dav \'e , R., Faucher-Gigu \`e re , C.-A., \"O zel , F., & Hopkins , P. F. 2017 a , , 464, 2840, 10.1093/mnras/stw2565

  6. [6]

    2017 b , , 470, 4698, 10.1093/mnras/stx1517

    Angl \'e s-Alc \'a zar , D., Faucher-Gigu \`e re , C.-A., Kere s , D., et al. 2017 b , , 470, 4698, 10.1093/mnras/stx1517

  7. [7]

    2017 c , , 472, L109, 10.1093/mnrasl/slx161

    Angl \'e s-Alc \'a zar , D., Faucher-Gigu \`e re , C.-A., Quataert , E., et al. 2017 c , , 472, L109, 10.1093/mnrasl/slx161

  8. [8]

    F., et al

    Angl \'e s-Alc \'a zar , D., Quataert , E., Hopkins , P. F., et al. 2021, , 917, 53, 10.3847/1538-4357/ac09e8

  9. [9]

    P., Tollerud , E

    Astropy Collaboration , Robitaille , T. P., Tollerud , E. J., et al. 2013, , 558, A33, 10.1051/0004-6361/201322068

  10. [10]

    M., Sip o cz , B

    Astropy Collaboration , Price-Whelan , A. M., Sip o cz , B. M., et al. 2018, , 156, 123, 10.3847/1538-3881/aabc4f

  11. [11]

    M., Lim , P

    Astropy Collaboration , Price-Whelan , A. M., Lim , P. L., et al. 2022, , 935, 167, 10.3847/1538-4357/ac7c74

  12. [12]

    2023, , 524, 5391, 10.1093/mnras/stad2046

    Ayromlou , M., Nelson , D., & Pillepich , A. 2023, , 524, 5391, 10.1093/mnras/stad2046

  13. [13]

    B., Pirzkal , N., Finkelstein , S

    Bagley , M. B., Pirzkal , N., Finkelstein , S. L., et al. 2024, , 965, L6, 10.3847/2041-8213/ad2f31

  14. [14]

    S., Devriendt , J., Slyz , A., et al

    Beckmann , R. S., Devriendt , J., Slyz , A., et al. 2017, , 472, 949, 10.1093/mnras/stx1831

  15. [15]

    S., Smith , M

    Bennett , J. S., Smith , M. C., Fielding , D. B., et al. 2024, arXiv e-prints, arXiv:2410.12909, 10.48550/arXiv.2410.12909

  16. [16]

    2022, , 512, 3703, 10.1093/mnras/stac648

    Bird , S., Ni , Y., Di Matteo , T., et al. 2022, , 512, 3703, 10.1093/mnras/stac648

  17. [17]

    2014, , 445, 2313, 10.1093/mnras/stu1923

    Bird , S., Vogelsberger , M., Haehnelt , M., et al. 2014, , 445, 2313, 10.1093/mnras/stu1923

  18. [18]

    Bluck , A. F. L., Maiolino , R., Brownson , S., et al. 2022, , 659, A160, 10.1051/0004-6361/202142643

  19. [19]

    Bluck , A. F. L., Piotrowska , J. M., & Maiolino , R. 2023, , 944, 108, 10.3847/1538-4357/acac7c

  20. [20]

    1952, , 112, 195, 10.1093/mnras/112.2.195

    Bondi , H. 1952, , 112, 195, 10.1093/mnras/112.2.195

  21. [21]

    M., et al

    Borrow , J., Schaller , M., Bah \'e , Y. M., et al. 2023, , 526, 2441, 10.1093/mnras/stad2928

  22. [22]

    2010, , 718, 1001, 10.1088/0004-637X/718/2/1001

    Bouch \'e , N., Dekel , A., Genzel , R., et al. 2010, , 718, 1001, 10.1088/0004-637X/718/2/1001

  23. [23]

    G., Schaye , J., Frenk , C

    Bower , R. G., Schaye , J., Frenk , C. S., et al. 2017, , 465, 32, 10.1093/mnras/stw2735

  24. [24]

    2023, arXiv e-prints, arXiv:2312.08277, 10.48550/arXiv.2312.08277

    Braspenning , J., Schaye , J., Schaller , M., et al. 2023, arXiv e-prints, arXiv:2312.08277, 10.48550/arXiv.2312.08277

  25. [25]

    S., et al

    Brennan , R., Choi , E., Somerville , R. S., et al. 2018, , 860, 14, 10.3847/1538-4357/aac2c4

  26. [26]

    2016, in Astrophysics and Space Science Library, Vol

    Brooks , A., & Christensen , C. 2016, in Astrophysics and Space Science Library, Vol. 418, Galactic Bulges, ed. E. Laurikainen , R. Peletier , & D. Gadotti , 317, 10.1007/978-3-319-19378-6_12

  27. [27]

    ARTEMIS emulator: exploring the effect of cosmology and galaxy formation physics on Milky Way-mass haloes and their satellites

    Brown , S. T., Fattahi , A., McCarthy , I. G., et al. 2024, arXiv e-prints, arXiv:2403.11692, 10.48550/arXiv.2403.11692

  28. [28]

    J., Owers , M

    Bryant , J. J., Owers , M. S., Robotham , A. S. G., et al. 2015, , 447, 2857, 10.1093/mnras/stu2635

  29. [29]

    A., Law , D

    Bundy , K., Bershady , M. A., Law , D. R., et al. 2015, , 798, 7, 10.1088/0004-637X/798/1/7

  30. [30]

    R., et al

    Busillo , V., Tortora , C., Napolitano , N. R., et al. 2023, , 525, 6191, 10.1093/mnras/stad2691

  31. [31]

    2023, , 520, 722, 10.1093/mnras/stad171

    Byrne , L., Faucher-Gigu \`e re , C.-A., Stern , J., et al. 2023, , 520, 722, 10.1093/mnras/stad171

  32. [32]

    C., Leja , J., Johnson , B

    Carnall , A. C., Leja , J., Johnson , B. D., et al. 2019, , 873, 44, 10.3847/1538-4357/ab04a2

  33. [33]

    L., Fielding , D

    Carr , C., Bryan , G. L., Fielding , D. B., Pandya , V., & Somerville , R. S. 2023, , 949, 21, 10.3847/1538-4357/acc4c7

  34. [34]

    M., Kartaltepe , J

    Casey , C. M., Kartaltepe , J. S., Drakos , N. E., et al. 2023, , 954, 31, 10.3847/1538-4357/acc2bc

  35. [35]

    C., Magnier , E

    Chambers , K. C., Magnier , E. A., Metcalfe , N., et al. 2016, arXiv e-prints, arXiv:1612.05560, 10.48550/arXiv.1612.05560

  36. [36]

    M., Koo , D

    Chen , Z., Faber , S. M., Koo , D. C., et al. 2020, , 897, 102, 10.3847/1538-4357/ab9633

  37. [37]

    P., Naab , T., Oser , L., & Moster , B

    Choi , E., Ostriker , J. P., Naab , T., Oser , L., & Moster , B. P. 2015, , 449, 4105, 10.1093/mnras/stv575

  38. [38]

    S., Ostriker , J

    Choi , E., Somerville , R. S., Ostriker , J. P., Hirschmann , M., & Naab , T. 2024, , 964, 54, 10.3847/1538-4357/ad245a

  39. [39]

    S., Ostriker , J

    Choi , E., Somerville , R. S., Ostriker , J. P., Naab , T., & Hirschmann , M. 2018, , 866, 91, 10.3847/1538-4357/aae076

  40. [40]

    A., Wyithe , J

    Correa , C. A., Wyithe , J. S. B., Schaye , J., & Duffy , A. R. 2015, , 450, 1514, 10.1093/mnras/stv689

  41. [41]

    A., Schaye , J., Bower , R

    Crain , R. A., Schaye , J., Bower , R. G., et al. 2015, , 450, 1937, 10.1093/mnras/stv725

  42. [42]

    2024, PySR: High-Performance Symbolic Regression in Python and Julia , Astrophysics Source Code Library, record ascl:2409.018

    Cranmer , M. 2024, PySR: High-Performance Symbolic Regression in Python and Julia , Astrophysics Source Code Library, record ascl:2409.018

  43. [43]

    2019, , 486, 2827, 10.1093/mnras/stz937

    Dav \'e , R., Angl \'e s-Alc \'a zar , D., Narayanan , D., et al. 2019, , 486, 2827, 10.1093/mnras/stz937

  44. [45]

    Dav \'e , R., Thompson , R., & Hopkins , P. F. 2016, , 462, 3265, 10.1093/mnras/stw1862

  45. [46]

    de Santi , N. S. M., Shao , H., Villaescusa-Navarro , F., et al. 2023, , 952, 69, 10.3847/1538-4357/acd1e2

  46. [47]

    2014, , 444, 2071, 10.1093/mnras/stu1427

    Dekel , A., & Mandelker , N. 2014, , 444, 2071, 10.1093/mnras/stu1427

  47. [48]

    M., Angl \'e s-Alc \'a zar , D., Thiele , L., et al

    Delgado , A. M., Angl \'e s-Alc \'a zar , D., Thiele , L., et al. 2023, , 526, 5306, 10.1093/mnras/stad2992

  48. [49]

    2016, arXiv e-prints, arXiv:1611.00036, 10.48550/arXiv.1611.00036

    DESI Collaboration , Aghamousa , A., Aguilar , J., et al. 2016, arXiv e-prints, arXiv:1611.00036, 10.48550/arXiv.1611.00036

  49. [50]

    B., Macci \`o , A

    Di Cintio , A., Brook , C. B., Macci \`o , A. V., et al. 2014, , 437, 415, 10.1093/mnras/stt1891

  50. [51]

    2008, , 676, 33, 10.1086/524921

    Di Matteo , T., Colberg , J., Springel , V., Hernquist , L., & Sijacki , D. 2008, , 676, 33, 10.1086/524921

  51. [52]

    M., Starkenburg , T

    Dickey , C. M., Starkenburg , T. K., Geha , M., et al. 2021, , 915, 53, 10.3847/1538-4357/abc014

  52. [53]

    E., & Torrey , P

    Diemer , B., Sparre , M., Abramson , L. E., & Torrey , P. 2017, , 839, 26, 10.3847/1538-4357/aa68e5

  53. [54]

    T., McLure , R

    Donnan , C. T., McLure , R. J., Dunlop , J. S., et al. 2024, , 533, 3222, 10.1093/mnras/stae2037

  54. [56]

    A., Macci \`o , A

    Dutton , A. A., Macci \`o , A. V., Buck , T., et al. 2019, , 486, 655, 10.1093/mnras/stz889

  55. [57]

    J., Willott , C., Alberts , S., et al

    Eisenstein , D. J., Willott , C., Alberts , S., et al. 2023, arXiv e-prints, arXiv:2306.02465, 10.48550/arXiv.2306.02465

  56. [58]

    S., Jenkins , A., et al

    Elbers , W., Frenk , C. S., Jenkins , A., et al. 2024, arXiv e-prints, arXiv:2403.12967, 10.48550/arXiv.2403.12967

  57. [59]

    Fabian , A. C. 2012, , 50, 455, 10.1146/annurev-astro-081811-125521

  58. [60]

    L., Bagley , M

    Finkelstein , S. L., Bagley , M. B., Ferguson , H. C., et al. 2023, , 946, L13, 10.3847/2041-8213/acade4

  59. [61]

    L., Bagley , M

    Finkelstein , S. L., Bagley , M. B., Arrabal Haro , P., et al. 2025, , 983, L4, 10.3847/2041-8213/adbbd3

  60. [62]

    2017, , 601, A143, 10.1051/0004-6361/201629478

    Fiore , F., Feruglio , C., Shankar , F., et al. 2017, , 601, A143, 10.1051/0004-6361/201629478

  61. [63]

    R., & Mirocha , J

    Furlanetto , S. R., & Mirocha , J. 2022, , 511, 3895, 10.1093/mnras/stac310

  62. [64]

    2024, , 529, 4896, 10.1093/mnras/stae817

    Gebhardt , M., Angl \'e s-Alc \'a zar , D., Borrow , J., et al. 2024, , 529, 4896, 10.1093/mnras/stae817

  63. [65]

    2023, , 954, L11, 10.3847/2041-8213/acee80

    Gelli , V., Salvadori , S., Ferrara , A., Pallottini , A., & Carniani , S. 2023, , 954, L11, 10.3847/2041-8213/acee80

  64. [66]

    2019, in International Conference on Machine Learning, PMLR, 2404--2414

    Greenberg, D., Nonnenmacher, M., & Macke, J. 2019, in International Conference on Machine Learning, PMLR, 2404--2414

  65. [67]

    A., Kocevski , D

    Grogin , N. A., Kocevski , D. D., Faber , S. M., et al. 2011, , 197, 35, 10.1088/0067-0049/197/2/35

  66. [68]

    2017, , 468, 3935, 10.1093/mnras/stx666

    Habouzit , M., Volonteri , M., & Dubois , Y. 2017, , 468, 3935, 10.1093/mnras/stx666

  67. [69]

    Harris el al., C. R. 2020, Nat, 585, 357, https://doi.org/10.1038/s41586-020-2649-2

  68. [70]

    J., & Tibshirani, R

    Hastie, T. J., & Tibshirani, R. J. 1990, Generalized additive models, Vol. 43 (CRC press)

  69. [71]

    C., & Hopkins , P

    Hayward , C. C., & Hopkins , P. F. 2017, , 465, 1682, 10.1093/mnras/stw2888

  70. [72]

    M., & Best , P

    Heckman , T. M., & Best , P. N. 2023, Galaxies, 11, 21, 10.3390/galaxies11010021

  71. [73]

    J., Chartier , N., et al

    Ho , M., Bartlett , D. J., Chartier , N., et al. 2024, arXiv e-prints, arXiv:2402.05137, 10.48550/arXiv.2402.05137

  72. [74]

    F., Cox , T

    Hopkins , P. F., Cox , T. J., Kere s , D., & Hernquist , L. 2008, , 175, 390, 10.1086/524363

  73. [75]

    F., Kere s , D., O \ n orbe , J., et al

    Hopkins , P. F., Kere s , D., O \ n orbe , J., et al. 2014, , 445, 581, 10.1093/mnras/stu1738

  74. [76]

    F., Wetzel , A., Kere s , D., et al

    Hopkins , P. F., Wetzel , A., Kere s , D., et al. 2018, , 480, 800, 10.1093/mnras/sty1690

  75. [77]

    F., Grudic , M

    Hopkins , P. F., Grudic , M. Y., Su , K.-Y., et al. 2024, The Open Journal of Astrophysics, 7, 18, 10.21105/astro.2309.13115

  76. [78]

    2021, , 502, 6010, 10.1093/mnras/stab357

    Huang , H.-J., Eifler , T., Mandelbaum , R., et al. 2021, , 502, 6010, 10.1093/mnras/stab357

  77. [79]

    Hunter, J. D. 2007, Computing in Science & Engineering, 9, 90, 10.1109/MCSE.2007.55

  78. [80]

    G., Gawiser , E., Faber , S

    Iyer , K. G., Gawiser , E., Faber , S. M., et al. 2019, , 879, 116, 10.3847/1538-4357/ab2052

  79. [81]

    G., Speagle , J

    Iyer , K. G., Speagle , J. S., Caplar , N., et al. 2024 a , , 961, 53, 10.3847/1538-4357/acff64

  80. [82]

    G., Tacchella , S., Genel , S., et al

    Iyer , K. G., Tacchella , S., Genel , S., et al. 2020, , 498, 430, 10.1093/mnras/staa2150

Showing first 80 references.

This paper was first reviewed by deepseek-v4-flash on August 5, 2026.