Search for: All records

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.

 

We measure the CO-to-H₂conversion factor (α_CO) in 37 galaxies at 2 kpc resolution, using the dust surface density inferred from far-infrared emission as a tracer of the gas surface density and assuming a constant dust-to-metal ratio. In total, we have ∼790 and ∼610 independent measurements ofα_COfor CO (2–1) and (1–0), respectively. The mean values forα_{CO (2–1)}andα_{CO (1–0)}are${9.3}_{-5.4}^{+4.6}$and${4.2}_{-2.0}^{+1.9}{M}_{\odot }{\mathrm{pc}}^{-2}{(\mathrm{K}\mathrm{km}{\mathrm{s}}^{-1})}^{-1}$, respectively. The CO-intensity-weighted mean is 5.69 forα_{CO (2–1)}and 3.33 forα_{CO (1–0)}. We examine howα_COscales with several physical quantities, e.g., the star formation rate (SFR), stellar mass, and dust-mass-weighted average interstellar radiation field strength ($\overline{U}$). Among them,$\overline{U}$, Σ_SFR, and the integrated CO intensity (W_CO) have the strongest anticorrelation with spatially resolvedα_CO. We provide linear regression results toα_COfor all quantities tested. At galaxy-integrated scales, we observe significant correlations betweenα_COandW_CO, metallicity,$\overline{U}$, and Σ_SFR. We also find thatα_COin each galaxy decreases with the stellar mass surface density (Σ_⋆) in high-surface-density regions (Σ_⋆≥ 100M_⊙pc⁻²), following the power-law relations${\alpha }_{\mathrm{CO}(2–1)}\propto {\mathrm{\Sigma }}_{\star }^{-0.5}$and${\alpha }_{\mathrm{CO}(1–0)}\propto {\mathrm{\Sigma }}_{\star }^{-0.2}$. The power-law index is insensitive to the assumed dust-to-metal ratio. We interpret the decrease inα_COwith increasing Σ_⋆as a result of higher velocity dispersion compared to isolated, self-gravitating clouds due to the additional gravitational force from stellar sources, which leads to the reduction inα_CO. The decrease inα_COat high Σ_⋆is important for accurately assessing molecular gas content and star formation efficiency in the centers of galaxies, which bridge “Milky Way–like” to “starburst-like” conversion factors.

 

Young stellar objects (YSOs) are the gold standard for tracing star formation in galaxies but have been unobservable beyond the Milky Way and Magellanic Clouds. But that all changed when the JWST was launched, which we use to identify YSOs in the Local Group galaxy M33, marking the first time that individual YSOs have been identified at these large distances. We present Mid-Infrared Instrument (MIRI) imaging mosaics at 5.6 and 21 $\mu$m that cover a significant portion of one of M33’s spiral arms that has existing panchromatic imaging from the Hubble Space Telescope and deep Atacama Large Millimeter/submillimeter Array CO measurements. Using these MIRI and Hubble Space Telescope images, we identify point sources using the new dolphot MIRI module. We identify 793 candidate YSOs from cuts based on colour, proximity to giant molecular clouds (GMCs), and visual inspection. Similar to Milky Way GMCs, we find that higher mass GMCs contain more YSOs and YSO emission, which further show YSOs identify star formation better than most tracers that cannot capture this relationship at cloud scales. We find evidence of enhanced star formation efficiency in the southern spiral arm by comparing the YSOs to the molecular gas mass.

 

Determining how the galactic environment, especially the high gas densities and complex dynamics in bar-fed galaxy centers, alters the star formation efficiency (SFE) of molecular gas is critical to understanding galaxy evolution. However, these same physical or dynamical effects also alter the emissivity properties of CO, leading to variations in the CO-to-H₂conversion factor (α_CO) that impact the assessment of the gas column densities and thus of the SFE. To address such issues, we investigate the dependence ofα_COon the local CO velocity dispersion at 150 pc scales using a new set of dust-basedα_COmeasurements and propose a newα_COprescription that accounts for CO emissivity variations across galaxies. Based on this prescription, we estimate the SFE in a sample of 65 galaxies from the PHANGS–Atacama Large Millimeter/submillimeter Array survey. We find increasing SFE toward high-surface-density regions like galaxy centers, while using a constant or metallicity-basedα_COresults in a more homogeneous SFE throughout the centers and disks. Our prescription further reveals a mean molecular gas depletion time of 700 Myr in the centers of barred galaxies, which is overall three to four times shorter than in nonbarred galaxy centers or the disks. Across the galaxy disks, the depletion time is consistently around 2–3 Gyr, regardless of the choice ofα_COprescription. All together, our results suggest that the high level of star formation activity in barred centers is not simply due to an increased amount of molecular gas, but also to an enhanced SFE compared to nonbarred centers or disk regions.

 

Abstract The CO-to-H 2 conversion factor ( α CO ) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower α CO in the centers of some barred galaxies on kiloparsec scales. To unveil the physical drivers of such variations, we obtained Atacama Large Millimeter/submillimeter Array bands (3), (6), and (7) observations toward the inner ∼2 kpc of NGC 3627 and NGC 4321 tracing 12 CO, 13 CO, and C 18 O lines on ∼100 pc scales. Our multiline modeling and Bayesian likelihood analysis of these data sets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of α CO . The central 300 pc nuclei in both galaxies show strong enhancement of temperature T k ≳ 100 K and density n H 2 > 10 3 cm −3 . Assuming a CO-to-H 2 abundance of 3 × 10 −4 , we derive 4–15 times lower α CO than the Galactic value across our maps, which agrees well with previous kiloparsec-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of α CO with low- J 12 CO optical depths ( τ CO ), as well as an anticorrelation with T k . The τ CO correlation explains most of the α CO variation in the three galaxy centers, whereas changes in T k influence α CO to second order. Overall, the observed line width and 12 CO/ 13 CO 2–1 line ratio correlate with τ CO variation in these centers, and thus they are useful observational indicators for α CO variation. We also test current simulation-based α CO prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. 

Recent strides have been made developing dust evolution models for galaxy formation simulations but these approaches vary in their assumptions and degree of complexity. Here, we introduce and compare two separate dust evolution models (labelled ‘Elemental’ and ‘Species’), based on recent approaches, incorporated into the gizmo code and coupled with fire-2 stellar feedback and interstellar medium physics. Both models account for turbulent dust diffusion, stellar production of dust, dust growth via gas-dust accretion, and dust destruction from time-resolved supernovae, thermal sputtering in hot gas, and astration. The ‘Elemental’ model tracks the evolution of generalized dust species and utilizes a simple, ‘tunable’ dust growth routine, while the ‘Species’ model tracks the evolution of specific dust species with set chemical compositions and incorporates a physically motivated, two-phase dust growth routine. We test and compare these models in an idealized Milky Way-mass galaxy and find that while both produce reasonable galaxy-integrated dust-to-metals (D/Z) ratios and predict gas-dust accretion as the main dust growth mechanism, a chemically motivated model is needed to reproduce the observed scaling relation between individual element depletions and D/Z with column density and local gas density. We also find the inclusion of theoretical metallic iron and O-bearing dust species are needed in the case of specific dust species in order to match observations of O and Fe depletions, and the integration of a sub-resolution dense molecular gas/CO scheme is needed to both match observed C depletions and ensure carbonaceous dust is not overproduced in dense environments.

 

Abstract The CO-to-H 2 conversion factor ( α CO ) is critical to studying molecular gas and star formation in galaxies. The value of α CO has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on ~kiloparsec scales revealed low values of α CO in the centers of some barred spiral galaxies, including NGC 3351. We present new Atacama Large Millimeter/submillimeter Array Band 3, 6, and 7 observations of 12 CO, 13 CO, and C 18 O lines on 100 pc scales in the inner ∼2 kpc of NGC 3351. Using multiline radiative transfer modeling and a Bayesian likelihood analysis, we infer the H 2 density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of 2–3 × 10 3 cm −3 in the central ∼1 kpc and a high temperature of 30–60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H 2 abundance of 3 × 10 −4 , our analysis yields α CO ∼ 0.5–2.0 M ⊙ (K km s −1 pc 2 ) −1 with a decreasing trend with galactocentric radius in the central ∼1 kpc. The inflows show a substantially lower α CO ≲ 0.1 M ⊙ (K km s −1 pc 2 ) −1 , likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted α CO of ∼1.5 M ⊙ (K km s −1 pc 2 ) −1 , which is similar to previous dust-modeling-based results at kiloparsec scales. This suggests that low α CO on kiloparsec scales in the centers of some barred galaxies may be due to the contribution of low-optical-depth CO emission in bar-driven inflows. 

We use PHANGS–James Webb Space Telescope (JWST) data to identify and classify 1271 compact 21μm sources in four nearby galaxies using MIRI F2100W data. We identify sources using a dendrogram-based algorithm, and we measure the background-subtracted flux densities for JWST bands from 2 to 21μm. Using the spectral energy distribution (SED) in JWST and HST bands plus ALMA and MUSE/VLT observations, we classify the sources by eye. Then we use this classification to define regions in color–color space and so establish a quantitative framework for classifying sources. We identify 1085 sources as belonging to the ISM of the target galaxies with the remainder being dusty stars or background galaxies. These 21μm sources are strongly spatially associated with Hiiregions (>92% of sources), while 74% of the sources are coincident with a stellar association defined in the HST data. Using SED fitting, we find that the stellar masses of the 21μm sources span a range of 10²–10⁴M_⊙with mass-weighted ages down to 2 Myr. There is a tight correlation between attenuation-corrected Hαand 21μm luminosity forL_ν,F2100W> 10¹⁹W Hz⁻¹. Young embedded source candidates selected at 21μm are found below this threshold and haveM_⋆< 10³M_⊙.

 

ABSTRACT Previous work has argued that atomic gas mass estimates of galaxies from 21-cm H i emission are systematically low due to a cold opaque atomic gas component. If true, this opaque component necessitates a $\sim 35{{\ \rm per\ cent}}$ correction factor relative to the mass from assuming optically thin H i emission. These mass corrections are based on fitting H i spectra with a single opaque component model that produces a distinct ‘top-hat’ shaped line profile. Here, we investigate this issue using deep, high spectral resolution H i VLA observations of M31 and M33 to test if these top-hat profiles are instead superpositions of multiple H i components along the line of sight. We fit both models and find that ${\gt}80{{\ \rm per\ cent}}$ of the spectra strongly prefer a multicomponent Gaussian model while ${\lt}2{{\ \rm per\ cent}}$ prefer the single opacity-corrected component model. This strong preference for multiple components argues against previous findings of lines of sight dominated by only cold H i. Our findings are enabled by the improved spectral resolution (0.42 ${\rm km\, s^{-1}}$), whereas coarser spectral resolution blends multiple components together. We also show that the inferred opaque atomic ISM mass strongly depends on the goodness-of-fit definition and is highly uncertain when the inferred spin temperature has a large uncertainty. Finally, we find that the relation of the H i surface density with the dust surface density and extinction has significantly more scatter when the inferred H i opacity correction is applied. These variations are difficult to explain without additionally requiring large variations in the dust properties. Based on these findings, we suggest that the opaque H i mass is best constrained by H i absorption studies.