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 Observations indicate dust populations vary between galaxies and within them, suggesting a complex life cycle and evolutionary history. Here we investigate the evolution of galactic dust populations across cosmic time using a suite of cosmological zoom-in simulations from the Feedback in Realistic Environments project, spanning $$M_{\rm vir}=10^{9-12}{M}_{\odot };\, M_{*}=10^{6-11}\, {M}_{\odot }$$. Our simulations incorporate a dust evolution model that accounts for the dominant sources of dust production, growth, and destruction and follows the evolution of specific dust species. All galactic dust populations in our suite exhibit similar evolutionary histories, with gas–dust accretion being the dominant producer of dust mass for all but the most metal-poor galaxies. Similar to previous works, we find the onset of efficient gas–dust accretion occurs above a ‘critical’ metallicity threshold (Zcrit). Due to this threshold, our simulations reproduce observed trends between galactic D/Z and metallicity and element depletion trends in the interstellar medium. However, we find Zcrit varies between dust species due to differences in key element abundances, dust physical properties, and life cycle processes resulting in $$Z_{\rm crit}\sim 0.05{\rm Z}_{\odot },\, 0.2{\rm Z}_{\odot },\, 0.5{\rm Z}_{\odot }$$ for metallic iron, silicates, and carbonaceous dust, respectively. These variations could explain the lack of small carbonaceous grains observed in the Magellanic Clouds. We also find a delay between the onset of gas–dust accretion and when a dust population reaches equilibrium, which we call the equilibrium time-scale (τequil). The relation between τequil and the metal enrichment time-scale of a galaxy, determined by its recent evolutionary history, can contribute to the scatter in the observed relation between galactic D/Z and metallicity. 

Context. Carbon monoxide (CO) is a poor tracer of H₂in the diffuse interstellar medium (ISM), where most of the carbon is not incorporated into CO molecules, unlike the situation at higher extinctions. Aims. We present a novel, indirect method for constraining H₂column densities (N_H2) without employing CO observations. We show that previously recognized nonlinearities in the relation between the extinction,A_V(H₂), derived from dust emission and the H Icolumn density (N_H I) are due to the presence of molecular gas. Methods. We employed archival (N_H2) data, obtained from the UV spectra of stars, and calculatedA_V(H₂) toward these sight lines using 3D extinction maps. The following relation fits the data: logN_H2= 1.38742 (logA_V(H₂)⁾³− 0.05359 (logA_V(H₂)⁾²+ 0.25722 logA_V(H₂) + 20.67191. This relation is useful for constrainingN_H2in the diffuse ISM as it requires onlyN_H Iand dust extinction data, which are both easily accessible. In 95% of the cases, the estimates produced by the fitted equation have deviations of less than a factor of 3.5. We constructed aN_H2map of our Galaxy and compared it to the CO integrated intensity (W_CO) distribution. Results. We find that the average ratio (X_CO) betweenN_H2andW_COis approximately equal to 2 × 10²⁰cm⁻²(K km s⁻¹)⁻¹, consistent with previous estimates. However, we find that theX_COfactor varies by orders of magnitude on arcminute scales between the outer and the central portions of molecular clouds. For regions withN_H2≳ 10²⁰cm⁻², we estimate that the average H₂fractional abundance,f_H2= 2N_H2/(2N_H2+N_H I), is 0.25. Multiple (distinct) largely atomic clouds are likely found along high-extinction sightlines (A_V≥ 1 mag), hence limitingf_H2in these directions. Conclusions. More than 50% of the lines of sight withN_H2≥ 10²⁰cm⁻²are untraceable by CO with aJ= 1−0 sensitivity limitW_CO= 1 K km s⁻¹. 

Recently, we demonstrated self-consistent formation of strongly-magnetized quasar accretion disks (QADs) from cosmological radiation-magnetohydrodynamic-thermochemical galaxy-star formation simulations, including the full STARFORGE physics shown previously to produce a reasonable IMF under typical ISM conditions. Here we study star formation and the stellar IMF in QADs, on scales from 100 au to 10 pc from the SMBH. We show it is critical to include physics often previously neglected, including magnetic fields, radiation, and (proto)stellar feedback. Closer to the SMBH, star formation is suppressed, but the (rare) stars that do form exhibit top-heavy IMFs. Stars can form only in special locations (e.g. magnetic field switches) in the outer QAD. Protostars accrete their natal cores rapidly but then dynamically decouple from the gas and ‘wander,’ ceasing accretion on timescales ~100 yr. Their jets control initial core accretion, but the ejecta are ‘swept up’ into the larger-scale QAD flow without much dynamical effect. The strong tidal environment strongly suppresses common-core multiplicity. The IMF shape depends sensitively on un-resolved dynamics of protostellar disks (PSDs), as the global dynamical times can become incredibly short (< yr) and tidal fields are incredibly strong, so whether PSDs can efficiently transport angular momentum or fragment catastrophically at <10 au scales requires novel PSD simulations to properly address. Most analytic IMF models and analogies with planet formation in PSDs fail qualitatively to explain the simulation IMFs, though we discuss a couple of viable models. 

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 Negative feedback from accreting supermassive black holes is considered crucial in suppressing star formation and quenching massive galaxies. However, several models and observations suggest that black hole feedback may have a positive effect, triggering star formation by compressing interstellar medium gas to higher densities. We investigate the dual role of black hole feedback using cosmological hydrodynamic simulations from the Feedback In Realistic Environment (FIRE) project, incorporating a novel implementation of hyper-refined accretion-disc winds. Focusing on a massive, star-forming galaxy at z ∼ 2 ($$M_{\rm halo} \sim 10^{12.5}\, {\rm M}_{\odot }$$), we demonstrate that strong quasar winds with a kinetic power of ∼1046 erg s−1, persisting for over 20 Myr, drive the formation of a central gas cavity and significantly reduce the surface density of star formation across the galaxy’s disc. The suppression of star formation primarily occurs by limiting the availability of gas for star formation rather than by evacuating the pre-existing star-forming gas reservoir (preventive feedback dominates over ejective feedback). Despite the overall negative impact of quasar winds, we identify several potential indicators of local positive feedback, including (1) the spatial anticorrelation between wind-dominated regions and star-forming clumps, (2) higher local star formation efficiency in compressed gas at the edge of the cavity, and (3) increased contribution of outflowing material to local star formation. Moreover, stars formed under the influence of quasar winds tend to be located at larger radial distances. Our findings suggest that both positive and negative AGN feedback can coexist within galaxies, although the local positive triggering of star formation has a minor influence on global galaxy growth. 

ABSTRACT Without active galactic nucleus (AGN) feedback, simulated massive, star-forming galaxies become too compact relative to observed galaxies at z ≲ 2. In this paper, we perform high-resolution re-simulations of a massive ($$M_{\star }\sim 10^{11}\, \rm {{\rm M}_{\odot }}$$) galaxy at z ∼ 2.3, drawn from the Feedback in Realistic Environments (FIRE) project. In the simulation without AGN feedback, the galaxy experiences a rapid starburst and shrinking of its half-mass radius. We experiment with driving mechanical AGN winds, using a state-of-the-art hyper-Lagrangian refinement technique to increase particle resolution. These winds reduce the gas surface density in the inner regions of the galaxy, suppressing the compact starburst and maintaining an approximately constant half-mass radius. Using radiative transfer, we study the impact of AGN feedback on the magnitude and extent of the multiwavelength continuum emission. When AGN winds are included, the suppression of the compact, dusty starburst results in lowered flux at FIR wavelengths (due to decreased star formation) but increased flux at optical-to-near-IR wavelengths (due to decreased dust attenuation, in spite of the lowered star formation rate), relative to the case without AGN winds. The FIR half-light radius decreases from ∼1 to $$\sim 0.1\, \rm {kpc}$$ in $$\lesssim 40\, \rm {Myr}$$ when AGN winds are not included, but increases to $$\sim 2\, \rm {kpc}$$ when they are. Interestingly, the half-light radius at optical-NIR wavelengths remains approximately constant over $$35\, \rm {Myr}$$, for simulations with and without AGN winds. In the case without winds, this occurs despite the rapid compaction, and is due to heavy dust obscuration in the inner regions of the galaxy. This work highlights the importance of forward-modelling when comparing simulated and observed galaxy populations. 

ABSTRACT Recent theoretical studies predict that the circumgalactic medium (CGM) around low-redshift, ∼L* galaxies could have substantial non-thermal pressure support in the form of cosmic rays. However, these predictions are sensitive to the specific model of cosmic ray transport employed, which is theoretically and observationally underconstrained. In this work, we propose a novel observational constraint for calculating the lower limit of the radially averaged, effective cosmic ray transport rate, $${\kappa _{\rm eff}^{\rm min}}$$. Under a wide range of assumptions (so long as cosmic rays do not lose a significant fraction of their energy in the galactic disc, regardless of whether the cosmic ray pressure is important or not in the CGM), we demonstrate a well-defined relationship between $${\kappa _{\rm eff}^{\rm min}}$$ and three observable galaxy properties: the total hydrogen column density, the average star formation rate, and the gas circular velocity. We use a suite of Feedback in Realistic Environments 2 galaxy simulations with a variety of cosmic ray transport physics to demonstrate that our analytical model of $${\kappa _{\rm eff}^{\rm min}}$$ is a robust lower limit of the true cosmic ray transport rate. We then apply our new model to calculate $${\kappa _{\rm eff}^{\rm min}}$$ for galaxies in the COS-Halos sample, and confirm this already reveals strong evidence for an effective transport rate that rises rapidly away from the interstellar medium to values $${\kappa _{\rm eff}^{\rm min}}\gtrsim 10^{30\!-\!31}\, {\rm cm}^2\, {\rm s}^{-1}$$ (corresponding to anisotropic streaming velocities of $$v^{\rm stream}_{\rm eff} \gtrsim 1000\, {\rm km}\, {\rm s}^{-1}$$) in the diffuse CGM, at impact parameters larger than 50–100 kpc. We discuss how future observations can provide qualitatively new constraints in our understanding of cosmic rays in the CGM and intergalactic medium. 

ABSTRACT Most observed stars are part of a multiple star system, but the formation of such systems and the role of environment and various physical processes is still poorly understood. We present a suite of radiation-magnetohydrodynamic simulations of star-forming molecular clouds from the STARFORGE project that include stellar feedback with varied initial surface density, magnetic fields, level of turbulence, metallicity, interstellar radiation field, simulation geometry and turbulent driving. In our fiducial cloud, the raw simulation data reproduces the observed multiplicity fractions for Solar-type and higher mass stars, similar to previous works. However, after correcting for observational incompleteness the simulation underpredicts these values. The discrepancy is likely due to the lack of disc fragmentation, as the simulation only resolves multiples that form either through capture or core fragmentation. The raw mass distribution of companions is consistent with randomly drawing from the initial mass function for the companions of $$\gt 1\, \mathrm{M}_{\rm \odot }$$ stars. However, accounting for observational incompleteness produces a flatter distribution similar to observations. We show that stellar multiplicity changes as the cloud evolves and anticorrelates with stellar density. This relationship also explains most multiplicity variations between runs, i.e. variations in the initial conditions that increase stellar density (increased surface density, reduced turbulence) also act to decrease multiplicity. While other parameters, such as metallicity, interstellar radiation, and geometry significantly affect the star formation history or the IMF, varying them produces no clear trend in stellar multiplicity properties.