Pillars of creation amongst destruction: Star formation in molecular clouds near R136 in 30 Doradus

New sensitive CO(2-1) observations of the 30 Doradus region in the Large Magellanic Cloud are presented. We identify a chain of three newly discovered molecular clouds we name KN1, KN2 and KN3 lying within 2--14 pc in projection from the young massive cluster R136 in 30 Doradus. Excited H$_2$ 2.12$\mu$m emission is spatially coincident with the molecular clouds, but ionized Br$\gamma$ emission is not. We interpret these observations as the tails of pillar-like structures whose ionized heads are pointing towards R136. Based on infrared photometry, we identify a new generation of stars forming within this structure.Comment: Accepted for publication in ApJ (includes 13 pages, 8 figures). For higher resolution figures please see http://www.das.uchile.cl/~vkalari/staplervk.pd

A 50 pc Scale View of Star Formation Efficiency across NGC 628

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.We thank the referee for helpful comments that improved this work. K.K. gratefully acknowledges support from grant KR 4598/1-2 from the German Research Foundation (DFG) Priority Program 1573. J.M.D.K. and M.C. gratefully acknowledge funding from the DFG in the form of an Emmy Noether Research Group (grant No. KR4801/1-1). J.M.D.K. gratefully acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 Research and Innovation Programme via the ERC Starting Grant MUSTANG (grant agreement No. 714907). B.G. gratefully acknowledges the support of the Australian Research Council as the recipient of a Future Fellowship (FT140101202). F.B. acknowledges funding from the European Union’s Horizon 2020 Research and Innovation Programme (grant agreement No. 726384—EMPIRE). G.B. is supported by CONICYT/ FONDECYT, Programa de Iniciación, Folio 11150220. A.H. acknowledges support from the Centre National d’Etudes Spatiales (CNES). E.R. acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), funding reference No. RGPIN-2017-03987. R.M. and E.S. acknowledge funding from the ERC under the European Union’s Horizon 2020 Research and Innovation Programme (grant agreement No. 694343). J.P. acknowledges support by the Programme National “Physique et Chimie du Milieu Interstellaire”(PCMI) of CNRS/INSU with INC/INP co-funded by CEA and CNES

PHANGS CO kinematics: disk orientations and rotation curves at 150 pc resolution

We present kinematic orientations and high resolution (150 pc) rotation curves for 67 main sequence star-forming galaxies surveyed in CO (2-1) emission by PHANGS-ALMA. Our measurements are based on the application of a new fitting method tailored to CO velocity fields. Our approach identifies an optimal global orientation as a way to reduce the impact of non-axisymmetric (bar and spiral) features and the uneven spatial sampling characteristic of CO emission in the inner regions of nearby galaxies. The method performs especially well when applied to the large number of independent lines-of-sight contained in the PHANGS CO velocity fields mapped at 1'' resolution. The high resolution rotation curves fitted to these data are sensitive probes of mass distribution in the inner regions of these galaxies. We use the inner slope as well as the amplitude of our fitted rotation curves to demonstrate that CO is a reliable global dynamical mass tracer. From the consistency between photometric orientations from the literature and kinematic orientations determined with our method, we infer that the shapes of stellar disks in the mass range of log($\rm M_{\star}(M_{\odot})$)=9.0-10.9 probed by our sample are very close to circular and have uniform thickness.Comment: 19 figures, 36 pages, accepted for publication in ApJ. Table of PHANGS rotation curves available from http://phangs.org/dat

The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk

Context. Cloud-scale surveys of molecular gas reveal the link between giant molecular cloud properties and star formation across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the normal star formation process, while galaxy centers tend to harbor denser gas that exhibits more extreme star formation. At high resolution, however, molecular clouds with exceptional gas properties and star formation activity may also be observed in normal disk environments. In this paper we study the brightest cloud traced in CO(2-1) emission in the disk of nearby spiral galaxy NGC 628. Aims. We characterize the properties of the molecular and ionized gas that is spatially coincident with an extremely bright H ii region in the context of the NGC 628 galactic environment. We investigate how feedback and large-scale processes influence the properties of the molecular gas in this region. Methods. High-resolution ALMA observations of CO(2-1) and CO(1−0) emission were used to characterize the mass and dynamical state of the 'headlight' molecular cloud. The characteristics of this cloud are compared to the typical properties of molecular clouds in NGC 628. A simple large velocity gradient (LVG) analysis incorporating additional ALMA observations of 13CO(1−0), HCO+(1−0), and HCN(1−0) emission was used to constrain the beam-diluted density and temperature of the molecular gas. We analyzed the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the H ii region. Results. The unusually bright headlight cloud is massive (1−2 x 107 M), with a beam-diluted density of nH2 = 5 x 104 cm−3 based on LVG modeling. It has a low virial parameter, suggesting that the CO emission associated with this cloud may be overluminous due to heating by the H ii region. A young (2−4 Myr) stellar population with mass 3 x105 M is associated. Conclusions. We argue that the headlight cloud is currently being destroyed by feedback from young massive stars. Due to the large mass of the cloud, this phase of the its evolution is long enough for the impact of feedback on the excitation of the gas to be observed. The high mass of the headlight cloud may be related to its location at a spiral co-rotation radius, where gas experiences reduced galactic shear compared to other regions of the disk and receives a sustained inflow of gas that can promote the mass growth of the cloud.CNH, AH and JP acknowledge support from the Programme National “Physique et Chimie du Milieu Interstellaire” (PCMI) of CNRS/INSU with INC/INP co-funded by CEA and CNES, and from the Programme National Cosmology and Galaxies (PNCG) of CNRS/INSU with INP and IN2P3, co-funded by CEA and CNES. AU acknowledges support from the Spanish funding grants AYA2016-79006-P (MINECO/FEDER) and PGC2018-094671-B-I00 (MCIU/AEI/FEDER). The work of AKL, JS, and DU is partially supported by the National Science Foundation under Grants No. 1615105, 1615109, and 1653300. FB acknowledges funding from the European Union’s Horizon 2020 research and innovation programme (grant agreement No 726384). APSH is a fellow of the International Max Planck Research School for Astronomy and Cosmic Physics at the University of Heidelberg (IMPRSHD). SCOG acknowledges support from the DFG via SFB 881 “The Milky Way System” (sub-projects B1, B2 and B8). JMDK gratefully acknowledges funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme via the ERC Starting Grant MUSTANG (grant agreement number 714907). JMDK and MC gratefully acknowledge funding from the German Research Foundation (DFG) in the form of an Emmy Noether Research Group (grant number KR4801/1-1). SEM acknowledges funding during part of this work from the Deutsche Forschungsgemeinschaft (DFG) via grant SCHI 536/7-2 as part of the priority program SPP 1573 “ISM-SPP: Physics of the Interstellar Medium”

The lifecycle of molecular clouds in nearby star-forming disc galaxies

It remains a major challenge to derive a theory of cloud-scale (⁠≲100 pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially resolved (∼100 pc) CO-to-H α flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10−30Myr⁠, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities ΣH2≥8M⊙pc−2⁠, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at ΣH2≤8M⊙pc−2 GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H α (75-90 per cent of the cloud lifetime), GMCs disperse within just 1−5Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4-10 per cent. These results show that galactic star formation is governed by cloud-scale, environmentally dependent, dynamical processes driving rapid evolutionary cycling. GMCs and H II regions are the fundamental units undergoing these lifecycles, with mean separations of 100−300pc in star-forming discs. Future work should characterize the multiscale physics and mass flows driving these lifecycles.MC and JMDK gratefully acknowledge funding from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) through an Emmy Noether Research Group (grant number KR4801/1-1) and the DFG Sachbeihilfe (grant number KR4801/2-1). JMDK, APSH, SMRJ, and DTH gratefully acknowledge funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme via the ERC Starting Grant MUSTANG (grant agreement number 714907). MC, JMDK, SMRJ, and DTH acknowledge support from the Australia-Germany Joint Research Cooperation Scheme (UA-DAAD, grant number 57387355). APSH, SMRJ, and DTH are fellows of the International Max Planck Research School for Astronomy and Cosmic Physics at the University of Heidelberg (IMPRS-HD). BG gratefully acknowledges the support of the Australian Research Council as the recipient of a Future Fellowship (FT140101202). CNC, AH, and JP acknowledge funding from the Programme National ‘Physique et Chimie du Milieu Interstellaire’ (PCMI) of the Centre national de la recherche scientifique/Institut national des sciences de l’Univers (CNRS/INSU) with the Institut de Chimie/Institut de Physique (INC/INP), co-funded by the Commissariat a l’ ` energie ´ atomique et aux energies alternatives (CEA) and the Centre ´ national d’etudes spatiales (CNES). AH acknowledges support ´ by the Programme National Cosmology et Galaxies (PNCG) of CNRS/INSU with the INP and the Institut national de physique nucleaire et de physique des particules (IN2P3), co-funded by ´ CEA and CNES. PL, ES, CF, DL, and TS acknowledge funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 694343). The work of AKL, JS, and DU is partially supported by the National Science Foundation (NSF) under Grants No. 1615105, 1615109, and 1653300. AKL also acknowledges partial support from the National Aeronautics and Space Administration (NASA) Astrophysics Data Analysis Program (ADAP) grants NNX16AF48G and NNX17AF39G. ER acknowledges the support of the Natural Sciences and Engineering Research Council of Canada (NSERC), funding reference number RGPIN-2017-03987. FB acknowledges funding from the ERC under the European Union’s Horizon 2020 research and innovation programme (grant agreement No. 726384). GB is supported by the Fondo de Fomento al Desarrollo Cient´ıfico y Tecnologico of the Comisi ´ on Nacional de ´ Investigacion Cient ´ ´ıfica y Tecnologica (CONICYT/FONDECYT), ´ Programa de Iniciacion, Folio 11150220. SCOG acknowledges ´ support from the DFG via SFB 881 ‘The Milky Way System’ (subprojects B1, B2, and B8) and also via Germany’s Excellence Strategy EXC-2181/1–390900948 (the Heidelberg STRUCTURES Excellence Cluster). KK gratefully acknowledges funding from the DFG in the form of an Emmy Noether Research Group (grant number KR4598/2-1, PI Kreckel). AU acknowledges support from the Spanish funding grants AYA2016-79006-P (MINECO/FEDER) and PGC2018-094671-B-I00 (MCIU/AEI/FEDER)

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

We measure the velocity dispersion, $\sigma$, and surface density, $\Sigma$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $\Sigma$ and $\sigma$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $\alpha_\mathrm{vir}\propto\sigma^2/\Sigma$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\propto\Sigma\,\sigma^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $\Sigma$, $\sigma$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $\sigma$-$\Sigma$ relation flattens at $\Sigma\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $\sigma$ for a given $\Sigma$ and high apparent $\alpha_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.Comment: Accepted for publication in ApJ. 45 pages, 11 figures, 8 tables, 4 Appendices; key results summarized in Figure 10. Machine-readable table can be downloaded at http://www.astronomy.ohio-state.edu/~sun.1608/datafile3.txt prior to publication. For a brief video describing the main results of this paper, please see https://www.youtube.com/watch?v=-_eL7t1PVq8&

Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-forming Galaxies

We compare the observed turbulent pressure in molecular gas, P_(turb), to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, P_(DE). To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multiwavelength data that trace the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that P_(turb) correlates with—but almost always exceeds—the estimated P_(DE) on kiloparsec scales. This indicates that the molecular gas is overpressurized relative to the large-scale environment. We show that this overpressurization can be explained by the clumpy nature of molecular gas; a revised estimate of P_(DE) on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed P_(turb) in galaxy disks. We also find that molecular gas with cloud-scale P_(turb) ≈ P_(DE) ≳ 10⁵ kB K cm⁻³ in our sample is more likely to be self-gravitating, whereas gas at lower pressure it appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between P_(turb) and the observed SFR surface density, Σ_(SFR), is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between Σ_(SFR) and kpc-scale P_(DE) in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale P_(DE) reported in previous works

The lifecycle of molecular clouds in nearby star-forming disc galaxies

It remains a major challenge to derive a theory of cloud-scale (⁠≲100 pc) star formation and feedback, describing how galaxies convert gas into stars as a function of the galactic environment. Progress has been hampered by a lack of robust empirical constraints on the giant molecular cloud (GMC) lifecycle. We address this problem by systematically applying a new statistical method for measuring the evolutionary timeline of the GMC lifecycle, star formation, and feedback to a sample of nine nearby disc galaxies, observed as part of the PHANGS-ALMA survey. We measure the spatially resolved (∼100 pc) CO-to-H α flux ratio and find a universal de-correlation between molecular gas and young stars on GMC scales, allowing us to quantify the underlying evolutionary timeline. GMC lifetimes are short, typically 10−30 Myr⁠, and exhibit environmental variation, between and within galaxies. At kpc-scale molecular gas surface densities Σ_(H₂) ≥ 8 M_⊙ pc⁻²⁠, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at Σ_(H₂) ≤ 8 M_⊙ pc⁻² GMCs decouple from galactic dynamics and live for an internal dynamical time-scale. After a long inert phase without massive star formation traced by H α (75–90 per cent of the cloud lifetime), GMCs disperse within just 1−5 Myr once massive stars emerge. The dispersal is most likely due to early stellar feedback, causing GMCs to achieve integrated star formation efficiencies of 4–10 per cent. These results show that galactic star formation is governed by cloud-scale, environmentally dependent, dynamical processes driving rapid evolutionary cycling. GMCs and H II regions are the fundamental units undergoing these lifecycles, with mean separations of 100−300 pc in star-forming discs. Future work should characterize the multiscale physics and mass flows driving these lifecycles

Molecular Gas Properties on Cloud Scales across the Local Star-forming Galaxy Population

Using the PHANGS–ALMA CO(2–1) survey, we characterize molecular gas properties on ~100 pc scales across 102,778 independent sightlines in 70 nearby galaxies. This yields the best synthetic view of molecular gas properties on cloud scales across the local star-forming galaxy population obtained to date. Consistent with previous studies, we observe a wide range of molecular gas surface densities (3.4 dex), velocity dispersions (1.7 dex), and turbulent pressures (6.5 dex) across the galaxies in our sample. Under simplifying assumptions about subresolution gas structure, the inferred virial parameters suggest that the kinetic energy of the molecular gas typically exceeds its self-gravitational binding energy at ~100 pc scales by a modest factor (1.3 on average). We find that the cloud-scale surface density, velocity dispersion, and turbulent pressure (1) increase toward the inner parts of galaxies, (2) are exceptionally high in the centers of barred galaxies (where the gas also appears less gravitationally bound), and (3) are moderately higher in spiral arms than in inter-arm regions. The galaxy-wide averages of these gas properties also correlate with the integrated stellar mass, star formation rate, and offset from the star-forming main sequence of the host galaxies. These correlations persist even when we exclude regions with extraordinary gas properties in galaxy centers, which contribute significantly to the inter-galaxy variations. Our results provide key empirical constraints on the physical link between molecular cloud populations and their galactic environment