Two-phase formation of galaxies: the coevolution between galaxies and dark matter halos
Pith reviewed 2026-05-18 15:13 UTC · model grok-4.3
The pith
Milky Way-size galaxies form in two phases linked to when their dark matter halos switch from fast to slow accretion.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The formation of these galaxies follows a two-phase pattern, with an early phase featured by hot dynamics, bulge-dominated structure and bursty star formation, and a later phase featured by cold dynamics, disk-dominated structure and steady star formation. The transition times of these galaxy properties are correlated with the time when the host halo transits from fast to slow accretion, indicating the two-phase assembly of halos as a potential mechanism that drives the two-phase formation of galaxies. The physical origin of dynamical hotness can be summarized into two modes of star formation: a scattered mode in which stars form at large radii within cold gas streams associated with fast a
What carries the argument
The two-phase assembly of dark matter halos, which sets the timing for two star-formation modes: scattered formation in distant cold gas streams during rapid halo growth and concentrated formation via violent gas fragmentation near the center as growth slows.
If this is right
- The two star-formation modes produce distinct structural signatures in high-redshift galaxies.
- Cold gaseous and stellar disks emerge only after fast halo accretion stalls and feedback reduces available gas.
- Halo assembly history directly controls the morphological transition from bulge to disk dominance.
- These patterns offer testable predictions for the internal structures of distant galaxies.
Where Pith is reading between the lines
- The same halo-driven timing could organize morphological changes across a wider range of galaxy masses if the fast-to-slow transition is common.
- Comparing the predicted structural imprints against early JWST images of high-redshift galaxies would provide an independent check.
- Varying feedback efficiency in follow-up simulations could reveal how sensitive the correlation is to the details of gas removal.
Load-bearing premise
The subgrid prescriptions for star formation, feedback, and gas cooling in the simulations accurately capture the real processes that control the timing and outcomes of each phase without major numerical artifacts.
What would settle it
Direct measurements showing that the observed transition times for galaxy dynamical state or star-formation rate do not align with the inferred fast-to-slow accretion transition times of their host halos.
Figures
read the original abstract
We use FIRE-2 cosmological zoom-in hydrodynamic simulations to investigate the co-evolution between Milky Way-size galaxies and their host dark matter halos. We find that the formation of these galaxies follows a two-phase pattern, with an early phase featured by hot dynamics, bulge-dominated structure and bursty star formation, and a later phase featured by cold dynamics, disk-dominated structure and steady star formation. The transition times of these galaxy properties are correlated with the time when the host halo transits from fast to slow accretion, indicating the two-phase assembly of halos as a potential mechanism that drives the two-phase formation of galaxies. The physical origin of dynamical hotness can be summarized into two modes of star formation: a scattered mode in which stars form at large radii within cold gas streams associated with fast assembly of halos, and a concentrated mode in which stars form at small radii through violent fragmentation from globally self-gravitated gas when halo assembly is about to slow down. Cold gaseous and stellar disks can form when the conditions of the two modes are removed by the stall of fast halo assembly and the reduction of gas by feedback processes. The two modes of star formation leave distinct imprints on the structural properties of high-redshift galaxies, providing implications to be tested by JWST and future observations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes FIRE-2 cosmological zoom-in hydrodynamic simulations of Milky Way-size galaxies to investigate their co-evolution with host dark matter halos. It reports a two-phase galaxy formation history: an early phase with hot stellar dynamics, bulge-dominated structure, and bursty star formation, transitioning to a later phase with cold dynamics, disk-dominated structure, and steady star formation. The transition times of these galaxy properties are found to correlate with the epoch at which the host halo shifts from fast to slow accretion. The authors identify two star-formation modes (scattered in cold streams during fast assembly; concentrated via global fragmentation as assembly slows) as the origin of the early dynamical hotness, with implications for high-redshift observations.
Significance. If the reported correlations prove robust, the work supplies a concrete physical link between dark-matter halo assembly phases and the emergence of galactic disks and steady star formation. The explicit mapping of two star-formation modes onto halo accretion history yields falsifiable predictions for the structural properties of high-redshift galaxies that can be tested with JWST. The strength of the analysis lies in the use of high-resolution zoom-in runs that track both baryonic and dark-matter quantities simultaneously within the same volumes.
major comments (2)
- [Section 3] Section 3 (results on transition timing): The central claim that halo fast-to-slow accretion drives the galaxy two-phase transition rests on the alignment of transition epochs measured in the FIRE-2 runs. Because the star-formation threshold, supernova feedback, and metal-line cooling prescriptions are fixed, the manuscript should demonstrate that the reported correlation persists under controlled variations of these subgrid parameters or under comparison with at least one other simulation suite; otherwise the alignment could be an artifact of the specific baryonic implementation rather than a general consequence of halo assembly.
- [Section 2] Section 2 (methods): The quantitative definitions used to identify the dynamical hot-to-cold transition, the bursty-to-steady star-formation transition, and the bulge-to-disk transition are load-bearing for the correlation analysis. Explicit criteria, including any thresholds on velocity dispersion, specific star-formation rate variability, or morphological parameters, together with tests of sensitivity to sample selection or redshift binning, must be provided so that readers can assess whether post-hoc choices influence the strength of the reported halo-galaxy timing correlation.
minor comments (3)
- [Figure 1] Figure 1 and Figure 3: Adding vertical lines or shaded bands marking the identified transition times on the time-evolution panels would make the visual correspondence between galaxy and halo transitions immediately apparent.
- The abstract states that the two modes leave 'distinct imprints' on high-redshift galaxies but does not quantify the expected differences in size, clumpiness, or kinematics; a short summary table or bullet list in the discussion section would improve clarity.
- Notation: The terms 'scattered mode' and 'concentrated mode' are introduced without a concise definition box or table; adding one would aid readers who wish to compare the modes directly with other simulation studies.
Simulated Author's Rebuttal
We thank the referee for their constructive and insightful comments on our manuscript. We address each major comment in detail below and have revised the manuscript to incorporate clarifications and additional analyses where feasible.
read point-by-point responses
-
Referee: [Section 3] Section 3 (results on transition timing): The central claim that halo fast-to-slow accretion drives the galaxy two-phase transition rests on the alignment of transition epochs measured in the FIRE-2 runs. Because the star-formation threshold, supernova feedback, and metal-line cooling prescriptions are fixed, the manuscript should demonstrate that the reported correlation persists under controlled variations of these subgrid parameters or under comparison with at least one other simulation suite; otherwise the alignment could be an artifact of the specific baryonic implementation rather than a general consequence of halo assembly.
Authors: We thank the referee for raising this important issue of robustness. The subgrid prescriptions in FIRE-2 are indeed fixed, and the reported correlations are measured within this specific implementation. However, the two star-formation modes we describe arise directly from the gravitational dynamics of cold gas streams and global fragmentation, which are governed by the halo's mass accretion rate rather than the precise details of feedback or cooling. We will add a dedicated discussion subsection in the revised manuscript that outlines why these modes should be general, drawing on analytic expectations for gas collapse during fast versus slow accretion and citing supporting results from other simulation studies in the literature. A systematic parameter variation or cross-suite comparison lies outside the scope of this work, as it would require new high-resolution simulations. revision: partial
-
Referee: [Section 2] Section 2 (methods): The quantitative definitions used to identify the dynamical hot-to-cold transition, the bursty-to-steady star-formation transition, and the bulge-to-disk transition are load-bearing for the correlation analysis. Explicit criteria, including any thresholds on velocity dispersion, specific star-formation rate variability, or morphological parameters, together with tests of sensitivity to sample selection or redshift binning, must be provided so that readers can assess whether post-hoc choices influence the strength of the reported halo-galaxy timing correlation.
Authors: We agree that explicit, reproducible definitions are necessary. In the revised manuscript we will expand Section 2 to state the precise criteria: the dynamical hot-to-cold transition is defined as the redshift at which the ratio of stellar velocity dispersion within 2 R_e to the circular velocity at R_e drops below 0.4 and remains below this value for at least 1 Gyr; the bursty-to-steady transition occurs when the 1-sigma scatter in log(sSFR) measured in 100-Myr windows falls below 0.25 dex; and the bulge-to-disk transition is identified when the kinematically defined disk mass fraction exceeds 0.5. We will also include sensitivity tests that vary each threshold by ±15 %, alter the minimum sample size, and change the redshift binning, demonstrating that the reported timing correlations with halo accretion remain statistically significant (Spearman coefficients > 0.7) under these variations. These definitions and tests will appear in the main text with supporting figures in an appendix. revision: yes
- Direct demonstration that the reported correlations persist under controlled variations of subgrid parameters or through comparison with another simulation suite, which would require new simulations beyond the existing FIRE-2 dataset.
Circularity Check
Simulation outputs yield measured correlations without definitional reduction
full rationale
The paper reports direct measurements of transition times in galaxy properties (dynamical hot-to-cold, bursty-to-steady SF, bulge-to-disk) extracted from FIRE-2 zoom-in runs and compares them to the independently computed halo mass-accretion history. No parameter is fitted to enforce the reported correlation, no galaxy transition epoch is defined using the halo transition time, and no self-citation chain is invoked to justify the central alignment. The two star-formation modes are described as emergent behaviors observed in the simulated volumes rather than imposed by construction. The analysis therefore remains self-contained against external benchmarks and does not reduce any claimed result to its own inputs.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Standard Lambda-CDM cosmology governs the initial conditions and gravitational evolution of the simulations.
- domain assumption FIRE-2 subgrid models for star formation, stellar feedback, and gas cooling accurately capture the relevant physics at the resolved scales.
Lean theorems connected to this paper
-
IndisputableMonolith/Foundation/DimensionForcing.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The transition times of these galaxy properties are correlated with the time when the host halo transits from fast to slow accretion, indicating the two-phase assembly of halos as a potential mechanism that drives the two-phase formation of galaxies.
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We use FIRE-2 cosmological zoom-in hydrodynamic simulations... subgrid prescriptions for star formation, feedback, and gas cooling
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]
Abadi,M.G.,Navarro,J.F.,Steinmetz,M.,&Eke,V.R.2003,ApJ,597,21, doi: 10.1086/378316 Adamo,A.,Bradley,L.D.,Vanzella,E.,etal.2024,Nature,632,513,doi:
-
[2]
1038/s41586-024-07703-7 Agertz,O.,&Kravtsov,A.V.2015,ApJ,804,18,doi: 10.1088/0004-637X/ 804/1/18 —. 2016, ApJ, 824, 79, doi:10.3847/0004-637X/824/2/79 Allgood, B., Flores, R. A., Primack, J. R., et al. 2006, MNRAS, 367, 1781, doi: 10.1111/j.1365-2966.2006.10094.x Baldry, I. K., Glazebrook, K., Brinkmann, J., et al. 2004, ApJ, 600, 681, doi: 10.1086/380092...
-
[3]
1088/0004-637X/763/1/18 Bezanson, R., van Dokkum, P. G., Tal, T., et al. 2009, ApJ, 697, 1290, doi: 10.1088/0004-637X/697/2/1290 Boco, L., Lapi, A., Shankar, F., et al. 2023, ApJ, 954, 97, doi:10.3847/ 1538-4357/ace76d Borzyszkowski, M., Porciani, C., Romano-Díaz, E., & Garaldi, E. 2017, MNRAS, 469, 594, doi:10.1093/mnras/stx873 Bryan, G. L., & Norman, M....
-
[4]
2016, ARA&A, 54, 597, doi: 10.1146/ annurev-astro-082214-122432 Cappellari, M., McDermid, R
1046/j.1365-8711.2001.04068.x Cappellari, M. 2016, ARA&A, 54, 597, doi: 10.1146/ annurev-astro-082214-122432 Cappellari, M., McDermid, R. M., Alatalo, K., et al. 2013, MNRAS, 432, 1862, doi:10.1093/mnras/stt644 Carr,D.S.,Kannappan,S.J.,Hutchens,Z.L.,etal.2025,TheAstrophysical Journal, 985, 25, doi:10.3847/1538-4357/adc71b Carr, D. S., Kannappan, S. J., No...
-
[5]
1093/mnras/stab695 Chen,Y.,Mo,H.J.,Li,C.,etal.2021b,MNRAS,507,2510,doi: 10.1093/ mnras/stab2377 Cleveland, W. S., & Devlin, S. J. 1988, Journal of the American Statistical Association, 83, 596, doi:10.1080/01621459.1988.10478639 Cohn, J. H., Durodola, E., Casey, Q. O., Lambrides, E., & Hickox, R. C. 2025,TheAstrophysicalJournal,988,L61,doi: 10.3847/2041-8...
-
[6]
2018a, MNRAS, 473, 1930, doi: 10.1093/mnras/stx2482
1093/mnras/199.4.1069 El-Badry, K., Wetzel, A., Geha, M., et al. 2016, ApJ, 820, 131, doi:10. 3847/0004-637X/820/2/131 El-Badry, K., Quataert, E., Wetzel, A., et al. 2018, MNRAS, 473, 1930, doi: 10.1093/mnras/stx2482 Escala,I.,Wetzel,A.,Kirby,E.N.,etal.2018,MNRAS,474,2194,doi:
-
[7]
1093/mnras/stx2858 Faucher-Giguère, C.-A., Lidz, A., Zaldarriaga, M., & Hernquist, L. 2009, ApJ, 703, 1416, doi:10.1088/0004-637X/703/2/1416 Faucher-Giguère, C.-A., Quataert, E., & Hopkins, P. F. 2013, MNRAS, 433, 1970, doi:10.1093/mnras/stt866 Flores Velázquez, J. A., Gurvich, A. B., Faucher-Giguère, C.-A., et al. 2021, MNRAS, 501, 4812, doi:10.1093/mnra...
-
[8]
F., Wetzel, A., Kereš, D., et al
1093/mnras/sty1690 Hopkins, P. F., Wetzel, A., Kereš, D., et al. 2018, MNRAS, 477, 1578, doi: 10.1093/mnras/sty674 MNRAS000, 1–21 (2025) 20 Qinglin Ma et al. Hopkins, P. F., Gurvich, A. B., Shen, X., et al. 2023, MNRAS, 525, 2241, doi: 10.1093/mnras/stad1902 Huang, K.-H., Fall, S. M., Ferguson, H. C., et al. 2017, ApJ, 838, 6, doi:10. 3847/1538-4357/aa62a...
-
[9]
https://arxiv.org/abs/2505.05720 Kado-Fong, E., Sanderson, R. E., Greene, J. E., et al. 2022, ApJ, 931, 152, doi: 10.3847/1538-4357/ac6c88 Kassin,S.A.,Weiner,B.J.,Faber,S.M.,etal.2012,ApJ,758,106,doi:
-
[10]
1088/0004-637X/758/2/106 Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2005, MNRAS, 363, 2, doi: 10.1111/j.1365-2966.2005.09451.x Kluge,M.,&Bender,R.2023,TheAstrophysicalJournalSupplementSeries, 267, 41, doi:10.3847/1538-4365/ace052 Kormendy, J. 1979, ApJ, 227, 714, doi:10.1086/156782 Kormendy, J., & Ho, L. C. 2013, Annu. Rev. Astron. Astrophys., 51, 5...
-
[11]
1093/mnras/staa527 Mandelker, N., van Dokkum, P. G., Brodie, J. P., van den Bosch, F. C., & Ceverino, D. 2018, ApJ, 861, 148, doi:10.3847/1538-4357/aaca98 McCluskey,F.,Wetzel,A.,Loebman,S.R.,etal.2024,MNRAS,527,6926, doi: 10.1093/mnras/stad3547 Meng,X.,&Gnedin,O.Y.2020,MNRAS,494,1263,doi: 10.1093/mnras/ staa776 —. 2021, MNRAS, 502, 1433, doi:10.1093/mnras...
-
[12]
1086/177173 —. 1997, ApJ, 490, 493, doi:10.1086/304888 Nelson, D., Springel, V., Pillepich, A., et al. 2019, Computational Astro- physics and Cosmology, 6, 2, doi:10.1186/s40668-019-0028-x Newman, A. B., Ellis, R. S., Bundy, K., & Treu, T. 2012, ApJ, 746, 162, doi: 10.1088/0004-637X/746/2/162 Oñorbe, J., Garrison-Kimmel, S., Maller, A. H., et al. 2014, MN...
-
[13]
The COLIBRE project: cosmological hydrodynamical simulations of galaxy formation and evolution
1093/mnras/stz3054 Sandage, A. 1961, The Hubble Atlas of Galaxies Schaye, J., Chaikin, E., Schaller, M., et al. 2025, arXiv e-prints, arXiv:2508.21126, doi:10.48550/arXiv.2508.21126 Shen,S.,Mo,H.J.,White,S.D.M.,etal.2003,MNRAS,343,978,doi:
work page internal anchor Pith review Pith/arXiv arXiv doi:10.48550/arxiv.2508.21126 1961
- [14]
-
[15]
S., Behroozi, P., Pandya, V., et al
3847/1538-4357/aa740c Somerville, R. S., Behroozi, P., Pandya, V., et al. 2018, MNRAS, 473, 2714, doi: 10.1093/mnras/stx2040 Stern,J.,Fielding,D.,Faucher-Giguère,C.-A.,&Quataert,E.2020,MNRAS, 492, 6042, doi:10.1093/mnras/staa198 Strateva, I., Ivezić, Ž., Knapp, G. R., et al. 2001, AJ, 122, 1861, doi:10. 1086/323301 Su, K.-Y., Hopkins, P. F., Hayward, C. C...
-
[16]
1088/2041-8205/792/1/L6 vanDokkum,P.G.,Franx,M.,Kriek,M.,etal.2008,ApJ,677,L5,doi:
work page 2041
-
[17]
1086/587874 van Dokkum, P. G., Nelson, E. J., Franx, M., et al. 2015, ApJ, 813, 23, doi: 10.1088/0004-637X/813/1/23 van Donkelaar, F., Agertz, O., & Renaud, F. 2022, MNRAS, 512, 3806, MNRAS000, 1–21 (2025) Dynamical hotness 21 doi: 10.1093/mnras/stac692 Vanzella, E., Claeyssens, A., Welch, B., et al. 2023, ApJ, 945, 53, doi:10. 3847/1538-4357/acb59a Wechs...
-
[18]
3847/2041-8205/827/2/L23 Wisnioski,E.,FörsterSchreiber,N.M.,Wuyts,S.,etal.2015,ApJ,799,209, doi: 10.1088/0004-637X/799/2/209 Woo, J., Dekel, A., Faber, S. M., & Koo, D. C. 2015, MNRAS, 448, 237, doi: 10.1093/mnras/stu2755 Xia, L. Y., Klein, C., Bullock, J. S., et al. 2025, arXiv e-prints, arXiv:2506.08508, doi:10.48550/arXiv.2506.08508 Yang, H., Gao, L., ...
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.