Search for: All records

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⁻³–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. 

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. 

Abstract Partial dust obscuration in active galactic nuclei (AGN) has been proposed as a potential explanation for some cases of AGN variability. The dust-gas mixture present in AGN tori is accelerated by radiation pressure, leading to the launching of an AGN wind. Dust under these conditions has been shown to be unstable to a generic class of fast-growing resonant drag instabilities (RDIs). In this work, we present the first numerical simulations of radiation-driven outflows that explicitly include dust dynamics in conditions resembling AGN winds. We investigate the implications of RDIs on the torus morphology, AGN variability, and the ability of radiation to effectively launch a wind. We find that the RDIs rapidly develop, reaching saturation at times much shorter than the global timescales of the outflows, resulting in the formation of filamentary structure on box-size scales with strong dust clumping and super-Alfvénic velocity dispersions. The instabilities lead to fluctuations in dust opacity and gas column density of 10-20% when integrated along mock observed lines-of-sight to the quasar accretion disk. These fluctuations occur over year to decade timescales and exhibit a red-noise power spectrum commonly observed for AGN. Additionally, we find that the radiation effectively couples with the dust-gas mixture, launching highly supersonic winds that entrain 70-90% of the gas, with a factor of ≲ 3 photon momentum loss relative to the predicted multiple-scattering momentum loading rate. Therefore, our findings suggest that RDIs play an important role in driving the clumpy nature of AGN tori and generating AGN variability consistent with observations. 

Abstract Radiation-dust driven outflows, where radiation pressure on dust grains accelerates gas, occur in many astrophysical environments. Almost all previous numerical studies of these systems have assumed that the dust was perfectly-coupled to the gas. However, it has recently been shown that the dust in these systems is unstable to a large class of ‘resonant drag instabilities’ (RDIs) which de-couple the dust and gas dynamics and could qualitatively change the nonlinear outcome of these outflows. We present the first simulations of radiation-dust driven outflows in stratified, inhomogeneous media, including explicit grain dynamics and a realistic spectrum of grain sizes and charge, magnetic fields and Lorentz forces on grains (which dramatically enhance the RDIs), Coulomb and Epstein drag forces, and explicit radiation transport allowing for different grain absorption and scattering properties. In this paper we consider conditions resembling giant molecular clouds (GMCs), H ii regions, and distributed starbursts, where optical depths are modest (≲ 1), single-scattering effects dominate radiation-dust coupling, Lorentz forces dominate over drag on grains, and the fastest-growing RDIs are similar, such as magnetosonic and fast-gyro RDIs. These RDIs generically produce strong size-dependent dust clustering, growing nonlinear on timescales that are much shorter than the characteristic times of the outflow. The instabilities produce filamentary and plume-like or ‘horsehead’ nebular morphologies that are remarkably similar to observed dust structures in GMCs and H ii regions. Additionally, in some cases they strongly alter the magnetic field structure and topology relative to filaments. Despite driving strong micro-scale dust clumping which leaves some gas ‘behind,’ an order-unity fraction of the gas is always efficiently entrained by dust.