pith. sign in

arxiv: 2606.10026 · v1 · pith:KEKQIDE5new · submitted 2026-06-08 · 🌌 astro-ph.GA · astro-ph.CO

Pressure-regulated feedback-modulated star formation as a subgrid model for galaxy formation simulations

Pith reviewed 2026-06-27 15:48 UTC · model grok-4.3

classification 🌌 astro-ph.GA astro-ph.CO
keywords star formationsubgrid modelgalaxy formation simulationsinterstellar mediumstellar feedbackgas depletion time
0
0 comments X

The pith

PRFM subgrid model ties star formation rate to the local balance between gravity and stellar feedback in galaxy disks.

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

The paper introduces two versions of a subgrid prescription for star formation in cosmological galaxy simulations, both derived from the pressure-regulated feedback-modulated theory. Instead of using fixed efficiency parameters, the model sets the star formation rate from the equilibrium between gravitational collapse and the momentum injected by stellar feedback. Calibrated directly on parsec-scale TIGRESS simulations, the prescriptions are tested in isolated Milky Way-like disks at resolutions from 10^5 to 10^7 solar masses. Both the volumetric and integrated implementations produce shorter gas depletion times than the IllustrisTNG approach, especially at high pressure and density, while the integrated version remains stable when the disk scale height is unresolved.

Core claim

The PRFM theory supplies a subgrid star formation model in which the local star formation rate is determined by the dynamic balance between gravity and the momentum supplied by stellar feedback; the two numerical implementations (volumetric when scale height is resolved, integrated when it is not) reproduce the star formation rates of resolved TIGRESS simulations more closely than empirical prescriptions and remain consistent across a wide range of mass resolutions.

What carries the argument

Pressure-regulated, feedback-modulated (PRFM) theory, which equates the star formation rate to the rate needed to maintain vertical equilibrium between gravitational weight and feedback momentum injection.

If this is right

  • Gas depletion times become shorter than in current empirical models, particularly inside dense, high-pressure regions.
  • At the highest tested resolution the two PRFM implementations agree with each other and with the underlying TIGRESS simulation for the global star formation rate.
  • The integrated formulation remains numerically stable and produces consistent results even when the gas disk scale height is not resolved.
  • Both versions supply a physically motivated effective equation of state that can be used directly in larger-volume cosmological simulations.

Where Pith is reading between the lines

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

  • Adoption of PRFM would likely reduce the need for separate tunable parameters that currently control star formation efficiency in cosmological codes.
  • Because the model predicts an explicit dependence of star formation on local pressure, it could alter the predicted distribution of star-forming gas across different galactic environments.
  • The robustness of the integrated version at low resolution suggests it could be applied immediately to existing large-volume runs without requiring changes to the numerical resolution.

Load-bearing premise

The calibration of PRFM on parsec-scale simulations can be applied unchanged as a subgrid rule at the much coarser resolutions used in cosmological galaxy formation runs.

What would settle it

A direct comparison, at the same mass resolution, of the star formation rate surface density versus mid-plane pressure relation in a full cosmological run using the PRFM subgrid model versus one using the IllustrisTNG prescription, checked against observed Kennicutt-Schmidt relations in nearby galaxies.

Figures

Figures reproduced from arXiv: 2606.10026 by Chang-Goo Kim, Eve C. Ostriker, Jan Burger, Sarah M. R. Jeffreson.

Figure 1
Figure 1. Figure 1: Equation of state comparison. Distribution of the gas effective pressure Peff (left and right) or weight W (center) as a function of the measured gas hydrogen number density nH,meas (left and right) or computed equilibrium density nH (center), in pixels of scale 1 × 1 kpc2 , and averaged over all simulation time-stamps with an interval of 10 Myr between 50 Myr and 600 Myr, inclusive. All results are from s… view at source ↗
Figure 2
Figure 2. Figure 2: Upper row: Effective gas velocity dispersion (σeff , orange, computed from the effective EoS), and measured stellar velocity dispersion (σ∗,z, blue) in the vertical direction as a function of galactocentric radius R from simulations adopting the PRFM-int implementation, at resolutions of 105 , 106 and 107 M⊙. Lower row: Ratio of these two velocity dispersions (blue lines), in comparison to the measured rat… view at source ↗
Figure 3
Figure 3. Figure 3: Ratio of dynamical time to depletion time, i.e. the star formation efficiency per dynamical time εdyn = tdyn/tdep (see Equation 1), as a function of gas pressure or weight. We show the distribution of gas cells in the plane spanned by the measured mid-plane effective pressure Peff (left and right columns) or computed weight W (center column) and the ratio of the dynamical time to depletion time tdyn/tdep. … view at source ↗
Figure 4
Figure 4. Figure 4: Depletion time as a function of gas density, from simulations with mass resolution 105 M⊙, at a simulation time of 50 Myr. The gas cell depletion time is tdep = Mg/M˙ ∗, and density is measured in the simulation as nH,meas (left and right) or computed as an equilibrium value nH (center). The black lines represent the least-squares best fit (log tdep ≈ −0.96 log nH −0.46) to the PRFM-vol simulation, includi… view at source ↗
Figure 5
Figure 5. Figure 5: Depletion and dynamical time as a function of gas pressure, from simulations with mass resolution 105 M⊙ at time 50 Myr. We show the distribution of gas cells in the plane spanned by the gas cell measured effective pressure Peff (left and right columns) or computed weight W (center column) and the gas cell depletion time tdep ≡ Mg/M˙ ∗ (top row) or the gas cell dynamical time tdyn (bottom row). For PRFM-vo… view at source ↗
Figure 6
Figure 6. Figure 6: Time evolution of the total galactic SFR in each set of simulations (PRFM-vol, PRFM-int, and TNG-model implementations) at each of the three resolutions tested (105 M⊙ in blue, 106 M⊙ in orange, and 107 M⊙ in green). The transparent blue line in the center panel is a copy of the 105 M⊙ evolution from the PRFM-vol case (blue line, left panel). SFRs obtained at high pressure (and density), relative to the im… view at source ↗
Figure 7
Figure 7. Figure 7: Distribution of SFR surface densities ΣSFR as a function of the mid-plane effective gas pressure Peff (left and right), or computed weight W (center), measured in projected pixels of scale 1 × 1 kpc2 , and averaged over all simulation time-stamps between 50 Myr and 600 Myr, inclusive. All three simulation resolutions of 105 M⊙ (top row), 106 M⊙ (center row) and 107 M⊙ (bottom row) are shown. The thick blac… view at source ↗
Figure 8
Figure 8. Figure 8: Distribution of gas cells in the plane spanned by the gas cell dynamical time tdyn and the gas cell depletion time tdep = Mg/M˙ ∗, at a simulation time of 50 Myr, at the three simulation resolutions of 105 M⊙ (top row), 106 M⊙ (center row) and 107 M⊙ (bottom row). Computation of timescales is described in [PITH_FULL_IMAGE:figures/full_fig_p018_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Maps of the gas surface density Σg (upper row), the stellar surface density Σ∗ (center row) and the SFR surface density ΣSFR (lower row) for simulations adopting the PRFM-vol implementation. Results are shown from simulations at three different mass resolutions (and corresponding typical cell size): 105 M⊙ (80 pc), 106 M⊙ (190 pc), 107 M⊙ (370 pc), at a simulation time of 50 Myr. In the top two rows, the m… view at source ↗
Figure 10
Figure 10. Figure 10: Upper panels: Histograms of the gas surface density Σg and the SFR surface density ΣSFR, summed over columns at a scale of 1 kpc. Results are from simulations adopting the PRFM-vol implementation at resolutions of 105 M⊙ (blue), 106 M⊙ (orange) and 107 M⊙ (green), at a time of 50 Myr. Lower panels: Histograms of the depletion time tdep computed from individual cells within the PRFM model for each simulati… view at source ↗
Figure 11
Figure 11. Figure 11: Histograms of the measured pressure (left), measured volume density (center), and measured scale-height (right) of star-forming gas from simulations adopting the PRFM-vol implementation at resolutions of 105 M⊙ (blue), 106 M⊙ (orange) and 107 M⊙ (green). All values are those obtained directly from the simulation in post-processing. −1 0 log (Peff/W) 0 1 2 3 PDF 105M 106M 107M −2 −1 0 1 log (nH,meas/nH) −1… view at source ↗
Figure 12
Figure 12. Figure 12: Comparison of measured values for ISM properties and equilibrium predictions from simulations using the PRFM￾vol implementation at a resolution of 105 M⊙ (blue), 106 M⊙ (orange) and 107 M⊙ (green). Left: Ratios of the effective pressure Peff at the midplane, computed via the EoS (Equation 15) from the measured gas density, to the dynamical equilibrium pressure computed via the ISM weight W (using Equation… view at source ↗
Figure 13
Figure 13. Figure 13: Maps of Σg, Σ∗, and ΣSFR as in [PITH_FULL_IMAGE:figures/full_fig_p024_13.png] view at source ↗
Figure 14
Figure 14. Figure 14: Histograms of Σg, ΣSFR, and tdep as in [PITH_FULL_IMAGE:figures/full_fig_p025_14.png] view at source ↗
Figure 15
Figure 15. Figure 15: Histograms of the computed ISM weight W from Equation 31 (left), computed equilibrium mid-plane volume density nH from Equation 29 (center), and computed equilibrium scale-height Hg from Equation 28 (right) for the PRFM-int implementation at resolutions of 105 M⊙ (blue), 106 M⊙ (orange) and 107 M⊙ (green) [PITH_FULL_IMAGE:figures/full_fig_p025_15.png] view at source ↗
Figure 16
Figure 16. Figure 16: Radial profiles of the gas surface densities Σg (top) and volume densities n (bottom) in the PRFM-vol (left) and PRFM-int (right) simulations at all numerical resolutions (see key). In the lower panels, solid lines show measured volume densities nH,meas, while the dotted lines in the lower right show calculated volume densities nH. The dashed horizontal black lines in lower panels indicate the star format… view at source ↗
read the original abstract

We present a new subgrid model for interstellar gas evolution in cosmological simulations of galaxy formation, based on the pressure-regulated, feedback-modulated (PRFM) theory of star formation. In contrast to the empirically pegged star formation prescriptions employed in current cosmological simulations, the PRFM model links the local star formation rate to the dynamic balance achieved in galactic interstellar gas between gravity and stellar feedback effects. With this formulation, both the star formation efficiency and the effective equation of state may be directly calibrated using numerical simulations, such as TIGRESS, which resolve physics of the interstellar medium and star formation at parsec scales. We develop, and implement in the Arepo moving-mesh code, two complementary classes of the subgrid model: a volumetric version (PRFM-vol) applicable when the gas disk scale height of a galaxy is numerically resolved in a simulation, and an integrated version (PRFM-int) that reconstructs the mid-plane density and pressure from vertical equilibrium considerations when the true gas scale height cannot be numerically resolved. Using isolated Milky-Way-like disk simulations across mass resolutions $10^5$-$10^7~{\rm M}_\odot$, we show that both implementations yield shorter gas depletion times than the IllustrisTNG prescription, especially in regions where pressure and density are large. At high resolution, PRFM-vol and PRFM-int agree closely with each other and with TIGRESS for the star formation rate; PRFM-int remains robust at all resolutions tested. These results demonstrate that PRFM-derived subgrid prescriptions provide a physically grounded and numerically stable framework for star formation across the dynamic range of galaxy formation simulations, paving the way for future cosmological applications.

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 / 1 minor

Summary. The paper presents a pressure-regulated, feedback-modulated (PRFM) subgrid model for interstellar gas and star formation in galaxy simulations. It develops two implementations in Arepo (PRFM-vol for resolved scale heights and PRFM-int for unresolved cases), calibrates star formation efficiency and effective EOS directly from TIGRESS parsec-scale simulations, and tests both on isolated Milky Way-like disk runs at mass resolutions 10^5–10^7 M_⊙. Results show shorter gas depletion times than the IllustrisTNG prescription (especially at high pressure/density), close agreement with TIGRESS at high resolution, and robustness of PRFM-int across resolutions. The work frames these as a physically grounded framework paving the way for cosmological applications.

Significance. If the direct transfer of the TIGRESS calibration holds without additional scale-dependent corrections, the model supplies a physically motivated alternative to empirical prescriptions, with explicit strengths in the TIGRESS-based calibration of both efficiency and EOS plus the dual volumetric/integrated implementations that address numerical resolution limits. The isolated-disk tests credibly demonstrate shorter depletion times and resolution robustness relative to TNG. However, because all quantitative results (depletion times, SFR agreement, convergence) are confined to isolated disks, the significance for cosmological runs with accretion and mergers remains prospective rather than demonstrated.

major comments (2)
  1. [Abstract] Abstract: The claim that PRFM-int 'remains robust at all resolutions tested' and that the prescriptions 'provide a physically grounded and numerically stable framework for star formation across the dynamic range of galaxy formation simulations' is load-bearing for the central contribution, yet rests exclusively on isolated MW-like disk simulations; no cosmological runs incorporating accretion, mergers, or varying environments are shown, leaving the assumption of direct transfer without scale-dependent corrections untested.
  2. [Abstract] Abstract: The statement that 'both the star formation efficiency and the effective equation of state may be directly calibrated using numerical simulations, such as TIGRESS' provides no quantitative error bars, exact fitting procedure, or data exclusion criteria for the TIGRESS-derived constants; this directly affects assessment of the reported agreement with TIGRESS SFRs and the shorter depletion times versus TNG.
minor comments (1)
  1. [Abstract] Abstract: Quantitative metrics for 'shorter gas depletion times' and 'agree closely' (e.g., specific depletion time ratios or SFR values with uncertainties) are not supplied, limiting evaluation of the strength of the comparison to TIGRESS and TNG.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their careful review and constructive comments. We address each major comment below, agreeing where revisions are warranted to better reflect the scope of our tests and the calibration details.

read point-by-point responses
  1. Referee: [Abstract] Abstract: The claim that PRFM-int 'remains robust at all resolutions tested' and that the prescriptions 'provide a physically grounded and numerically stable framework for star formation across the dynamic range of galaxy formation simulations' is load-bearing for the central contribution, yet rests exclusively on isolated MW-like disk simulations; no cosmological runs incorporating accretion, mergers, or varying environments are shown, leaving the assumption of direct transfer without scale-dependent corrections untested.

    Authors: We agree that all quantitative results, including resolution robustness and depletion times, are demonstrated exclusively in isolated Milky Way-like disk simulations. The manuscript presents the PRFM model as a framework that paves the way for cosmological applications rather than claiming validation in those regimes. We will revise the abstract to explicitly state that the current tests are limited to isolated disks, that robustness is shown with respect to numerical resolution in these setups, and that cosmological runs with accretion and mergers are left for future work. This removes any implication that direct transfer without scale-dependent corrections has been tested in complex environments. revision: yes

  2. Referee: [Abstract] Abstract: The statement that 'both the star formation efficiency and the effective equation of state may be directly calibrated using numerical simulations, such as TIGRESS' provides no quantitative error bars, exact fitting procedure, or data exclusion criteria for the TIGRESS-derived constants; this directly affects assessment of the reported agreement with TIGRESS SFRs and the shorter depletion times versus TNG.

    Authors: The calibration of the star formation efficiency and effective equation of state from TIGRESS, including the fitting procedure and associated constants, is described in detail in Section 3 of the manuscript. To improve clarity in the abstract, we will add a concise reference to the calibration method and note that quantitative uncertainties and fitting details are provided in the main text. This will allow readers to assess the agreement with TIGRESS and comparisons to TNG more directly from the abstract. revision: yes

Circularity Check

1 steps flagged

PRFM calibration from TIGRESS makes reported agreement with TIGRESS a fitted result by construction

specific steps
  1. fitted input called prediction [Abstract]
    "With this formulation, both the star formation efficiency and the effective equation of state may be directly calibrated using numerical simulations, such as TIGRESS, which resolve physics of the interstellar medium and star formation at parsec scales. [...] At high resolution, PRFM-vol and PRFM-int agree closely with each other and with TIGRESS for the star formation rate"

    The efficiency and EOS parameters are fitted/calibrated from TIGRESS; therefore the reported close agreement with TIGRESS SFRs at high resolution in the Arepo runs is a direct consequence of that calibration rather than an independent prediction.

full rationale

The paper explicitly states that star formation efficiency and effective EOS are calibrated directly from TIGRESS simulations. It then presents close agreement between the Arepo implementations and TIGRESS SFRs at high resolution as a validation result. This agreement is forced by the calibration step rather than an independent test. The comparison to IllustrisTNG and resolution robustness tests retain independent content, and no load-bearing self-citation chain or definitional reduction appears in the derivation of the PRFM framework itself. The assumption of direct transfer to cosmological runs without scale corrections is noted as future work rather than a derived claim.

Axiom & Free-Parameter Ledger

1 free parameters · 1 axioms · 0 invented entities

Central claim rests on transferability of TIGRESS-calibrated PRFM relations to unresolved scales and on the validity of vertical-equilibrium reconstruction for the integrated version; no new entities postulated.

free parameters (1)
  • PRFM calibration constants from TIGRESS
    Star formation efficiency and effective EOS stated to be directly calibrated from TIGRESS; these act as fitted inputs to the subgrid model.
axioms (1)
  • domain assumption Vertical hydrostatic equilibrium allows reconstruction of mid-plane density and pressure when scale height is unresolved
    Invoked for PRFM-int version in the abstract description of the integrated model.

pith-pipeline@v0.9.1-grok · 5851 in / 1408 out tokens · 22104 ms · 2026-06-27T15:48:58.791990+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

61 extracted references · 61 canonical work pages · 6 internal anchors

  1. [1]

    L., Teyssier, R., & Dekel, A

    Andalman, Z. L., Teyssier, R., & Dekel, A. 2025, MNRAS, 540, 3350, doi: 10.1093/mnras/staf930

  2. [2]

    K., S´ anchez, S

    Barrera-Ballesteros, J. K., S´ anchez, S. F., Heckman, T., et al. 2021, MNRAS, 503, 3643, doi: 10.1093/mnras/stab755

  3. [3]

    , keywords =

    Burger, J. D., Springel, V., Ostriker, E. C., et al. 2025, MNRAS, 544, 1390, doi: 10.1093/mnras/staf1720

  4. [4]

    , keywords =

    Casey, C. M., Akins, H. B., Shuntov, M., et al. 2024, ApJ, 965, 98, doi: 10.3847/1538-4357/ad2075

  5. [5]

    Cen, R., & Ostriker, J. P. 1992, ApJL, 399, L113, doi: 10.1086/186620 Dav´ e, R., Angl´ es-Alc´ azar, D., Narayanan, D., et al. 2019, MNRAS, 486, 2827, doi: 10.1093/mnras/stz937

  6. [6]

    , keywords =

    Li, Z. 2023, MNRAS, 523, 3201, doi: 10.1093/mnras/stad1557

  7. [7]

    L., Pan, H.-A., Bluck, A

    Ellison, S. L., Pan, H.-A., Bluck, A. F. L., et al. 2024, MNRAS, 527, 10201, doi: 10.1093/mnras/stad3778

  8. [8]

    Elmegreen, B. G. 1989, ApJ, 338, 178, doi: 10.1086/167192

  9. [9]

    M., Hernquist, L., et al

    Genel, S., Fall, S. M., Hernquist, L., et al. 2015, ApJL, 804, L40, doi: 10.1088/2041-8205/804/2/L40

  10. [10]

    2018, MNRAS, 479, 3042, doi: 10.1093/mnras/sty1653

    Girichidis, P., Naab, T., Hanasz, M., & Walch, S. 2018, MNRAS, 479, 3042, doi: 10.1093/mnras/sty1653

  11. [11]

    , keywords =

    Gurvich, A. B., Faucher-Gigu` ere, C.-A., Richings, A. J., et al. 2020, MNRAS, 498, 3664, doi: 10.1093/mnras/staa2578

  12. [12]

    , keywords =

    Hassan, S., Ostriker, E. C., Kim, C.-G., et al. 2023, MNRAS to be submitted —. 2024, ApJ, 975, 151, doi: 10.3847/1538-4357/ad73a4

  13. [13]

    2019, MNRAS, 483, 3363, doi: 10.1093/mnras/sty3252

    Hu, C.-Y. 2019, MNRAS, 483, 3363, doi: 10.1093/mnras/sty3252

  14. [14]

    Jeffreson, S. M. R., Ostriker, E. C., Kim, C.-G., et al. 2024, ApJ, 975, 113, doi: 10.3847/1538-4357/ad793f

  15. [15]

    Jeffreson, S. M. R., Sun, J., & Wilson, C. D. 2022, MNRAS, 515, 1663, doi: 10.1093/mnras/stac1874

  16. [16]

    1992, ApJ, 391, 502, doi: 10.1086/171366

    Katz, N. 1992, ApJ, 391, 502, doi: 10.1086/171366

  17. [17]

    C., & Evans, N

    Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531, doi: 10.1146/annurev-astro-081811-125610

  18. [18]

    Kennicutt, Jr., R. C. 1998, ApJ, 498, 541, doi: 10.1086/305588

  19. [19]

    Kim, C.-G., Kim, J.-G., Gong, M., & Ostriker, E. C. 2023a, ApJ, 946, 3, doi: 10.3847/1538-4357/acbd3a

  20. [20]

    Kim, C.-G., & Ostriker, E. C. 2015, ApJ, 802, 99, doi: 10.1088/0004-637X/802/2/99 —. 2017, ApJ, 846, 133, doi: 10.3847/1538-4357/aa8599 —. 2018, ApJ, 853, 173, doi: 10.3847/1538-4357/aaa5ff

  21. [21]

    C., & Kim, W.-T

    Kim, C.-G., Ostriker, E. C., & Kim, W.-T. 2013, ApJ, 776, 1, doi: 10.1088/0004-637X/776/1/1

  22. [22]

    , keywords =

    Kim, C.-G., Ostriker, E. C., Somerville, R. S., et al. 2020a, ApJ, 900, 61, doi: 10.3847/1538-4357/aba962

  23. [23]

    , keywords =

    Kim, C.-G., Ostriker, E. C., Fielding, D. B., et al. 2020b, ApJL, 903, L34, doi: 10.3847/2041-8213/abc252

  24. [24]

    , keywords =

    Kim, C.-G., Ostriker, E. C., Kim, J.-G., et al. 2024, ApJ, 972, 67, doi: 10.3847/1538-4357/ad59ab

  25. [25]

    Kim, J.-G., Gong, M., Kim, C.-G., & Ostriker, E. C. 2023b, ApJS, 264, 10, doi: 10.3847/1538-4365/ac9b1d

  26. [26]

    2014, ApJS, 210, 14, doi: 10.1088/0067-0049/210/1/14

    Kim, J.-h., Abel, T., Agertz, O., et al. 2014, ApJS, 210, 14, doi: 10.1088/0067-0049/210/1/14

  27. [27]

    Kim, W.-T., Kim, C.-G., & Ostriker, E. C. 2020c, ApJ, 898, 35, doi: 10.3847/1538-4357/ab9b87 Labb´ e, I., van Dokkum, P., Nelson, E., et al. 2023, Nature, 616, 266, doi: 10.1038/s41586-023-05786-2

  28. [28]

    B., Kim, J.-G., Kim, C.-G., & Ostriker, E

    Linzer, N. B., Kim, J.-G., Kim, C.-G., & Ostriker, E. C. 2024, ApJ, 975, 173, doi: 10.3847/1538-4357/ad7733

  29. [29]

    1994, ApJ, 433, 687, doi: 10.1086/174677 —

    Malhotra, S. 1994, ApJ, 433, 687, doi: 10.1086/174677 —. 1995, ApJ, 448, 138, doi: 10.1086/175946

  30. [30]

    , keywords =

    Mason, C. A., Trenti, M., & Treu, T. 2023, MNRAS, 521, 497, doi: 10.1093/mnras/stad035

  31. [31]

    M., Stanimirovi´ c, S., & Rybarczyk, D

    McClure-Griffiths, N. M., Stanimirovi´ c, S., & Rybarczyk, D. R. 2023, ARA&A, 61, 19, doi: 10.1146/annurev-astro-052920-104851

  32. [32]

    2020, A&A, 643, A78, doi: 10.1051/0004-6361/202039008

    Molina, J., Ibar, E., Godoy, N., et al. 2020, A&A, 643, A78, doi: 10.1051/0004-6361/202039008

  33. [33]

    Naab, T., & Ostriker, J. P. 2017, ARA&A, 55, 59, doi: 10.1146/annurev-astro-081913-040019

  34. [34]

    Monthly Notices of the Royal Astronomical Society , author =

    Hernquist, L. 2012, MNRAS, 421, 3127, doi: 10.1111/j.1365-2966.2012.20536.x

  35. [35]

    2018 MNRAS, 475, 624, doi:10.1093/mnras/stx3040

    Nelson, D., Pillepich, A., Springel, V., et al. 2018, MNRAS, 475, 624, doi: 10.1093/mnras/stx3040

  36. [36]

    , keywords =

    Ostriker, E. C., & Kim, C.-G. 2022, ApJ, 936, 137, doi: 10.3847/1538-4357/ac7de2

  37. [37]

    C., McKee, C

    Ostriker, E. C., McKee, C. F., & Leroy, A. K. 2010, ApJ, 721, 975, doi: 10.1088/0004-637X/721/2/975

  38. [38]

    C., & Shetty, R

    Ostriker, E. C., & Shetty, R. 2011, ApJ, 731, 41, doi: 10.1088/0004-637X/731/1/41

  39. [39]

    , keywords =

    Pandya, V., Fielding, D. B., Angl´ es-Alc´ azar, D., et al. 2021, MNRAS, 508, 2979, doi: 10.1093/mnras/stab2714

  40. [40]
  41. [41]

    , keywords =

    Pillepich, A., Nelson, D., Springel, V., et al. 2019, MNRAS, 490, 3196, doi: 10.1093/mnras/stz2338

  42. [42]

    MNRAS , archivePrefix = "arXiv", eprint =

    Schaye, J., & Dalla Vecchia, C. 2008, MNRAS, 383, 1210, doi: 10.1111/j.1365-2966.2007.12639.x

  43. [43]

    , archivePrefix = "arXiv", eprint =

    Schaye, J., Crain, R. A., Bower, R. G., et al. 2015, MNRAS, 446, 521, doi: 10.1093/mnras/stu2058 31

  44. [44]

    2026, MNRAS, 545, staf2119, doi: 10.1093/mnras/staf2119

    Shen, X., Kannan, R., Puchwein, E., et al. 2026, MNRAS, 545, staf2119, doi: 10.1093/mnras/staf2119

  45. [45]

    2025, ApJ, 987, 204, doi: 10.3847/1538-4357/adda3d

    Weber, M. 2025, ApJ, 987, 204, doi: 10.3847/1538-4357/adda3d

  46. [46]

    A model for unresolved cool clouds entrained in galactic winds in cosmological simulations

    Smith, M. C., Fielding, D. B., Bryan, G. L., et al. 2024a, MNRAS, 535, 3550, doi: 10.1093/mnras/stae2589 —. 2024b, MNRAS, 527, 1216, doi: 10.1093/mnras/stad3168

  47. [47]

    Physical Models of Galaxy Formation in a Cosmological Framework

    Somerville, R. S., & Dav´ e, R. 2015, ARA&A, 53, 51, doi: 10.1146/annurev-astro-082812-140951

  48. [48]

    2009, MNRAS, 400, 1181, doi: 10.1111/j.1365-2966.2009.15528.x

    Springel, V. 2010, MNRAS, 401, 791, doi: 10.1111/j.1365-2966.2009.15715.x

  49. [49]

    , keywords =

    Springel, V., & Hernquist, L. 2003, MNRAS, 339, 289, doi: 10.1046/j.1365-8711.2003.06206.x

  50. [50]

    2018, title First results from the IllustrisTNG simulations: matter and galaxy clustering , , 475, 676, 10.1093/mnras/stx3304

    Springel, V., Pakmor, R., Pillepich, A., et al. 2018, MNRAS, 475, 676, doi: 10.1093/mnras/stx3304

  51. [51]

    P., Kim, C.-G., Bryan, G

    Steinwandel, U. P., Kim, C.-G., Bryan, G. L., et al. 2024, ApJ, 960, 100, doi: 10.3847/1538-4357/ad09e1

  52. [52]

    2021, ApJ, 910, 131, doi: 10.3847/1538-4357/abe5ab

    Su, Y., Yang, J., Yan, Q.-Z., et al. 2021, ApJ, 910, 131, doi: 10.3847/1538-4357/abe5ab

  53. [53]

    K., Ostriker, E

    Sun, J., Leroy, A. K., Ostriker, E. C., et al. 2020, ApJ, 892, 148, doi: 10.3847/1538-4357/ab781c —. 2023, ApJL, 945, L19, doi: 10.3847/2041-8213/acbd9c

  54. [54]

    , keywords =

    Tacconi, L. J., Genzel, R., & Sternberg, A. 2020, ARA&A, 58, 157, doi: 10.1146/annurev-astro-082812-141034

  55. [55]

    M., et al

    Teng, Y.-H., Chiang, I.-D., Sandstrom, K. M., et al. 2024, ApJ, 961, 42, doi: 10.3847/1538-4357/ad10ae

  56. [56]

    C., et al

    Vijayakumar, V., Sun, J., Ostriker, E. C., et al. 2025, ApJ, 989, 66, doi: 10.3847/1538-4357/ade800

  57. [57]

    2013, title A model for cosmological simulations of galaxy formation physics , , 436, 3031, 10.1093/mnras/stt1789

    Vogelsberger, M., Genel, S., Sijacki, D., et al. 2013, MNRAS, 436, 3031, doi: 10.1093/mnras/stt1789

  58. [58]

    doi:10.1038/s42254-019-0127-2 , eprint =

    Vogelsberger, M., Marinacci, F., Torrey, P., & Puchwein, E. 2020, Nature Reviews Physics, 2, 42, doi: 10.1038/s42254-019-0127-2

  59. [59]

    Introducing the

    Vogelsberger, M., Genel, S., Springel, V., et al. 2014, MNRAS, 444, 1518, doi: 10.1093/mnras/stu1536

  60. [60]

    A., Elbaz, D., et al

    Xiao, M., Oesch, P. A., Elbaz, D., et al. 2024, Nature, 635, 311, doi: 10.1038/s41586-024-08094-5

  61. [61]

    2026, ApJ, 999, 41, doi: 10.3847/1538-4357/ae3c08

    Zhou, Y., Di Matteo, T., Bird, S., et al. 2026, ApJ, 999, 41, doi: 10.3847/1538-4357/ae3c08