Dust and gas modelling in radiative transfer simulations of disc-dominated galaxies with RADMC-3D
Pith reviewed 2026-05-23 07:11 UTC · model grok-4.3
The pith
A radiative transfer pipeline for galaxy simulations demonstrates that dust grain composition and size must be modeled adequately to achieve converged observables at the tens-of-percent level.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The central discovery is that the RTGen pipeline, when using an appropriate model for dust grains' composition and size, produces spectral energy distributions, continuum images, and CO luminosity maps for the simulated galaxies that match literature results from observations and theory, reaching convergence at the few tens of percent level.
What carries the argument
The RTGen pipeline, which performs Monte Carlo radiative transfer a posteriori on hydrodynamic galaxy simulations to compute dust temperatures followed by ray tracing for SEDs, images, and line profiles, with explicit modeling of dust and the atomic-to-molecular transition.
If this is right
- The pipeline predicts accurate spectral energy distributions for the studied galaxies.
- Continuum and CO luminosity images agree with results from observations and theoretical studies.
- Dust modelling has an important impact on the convergence of predicted galaxy observables.
- Adequate modelling of dust grains composition and size is required for convergence.
- The framework enables high-accuracy studies of the interstellar medium in galaxies.
Where Pith is reading between the lines
- Mock images generated this way could help interpret observations of high-redshift galaxies with JWST and ALMA.
- Similar pipelines might be tested on simulations with different resolutions or feedback prescriptions to assess robustness.
- The method could extend to studying the impact on other emission lines beyond CO.
- Applying the pipeline to cosmological simulations rather than isolated galaxies could reveal environmental effects on observables.
Load-bearing premise
That the selected dust abundance, composition, grain size distribution, and atomic-to-molecular transition model accurately represent the physical conditions in the simulated galaxies.
What would settle it
Direct comparison of the pipeline's predicted SEDs and images against multi-wavelength observations of real disc galaxies with similar properties, checking if residuals exceed a few tens of percent.
Figures
read the original abstract
Bridging theory and observations is a key task to understand galaxy formation and evolution. With the advent of state-of-the-art observational facilities, an accurate modelling of galaxy observables through radiative transfer simulations coupled to hydrodynamic simulations of galaxy formation must be performed. We present a novel pipeline, dubbed RTGen, based on the Monte Carlo radiative transfer code RADMC-3D , and explore the impact of the physical assumptions and modelling of dust and gas phases on the resulting galaxy observables. In particular, we address the impact of the dust abundance, composition, and grain size, as well as model the atomic-to-molecular transition and study the resulting emission from molecular gas. We apply Monte Carlo radiative transfer a posteriori to determine the dust temperature in six different hydrodynamic simulations of isolated galaxies. Afterwards, we apply ray tracing to compute the spectral energy distribution, as well as continuum images and spectral line profiles. We find our pipeline to predict accurate spectral energy distribution distributions of the studied galaxies, as well as continuum and CO luminosity images, in good agreement with literature results from both observations and theoretical studies. In particular, we find the dust modelling to have an important impact on the convergence of the resulting predicted galaxy observables, and that an adequate modelling of dust grains composition and size is required. We conclude that our novel framework is ready to perform high-accuracy studies of the observables of the ISM, reaching few tens percent convergence under the studied baseline configuration. This will enable robust studies of galaxy formation, and in particular of the nature of massive clumps in high-redshift galaxies, through the generation of mock images mimicking observations from state-of-the-art facilities such as JWST and ALMA.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces the RTGen pipeline, which couples the RADMC-3D Monte Carlo radiative transfer code to hydrodynamic simulations of six isolated disc-dominated galaxies. It computes dust temperatures, spectral energy distributions (SEDs), continuum images, and CO luminosity images and line profiles. The authors explore the effects of dust abundance, composition, grain size distribution, and the atomic-to-molecular gas transition, concluding that their pipeline produces observables in good agreement with literature and achieves convergence to within a few tens of percent under a baseline dust and gas modeling configuration.
Significance. If the reported convergence and agreement with observations and theory are robust, this work would offer a practical framework for generating realistic mock observations from galaxy formation simulations. This is particularly relevant for interpreting data from facilities like JWST and ALMA, and for studying the interstellar medium in high-redshift galaxies. The focus on the sensitivity to dust modeling choices is a positive aspect, provided it is backed by quantitative tests.
major comments (2)
- [Abstract] Abstract: The central claim that the pipeline 'reaches few tens percent convergence under the studied baseline configuration' is not accompanied by any quantitative metrics, error bars, or explicit description of the convergence measurement procedure across the six simulations. This is load-bearing for the headline result regarding the pipeline's readiness for high-accuracy studies.
- [Abstract] Abstract: The assertion that 'an adequate modelling of dust grains composition and size is required' for convergence is stated without reporting sensitivity tests to alternative grain size distributions, compositions, dust abundances, or atomic-to-molecular transition thresholds. Without such tests, it remains unclear whether the reported convergence is general or specific to the chosen baseline parameters.
Simulated Author's Rebuttal
We thank the referee for their constructive feedback on the abstract. We agree that the abstract requires strengthening with explicit references to quantitative results and sensitivity tests already present in the main text. We will revise the abstract accordingly while preserving its concise nature.
read point-by-point responses
-
Referee: [Abstract] The central claim that the pipeline 'reaches few tens percent convergence under the studied baseline configuration' is not accompanied by any quantitative metrics, error bars, or explicit description of the convergence measurement procedure across the six simulations.
Authors: We acknowledge the abstract does not embed the full quantitative details. The main manuscript (Sections 4 and 5) reports direct comparisons of SEDs, continuum images, and CO line profiles across the six simulations, showing typical differences of 20-40% between baseline and varied configurations, with explicit convergence metrics derived from pixel-by-pixel and integrated flux ratios. We will revise the abstract to state 'reaching convergence to within ~30% for key observables under the baseline configuration, as quantified in Sections 4-5' and add a parenthetical reference to the convergence procedure. revision: yes
-
Referee: [Abstract] The assertion that 'an adequate modelling of dust grains composition and size is required' for convergence is stated without reporting sensitivity tests to alternative grain size distributions, compositions, dust abundances, or atomic-to-molecular transition thresholds.
Authors: The manuscript does perform and report these sensitivity tests: Section 3.2 varies dust-to-gas ratio, composition (silicate vs. graphite fractions), grain size distributions (MRN vs. alternative power laws), and H2 formation thresholds, with results in Figures 6-9 showing that non-baseline choices increase discrepancies by factors of 1.5-3. The abstract will be revised to read 'an adequate modelling of dust grain composition and size is required, as demonstrated by sensitivity tests in Section 3.2' to make this explicit. revision: yes
Circularity Check
No significant circularity; pipeline applies external RT code to independent hydro simulations
full rationale
The paper describes applying the external RADMC-3D Monte Carlo code a posteriori to six independent hydrodynamic simulations of isolated galaxies, then computing SEDs, images, and line profiles via ray tracing. Outputs are compared to external literature (observations and other theoretical studies). No equations or claims reduce by construction to self-defined quantities, fitted inputs renamed as predictions, or load-bearing self-citations. The statement that dust modeling impacts convergence is an empirical observation from the runs under the chosen baseline parameters, not a definitional tautology. The physical adequacy of the dust and transition models is an assumption, but the derivation chain itself does not collapse into its inputs.
Axiom & Free-Parameter Ledger
free parameters (2)
- dust abundance
- dust composition and grain size distribution
axioms (2)
- domain assumption Standard Monte Carlo radiative transfer assumptions in RADMC-3D are sufficient to compute dust temperatures and emission from the hydrodynamic outputs.
- domain assumption The atomic-to-molecular transition model used produces realistic molecular gas distributions for the line emission calculation.
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We adopt the scaling relation log10(DGR)=2.445 log10(Z/Z⊙)−2.029 … C/O=0.8 … n(a)∼a^−3.5 … β_CO fitted to ~1 % tolerance
-
IndisputableMonolith/Foundation/AlexanderDuality.leanalexander_duality_circle_linking unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We apply Monte Carlo radiative transfer a posteriori … reaching few tens percent convergence under the studied baseline configuration
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]
, " * write output.state after.block = add.period write newline
ENTRY address archiveprefix author booktitle chapter edition editor howpublished institution eprint journal key month note number organization pages publisher school series title type volume year label extra.label sort.label short.list INTEGERS output.state before.all mid.sentence after.sentence after.block FUNCTION init.state.consts #0 'before.all := #1 ...
-
[2]
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 global.max substring 't := if while FUNCTION word.in bbl.in " " * FUNCTION format....
- [3]
- [4]
- [5]
- [6]
-
[7]
Bjorkman , J. E. & Wood , K. 2001, , 554, 615
work page 2001
- [8]
- [9]
-
[10]
Bohren , C. F. & Huffman , D. R. 1983, Absorption and scattering of light by small particles , Research supported by the University of Arizona and Institute of Occupational and Environmental Health. New York, Wiley-Interscience, 1983, 541 p
work page 1983
-
[11]
Bolatto , A. D., Wolfire , M., & Leroy , A. K. 2013, , 51, 207
work page 2013
-
[12]
A., van der Werf , P., Weiss , A., et al
Boogaard , L. A., van der Werf , P., Weiss , A., et al. 2020, , 902, 109
work page 2020
-
[13]
Bothwell , M. S., Smail , I., Chapman , S. C., et al. 2013, , 429, 3047
work page 2013
-
[14]
Bressan , A., Fagotto , F., Bertelli , G., & Chiosi , C. 1993, , 100, 647
work page 1993
- [15]
- [16]
- [17]
- [18]
-
[19]
R., Bovino , S., Lupi , A., Schleicher , D
Capelo , P. R., Bovino , S., Lupi , A., Schleicher , D. R. G., & Grassi , T. 2018, , 475, 3283
work page 2018
-
[20]
CASA Team , Bean , B., Bhatnagar , S., et al. 2022, , 134, 114501
work page 2022
-
[21]
Casey , C. M. 2012, , 425, 3094
work page 2012
-
[22]
M., Narayanan , D., & Cooray , A
Casey , C. M., Narayanan , D., & Cooray , A. 2014, , 541, 45
work page 2014
- [23]
- [24]
-
[25]
Dav \'e , R., Angl \'e s-Alc \'a zar , D., Narayanan , D., et al. 2019, , 486, 2827
work page 2019
-
[26]
Dav \'e , R., Thompson , R., & Hopkins , P. F. 2016, , 462, 3265
work page 2016
- [27]
- [28]
-
[29]
Dessauges-Zavadsky , M., Richard , J., Combes , F., et al. 2023, , 519, 6222
work page 2023
-
[30]
2019, Nature Astronomy, 3, 1115
Dessauges-Zavadsky , M., Richard , J., Combes , F., et al. 2019, Nature Astronomy, 3, 1115
work page 2019
-
[31]
Dessauges-Zavadsky , M., Zamojski , M., Schaerer , D., et al. 2015, , 577, A50
work page 2015
-
[32]
Dewdney , P. E., Hall , P. J., Schilizzi , R. T., & Lazio , T. J. L. W. 2009, IEEE Proceedings, 97, 1482
work page 2009
-
[33]
Diemer , B., Stevens , A. R. H., Forbes , J. C., et al. 2018, , 238, 33
work page 2018
-
[34]
Draine , B. T. 1978, , 36, 595
work page 1978
-
[35]
Draine , B. T. 2003, , 41, 241
work page 2003
-
[36]
Draine , B. T. & Bertoldi , F. 1996, , 468, 269
work page 1996
-
[37]
Draine , B. T. & Lee , H. M. 1984, , 285, 89
work page 1984
-
[38]
P., Juhasz , A., Pohl , A., et al
Dullemond , C. P., Juhasz , A., Pohl , A., et al. 2012, RADMC-3D: A multi-purpose radiative transfer tool , Astrophysics Source Code Library, record ascl:1202.015
work page 2012
-
[39]
Esmerian , C. J. & Gnedin , N. Y. 2022, , 940, 74
work page 2022
- [40]
- [41]
-
[42]
J., Chatzikos , M., Guzm \'a n , F., et al
Ferland , G. J., Chatzikos , M., Guzm \'a n , F., et al. 2017, , 53, 385
work page 2017
-
[43]
Ferland , G. J., Korista , K. T., Verner , D. A., et al. 1998, , 110, 761
work page 1998
-
[44]
Ferrarotti , A. S. & Gail , H. P. 2006, , 447, 553
work page 2006
- [45]
- [46]
-
[47]
Gardner , J. P., Mather , J. C., Clampin , M., et al. 2006, , 123, 485
work page 2006
-
[48]
Gnedin , N. Y. & Abel , T. 2001, , 6, 437
work page 2001
-
[49]
Gnedin , N. Y. & Draine , B. T. 2014, , 795, 37
work page 2014
-
[50]
Gnedin , N. Y. & Kravtsov , A. V. 2010, , 714, 287
work page 2010
-
[51]
Gnedin , N. Y. & Kravtsov , A. V. 2011, , 728, 88
work page 2011
- [52]
-
[53]
Y., Tassis , K., & Kravtsov , A
Gnedin , N. Y., Tassis , K., & Kravtsov , A. V. 2009, , 697, 55
work page 2009
- [54]
-
[55]
J., Abergel , A., Abreu , A., et al
Griffin , M. J., Abergel , A., Abreu , A., et al. 2010, , 518, L3
work page 2010
-
[56]
Henyey , L. G. & Greenstein , J. L. 1941, , 93, 70
work page 1941
- [57]
-
[58]
Hirashita , H. & Kobayashi , H. 2013, Earth, Planets and Space, 65, 1083
work page 2013
-
[59]
Hopkins , P. F. & Lee , H. 2016, , 456, 4174
work page 2016
-
[60]
Khatri , P., Porciani , C., Romano-D \' az , E., Seifried , D., & Sch \"a be , A. 2024, , 688, A194
work page 2024
- [61]
-
[62]
Krumholz , M. R. 2013, , 436, 2747
work page 2013
-
[63]
Krumholz , M. R. 2014, , 437, 1662
work page 2014
- [64]
-
[65]
K., Walter , F., Brinks , E., et al
Leroy , A. K., Walter , F., Brinks , E., et al. 2008, , 136, 2782
work page 2008
- [66]
-
[67]
F., Yajima , H., Zhu , Q., & Maji , M
Li , Y., Gu , M. F., Yajima , H., Zhu , Q., & Maji , M. 2020, , 494, 1919
work page 2020
- [68]
-
[69]
Lucy , L. B. 1999, , 344, 282
work page 1999
-
[70]
Ludlow , A. D., Navarro , J. F., Angulo , R. E., et al. 2014, , 441, 378
work page 2014
-
[71]
R., Volonteri , M., & Silk , J
Lupi , A., Bovino , S., Capelo , P. R., Volonteri , M., & Silk , J. 2018, , 474, 2884
work page 2018
- [72]
- [73]
-
[74]
S., Rumpl , W., & Nordsieck , K
Mathis , J. S., Rumpl , W., & Nordsieck , K. H. 1977, , 217, 425
work page 1977
- [75]
-
[76]
Mayer , L., Moore , B., Quinn , T., Governato , F., & Stadel , J. 2002, , 336, 119
work page 2002
- [77]
-
[78]
Nagy , D., Dessauges-Zavadsky , M., Richard , J., et al. 2022, , 657, A25
work page 2022
-
[79]
Narayanan , D., Turk , M. J., Robitaille , T., et al. 2021, , 252, 12
work page 2021
- [80]
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.