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 non-linear 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 non-linear on time-scales 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.
In dusty cool-star outflow or ejection events around asymptotic giant branch (AGB) or R Coronae Borealis or RCB-like stars, dust is accelerated by radiation from the star and coupled to the gas via collisional drag forces. It has recently been shown that such dust-gas mixtures are unstable to a super-class of instabilities called the resonant drag instabilities (RDIs), which promote dust clustering. We therefore consider idealized simulations of the RDIs operating on a spectrum of dust grain sizes subject to radiative acceleration (allowing for different grain optical properties), coupled to the gas with a realistic drag law, including or excluding the effects of magnetic fields and charged grains, and calculate for the first time how the RDIs could contribute to observed variability. We show that the RDIs naturally produce significant variations (spatially and temporally) ($\sim 10\!-\!20{{\ \rm per\ cent}}$ 1 σ-level) in the extinction, corresponding to $\sim 0.1\!-\!1\,$mag level in the stellar types above, on time-scales of order months to a year. The fluctuations are surprisingly robust to the assumed size of the source as they are dominated by large-scale modes, which also means their spatial structure could be resolved in some nearby systems. We also quantify how this produces variations in the line-of-sight grain size-distribution. All of these variations are similar to those observed, suggesting that the RDIs may play a key role driving observed spatial and temporal variability in dust extinction within dusty outflow/ejection events around cool stars. We further propose that the measured variations in grain sizes could directly be used to identify the presence of the RDIs in close by systems with observations.
more » « less- NSF-PAR ID:
- 10369903
- Publisher / Repository:
- Oxford University Press
- Date Published:
- Journal Name:
- Monthly Notices of the Royal Astronomical Society
- Volume:
- 515
- Issue:
- 4
- ISSN:
- 0035-8711
- Page Range / eLocation ID:
- p. 4797-4809
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
ABSTRACT -
ABSTRACT We study the non-linear evolution of the acoustic ‘resonant drag instability’ (RDI) using numerical simulations. The acoustic RDI is excited in a dust–gas mixture when dust grains stream through gas, interacting with sound waves to cause a linear instability. We study this process in a periodic box by accelerating neutral dust with an external driving force. The instability grows as predicted by linear theory, eventually breaking into turbulence and saturating. As in linear theory, the non-linear behaviour is characterized by three regimes – high, intermediate, and low wavenumbers – the boundary between which is determined by the dust–gas coupling strength and the dust-to-gas mass ratio. The high and intermediate wavenumber regimes behave similarly to one another, with large dust-to-gas ratio fluctuations while the gas remains largely incompressible. The saturated state is highly anisotropic: dust is concentrated in filaments, jets, or plumes along the direction of acceleration, with turbulent vortex-like structures rapidly forming and dissipating in the perpendicular directions. The low-wavenumber regime exhibits large fluctuations in gas and dust density, but the dust and gas remain more strongly coupled in coherent ‘fronts’ perpendicular to the acceleration. These behaviours are qualitatively different from those of dust ‘passively’ driven by external hydrodynamic turbulence, with no back-reaction force from dust on to gas. The virulent nature of these instabilities has interesting implications for dust-driven winds in a variety of astrophysical systems, including around cool stars, in dusty torii around active-galactic-nuclei, and in and around giant molecular clouds.more » « less
-
more » « less
ABSTRACT We study the linear growth and non-linear saturation of the ‘acoustic Resonant Drag Instability’ (RDI) when the dust grains, which drive the instability, have a wide, continuous spectrum of different sizes. This physics is generally applicable to dusty winds driven by radiation pressure, such as occurs around red-giant stars, star-forming regions, or active galactic nuclei. Depending on the physical size of the grains compared to the wavelength of the radiation field that drives the wind, two qualitatively different regimes emerge. In the case of grains that are larger than the radiation’s wavelength – termed the constant-drift regime – the grain’s equilibrium drift velocity through the gas is approximately independent of grain size, leading to strong correlations between differently sized grains that persist well into the saturated non-linear turbulence. For grains that are smaller than the radiation’s wavelength – termed the non-constant-drift regime – the linear instability grows more slowly than the single-grain-size RDI and only the larger grains exhibit RDI-like behaviour in the saturated state. A detailed study of grain clumping and grain–grain collisions shows that outflows in the constant-drift regime may be effective sites for grain growth through collisions, with large collision rates but low collision velocities.
-
more » « less
ABSTRACT Radiation pressure on dust is an important feedback process around star clusters and may eject gas from bright subregions in star-forming galaxies. The Eddington ratio has previously been constructed for galaxy-averaged observations, individual star clusters, and Galactic H ii regions. Here we assess the role of radiation pressure in thousands of subregions across two local star-forming galaxies, NGC 6946 and NGC 5194. Using a model for the spectral energy distribution from stellar population synthesis and realistic dust grain scattering and absorption, we compute flux and radiation pressure-mean opacities and population-averaged optical depth 〈τRP〉. Using Monte-Carlo calculations, we assess the momentum coupling through a dusty column to the stellar continuum. Optically-thin regions around young stellar populations are 30–50 times super-Eddington. We calculate the Eddington ratio for the subregions including the local mass of young and old stars and cool atomic and molecular gas. We compute the fraction of the total star formation that is currently super-Eddington, and provide an assessment of the role of radiation pressure in the dusty gas dynamics. Depending on the assumed height of the dusty gas and the age of the stellar population, we find that ∼0–10 per cent of the sightlines are super-Eddington. These regions may be accelerated to ∼5–15 km s−1 by radiation pressure alone. Additionally, our results show that for beamed radiation, the function 1 − exp (−〈τRP〉) is an excellent approximation to the momentum transfer. Opacities and optical depths are tabulated for SEDs of different stellar ages and for continuous star formation.
-
null (Ed.)Context. The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: star birth and death, grain formation and destruction, and galactic inflows and outflows. Understanding such processes and their relative importance is essential to following galaxy evolution and the chemical enrichment through the cosmic epochs, and to interpreting current and future observations. Despite the importance of such topics, the contribution of different stellar sources to the chemical enrichment of galaxies, for example massive stars exploding as Type II supernovae (SNe) and low-mass stars, as well as the mechanisms driving the evolution of dust grains, such as for example grain growth in the ISM and destruction by SN shocks, remain controversial from both observational and theoretical viewpoints. Aims. In this work, we revise the current description of metal and dust evolution in the ISM of local low-metallicity dwarf galaxies and develop a new description of Lyman-break galaxies (LBGs) which are considered to be their high-redshift counterparts in terms of star formation, stellar mass, and metallicity. Our goal is to reproduce the observed properties of such galaxies, in particular (i) the peak in dust mass over total stellar mass (sMdust) observed within a few hundred million years; and (ii) the decrease in sMdust at a later time. Methods. We fitted spectral energy distribution of dwarf galaxies and LBGs with the “Code Investigating GALaxies Emission”, through which the total stellar mass, dust mass, and star formation rate are estimated. For some of the dwarf galaxies considered, the metal and gas content are available from the literature. We computed different prescriptions for metal and dust evolution in these systems (e.g. different initial mass functions for stars, dust condensation fractions, SN destruction, dust accretion in the ISM, and inflow and outflow efficiency), and we fitted the properties of the observed galaxies through the predictions of the models. Results. Only some combinations of models are able to reproduce the observed trend and simultaneously fit the observed properties of the galaxies considered. In particular, we show that (i) a top-heavy initial mass function that favours the formation of massive stars and a dust condensation fraction for Type II SNe of around 50% or more help to reproduce the peak of sMdust observed after ≈100 Myr from the beginning of the baryon cycle for both dwarf galaxies and LBGs; (ii) galactic outflows play a crucial role in reproducing the observed decline in sMdust with age and are more efficient than grain destruction from Type II SNe both in local galaxies and at high-redshift; (iii) a star formation efficiency (mass of gas converted into stars) of a few percent is required to explain the observed metallicity of local dwarf galaxies; and (iv) dust growth in the ISM is not necessary in order to reproduce the values of sMdust derived for the galaxies under study, and, if present, the effect of this process would be erased by galactic outflows.more » « less
