Source author record

Erik Rosolowsky

Erik Rosolowsky 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.SR astro-ph.IM astro-ph astro-ph.HE gr-qc nucl-th

Catalog footprint

What is connected

85works

8topics

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.

preprint2022arXiv

12CO (3-2) High-Resolution Survey (COHRS) of the Galactic Plane: Complete Data Release

We present the full data release of 12CO (3-2) High-Resolution Survey (COHRS), which has mapped the inner Galactic plane over the range of 9.5$^{\circ}$ $\le$ l $\le$ 62.3$^{\circ}$ and $|b| \le 0.5^{\circ}$. The COHRS has been carried out using the Heterodyne Array Receiver Program (HARP) on the 15 m James Clerk Maxwell Telescope (JCMT) in Hawaii. The released data are smoothed to have a spatial resolution of 16.6 arcsec and a velocity resolution of 0.635 km/s, achieving a mean root-mean-square of $\sim 0.6$ K on $T_\mathrm{A}^*$. The COHRS data are useful for investigating detailed three-dimensional structures of individual molecular clouds and large-scale structures such as spiral arms in the Galactic plane. Furthermore, data from other available public surveys of different CO isotopologues and transitions with similar angular resolutions to this survey, such as FUGIN, SEDIGISM, and CHIMPS/CHIMPS2, allow studying the physical properties of molecular clouds and comparing their states with each other. In this paper, we report further observations on R2 and improved data reduction since the original COHRS release. We discuss the characteristics of the COHRS data and present integrated-emission images and a position-velocity (PV) map of the region covered. The PV map shows a good match with the spiral-arm traces from the existing CO and HI surveys. We also obtain and compare integrated one-dimensional distributions of 12CO (1-0) and (3-2) and those of star-forming populations to each other.

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

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores.

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

pyspeckit: A spectroscopic analysis and plotting package

pyspeckit is a toolkit and library for spectroscopic analysis in Python. We describe the pyspeckit package and highlight some of its capabilities, such as interactively fitting a model to data, akin to the historically widely-used splot function in IRAF. pyspeckit employs the Levenberg-Marquardt optimization method via the mpfit and lmfit implementations, and important assumptions regarding error estimation are described here. Wrappers to use pymc and emcee as optimizers are provided. A parallelized wrapper to fit lines in spectral cubes is included. As part of the astropy affiliated package ecosystem, pyspeckit is open source and open development and welcomes input and collaboration from the community.

preprint2022arXiv

Stellar feedback in M83 as observed with MUSE -- I. Overview, an unprecedented view of the stellar and gas kinematics and evidence of outflowing gas

We present a large VLT/MUSE mosaic (3.8 x 3.8 kpc) of the nearby spiral galaxy M83, with a spatial resolution ~20 pc. We obtained the kinematics of the stars and ionised gas, and compared them with molecular gas kinematics from ALMA CO(2-1). We separated the ionised gas into HII regions and diffuse ionised gas (DIG) and determined the fraction of Ha luminosity originating from the DIG (f_DIG). We observe that both stars and gas trace the galactic disk rotation, as well as a fast-rotating nuclear component, likely connected to secular processes driven by the galactic bar. In the gas kinematics, we observe a stream east of the nucleus, redshifted with respect to the disk. The stream is surrounded by an extended ionised gas region with enhanced velocity dispersion and a high ionisation state, which is largely consistent with being ionised by slow shocks. We interpret this feature as either the superposition of the disk and an extraplanar layer of DIG, or as a bar-driven inflow of shocked gas. A double Gaussian component fit to the Ha line also reveals the presence of a nuclear biconic structure whose axis of symmetry is perpendicular to the bar. The two cones appear blue- and redshifted along the line of sight and stand out for having an Ha emission separated by up to 200 km s-1 from that of the disk, and a high velocity dispersion ~80-200 km s-1. At the far end of the cones, we observe that the gas is consistent with being ionised by shocks. These features had never been observed before in M83; we postulate that they are tracing a starburst-driven outflow shocking into the surrounding ISM. Finally, we obtain f_DIG ~ 13% in our field of view. We inspect the emission of the HII regions and DIG in `BPT' diagrams, finding that in HII regions photoionisation accounts for 99.8% of the Ha flux, whereas the DIG has a mixed contribution from photoionisation (94.9%) and shocks (5.1%). [abridged]

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 Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) II. The Spatially Resolved Recent Star Formation History of M33

We measure the spatially resolved recent star formation history (SFH) of M33 using optical images taken with the Hubble Space Telescope as part of the Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) survey. The area covered by the observations used in this analysis covers a de-projected area of $\sim$38 kpc$^{2}$ and extends to $\sim$3.5 and $\sim$2 kpc from the center of M33 along the major and semi-major axes, respectively. We divide the PHATTER optical survey into 2005 regions that measure 24 arcsec, $\sim$100 pc, on a side and fit color magnitude diagrams for each region individually to measure the spatially resolved SFH of M33 within the PHATTER footprint. There are significant fluctuations in the SFH on small spatial scales and also galaxy-wide scales that we measure back to about 630 Myr ago. We observe a more flocculent spiral structure in stellar populations younger than about 80 Myr, while the structure of the older stellar populations is dominated by two spiral arms. We also observe a bar in the center of M33, which dominates at ages older than about 80 Myr. Finally, we find that the mean star formation rate (SFR) over the last 100 Myr within the PHATTER footprint is 0.32$\pm$0.02 M$_{\odot}$ yr$^{-1}$. We measure a current SFR (over the last 10 Myr) of 0.20$\pm$0.03 M$_{\odot}$ yr$^{-1}$. This SFR is slightly higher than previous measurements from broadband estimates, when scaled to account for the fraction of the D25 area covered by the PHATTER survey footprint.

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.

preprint2020arXiv

[CI](1-0) and [CI](2-1) in resolved local galaxies

We present resolved [CI] line intensities of 18 nearby galaxies observed with the SPIRE FTS spectrometer on the Herschel Space Observatory. We use these data along with resolved CO line intensities from $J_\mathrm{up} = 1$ to 7 to interpret what phase of the interstellar medium the [CI] lines trace within typical local galaxies. A tight, linear relation is found between the intensities of the CO(4-3) and [CI](2-1) lines; we hypothesize this is due to the similar upper level temperature of these two lines. We modeled the [CI] and CO line emission using large velocity gradient models combined with an empirical template. According to this modeling, the [CI](1-0) line is clearly dominated by the low-excitation component. We determine [CI] to molecular mass conversion factors for both the [CI](1-0) and [CI](2-1) lines, with mean values of $α_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $α_{\mathrm{[CI](2-1)}} = 34 $ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ with logarithmic root-mean-square spreads of 0.20 and 0.32 dex, respectively. The similar spread of $α_{\mathrm{[CI](1-0)}}$ to $α_{\mathrm{CO}}$ (derived using the CO(2-1) line) suggests that [CI](1-0) may be just as good a tracer of cold molecular gas as CO(2-1) in galaxies of this type. On the other hand, the wider spread of $α_{\mathrm{[CI](2-1)}}$ and the tight relation found between [CI](2-1) and CO(4-3) suggest that much of the [CI](2-1) emission may originate in warmer molecular gas.

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

Efficient methanol production on the dark side of a prestellar core

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

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

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

The role played by magnetic field during star formation is an important topic in astrophysics. We investigate the correlation between the orientation of star-forming cores (as defined by the core major axes) and ambient magnetic field directions in 1) a 3D MHD simulation, 2) synthetic observations generated from the simulation at different viewing angles, and 3) observations of nearby molecular clouds. We find that the results on relative alignment between cores and background magnetic field in synthetic observations slightly disagree with those measured in fully 3D simulation data, which is partly because cores identified in projected 2D maps tend to coexist within filamentary structures, while 3D cores are generally more rounded. In addition, we examine the progression of magnetic field from pc- to core-scale in the simulation, which is consistent with the anisotropic core formation model that gas preferably flow along the magnetic field toward dense cores. When comparing the observed cores identified from the GBT Ammonia Survey (GAS) and Planck polarization-inferred magnetic field orientations, we find that the relative core-field alignment has a regional dependence among different clouds. More specifically, we find that dense cores in the Taurus molecular cloud tend to align perpendicular to the background magnetic field, while those in Perseus and Ophiuchus tend to have random (Perseus) or slightly parallel (Ophiuchus) orientations with respect to the field. We argue that this feature of relative core-field orientation could be used to probe the relative significance of the magnetic field within the cloud.

preprint2020arXiv

The MUSTANG-2 Galactic Plane Survey (MGPS90) pilot

We report the results of a pilot program for a Green Bank Telescope (GBT) MUSTANG Galactic Plane survey at 3 mm (90 GHz), MGPS90. The survey achieves a typical $1σ$ depth of $1-2$ mJy beam$^{-1}$ with a 9" beam. We describe the survey parameters, quality assessment process, cataloging, and comparison with other data sets. We have identified 709 sources over seven observed fields selecting some of the most prominent millimeter-bright regions between $0°< \ell < 50°$ (total area $\approx 7.5 °^2$). The majority of these sources have counterparts at other wavelengths. By applying flux selection criteria to these sources, we successfully recovered several known hypercompact HII (HCHII) regions, but did not confirm any new ones. We identify 126 sources that have mm-wavelength counterparts but do not have cm-wavelength counterparts and are therefore candidate HCHII regions; of these, 10 are morphologically compact and are strong candidates for new HCHII regions. Given the limited number of candidates in the extended area in this survey compared to the relatively large numbers seen in protoclusters W51 and W49, it appears that most HCHII regions exist within dense protoclusters. Comparing the counts of HCHII to ultracompact HII (UCHII) regions, we infer the HCHII region lifetime is 16-46% that of the UCHII region lifetime. We additionally separated the 3 mm emission into dust and free-free emission by comparing with archival 870 $μ$m and 20 cm data. In the selected pilot fields, most ($\gtrsim80$%) of the 3 mm emission comes from plasma, either through free-free or synchrotron emission.

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.

preprint2020arXiv

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Context : Star formation takes place in cold dense cores in molecular clouds. Earlier observations have found that dense cores exhibit subsonic non-thermal velocity dispersions. In contrast, CO observations show that the ambient large-scale cloud is warmer and has supersonic velocity dispersions. Aims : We aim to study the ammonia ($\rm NH_3$) molecular line profiles with exquisite sensitivity towards the coherent cores in L1688 in order to study their kinematical properties in unprecedented detail. Methods : We used $\rm NH_3$ (1,1) and (2,2) data from the first data release (DR1) in the Green Bank Ammonia Survey (GAS). We first smoothed the data to a larger beam of 1' to obtain substantially more extended maps of velocity dispersion and kinetic temperature, compared to the DR1 maps. We then identified the coherent cores in the cloud and analysed the averaged line profiles towards the cores. Results : For the first time, we detected a faint (mean $\rm NH_3$(1,1) peak brightness $<$0.25 K in $T_{MB}$), supersonic component towards all the coherent cores in L1688. We fitted two components, one broad and one narrow, and derived the kinetic temperature and velocity dispersion of each component. The broad components towards all cores have supersonic linewidths ($\mathcal{M}_S \ge 1$). This component biases the estimate of the narrow dense core component's velocity dispersion by $\approx$28% and the kinetic temperature by $\approx$10%, on average, as compared to the results from single-component fits. Conclusions : Neglecting this ubiquitous presence of a broad component towards all coherent cores causes the typical single-component fit to overestimate the temperature and velocity dispersion. This affects the derived detailed physical structure and stability of the cores estimated from $\rm NH_3$ observations.

preprint2020arXiv

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

Recent observations of global velocity gradients across and along molecular filaments have been interpreted as signs of gas accreting onto and along these filaments, potentially feeding star-forming cores and proto-clusters. The behavior of velocity gradients in filaments, however, has not been studied in detail, particularly on small scales (< 0.1 pc). In this paper, we present MUFASA, an efficient, robust, and automatic method to fit ammonia lines with multiple velocity components, generalizable to other molecular species. We also present CRISPy, a Python package to identify filament spines in 3D images (e.g., position-position-velocity cubes), along with a complementary technique to sort fitted velocity components into velocity-coherent filaments. In NGC 1333, we find a wealth of velocity gradient structures on a beam-resolved scale of ~0.05 pc. Interestingly, these local velocity gradients are not randomly oriented with respect to filament spines and their perpendicular, i.e., radial, component decreases in magnitude towards the spine for many filaments. Together with remarkably constant velocity gradients on larger scales along many filaments, these results suggest a scenario in which gas falling onto filaments is progressively damped and redirected to flow along these filaments.

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

Canada and the SKA from 2020-2030

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

preprint2019arXiv

Canadian Investigations of the Interstellar Medium

The interstellar medium mediates galactic evolution as the reservoir of material for future star formation and the repository of energy and matter output by stellar processes. Canadians have played leading roles in ISM science for decades. The Canadian Galactic Plane Survey identified a wealth of small-scale structure in H I emission as well as self-absorption and in the structure of polarized emission. These observations demonstrated that no phase of the ISM, including the transition from atomic gas to star formation, can be understood in isolation. Canadians have also played leading roles in the characterization of dust with Planck and balloon-borne telescopes. Canadians have also used pulsar scintillometry to measure structure in the ISM at the smallest scales, below 1 AU. The 2020s offer many opportunities for ISM science in Canada. A major but cost-effective upgrade to the Synthesis Telescope with broadband (400-1800 MHz) single-pixel feeds would enable broadband polarimetry as well as wide-area, arcminute surveys of radio recombination lines. The next generation of balloon-borne telescopes will investigate magnetic fields and dust properties. Large single dishes, particularly the Green Bank Telescope, remain essential for our understanding of the diffuse structure which characterizes the ISM. Very long baseline interferometry capability enables parallax measurements of pulsars and masers and for further progress in scintillometry. ISM astronomers will continue to participate in cosmological experiments including CHIME and CHORD. Protecting quiet radio frequency interference environments will be ever more important as broadband observations are ever more central to ISM science. Computational capability is essential both for numerical work and for handling the observational data.

preprint2019arXiv

CLOVER: Convnet Line-fitting Of Velocities in Emission-line Regions

When multiple star-forming gas structures overlap along the line-of-sight and emit optically thin emission at significantly different radial velocities, the emission can become non-Gaussian and often exhibits two distinct peaks. Traditional line-fitting techniques can fail to account adequately for these double-peaked profiles, providing inaccurate cloud kinematics measurements. We present a new method called Convnet Line-fitting Of Velocities in Emission-line Regions (CLOVER) for distinguishing between one-component, two-component, and noise-only emission lines using 1D convolutional neural networks trained with synthetic spectral cubes. CLOVER utilizes spatial information in spectral cubes by predicting on $3\times3$ pixel sub-cubes, using both the central pixel's spectrum and the average spectrum over the $3\times3$ grid as input. On an unseen set of 10,000 synthetic spectral cubes in each predicted class, CLOVER has classification accuracies of $\sim99\%$ for the one-component class and $\sim97\%$ for the two-component class. For the noise-only class, which is analogous to a signal-to-noise cutoff of four for traditional line-fitting methods, CLOVER has classification accuracy of $100\%$. CLOVER also has exceptional performance on real observations, correctly distinguishing between the three classes across a variety of star-forming regions. In addition, CLOVER quickly and accurately extracts kinematics directly from spectra identified as two-component class members. Moreover, we show that CLOVER is easily scalable to emission lines with hyperfine splitting, making it an attractive tool in the new era of large-scale NH$_3$ and N$_2$H$^+$ mapping surveys.

preprint2019arXiv

Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

preprint2019arXiv

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $β$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $β\sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

preprint2019arXiv

LRP2020: The cosmic origin and evolution of the elements

The origin of many elements of the periodic table remains an unsolved problem. While many nucleosynthetic channels are broadly understood, significant uncertainties remain regarding certain groups of elements such as the intermediate and rapid neutron-capture processes, the p-process, or the origin of odd-Z elements in the most metal-poor stars. Canada has a long tradition of leadership in nuclear astrophysics, dating back to the work of Alastair Cameron in the 1950s. Recent faculty hires have further boosted activity in the field, including transient observation and theory, survey science on galactic nucleosynthesis, and nuclear experiments. This white paper contains a brief overview of recent activity in the community, highlighting strengths in each sub-field, and provides recommendations to improve interdisciplinary collaboration. Sustaining Canadian leadership in the next decade will require, on the observational side, access to transient and non-transient surveys like LSST, SKA, or MSE, support for target-of-opportunity observing in current and future Canadian telescopes, and participation in next-generation X-ray telescopes such as ATHENA. State-of-the-art theoretical predictions will require an ambitious succession plan for the Niagara supercomputer to support large parallel jobs. We propose a funding instrument for postdoctoral training that reflects the interdisciplinary nature of nuclear astrophysics research, and the creation of a national collaborative funding program that allows for joint projects and workshop organization.

preprint2019arXiv

The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

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.

preprint2019arXiv

The Next Generation Very Large Array

The next generation Very Large Array (ngVLA) is a transformational radio observatory being designed by the U.S. National Radio Astronomy Observatory (NRAO). It will provide order of magnitude improvements in sensitivity, resolution, and uv coverage over the current Jansky Very Large Array (VLA) at ~1.2-50 GHz and extend the frequency range up to 70-115 GHz. This document is a white paper written by members of the Canadian community for the 2020 Long Range Plan panel, which will be making recommendations on Canada's future directions in astronomy. Since Canadians have been historically major users of the VLA and have been valued partners with NRAO for ALMA, Canada's participation in ngVLA is welcome. Canadians have been actually involved in ngVLA discussions for the past five years, and have played leadership roles in the ngVLA Science and Technical Advisory Councils. Canadian technologies are also very attractive for the ngVLA, in particular our designs for radio antennas, receivers, correlates, and data archives, and our industrial capacities to realize them. Indeed, the Canadian designs for the ngVLA antennas and correlator/beamformer are presently the baseline models for the project. Given the size of Canada's radio community and earlier use of the VLA (and ALMA), we recommend Canadian participation in the ngVLA at the 7% level. Such participation would be significant enough to allow Canadian leadership in gVLA's construction and usage. Canada's participation in ngVLA should not preclude its participation in SKA; access to both facilities is necessary to meet Canada's radio astronomy needs. Indeed, ngVLA will fill the gap between those radio frequencies observable with the SKA and ALMA at high sensitivities and resolutions. Canada's partnership in ngVLA will give it access to cutting-edge facilities together covering approximately three orders of magnitude in frequency.

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

Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations

We have developed an iterative procedure to systematically combine the millimeter and submillimeter images of OB cluster-forming molecular clouds, which were taken by ground based (CSO, JCMT, APEX, IRAM-30m) and space telescopes (Herschel, Planck). For the seven luminous ($L$$>$10$^{6}$ $L_{\odot}$) Galactic OB cluster-forming molecular clouds selected for our analyses, namely W49A, W43-Main, W43-South, W33, G10.6-0.4, G10.2-0.3, G10.3-0.1, we have performed single-component, modified black-body fits to each pixel of the combined (sub)millimeter images, and the Herschel PACS and SPIRE images at shorter wavelengths. The $\sim$10$"$ resolution dust column density and temperature maps of these sources revealed dramatically different morphologies, indicating very different modes of OB cluster-formation, or parent molecular cloud structures in different evolutionary stages. The molecular clouds W49A, W33, and G10.6-0.4 show centrally concentrated massive molecular clumps that are connected with approximately radially orientated molecular gas filaments. The W43-Main and W43-South molecular cloud complexes, which are located at the intersection of the Galactic near 3-kpc (or Scutum) arm and the Galactic bar, show a widely scattered distribution of dense molecular clumps/cores over the observed $\sim$10 pc spatial scale. The relatively evolved sources G10.2-0.3 and G10.3-0.1 appear to be affected by stellar feedback, and show a complicated cloud morphology embedded with abundant dense molecular clumps/cores. We find that with the high angular resolution we achieved, our visual classification of cloud morphology can be linked to the systematically derived statistical quantities (i.e., the enclosed mass profile, the column density probability distribution function, the two-point correlation function of column density, and the probability distribution function of clump/core separations).

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

Effects of Galactic Disc Inclination and Resolution on Observed GMC Properties and Larson's Scaling Relations

With ALMA making it possible to resolve giant molecular clouds (GMCs) in other galaxies, it is becoming necessary to quantify the observational bias on measured GMC properties. Using a hydrodynamical simulation of a barred spiral galaxy, we compared the physical properties of GMCs formed in position-position-position space (PPP) to the observational position-position-velocity space (PPV). We assessed the effect of disc inclination: face-on (PPV_face) and edge-on (PPV_edge), and resolution: 1.5 pc versus 24 pc, on GMC properties and the further implications of using Larson's scaling relations for mass-radius and velocity dispersion-radius. The low-resolution PPV data are generated by simulating ALMA Cycle 3 observations using the CASA package. Results show that the median properties do not differ strongly between PPP and PPV_face under both resolutions, but PPV_edge clouds deviate from these two. The differences become magnified when switching to the lower, but more realistic resolution. The discrepancy can lead to opposite results for the virial parameter's measure of gravitational binding, and therefore the dynamical state of the clouds. The power-law indices for the two Larson's scaling relations decrease going from PPP, PPV_face to PPV_edge and decrease from high to low resolutions. We conclude that the relations are not entirely driven by the underlying physical origin and therefore have to be used with caution when considering the environmental dependence, dynamical state, and the extragalactic CO-to-H2 conversion factor of GMCs.

preprint2016arXiv

The inner mass distribution of late-type spiral galaxies from SAURON stellar kinematic maps

We infer the central mass distributions within 0.4-1.2 disc scale lengths of 18 late-type spiral galaxies using two different dynamical modelling approaches - the Asymmetric Drift Correction (ADC) and axisymmetric Jeans Anisotropic Multi-gaussian expansion (JAM) model. ADC adopts a thin disc assumption, whereas JAM does a full line-of-sight velocity integration. We use stellar kinematics maps obtained with the integral-field spectrograph SAURON to derive the corresponding circular velocity curves from the two models. To find their best-fit values, we apply Markov Chain Monte Carlo (MCMC) method. ADC and JAM modelling approaches are consistent within 5% uncertainty when the ordered motions are significant comparable to the random motions, i.e, $\overline{v_ϕ}/σ_R$ is locally greater than 1.5. Below this value, the ratio $v_\mathrm{c,JAM}/v_\mathrm{c,ADC}$ gradually increases with decreasing $\overline{v_ϕ}/σ_R$, reaching $v_\mathrm{c,JAM}\approx 2 \times v_\mathrm{c,ADC}$. Such conditions indicate that the stellar masses of the galaxies in our sample are not confined to their disk planes and likely have a non-negligible contribution from their bulges and thick disks.

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

Giant Molecular Clouds in the Early-Type Galaxy NGC4526

We present a high spatial resolution ($\approx 20$ pc) of $^{12}$CO($2-1$) observations of the lenticular galaxy NGC4526. We identify 103 resolved Giant Molecular Clouds (GMCs) and measure their properties: size $R$, velocity dispersion $σ_v$, and luminosity $L$. This is the first GMC catalog of an early-type galaxy. We find that the GMC population in NGC4526 is gravitationally bound, with a virial parameter $α\sim 1$. The mass distribution, $dN/dM \propto M^{-2.39 \pm 0.03}$, is steeper than that for GMCs in the inner Milky Way, but comparable to that found in some late-type galaxies. We find no size-linewidth correlation for the NGC4526 clouds, in contradiction to the expectation from Larson's relation. In general, the GMCs in NGC4526 are more luminous, denser, and have a higher velocity dispersion than equal size GMCs in the Milky Way and other galaxies in the Local Group. These may be due to higher interstellar radiation field than in the Milky Way disk and weaker external pressure than in the Galactic center. In addition, a kinematic measurement of cloud rotation shows that the rotation is driven by the galactic shear. For the vast majority of the clouds, the rotational energy is less than the turbulent and gravitational energy, while the four innermost clouds are unbound and will likely be torn apart by the strong shear at the galactic center. We combine our data with the archival data of other galaxies to show that the surface density $Σ$ of GMCs is not approximately constant as previously believed, but varies by $\sim 3$ orders of magnitude. We also show that the size and velocity dispersion of GMC population across galaxies are related to the surface density, as expected from the gravitational and pressure equilibrium, i.e. $σ_v R^{-1/2} \propto Σ^{1/2}$.

preprint2015arXiv

Graph-based interpretation of the Molecular Interstellar Medium Segmentation

We present a generalization of the Giant Molecular Cloud (GMC) identification problem based on cluster analysis. The method we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation) considers the dendrogram of emission in the broader framework of graph theory and utilizes spectral clustering to find discrete regions with similar emission properties. For Galactic molecular cloud structures, we show that the characteristic volume and/or integrated CO luminosity are useful criteria to define the clustering, yielding emission structures that closely reproduce "by-eye" identification results. SCIMES performs best on well-resolved, high-resolution data, making it complementary to other available algorithms. Using 12CO(1-0) data for the Orion-Monoceros complex, we demonstrate that SCIMES provides robust results against changes of the dendrogram-construction parameters, noise realizations and degraded resolution. By comparing SCIMES with other cloud decomposition approaches, we show that our method is able to identify all canonical clouds of the Orion-Monoceros region, avoiding the over-division within high resolution survey data that represents a common limitation of several decomposition algorithms. The Orion-Monoceros objects exhibit hierarchies and size-line width relationships typical to the turbulent gas in molecular clouds, although "the Scissors" region deviates from this common description. SCIMES represents a significant step forward in moving away from pixel-based cloud segmentation towards a more physical-oriented approach, where virtually all properties of the ISM can be used for the segmentation of discrete objects.

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

The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations

We present 107 maps of continuum emission at 350 microns from Galactic molecular clumps. Observed sources were mainly selected from the Bolocam Galactic Plane Survey (BGPS) catalog, with 3 additional maps covering star forming regions in the outer Galaxy. The higher resolution of the SHARC-II images (8.5" beam) compared with the 1.1 mm images from BGPS (33" beam) allowed us to identify a large population of smaller substructures within the clumps. A catalog is presented for the 1386 sources extracted from the 350 micron maps. The color temperature distribution of clumps based on the two wavelengths has a median of 13.3 K and mean of 16.3 +- 0.4 K, assuming an opacity law index of 1.7. For the structures with good determination of color temperatures, the mean ratio of gas temperature, determined from NH3 observations, to dust color temperature is 0.88 and the median ratio is 0.76. About half the clumps have more than two substructures and 22 clumps have more than 10. The fraction of the mass in dense substructures seen at 350 microns compared to the mass of their parental clump is ~0.19, and the surface densities of these substructures are, on average, 2.2 times those seen in the clumps identified at 1.1 mm. For a well-characterized sample, 88 structures (31%) exceed a surface density of 0.2 g cm^(-2), and 18 (6%) exceed 1.0 g cm^(-2), thresholds for massive star formation suggested by theorists.

preprint2015arXiv

The Bolocam Galactic Plane Survey. XIII. Physical Properties and Mass Functions of Dense Molecular Cloud Structures

We use the distance probability density function (DPDF) formalism of Ellsworth-Bowers et al. (2013, 2015) to derive physical properties for the collection of 1,710 Bolocam Galactic Plane Survey (BGPS) version 2 sources with well-constrained distance estimates. To account for Malmquist bias, we estimate that the present sample of BGPS sources is 90% complete above 400 $M_\odot$ and 50% complete above 70 $M_\odot$. The mass distributions for the entire sample and astrophysically motivated subsets are generally fitted well by a lognormal function, with approximately power-law distributions at high mass. Power-law behavior emerges more clearly when the sample population is narrowed in heliocentric distance (power-law index $α= 2.0\pm0.1$ for sources nearer than 6.5 kpc and $α= 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq α\leq 2.05$ for various subsamples of sources, intermediate between that of giant molecular clouds and the stellar initial mass function. The fit to the entire sample yields a high-mass power-law $\hatα = 1.94_{-0.10}^{+0.34}$. Physical properties of BGPS sources are consistent with large molecular cloud clumps or small molecular clouds, but the fractal nature of the dense interstellar medium makes difficult the mapping of observational categories to the dominant physical processes driving the observed structure. The face-on map of the Galactic disk's mass surface density based on BGPS dense molecular cloud structures reveals the high-mass star-forming regions W43, W49, and W51 as prominent mass concentrations in the first quadrant. Furthermore, we present a 0.25-kpc resolution map of the dense gas mass fraction across the Galactic disk that peaks around 5%.

preprint2015arXiv

The dense gas mass fraction in the W51 cloud and its protoclusters

We present new 2 cm and 6 cm maps of H2CO, radio recombination lines, and the radio continuum in the W51 star forming complex acquired with Arecibo and the Green Bank Telescope at ~50" resolution. We use H2CO absorption to determine the relative line-of-sight positions of molecular and ionized gas. We measure gas densities using the H2CO densitometer, including continuous measurements of the dense gas mass fraction (DGMF) over the range $10^4$ cm$^{-3}$ < n(H$_2$) < $10^6$ cm$^{-3}$- this is the first time a dense gas mass fraction has been measured over a range of densities with a single data set. The DGMF in W51A is high,f >~70% above $n>10^4$ cm$^{-3}$, while it is low, f<20%, in W51 B. We did not detect any H2CO emission throughout the W51 GMC; all gas dense enough to emit under normal conditions is in front of bright continuum sources and therefore is seen in absorption instead. The data set has been made public at http://dx.doi.org/10.7910/DVN/26818. Conclusions. (1) The dense gas fraction in the W51 A and B clouds shows that W51 A will continue to form stars vigorously, while star formation has mostly ended in W51 B. The lack of dense, star-forming gas around W51 C indicates that collect-and-collapse is not acting or is inefficient in W51. (2) Ongoing high-mass star formation is correlated with n ~ 1x10$^5$ cm$^{-3}$ gas. Gas with n > 10$^4$ cm$^{-3}$ is weakly correlated with low and moderate mass star formation, but does not strongly correlate with high-mass star formation. (3) The nondetection of H$_2$CO emission implies that the emission detected in other galaxies, e.g. Arp 220, comes from high-density gas that is not directly affiliated with already-formed massive stars. Either the non-star-forming ISM of these galaxies is very dense, implying the star formation density threshold is higher, or H II regions have their emission suppressed.

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).

preprint2015arXiv

What is a GMC? Are Observers and Simulators Discussing the Same Star-forming Clouds?

As both simulations and observations reach the resolution of the star-forming molecular clouds, it becomes important to clarify if these two techniques are discussing the same objects in galaxies. We compare clouds formed in a high resolution galaxy simulation identified as continuous structures within a contour, in the simulator's position-position-position (PPP) co-ordinate space and the observer's position-position-velocity space (PPV). Results indicate that the properties of the cloud populations are similar in both methods and up to 70% of clouds have a single counterpart in the opposite data structure. Comparing individual clouds in a one-to-one match reveals a scatter in properties mostly within a factor of two. However, the small variations in mass, radius and velocity dispersion produce significant differences in derived quantities such as the virial parameter. This makes it difficult to determine if a structure is truely gravitationally bound. The three cloud types originally found in the simulation in Fujimoto et al. (2014) are identified in both data sets, with around 80% of the clouds retaining their type between identification methods. We also compared our results when using a peak decomposition method to identify clouds in both PPP and PPV space. The number of clouds increased with this technique, but the overall cloud properties remained similar. However, the more crowded environment lowered the ability to match clouds between techniques to 40%. The three cloud types also became harder to separate, especially in the PPV data set. The method used for cloud identification therefore plays a critical role in determining cloud properties, but both PPP and PPV can potentially identify the same structures.

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

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

preprint2014arXiv

Comparing Simulated Emission from Molecular Clouds Using Experimental Design

We propose a new approach to comparing simulated observations that enables us to determine the significance of the underlying physical effects. We utilize the methodology of experimental design, a subfield of statistical analysis, to establish a framework for comparing simulated position-position-velocity data cubes to each other. We propose three similarity metrics based on methods described in the literature: principal component analysis, the spectral correlation function, and the Cramer multi-variate two sample similarity statistic. Using these metrics, we intercompare a suite of mock observational data of molecular clouds generated from magnetohydrodynamic simulations with varying physical conditions. Using this framework, we show that all three metrics are sensitive to changing Mach number and temperature in the simulation sets, but cannot detect changes in magnetic field strength and initial velocity spectrum. We highlight the shortcomings of one-factor-at-a-time designs commonly used in astrophysics and propose fractional factorial designs as a means to rigorously examine the effects of changing physical properties while minimizing the investment of computational resources.

preprint2014arXiv

The Bolocam Galactic Plane Survey. XII. Distance Catalog Expansion Using Kinematic Isolation of Dense Molecular Cloud Structures With 13CO(1-0)

We present an expanded distance catalog for 1,710 molecular cloud structures identified in the Bolocam Galactic Plane Survey (BGPS) version 2, representing a nearly threefold increase over the previous BGPS distance catalog. We additionally present a new method for incorporating extant data sets into our Bayesian distance probability density function (DPDF) methodology. To augment the dense-gas tracers (e.g., HCO+(3-2), NH3(1,1)) used to derive line-of-sight velocities for kinematic distances, we utilize the Galactic Ring Survey 13CO(1-0) data to morphologically extract velocities for BGPS sources. The outline of a BGPS source is used to select a region of the GRS 13CO data, along with a reference region to subtract enveloping diffuse emission, to produce a line profile of 13CO matched to the BGPS source. For objects with a HCO+(3-2) velocity, \approx 95% of the new 13CO(1-0) velocities agree with that of the dense gas. A new prior DPDF for kinematic distance ambiguity (KDA) resolution, based on a validated formalism for associating molecular cloud structures with known objects from the literature, is presented. We demonstrate this prior using catalogs of masers with trigonometric parallaxes and HII regions with robust KDA resolutions. The distance catalog presented here contains well-constrained distance estimates for 20% of BGPS V2 sources, with typical distance uncertainties \lesssim 0.5 kpc. Approximately 75% of the well-constrained sources lie within 6 kpc of the Sun, concentrated in the Scutum-Centarus arm. Galactocentric positions of objects additionally trace out portions of the Sagittarius, Perseus, and Outer arms in the first and second Galactic quadrants, and we also find evidence for significant regions of interarm dense gas.

preprint2014arXiv

The Interstellar Medium and Star Formation in Edge-On Galaxies. II. NGC 4157, 4565, and 5907

We present a study of the vertical structure of the gaseous and stellar disks in a sample of edge-on galaxies (NGC 4157, 4565, and 5907) using BIMA/CARMA 12CO (J = 1 --> 0), VLA H I, and Spitzer 3.6 micron data. In order to take into account projection effects when we measure the disk thickness as a function of radius, we first obtain the inclination by modeling the radio data. Using the measurement of the disk thicknesses and the derived radial profiles of gas and stars, we estimate the corresponding volume densities and vertical velocity dispersions. Both stellar and gas disks have smoothly varying scale heights and velocity dispersions, contrary to assumptions of previous studies. Using the velocity dispersions, we find that the gravitational instability parameter Q follows a fairly uniform profile with radius and is greater than or equal to 1 across the star forming disk. The star formation law has a slope that is significantly different from those found in more face-on galaxy studies, both in deprojected and pixel-by-pixel plots. Midplane gas pressure based on the varying scale heights and velocity dispersions appears to roughly hold a power-law correlation with the midplane volume density ratio.

preprint2014arXiv

The JCMT Legacy Survey of the Gould Belt: a molecular line study of the Ophiuchus molecular cloud

CO, $^{13}$CO and C$^{18}$O ${\it J}$ = 3--2 observations are presented of the Ophiuchus molecular cloud. The $^{13}$CO and C$^{18}$O emission is dominated by the Oph A clump, and the Oph B1, B2, C, E, F and J regions. The optically thin(ner) C$^{18}$O line is used as a column density tracer, from which the gravitational binding energy is estimated to be $4.5 \times 10^{39}$ J (2282 $M_\odot$ km$^2$ s$^{-2}$). The turbulent kinetic energy is $6.3 \times 10^{38}$ J (320 $M_\odot$ km$^2$ s$^{-2}$), or 7 times less than this, and therefore the Oph cloud as a whole is gravitationally bound. Thirty protostars were searched for high velocity gas, with eight showing outflows, and twenty more having evidence of high velocity gas along their lines-of-sight. The total outflow kinetic energy is $1.3 \times 10^{38}$ J (67 $M_\odot$ km$^2$ s$^{-2}$), corresponding to 21$\%$ of the cloud's turbulent kinetic energy. Although turbulent injection by outflows is significant, but does ${\it not}$ appear to be the dominant source of turbulence in the cloud. 105 dense molecular clumplets were identified, which had radii $\sim$ 0.01--0.05 pc, virial masses $\sim$ 0.1--12 $M_\odot$, luminosities $\sim$ 0.001--0.1 K~km s$^{-1}$ pc$^{-2}$, and excitation temperatures $\sim$ 10--50K. These are consistent with the standard GMC based size-line width relationships, showing that the scaling laws extend down to size scales of hundredths of a parsec, and to sub solar-mass condensations. There is however no compelling evidence that the majority of clumplets are undergoing free-fall collapse, nor that they are pressure confined.

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

Hierarchical Structure of Magnetohydrodynamic Turbulence In Position-Position-Velocity Space

Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic Position-Position-Velocity (PPV) emission cubes of optically thin isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the physics of the gas, i.e. self-gravity and the global sonic and Alfvenic Mach number. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observed a strong sonic and Alfvenic dependency when we apply the the statistical moments (i.e. mean, variance, skewness, kurtosis) to the dendrogram distribution. Larger magnetic field and sonic Mach number correspond to larger values of the moments. Application of the dendrogram to 3D density cubes, also known as Position-Position-Position cubes (PPP), reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions in the direction of obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP.

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 Bolocam Galactic Plane Survey IX: Data Release 2 and Outer Galaxy Extension

We present a re-reduction and expansion of the Bolocam Galactic Plane Survey, first presented by Aguirre et al. (2011) and Rosolowsky et al. (2010). The BGPS is a 1.1 mm survey of dust emission in the Northern galactic plane, covering longitudes -10 < \ell < 90 and latitudes |b| < 0.5 with a typical 1-σRMS sensitivity of 30-100 mJy in a 33" beam. Version 2 of the survey includes an additional 20 square degrees of coverage in the 3rd and 4th quadrants and 2 square degrees in the 1st quadrant. The new data release has improved angular recovery, with complete recovery out to 80" and partial recovery to 300", and reduced negative bowls around bright sources resulting from the atmospheric subtraction process. We resolve the factor of 1.5 flux calibration offset between the v1.0 data release and other data sets and determine that there is no offset between v2.0 and other data sets. The v2.0 pointing accuracy is tested against other surveys and demonstrated to be accurate and an improvement over v1.0. We present simulations and tests of the pipeline and its properties, including measurements of the pipeline's angular transfer function. The Bolocat cataloging tool was used to extract a new catalog, which includes 8594 sources, with 591 in the expanded regions. We have demonstrated that the Bolocat 40" and 80" apertures are accurate even in the presence of strong extended background emission. The number of sources is lower than in v1.0, but the amount of flux and area included in identified sources is larger.

preprint2013arXiv

The Bolocam Galactic Plane Survey. VIII. A Mid-Infrared Kinematic Distance Discrimination Method

We present a new distance estimation method for dust-continuum-identified molecular cloud clumps. Recent (sub-)millimeter Galactic plane surveys have cataloged tens of thousands of these objects, but detailed study of their physical properties requires robust distance determinations. We derive Bayesian distance probability density functions (DPDFs) for 770 objects from the Bolocam Galactic Plane Survey in the longitude range 7.5 < l < 65. The DPDF formalism is based on kinematic distances, and uses external data sets to place prior distance probabilities to resolve the kinematic distance ambiguity (KDA) for objects in the inner Galaxy. We present priors related to the mid-infrared absorption of dust in dense molecular regions and the Galactic distribution of molecular gas. By assuming a numerical model of Galactic mid-infrared emission and simple radiative transfer, we match the morphology of millimeter thermal dust emission with mid-infrared absorption to compute a prior DPDF for distance discrimination. Selecting objects first from (sub-)millimeter source catalogs avoids a bias towards the darkest infrared dark clouds (IRDCs) and extends the range of heliocentric distance and contrast range probed by mid-infrared extinction. We derive well-constrained KDA resolutions for 618 molecular cloud clumps, with approximately 15% placed at or beyond the tangent distance. Objects sufficiently dark to be cataloged as IRDCs are generally placed at the near kinematic distance. Distance comparisons with Galactic Ring Survey KDA resolutions yield a 92% agreement, and a Galactic face-on view reveals sections of the Sagittarius and Scutum-Centaurus Arms. This KDA-resolution method for large catalogs of sources through the combination of (sub-)millimeter and mid-infrared observations of molecular cloud clumps is generally applicable to other dust-continuum Galactic plane surveys.

preprint2013arXiv

The Bolocam Galactic Plane Survey. X. A Complete Spectroscopic Catalog of Dense Molecular Gas Observed toward 1.1 mm Dust Continuum Sources with 7.5 <= l <= 194 degrees

The Bolocam Galactic Plane Survey (BGPS) is a 1.1 mm continuum survey of dense clumps of dust throughout the Galaxy covering 170 square degrees. We present spectroscopic observations using the Heinrich Hertz Submillimeter Telescope of the dense gas tracers, HCO+ and N2H+ 3-2, for all 6194 sources in the Bolocam Galactic Plane Survey v1.0.1 catalog between 7.5 <= l <= 194 degrees. This is the largest targeted spectroscopic survey of dense molecular gas in the Milky Way to date. We find unique velocities for 3126 (50.5%) of the BGPS v1.0.1 sources observed. Strong N2H+ 3-2 emission (T_{mb} > 0.5 K) without HCO+ 3-2 emission does not occur in this catalog. We characterize the properties of the dense molecular gas emission toward the entire sample. HCO+ is very sub-thermally populated and the 3-2 transitions are optically thick toward most BGPS clumps. The median observed line width is 3.3 km/s consistent with supersonic turbulence within BGPS clumps. We find strong correlations between dense molecular gas integrated intensities and 1.1 mm peak flux and the gas kinetic temperature derived from previously published NH3 observations. These intensity correlations are driven by the sensitivity of the 3-2 transitions to excitation conditions rather than by variations in molecular column density or abundance. We identify a subset of 113 sources with stronger N2H+ than HCO+ integrated intensity, but we find no correlations between the N2H+ / HCO+ ratio and 1.1 mm continuum flux density, gas kinetic temperature, or line width. Self-absorbed profiles are rare (1.3%).

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

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.

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

Minimal HCN emission from Molecular Clouds in M33

Since HCN emission has been shown to be a linear tracer of ongoing star formation activity, we have searched for HCN (J = 1->0) emission from known GMCs in the nearby galaxy M33. No significant HCN emission has been found along any of the lines of sight. We find two lines of sight where CO-to-HCN integrated intensity ratios up to 280, nearly a factor of 6 above what is found in comparable regions of the Milky Way. Star formation tracers suggest that the HCN-to-star formation rate ratio (L_HCN/M_SFR) is a factor of six lower than what is observed in the Milky Way (on average) and local extragalactic systems. Simple chemical models accounting for the sub-solar N/O ratio suggest that depletion cannot account for the high CO-to-HCN ratios. Given HCN formation requires high extinction (A_V > 4), low metallicity may yield reduced dust shielding and thus a high CO/HCN ratio. The turbulence and structure of GMCs in M33 are comparable to those found in other systems, so the differences are unlikely to result from different GMC properties. Since lower CO-to-HCN ratios are associated with the highest rates of star formation, we attribute the deficits in part to evolutionary effects within GMCs.

preprint2011arXiv

The Bolocam Galactic Plane Survey V: HCO+ and N2H+ Spectroscopy of 1.1 mm Dust Continuum Sources

We present the results of observations of 1882 sources in the Bolocam Galactic Plane Survey (BGPS) at 1.1 mm with the 10m Heinrich Hertz Telescope simultaneously in HCO+ J=3-2 and N2H+ J=3-2. We detect 77% of these sources in HCO^+ and 51% in N2H+ at greater than 3$σ$. We find a strong correlation between the integrated intensity of both dense gas tracers and the 1.1 mm dust emission of BGPS sources. We determine kinematic distances for 529 sources (440 in the first quadrant breaking the distance ambiguity and 89 in the second quadrant) We derive the size, mass, and average density for this subset of clumps. The median size of BGPS clumps is 0.75 pc with a median mass of 330 M$_{\odot}$ (assuming T_{Dust}=20 K). The median HCO+ linewidth is 2.9 km s$^{-1}$ indicating that BGPS clumps are dominated by supersonic turbulence or unresolved kinematic motions. We find no evidence for a size-linewidth relationship for BGPS clumps. We analyze the effects of the assumed dust temperature on the derived clump properties with a Monte Carlo simulation and we find that changing the temperature distribution will change the median source properties (mass, volume-averaged number density, surface density) by factors of a few. The observed differential mass distribution has a power-law slope that is intermediate between that observed for diffuse CO clouds and the stellar IMF. BGPS clumps represent a wide range of objects (from dense cores to more diffuse clumps) and are typically characterized by larger sizes and lower densities than previously published surveys of high-mass star-forming regions. This collection of objects is a less-biased sample of star-forming regions in the Milky Way that likely span a wide range of evolutionary states.

preprint2011arXiv

The Bolocam Galactic Plane Survey. VII. Characterizing the Properties of Massive Star-Forming Regions

We present the results of a GBT survey of NH3(1,1), (2,2), (3,3) lines towards 631 Bolocam Galactic Plane Survey (BGPS) sources at a range of Galactic longitudes in the inner Galaxy. We have detected the NH3(1,1) line towards 72% of our targets (456), demonstrating that the high column density features identified in the BGPS and other continuum surveys accurately predict the presence of dense gas. We have determined kinematic distances and resolved the distance ambiguity for all BGPS sources detected in NH3. The BGPS sources trace the locations of the Scutum and Sagittarius spiral arms, with the number of sources peaking between Galactocentric radii of 4-5 kpc. We measure the physical properties of each source and find that depending on the distance, BGPS sources are primarily clumps, with some cores and clouds. We have examined the physical properties as a function of Galactocentric distance, and find a mean gas kinetic temperature of 15.6 K, and that the NH3 column density and abundance decrease by nearly an order of magnitude between Galactocentric radii of 3 - 11 kpc. Comparing sources at similar distances demonstrates that the physical properties are indistinguishable, which suggests a similarity in clump structure across the Galactic disk. We have also compared the BGPS sources to criteria for efficient star formation presented independently by Heiderman et al. and Lada et al., and massive star formation presented by Kauffmann et al. 48% of our sample should be forming stars (including massive stars) with high efficiency, and 87% contain subregions that should be efficiently forming stars. Indeed, we find that 67% of the sample exhibit signs of star formation activity based on an association with a mid-IR source.

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.

preprint2010arXiv

Dense Gas Tracers in Perseus: Relating the N2H+, NH3, and Dust Continuum Properties of Pre- and Proto-Stellar Cores

We investigate 35 pre-stellar cores and 36 proto-stellar cores in the Perseus molecular cloud. We find a very tight correlation between the physical parameters describing the N2H+ and NH3 gas. Both the velocity centroids and the line widths of N2H+ and NH3 correlate much better than either species correlates with CO, as expected if the nitrogen-bearing species are probing primarily the dense core gas where the CO has been depleted. We also find a tight correlation in the inferred abundance ratio between N2H+ and para-NH3 across all cores, with N(p-NH3)/N(N2H+)= 22 +/- 10. We find a mild correlation between NH3 (and N2H+) column density and the (sub)millimeter dust continuum derived H2 column density for pre-stellar cores, N(p-NH3)/N(H2) ~ 10e-8, but do not find a fixed ratio for proto-stellar cores. The observations suggest that in the Perseus molecular cloud the formation and destruction mechanisms for the two nitrogen-bearing species are similar, regardless of the physical conditions in the dense core gas. While the equivalence of N2H+ and NH3 as powerful tracers of dense gas is validated, the lack of correspondence between these species and the (sub)millimeter dust continuum observations for proto-stellar cores is disconcerting and presently unexplained.

preprint2010arXiv

Planetary nebulae in M33: probes of AGB nucleosynthesis and ISM abundances

We have obtained deep optical spectrophotometry of 16 planetary nebulae in M33, mostly located in the central two kpc of the galaxy, with the Subaru and Keck telescopes. We have derived electron temperatures and chemical abundances from the detection of the [OIII]4363 line for the whole sample. We have found one object with an extreme nitrogen abundance, 12+log(N/H)=9.20, accompanied by a large helium content. After combining our data with those available in the literature for PNe and HII regions, we have examined the behavior of nitrogen, neon, oxygen and argon in relation to each other, and as a function of galactocentric distance. We confirm the good correlation between Ne/H and O/H for PNe in M33. Ar/H is also found to correlate with O/H. This strengthens the idea that at the metallicity of the bright PNe analyzed in M33, which is similar to that found in the LMC, these elements have not been significantly modified during the dredge-up processes that take place during the AGB phase of their progenitor stars. We find no significant oxygen abundance offset between PNe and HII regions at any given galactocentric distance, despite the fact that these objects represent different age groups in the evolution of the galaxy. Combining the results from PNe and HII regions, we obtain a representative slope of the ISM alpha-element (O, Ar, Ne) abundance gradient in M33 of -0.025 +/- 0.006 dex/kpc. Both PNe and HII regions display a large abundance dispersion at any given distance from the galactic center. We find that the N/O ratio in PNe is enhanced, relative to the HII regions, by approximately 0.8 dex.

preprint2010arXiv

The Bolocam Galactic Plane Survey -- III. Characterizing Physical Properties of Massive Star-Forming Regions in the Gemini OB1 Molecular Cloud

We present the 1.1 millimeter Bolocam Galactic Plane Survey (BGPS) observations of the Gemini OB1 molecular cloud complex, and targeted ammonia observations of the BGPS sources. When paired with molecular spectroscopy of a dense gas tracer, millimeter observations yield physical properties such as masses, radii, mean densities, kinetic temperatures and line widths. We detect 34 distinct BGPS sources above 5-sigma=0.37 Jy/beam with corresponding 5-sigma detections in the ammonia (1,1) transition. Eight of the objects show water maser emission (20%). We find a mean millimeter source FWHM of 1.12 pc, and a mean kinetic temperature of 20 K for the sample of 34 BGPS sources. The observed ammonia line widths are dominated by non-thermal motions, typically found to be a few times the thermal sound speed expected for the derived kinetic temperature. We calculate the mass for each source from the millimeter flux assuming the sources are isothermal and find a mean isothermal mass within a 120" aperture of 230 +/- 180 solar masses. We find a total mass of 8,400 solar masses for all BGPS sources in the Gemini OB1 molecular cloud, representing 6.5% of the cloud mass. By comparing the millimeter isothermal mass to the virial mass within a radius equal to the mm source size calculated from the ammonia line widths, we find a mean virial parameter (M_vir/M_iso) of 1.0 +/- 0.9 for the sample. We find mean values for the distributions of column densities of 10^22 cm^-2 for H_2, and 3.0x10^14 cm^-2 for ammonia, giving a mean ammonia abundance of 3.0x10^-8 relative to H_2. We find volume-averaged densities on the order of 10^3-10^4 cm^-3. The sizes and densities suggest that in the Gem OB1 region the BGPS is detecting the clumps from which stellar clusters form, rather than smaller, higher density cores where single stars or small multiple systems form.

preprint2010arXiv

The Bolocam Galactic Plane Survey IV: 1.1 and 0.35 mm Dust Continuum Emission in the Galactic Center Region

The Bolocam Galactic Plane Survey (BGPS) data for a six square degree region of the Galactic plane containing the Galactic center is analyzed and compared to infrared and radio continuum data. The BGPS 1.1 mm emission consists of clumps interconnected by a network of fainter filaments surrounding cavities, a few of which are filled with diffuse near-IR emission indicating the presence of warm dust or with radio continuum characteristic of HII regions or supernova remnants. New 350 μm images of the environments of the two brightest regions, Sgr A and B, are presented. Sgr B2 is the brightest mm-emitting clump in the Central Molecular Zone and may be forming the closest analog to a super star cluster in the Galaxy. The Central Molecular Zone (CMZ) contains the highest concentration of mm and sub-mm emitting dense clumps in the Galaxy. Most 1.1 mm features at positive longitudes are seen in silhouette against the 3.6 to 24 μm background observed by the Spitzer Space Telescope. However, only a few clumps at negative longitudes are seen in absorption, confirming the hypothesis that positive longitude clumps in the CMZ tend to be on the near-side of the Galactic center, consistent with the suspected orientation of the central bar in our Galaxy. Some 1.1 mm cloud surfaces are seen in emission at 8 μm, presumably due to polycyclic aromatic hydrocarbons (PAHs). A ~0.2\degree (~30 pc) diameter cavity and infrared bubble between l \approx 0.0\degree and 0.2\degree surrounds the Arches and Quintuplet clusters and Sgr A. The bubble contains several clumpy dust filaments that point toward Sgr A\ast; its potential role in their formation is explored. [abstract truncated]

preprint2010arXiv

The Bolocam Galactic Plane Survey: Survey Description and Data Reduction

We present the Bolocam Galactic Plane Survey (BGPS), a 1.1 mm continuum survey at 33" effective resolution of 170 square degrees of the Galactic Plane visible from the northern hemisphere. The survey is contiguous over the range -10.5 < l < 90.5, |b| < 0.5 and encompasses 133 square degrees, including some extended regions |b| < 1.5. In addition to the contiguous region, four targeted regions in the outer Galaxy were observed: IC1396, a region towards the Perseus Arm, W3/4/5, and Gem OB1. The BGPS has detected approximately 8400 clumps over the entire area to a limiting non-uniform 1-sigma noise level in the range 11 to 53 mJy/beam in the inner Galaxy. The BGPS source catalog is presented in a companion paper (Rosolowsky et al. 2010). This paper details the survey observations and data reduction methods for the images. We discuss in detail the determination of astrometric and flux density calibration uncertainties and compare our results to the literature. Data processing algorithms that separate astronomical signals from time-variable atmospheric fluctuations in the data time-stream are presented. These algorithms reproduce the structure of the astronomical sky over a limited range of angular scales and produce artifacts in the vicinity of bright sources. Based on simulations, we find that extended emission on scales larger than about 5.9' is nearly completely attenuated (> 90%) and the linear scale at which the attenuation reaches 50% is 3.8'. Comparison with other millimeter-wave data sets implies a possible systematic offset in flux calibration, for which no cause has been discovered. This presentation serves as a companion and guide to the public data release through NASA's Infrared Processing and Analysis Center (IPAC) Infrared Science Archive (IRSA). New data releases will be provided through IPAC IRSA with any future improvements in the reduction.

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).

preprint2009arXiv

Ionized Gas Towards Molecular Clumps: Physical Properties of Massive Star Forming Regions

We have conducted a search for ionized gas at 3.6 cm, using the Very Large Array, towards 31 Galactic intermediate- and high-mass clumps detected in previous millimeter continuum observations. In the 10 observed fields, 35 HII regions are identified, of which 20 are newly discovered. Many of the HII regions are multiply peaked indicating the presence of a cluster of massive stars. We find that the ionized gas tends to be associated towards the millimeter clumps; of the 31 millimeter clumps observed, 9 of these appear to be physically related to ionized gas, and a further 6 have ionized gas emission within 1'. For clumps with associated ionized gas, the combined mass of the ionizing massive stars is compared to the clump masses to provide an estimate of the instantaneous star formation efficiency. These values range from a few percent to 25%, and have an average of 7 +/- 8%. We also find a correlation between the clump mass and the mass of the ionizing massive stars within it, which is consistent with a power law. This result is comparable to the prediction of star formation by competitive accretion that a power law relationship exists between the mass of the most massive star in a cluster and the total mass of the remaining stars.

preprint2009arXiv

The Bolocam Galactic Plane Survey -- II. Catalog of the Image Data

We present a catalog of 8358 sources extracted from images produced by the Bolocam Galactic Plane Survey (BGPS). The BGPS is a survey of the millimeter dust continuum emission from the northern Galactic plane. The catalog sources are extracted using a custom algorithm, Bolocat, which was designed specifically to identify and characterize objects in the large-area maps generated from the Bolocam instrument. The catalog products are designed to facilitate follow-up observations of these relatively unstudied objects. The catalog is 98% complete from 0.4 Jy to 60 Jy over all object sizes for which the survey is sensitive (<3.5'). We find that the sources extracted can best be described as molecular clumps -- large dense regions in molecular clouds linked to cluster formation. We find the flux density distribution of sources follows a power law with dN/dS ~S^(-2.4 +/- 0.1) and that the mean Galactic latitude for sources is significantly below the midplane: <b>=(-0.095 +/- 0.001) deg.

preprint2008arXiv

The Gas Temperature of Starless Cores in Perseus

In this paper we study the determinants of starless core temperatures in the Perseus molecular cloud. We use NH3 (1,1) and (2,2) observations to derive core temperatures (T_kin) and data from the COMPLETE Survey of Star Forming Regions and the c2d Spitzer Legacy Survey for observations of the other core and molecular cloud properties. The kinetic temperature distribution probed by NH3 is in the fairly narrow range of 9 - 15 K. We find that cores within the clusters IC348 and NGC1333 are significantly warmer than "field" starless cores, and T_kin is higher within regions of larger extinction-derived column density. Starless cores in the field are warmer when they are closer to class O/I protostars, but this effect is not seen for those cores in clusters. For field starless cores, T_kin is higher in regions in which the 13CO linewidth and the 1.1mm flux from the core are larger, and T_kin is lower when the the peak column density within the core and average volume density of the core are larger. There is no correlation between T_kin and 13CO linewidth, 1.1mm flux, density or peak column density for those cores in clusters. The temperature of the cloud material along the line of sight to the core, as measured by CO or far-infrared emission from dust, is positively correlated with core temperature when considering the collection of cores in the field and in clusters, but this effect is not apparent when the two subsamples of cores are considered separately.

Erik Rosolowsky

What is connected

Connect this record

See the researcher in context

Building this map preview

85 published item(s)

12CO (3-2) High-Resolution Survey (COHRS) of the Galactic Plane: Complete Data Release

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)

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

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

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

pyspeckit: A spectroscopic analysis and plotting package

Stellar feedback in M83 as observed with MUSE -- I. Overview, an unprecedented view of the stellar and gas kinematics and evidence of outflowing gas

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

The Panchromatic Hubble Andromeda Treasury: Triangulum Extended Region (PHATTER) II. The Spatially Resolved Recent Star Formation History of M33

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

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

[CI](1-0) and [CI](2-1) in resolved local galaxies

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

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

Efficient methanol production on the dark side of a prestellar core

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

Relative Alignment between Dense Molecular Cores and Ambient Magnetic Field: The Synergy of Numerical Models and Observations

The MUSTANG-2 Galactic Plane Survey (MGPS90) pilot

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

Ubiquitous $\rm NH_3$ supersonic component in L1688 coherent cores

Velocity-coherent Filaments in NGC 1333: Evidence for Accretion Flow?

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

Canada and the SKA from 2020-2030

Canadian Investigations of the Interstellar Medium

CLOVER: Convnet Line-fitting Of Velocities in Emission-line Regions

Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

LRP2020: The cosmic origin and evolution of the elements

The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks

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

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

The Next Generation Very Large Array

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

Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations

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

Effects of Galactic Disc Inclination and Resolution on Observed GMC Properties and Larson's Scaling Relations

The inner mass distribution of late-type spiral galaxies from SAURON stellar kinematic maps

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

Giant Molecular Clouds in the Early-Type Galaxy NGC4526

Graph-based interpretation of the Molecular Interstellar Medium Segmentation

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

The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations

The Bolocam Galactic Plane Survey. XIII. Physical Properties and Mass Functions of Dense Molecular Cloud Structures

The dense gas mass fraction in the W51 cloud and its protoclusters

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

What is a GMC? Are Observers and Simulators Discussing the Same Star-forming Clouds?

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main

Comparing Simulated Emission from Molecular Clouds Using Experimental Design

The Bolocam Galactic Plane Survey. XII. Distance Catalog Expansion Using Kinematic Isolation of Dense Molecular Cloud Structures With 13CO(1-0)

The Interstellar Medium and Star Formation in Edge-On Galaxies. II. NGC 4157, 4565, and 5907

The JCMT Legacy Survey of the Gould Belt: a molecular line study of the Ophiuchus molecular cloud

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

Hierarchical Structure of Magnetohydrodynamic Turbulence In Position-Position-Velocity Space

Molecular Gas and Star Formation in Nearby Disk Galaxies

The Bolocam Galactic Plane Survey IX: Data Release 2 and Outer Galaxy Extension

The Bolocam Galactic Plane Survey. VIII. A Mid-Infrared Kinematic Distance Discrimination Method

The Bolocam Galactic Plane Survey. X. A Complete Spectroscopic Catalog of Dense Molecular Gas Observed toward 1.1 mm Dust Continuum Sources with 7.5 <= l <= 194 degrees

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

Low CO Luminosities in Dwarf Galaxies

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

Minimal HCN emission from Molecular Clouds in M33

The Bolocam Galactic Plane Survey V: HCO+ and N2H+ Spectroscopy of 1.1 mm Dust Continuum Sources