Source author record

Frank Bigiel

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

45works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

45 published item(s)

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

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

Planetary Nebula Luminosity Function distances for 19 galaxies observed by PHANGS-MUSE

We provide new planetary nebula luminosity function (PNLF) distances to 19 nearby spiral galaxies that were observed with VLT/MUSE by the PHANGS collaboration. Emission line ratios are used to separate planetary nebulae (PNe) from other bright [OIII] emitting sources like compact supernovae remnants (SNRs) or HII regions. While many studies have used narrowband imaging for this purpose, the detailed spectral line information provided by integral field unit (IFU) spectroscopy grants a more robust way of categorising different [OIII] emitters. We investigate the effects of SNR contamination on the PNLF and find that we would fail to classify all objects correctly, when limited to the same data narrowband imaging provides. However, the few misclassified objects usually do not fall on the bright end of the luminosity function, and only in three cases does the distance change by more than $1σ$. We find generally good agreement with literature values from other methods. Using metallicity constraints that have also been derived from the same IFU data, we revisit the PNLF zero point calibration. Over a range of $8.34 < 12 + \log(\mathrm{O}/\mathrm{H}) < 8.59$, our sample is consistent with a constant zero point and yields $M^*=-4.542^{+0.103}_{-0.059}\, \mathrm{mag}$, within $1σ$ of other literature values. MUSE pushes the limits of PNLF studies and makes galaxies beyond $20\, \mathrm{Mpc}$ accessible for this kind of analysis. This approach to the PNLF shows great promise for leveraging existing archival IFU data on nearby galaxies.

preprint2022arXiv

The Gas-Star Formation Cycle in Nearby Star-forming Galaxies II. Resolved Distributions of CO and H$α$ Emission for 49 PHANGS Galaxies

The relative distribution of molecular gas and star formation in galaxies gives insight into the physical processes and timescales of the cycle between gas and stars. In this work, we track the relative spatial configuration of CO and H$α$ emission at high resolution in each of our galaxy targets, and use these measurements to quantify the distributions of regions in different evolutionary stages of star formation: from molecular gas without star formation traced by H$α$ to star-forming gas, and to HII regions. The large sample, drawn from the Physics at High Angular resolution in Nearby GalaxieS ALMA and narrowband H$α$ (PHANGS-ALMA and PHANGS-H$α$) surveys, spans a wide range of stellar mass and morphological types, allowing us to investigate the dependencies of the gas-star formation cycle on global galaxy properties. At a resolution of 150 pc, the incidence of regions in different stages shows a dependence on stellar mass and Hubble type of galaxies over the radial range probed. Massive and/or earlier-type galaxies exhibit a significant reservoir of molecular gas without star formation traced by H$α$, while lower-mass galaxies harbor substantial HII regions that may have dispersed their birth clouds or formed from low-mass, more isolated clouds. Galactic structures add a further layer of complexity to relative distribution of CO and H$α$ emission. Trends between galaxy properties and distributions of gas traced by CO and H$α$ are visible only when the observed spatial scale is $\ll$ 500 pc, reflecting the critical resolution requirement to distinguish stages of star formation process.

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

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.

preprint2021arXiv

WALLABY Pre-Pilot Survey: The effects of tidal interaction on radial distribution of color in galaxies of the Eridanus supergroup

We study the tidal interaction of galaxies in the Eridanus supergroup, using HI data from the pre-pilot survey of WALLABY (Widefield ASKAP L-band Legacy All-sky Blind surveY). We obtain optical photometric measurements and quantify the strength of tidal perturbation using a tidal parameter $S_{sum}$. For low-mass galaxies of $M_* \lesssim 10^9 M_\odot$, we find a dependence of decreasing HI-to-optical disk size ratio with increasing $S_{sum}$, but no dependence of HI spectral line asymmetry with $S_{sum}$. This is consistent with the behavior expected under tidal stripping. We confirm that the color profile shape and color gradient depend on the stellar mass, but there is additional correlation of low-mass galaxies having their color gradients within $2R_{50}$ increasing with higher $S_{sum}$. For these low-mass galaxies, the dependence of color gradients on $S_{sum}$ is driven by color becoming progressively redder in the inner disk when tidal perturbations are stronger. For high-mass galaxies, there is no dependence of color gradients on $S_{sum}$, and we find a marginal reddening throughout the disks with increasing $S_{sum}$. Our result highlights tidal interaction as an important environmental effect in producing the faint end of the star formation suppressed sequence in galaxy groups.

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

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

When Gas Dynamics Decouples from Galactic Rotation: Characterizing ISM Circulation in Disk Galaxies

In galactic disks, galactic rotation sets the bulk motion of gas, and its energy and momentum can be transferred toward small scales. Additionally, in the interstellar medium, random and noncircular motions arise from stellar feedback, cloud-cloud interactions, and instabilities, among other processes. Our aim is to comprehend to which extent small-scale gas dynamics is decoupled from galactic rotation. We study the relative contributions of galactic rotation and local noncircular motions to the circulation of gas, $Γ$, a macroscopic measure of local rotation, defined as the line integral of the velocity field around a closed path. We measure the circulation distribution as a function of spatial scale in a set of simulated disk galaxies and we model the velocity field as the sum of galactic rotation and a Gaussian random field. The random field is parameterized by a broken power law in Fourier space, with a break at the scale $λ_c$. We define the spatial scale $λ_{\rm eq}$ at which galactic rotation and non-circular motions contribute equally to $Γ$. For our simulated galaxies, the gas dynamics at the scale of molecular clouds is usually dominated by noncircular motions, but in the center of galactic disks galactic rotation is still relevant. Our model shows that the transfer of rotation from large scales breaks at the scale $λ_c$ and this transition is necessary to reproduce the circulation distribution. We find that $λ_{\rm eq}$, and therefore the structure of the gas velocity field, is set by the local conditions of gravitational stability and stellar feedback.

preprint2019arXiv

Mapping electron temperature variations across a spiral arm in NGC 1672

We report one of the first extragalactic observations of electron temperature variations across a spiral arm. Using MUSE mosaic observations of the nearby galaxy NGC 1672, we measure the [N II]5755 auroral line in a sample of 80 HII regions in the eastern spiral arm of NGC1672. We discover systematic temperature variations as a function of distance perpendicular to the spiral arm. The electron temperature is lowest on the spiral arm itself and highest on the downstream side. Photoionization models of different metallicity, pressure, and age of the ionizing source are explored to understand what properties of the interstellar medium drive the observed temperature variations. An azimuthally varying metallicity appears to be the most likely cause of the temperature variations. The electron temperature measurements solidify recent discoveries of azimuthal variations of oxygen abundance based on strong lines, and rule out the possibility that the abundance variations are artefacts of the strong-line calibrations.

preprint2019arXiv

The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations

The processes regulating star formation in galaxies are thought to act across a hierarchy of spatial scales. To connect extragalactic star formation relations from global and kpc-scale measurements to recent cloud-scale resolution studies, we have developed a simple, robust method that quantifies the scale dependence of the relative spatial distributions of molecular gas and recent star formation. In this paper, we apply this method to eight galaxies with roughly 1 arcsec resolution molecular gas imaging from the PHANGS-ALMA and PAWS surveys that have matched resolution, high quality narrowband Halpha imaging. At a common scale of 140pc, our massive (log(Mstar/Msun)=9.3-10.7), normally star-forming (SFR/Msun/yr=0.3-5.9) galaxies exhibit a significant reservoir of quiescent molecular gas not associated with star formation as traced by Halpha emission. Galactic structures act as backbones for both molecular and HII region distributions. As we degrade the spatial resolution, the quiescent molecular gas disappears, with the most rapid changes occurring for resolutions up to about 0.5kpc. As the resolution becomes poorer, the morphological features become indistinct for spatial scales larger than about 1kpc. The method is a promising tool to search for relationships between the quiescent or star-forming molecular reservoir and galaxy properties, but requires a larger sample size to identify robust correlations between the star-forming molecular gas fraction and global galaxy parameters.

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

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

Radial gas motions in The HI Nearby Galaxy Survey (THINGS)

The study of 21cm line observations of atomic hydrogen allows detailed insight into the kinematics of spiral galaxies. We use sensitive high-resolution VLA data from The HI Nearby Galaxy Survey (THINGS) to search for radial gas flows primarily in the outer parts (up to $3\times r_{25}$) of ten nearby spiral galaxies. Inflows are expected to replenish the gas reservoir and fuel star formation under the assumption that galaxies evolve approximately in steady state. We carry out a detailed investigation of existing tilted ring fitting schemes and discover systematics that can hamper their ability to detect signatures of radial flows. We develop a new Fourier decomposition scheme that fits for rotational and radial velocities and simultaneously determines position angle and inclination as a function of radius. Using synthetic velocity fields we show that our novel fitting scheme is less prone to such systematic errors and that it is well suited to detect radial inflows in disks. We apply our fitting scheme to ten THINGS galaxies and find clear indications of, at least partly previously unidentified, radial gas flows, in particular for NGC 2403 and NGC 3198 and to a lesser degree for NGC 7331, NGC 2903 and NGC 6946. The mass flow rates are of the same order but usually larger than the star formation rates. At least for these galaxies a scenario in which continuous mass accretion feeds star formation seems plausible. The other galaxies show a more complicated picture with either no clear inflow, outward motions or complex kinematic signatures.

preprint2015arXiv

A Survey for Hydroxyl in the THOR Pilot Region around W43

We report on observations of the hydroxyl radical (OH) within The H{\sc I}, OH Recombination line survey (THOR) pilot region. The region is bounded approximately between Galactic coordinates l=29.2 to 31.5$^\circ$ and b=-1.0 to +1.0$^\circ$ and includes the high-mass star forming region W43. We identify 103 maser sites, including 72 with 1612\,MHz masers, 42 showing masers in either of the main line transitions at 1665 and 1667\,MHz and four showing 1720\,MHz masers. Most maser sites with either main-line or 1720\,MHz emission are associated with star formation, whereas most of the 1612\,MHz masers are associated with evolved stars. We find that nearly all of the main-line maser sites are co-spatial with an infrared source, detected by GLIMPSE. We also find diffuse OH emission, as well as OH in absorption towards selected unresolved or partially resolved sites. Extended OH absorption is found towards the well known star forming complex W43 Main.

preprint2015arXiv

Gas-phase metallicity profiles of the Bluedisk galaxies: Is metallicity in a local star-formation regulated equilibrium?

As part of the Bluedisk survey we analyse the radial gas-phase metallicity profiles of 50 late-type galaxies We compare the metallicity profiles of a sample of HI-rich galaxies against a control sample of HI-'normal' galaxies. We find the metallicity gradient of a galaxy to be strongly correlated with its HI mass fraction (M(HI) / Mstar). We note that some galaxies exhibit a steeper metallicity profile in the outer disc than in the inner disc. These galaxies are found in both the HI-rich and control samples. This contradicts a previous indication that these outer drops are exclusive to HI-rich galaxies. These effects are not driven by bars, although we do find some indication that barred galaxies have flatter metallicity profiles. By applying a simple analytical model we are able to account for the variety of metallicity profiles that the two samples present. The success of this model implies that the metallicity in these isolated galaxies may be in a local equilibrium, regulated by star formation. This insight could provide an explanation of the observed local mass-metallicity relation.

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

An observational and theoretical view of the radial distribution of HI gas in galaxies

We analyze the radial distribution of HI gas for 23 disk galaxies with unusually high HI content from the Bluedisk sample, along with a similar-sized sample of "normal" galaxies. We propose an empirical model to fit the radial profile of the HI surface density, an exponential function with a depression near the center. The radial HI surface density profiles are very homogeneous in the outer regions of the galaxy; the exponentially declining part of the profile has a scale-length of $\sim 0.18$ R1, where R1 is the radius where the column density of the HI is 1 M$_{\odot}$ pc$^{-2}$. This holds for all galaxies, independent of their stellar or HI mass. The homogenous outer profiles, combined with the limited range in HI surface density in the non-exponential inner disk, results in the well-known tight relation between HI size and HI mass. By comparing the radial profiles of the HI-rich galaxies with those of the control systems, we deduce that in about half the galaxies, most of the excess gas lies outside the stellar disk, in the exponentially declining outer regions of the HI disk. In the other half, the excess is more centrally peaked. We compare our results with existing smoothed-particle hydrodynamical simulations and semi-analytic models of disk galaxy formation in a $Λ$ Cold Dark Matter universe. Both the hydro simulations and the semi-analytic models reproduce the HI surface density profiles and the HI size-mass relation without further tuning of the simulation and model inputs. In the semi-analytic models, the universal shape of the outer HI radial profiles is a consequence of the {\em assumption} that infalling gas is always distributed exponentially. The conversion of atomic gas to molecular form explains the limited range of HI surface densities in the inner disk. These two factors produce the tight HI mass-size relation.

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

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

Indications of a sub-linear and non-universal Kennicutt-Schmidt relationship

We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population. Assuming standard conversion factors to estimate Sigma_SFR and Sigma_mol from the observations, we find that the KS parameters vary between galaxies, indicating that no universal relationship holds for all galaxies. The KS slope of the whole population is 0.76, with the 2sigma range extending from 0.58 to 0.94. These results imply that the molecular gas depletion time is not constant, but varies from galaxy to galaxy, and increases with the molecular gas surface density. Therefore, other galactic properties besides just Sigma_mol affect Sigma_SFR, such as the gas fraction or stellar mass. The non-universality of the KS relationship indicates that a comprehensive theory of star formation must take into account additional physical processes that may vary from galaxy to galaxy.

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

The Bluedisks project, a study of unusually HI-rich galaxies: I. HI Sizes and Morphology

We introduce the "Bluedisk" project, a large program at the Westerbork Synthesis Radio Telescope (WSRT) that has mapped the HI in a sample of 23 nearby galaxies with unusually high HI mass fractions, along with a similar-sized sample of control galaxies matched in stellar mass, size, inclination and redshift. This paper presents the sample selection, observational set-up, data reduction strategy, and a first analysis of the sizes and structural properties of the HI disks. We find that the HI-rich galaxies lie on the same HI mass versus HI size relation as normal spiral galaxies, extending it to total HI masses of $2 \times 10^{10} M_{\odot}$ and radii R1 of $\sim 100$ kpc (where R1 is defined as the radius where the HI column density reaches 1 $M_{\odot}$ pc$^{-2}$). HI-rich galaxies have significantly larger values of HI-to-optical size ratio at fixed stellar mass, concentration index, stellar and star formation rate surface density compared to the control sample. The disks of HI-rich galaxies are also significantly more clumpy (i.e. have higher HI Gini and $Δ$Area coefficient) than those of normal spirals. There is no evidence that the disks of HI-rich galaxies are more disturbed: HI-rich galaxies exhibit no difference with respect to control samples in their distributions of HI asymmetry indices or optical/HI disk position angle differences. In fact, the center of the HI distribution corresponds more closely with the center of the optical light in the HI-rich galaxies than in the controls. All these results argue against a scenario in which new gas has been brought in by mergers. It is possible that they may be more consistent with cooling from a surrounding quasi-static halo of warm/hot gas.

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

Evidence for a non-universal Kennicutt-Schmidt relationship using hierarchical Bayesian linear regression

For investigating the relationship between the star formation rate and gas surface density, we develop a Bayesian linear regression method that rigorously treats measurement uncertainties and accounts for hierarchical data structure. The hierarchical method simultaneously estimates the intercept, slope, and scatter about the regression line of each individual subject and the population. Using synthetic datasets, we demonstrate that the method recovers the underlying parameters of both the individuals and the population, especially when compared to commonly employed ordinary least squares techniques, such as the bisector fit. We apply the hierarchical method to estimate the Kennicutt-Schmidt (KS) parameters of a sample of spiral galaxies compiled by Bigiel et al. (2008). We find significant variation in the KS parameters, indicating that no single relationship holds for all galaxies. This suggests that the relationship between molecular gas and star formation differs between galaxies, possibly due to the influence of other properties within a given galaxy, such as metallicity, molecular gas fraction, or stellar mass. In four galaxies the slope estimates are sub-linear (at the 95% level). We estimate the mean index of the KS relationship for the population to be 0.84, with 95% range [0.63, 1.0]. For the galaxies with sub-linear KS relationships, a possible interpretation is that CO emission is tracing some molecular gas that is not directly associated with star formation. Equivalently, a sub-linear KS relationship may be indicative of an increasing gas depletion time at higher surface densities, as traced by CO emission. The Bayesian method can account for all sources of uncertainties, including variations in the conversion of observed intensities to star formation rates and gas surface densities, and is therefore well suited for a thorough statistical analysis of the KS law.

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

On the Origin of the Supergiant HI Shell and Putative Companion in NGC 6822

We present new Hubble Space Telescope Advanced Camera for Surveys imaging of six positions spanning 5.8 kpc of the HI major axis of the Local Group dIrr NGC 6822, including both the putative companion galaxy and the large HI hole. The resulting deep color magnitude diagrams show that NGC 6822 has formed >50% of its stars in the last ~5 Gyr. The star formation histories of all six positions are similar over the most recent 500 Myr, including low-level star formation throughout this interval and a weak increase in star formation rate during the most recent 50 Myr. Stellar feedback can create the giant HI hole, assuming that the lifetime of the structure is longer than 500 Myr; such long-lived structures have now been observed in multiple systems and may be the norm in galaxies with solid-body rotation. The old stellar populations (red giants and red clump stars) of the putative companion are consistent with those of the extended halo of NGC 6822; this argues against the interpretation of this structure as a bona fide interacting companion galaxy and against its being linked to the formation of the HI hole via an interaction. Since there is no evidence in the stellar population of a companion galaxy, the most likely explanation of the extended HI structure in NGC 6822 is a warped disk inclined to the line of sight.

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

Finding Dwarf Galaxies From Their Tidal Imprints

We describe ongoing work on a new method that allows one to determine the mass and relative position (in galactocentric radius and azimuth) of galactic companions purely from analysis of observed disturbances in gas disks. Here, we demonstrate the validity of this method, which we call Tidal Analysis, by applying it to local spirals with known optical companions, namely M51 and NGC 1512. These galaxies span the range from having a very low mass companion ($\sim$ one-hundredth the mass of the primary galaxy) to a fairly massive companion ($\sim$ one-third the mass of the primary galaxy). This approach has broad implications for many areas of astrophysics -- for the indirect detection of dark matter (or dark-matter dominated dwarf galaxies), and for galaxy evolution in its use as a decipher for the dynamical impact of satellites on galactic disks. Here, we provide a proof of principle of the method by applying it to infer and quantitatively characterize optically visible galactic companions of local spirals, from the analysis of observed disturbances in outer gas disks.

preprint2010arXiv

SLUG: A new way to Stochastically Light Up Galaxies

We present SLUG, a new code to "Stochastically Light Up Galaxies". SLUG populates star clusters by randomly drawing stars from an initial mass function (IMF) and then following their time evolution with stellar models and an observationally-motivated prescription for cluster disruption. For a choice of star formation history, metallicity, and IMF, SLUG outputs synthetic photometry for clusters and field stars with a proper treatment of stochastic star formation. SLUG generates realistic distributions of star clusters, demonstrating the range of properties that result from finite sampling of an IMF and a random distribution of ages. The simulated data sets provide a quantitative means to address open problems in studies of star formation in galaxies and clusters, such as a test for IMF variations that are suggested by the systematic deficiency in the H-alpha/UV ratio in outer disks or in dwarf galaxies. SLUG will be made publicly available through the website http://sites.google.com/site/runslug/.

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.

preprint2008arXiv

The Star Formation Law in Nearby Galaxies on Sub-Kpc Scales

(Abridged) We present a comprehensive analysis of the relationship between star formation rate surface density (SFR SD) and gas surface density (gas SD) at sub-kpc resolution in a sample of 18 nearby galaxies. We use high resolution HI data from THINGS, CO data from HERACLES and BIMA SONG, 24 micron data from the Spitzer Space Telescope, and UV data from GALEX. We target 7 spiral galaxies and 11 late-type/dwarf galaxies and investigate how the star formation law differs between the H2-dominated centers of spiral galaxies, their HI-dominated outskirts and the HI-rich late-type/dwarf galaxies. We find that a Schmidt-type power law with index N=1.0+-0.2 relates the SFR SD and the H2 SD across our sample of spiral galaxies, i.e., that H2 forms stars at a constant efficiency in spirals. The average molecular gas depletion time is ~2*10^9 yrs. We interpret the linear relation and constant depletion time as evidence that stars are forming in GMCs with approximately uniform properties and that the H2 SD may be more a measure of the filling fraction of giant molecular clouds than changing conditions in the molecular gas. The relationship between total gas SD and SFR SD varies dramatically among and within spiral galaxies. Most galaxies show little or no correlation between the HI SD and the SFR SD. As a result, the star formation efficiency (SFE = SFR SD / gas SD) varies strongly across our sample and within individual galaxies. We show that in spirals the SFE is a clear function of radius, while the dwarf galaxies in our sample display SFEs similar to those found in the outer optical disks of the spirals. Another general feature of our sample is a sharp saturation of the HI SD at ~9 M_sol/pc^2 in both the spiral and dwarf galaxies.