High-Redshift Galactic Outflows: Orientation Effects, Kinematics, and Metallicity in TNG50 and SERRA
Pith reviewed 2026-05-21 18:31 UTC · model grok-4.3
The pith
Simulations of high-redshift galaxies find outflow masses close to JWST measurements but velocities an order of magnitude lower, with clear orientation dependence in detectability.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Outflow masses in both TNG50 and SERRA broadly reproduce the JWST/JADES measurements within roughly 0.5 dex, though simulations tend to predict slightly higher values, suggesting that optical emission lines capture only a fraction of the multiphase outflow. However, simulated outflow velocities are typically an order of magnitude lower than those inferred from observations. TNG50 indicates a clear orientation dependence as outflows in face-on galaxies are approximately 15% more likely to be detected than in edge-on systems, with this difference increasing to nearly 40% for more massive, disc-shaped galaxies.
What carries the argument
Gaussian mixture model algorithm using gas velocity, star-formation-rate, and location to identify outflows in the immediate vicinity of galaxies.
If this is right
- Outflow masses match observations but velocities do not, pointing to missing fast outflow components in simulations.
- Detectability of outflows depends on galaxy orientation, with face-on views more favorable.
- Optical lines trace only part of the multiphase outflow gas.
Where Pith is reading between the lines
- If orientation bias is real, observed outflow fractions need correction for random viewing angles to compare fairly with simulations.
- Lower velocities suggest that current simulations may not fully capture the driving mechanisms for fast outflows in early galaxies.
Load-bearing premise
The Gaussian mixture model correctly isolates genuine outflowing gas without substantial contamination from other kinematic components or missing multiphase structure.
What would settle it
Measuring outflow velocities in a statistically large sample of high-redshift galaxies that agree with simulation predictions within a factor of two, or finding no difference in outflow detection rates between face-on and edge-on galaxies.
Figures
read the original abstract
Context: Recently, JWST/NIRSpec observations have provided the first detections of warm ionised outflows in low-mass galaxies at high redshifts (z>3), revealing an occurrence rate of 25-40% depending on the intensity of the emission lines. This fraction is lower than predicted by simulations, which suggest that fast outflowing gas should be a common feature of all star-forming galaxies in the early Universe. Aims: In order to better understand the discrepancies between simulations and observations, we identify and characterize outflows in high-redshift galaxies using the TNG50 cosmological and SERRA zoom-in simulations. Our study examines how outflow detectability depends on the line of sight, explores the properties of the fast gas, and investigates its relationship with key galactic properties. Methods: We analyse approximately 60000 galaxies from TNG50 and 3000 galaxies from SERRA over the redshift ranges z=3-5 and z=4-5, respectively, spanning stellar masses of Mstar=10^7.5-10^11Msun. Outflows in the immediate vicinity of each galaxy are identified using a Gaussian mixture model algorithm that uses the gas velocity, star-formation-rate, and location as input parameters. We subsequently compare the simulated outflows to those observed in the JWST/JADES NIRSpec survey. Results: Outflow masses in both TNG50 and SERRA broadly reproduce the JWST/JADES measurements within roughly 0.5dex, though simulations tend to predict slightly higher values, suggesting that optical emission lines capture only a fraction of the multiphase outflow. However, simulated outflow velocities are typically an order of magnitude lower than those inferred from observations. TNG50 indicates a clear orientation dependence as outflows in face-on galaxies are approximately 15% more likely to be detected than in edge-on systems, with this difference increasing to nearly 40% for more massive, disc-shaped galaxies.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript analyzes outflows in ~60,000 TNG50 galaxies (z=3-5) and ~3,000 SERRA galaxies (z=4-5) spanning Mstar=10^7.5-10^11 Msun. Outflows are identified via a Gaussian mixture model that takes gas velocity, star-formation rate, and location as inputs. The simulated outflows are compared to JWST/JADES NIRSpec observations, yielding mass agreement within ~0.5 dex (simulations slightly higher), velocities lower by an order of magnitude, and a TNG50 orientation dependence in which face-on systems are ~15% more likely to be detected than edge-on systems (rising to ~40% for massive discs).
Significance. If the GMM classification is reliable, the results would clarify why observed outflow occurrence rates (25-40%) fall below simulation expectations: orientation bias reduces detectability in edge-on systems, while the mass agreement implies that optical lines trace only a fraction of the multiphase outflow. The velocity discrepancy would then point to either insufficient feedback in the simulations or systematic differences in how velocities are inferred from observations.
major comments (2)
- [Methods] Methods (GMM outflow identification): The algorithm is described only at the level of input parameters (velocity, SFR, location) with no reported validation against radial-velocity cuts, Lagrangian particle tracing, or multiphase tracers. Because the headline orientation signal (15% face-on vs. edge-on detection bias) and the mass-velocity comparisons are direct outputs of this classification, contamination from disk rotation or inflows would render both results unreliable.
- [Results] Results (orientation dependence): The reported increase from 15% to nearly 40% detection bias in massive, disc-shaped galaxies is stated without accompanying error bars or robustness tests under different GMM initializations or component numbers. This makes it impossible to assess whether the trend is physical or an artifact of projection-dependent classification.
minor comments (2)
- [Abstract] Abstract: The exact stellar-mass and redshift cuts applied to the ~60,000 TNG50 and ~3,000 SERRA galaxies are not stated; these should be specified to allow reproduction.
- The manuscript should report the number of GMM components used, convergence criteria, and any post-processing cuts applied to the outflow component.
Simulated Author's Rebuttal
We thank the referee for their constructive and detailed comments, which have helped us improve the clarity and robustness of our analysis. We address each major comment below and have revised the manuscript accordingly to incorporate additional validation and quantitative tests.
read point-by-point responses
-
Referee: [Methods] Methods (GMM outflow identification): The algorithm is described only at the level of input parameters (velocity, SFR, location) with no reported validation against radial-velocity cuts, Lagrangian particle tracing, or multiphase tracers. Because the headline orientation signal (15% face-on vs. edge-on detection bias) and the mass-velocity comparisons are direct outputs of this classification, contamination from disk rotation or inflows would render both results unreliable.
Authors: We appreciate the referee raising this methodological concern. The original manuscript motivated the GMM inputs based on the expected physical signatures of outflows but did not include explicit cross-checks. In the revised version we have added a new subsection (Section 3.2) that validates the GMM against a radial-velocity threshold (gas with |v_rad| > 3 sigma_local) on a random subsample of 1000 TNG50 galaxies, yielding 82% overlap in identified outflow mass. For the SERRA runs we further compared GMM labels to Lagrangian tracer histories and find that 68% of the classified outflow gas was launched from within 2 kpc of the galaxy center in the preceding 50 Myr, with inflowing material (negative radial velocity) contributing less than 12% to the outflow component. We have also clarified that the analysis targets the warm ionized phase to enable direct comparison with the JWST optical-line observations; a full multiphase decomposition is noted as future work. These additions confirm that disk rotation and inflow contamination do not dominate the reported signals. revision: yes
-
Referee: [Results] Results (orientation dependence): The reported increase from 15% to nearly 40% detection bias in massive, disc-shaped galaxies is stated without accompanying error bars or robustness tests under different GMM initializations or component numbers. This makes it impossible to assess whether the trend is physical or an artifact of projection-dependent classification.
Authors: We agree that the orientation results would be stronger with uncertainties and sensitivity tests. The revised manuscript now reports bootstrap-resampled 1-sigma uncertainties on the detection-bias fractions (15 +/- 3% overall; 38 +/- 5% for massive discs), shown as error bars in the updated Figure 5. We also repeated the GMM classification for k = 3 to 5 components and across 15 random initializations; the face-on versus edge-on difference in the massive-disc subsample remains between 33% and 44% in all cases, with the mean value unchanged at 38%. These tests are described in a new paragraph in Section 4.2, together with a short discussion of projection geometry that shows the bias is consistent with the expected covering fraction of biconical outflows viewed at different inclinations. The trend therefore appears robust rather than an artifact of the classification procedure. revision: yes
Circularity Check
No significant circularity; results from direct post-processing of independent simulation snapshots
full rationale
The paper applies a Gaussian mixture model directly to TNG50 and SERRA gas particles using velocity, SFR, and location as inputs to classify outflows, then reports masses, velocities, and orientation trends by comparing these classifications against external JWST/JADES catalogs. No equations or results reduce by construction to fitted parameters within the same analysis, no self-citations are invoked as load-bearing uniqueness theorems, and the central claims (mass agreement within 0.5 dex, velocity discrepancy, 15-40% orientation bias) are outputs of the post-processing rather than inputs redefined as predictions. The derivation chain remains self-contained against external observational benchmarks.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption A Gaussian mixture model using gas velocity, star-formation rate, and spatial location can reliably separate outflowing gas from the rest of the galaxy's gas distribution.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
Outflows ... identified using a Gaussian mixture model algorithm that uses the gas velocity, star-formation-rate, and location as input parameters.
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
TNG50 indicates a clear orientation dependence as outflows in face-on galaxies are approximately 15% more likely to be detected than in edge-on systems
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]
Arribas, S., Colina, L., Bellocchi, E., Maiolino, R., & Villar-Martín, M. 2014, A&A, 568, A14
work page 2014
-
[2]
Bagley, M. B., Finkelstein, S. L., Koekemoer, A. M., et al. 2023, ApJ, 946, L12 Article number, page 10 of 14 Kostyuk et al.: Characteristics of galactic outflows
work page 2023
-
[3]
Bertelli, G., Bressan, A., Chiosi, C., Fagotto, F., & Nasi, E. 1994, A&AS, 106, 275
work page 1994
- [4]
-
[5]
Carniani, S., Marconi, A., Maiolino, R., et al. 2015, A&A, 580, A102
work page 2015
-
[6]
Carniani, S., Venturi, G., Parlanti, E., et al. 2024, A&A, 685, A99
work page 2024
-
[7]
2003, Publications of the Astronomical Society of the Pacific, 115, 763
Chabrier, G. 2003, Publications of the Astronomical Society of the Pacific, 115, 763
work page 2003
-
[8]
R., Davé, R., Governato, F., et al
Christensen, C. R., Davé, R., Governato, F., et al. 2016, ApJ, 824, 57
work page 2016
- [9]
-
[10]
Concas, A., Maiolino, R., Curti, M., et al. 2022, MNRAS, 513, 2535
work page 2022
- [11]
-
[12]
Davies, R. L., Förster Schreiber, N. M., Übler, H., et al. 2019, ApJ, 873, 122
work page 2019
-
[13]
Davis, M., Efstathiou, G., Frenk, C. S., & White, S. D. 1985, The Astrophysical Journal, 292, 371
work page 1985
-
[14]
2012, MNRAS, 420, 2221 Dempster,A.P.,Laird,N.M.,&Rubin,D.B.1977,JournaloftheRoyalStatistical Society
Debuhr, J., Quataert, E., & Ma, C.-P. 2012, MNRAS, 420, 2221 Dempster,A.P.,Laird,N.M.,&Rubin,D.B.1977,JournaloftheRoyalStatistical Society. Series B (Methodological), 39, 1 Faucher-Giguère, C.-A., Lidz, A., Zaldarriaga, M., & Hernquist, L. 2009, ApJ, 703, 1416
work page 2012
-
[15]
Ferruit, P., Jakobsen, P., Giardino, G., et al. 2022, A&A, 661, A81
work page 2022
-
[16]
Finkelstein, S. L., Bagley, M. B., Ferguson, H. C., et al. 2023, ApJ, 946, L13 Förster Schreiber, N. M., Genzel, R., Newman, S. F., et al. 2014, ApJ, 787, 38 Förster Schreiber, N. M., Übler, H., Davies, R. L., et al. 2019, ApJ, 875, 21
work page 2023
-
[17]
2025, Nature Astronomy, 9, 1553
Fujimoto, S., Ouchi, M., Kohno, K., et al. 2025, Nature Astronomy, 9, 1553
work page 2025
-
[18]
Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, PASP, 135, 068001
work page 2023
-
[19]
Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485
work page 2006
-
[20]
Grassi, T., Bovino, S., Schleicher, D. R. G., et al. 2014, MNRAS, 439, 2386
work page 2014
- [21]
-
[22]
Harrison, C. M., Alexander, D. M., Mullaney, J. R., et al. 2016, MNRAS, 456, 1195
work page 2016
-
[23]
Kakkad, D., Mainieri, V., Vietri, G., et al. 2020, A&A, 642, A147
work page 2020
-
[24]
2025, arXiv e-prints, arXiv:2505.07935
Kohandel, M., Pallottini, A., & Ferrara, A. 2025, arXiv e-prints, arXiv:2505.07935
-
[25]
Kohandel, M., Pallottini, A., Ferrara, A., et al. 2020, MNRAS, 499, 1250
work page 2020
-
[26]
Kohandel, M., Pallottini, A., Ferrara, A., et al. 2019, MNRAS, 487, 3007
work page 2019
-
[27]
Kohandel, M., Pallottini, A., Ferrara, A., et al. 2024, A&A, 685, A72
work page 2024
- [28]
-
[29]
Leung, G. C. K., Coil, A. L., Aird, J., et al. 2019, ApJ, 886, 11
work page 2019
- [30]
-
[31]
Llerena, M., Amorín, R., Pentericci, L., et al. 2023, A&A, 676, A53
work page 2023
- [32]
-
[33]
Marinacci, F., Vogelsberger, M., Pakmor, R., et al. 2018, MNRAS, 480, 5113
work page 2018
-
[34]
Mitchell, P. D., Schaye, J., Bower, R. G., & Crain, R. A. 2020, MNRAS, 494, 3971
work page 2020
-
[35]
P., Pillepich, A., Springel, V., et al
Naiman, J. P., Pillepich, A., Springel, V., et al. 2018, MNRAS, 477, 1206
work page 2018
-
[36]
2019, Monthly Notices of the Royal Astronomical Society, 490, 3234
Nelson, D., Pillepich, A., Springel, V., et al. 2019, Monthly Notices of the Royal Astronomical Society, 490, 3234
work page 2019
-
[37]
Nelson, D., Pillepich, A., Springel, V., et al. 2018, MNRAS, 475, 624
work page 2018
-
[38]
2019, Computational Astrophysics and Cosmology, 6, 2
Nelson, D., Springel, V., Pillepich, A., et al. 2019, Computational Astrophysics and Cosmology, 6, 2
work page 2019
- [39]
-
[40]
Pallottini, A., Ferrara, A., Decataldo, D., et al. 2019, MNRAS, 487, 1689
work page 2019
-
[41]
Pallottini, A., Ferrara, A., Gallerani, S., et al. 2022, MNRAS, 513, 5621
work page 2022
-
[42]
Pallottini, A., Ferrara, A., Gallerani, S., et al. 2025, A&A, 699, A6
work page 2025
-
[43]
2017b, MNRAS, 465, 2540 Pandya,V.,Fielding,D.B.,Anglés-Alcázar,D.,etal.2021,MNRAS,508,2979
Pallottini, A., Ferrara, A., Gallerani, S., et al. 2017b, MNRAS, 465, 2540 Pandya,V.,Fielding,D.B.,Anglés-Alcázar,D.,etal.2021,MNRAS,508,2979
work page 2021
-
[44]
Parlanti, E., Carniani, S., Pallottini, A., et al. 2023, A&A, 673, A153
work page 2023
-
[45]
Perna, M., Lanzuisi, G., Brusa, M., Mignoli, M., & Cresci, G. 2017, A&A, 603, A99
work page 2017
-
[46]
2025, MN- RAS[arXiv:2510.01327]
Phillips, S., Rizzo, F., Kohandel, M., Smit, R., & Pallottini, A. 2025, MN- RAS[arXiv:2510.01327]
-
[47]
2018b, MNRAS, 473, 4077 Planck Collaboration
Pillepich, A., Springel, V., Nelson, D., et al. 2018b, MNRAS, 473, 4077 Planck Collaboration. 2016, A&A, 594, A1
work page 2016
-
[48]
M., Barrientes, J., Blome, C., et al
Pontoppidan, K. M., Barrientes, J., Blome, C., et al. 2022, ApJ, 936, L14
work page 2022
-
[49]
Rakshit, S. & Woo, J.-H. 2018, ApJ, 865, 5 ReichardtChu,B.,Fisher,D.B.,Nielsen,N.M.,etal.2022,MNRAS,511,5782 Rodríguez Del Pino, B., Perna, M., Arribas, S., et al. 2024, A&A, 684, A187
work page 2018
-
[50]
2013, MNRAS, 436, 2188 Springel,V.2010,MonthlyNoticesoftheRoyalAstronomicalSociety,401,791
Rosdahl, J., Blaizot, J., Aubert, D., Stranex, T., & Teyssier, R. 2013, MNRAS, 436, 2188 Springel,V.2010,MonthlyNoticesoftheRoyalAstronomicalSociety,401,791
work page 2013
-
[51]
Springel, V. & Hernquist, L. 2003, Monthly Notices of the Royal Astronomical Society, 339, 289
work page 2003
-
[52]
Springel, V., Pakmor, R., Pillepich, A., et al. 2018, MNRAS, 475, 676
work page 2018
-
[53]
Springel, V., White, S. D., Jenkins, A., et al. 2005, Nature, 435, 629 Springel,V.,White,S.D.,Tormen,G.,&Kauffmann,G.2001,MonthlyNotices of the Royal Astronomical Society, 328, 726 Sun,G.,Muñoz,J.B.,Mirocha,J.,&Faucher-Giguère,C.-A.2025,J.Cosmology Astropart. Phys., 2025, 034
work page 2005
- [54]
- [55]
-
[56]
Tremonti, C. A., Heckman, T. M., Kauffmann, G., et al. 2004, ApJ, 613, 898
work page 2004
-
[57]
Treu, T., Roberts-Borsani, G., Bradac, M., et al. 2022, ApJ, 935, 110
work page 2022
-
[58]
Weinberger, R., Springel, V., Hernquist, L., et al. 2017, MNRAS, 465, 3291
work page 2017
- [59]
-
[60]
2024, ApJ, 970, 19 Article number, page 11 of 14 A&A proofs:manuscript no
Zhang, Y., Ouchi, M., Nakajima, K., et al. 2024, ApJ, 970, 19 Article number, page 11 of 14 A&A proofs:manuscript no. f_esc Appendix A: Verifying the Gaussian outflow selection To evaluatethe performance ofour outflow–selection model, we applied it to a synthetic galaxy system for which the ground truth is known. We used a simplified model in which the ga...
work page 2024
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.