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.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high far-infrared surface brightness, and combine these with local universe ...CO outflow results from the literature. The CO line ratios and spatial outflow structure of our targets provide some constraints on the conversion steps from observables to physical quantities such as molecular mass outflow rates. Where available, ratios between outflow emission in higher J CO transitions and in CO(1-0) are typically consistent with excitation
R
i
1
≲ 1. However, for IRAS 13120−5453,
R
31
= 2.10 ± 0.29 indicates optically thin CO in the outflow. Like much of the outflow literature, we use
α
CO(1 − 0)
= 0.8, and we present arguments for using
C
= 1 in deriving molecular mass outflow rates
Ṁ
out
=
C
M
out
v
out
/
R
out
. We compare the two main methods for molecular outflow detection: CO millimeter interferometry and
Herschel
OH-based spectroscopic outflow searches. For 26 sources studied with both methods, we find an 80% agreement in detecting
v
out
≳ 150 km s
−1
outflows, and non-matches can be plausibly ascribed to outflow geometry and signal-to-noise ratio. For a published sample of 12 bright ultraluminous infrared galaxies with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with active galactic nucleus (AGN) luminosity and with bolometric luminosity than with far-infrared surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in nonAGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short flow times of ∼10
6
yr and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGNs.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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
.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
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.
ABSTRACT Distance uncertainties plague our understanding of the physical scales relevant to the physics of star formation in extragalactic studies. The planetary nebulae luminosity function (PNLF) is ...one of very few techniques that can provide distance estimates to within ∼10%; however, it requires a planetary nebula (PN) sample that is uncontaminated by other ionizing sources. We employ optical integral field unit spectroscopy using the Multi-Unit Spectroscopic Explorer on the Very Large Telescope to measure O iii line fluxes for sources unresolved on 50 pc scales within the central star-forming galaxy disk of NGC 628. We use diagnostic line ratios to identify 62 PNe, 30 supernova remnants, and 87 H ii regions within our fields. Using the 36 brightest PNe, we determine a new PNLF distance modulus of 29.91 − 0.13 + 0.08 mag (9.59 − 0.57 + 0.35 Mpc), which is in good agreement with literature values, but significantly larger than the previously reported PNLF distance. We are able to explain the discrepancy and recover the previous result when we reintroduce SNR contaminants to our sample. This demonstrates the power of full spectral information over narrowband imaging in isolating PNe. Given our limited spatial coverage within the Galaxy, we show that this technique can be used to refine distance estimates, even when IFU observations cover only a fraction of a galaxy disk.
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 a follow-up analysis examining the dynamics and structures of 41 massive, large star-forming galaxies at
z
∼ 0.67 − 2.45 using both ionized and molecular gas kinematics. We fit ...the galaxy dynamics with models consisting of a bulge, a thick, turbulent disk, and an NFW dark matter halo, using code that fully forward-models the kinematics, including all observational and instrumental effects. We explore the parameter space using Markov Chain Monte Carlo (MCMC) sampling, including priors based on stellar and gas masses and disk sizes. We fit the full sample using extracted 1D kinematic profiles. For a subset of 14 well-resolved galaxies, we also fit the 2D kinematics. The MCMC approach robustly confirms the results from least-squares fitting presented in Paper I: the sample galaxies tend to be baryon-rich on galactic scales (within one effective radius). The 1D and 2D MCMC results are also in good agreement for the subset, demonstrating that much of the galaxy dynamical information is captured along the major axis. The 2D kinematics are more affected by the presence of noncircular motions, which we illustrate by constructing a toy model with constant inflow for one galaxy that exhibits residual signatures consistent with radial motions. This analysis, together with results from Paper I and other studies, strengthens the finding that massive, star-forming galaxies at
z
∼ 1 − 2 are baryon-dominated on galactic scales, with lower dark matter fractions toward higher baryonic surface densities. Finally, we present details of the kinematic fitting code used in this analysis.
We use integral field spectroscopy from the PHANGS–MUSE survey, which resolves the ionised interstellar medium structure at ∼50 pc resolution in 19 nearby spiral galaxies, to study the origin of the ...diffuse ionised gas (DIG). We examine the physical conditions of the diffuse gas by first removing morphologically defined H
II
regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines in the DIG. A simple model for the leakage and propagation of ionising radiation from H
II
regions is able to reproduce the observed distribution of H
α
in the DIG. This model infers a typical mean free path for the ionising radiation of 1.9 kpc for photons propagating within the disc plane. Leaking radiation from H
II
regions also explains the observed decrease in line ratios of low-ionisation species (S
II
/H
α
, N
II
/H
α
, and O
I
/H
α
) with increasing H
α
surface brightness (Σ
H
α
). Emission from hot low-mass evolved stars, however, is required to explain: (1) the enhanced low-ionisation line ratios observed in the central regions of some of the galaxies in our sample; (2) the observed trends of a flat or decreasing O
III
/H
β
with Σ
H
α
; and (3) the offset of some DIG regions from the typical locus of H
II
regions in the Baldwin–Phillips–Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LINERs). Hot low-mass evolved stars make a small contribution to the energy budget of the DIG (2% of the galaxy-integrated H
α
emission), but their harder spectra make them fundamental contributors to O
III
emission. The DIG might result from a superposition of two components, an energetically dominant contribution from young stars and a more diffuse background of harder ionising photons from old stars. This unified framework bridges observations of the Milky Way DIG with LI(N)ER-like emission observed in nearby galaxy bulges.
Full text
Available for:
FMFMET, NUK, UL, UM, UPUK
We combine new sensitive, wide-field CO data from the HERACLES survey with ultraviolet and infrared data from GALEX and Spitzer to compare the surface densities of H2, Delta *SH2, and the recent star ...formation rate, Delta *SSFR, over many thousands of positions in 30 nearby disk galaxies. We more than quadruple the size of the galaxy sample compared to previous work and include targets with a wide range of galaxy properties. Even though the disk galaxies in this study span a wide range of properties, we find a strong, and approximately linear correlation between Delta *SSFR and Delta *SH2 at our common resolution of 1 kpc. This implies a roughly constant median H2 consumption time, Delta *tH2 Dep = Delta *SH2/ Delta *SSFR, of ~2.35 Gyr (including heavy elements) across our sample. At 1 kpc resolution, there is only a weak correlation between Delta *SH2 and Delta *tH2 Dep over the range Delta *SH2 5-100 M pc--2, which is probed by our data. We compile a broad set of literature measurements that have been obtained using a variety of star formation tracers, sampling schemes, and physical scales and show that overall, these data yield almost exactly the same results, although with more scatter. We interpret these results as strong, albeit indirect evidence that star formation proceeds in a uniform way in giant molecular clouds in the disks of spiral galaxies.