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
A Framework for Modeling Polycyclic Aromatic Hydrocarbon Emission in Galaxy Evolution Simulations
We present a new methodology for simulating mid-infrared emission from polycyclic aromatic hydrocarbons (PAHs) in galaxy evolution simulations. To do this, we combine theoretical models of PAH emission features as they respond to varying interstellar radiation fields, grain-size distributions, and ionization states with a new model for dust evolution in galaxy simulations. We apply these models to three idealized arepo galaxy evolution simulations within the smuggle physics framework. We use these simulations to develop numerical experiments investigating the buildup of PAH masses and luminosities in galaxies in idealized analogs of the Milky Way, a dwarf galaxy, and a starburst disk. Our main results are as follows. Galaxies with high specific star formation rates have increased feedback energy per unit mass, and are able to shatter grains efficiently, driving up the fraction of ultrasmall grains. At the same time, in our model large radiation fields per unit gas density convert aliphatic grains into aromatics. The fraction of dust grains in the form of PAHs (q_PAH) can be understood as a consequence of these processes, and in our model PAHs form primarily from interstellar processing (shattering) of larger grains rather than from the growth of smaller grains. We find that the hardness of the radiation field plays a larger role than variations in the grain-size distribution in setting the total integrated PAH luminosities, though cosmological simulations are necessary to investigate fully the complex interplay of processes that drive PAH band luminosities in galaxies.
more »
« less
- Award ID(s):
- 2008490
- PAR ID:
- 10538549
- Publisher / Repository:
- ApJ
- Date Published:
- Journal Name:
- The Astrophysical Journal
- Volume:
- 951
- Issue:
- 2
- ISSN:
- 0004-637X
- Page Range / eLocation ID:
- 100
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
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
-
ABSTRACT Interstellar chemistry is important for galaxy formation, as it determines the rate at which gas can cool, and enables us to make predictions for observable spectroscopic lines from ions and molecules. We explore two central aspects of modelling the chemistry of the interstellar medium (ISM): (1) the effects of local stellar radiation, which ionizes and heats the gas, and (2) the depletion of metals on to dust grains, which reduces the abundance of metals in the gas phase. We run high-resolution (400 M⊙ per baryonic particle) simulations of isolated disc galaxies, from dwarfs to Milky Way-mass, using the fire galaxy formation models together with the chimes non-equilibrium chemistry and cooling module. In our fiducial model, we couple the chemistry to the stellar fluxes calculated from star particles using an approximate radiative transfer scheme; and we implement an empirical density-dependent prescription for metal depletion. For comparison, we also run simulations with a spatially uniform radiation field, and without metal depletion. Our fiducial model broadly reproduces observed trends in H i and H2 mass with stellar mass, and in line luminosity versus star formation rate for [C ii]$$_{158 \rm {\mu m}}$$, [O i]$$_{63 \rm {\mu m}}$$, [O iii]$$_{88 \rm {\mu m}}$$, [N ii]$$_{122 \rm {\mu m}}$$, and H α6563Å. Our simulations with a uniform radiation field predict fainter luminosities, by up to an order of magnitude for [O iii]$$_{88 \rm {\mu m}}$$ and H α6563Å, while ignoring metal depletion increases the luminosity of carbon and oxygen lines by a factor ≈ 2. However, the overall evolution of the galaxy is not strongly affected by local stellar fluxes or metal depletion, except in dwarf galaxies where the inclusion of local fluxes leads to weaker outflows and hence higher gas fractions.more » « less
-
Abstract We present James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) integral-field spectroscopy of the nearby merging, luminous infrared galaxy, NGC 7469. This galaxy hosts a Seyfert type-1.5 nucleus, a highly ionized outflow, and a bright, circumnuclear star-forming ring, making it an ideal target to study active galactic nucleus (AGN) feedback in the local universe. We take advantage of the high spatial/spectral resolution of JWST/MIRI to isolate the star-forming regions surrounding the central active nucleus and study the properties of the dust and warm molecular gas on ∼100 pc scales. The starburst ring exhibits prominent polycyclic aromatic hydrocarbon (PAH) emission, with grain sizes and ionization states varying by only ∼30%, and a total star formation rate of 10–30 M ⊙ yr −1 derived from fine structure and recombination emission lines. Using pure rotational lines of H 2 we detect 1.2 × 10 7 M ⊙ of warm molecular gas at a temperature higher than 200 K in the ring. All PAH bands get significantly weaker toward the central source, where larger and possibly more ionized grains dominate the emission, likely the result of the ionizing radiation and/or the fast wind emerging from the AGN. The small grains and warm molecular gas in the bright regions of the ring however display properties consistent with normal star-forming regions. These observations highlight the power of JWST to probe the inner regions of dusty, rapidly evolving galaxies for signatures of feedback and inform models that seek to explain the coevolution of supermassive black holes and their hosts.more » « less
-
Abstract We present James Webb Space Telescope (JWST) Near Infrared Spectrograph (NIRSpec) integral field spectroscopy of the nearby luminous infrared galaxy NGC 7469. We take advantage of the high spatial/spectral resolution and wavelength coverage of JWST/NIRSpec to study the 3.3μm neutral polycyclic aromatic hydrocarbon (PAH) grain emission on ∼200 pc scales. A clear change in the average grain properties between the star-forming ring and the central AGN is found. Regions in the vicinity of the AGN, with [Neiii]/[Neii] > 0.25, tend to have larger grain sizes and lower aliphatic-to-aromatic (3.4/3.3) ratios, indicating that smaller grains are preferentially removed by photodestruction in the vicinity of the AGN. PAH emission at the nucleus is weak and shows a low 11.3/3.3 PAH ratio. We find an overall suppression of the total PAH emission relative to the ionized gas in the central 1 kpc region of the AGN in NGC 7469 compared to what has been observed with Spitzer on 3 kpc scales. However, the fractional 3.3μm–to–total PAH power is enhanced in the starburst ring, possibly due to a variety of physical effects on subkiloparsec scales, including recurrent fluorescence of small grains or multiple photon absorption by large grains. Finally, the IFU data show that while the 3.3μm PAH-derived star formation rate (SFR) in the ring is 27% higher than that inferred from the [Neii] and [Neiii] emission lines, the integrated SFR derived from the 3.3μm feature would be underestimated by a factor of 2 due to the deficit of PAHs around the AGN, as might occur if a composite system like NGC 7469 were to be observed at high redshift.more » « less