We present EMPIRE, an IRAM 30 m large program that mapped λ = 3-4 mm dense gas tracers at ∼1-2 kpc resolution across the whole star-forming disk of nine nearby massive spiral galaxies. We describe ...the EMPIRE observing and reduction strategies and show new whole-galaxy maps of HCN(1−0), HCO+(1−0), HNC(1−0), and CO(1−0). We explore how the HCN-to-CO and IR-to-HCN ratios, observational proxies for the dense gas fraction and dense gas star formation efficiency, depend on host galaxy and local environment. We find that the fraction of dense gas correlates with stellar surface density, gas surface density, molecular-to-atomic gas ratio, and dynamical equilibrium pressure. In EMPIRE, the star formation rate per unit dense gas is anticorrelated with these same environmental parameters. Thus, although dense gas appears abundant in the central regions of many spiral galaxies, this gas appears relatively inefficient at forming stars. These results qualitatively agree with previous work on nearby galaxies and the Milky Way's Central Molecular Zone. To first order, EMPIRE demonstrates that the conditions in a galaxy disk set the gas density distribution and that the dense gas traced by HCN shows an environment-dependent relation to star formation. However, our results also show significant ( 0.2 dex) galaxy-to-galaxy variations. We suggest that gas structure below the scale of our observations and dynamical effects likely also play an important role.
We identify stellar structures in the PHANGS sample of 74 nearby galaxies and construct morphological masks of sub-galactic environments based on
Spitzer
3.6
μ
m images. At the simplest level, we ...distinguish five environments: centres, bars, spiral arms, interarm regions, and discs without strong spirals. Slightly more sophisticated masks include rings and lenses, which are publicly released but not explicitly used in this paper. We examine trends with environment in the molecular gas content, star formation rate, and depletion time using PHANGS–ALMA CO(2–1) intensity maps and tracers of star formation. The interarm regions and discs without strong spirals clearly dominate in area, whereas molecular gas and star formation are quite evenly distributed among the five basic environments. We reproduce the molecular Kennicutt–Schmidt relation with a slope compatible with unity within the uncertainties and without significant slope differences among environments. In contrast to what has been suggested by early studies, we find that bars are not always deserts devoid of gas and star formation, but instead they show large diversity. Similarly, spiral arms do not account for most of the gas and star formation in disc galaxies, and they do not have shorter depletion times than the interarm regions. Spiral arms accumulate gas and star formation, without systematically boosting the star formation efficiency. Centres harbour remarkably high surface densities and on average shorter depletion times than other environments. Centres of barred galaxies show higher surface densities and wider distributions compared to the outer disc; yet, depletion times are similar to unbarred galaxies, suggesting highly intermittent periods of star formation when bars episodically drive gas inflow, without enhancing the central star formation efficiency permanently. In conclusion, we provide quantitative evidence that stellar structures in galaxies strongly affect the organisation of molecular gas and star formation, but their impact on star formation efficiency is more subtle.
Star formation is a multi-scale process that requires tracing cloud formation and stellar feedback within the local ( kpc) and global galaxy environment. We present first results from two large ...observing programs on the Atacama Large Millimeter/submillimeter Array (ALMA)and the Very Large Telescope/Multi Unit Spectroscopic Explorer(VLT/MUSE), mapping cloud scales (1″ = 47 pc) in both molecular gas and star-forming tracers across 90 kpc2 of the central disk of NGC 628 to probe the physics of star formation. Systematic spatial offsets between molecular clouds and H ii regions illustrate the time evolution of star-forming regions. Using uniform sampling of both maps on 50-500 pc scales, we infer molecular gas depletion times of 1-3 Gyr, but also find that the increase of scatter in the star formation relation on small scales is consistent with gas and H ii regions being only weakly correlated at the cloud (50 pc) scale. This implies a short overlap phase for molecular clouds and H ii regions, which we test by directly matching our catalog of 1502 H ii regions and 738 GMCs. We uncover only 74 objects in the overlap phase, and we find depletion times >1 Gyr, significantly longer than previously reported for individual star-forming clouds in the Milky Way. Finally, we find no clear trends that relate variations in the depletion time observed on 500 pc scales to physical drivers (metallicity, molecular and stellar-mass surface density, molecular gas boundedness) on 50 pc scales.
Interstellar dust and starlight are modeled for the galaxies of the project "Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel." The galaxies were observed by the Infrared Array ...Camera and the Multiband Imaging Photometer for Spitzer on Spitzer Space Telescope, and the Photodetector Array Camera and Spectrometer and the Spectral and Photometric Imaging Receiver on Herschel Space Observatory. With data from 3.6 to 500 m, dust models are strongly constrained. Using a physical dust model, for each pixel in each galaxy we estimate (1) dust surface density, (2) dust mass fraction in polycyclic aromatic hydrocarbons (PAHs), (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in subregions with high starlight intensity. The dust models successfully reproduce the observed global and resolved spectral energy distributions. With the angular resolution of Herschel, we obtain well-resolved maps (available online) for the dust properties. As in previous studies, we find the PAH fraction to be an increasing function of metallicity, with a threshold oxygen abundance Z/Z 0.1, but we find the data to be fitted best with increasing linearly with above a threshold value of 0.15(O/H) . We obtain total dust masses for each galaxy by summing the dust mass over the individual map pixels; these "resolved" dust masses are consistent with the masses inferred from a model fit to the global photometry. The global dust-to-gas ratios obtained from this study are found to correlate with galaxy metallicities. Systems with Z/Z 0.5 have most of their refractory elements locked up in dust, whereas in systems with Z/Z 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that is suppressed in regions with unusually warm dust with . With knowledge of one long-wavelength flux density ratio (e.g., f160/f500), the minimum starlight intensity heating the dust ( ) can be estimated to within ∼50%, despite a variation in of more than two orders of magnitude. For the adopted dust model, dust masses can be estimated to within ∼0.2 dex accuracy using the f160/f500 flux ratio and the integrated dust luminosity, and to ∼0.07 dex accuracy using the 500 m luminosity alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions for estimating starlight heating intensity and dust mass may be useful for studies of high-redshift galaxies.
Aims.
The complexity of star formation at the physical scale of molecular clouds is not yet fully understood. We investigate the mechanisms regulating the formation of stars in different environments ...within nearby star-forming galaxies from the Physics at High Angular resolution in Nearby GalaxieS (PHANGS) sample.
Methods.
Integral field spectroscopic data and radio-interferometric observations of 18 galaxies were combined to explore the existence of the resolved star formation main sequence (Σ
stellar
versus Σ
SFR
), resolved Kennicutt–Schmidt relation (Σ
mol. gas
versus Σ
SFR
), and resolved molecular gas main sequence (Σ
stellar
versus Σ
mol. gas
), and we derived their slope and scatter at spatial resolutions from 100 pc to 1 kpc (under various assumptions).
Results.
All three relations were recovered at the highest spatial resolution (100 pc). Furthermore, significant variations in these scaling relations were observed across different galactic environments. The exclusion of non-detections has a systematic impact on the inferred slope as a function of the spatial scale. Finally, the scatter of the Σ
mol. gas + stellar
versus Σ
SFR
correlation is smaller than that of the resolved star formation main sequence, but higher than that found for the resolved Kennicutt–Schmidt relation.
Conclusions.
The resolved molecular gas main sequence has the tightest relation at a spatial scale of 100 pc (scatter of 0.34 dex), followed by the resolved Kennicutt–Schmidt relation (0.41 dex) and then the resolved star formation main sequence (0.51 dex). This is consistent with expectations from the timescales involved in the evolutionary cycle of molecular clouds. Surprisingly, the resolved Kennicutt–Schmidt relation shows the least variation across galaxies and environments, suggesting a tight link between molecular gas and subsequent star formation. The scatter of the three relations decreases at lower spatial resolutions, with the resolved Kennicutt–Schmidt relation being the tightest (0.27 dex) at a spatial scale of 1 kpc. Variation in the slope of the resolved star formation main sequence among galaxies is partially due to different detection fractions of Σ
SFR
with respect to Σ
stellar
.
ABSTRACT We analyze the velocity dispersions of individual H i and CO profiles in a number of nearby galaxies from the high-resolution HERACLES CO and THINGS H i surveys. Focusing on regions with ...bright CO emission, we find a CO dispersion value = 7.3 1.7 km s−1. The corresponding H i dispersion = 11.7 2.3 km s−1, yielding a mean dispersion ratio = 1.4 0.2, independent of radius. We find that the CO velocity dispersion increases toward lower peak fluxes. This is consistent with previous work where we showed that when using spectra averaged ("stacked") over large areas, larger values for the CO dispersion are found, and a lower ratio = 1.0 0.2. The stacking method is more sensitive to low-level diffuse emission, whereas individual profiles trace narrow-line, GMC-dominated, bright emission. These results provide further evidence that disk galaxies contain not only a thin, low velocity dispersion, high density CO disk that is dominated by GMCs, but also a fainter, higher dispersion, diffuse disk component.
The distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined ...the gas phase oxygen abundances for 7138 H ii regions across the disks of eight nearby galaxies using Very Large Telescope/Multi Unit Spectroscopic Explorer (MUSE) optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first-order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (ΔO/H) and explore azimuthal variations. Across each galaxy, we find low ( = 0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those H ii regions that are 1 outliers toward both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher H luminosity, lower H velocity dispersion, younger star clusters, and associated molecular gas clouds showing higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased H velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but cannot always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the interstellar medium.
Abstract
We present new 3 mm observations of the ionized gas toward the nuclear starburst in the nearby (
D
∼ 3.5 Mpc) galaxy NGC 253. With ALMA, we detect emission from the H40
α
and He40
α
lines in ...the central 200 pc of this galaxy on spatial scales of ∼4 pc. The recombination line emission primarily originates from a population of approximately a dozen embedded super star clusters in the early stages of formation. We find that emission from these clusters is characterized by electron temperatures ranging from 7000 to 10,000 K and measures an average singly ionized helium abundance 〈
Y
+
〉 = 0.25 ± 0.06, both of which are consistent with values measured for H
ii
regions in the center of the Milky Way. We also report the discovery of unusually broad line width recombination line emission originating from seven of the embedded clusters. We suggest that these clusters contribute to the launching of the large-scale hot wind observed to emanate from the central starburst. Finally, we use the measured recombination line fluxes to improve the characterization of overall embedded cluster properties, including the distribution of cluster masses and the fractional contribution of the clustered star formation to the total starburst, which we estimate is at least 50%.
We present a study correlating the spatial locations of young star clusters with those of molecular clouds in NGC 5194, in order to investigate the time-scale over which clusters separate from their ...birth clouds. The star cluster catalogues are from the Legacy ExtraGalactic UV Survey (LEGUS) and the molecular clouds from the Plateau de Bure Interefrometer Arcsecond Whirpool Survey (PAWS). We find that younger star clusters are spatially closer to molecular clouds than older star clusters. The median age for clusters associated with clouds is 4 Myr, whereas it is 50 Myr for clusters that are sufficiently separated from a molecular cloud to be considered unassociated. After ∼6 Myr, the majority of the star clusters lose association with their molecular gas. Younger star clusters are also preferentially located in stellar spiral arms where they are hierarchically distributed in kpc-size regions for 50-100 Myr before dispersing. The youngest star clusters are more strongly clustered, yielding a two-point correlation function with α = -0.28 ± 0.04, than the giant molecular cloud (GMCs) (α = -0.09 ± 0.03) within the same PAWS field. However, the clustering strength of the most massive GMCs, supposedly the progenitors of the young clusters for a star formation efficiency of a few per cent, is comparable (α = -0.35 ± 0.05) to that of the clusters. We find a galactocentric dependence for the coherence of star formation, in which clusters located in the inner region of the galaxy reside in smaller star-forming complexes and display more homogeneous distributions than clusters further from the centre. This result suggests a correlation between the survival of a cluster complex and its environment.
Abstract
High critical density molecular lines like HCN (1-0) or HCO+ (1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth ...of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high-density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC (1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C (1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large programme EMPIRE and new Atacama Large Millimetre/submillimetre Array observations, which together target six nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN (1-0) and HCO+ (1-0) to be τ ∼ 2–11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5 and 20 × 105 cm−3, 1 and 3 × 105 cm−3 and 9 × 104 cm−3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, η, on these estimates and find that the effective critical densities decrease by a factor of
$\frac{\eta _{12}}{\eta _{13}}\,\tau _{12}$
. A comparison to existing work in NGC 5194 and NGC 253 shows the HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of active galactic nucleus dominated nearby galaxies.
PDF
