An official website of the United States government
Abstract Most simulations of galaxies and massive giant molecular clouds (GMCs) cannot explicitly resolve the formation (or predict the main-sequence masses) of individual stars. So they must use some prescription for the amount of feedback from an assumed population of massive stars (e.g. sampling the initial mass function, IMF). We perform a methods study of simulations of a star-forming GMC with stellar feedback from UV radiation, varying only the prescription for determining the luminosity of each stellar mass element formed (according to different IMF sampling schemes). We show that different prescriptions can lead to widely varying (factor of ∼3) star formation efficiencies (on GMC scales) even though the average mass-to-light ratios agree. Discreteness of sources is important: radiative feedback from fewer, more-luminous sources has a greater effect for a given total luminosity. These differences can dominate over other, more widely recognized differences between similar literature GMC-scale studies (e.g. numerical methods, cloud initial conditions, presence of magnetic fields). Moreover the differences in these methods are not purely numerical: some make different implicit assumptions about the nature of massive star formation, and this remains deeply uncertain in star formation theory.
more »
« less
Predicting the Radiation Field of Molecular Clouds Using Denoising Diffusion Probabilistic Models
Abstract Accurately quantifying the impact of radiation feedback in star formation is challenging. To address this complex problem, we employ deep-learning techniques known as denoising diffusion probabilistic models (DDPMs) to predict the interstellar radiation field (ISRF) strength based on three-band dust emission at 4.5, 24, and 250μm. We adopt magnetohydrodynamic simulations from the STARFORGE project that model star formation and giant molecular cloud (GMC) evolution. We generate synthetic dust emission maps matching observed spectral energy distributions in the Monoceros R2 (MonR2) GMC. We train DDPMs to estimate the ISRF using synthetic three-band dust emission. The dispersion between the predictions and true values is within a factor of 0.1 for the test set. We extended our assessment of the diffusion model to include new simulations with varying physical parameters. While there is a consistent offset observed in these out-of-distribution simulations, the model effectively constrains the relative intensity to within a factor of 2. Meanwhile, our analysis reveals a weak correlation between the ISRF solely derived from dust temperature and the actual ISRF. We apply our trained model to predict the ISRF in MonR2, revealing a correspondence between intense ISRF, bright sources, and high dust emission, confirming the model’s ability to capture ISRF variations. Our model robustly predicts radiation feedback distribution, even in complex, poorly constrained ISRF environments like those influenced by nearby star clusters. However, precise ISRF predictions require an accurate training data set mirroring the target molecular cloud’s unique physical conditions.
more »
« less
- PAR ID:
- 10533107
- Publisher / Repository:
- ApJ
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 958
- Issue:
- 1
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 97
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract The dust grain size distribution (GSD) likely varies significantly across star-forming environments in the Universe, but its impact on star formation remains unclear. This ambiguity arises because the GSD interacts nonlinearly with processes like heating, cooling, radiation, and chemistry, which have competing effects and varying environmental dependencies. Processes such as grain coagulation, expected to be efficient in dense star-forming regions, reduce the abundance of small grains and increase that of larger grains. Motivated by this, we investigate the effects of similar GSD variations on the thermochemistry and evolution of giant molecular clouds (GMCs) using magnetohydrodynamic simulations spanning a range of cloud masses and grain sizes, which explicitly incorporate the dynamics of dust grains within the full-physics framework of the STARFORGE project. We find that grain size variations significantly alter GMC thermochemistry: the leading-order effect is that larger grains, under fixed dust mass, GSD dynamic range, and dust-to-gas ratio, result in lower dust opacities. This reduced opacity permits interstellar radiation field and internal radiation photons to penetrate more deeply. This leads to rapid gas heating and inhibited star formation. Star formation efficiency is highly sensitive to grain size, with an order-of-magnitude reduction when grain size dynamic range increases from 10−3–0.1μm to 0.1–10μm. Additionally, warmer gas suppresses low-mass star formation, and decreased opacities result in a greater proportion of gas in diffuse ionized structures.more » « less
-
Abstract Stars form within dense cores composed of both gas and dust within molecular clouds. However, despite the crucial role that dust plays in the star formation process, its dynamics is frequently overlooked, with the common assumption being a constant, spatially uniform dust-to-gas ratio and grain size spectrum. In this study, we introduce a set of radiation-dust-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project. These simulations expand upon the earlier radiation MHD models, which included cooling, individual star formation, and feedback. Notably, they explicitly address the dynamics of dust grains, considering radiation, drag, and Lorentz forces acting on a diverse size spectrum of live dust grains. We find that once stars exceed a certain mass threshold (∼2M⊙), their emitted radiation can evacuate dust grains from their vicinity, giving rise to a dust-suppressed zone of size ∼100 au. This removal of dust, which interacts with gas through cooling, chemistry, drag, and radiative transfer, alters the gas properties in the region. Commencing during the early accretion stages and preceding the main-sequence phase, this process results in a mass-dependent depletion in the accreted dust-to-gas (ADG) mass ratio within both the circumstellar disk and the star. We predict that massive stars (≳10M⊙) would exhibit ADG ratios that are approximately 1 order of magnitude lower than that of their parent clouds. Consequently, stars, their disks, and circumstellar environments would display notable deviations in the abundances of elements commonly associated with dust grains, such as carbon and oxygen.more » « less
-
null (Ed.)Abstract We present a large suite of MHD simulations of turbulent, star-forming giant molecular clouds (GMCs) with stellar feedback, extending previous work by simulating 10 different random realizations for each point in the parameter space of cloud mass and size. It is found that once the clouds disperse due to stellar feedback, both self-gravitating star clusters and unbound stars generally remain, which arise from the same underlying continuum of substructured stellar density, ie. the hierarchical cluster formation scenario. The fraction of stars that are born within gravitationally-bound star clusters is related to the overall cloud star formation efficiency set by stellar feedback, but has significant scatter due to stochastic variations in the small-scale details of the star-forming gas flow. We use our numerical results to calibrate a model for mapping the bulk properties (mass, size, and metallicity) of self-gravitating GMCs onto the star cluster populations they form, expressed statistically in terms of cloud-level distributions. Synthesizing cluster catalogues from an observed GMC catalogue in M83, we find that this model predicts initial star cluster masses and sizes that are in good agreement with observations, using only standard IMF and stellar evolution models as inputs for feedback. Within our model, the ratio of the strength of gravity to stellar feedback is the key parameter setting the masses of star clusters, and of the various feedback channels direct stellar radiation (photon momentum and photoionization) is the most important on GMC scales.more » « less
-
ABSTRACT Hydrogen emission lines can provide extensive information about star-forming galaxies in both the local and high-redshift Universe. We present a detailed Lyman continuum (LyC), Lyman-α (Lyα), and Balmer line (Hα and Hβ) radiative transfer study of a high-resolution isolated Milky Way simulation using the state-of-the-art Arepo-RT radiation hydrodynamics code with the SMUGGLE galaxy formation model. The realistic framework includes stellar feedback, non-equilibrium thermochemistry accounting for molecular hydrogen, and dust grain evolution in the interstellar medium (ISM). We extend our publicly available Cosmic Lyα Transfer (COLT) code with photoionization equilibrium Monte Carlo radiative transfer and various methodology improvements for self-consistent end-to-end (non-)resonant line predictions. Accurate LyC reprocessing to recombination emission requires modelling pre-absorption by dust ($$f_\text{abs} \approx 27.5\,\rm{per\,\,cent}$$), helium ionization ($$f_\text{He} \approx 8.7\,\rm{per\,\,cent}$$), and anisotropic escape fractions ($$f_\text{esc} \approx 7.9\,\rm{per\,\,cent}$$), as these reduce the available budget for hydrogen line emission ($$f_\text{H} \approx 55.9\,\rm{per\,\,cent}$$). We investigate the role of the multiphase dusty ISM, disc geometry, gas kinematics, and star formation activity in governing the physics of emission and escape, focusing on the time variability, gas-phase structure, and spatial spectral, and viewing angle dependence of the emergent photons. Isolated disc simulations are well-suited for comprehensive observational comparisons with local Hα surveys, but would require a proper cosmological circumgalactic medium (CGM) environment as well as less dust absorption and rotational broadening to serve as analogs for high-redshift Lyα emitting galaxies. Future applications of our framework to next-generation cosmological simulations of galaxy formation including radiation-hydrodynamics that resolve ≲10 pc multiphase ISM and ≲1 kpc CGM structures will provide crucial insights and predictions for current and upcoming Lyα observations.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript
- Journal Article:
- https://doi.org/10.3847/1538-4357/acfedc
