-
Thermodynamics of Giant Molecular Clouds: The Effects of Dust Grain Size
Authors:
Nadine H. Soliman,
Philip F. Hopkins,
Michael Y. Grudić
Abstract:
The dust grain size distribution (GSD) likely varies significantly across different star-forming environments in the Universe, but the overall impact of this variation on star formation remains unclear. This ambiguity arises because the GSD interacts non-linearly with processes like heating/cooling, radiation, and chemistry, which have competing effects and different environmental dependencies. In…
▽ More
The dust grain size distribution (GSD) likely varies significantly across different star-forming environments in the Universe, but the overall impact of this variation on star formation remains unclear. This ambiguity arises because the GSD interacts non-linearly with processes like heating/cooling, radiation, and chemistry, which have competing effects and different environmental dependencies. In this study, we investigate the effects of GSD variation on the thermochemistry and evolution of giant molecular clouds (GMCs). To achieve this, we conducted radiation-dust-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 differences in grain size significantly alter the thermochemistry of GMCs. Specifically, we show that the leading-order effect is that larger grains, under fixed dust mass and dust-to-gas ratio conditions, result in lower dust opacities. This reduced opacity permits ISRF photons to penetrate more deeply and allows internal radiation field photons to permeate more extensively into the cloud, resulting in rapid gas heating and the inhibition of star formation. We find that star formation efficiency is highly sensitive to grain size, with an order of magnitude reduction in efficiency when grain size increases from 0.1 $\rmμm$ to 10 $\rmμm$. Additionally, we note that warmer gas suppresses the formation of low-mass stars. Moreover, as a consequence of the decreased opacities, we observe a greater proportion of gas residing in diffuse ionized structures.
△ Less
Submitted 12 July, 2024;
originally announced July 2024.
-
Dust-Evacuated Zones Near Massive Stars: Consequences of Dust Dynamics on Star-forming Regions
Authors:
Nadine H. Soliman,
Philip F. Hopkins,
Michael Y. Grudić
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 st…
▽ More
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 {\small 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 interactions between radiation and dust significantly influence the properties of gas surrounding and accreting onto massive stars. Specifically, we find that once stars exceed a certain mass threshold ($\sim 2 M_{\odot}$), their emitted radiation can evacuate dust grains from their vicinity, giving rise to a dust-suppressed zone of size $\sim 100$ AU. 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 disc and the star. We predict massive stars ($\gtrsim 10 M_{\odot}$) would exhibit ADG ratios that are approximately one order of magnitude lower than that of their parent clouds. Consequently, stars, their discs, and circumstellar environments would display notable deviations in the abundances of elements commonly associated with dust grains, such as carbon and oxygen.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Microphysical Regulation of Non-Ideal MHD in Weakly-Ionized Systems: Does the Hall Effect Matter?
Authors:
Philip F. Hopkins,
Jonathan Squire,
Raphael Skalidis,
Nadine H. Soliman
Abstract:
The magnetohydrodynamics (MHD) equations plus 'non-ideal' (Ohmic, Hall, ambipolar) resistivities are widely used to model weakly-ionized astrophysical systems. We show that if gradients in the magnetic field become too steep, the implied charge drift speeds become much faster than microphysical signal speeds, invalidating the assumptions used to derive both the resistivities and MHD equations them…
▽ More
The magnetohydrodynamics (MHD) equations plus 'non-ideal' (Ohmic, Hall, ambipolar) resistivities are widely used to model weakly-ionized astrophysical systems. We show that if gradients in the magnetic field become too steep, the implied charge drift speeds become much faster than microphysical signal speeds, invalidating the assumptions used to derive both the resistivities and MHD equations themselves. Generically this situation will excite microscale instabilities that suppress the drift and current. We show this could be relevant at low ionization fractions especially if Hall terms appear significant, external forces induce supersonic motions, or dust grains become a dominant charge carrier. Considering well-established treatments of super-thermal drifts in laboratory, terrestrial, and Solar plasmas as well as conduction and viscosity models, we generalize a simple prescription to rectify these issues, where the resistivities are multiplied by a correction factor that depends only on already-known macroscopic quantities. This is generalized for multi-species and weakly-ionized systems, and leaves the equations unchanged in the drift limits for which they are derived, but restores physical behavior (driving the system back towards slow drift by diffusing away small-scale gradients in the magnetic field) if the limits are violated. This has important consequences: restoring intuitive behaviors such as the system becoming hydrodynamic in the limit of zero ionization; suppressing magnetic structure on scales below a critical length which can comparable to circumstellar disk sizes; limiting the maximum magnetic amplification; and suppressing the effects of the Hall term in particular. This likely implies that the Hall term does not become dynamically important under most conditions of interest in these systems.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Co-Evolution vs. Co-existence: The Effect of Accretion Modelling on the Evolution of Black Holes and Host Galaxies
Authors:
Nadine H. Soliman,
Andrea V. Macciò,
Marvin Blank
Abstract:
We append two additional black hole (BH) accretion models, namely viscous disc and gravitational torque-driven accretion, into the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project of galaxy simulations. We show that these accretion models, characterized by a weaker dependence on the BH mass compared to the commonly used Bondi-Hoyle accretion, naturally create a common evo…
▽ More
We append two additional black hole (BH) accretion models, namely viscous disc and gravitational torque-driven accretion, into the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project of galaxy simulations. We show that these accretion models, characterized by a weaker dependence on the BH mass compared to the commonly used Bondi-Hoyle accretion, naturally create a common evolutionary track (co-existence) between the mass of the BH and the stellar mass of the galaxy, even without any direct coupling via feedback (FB). While FB is indeed required to control the final BH and stellar mass of the galaxies, our results suggest that FB might not be the leading driver of the cosmic co-evolution between these two quantities; in these models, co-evolution is simply determined by the shared central gas supply. Conversely, simulations using Bondi-Hoyle accretion show a two-step evolution, with an early growth of stellar mass followed by exponential growth of the central supermassive black hole (SMBH). Our results show that the modelling of BH accretion (sometimes overlooked) is an extremely important part of BH evolution and can improve our understanding of how scaling relations emerge and evolve, and whether SMBH and stellar mass co-exist or co-evolve through cosmic time.
△ Less
Submitted 25 July, 2023;
originally announced July 2023.
-
Dust Dynamics in AGN Winds: A New Mechanism For Multiwavelength AGN Variability
Authors:
Nadine H. Soliman,
Philip F. Hopkins
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 torii 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). We present the…
▽ More
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 torii 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). We present the first set of numerical simulations of radiation driven outflows that include explicit dust dynamics in conditions resembling AGN winds and discuss the implications of the RDIs on the morphology of the AGN torus, AGN variability, and the ability of the 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 on micro-scales. This results in 10-20% fluctuations in dust opacity and gas column density, integrated along mock observed lines-of-sight to the quasar accretion disk, over year to decade timescales with a red-noise power spectrum which is commonly observed for AGN. Additionally, all our simulations show that the radiation is sufficiently coupled to the dust-gas mixture launching highly super-sonic winds, which entrain 70-90% of gas, with a factor of $\lesssim 3$ photon momentum loss relative to the ideal case. Therefore, the RDIs could play an important role in generating the clumpy nature of AGN torii and driving AGN variability consistent with observations.
△ Less
Submitted 24 October, 2022;
originally announced October 2022.
-
NIHAO XXVII: Crossing the green valley
Authors:
Marvin Blank,
Andrea V. Macciò,
Xi Kang,
Keri L. Dixon,
Nadine H. Soliman
Abstract:
The transition of high-mass galaxies from being blue and star forming to being red and dead is a crucial step in galaxy evolution, yet not fully understood. In this work, we use the NIHAO suite of galaxy simulations to investigate the relation between the transition time through the green valley and other galaxy properties. The typical green valley crossing time of our galaxies is approximately 40…
▽ More
The transition of high-mass galaxies from being blue and star forming to being red and dead is a crucial step in galaxy evolution, yet not fully understood. In this work, we use the NIHAO suite of galaxy simulations to investigate the relation between the transition time through the green valley and other galaxy properties. The typical green valley crossing time of our galaxies is approximately 400 Myr, somewhat shorter than observational estimates. The crossing of the green valley is triggered by the onset of AGN feedback and the subsequent shut down of star formation. Interestingly the time spent in the green valley is not related to any other galaxy properties, such as stellar age or metallicity, or the time at which the star formation quenching takes place. The crossing time is set by two main contributions: the ageing of the current stellar population and the residual star formation in the green valley. These effects are of comparable magnitude, while major and minor mergers have a negligible contribution. Most interestingly, we find the time that a galaxy spends to travel through the green valley is twice the $e$-folding time of the star formation quenching. This result is stable against galaxy properties and the exact numerical implementation of AGN feedback in the simulation. Assuming a typical crossing time of about one Gyr inferred from observations, our results imply that any mechanism or process aiming to quench star formation, must do it on a typical timescale of 500 Myr.
△ Less
Submitted 25 April, 2022;
originally announced April 2022.
-
Dust in the Wind with Resonant Drag Instabilities: I. The Dynamics of Dust-Driven Outflows in GMCs and HII Regions
Authors:
Philip F. Hopkins,
Anna L. Rosen,
Jonathan Squire,
Georgia V. Panopoulou,
Nadine H. Soliman,
Darryl Seligman,
Ulrich P. Steinwandel
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…
▽ More
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), HII regions, and distributed starbursts, where optical depths are modest, 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 HII 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.
△ Less
Submitted 8 November, 2022; v1 submitted 9 July, 2021;
originally announced July 2021.
-
NIHAO XXVI: Nature versus nurture, the Star Formation Main Sequence and the origin of its scatter
Authors:
Marvin Blank,
Liam E. Meier,
Andrea V. Macciò,
Aaron A. Dutton,
Keri L. Dixon,
Nadine H. Soliman,
Xi Kang
Abstract:
We investigate how the NIHAO galaxies match the observed star formation main sequence (SFMS) and what the origin of its scatter is. The NIHAO galaxies reproduce the SFMS and generally agree with observations, but the slope is about unity and thus significantly larger than observed values. This is because observed galaxies at large stellar masses, although still being part of the SFMS, are already…
▽ More
We investigate how the NIHAO galaxies match the observed star formation main sequence (SFMS) and what the origin of its scatter is. The NIHAO galaxies reproduce the SFMS and generally agree with observations, but the slope is about unity and thus significantly larger than observed values. This is because observed galaxies at large stellar masses, although still being part of the SFMS, are already influenced by quenching. This partial suppression of star formation by AGN feedback leads to lower star formation rates and therefore to lower observed slopes. We confirm that including the effects of AGN in our galaxies leads to slopes in agreement with observations. We find the deviation of a galaxy from the SFMS is correlated with its $z=0$ dark matter halo concentration and thus with its halo formation time. This means galaxies with a higher-than-average star formation rate (SFR) form later and vice versa. We explain this apparent correlation with the SFR by re-interpreting galaxies that lie above the SFMS (higher-than-average SFR) as lying to the left of the SFMS (lower-than-average stellar mass) and vice versa. Thus later forming haloes have a lower-than-average stellar mass, this is simply because they have had less-than-average time to form stars, and vice versa. It is thus the nature, i.e. how and when these galaxies form, that sets the path of a galaxy in the SFR versus stellar mass plane.
△ Less
Submitted 6 September, 2020; v1 submitted 31 August, 2020;
originally announced August 2020.