Source author record

Adam K. Leroy

Adam K. Leroy appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph astro-ph.IM

Catalog footprint

What is connected

76works

4topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

Evolved Supergiants in PHANGS I: Red Supergiants in 19 Galaxies between 5-20 Mpc with HST and JWST

Red supergiants (RSGs) are important for our understanding of supernova progenitors, stellar populations, stellar evolution, mass loss and dust production. Extragalactic surveys of RSGs have a long history in the Local Group, but few studies exist beyond that due to the limited resolution and sensitivity of ground-based and previous space-based infrared observatories. Here we demonstrate the combined power of HST and JWST to push systematic searches of RSGs out to $\sim$20 Mpc. We introduce a catalog of 97057 RSGs -- the largest single-survey release of RSGs -- with masses $\gtrsim$10 M$_{\odot}$ in 19 galaxies from the PHANGS HST+JWST Treasury program. We use HST F814W and JWST F200W photometry to select stars as RSGs based on predicted colors and magnitudes from PARSEC isochrones. The spatial distribution of our recovered RSGs follow the familiar pattern of mostly being concentrated in active star-forming regions such as spiral arms and central starburst rings. The RSG number density on kpc-scales is strongly correlated ($r_s$$\sim$0.82) with local star-formation rate density ($Σ_{SFR}$) traced by extinction-corrected far-ultraviolet (FUV) from GALEX+WISE, and weakly correlated ($r_s$$\sim$0.57) with the total stellar mass density ($Σ_*$), traced by near-infrared emission from WISE+Spitzer. The number of RSGs per mass of stellar populations with ages 6-30 Myr (the likely age range of RSGs $>$10 M$_{\odot}$) is $\sim$1 per 10$^{3.77\pm0.27}$ M$_{\odot}$, assuming constant star-formation rates from FUV+W4. Our sample will be a useful resource for tracking progenitors and feedback sites of future supernovae in PHANGS, age-dating stellar populations, and more.

preprint2022arXiv

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

The complex physical, kinematic, and chemical properties of galaxy centres make them interesting environments to examine with molecular line emission. We present new $2-4$" (${\sim}75{-}150$ pc at $7.7$ Mpc) observations at 2 and 3 mm covering the central $50$" (${\sim}1.9$ kpc) of the nearby double-barred spiral galaxy NGC 6946 obtained with the IRAM Plateau de Bure Interferometer. We detect spectral lines from ten molecules: CO, HCN, HCO$^+$, HNC, CS, HC$_3$N, N$_2$H$^+$, C$_2$H, CH$_3$OH, and H$_2$CO. We complemented these with published 1mm CO observations and 33 GHz continuum observations to explore the star formation rate surface density ${Σ_{\mathrm{SFR}}}$ on 150 pc scales. In this paper, we analyse regions associated with the inner bar of NGC 6946 $-$ the nuclear region (NUC), the northern (NBE), and southern inner bar end (SBE) and we focus on short-spacing corrected bulk (CO) and dense gas tracers (HCN, HCO$^+$, and HNC). We find that HCO$^+$ correlates best with ${Σ_{\mathrm{SFR}}}$, but the dense gas fraction ($f_{\mathrm{dense}}$) and star formation efficiency of the dense gas (${\mathrm{SFE_{dense}}}$) fits show different behaviours than expected from large-scale disc observations.The SBE has a higher ${Σ_{\mathrm{SFR}}}$, $f_{\mathrm{dense}}$, and shocked gas fraction than the NBE. We examine line ratio diagnostics and find a higher CO(2-1)/CO(1-0) ratio towards NBE than for the NUC. Moreover, comparison with existing extragalactic datasets suggests that using the HCN/HNC ratio to probe kinetic temperatures is not suitable on kiloparsec and sub-kiloparsec scales in extragalactic regions. Lastly, our study shows that the HCO$^+$/HCN ratio might not be a unique indicator to diagnose AGN activity in galaxies.

preprint2022arXiv

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

We present the CO isotopologue Line Atlas within the Whirpool galaxy Survey (CLAWS) based on an IRAM 30-m large programme which provides a benchmark study of numerous, faint CO isotopologues in the mm-wavelength regime across the full disc of M51 (NGC 5194). The survey's core goal is to use the low-J CO isotopologue lines to constrain CO excitation and chemistry, and therefrom the local physical conditions of the gas. In this survey paper, we describe the CLAWS observing and data reduction strategies. We map the J=1-0 and 2-1 transitions of the CO isotopologues $^{12}$CO,$^{13}$CO, C$^{18}$O and C$^{17}$O, as well as several supplementary lines within the 1 mm and 3 mm window (CN(1-0), CS(2-1), CH$_3$OH(2-1), N$_2$H$^+$(1-0), HC$_3$N(10-9)) at ~1 kpc resolution. A total observation time of 149 h offers unprecedented sensitivity. We use these data to explore several CO isotopologue line ratios in detail, study their radial (and azimuthal) trends and investigate whether changes in line ratios stem from changes in ISM properties such as gas temperatures, densities or chemical abundances. For example, we find negative radial trends for the $^{13}$CO}/$^{12}$CO, C$^{18}$O/$^{12}$CO and C$^{18}$O/$^{13}$CO line ratios in their J=1-0 transitions. We also find variations with local environment, such as higher $^{12}$CO(2-1)/(1-0) or $^{13}$CO/$^{12}$CO(1-0) line ratios in interarm regions compared to spiral arm regions. We propose that these aforementioned variations of CO line ratios are most likely due to a variation of the optical depth, while abundance variations due to selective nucleosynthesis on a galaxy-wide scale could also play a role. We also study the CO spectral line energy distribution (SLED) using archival JCMT $^{12}$CO(3-2) data and find a variation of the SLED shape with local environmental parameters further underlying changes in optical depth, gas temperatures or densities.

preprint2022arXiv

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

We use integral field spectroscopy from the PHANGS-MUSE survey, which resolves the ionised interstellar medium at ${\sim}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 HII regions and then binning the low-surface-brightness areas to achieve significant detections of the key nebular lines. A simple model for the leakage and propagation of ionising radiation from HII regions is able to reproduce the observed distribution of H$α$ in the DIG. Leaking radiation from HII regions also explains the observed decrease in line ratios of low-ionisation species ([SII]/H$α$, [NII]/H$α$ and [OI]/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 galaxies; (2) the observed trends of a flat or decreasing [OIII]/H$β$ with $Σ_{Hα}$; and (3) the offset of some DIG regions from the locus of HII regions in the Baldwin-Phillips-Terlevich (BPT) diagram, extending into the area of low-ionisation (nuclear) emission-line regions (LI[N]ERs). 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 [OIII] 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.

preprint2022arXiv

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

The processes of star formation and feedback, regulating the cycle of matter between gas and stars on the scales of giant molecular clouds (GMCs; $\sim$100pc), play a major role in governing galaxy evolution. Measuring the time-scales of GMC evolution is important to identify and characterise the specific physical mechanisms that drive this transition. By applying a robust statistical method to high-resolution CO and narrow-band H$α$ imaging from the PHANGS survey, we systematically measure the evolutionary timeline from molecular clouds to exposed young stellar regions on GMC scales, across the discs of an unprecedented sample of 54 star-forming main-sequence galaxies (excluding their unresolved centres). We find that clouds live for about $1{-}3$ GMC turbulence crossing times ($5{-}30$Myr) and are efficiently dispersed by stellar feedback within $1{-}5$Myr once the star-forming region becomes partially exposed, resulting in integrated star formation efficiencies of $1{-}8$%. These ranges reflect physical galaxy-to-galaxy variation. In order to evaluate whether galactic environment influences GMC evolution, we correlate our measurements with average properties of the GMCs and their local galactic environment. We find several strong correlations that can be physically understood, revealing a quantitative link between galactic-scale environmental properties and the small-scale GMC evolution. Notably, the measured CO-visible cloud lifetimes become shorter with decreasing galaxy mass, mostly due to the increasing presence of CO-dark molecular gas in such environment. Our results represent a first step towards a comprehensive picture of cloud assembly and dispersal, which requires further extension and refinement with tracers of the atomic gas, dust, and deeply-embedded stars.

preprint2022arXiv

The Morpho-Kinematic Architecture of Super Star Clusters in the Center of NGC253

The center of the nearby galaxy NGC\,253 hosts a population of more than a dozen super star clusters (SSCs) which are still in the process of forming. The majority of the star formation of the burst is concentrated in these SSCs, and the starburst is powering a multiphase outflow from the galaxy. In this work, we measure the 350~GHz dust continuum emission towards the center of NGC\,253 at 47~milliarcsecond (0.8 pc) resolution using data from the Atacama Large Millimeter/submillimeter Array (ALMA). We report the detection of 350~GHz (dust) continuum emission in the outflow for the first time, associated with the prominent South-West streamer. In this feature, the dust emission has a width of $\approx$~8~pc, is located at the outer edge of the CO emission, and corresponds to a molecular gas mass of $\sim~(8-17)\times10^6$~M$_\odot$. In the starburst nucleus, we measure the resolved radial profiles, sizes, and molecular gas masses of the SSCs. Compared to previous work at somewhat lower spatial resolution, the SSCs here break apart into smaller substructures with radii $0.4-0.7$~pc. In projection, the SSCs, dust, and dense molecular gas appear to be arranged as a thin, almost linear, structure roughly 155~pc in length. The morphology and kinematics of this structure can be well explained as gas following $x_2$ orbits at the center of a barred potential. We constrain the morpho-kinematic arrangement of the SSCs themselves, finding that an elliptical, angular momentum-conserving ring is a good description of the both morphology and kinematics of the SSCs.

preprint2022arXiv

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

We present the PHANGS-MUSE survey, a programme using the MUSE IFS at the ESO VLT to map 19 massive $(9.4 < \log(M_{*}/M_\odot) < 11.0)$ nearby (D < 20 Mpc) star-forming disc galaxies. The survey consists of 168 MUSE pointings (1'x1' each), a total of nearly 15 Million spectra, covering ~1.5 Million independent spectra. PHANGS-MUSE provides the first IFS view of star formation across different local environments (including galaxy centres, bars, spiral arms) in external galaxies at a median resolution of 50~pc, better than the mean inter-cloud distance in the ionised interstellar medium. This `cloud-scale' resolution allows detailed demographics and characterisations of HII regions and other ionised nebulae. PHANGS-MUSE further delivers a unique view on the associated gas and stellar kinematics, and provides constraints on the star formation history. The PHANGS-MUSE survey is complemented by dedicated ALMA CO(2-1) and multi-band HST observations, therefore allowing us to probe the key stages of the star formation process from molecular clouds to HII regions and star clusters. This paper describes the scientific motivation, sample selection, observational strategy, data reduction and analysis process of the PHANGS-MUSE survey. We present our bespoke automated data-reduction framework, which is built on the reduction recipes provided by ESO, but additionally allows for mosaicking and homogenisation of the point spread function. We further present a detailed quality assessment and a brief illustration of the potential scientific applications of the large set of PHANGS-MUSE data products generated by our data analysis framework. The data cubes and analysis data products described in this paper represent the basis for the first PHANGS-MUSE public data release and are available in the ESO archive and via the Canadian Astronomy Data Centre.

preprint2021arXiv

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

We present improved methods for segmenting CO emission from galaxies into individual molecular clouds, providing an update to the CPROPS algorithms presented by Rosolowsky & Leroy (2006; arXiv:astro-ph/0601706 ). The new code enables both homogenization of the noise and spatial resolution among data, which allows for rigorous comparative analysis. The code also models the completeness of the data via false source injection and includes an updated segmentation approach to better deal with blended emission. These improved algorithms are implemented in a publicly available python package, PYCPROPS. We apply these methods to ten of the nearest galaxies in the PHANGS-ALMA survey, cataloguing CO emission at a common 90 pc resolution and a matched noise level. We measure the properties of 4986 individual clouds identified in these targets. We investigate the scaling relations among cloud properties and the cloud mass distributions in each galaxy. The physical properties of clouds vary among galaxies, both as a function of galactocentric radius and as a function of dynamical environment. Overall, the clouds in our target galaxies are well-described by approximate energy equipartition, although clouds in stellar bars and galaxy centres show elevated line widths and virial parameters. The mass distribution of clouds in spiral arms has a typical mass scale that is 2.5x larger than interarm clouds and spiral arms clouds show slightly lower median virial parameters compared to interarm clouds (1.2 versus 1.4).

preprint2021arXiv

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

We measure the low-J CO line ratio R21=CO(2-1)/CO(1-0), R32=CO(3-2)/CO(2-1), and R31 = CO(3-2)/CO(1-0) using whole-disk CO maps of nearby galaxies. We draw CO(2-1) from PHANGS--ALMA, HERACLES, and follow-up IRAM surveys; CO(1-0) from COMING and the Nobeyama CO Atlas of Nearby Spiral Galaxies; and CO(3-2) from the JCMT NGLS and APEX LASMA mapping. Altogether this yields 76, 47, and 29 maps of R21, R32, and R31 at 20" \sim 1.3 kpc resolution, covering 43, 34, and 20 galaxies. Disk galaxies with high stellar mass, log10 M_* [Msun]=10.25-11 and star formation rate, SFR=1-5 Msun/yr, dominate the sample. We find galaxy-integrated mean values and 16%-84% range of R21 = 0.65 (0.50-0.83), R32=0.50 (0.23-0.59), and R31=0.31 (0.20-0.42). We identify weak trends relating galaxy-integrated line ratios to properties expected to correlate with excitation, including SFR/M_* and SFR/L_CO. Within galaxies, we measure central enhancements with respect to the galaxy-averaged value of \sim 0.18^{+0.09}_{-0.14} dex for R21, 0.27^{+0.13}_{-0.15} dex for R31, and 0.08^{+0.11}_{-0.09} dex for R32. All three line ratios anti-correlate with galactocentric radius and positively correlate with the local star formation rate surface density and specific star formation rate, and we provide approximate fits to these relations. The observed ratios can be reasonably reproduced by models with low temperature, moderate opacity, and moderate densities, in good agreement with expectations for the cold ISM. Because the line ratios are expected to anti-correlate with the CO(1-0)-to-H_2 conversion factor, alphaCO^(1-0), these results have general implications for the interpretation of CO emission from galaxies.

preprint2021arXiv

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{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 $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line 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\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$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\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(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 $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.

preprint2021arXiv

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

The sensitivity and angular resolution of photometric surveys executed by the Hubble Space Telescope (HST) enable studies of individual star clusters in galaxies out to a few tens of megaparsecs. The fitting of spectral energy distributions (SEDs) of star clusters is essential for measuring their physical properties and studying their evolution. We report on the use of the publicly available Code Investigating GALaxy Emission (CIGALE) SED fitting package to derive ages, stellar masses, and reddenings for star clusters identified in the Physics at High Angular resolution in Nearby GalaxieS-HST (PHANGS-HST) survey. Using samples of star clusters in the galaxy NGC 3351, we present results of benchmark analyses performed to validate the code and a comparison to SED fitting results from the Legacy ExtraGalactic Ultraviolet Survey (LEGUS). We consider procedures for the PHANGS-HST SED fitting pipeline, e.g., the choice of single stellar population models, the treatment of nebular emission and dust, and the use of fluxes versus magnitudes for the SED fitting. We report on the properties of clusters in NGC 3351 and find, on average, the clusters residing in the inner star-forming ring of NGC 3351 are young ($< 10$ Myr) and massive ($10^{5} M_{\odot}$) while clusters in the stellar bulge are significantly older. Cluster mass function fits yield $β$ values around -2, consistent with prior results with a tendency to be shallower at the youngest ages. Finally, we explore a Bayesian analysis with additional physically-motivated priors for the distribution of ages and masses and analyze the resulting cluster distributions.

preprint2021arXiv

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

We use an unprecedented sample of about 23,000 HII regions detected at an average physical resolution of 67pc in the PHANGS-MUSE sample to study the extragalactic HII region Ha luminosity function (LF). Our observations probe the star-forming disk of 19 nearby spiral galaxies with low inclination and located close to the star formation main sequence at z=0. The mean LF slope $α$ in our sample is =1.73 with a $σ$ of 0.15. We find that $α$ decreases with the galaxy's star formation rate surface density and argue that this is driven by an enhanced clustering of young stars at high gas surface densities. Looking at the HII regions within single galaxies we find that no significant variations occur between the LF of the inner and outer part of the star-forming disk, whereas the LF in the spiral arm areas is shallower than in the inter-arm areas for six out of the 13 galaxies with clearly visible spiral arms. We attribute these variations to the spiral arms increasing the molecular clouds arm--inter-arm mass contrast and find suggestive evidence that they are more evident for galaxies with stronger spiral arms. Furthermore, we find systematic variations in $α$ between samples of HII regions with high and low ionization parameter q and argue that they are driven by the aging of HII regions.

preprint2021arXiv

Resolving the Dust-to-Metals Ratio and CO-to-H$_2$ Conversion Factor in the Nearby Universe

We investigate the relationship between the dust-to-metals ratio (D/M) and the local interstellar medium environment at ~2 kpc resolution in five nearby galaxies: IC342, M31, M33, M101, and NGC628. A modified blackbody model with a broken power-law emissivity is used to model the dust emission from 100 to 500 um observed by Herschel. We utilize the metallicity gradient derived from auroral line measurements in HII regions whenever possible. Both archival and new CO rotational line and HI 21 cm maps are adopted to calculate gas surface density, including new wide field CO and HI maps for IC342 from IRAM and the VLA, respectively. We experiment with several prescriptions of CO-to-H$_2$ conversion factor, and compare the resulting D/M-metallicity and D/M-density correlations, both of which are expected to be non-negative from depletion studies. The D/M is sensitive to the choice of the conversion factor. The conversion factor prescriptions based on metallicity only yield too much molecular gas in the center of IC342 to obtain the expected correlations. Among the prescriptions tested, the one that yields the expected correlations depends on both metallicity and surface density. The 1-$σ$ range of the derived D/M spans 0.40-0.58. Compared to chemical evolution models, our measurements suggest that the dust growth time scale is much shorter than the dust destruction time scale. The measured D/M is consistent with D/M in galaxy-integrated studies derived from infrared dust emission. Meanwhile, the measured D/M is systematically higher than the D/M derived from absorption, which likely indicates a systematic offset between the two methods.

preprint2021arXiv

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

The PHANGS program is building the first dataset to enable the multi-phase, multi-scale study of star formation across the nearby spiral galaxy population. This effort is enabled by large survey programs with ALMA, VLT/MUSE, and HST, with which we have obtained CO(2-1) imaging, optical spectroscopic mapping, and high resolution UV-optical imaging, respectively. Here, we present PHANGS-HST, which has obtained five band NUV-U-B-V-I imaging of the disks of 38 spiral galaxies at distances of 4-23 Mpc, and parallel V and I band imaging of their halos, to provide a census of tens of thousands of compact star clusters, and multi-scale stellar associations. The combination of HST, ALMA, and VLT/MUSE observations will yield an unprecedented joint catalog of the observed and physical properties of ~100,000 star clusters, associations, HII regions, and molecular clouds. With these basic units of star formation, PHANGS will systematically chart the evolutionary cycling between gas and stars, across a diversity of galactic environments found in nearby galaxies. We discuss the design of the PHANGS-HST survey, and provide an overview of the HST data processing pipeline and first results, highlighting new methods for selecting star cluster candidates, morphological classification of candidates with convolutional neural networks, and identification of stellar associations over a range of physical scales with a watershed algorithm. We describe the cross-observatory imaging, catalogs, and software products to be released. These high-level science products will seed a broad range of investigations, in particular, the study of embedded stellar populations and dust with JWST, for which a PHANGS Cycle 1 Treasury program to obtain eight band 2-21 $μ$m imaging has been approved.

preprint2020arXiv

A model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling

In Meidt et al. (2018), we showed that gas kinematics on the scale of individual molecular clouds are not dominated by self-gravity but also track a component that originates with orbital motion in the potential of the host galaxy. This agrees with observed cloud line widths, which show systematic variations from virial motions with environment, pointing at the influence of the galaxy potential. In this paper, we hypothesize that these motions act to slow down the collapse of gas and so help regulate star formation. Extending the results of Meidt et al. (2018), we derive a dynamical collapse timescale that approaches the free-fall time only once the gas has fully decoupled from the galactic potential. Using this timescale we make predictions for how the fraction of free-falling, strongly self-gravitating gas varies throughout the disks of star-forming galaxies. We also use this collapse timescale to predict variations in the molecular gas star formation efficiency, which is lowered from a maximum, feedback-regulated level in the presence of strong coupling to the galactic potential. Our model implies that gas can only decouple from the galaxy to collapse and efficiently form stars deep within clouds. We show that this naturally explains the observed drop in star formation rate per unit gas mass in the Milky Way's CMZ and other galaxy centers. The model for a galactic bottleneck to star formation also agrees well with resolved observations of dense gas and star formation in galaxy disks and the properties of local clouds.

preprint2020arXiv

Detection of the diffuse HI emission in the Circumgalactic Medium of NGC 891 and NGC 4565

We present detections of 21-cm emission from neutral hydrogen (HI) in the circumgalactic medium (CGM) of the local edge-on galaxies NGC 891 and NGC 4565 using the Robert C. Byrd Green Bank Telescope (GBT). With our 5$σ$ sensitivity of $8.2 \times 10^{16}$ cm$^{-2}$ calculated over a 20 km s$^{-1}$ channel, we achieve $>5σ$ detections out to $90-120$ kpc along the minor axes. The velocity width of the CGM emission is as large as that of the disk $\approx 500$ km s$^{-1}$, indicating the existence of a diffuse component permeating the halo. We compare our GBT measurements with interferometric data from the Westerbork Synthesis Radio Telescope (WSRT). The WSRT maps the HI emission from the disk at high S/N but has limited surface brightness sensitivity at the angular scales probed with the GBT. After convolving the WSRT data to the spatial resolution of the GBT (FWHM = 9.1$'$), we find that the emission detected by the WSRT accounts for $48^{+15}_{-25}$% ($58^{+4}_{-18}$%) of the total flux recovered by the GBT from the CGM of NGC 891(NGC 4565). The existence of significant GBT-only flux suggests the presence of a large amount of diffuse, low column density HI emission in the CGM. For reasonable assumptions, the extended diffuse HI could account for $5.2\pm0.9$% and $2.0\pm0.8$% of the total HI emission of NGC 891 and NGC 4565.

preprint2020arXiv

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

PHANGS-HST is an ultraviolet-optical imaging survey of 38 spiral galaxies within ~20 Mpc. Combined with the PHANGS-ALMA, PHANGS-MUSE surveys and other multiwavelength data, the dataset will provide an unprecedented look into the connections between young stars, HII regions, and cold molecular gas in these nearby star-forming galaxies. Accurate distances are needed to transform measured observables into physical parameters (e.g., brightness to luminosity, angular to physical sizes of molecular clouds, star clusters and associations). PHANGS-HST has obtained parallel ACS imaging of the galaxy halos in the F606W and F814W bands. Where possible, we use these parallel fields to derive tip of the red giant branch (TRGB) distances to these galaxies. In this paper, we present TRGB distances for 11 galaxies from ~4 to ~15 Mpc, based on the first year of PHANGS-HST observations. Five of these represent the first published TRGB distance measurements (IC 5332, NGC 2835, NGC 4298, NGC 4321, and NGC 4328), and eight of which are the best available distances to these targets. We also provide a compilation of distances for the 118 galaxies in the full PHANGS sample, which have been adopted for the first PHANGS-ALMA public data release.

preprint2020arXiv

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

We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $Σ_\mathrm{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 $Σ_\mathrm{SFR}$ and kpc-scale $P_\mathrm{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_\mathrm{DE}$ reported in previous works.

preprint2020arXiv

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.

preprint2020arXiv

Spatial Power Spectra of Dust across the Local Group: No Constraint on Disc Scale Height

We analyze the 1D spatial power spectra of dust surface density and mid to far-infrared emission at $24-500\,μ$m in the LMC, SMC, M31, and M33. By forward-modelling the point-spread-function (PSF) on the power spectrum, we find that nearly all power spectra have a single power-law and point source component. A broken power-law model is only favoured for the LMC 24 $μ$m MIPS power spectrum and is due to intense dust heating in 30 Doradus. We also test for local power spectrum variations by splitting the LMC and SMC maps into $820$ pc boxes. We find significant variations in the power-law index with no strong evidence for breaks. The lack of a ubiquitous break suggests that the spatial power spectrum does not constrain the disc scale height. This contradicts claims of a break where the turbulent motion changes from 3D to 2D. The power spectrum indices in the LMC, SMC, and M31 are similar ($2.0-2.5$). M33 has a flatter power spectrum ($1.3$), similar to more distant spiral galaxies with a centrally-concentrated H$_{2}$ distribution. We compare the power spectra of HI, CO, and dust in M31 and M33 and find that HI power spectra are consistently flatter than CO power spectra. These results cast doubt on the idea that the spatial power spectrum traces large scale turbulent motion in nearby galaxies. Instead, we find that the spatial power spectrum is influenced by (1) the PSF on scales below $\sim3$ times the FWHM, (2) bright compact regions (30 Doradus), and (3) the global morphology of the tracer (an exponential CO disc).

preprint2020arXiv

The evolution of neutral hydrogen over the past 11 Gyr via HI 21 cm absorption

We present the results of a blind search for intervening HI 21 cm absorption toward 260 radio sources in the redshift range 0<z<2.74 with the Green Bank Telescope. The survey has the sensitivity to detect sub-damped Lyman-alpha (DLA) systems for HI spin temperatures $T_s/f$ = 100 K, and despite the successful re-detection of ten known 21 cm absorbers in the sample, we detect no new absorption lines in the full survey. Sources detected in 21 cm absorption were also searched for hydroxyl (OH) 18 cm absorption and we re-detect 1667 MHz OH absorption toward PKS 1830-211. We searched for intervening HI 21 cm absorption along the line of sight in each source achieving a total redshift coverage of $Δz$ = 88.64 (comoving absorption path of $ΔX$ = 159.5) after removing regions affected by radio frequency interference. We compute a 95% confidence upper limit on the column density frequency distribution $f(N_{\rm HI})$ and set a statistical constraint on the spin temperature $T_s$ in the range 100-1000 K, consistent with prior redshifted optical DLA surveys and HI 21 cm emission observations at the same redshifts. We infer a value for the cosmological mass density of neutral gas, $Ω_{\rm HI}$. Through comparison with prior $Ω_{\rm HI}$ measurements, we place a statistical constraint on the mean spin temperature of $T_s/f$ = 175 K. Our derived $Ω_{\textrm {HI}}$ values support a relative mild evolution in $Ω_{\textrm {HI}}$ over the last 11 Gyr and are consistent with other methods that measure $Ω_{\textrm {HI}}$.

preprint2020arXiv

The Molecular ISM in the Super Star Clusters of the Starburst NGC 253

We present submillimeter spectra of the (proto-)super star cluster (SSC) candidates in the starbursting center of the nearby galaxy NGC 253 identified by Leroy et al. (2018). The 2.5pc resolution of our ALMA cycle 3 observations approach the size of the SSCs and allows the study of physical and chemical properties of the molecular gas in these sources. In the 14 SSC sources and in the frequency ranges 342.0-345.8 GHz and 353.9-357.7 GHz we detect 55 lines belonging to 19 different chemical species. The SSCs differ significantly in chemical complexity, with the richest clusters showing 19 species and the least complex showing 4 species. We detect HCN isotopologues and isomers (H$^{13}$CN, HC$^{15}$N, H$^{15}$NC), abundant HC$_3$N, SO and S$^{18}$O, SO$_2$, and H$_2$CS. The gas ratios CO/HCN, CO/HCO$^+$ are low, ~1-10, implying high dense gas fractions in the SSCs. Line ratio analyses suggests chemistry consistent with photon-dominated regions and mechanical heating. None of the SSCs near the galaxy center show line ratios that imply an X-ray dominated region, suggesting that heating by any (still unknown) AGN does not play a major role. The gas temperatures are high in most sources, with an average rotational temperature of ~130 K in SO$_2$. The widespread existence of vibrationally excited HCN and HC$_3$N transitions implies strong IR radiation fields, potentially trapped by a greenhouse effect due to high continuum opacities.

preprint2020arXiv

The turbulent gas structure in the centers of NGC253 and the Milky Way

We compare molecular gas properties in the starbursting center of NGC253 and the Milky Way Galactic Center (GC) on scales of ~1-100 pc using dendograms and resolution-, area- and noise-matched datasets in CO (1-0) and CO (3-2). We find that the size-line width relations in NGC253 and the GC have similar slope, but NGC253 has larger line widths by factors of ~2-3. The $σ^2/R$ dependency on column density shows that, in the GC, on scales of 10-100 pc the kinematics of gas over $N>3\times10^{21}$ cm$^{-2}$ are compatible with gravitationally bound structures. In NGC253 this is only the case for column densities $N>3\times10^{22}$ cm$^{-2}$. The increased line widths in NGC253 originate in the lower column density gas. This high-velocity dispersion, not gravitationally self-bound gas is likely in transient structures created by the combination of high average densities and feedback in the starburst. The high densities turns the gas molecular throughout the volume of the starburst, and the injection of energy and momentum by feedback significantly increases the velocity dispersion at a given spatial scale over what is observed in the GC.

preprint2019arXiv

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

Cloud-scale surveys of molecular gas reveal the link between molecular clouds properties and star formation (SF) across a range of galactic environments. Cloud populations in galaxy disks are considered to be representative of the `normal' SF. At high resolution, however, clouds with exceptional gas properties and SF activity may also be observed in normal disk environments. In this paper, we study the brightest cloud traced in CO emission in the disk of NGC628. The cloud is spatially coincident with an extremely bright HII region. We characterize its molecular gas properties and investigate how feedback and large-scale processes influence the properties of the molecular gas. High resolution CO ALMA observations are used to characterize its mass and dynamical state, which are compared to other clouds in NGC628. A LVG analysis is used to constrain the beam-diluted density and temperature of the molecular gas. We analyze the MUSE spectrum using Starburst99 to characterize the young stellar population associated with the HII region. The cloud is massive ($1-2\times10^7$M$_{\odot}$), with a beam-diluted density of $n_{\rm H_2}=5\times10^4$ cm$^{-3}$. It has a low virial parameter, suggesting that its CO emission may be overluminous due to heating by the HII region. A young ($2-4$ Myr), massive $3\times10^{5}$ M$_{\odot}$ stellar population is associated. We argue that the cloud is currently being destroyed by feedback from young massive stars. Due to the cloud's large mass, this phase of the cloud's evolution is long enough for the impact of feedback on the excitation of the gas to be observed. Its high mass 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 cloud's mass growth.

preprint2019arXiv

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

It remains a major challenge to derive a theory of cloud-scale ($\lesssim100$ 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 ($\sim100$ 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 $Σ_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $Σ_{\rm H_2}\leqslant8$M$_{\odot}$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% 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% These results show that galactic star formation is governed by cloud-scale, environmentally-dependent, dynamical processes driving rapid evolutionary cycling. GMCs and HII regions are the fundamental units undergoing these lifecycles, with mean separations of 100-300 pc in star-forming discs. Future work should characterise the multi-scale physics and mass flows driving these lifecycles.

preprint2018arXiv

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

We measure the velocity dispersion, $σ$, and surface density, $Σ$, 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. $Σ$ and $σ$ 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 $α_\mathrm{vir}\proptoσ^2/Σ$, 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Σ\,σ^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $Σ$, $σ$, 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 $σ$-$Σ$ relation flattens at $Σ\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $σ$ for a given $Σ$ and high apparent $α_\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.

preprint2016arXiv

A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes

The cloud-scale density, velocity dispersion, and gravitational boundedness of the interstellar medium (ISM) vary within and among galaxies. In turbulent models, these properties play key roles in the ability of gas to form stars. New high fidelity, high resolution surveys offer the prospect to measure these quantities across galaxies. We present a simple approach to make such measurements and to test hypotheses that link small-scale gas structure to star formation and galactic environment. Our calculations capture the key physics of the Larson scaling relations, and we show good correspondence between our approach and a traditional "cloud properties" treatment. However, we argue that our method is preferable in many cases because of its simple, reproducible characterization of all emission. Using, low-J 12CO data from recent surveys, we characterize the molecular ISM at 60pc resolution in the Antennae, the Large Magellanic Cloud, M31, M33, M51, and M74. We report the distributions of surface density, velocity dispersion, and gravitational boundedness at 60pc scales and show galaxy-to-galaxy and intra-galaxy variations in each. The distribution of flux as a function of surface density appears roughly lognormal with a 1sigma width of ~0.3 dex, though the center of this distribution varies from galaxy to galaxy. The 60pc resolution line width and molecular gas surface density correlate well, which is a fundamental behavior expected for virialized or free-falling gas. Varying the measurement scale for the LMC and M31, we show that the molecular ISM has higher surface densities, lower line widths, and more self-gravity at smaller scales.

preprint2016arXiv

Dense Gas Fraction and Star Formation Efficiency Variations in the Antennae Galaxies

We use the CARMA millimeter interferometer to map the Antennae Galaxies (NGC4038/39), tracing the bulk of the molecular gas via the 12CO(1-0) line and denser molecular gas via the high density transitions HCN(1-0), HCO+(1-0), CS(2-1), and HNC(1-0). We detect bright emission from all tracers in both the two nuclei and three locales in the overlap region between the two nuclei. These three overlap region peaks correspond to previously identified "supergiant molecular clouds". We combine the CARMA data with Herschel infrared (IR) data to compare observational indicators of the star formation efficiency (SFR/H2~IR/CO), dense gas fraction (HCN/CO), and dense gas star formation efficiency (IR/HCN). Regions within the Antennae show ratios consistent with those seen for entire galaxies, but these ratios vary by up to a factor of 6 within the galaxy. The five detected regions vary strongly in both their integrated intensities and these ratios. The northern nucleus is the brightest region in mm-wave line emission, while the overlap region is the brightest part of the system in the IR. We combine the CARMA and Herschel data with ALMA CO data to report line ratio patterns for each bright point. CO shows a declining spectral line energy distribution, consistent with previous studies. HCO+(1-0) emission is stronger than HCN(1-0) emission, perhaps indicating either more gas at moderate densities or higher optical depth than is commonly seen in more advanced mergers.

preprint2016arXiv

The Molecular Clouds Fueling a 1/5 Solar Metallicity Starburst

Using the Atacama Large Millimeter/submillimeter Array, we have made the first high spatial and spectral resolution observations of the molecular gas and dust in the prototypical blue compact dwarf galaxy II Zw 40. The CO(2-1) and CO(3-2) emission is clumpy and distributed throughout the central star-forming region. Only one of eight molecular clouds has associated star formation. The continuum spectral energy distribution is dominated by free-free and synchrotron: at 870$μm$, only 50% of the emission is from dust. We derive a CO-to-H$_2$ conversion factor by several methods including a new method that uses simple photodissocation models and resolved CO line intensity measurements to derive a relationship that uniquely predicts $α_{CO}$ for a given metallicity. We find that the CO-to-H$_2$ conversion factor is 4 to 35 times that of the Milky Way (18.1 to 150.5 M$_\odot$ / (K km/s pc$^2$)). The star formation efficiency of the molecular gas at least 10 times higher than that found in normal spiral galaxies, which is likely due to the burst-dominated star formation history of II Zw 40 rather than an intrinsically higher efficiency. The molecular clouds within II Zw 40 resemble those in other strongly interacting systems like the Antennae: overall they have high size-linewidth coefficients and molecular gas surface densities. These properties appear to be due to the high molecular gas surface densities produced in this merging system rather than to increased external pressure. Overall, these results paint a picture of II Zw 40 as a complex, rapidly evolving system whose molecular gas properties are dominated by the large-scale gas shocks from its on-going merger.

preprint2016arXiv

The Relationship Between Molecular Gas, HI, and Star Formation in the Low-Mass, Low-Metallicity Magellanic Clouds

The Magellanic Clouds provide the only laboratory to study the effect of metallicity and galaxy mass on molecular gas and star formation at high (~20 pc) resolution. We use the dust emission from HERITAGE Herschel data to map the molecular gas in the Magellanic Clouds, avoiding the known biases of CO emission as a tracer of H$_{2}$. Using our dust-based molecular gas estimates, we find molecular gas depletion times of ~0.4 Gyr in the LMC and ~0.6 SMC at 1 kpc scales. These depletion times fall within the range found for normal disk galaxies, but are shorter than the average value, which could be due to recent bursts in star formation. We find no evidence for a strong intrinsic dependence of the molecular gas depletion time on metallicity. We study the relationship between gas and star formation rate across a range in size scales from 20 pc to ~1 kpc, including how the scatter in molecular gas depletion time changes with size scale, and discuss the physical mechanisms driving the relationships. We compare the metallicity-dependent star formation models of Ostriker, McKee, and Leroy (2010) and Krumholz (2013) to our observations and find that they both predict the trend in the data, suggesting that the inclusion of a diffuse neutral medium is important at lower metallicity.

preprint2015arXiv

ALMA Multi-line Imaging of the Nearby Starburst Galaxy NGC 253

We present spatially resolved ($\sim$50 pc) imaging of molecular gas species in the central kiloparsec of the nearby starburst galaxy NGC 253, based on observations taken with the Atacama Large Millimeter/submillimeter Array (ALMA). A total of 50 molecular lines are detected over a 13 GHz bandwidth imaged in the 3 mm band. Unambiguous identifications are assigned for 27 lines. Based on the measured high CO/C$^{17}$O isotopic line ratio ($\gtrsim$350), we show that $^{12}$CO(1-0) has moderate optical depths. A comparison of the HCN and HCO$^{+}$ with their $^{13}$C-substituted isotopologues shows that the HCN(1-0) and HCO$^{+}$(1-0) lines have optical depths at least comparable to CO(1-0). H$^{13}$CN/H$^{13}$CO$^{+}$ (and H$^{13}$CN/HN$^{13}$C) line ratios provide tighter constraints on dense gas properties in this starburst. SiO has elevated abundances across the nucleus. HNCO has the most distinctive morphology of all the bright lines, with its global luminosity dominated by the outer parts of the central region. The dramatic variation seen in the HNCO/SiO line ratio suggests that some of the chemical signatures of shocked gas are being erased in the presence of dominating central radiation fields (traced by C$_{2}$H and CN). High density molecular gas tracers (including HCN, HCO$^+$, and CN) are detected at the base of the molecular outflow. We also detect hydrogen $β$ recombination lines that, like their $α$ counterparts, show compact, centrally peaked morphologies, distinct from the molecular gas tracers. A number of sulfur based species are mapped (CS, SO, NS, C$_{2}$S, H$_{2}$CS and CH$_{3}$SH) and have morphologies similar to SiO.

preprint2015arXiv

Faint CO Line Wings in Four Star-Forming (Ultra)luminous Infrared Galaxies

We report the results of a search for large velocity width, low-intensity line wings - a commonly used signature of molecular outflows - in four low redshift (ultra)luminous infrared galaxies (U/LIRGs) that appear to be dominated by star formation. The targets were drawn from a sample of fourteen such galaxies presented in Chung et al. (2011), who showed the stacked CO spectrum of the sample to exhibit 1000 km/s-wide line wings. We obtained sensitive, wide bandwidth imaging of our targets using the IRAM Plateau de Bure Interferometer. We detect each target at very high significance but do not find the claimed line wings in these four targets. Instead, we constrain the flux in the line wings to be only a few percent. Casting our results as mass outflow rates following Cicone et al. (2014) we show them to be consistent with a picture in which very high mass loading factors preferentially occur in systems with high AGN contributions to their bolometric luminosity. We identify one of our targets, IRAS05083 (VII Zw 31), as a candidate molecular outflow.

preprint2015arXiv

Microwave Continuum Emission and Dense Gas Tracers in NGC 3627: Combining Jansky VLA and ALMA Observations

We present Karl G. Jansky Very Large Array (VLA) Ka band (33 GHz) and Atacama Large Millimeter Array (ALMA) Band 3 (94.5 GHz) continuum images covering the nucleus and two extranuclear star-forming regions within the nearby galaxy NGC 3627 (M 66), observed as part of the Star Formation in Radio Survey (SFRS). Both images achieve an angular resolution of $\lesssim$2\arcsec, allowing us to map the radio spectral indices and estimate thermal radio fractions at a linear resolution of $\lesssim$90 pc at the distance of NGC 3627. The thermal fraction at 33 GHz reaches unity at and around the peaks of each HII region; we additionally observed the spectral index between 33 and 94.5 GHz to become both increasingly negative and positive away from the peaks of the HII regions, indicating an increase of non-thermal extended emission from diffusing cosmic-ray electrons and the possible presence of cold dust, respectively. While the ALMA observations were optimized for collecting continuum data, they also detected line emission from the $J=1\rightarrow0$ transitions of HCN and HCO$^{+}$. The peaks of dense molecular gas traced by these two spectral lines are spatially offset from the peaks of the 33 and 94.5 GHz continuum emission for the case of the extranuclear star-forming regions, indicating that our data reach an angular resolution at which one can spatially distinguish sites of recent star formation from the sites of future star formation. Finally, we find trends of decreasing dense gas fraction and velocity dispersion with increasing star formation efficiency among the three regions observed, indicating that the dynamical state of the dense gas, rather than its abundance, plays a more significant role in the star formation process.

preprint2015arXiv

Next Generation Very Large Array Memo No. 7 Science Working Group 2: "Galaxy Ecosystems": The Matter Cycle in and Around Galaxies

This white paper discusses how a "next-generation" Very Large Array (ngVLA) operating in the frequency range 1-116 GHz could be a groundbreaking tool to study the detailed astrophysics of the "matter cycle" in the Milky Way and other galaxies. If optimized for high brightness sensitivity, the ngVLA would bring detailed microwave spectroscopy and modeling of the full radio spectral energy distribution into regular use as survey tools at resolutions of 0.1- 1 arcseconds. This wavelength range includes powerful diagnostics of density, excitation, and chemistry in the cold ISM, as well as multiple tracers of the rate of recent star formation, the magnetic field, shocks, and properties of the ionized ISM. We highlight design considerations that would make this facility revolutionary in this area, the foremost of which is a large amount of collecting area on ~km-length baselines. We also emphasize the strong case for harnessing the large proposed collecting area of the ngVLA for very long baseline applications as part of the core design. This would allow measurements of the three dimensional space motions of galaxies to beyond the Local Group and mapping of the Milky Way out to the far side of the disk. Finally, we discuss the gains from the proposed combination of very high resolution and sensitivity to thermal emission, which include observing the feeding of black holes and resolving forming protoclusters.

preprint2015arXiv

Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)

We describe and execute a novel approach to observationally estimate the lifetimes of giant molecular clouds (GMCs). We focus on the cloud population between the two main spiral arms in M51 (the inter-arm region) where cloud destruction via shear and star formation feedback dominates over formation processes. By monitoring the change in GMC number densities and properties from one side of the inter-arm to the other, we estimate the lifetime as a fraction of the inter-arm travel time. We find that GMC lifetimes in M51's inter-arm are finite and short, 20 to 30 Myr. Such short lifetimes suggest that cloud evolution is influenced by environment, in which processes can disrupt GMCs after a few free-fall times. Over most of the region under investigation shear appears to regulate the lifetime. As the shear timescale increases with galactocentric radius, we expect cloud destruction to switch primarily to star formation feedback at larger radii. We identify a transition from shear- to feedback-dominated disruption through a change in the behavior of the GMC number density. The signature suggests that shear is more efficient at completely dispersing clouds, whereas feedback transforms the population, e.g. by fragmenting high mass clouds into lower mass pieces. Compared to the characteristic timescale for molecular hydrogen in M51, our short lifetimes suggest that gas can remain molecular while clouds disperse and reassemble. We propose that galaxy dynamics regulates the cycling of molecular material from diffuse to bound (and ultimately star-forming) objects, contributing to long observed molecular depletion times in normal disk galaxies. We also speculate that, in more extreme environments such as elliptical galaxies and concentrated galaxy centers, star formation can be suppressed when the shear timescale becomes so short that some clouds can not survive to collapse and form stars.

preprint2015arXiv

The Multi-Phase Cold Fountain in M82 Revealed by a Wide, Sensitive Map of the Molecular ISM

We present a wide area (~ 8 x 8 kpc), sensitive map of CO (2-1) emission around the nearby starburst galaxy M82. Molecular gas extends far beyond the stellar disk, including emission associated with the well-known outflow as far as 3 kpc from M82's midplane. Kinematic signatures of the outflow are visible in both the CO and HI emission: both tracers show a minor axis velocity gradient and together they show double peaked profiles, consistent with a hot outflow bounded by a cone made of a mix of atomic and molecular gas. Combining our CO and HI data with observations of the dust continuum, we study the changing properties of the cold outflow as it leaves the disk. While H_2 dominates the ISM near the disk, the dominant phase of the cool medium changes as it leaves the galaxy and becomes mostly atomic after about a kpc. Several arguments suggest that regardless of phase, the mass in the cold outflow does not make it far from the disk; the mass flux through surfaces above the disk appears to decline with a projected scale length of ~ 1-2 kpc. The cool material must also end up distributed over a much wider angle than the hot outflow based on the nearly circular isophotes of dust and CO at low intensity and the declining rotation velocities as a function of height from the plane. The minor axis of M82 appears so striking at many wavelengths because the interface between the hot wind cavity and the cool gas produces Halpha, hot dust, PAH emission, and scattered UV light. We also show the level at which a face-on version of M82 would be detectable as an outflow based on unresolved spectroscopy. Finally, we consider multiple constraints on the CO-to-H$_2$ conversion factor, which must change across the galaxy but appears to be only a factor of ~ 2 lower than the Galactic value in the outflow.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

We map the distribution of dust in M31 at 25pc resolution, using stellar photometry from the Panchromatic Hubble Andromeda Treasury. We develop a new mapping technique that models the NIR color-magnitude diagram (CMD) of red giant branch (RGB) stars. The model CMDs combine an unreddened foreground of RGB stars with a reddened background population viewed through a log-normal column density distribution of dust. Fits to the model constrain the median extinction, the width of the extinction distribution, and the fraction of reddened stars. The resulting extinction map has >4 times better resolution than maps of dust emission, while providing a more direct measurement of the dust column. There is superb morphological agreement between the new map and maps of the extinction inferred from dust emission by Draine et al. 2014. However, the widely-used Draine & Li (2007) dust models overpredict the observed extinction by a factor of ~2.5, suggesting that M31's true dust mass is lower and that dust grains are significantly more emissive than assumed in Draine et al. (2014). The discrepancy we identify is consistent with similar findings in the Milky Way by the Planck Collaboration (2015), but has a more complex dependence on parameters from the Draine & Li (2007) dust models. We also show that the discrepancy with the Draine et al. (2014) map is lowest where the interstellar radiation field has a harder spectrum than average. We discuss possible improvements to the CMD dust mapping technique, and explore further applications.

preprint2015arXiv

The Panchromatic Hubble Andromeda Treasury XI: The Spatially-Resolved Recent Star Formation History of M31

We measure the recent star formation history (SFH) across M31 using optical images taken with the \texit{Hubble Space Telescope} as part of the Panchromatic Hubble Andromeda Treasury (PHAT). We fit the color-magnitude diagrams in ~9000 regions that are ~100 pc $\times$ 100 pc in projected size, covering a 0.5 square degree area (~380 kpc$^2$, deprojected) in the NE quadrant of M31. We show that the SFHs vary significantly on these small spatial scales but that there are also coherent galaxy-wide fluctuations in the SFH back to ~500 Myr, most notably in M31's 10-kpc star-forming ring. We find that the 10-kpc ring is at least 400 Myr old, showing ongoing star formation over the past ~500 Myr. This indicates the presence of molecular gas in the ring over at least 2 dynamical times at this radius. We also find that the ring's position is constant throughout this time, and is stationary at the level of 1 km/s, although there is evidence for broadening of the ring due to diffusion of stars into the disk. Based on existing models of M31's ring features, the lack of evolution in the ring's position makes a purely collisional ring origin highly unlikely. We find that the global SFR has been fairly constant over the last ~500 Myr, though it does show a small increase at 50 Myr that is 1.3 times the average SFR over the past 100 Myr. During the last ~500 Myr, ~60% of all SF occurs in the 10-kpc ring. Finally, we find that in the past 100 Myr, the average SFR over the PHAT survey area is $0.28\pm0.03$ M$_\odot$ yr$^{-1}$ with an average deprojected intensity of $7.3 \times 10^{-4}$ M$_\odot$ yr$^{-1}$ kpc$^{-2}$, which yields a total SFR of ~0.7 M$_\odot$ yr$^{-1}$ when extrapolated to the entire area of M31's disk. This SFR is consistent with measurements from broadband estimates. [abridged]

preprint2015arXiv

The Relationship between CO Emission and Visual Extinction Traced by Dust Emission in the Magellanic Clouds

To test the theoretical understanding that finding bright CO emission depends primarily on dust shielding, we investigate the relationship between CO emission ($I_{\rm CO}$) and the amount of dust (estimated from IR emission and expressed as "$A_V$") across the Large Magellanic Cloud, the Small Magellanic Cloud, and the Milky Way. We show that at our common resolution of 10 pc scales, $I_{\rm CO}$ given a fixed line-of-sight $A_V$ is similar across all three systems despite the difference in metallicity. We find some evidence for a secondary dependence of $I_{\rm CO}$ on radiation field; in the LMC, $I_{\rm CO}$ at a given $A_V$ is smaller in regions of high $T_{\rm dust}$, perhaps because of an increased photodissociating radiation field. We suggest a simple but useful picture in which the CO-to-H$_2$ conversion factor (\xco) depends on two separable factors: (1) the distribution of gas column densities, which maps to an extinction distribution via a dust-to-gas ratio; and (2) the dependence of $I_{\rm CO}$ on $A_V$. Assuming that the probability distribution function (PDF) of local Milky Way clouds is universal, this approach predicts a dependence of \xco\ on $Z$ between $Z^{-1}$ and $Z^{-2}$ above about a third solar metallicity. Below this metallicity, CO emerges from only the high column density parts of the cloud and so depends very sensitively on the adopted PDF and the H$_2$/{\sc Hi} prescription. The PDF of low metallicity clouds is thus of considerable interest and the uncertainty associated with even an ideal prescription for \xco\ at very low metallicity will be large.

preprint2015arXiv

Variations in the Star Formation Efficiency of the Dense Molecular Gas across the Disks of Star-Forming Galaxies

We present a new survey of HCN(1-0) emission, a tracer of dense molecular gas, focused on the little-explored regime of normal star-forming galaxy disks. Combining HCN, CO, and infrared (IR) emission, we investigate the role of dense gas in Star Formation (SF), finding systematic variations in both the apparent dense gas fraction and the apparent SF efficiency (SFE) of dense gas. The latter may be unexpected, given the popularity of gas density threshold models to explain SF scaling relations. We used the IRAM 30-m telescope to observe HCN(1-0) across 29 nearby disk galaxies whose CO(2-1) emission has previously been mapped by the HERACLES survey. Because our observations span a range of galactocentric radii, we are able to investigate the properties of the dense gas as a function of local conditions. We focus on how the IR/CO, HCN/CO, and IR/HCN ratios (observational cognates of the SFE, dense gas fraction, and dense gas SFE) depend on the stellar surface density and the molecular/atomic ratio. The HCN/CO ratio correlates tightly with these two parameters across a range of 2.1 dex and increases in the high surface density parts of galaxies. Simultaneously, the IR/HCN ratio decreases systematically with these same parameters and is ~6-8 times lower near galaxy centers than in the outer regions. For fixed line-mass conversion factors, these results are incompatible with a simple model in which SF depends only on the gas mass above some density threshold. Only a specific set of environment-dependent conversion factors can render our observations compatible with such a model. Whole cloud models, such as the theory of turbulence regulated SF, do a better job of matching our data. We explore one such model in which variations in the Mach number and in the mean density would respectively drive the trends within galaxy disks and the differences between disk and merging galaxies (abridged).

preprint2014arXiv

A Test of Star Formation Laws in Disk Galaxies. II. Dependence on dynamical properties

We use observed radial profiles of mass surface densities of total, $Σ_g$, & molecular, $Σ_{\rm H2}$, gas, rotation velocity & star formation rate (SFR) surface density, $Σ_{\rm sfr}$, of the molecular-rich ($Σ_{\rm H2}\geΣ_{\rm HI}/2$) regions of 16 nearby disk galaxies to test several star formation laws: a Kennicutt-Schmidt law, $Σ_{\rm sfr}=A_gΣ_{g,2}^{1.5}$; a Constant Molecular law, $Σ_{\rm sfr}=A_{\rm H2}Σ_{\rm H2,2}$; the turbulence-regulated laws of Krumholz & McKee (KM05) and Krumholz et al. (KMT09), a Gas-$Ω$ law, $Σ_{\rm sfr}=B_ΩΣ_gΩ$; and a shear-driven GMC Collision law, $Σ_{\rm sfr}=B_{\rm CC}Σ_gΩ(1-0.7β)$, where $β\equiv d {\rm ln} v_{\rm circ}/d {\rm ln} r$. If allowed one free normalization parameter for each galaxy, these laws predict the SFR with rms errors of factors of 1.4 - 1.8. If a single normalization parameter is used by each law for the entire galaxy sample, then rms errors range from factors of 1.5 - 2.1. Although the Constant Molecular law gives the smallest errors, the improvement over KMT, Kennicutt-Schmidt & GMC Collision laws is not especially significant, particularly given the different observational inputs that the laws utilize and the scope of included physics, which ranges from empirical relations to detailed treatment of interstellar medium processes. We next search for variation of star formation law parameters with local & global galactic dynamical properties of disk shear rate (related to $β$), rotation speed & presence of a bar. We demonstrate with high significance that higher shear rates enhance star formation efficiency per local orbital time. Such a trend is expected if GMC collisions play an important role in star formation, while an opposite trend would be expected if development of disk gravitational instabilities is the controlling physics.

preprint2014arXiv

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

We use ALMA to derive the mass, length, and time scales associated with the nuclear starburst in NGC 253. This region forms ~2 M_sun/yr of stars and resembles other starbursts in scaling relations, with star formation consuming the gas reservoir 10 times faster than in galaxy disks. We present observations of CO, the high effective density transitions HCN(1-0), HCO+(1-0), CS(2-1), and their isotopologues. We identify ten clouds that appear as peaks in line emission and enhancements in the HCN-to-CO ratio. These clouds are massive (~10^7 M_sun) structures with sizes (~30 pc) similar to GMCs in other systems. Compared to disk galaxy GMCs, they show high line widths (~20-40 km/s) given their size, with implied Mach numbers ~90. The clouds also show high surface (~6,000 M_sun/pc^2) and volume densities (n_H2~2,000 cm^-3). Given these, self-gravity can explain the line widths. This short free fall time (~0.7 Myr) helps explain the more efficient star formation in NGC 253. We also consider the starburst region as a whole. The geometry is confused by the high inclination, but simple models support a non-axisymmetric, bar-like geometry with a compact, clumpy region of high gas density embedded in an extended CO distribution. Even for the whole region, the surface density still exceeds that of a disk galaxy GMC. The orbital time (~10 Myr), disk free fall time (<~ 3 Myr), and disk crossing time (<~ 3 Myr) are each much shorter than in a normal spiral galaxy disk. Some but not all aspects of the structure correspond to predictions from assuming vertical dynamical equilibrium or a marginally stable rotating disk. Finally, the CO-to-H2 conversion factor implied by our cloud calculations is approximately Galactic, contrasting with results showing a low value for the whole starburst region. The contrast provides resolved support for the idea of mixed molecular ISM phases in starburst galaxies.

preprint2014arXiv

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

We present new Karl G. Jansky Very Large Array radio continuum images of the nuclei of Arp 220, the nearest ultra-luminous infrared galaxy. These images have both the angular resolution to study detailed morphologies of the two nuclei that power the system and sensitivity to a wide range of spatial scales. At 33 GHz, and with a resolution of 0".081 x 0".063 (29.9 x 23.3 pc), we resolve the emission surrounding both nuclei and conclude that is mostly synchrotron in nature. The spatial distributions of radio emission in both nuclei are well described by exponential profiles. These have deconvolved half-light radii of 51 and 35 pc for the eastern and western nuclei, and they match the number density profile of radio supernovae observed with very long baseline interferometry. This similarity might be due to the fast cooling of cosmic rays electrons caused by the presence of a strong (~ mG) magnetic field in this system. We estimate high luminosity surface densities of $\mathrm{Σ_{IR} \sim 4.2^{+1.6}_{-0.7} \times 10^{13}}$ (east) and $\mathrm{\sim 9.7^{+3.7}_{-2.4} \times 10^{13}~(west)~L_{\odot}~kpc^{-2}}$, and star formation rate surface densities of $\mathrm{Σ_{SFR} \sim 10^{3.7\pm0.1}}$ (east) and $\mathrm{\sim 10^{4.1\pm0.1}~(west)~M_{\odot}~yr^{-1}~kpc^{-2}}$. These values, especially for the western nucleus are, to our knowledge, the highest luminosity and star formation rate surface densities measured for any star-forming system. Despite these high values, the nuclei lie below the dusty Eddington limit in which radiation pressure is balanced only by self-gravity. The small measured sizes also imply that the nuclei of Arp 220 are only transparent in the frequency range ~ 5 to 350 GHz. Our results offer no clear evidence that an active galactic nucleus dominates the emission from either nucleus at 33 GHz.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51

Using data from the PdBI Arcsecond Whirlpool Survey (PAWS), we have generated the largest extragalactic Giant Molecular Cloud (GMC) catalog to date, containing 1,507 individual objects. GMCs in the inner M51 disk account for only 54% of the total 12CO(1-0) luminosity of the survey, but on average they exhibit physical properties similar to Galactic GMCs. We do not find a strong correlation between the GMC size and velocity dispersion, and a simple virial analysis suggests that 30% of GMCs in M51 are unbound. We have analyzed the GMC properties within seven dynamically-motivated galactic environments, finding that GMCs in the spiral arms and in the central region are brighter and have higher velocity dispersions than inter-arm clouds. Globally, the GMC mass distribution does not follow a simple power law shape. Instead, we find that the shape of the mass distribution varies with galactic environment: the distribution is steeper in inter-arm region than in the spiral arms, and exhibits a sharp truncation at high masses for the nuclear bar region. We propose that the observed environmental variations in the GMC properties and mass distributions are a consequence of the combined action of large-scale dynamical processes and feedback from high mass star formation. We describe some challenges of using existing GMC identification techniques for decomposing the 12CO(1-0) emission in molecule-rich environments, such as M51's inner disk.

preprint2014arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51

The kinematic complexity and the favorable position of M51 on the sky make this galaxy an ideal target to test different theories of spiral arm dynamics. Taking advantage of the new high resolution PdBI Arcsecond Whirlpool Survey (PAWS) data, we undertake a detailed kinematic study of M51 to characterize and quantify the origin and nature of the non-circular motions. Using a tilted-ring analysis supported by several other archival datasets we update the estimation of M51's position angle (PA=(173 +/- 3) deg) and inclination (i=(22 +/- 5) deg). Harmonic decomposition of the high resolution (40 pc) CO velocity field shows the first kinematic evidence of an m=3 wave in the inner disk of M51 with a corotation at R(CR,m=3)=1.1 +/- 0.1 kpc and a pattern speed of Omega_p(m=3) = 140 km/(s kpc). This mode seems to be excited by the nuclear bar, while the beat frequencies generated by the coupling between the m=3 mode and the main spiral structure confirm its density-wave nature. We observe also a signature of an m=1 mode that is likely responsible for the lopsidedness of M51 at small and large radii. We provide a simple method to estimate the radial variation of the amplitude of the spiral perturbation (Vsp) attributed to the different modes. The main spiral arm structure has <Vsp>=50-70 km/s, while the streaming velocity associated with the m=1 and m=3 modes is, in general, 2 times lower. Our joint analysis of HI and CO velocity fields at low and high spatial resolution reveals that the atomic and molecular gas phases respond differently to the spiral perturbation due to their different vertical distribution and emission morphology.

preprint2014arXiv

The Radio Continuum-Star Formation Rate Relation in WSRT SINGS Galaxies

We study the spatially resolved Radio Continuum-Star Formation Rate (RC-SFR) relation using state-of-the-art star-formation (SF) tracers in a sample of 17 THINGS galaxies. We use hybrid Sigma_SFR maps (GALEX FUV plus Spitzer 24 mu), RC maps at 22/18 cm from the WSRT SINGS survey, and H-alpha maps to correct for thermal RC emission. We compare azimuthally averaged radial profiles of the RC and FUV/MIR-based Sigma_SFR maps and study pixel-by-pixel correlations at fixed linear scales of 1.2 and 0.7 kpc. The ratio of the integrated SFRs from the RC emission to that of the FUV/MIR-based SF tracers is R_int = 0.78 +/- 0.38, consistent with Condon's relation. We find a tight correlation between the radial profiles of the radio and FUV/MIR-based Sigma_SFR for the entire extent of the disk. The ratio R of the azimuthally averaged radio to FUV/MIR-based Sigma_SFR agrees with the integrated ratio with only small quasi-random fluctuations as function of radius. Pixel-by-pixel plots show a tight correlation in log-log diagrams of radio to FUV/MIR-based Sigma_SFR, with a typical standard deviation of a factor of two. Averaged over our sample we find (Sigma_SFR)_RC ~ (Sigma_SFR)_hyb^{0.63+/-0.25} implying that data points with high Sigma_SFR are relatively radio dim, whereas the reverse is true for low Sigma_SFR. We interpret this as a result of spectral ageing of CRe, which is supported by the radio spectral index: data points dominated by young CRe are relatively radio dim, those dominated by old CRe are relatively radio bright. The ratio of radio to FUV/MIR-based integrated SFR is independent of global galaxy parameters, suggesting that we can use RC emission as a universal SF tracer for galaxies, if we restrict ourselves to global or azimuthally averaged measurements. A magnetic field-SFR relation, B ~ SFR_hyb^{0.30+/-0.02}, holding both globally and locally, can explain our results. (abridged)

preprint2013arXiv

A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud

We compare the properties of giant molecular clouds (GMCs) in M51 identified by the Plateau de Bure Interferometer Whirlpool Arcsecond Survey (PAWS) with GMCs identified in wide-field, high resolution surveys of CO emission in M33 and the Large Magellanic Cloud (LMC). We find that GMCs in M51 are larger, brighter and have higher velocity dispersions relative to their size than equivalent structures in M33 and the LMC. These differences imply that there are genuine variations in the average mass surface density of the different GMC populations. To explain this, we propose that the pressure in the interstellar medium surrounding the GMCs plays a role in regulating their density and velocity dispersion. We find no evidence for a correlation between size and linewidth in any of M51, M33 or the LMC when the CO emission is decomposed into GMCs, although moderately robust correlations are apparent when regions of contiguous CO emission (with no size limitation) are used. Our work demonstrates that observational bias remains an important obstacle to the identification and study of extragalactic GMC populations using CO emission, especially in molecule-rich galactic environments.

preprint2013arXiv

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING). III. The Dependence of Atomic and Molecular Gas Surface Densities on Galaxy Properties

We investigate the correlation between CO and HI emission in 18 nearby galaxies from the CARMA Survey Toward IR-Bright Nearby Galaxies (STING) at sub-kpc and kpc scales. Our sample, spanning a wide range in stellar mass and metallicity, reveals evidence for a metallicity dependence of the HI column density measured in regions exhibiting CO emission. Such a dependence is predicted by the equilibrium model of McKee & Krumholz, which balances H_2 formation and dissociation. The observed HI column density is often smaller than predicted by the model, an effect we attribute to unresolved clumping, although values close to the model prediction are also seen. We do not observe HI column densities much larger than predicted, as might be expected were there a diffuse HI component that did not contribute to H_2 shielding. We also find that the H_2 column density inferred from CO correlates strongly with the stellar surface density, suggesting that the local supply of molecular gas is tightly regulated by the stellar disk.

preprint2013arXiv

Clumping and the Interpretation of kpc-Scale Maps of the Interstellar Medium: Smooth HI and Clumpy, Variable H2 Surface Density

Many recent models consider the structure of individual interstellar medium (ISM) clouds as a way to explain observations of large parts of galaxies. To compare such models to observations, one must understand how to translate between surface densities observed averaging over large (~kpc) scales and surface densities on the scale of individual clouds (~pc scale), which are treated by models. We define a "clumping factor" that captures this translation as the ratio of the mass-weighted surface density, which is often the quantity of physical interest, to the area-weighted surface density, which is observed. We use high spatial resolution (sub-kpc) maps of CO and HI emission from nearby galaxies to measure the clumping factor of both atomic and molecular gas. The molecular and atomic ISM exhibit dramatically different degrees of clumping. As a result, the ratio H2/HI measured at ~kpc resolution cannot be trivially interpreted as a cloud-scale ratio of surface densities. HI emission appears very smooth, with a clumping factor of only ~1.3. Based on the scarce and heterogeneous high resolution data available, CO emission is far more clumped with a widely variable clumping factor, median ~7 for our heterogeneous data. Our measurements do not provide evidence for a universal mass-weighted surface density of molecular gas, but also cannot conclusively rule out such a scenario. We suggest that a more sophisticated treatment of molecular ISM structure, one informed by high spatial resolution CO maps, is needed to link cloud-scale models to kpc-scale observations of galaxies.

preprint2013arXiv

Dust May Be More Rare Than Expected in Metal Poor Galaxies

'Normal' galaxies observed at z>6, when the Universe was <1 billion years old, thus far show no evidence of the cold dust that accompanies star formation in the local Universe, where the dust-to-gas mass ratio is 1%. A prototypical example is 'Himiko' (z=6.6), which a mere 840 Myr after the Big Bang is forming stars at a rate of 30-100 Msun/yr, yielding a mass assembly time M^{star}/SFR 150x10^6 yr. Himiko is estimated to have a low fraction (2-3% of the Solar value) of elements heavier than helium (metallicity), and although its gas mass cannot be asserted at this time its dust-to-stellar mass ratio is constrained to be <0.05%. The local galaxy I Zw 18, with a metallicity 4% solar and forming stars less rapidly than Himiko but still vigorously for its mass (M^{star}/SFR 1.6x10^9 yr), is also very dust deficient and perhaps one of the best analogues of primitive galaxies accessible to detailed study. Here we report observations of dust emission from I Zw 18 from which we determine its dust mass to be 450-1800 Msun, yielding a dust-to-stellar mass ratio \approx 10^{-6}-10^{-5} and a dust-to-gas mass ratio 3.2-13x10^{-6}. If I Zw 18 is a reasonable analog of Himiko, then Himiko's dust mass is \approx 50,000 Msun, a factor of 100 below the current upper limit. These numbers are considerably uncertain, but if most high-z galaxies are more like Himiko than like the quasar host SDSS J114816.64+525150.3, then the prospects for detecting the gas and dust in them are much poorer than hitherto anticipated.

preprint2013arXiv

Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure

We use the high spatial and spectral resolution of the PAWS CO(1-0) survey of the inner 9 kpc of the iconic spiral galaxy M51 to examine the effect of gas streaming motions on the star-forming properties of individual GMCs. We compare our view of gas flows in M51 -- which arise due to departures from axi-symmetry in the gravitational potential (i.e. the nuclear bar and spiral arms) -- with the global pattern of star formation as traced by Halpha and 24μm emission. We find that the dynamical environment of GMCs strongly affects their ability to form stars, in the sense that GMCs situated in regions with large streaming motions can be stabilized, while similarly massive GMCs in regions without streaming go on to efficiently form stars. We argue that this is the result of reduced surface pressure felt by clouds embedded in an ambient medium undergoing large streaming motions, which prevents collapse. Indeed, the variation in gas depletion time expected based on the observed streaming motions throughout the disk of M51 quantitatively agrees with the variation in observed gas depletion time scale. The example of M51 shows that streaming motions, triggered by gravitational instabilities in the form of bars and spiral arms, can alter the star formation law; this can explain the variation in gas depletion time among galaxies with different masses and morphologies. In particular, we can explain the long gas depletion times in spiral galaxies compared to dwarf galaxies and starbursts. We suggest that adding a dynamical pressure term to the canonical free-fall time produces a single star formation law that can be applied to all star-forming regions and galaxies, across cosmic time.

preprint2013arXiv

Molecular Gas and Star Formation in Nearby Disk Galaxies

We compare molecular gas traced by 12CO(2-1) maps from the HERACLES survey, with tracers of the recent star formation rate (SFR) across 30 nearby disk galaxies. We demonstrate a first-order linear correspondence between Sig_mol and Sig_SFR but also find important second-order systematic variations in the apparent molecular gas depletion time, t_dep^mol = Sig_mol / Sig_SFR. At our 1 kpc common resolution, CO correlates closely with many tracers of the recent SFR. Weighting each line of sight equally and using a fixed, Milky Way alpha_CO, our data yield a molecular gas depletion time, t_dep^mol=Sig_mol/Sig_SFR ~ 2.2 Gyr with 0.3 dex scatter, in good agreement with literature data. We apply a forward-modeling approach to constrain the power-law index, N, that relates the SFR surface density and the molecular gas surface density and find N=1+/-0.15 for our full data set with some variation from galaxy to galaxy. However, we caution that a power law treatment oversimplifies the topic given that we observe correlations between t_dep^mol and other local and global quantities. The strongest of these are a decreased t_dep^mol in low-mass, low-metallicity galaxies and a correlation of the kpc-scale t_dep^mol with dust-to-gas ratio, D/G. These correlations can be explained by a CO-to-H2 conversion factor that depends on D/G in the theoretically expected way. This is not a unique interpretation, but external evidence of conversion factor variations makes it a conservative one. After applying a D/G-dependent alpha_CO, some weak correlations between t_dep^mol and local conditions persist. In particular, we observe lower t_dep^mol and enhanced CO excitation associated with some nuclear gas concentrations. These appear to reflect real enhancements in the SFR/H2 and t_dep appears multivalued at fixed Sig_mol, supporting the the idea of "disk" and "starburst" modes driven by environmental factors.

preprint2013arXiv

Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View

We analyse the distribution of CO brightness temperature and integrated intensity in M51 at ~40 pc resolution using new CO data from the Plateau de Bure Arcsecond Whirlpool Survey (PAWS). We present probability distribution functions (PDFs) of the CO emission within the PAWS field, which covers the inner 11 x 7 kpc of M51. We find variations in the shape of CO PDFs within different M51 environments, and between M51 and M33 and the Large Magellanic Cloud (LMC). Globally, the PDFs for the inner disk of M51 can be represented by narrow lognormal functions that cover 1 to 2 orders of magnitude in CO brightness and integrated intensity. The PDFs for M33 and the LMC are narrower and peak at lower CO intensities. However, the CO PDFs for different dynamical environments within the PAWS field depart from the shape of the global distribution. The PDFs for the interarm region are approximately lognormal, but in the spiral arms and central region of M51, they exhibit diverse shapes with a significant excess of bright CO emission. The observed environmental dependence of the shape of the CO PDFs is qualitatively consistent with changes that would be expected if molecular gas in the spiral arms has a larger range of average densities, gas temperatures and velocity fluctuations, though further work is required to disentangle the importance of large-scale dynamical effects versus star formation feedback in regulating these properties. We show that the shape of the CO PDFs for different M51 environments is only weakly related to global properties of the CO emission, but is strongly correlated with some properties of the local giant molecular cloud (GMC) and young stellar cluster populations. For galaxies with strong spiral structure such as M51, our results indicate that galactic-scale dynamical processes play a significant role in the formation and evolution of GMCs and stellar clusters.(abridged)

preprint2013arXiv

The CO-to-H2 Conversion Factor

CO line emission represents the most accessible and widely used tracer of the molecular interstellar medium. This renders the translation of observed CO intensity into total H2 gas mass critical to understand star formation and the interstellar medium in our Galaxy and beyond. We review the theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H2 "conversion factor," Xco, in different environments. In the Milky Way disk, we recommend a conversion factor Xco = 2x10^{20} cm^-2/(K km/s)^-1 with +/-30% uncertainty. Studies of other "normal galaxies" return similar values in Milky Way-like disks, but with greater scatter and systematic uncertainty. Departures from this Galactic conversion factor are both observed and expected. Dust-based determinations, theoretical arguments, and scaling relations all suggest that Xco increases with decreasing metallicity, turning up sharply below metallicity ~1/3-1/2 solar in a manner consistent with model predictions that identify shielding as a key parameter. Based on spectral line modeling and dust observations, Xco appears to drop in the central, bright regions of some but not all galaxies, often coincident with regions of bright CO emission and high stellar surface density. This lower Xco is also present in the overwhelmingly molecular interstellar medium of starburst galaxies, where several lines of evidence point to a lower CO-to-H2 conversion factor. At high redshift, direct evidence regarding the conversion factor remains scarce; we review what is known based on dynamical modeling and other arguments.

preprint2013arXiv

The Fueling Diagram: Linking Galaxy Molecular-to-Atomic Gas Ratios to Interactions and Accretion

To assess how external factors such as local interactions and fresh gas accretion influence the global ISM of galaxies, we analyze the relationship between recent enhancements of central star formation and total molecular-to-atomic (H2/HI) gas ratios, using a broad sample of field galaxies spanning early-to-late type morphologies, stellar masses of 10^(7.2-11.2) Msun, and diverse stages of evolution. We find that galaxies occupy several loci in a "fueling diagram" that plots H2/HI vs. mass-corrected blue-centeredness, a metric tracing the degree to which galaxies have bluer centers than the average galaxy at their stellar mass. Spiral galaxies show a positive correlation between H2/HI and mass-corrected blue-centeredness. When combined with previous results linking mass-corrected blue-centeredness to external perturbations, this correlation suggests a link between local galaxy interactions and molecular gas inflow/replenishment. Intriguingly, E/S0 galaxies show a more complex picture: some follow the same correlation, some are quenched, and a distinct population of blue-sequence E/S0 galaxies (with masses below key transitions in gas richness) defines a separate loop in the fueling diagram. This population appears to be composed of low-mass merger remnants currently in late- or post-starburst states, in which the burst first consumes the H2 while the galaxy center keeps getting bluer, then exhausts the H2, at which point the burst population reddens as it ages. Multiple lines of evidence suggest connected evolutionary sequences in the fueling diagram. In particular, tracking total gas-to-stellar mass ratios within the diagram provides evidence of fresh gas accretion onto low-mass E/S0s emerging from central starbursts. Drawing on a comprehensive literature search, we suggest that virtually all galaxies follow the same evolutionary patterns found in our broad sample.

preprint2013arXiv

The Green Bank Telescope Maps the Dense, Star-Forming Gas in the Nearby Starburst Galaxy M82

Observations of the Milky Way and nearby galaxies show that dense molecular gas correlates with recent star formation, suggesting that the formation of this gas phase may help regulate star formation. A key test of this idea requires wide-area, high-resolution maps of dense molecular gas in galaxies to explore how local physical conditions drive dense gas formation, but these observations have been limited because of the faintness of dense gas tracers like HCN and HCO+. Here we demonstrate the power of the Robert C. Byrd Green Bank Telescope -- the largest single-dish millimeter radio telescope -- for mapping dense gas in galaxies by presenting the most sensitive maps yet of HCN and HCO+ in the starburst galaxy M82. The HCN and HCO+ in the disk of this galaxy correlates with both recent star formation and more diffuse molecular gas and shows kinematics consistent with a rotating torus. The HCO+ emission extending to the north and south of the disk is coincident with the outflow previously identified in CO and traces the eastern edge of the hot outflowing gas. The central starburst region has a higher ratio of star formation to dense gas than the outer regions, pointing to the starburst as a key driver of this relationship. These results establish that the GBT can efficiently map the dense molecular gas at 90 GHz in nearby galaxies, a capability that will increase further with the 16 element feed array under construction.

preprint2013arXiv

The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

The PdBI (Plateau de Bure Interferometer) Arcsecond Whirlpool Survey (PAWS) has mapped the molecular gas in the central ~9kpc of M51 in its 12CO(1-0) line emission at cloud-scale resolution of ~40pc using both IRAM telescopes. We utilize this dataset to quantitatively characterize the relation of molecular gas (or CO emission) to other tracers of the interstellar medium (ISM), star formation and stellar populations of varying ages. Using 2-dimensional maps, a polar cross-correlation technique and pixel-by-pixel diagrams, we find: (a) that (as expected) the distribution of the molecular gas can be linked to different components of the gravitational potential, (b) evidence for a physical link between CO line emission and radio continuum that seems not to be caused by massive stars, but rather depend on the gas density, (c) a close spatial relation between the PAH and molecular gas emission, but no predictive power of PAH emission for the molecular gas mass,(d) that the I-H color map is an excellent predictor of the distribution (and to a lesser degree the brightness) of CO emission, and (e) that the impact of massive (UV-intense) young star-forming regions on the bulk of the molecular gas in central ~9kpc can not be significant due to a complex spatial relation between molecular gas and star-forming regions that ranges from co-spatial to spatially offset to absent. The last point, in particular, highlights the importance of galactic environment -- and thus the underlying gravitational potential -- for the distribution of molecular gas and star formation.

preprint2013arXiv

The Starburst-Driven Molecular Wind in NGC 253 and the Suppression of Star Formation

The under-abundance of very massive galaxies in the universe is frequently attributed to the effect of galactic winds. Although ionized galactic winds are readily observable most of the expelled mass is likely in cooler atomic and molecular phases. Expanding molecular shells observed in starburst systems such as NGC 253 and M 82 may facilitate the entrainment of molecular gas in the wind. While shell properties are well constrained, determining the amount of outflowing gas emerging from such shells and the connection between this gas and the ionized wind requires spatial resolution <100 pc coupled with sensitivity to a wide range of spatial scales, hitherto not available. Here we report observations of NGC 253, a nearby starburst galaxy (D~3.4 Mpc) known to possess a wind, which trace the cool molecular wind at 50 pc resolution. At this resolution the extraplanar molecular gas closely tracks the Hα filaments, and it appears connected to molecular expanding shells located in the starburst region. These observations allow us to directly measure the molecular outflow rate to be > 3 Msun/yr and likely ~9 Msun/yr. This implies a ratio of mass-outflow rate to star formation rate of at least η~1-3, establishing the importance of the starburst-driven wind in limiting the star formation activity and the final stellar content.

preprint2012arXiv

A High Resolution Study of the HI-H2 Transition across the Perseus Molecular Cloud

To investigate the fundamental principles of H2 formation in a giant molecular cloud (GMC), we derive the HI and H2 surface density (Sigma_HI and Sigma_H2) images of the Perseus molecular cloud on sub-pc scales (~0.4 pc). We use the far-infrared data from the Improved Reprocessing of the IRAS Survey and the V-band extinction image provided by the COMPLETE Survey to estimate the dust column density image of Perseus. In combination with the HI data from the Galactic Arecibo L-band Feed Array HI Survey and an estimate of the local dust-to-gas ratio, we then derive the Sigma_H2 distribution across Perseus. We find a relatively uniform Sigma_HI ~ 6-8 Msun pc^-2 for both dark and star-forming regions, suggesting a minimum HI surface density required to shield H2 against photodissociation. As a result, a remarkably tight and consistent relation is found between Sigma_H2/Sigma_HI and Sigma_HI+Sigma_H2. The transition between the HI- and H2-dominated regions occurs at N(HI)+2N(H2) ~ (8-14) x 10^20 cm^-2. Our findings are consistent with predictions for H2 formation in equilibrium, suggesting that turbulence may not be of primary importance for H2 formation. However, the importance of a warm neutral medium for H2 shielding, an internal radiation field, and the timescale of H2 formation still remain as open questions. We also compare H2 and CO distributions and estimate the fraction of "CO-dark" gas, f_DG ~ 0.3. While significant spatial variations of f_DG are found, we do not find a clear correlation with the mean V-band extinction.

preprint2012arXiv

Dust-to-Gas Ratio in the Extremely Metal Poor Galaxy I ZW 18

The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3σ$ limit). This upper limit is driven by the non-detection at 160 $μ$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.

preprint2012arXiv

Estimating the Star Formation Rate at 1 kpc Scales in Nearby Galaxies

Using combinations of Hα, ultraviolet (UV), and infrared (IR) emission, we estimate the star formation rate (SFR) surface density, Σ_SFR, at 1 kpc resolution for 30 disk galaxies that are targets of the IRAM HERACLES CO survey. We present a new physically-motivated IR spectral energy distribution-based approach to account for possible contributions to 24\mum emission not associated with recent star formation. Considering a variety of "reference" SFRs from the literature, we revisit the calibration of the 24\mum term in hybrid (UV+IR or Hα+IR) tracers. We show that the overall calibration of this term remains uncertain at the factor of two level because of the lack of wide-field, robust reference SFR estimates. Within this uncertainty, published calibrations represent a reasonable starting point for 1 kpc-wide areas of star-forming disk galaxies but we re-derive and refine the calibration of the IR term in these tracers to match our resolution and approach to 24\mum emission. We compare a large suite of Σ_SFR estimates and find that above Σ_SFR \sim 10^-3 M_\odot yr^-1 kpc^-2 the systematic differences among tracers are less than a factor of two across two orders of magnitude dynamic range. We caution that methodology and data both become serious issues below this level. We note from simple model considerations that focusing on a part of a galaxy dominated by a single stellar population the intrinsic uncertainty in Hαand FUV-based SFRs are \sim 0.3 and \sim 0.5 dex.

preprint2012arXiv

Low CO Luminosities in Dwarf Galaxies

[Abridged] We present maps of CO 2-1 emission covering the entire star-forming disks of 16 nearby dwarf galaxies observed by the IRAM HERACLES survey. The data have 13 arcsec angular resolution, ~250 pc at our average distance of 4 Mpc, and sample the galaxies by 10-1000 resolution elements. We apply stacking techniques to perform the first sensitive search for CO emission in dwarfs outside the Local Group ranging from single lines-of-sight, stacked over IR-bright regions of embedded star formation, and stacked over the entire galaxy. We detect 5 dwarfs in CO with total luminosities of L_CO = 3-28 1e6 Kkmspc2. The other 11 dwarfs remain undetected in CO even in the stacked data and have L_CO < 0.4-8 1e6 Kkmspc2. We combine our sample of dwarfs with a large literature sample of spirals to study scaling relations of L_CO with M_B and metallicity. We find that dwarfs with metallicities of Z ~ 1/2-1/10 Z_sun have L_CO about 1e2-1e4x smaller than spirals and that their L_CO per unit L_B is 10-100x smaller. A comparison with tracers of star formation (FUV and 24 micron) shows that L_CO per unit SFR is 10-100x smaller in dwarfs. One possible interpretation is that dwarfs form stars much more efficiently, however we argue that the low L_CO/SFR ratio is due to significant changes of the CO-to-H2 conversion factor, alpha_CO, in low metallicity environments. Assuming a constant H2 depletion time of 1.8 Gyr (as found for nearby spirals) implies alpha_CO values for dwarfs with Z ~ 1/2-1/10 Z_sun that are more than 10x higher than those found in solar metallicity spirals. This significant increase of alpha_CO at low metallicity is consistent with previous studies, in particular those which model dust emission to constrain H2 masses. Even though it is difficult to parameterize the metallicity dependence of alpha_CO, our results suggest that CO is increasingly difficult to detect at lower metallicities.

preprint2012arXiv

The VIRUS-P Exploration of Nearby Galaxies (VENGA): The Xco Gradient in NGC 628

We measure the radial profile of the 12CO(1-0) to H_2 conversion factor (Xco) in NGC 628. The Hαemission from the VENGA integral field spectroscopy is used to map the star formation rate surface density (Σ_{SFR}). We estimate the molecular gas surface density (Σ_{H2}) from Σ_{SFR} by inverting the molecular star formation law (SFL), and compare it to the CO intensity to measure Xco. We study the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers (Hαvs. far-UV plus 24μm), and CO maps from different telescopes (single-dish and interferometers). The observed Xco profile is robust against these systematics, drops by a factor of 2 from R~7 kpc to the center of the galaxy, and is well fit by a gradient Δlog(Xco)=0.06\pm0.02 dex kpc^-1. We study how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical models show that the gradient in Xco can be explained by a combination of decreasing metallicity, and decreasing Σ_{H2} with radius. Photoelectric heating from the local UV radiation field appears to contribute to the decrease of Xco in higher density regions. Our results show that galactic environment plays an important role at setting the physical conditions in star forming regions, in particular the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We caution against adopting a single Xco value when large changes in gas surface density or metallicity are present.

preprint2011arXiv

A Molecular Star Formation Law in the Atomic Gas Dominated Regime in Nearby Galaxies

We use the IRAM HERACLES survey to study CO emission from 33 nearby spiral galaxies down to very low intensities. Using atomic hydrogen (HI) data, mostly from THINGS, we predict the local mean CO velocity from the mean HI velocity. By renormalizing the CO velocity axis so that zero corresponds to the local mean HI velocity we are able to stack spectra coherently over large regions as function of radius. This enables us to measure CO intensities with high significance as low as Ico = 0.3 K km/s (H2_SD = 1 Msun/pc2), an improvement of about one order of magnitude over previous studies. We detect CO out to radii Rgal = R25 and find the CO radial profile to follow a uniform exponential decline with scale length of 0.2 R25. Comparing our sensitive CO profiles to matched profiles of HI, Halpha, FUV, and IR emission at 24um and 70um, we observe a tight, roughly linear relation between CO and IR intensity that does not show any notable break between regions that are dominated by molecular (H2) gas (H2_SD > HI_SD) and those dominated by atomic gas (H2_SD < HI_SD). We use combinations of FUV+24um and Halpha+24um to estimate the recent star formation rate (SFR) surface density, SFR_SD, and find approximately linear relations between SFR_SD and H2_SD. We interpret this as evidence for stars forming in molecular gas with little dependence on the local total gas surface density. While galaxies display small internal variations in the SFR-to-H2 ratio, we do observe systematic galaxy-to-galaxy variations. These galaxy-to-galaxy variations dominate the scatter in relations between CO and SFR tracers measured at large scales. The variations have the sense that less massive galaxies exhibit larger ratios of SFR-to-CO than massive galaxies. Unlike the SFR-to-CO ratio, the balance between HI and H2 depends strongly on the total gas surface density and radius. It must also depend on additional parameters.

preprint2011arXiv

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING) II: Molecular Gas Star Formation Law and Depletion Time Across the Blue Sequence

We present an analysis of the relationship between molecular gas and current star formation rate surface density at sub-kpc and kpc scales in a sample of 14 nearby star-forming galaxies. Measuring the relationship in the bright, high molecular gas surface density ($\Shtwo\gtrsim$20 \msunpc) regions of the disks to minimize the contribution from diffuse extended emission, we find an approximately linear relation between molecular gas and star formation rate surface density, $\nmol\sim0.96\pm0.16$, with a molecular gas depletion time $\tdep\sim2.30\pm1.32$ Gyr. We show that, in the molecular regions of our galaxies there are no clear correlations between \tdep\ and the free-fall and effective Jeans dynamical times throughout the sample. We do not find strong trends in the power-law index of the spatially resolved molecular gas star formation law or the molecular gas depletion time across the range of galactic stellar masses sampled (\mstar $\sim$$10^{9.7}-10^{11.5}$ \msun). There is a trend, however, in global measurements that is particularly marked for low mass galaxies. We suggest this trend is probably due to the low surface brightness CO, and it is likely associated with changes in CO-to-H2 conversion factor.

preprint2011arXiv

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING): Molecular Gas Star Formation Law in NGC4254

This study explores the effects of different assumptions and systematics on the determination of the local, spatially resolved star formation law. Using four star formation rate (SFR) tracers (Hαwith azimuthally averaged extinction correction, mid-infrared 24 micron, combined Hαand mid-infrared 24 micron, and combined far-ultraviolet and mid-infrared 24 micron), several fitting procedures, and different sampling strategies we probe the relation between SFR and molecular gas at various spatial resolutions and surface densities within the central 6.5 kpc in the disk of NGC4254. We find that in the high surface brightness regions of NGC4254 the form of the molecular gas star formation law is robustly determined and approximately linear and independent of the assumed fraction of diffuse emission and the SFR tracer employed. When the low surface brightness regions are included, the slope of the star formation law depends primarily on the assumed fraction of diffuse emission. In such case, results range from linear when the fraction of diffuse emission in the SFR tracer is ~30% or less (or when diffuse emission is removed in both the star formation and the molecular gas tracer), to super-linear when the diffuse fraction is ~50% and above. We find that the tightness of the correlation between gas and star formation varies with the choice of star formation tracer. The 24 micron SFR tracer by itself shows the tightest correlation with the molecular gas surface density, whereas the Hαcorrected for extinction using an azimuthally-averaged correction shows the highest dispersion. We find that for R<0.5R_25 the local star formation efficiency is constant and similar to that observed in other large spirals, with a molecular gas depletion time ~2 Gyr.

preprint2011arXiv

Complex Radio Spectral Energy Distributions in Luminous and Ultraluminous Infrared Galaxies

We use the Expanded Very Large Array to image radio continuum emission from local luminous and ultraluminous infrared galaxies (LIRGs and ULIRGs) in 1 GHz windows centered at 4.7, 7.2, 29, and 36 GHz. This allows us to probe the integrated radio spectral energy distribution (SED) of the most energetic galaxies in the local universe. The 4-8 GHz flux densities agree well with previous measurements. They yield spectral indices α\approx -0.67 (where F_ν\propto ν^α) with \pm 0.15 (1σ) scatter, typical of nonthermal (synchrotron) emission from star-forming galaxies. The contrast of our 4-8 GHz data with literature 1.5 and 8.4 GHz flux densities gives further evidence for curvature of the radio SED of U/LIRGs. The SED appears flatter near \sim 1 GHz than near \sim 6 GHz, suggesting significant optical depth effects at the lower frequencies. The high frequency (28-37 GHz) flux densities are low compared to extrapolations from the 4-8 GHz data. We confirm and extend to higher frequency a previously observed deficit of high frequency radio emission for luminous starburst galaxies.

preprint2011arXiv

Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview

The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.

preprint2011arXiv

The CO-to-H2 Conversion Factor From Infrared Dust Emission Across the Local Group

We estimate the conversion factor relating CO emission to H2 mass, alpha_CO, in five Local Group galaxies that span approximately an order of magnitude in metallicity - M31, M 33, the Large Magellanic Cloud (LMC), NGC 6822, and the Small Magellanic Cloud (SMC). We model the dust mass along the line of sight from infrared (IR) emission and then solve for the alpha_CO that best allows a single gas-to-dust ratio (delta_GDR) to describe each system. This approach remains sensitive to CO-dark envelopes of H2 surrounding molecular clouds. In M 31, M 33, and the LMC we find alpha_CO \approx 3-9 M_sun pc^-2 (K km s^-1)^-1, consistent with the Milky Way value within the uncertainties. The two lowest metallicity galaxies in our sample, NGC 6822 and the SMC (12 + log(O/H) \approx 8.2 and 8.0), exhibit a much higher alpha_CO. Our best estimates are α_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and α_SMC \approx 70 M_sun/pc^-2 (K km s-1)-1. These results are consistent with the conversion factor becoming CO a strong function of metallicity around 12 + log(O/H) \sim 8.4 - 8.2. We favor an interpretation where decreased dust-shielding leads to the dominance of CO-free envelopes around molecular clouds below this metallicity.

preprint2011arXiv

The Spitzer Spectroscopic Survey of the Small Magellanic Cloud (S4MC): Probing the Physical State of Polycyclic Aromatic Hydrocarbons in a Low-Metallicity Environment

We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.

preprint2011arXiv

The State of the Gas and the Relation Between Gas and Star Formation at Low Metallicity: the Small Magellanic Cloud

We compare atomic gas, molecular gas, and the recent star formation rate (SFR) inferred from H-alpha in the Small Magellanic Cloud (SMC). By using infrared dust emission and local dust-to-gas ratios, we construct a map of molecular gas that is independent of CO emission. This allows us to disentangle conversion factor effects from the impact of metallicity on the formation and star formation efficiency of molecular gas. On scales of 200 pc to 1 kpc we find a characteristic molecular gas depletion time of ~1.6 Gyr, similar to that observed in the molecule-rich parts of large spiral galaxies on similar spatial scales. This depletion time shortens on much larger scales to ~0.6 Gyr because of the presence of a diffuse H-alpha component, and lengthens on much smaller scales to ~7.5 Gyr because the H-alpha and H2 distributions differ in detail. We estimate the systematic uncertainties in our measurement to be a factor of 2-3. We suggest that the impact of metallicity on the physics of star formation in molecular gas has at most this magnitude. The relation between SFR and neutral (H2+HI) gas surface density is steep, with a power-law index ~2.2+/-0.1, similar to that observed in the outer disks of large spiral galaxies. At a fixed total gas surface density the SMC has a 5-10 times lower molecular gas fraction (and star formation rate) than large spiral galaxies. We explore the ability of the recent models by Krumholz et al. (2009) and Ostriker et al. (2010) to reproduce our observations. We find that to explain our data at all spatial scales requires a low fraction of cold, gravitationally-bound gas in the SMC. We explore a combined model that incorporates both large scale thermal and dynamical equilibrium and cloud-scale photodissociation region structure and find that it reproduces our data well, as well as predicting a fraction of cold atomic gas very similar to that observed in the SMC.

preprint2011arXiv

What Phase of the Interstellar Medium Correlates with the Star Formation Rate?

Nearby spiral galaxies show an extremely tight correlation between tracers of molecular hydrogen (H_2) in the interstellar medium (ISM) and tracers of recent star formation, but it is unclear whether this correlation is fundamental or accidental. In the galaxies that have been surveyed to date, H_2 resides predominantly in gravitationally bound clouds cooled by carbon monoxide (CO) molecules, but in galaxies of low metal content the correlations between bound clouds, CO, and H_2 break down, and it is unclear if the star formation rate will then correlate with H_2 or with some other quantity. Here we show that star formation will continue to follow H_2 independent of metallicity. This is not because H_2 is directly important for cooling, but instead because the transition from predominantly atomic hydrogen (HI) to H_2 occurs under the same conditions as a dramatic drop in gas temperature and Bonnor-Ebert mass that destabilizes clouds and initiates collapse. We use this model to compute how star formation rate will correlate with total gas mass, with mass of gas where the hydrogen is H_2, and with mass of gas where the carbon is CO in galaxies of varying metallicity, and show that preliminary observations match the trend we predict.

preprint2010arXiv

Regulation of Star Formation Rates in Multiphase Galactic Disks: a Thermal/Dynamical Equilibrium Model

We develop a model for regulation of galactic star formation rates Sigma_SFR in disk galaxies, in which ISM heating by stellar UV plays a key role. By requiring simultaneous thermal and (vertical) dynamical equilibrium in the diffuse gas, and star formation at a rate proportional to the mass of the self-gravitating component, we obtain a prediction for Sigma_SFR as a function of the total gaseous surface density Sigma and the density of stars + dark matter, rho_sd. The physical basis of this relationship is that thermal pressure in the diffuse ISM, which is proportional to the UV heating rate and therefore to Sigma_SFR, must adjust to match the midplane pressure set by the vertical gravitational field. Our model applies to regions where Sigma < 100 Msun/pc^2. In low-Sigma_SFR (outer-galaxy) regions where diffuse gas dominates, the theory predicts Sigma_SFR \propto Sigma (rho_sd)^1/2. The decrease of thermal equilibrium pressure when Sigma_SFR is low implies, consistent with observations, that star formation can extend (with declining efficiency) to large radii in galaxies, rather than having a sharp cutoff. The main parameters entering our model are the ratio of thermal pressure to total pressure in the diffuse ISM, the fraction of diffuse gas that is in the warm phase, and the star formation timescale in self-gravitating clouds; all of these are (in principle) direct observables. At low surface density, our model depends on the ratio of the mean midplane FUV intensity (or thermal pressure in the diffuse gas) to the star formation rate, which we set based on Solar neighborhood values. We compare our results to recent observations, showing good agreement overall for azimuthally-averaged data in a set of spiral galaxies. For the large flocculent spiral galaxies NGC 7331 and NGC 5055, the correspondence between theory and observation is remarkably close.

preprint2010arXiv

The Scale Dependence of the Molecular Gas Depletion Time in M33

We study the Local Group spiral galaxy M33 to investigate how the observed scaling between the (kpc-averaged) surface density of molecular gas (Σ_H2) and recent star formation rate (Σ_SFR) relates to individual star-forming regions. To do this, we measure the ratio of CO emission to extinction-corrected Halpha emission in apertures of varying sizes centered both on peaks of CO and Halpha emission. We parameterize this ratio as a molecular gas (H_2) depletion time (τ_dep). On large (kpc) scales, our results are consistent with a molecular star formation law (Sigma_SFR \sim Sigma_H2^b) with b \sim 1.1 - 1.5 and a median τ_dep \sim 1 Gyr, with no dependence on type of region targeted. Below these scales, τ_dep is a strong function of adopted angular scale and the type of region that is targeted. Small (\lesssim 300pc) apertures centered on CO peaks have very long τ_dep (i.e., high CO-to-Halpha flux ratio) and small apertures targeted toward Halpha peaks have very short τ_dep. This implies that the star formation law observed on kpc scales breaks down once one reaches aperture sizes of \lesssim 300pc. For our smallest apertures (75pc), the difference in τ_dep between the two types of regions is more than one order of magnitude. This scale behavior emerges from averaging over star-forming regions with a wide range of CO-to-Halpha ratios with the natural consequence that the breakdown in the star formation law is a function of the surface density of the regions studied. We consider the evolution of individual regions the most likely driver for region-to-region differences in τ_dep (and thus the CO-to-Halpha ratio).

preprint2008arXiv

The Star Formation Efficiency in Nearby Galaxies: Measuring Where Gas Forms Stars Effectively

We measure the star formation efficiency (SFE), the star formation rate per unit gas, in 23 nearby galaxies and compare it to expectations from proposed star formation laws and thresholds. We use HI maps from THINGS and derive H2 maps from HERACLES and BIMA SONG CO. We estimate the star formation rate by combining GALEX FUV maps and SINGS 24 micron maps, infer stellar surface density profiles from SINGS 3.6 micron data, and use kinematics from THINGS. We measure the SFE as a function of: the free-fall and orbital timescales; midplane gas pressure; stability of the gas disk to collapse (including the effects of stars); the ability of perturbations to grow despite shear; and the ability of a cold phase to form. In spirals, the SFE of H2 alone is nearly constant at 5.25 +/- 2.5 x 10^(-10) yr^(-1) (equivalent to an H2 depletion time of 1.9x10^9 yr) as a function of all of these variables at our 800 pc resolution. Where the ISM is mostly HI, on the other hand, the SFE decreases with increasing radius in both spiral and dwarf galaxies, a decline reasonably described by an exponential with scale length 0.2-0.25 r_25. We interpret this decline as a strong dependence of GMC formation on environment. The ratio of H2 to HI appears to be a smooth function of radius, stellar surface density, and pressure spanning from the H2-dominated to HI-dominated ISM. The radial decline in SFE is too steep to be reproduced only by increases in the free-fall time or orbital time. Thresholds for large-scale instability suggest that our disks are stable or marginally stable and do not show a clear link to the declining SFE. We suggest that ISM physics below the scales that we observe - phase balance in the HI, H2 formation and destruction, and stellar feedback - governs the formation of GMCs from HI.

preprint2006arXiv

The Spitzer Survey of the Small Magellanic Cloud: S3MC Imaging and Photometry in the Mid- and Far-Infrared Wavebands

We present the initial results from the Spitzer Survey of the Small Magellanic Cloud (S3MC), which imaged the star-forming body of the Small Magellanic Cloud (SMC) in all seven MIPS and IRAC wavebands. We find that the F_8/F_24 ratio (an estimate of PAH abundance) has large spatial variations and takes a wide range of values that are unrelated to metallicity but anticorrelated with 24 um brightness and F_24/F_70 ratio. This suggests that photodestruction is primarily responsible for the low abundance of PAHs observed in star-forming low-metallicity galaxies. We use the S3MC images to compile a photometric catalog of ~400,000 mid- and far-infrared point sources in the SMC. The sources detected at the longest wavelengths fall into four main categories: 1) bright 5.8 um sources with very faint optical counterparts and very red mid-infrared colors ([5.8]-[8.0]>1.2), which we identify as YSOs. 2) Bright mid-infrared sources with mildly red colors (0.16<[5.8]-[8.0]<0.6), identified as carbon stars. 3) Bright mid-infrared sources with neutral colors and bright optical counterparts, corresponding to oxygen-rich evolved stars. And, 4) unreddened early B stars (B3 to O9) with a large 24 um excess. This excess is reminiscent of debris disks, and is detected in only a small fraction of these stars (<5%). The majority of the brightest infrared point sources in the SMC fall into groups one to three. We use this photometric information to produce a catalog of 282 bright YSOs in the SMC with a very low level of contamination (~7%).

Adam K. Leroy

What is connected

Connect this record

See the researcher in context

Building this map preview

76 published item(s)

Evolved Supergiants in PHANGS I: Red Supergiants in 19 Galaxies between 5-20 Mpc with HST and JWST

A 2-3 mm high-resolution molecular line survey towards the centre of the nearby spiral galaxy NGC 6946

A CO isotopologue Line Atlas within the Whirlpool galaxy Survey (CLAWS)

A tale of two DIGs: The relative role of HII regions and low-mass hot evolved stars in powering the diffuse ionised gas (DIG) in PHANGS-MUSE galaxies

Environmental dependence of the molecular cloud lifecycle in 54 main sequence galaxies

The Morpho-Kinematic Architecture of Super Star Clusters in the Center of NGC253

The PHANGS-MUSE survey -- Probing the chemo-dynamical evolution of disc galaxies

Giant Molecular Cloud Catalogues for PHANGS-ALMA: Methods and Initial Results

Low-J CO Line Ratios From Single Dish CO Mapping Surveys and PHANGS-ALMA

Molecular Gas Properties and CO-to-H2 Conversion Factors in the Central Kiloparsec of NGC 3351

PHANGS-HST: Star Cluster Spectral Energy Distribution Fitting with CIGALE

PHANGS-MUSE: the HII region luminosity function of local star-forming galaxies

Resolving the Dust-to-Metals Ratio and CO-to-H$_2$ Conversion Factor in the Nearby Universe

The PHANGS-HST Survey: Physics at High Angular resolution in Nearby GalaxieS with the Hubble Space Telescope

A model for the onset of self-gravitation and star formation in molecular gas governed by galactic forces: II. the bottleneck to collapse set by cloud-environment decoupling

Detection of the diffuse HI emission in the Circumgalactic Medium of NGC 891 and NGC 4565

Distances to PHANGS Galaxies: New Tip of the Red Giant Branch Measurements and Adopted Distances

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

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

Spatial Power Spectra of Dust across the Local Group: No Constraint on Disc Scale Height

The evolution of neutral hydrogen over the past 11 Gyr via HI 21 cm absorption

The Molecular ISM in the Super Star Clusters of the Starburst NGC 253

The turbulent gas structure in the centers of NGC253 and the Milky Way

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

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

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

A Portrait of Cold Gas in Galaxies at 60pc Resolution and a Simple Method to Test Hypotheses That Link Small-Scale ISM Structure to Galaxy-Scale Processes

Dense Gas Fraction and Star Formation Efficiency Variations in the Antennae Galaxies

The Molecular Clouds Fueling a 1/5 Solar Metallicity Starburst

The Relationship Between Molecular Gas, HI, and Star Formation in the Low-Mass, Low-Metallicity Magellanic Clouds

ALMA Multi-line Imaging of the Nearby Starburst Galaxy NGC 253

Faint CO Line Wings in Four Star-Forming (Ultra)luminous Infrared Galaxies

Microwave Continuum Emission and Dense Gas Tracers in NGC 3627: Combining Jansky VLA and ALMA Observations

Next Generation Very Large Array Memo No. 7 Science Working Group 2: "Galaxy Ecosystems": The Matter Cycle in and Around Galaxies

Short GMC lifetimes: an observational estimate with the PdBI Arcsecond Whirlpool Survey (PAWS)

The Multi-Phase Cold Fountain in M82 Revealed by a Wide, Sensitive Map of the Molecular ISM

The Panchromatic Hubble Andromeda Treasury VIII: A Wide-Area, High-Resolution Map of Dust Extinction in M31

The Panchromatic Hubble Andromeda Treasury XI: The Spatially-Resolved Recent Star Formation History of M31

The Relationship between CO Emission and Visual Extinction Traced by Dust Emission in the Magellanic Clouds

Variations in the Star Formation Efficiency of the Dense Molecular Gas across the Disks of Star-Forming Galaxies

A Test of Star Formation Laws in Disk Galaxies. II. Dependence on dynamical properties

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

High-Resolution Radio Continuum Measurements of the Nuclear Disks of Arp 220

The PdBI Arcsecond Whirlpool Survey (PAWS): Environmental Dependence of Giant Molecular Cloud Properties in M51

The PdBI Arcsecond Whirlpool Survey (PAWS): Multi-phase cold gas kinematic of M51

The Radio Continuum-Star Formation Rate Relation in WSRT SINGS Galaxies

A Comparative Study of Giant Molecular Clouds in M51, M33 and the Large Magellanic Cloud

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING). III. The Dependence of Atomic and Molecular Gas Surface Densities on Galaxy Properties

Clumping and the Interpretation of kpc-Scale Maps of the Interstellar Medium: Smooth HI and Clumpy, Variable H2 Surface Density

Dust May Be More Rare Than Expected in Metal Poor Galaxies

Gas Kinematics on GMC scales in M51 with PAWS: cloud stabilization through dynamical pressure

Molecular Gas and Star Formation in Nearby Disk Galaxies

Probability Distribution Functions OF 12CO(J = 1-0) Brightness and Integrated Intensity in M51: The PAWS View

The CO-to-H2 Conversion Factor

The Fueling Diagram: Linking Galaxy Molecular-to-Atomic Gas Ratios to Interactions and Accretion

The Green Bank Telescope Maps the Dense, Star-Forming Gas in the Nearby Starburst Galaxy M82

The PdBI Arcsecond Whirlpool Survey (PAWS). I. A Cloud-Scale/Multi-Wavelength View of the Interstellar Medium in a Grand-Design Spiral Galaxy

The Starburst-Driven Molecular Wind in NGC 253 and the Suppression of Star Formation

A High Resolution Study of the HI-H2 Transition across the Perseus Molecular Cloud

Dust-to-Gas Ratio in the Extremely Metal Poor Galaxy I ZW 18

Estimating the Star Formation Rate at 1 kpc Scales in Nearby Galaxies

Low CO Luminosities in Dwarf Galaxies

The VIRUS-P Exploration of Nearby Galaxies (VENGA): The Xco Gradient in NGC 628

A Molecular Star Formation Law in the Atomic Gas Dominated Regime in Nearby Galaxies

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING) II: Molecular Gas Star Formation Law and Depletion Time Across the Blue Sequence

CARMA Survey Toward Infrared-bright Nearby Galaxies (STING): Molecular Gas Star Formation Law in NGC4254

Complex Radio Spectral Energy Distributions in Luminous and Ultraluminous Infrared Galaxies

Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview

The CO-to-H2 Conversion Factor From Infrared Dust Emission Across the Local Group

The Spitzer Spectroscopic Survey of the Small Magellanic Cloud (S4MC): Probing the Physical State of Polycyclic Aromatic Hydrocarbons in a Low-Metallicity Environment

The State of the Gas and the Relation Between Gas and Star Formation at Low Metallicity: the Small Magellanic Cloud

What Phase of the Interstellar Medium Correlates with the Star Formation Rate?

Regulation of Star Formation Rates in Multiphase Galactic Disks: a Thermal/Dynamical Equilibrium Model

The Scale Dependence of the Molecular Gas Depletion Time in M33