sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Rosolowsky, E.

Normalized to: Rosolowsky, E.

175 article(s) in total. 830 co-authors, from 1 to 48 common article(s). Median position in authors list is 5,0.

[1]  [pdf] - 2042415
The X-ray Emissivity of Low-Density Stellar Populations
Comments: 28 pages, 23 figures. MNRAS, in press
Submitted: 2020-01-26
The dynamical production of low-mass X-ray binaries and brighter cataclysmic variables (CVs) in dense globular clusters is well-established. We investigate how the X-ray emissivity of fainter X-ray binaries (principally CVs and coronally active binaries) varies between different environments. We compile calculations (largely from the literature) of the X-ray emissivity of old stellar populations, including open and globular clusters and several galaxies. We investigate three literature claims of unusual X-ray sources in low-density stellar populations. We show that a suggested quiescent neutron star in the open cluster NGC 6819 is a foreground M dwarf. We show that the suggested diffuse X-ray emission from an old nova shell in the globular cluster NGC 6366 is actually a background galaxy cluster. And we show that a suggested population of quiescent X-ray binaries in the Sculptor Dwarf Galaxy is mostly (perhaps entirely) background galaxies. We find that above densities of $10^4$ M$_{\odot}$/pc$^3$, the X-ray emissivity of globular clusters increases, due to dynamical production of X-ray emitting systems. Below this density, globular clusters have lower X-ray emissivity than the other populations, and we do not see a strong dependence of X-ray emissivity due to density effects. We find significant correlations between X-ray emissivity and binary fraction, metallicity, and density. Sampling these fits via bootstrap techniques gives less significant correlations, but confirms the effect of metallicity on low-density populations, and that of density on the full globular cluster sample.
[2]  [pdf] - 2034714
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
Comments: Accepted for publication in ApJ, 32 pages, 11 figures
Submitted: 2020-01-21
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.
[3]  [pdf] - 2032037
[CI](1-0) and [CI](2-1) in resolved local galaxies
Comments: 15 pages, 10 figures, published in ApJ
Submitted: 2020-01-14
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 $\alpha_{\mathrm{[CI](1-0)}} = 7.3$ M$_{\mathrm{sun}}$ K$^{-1}$ km$^{-1}$ s pc$^{-2}$ and $\alpha_{\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 $\alpha_{\mathrm{[CI](1-0)}}$ to $\alpha_{\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 $\alpha_{\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.
[4]  [pdf] - 2029759
Spatial Power Spectra of Dust across the Local Group: No Constraint on Disc Scale Height
Comments: Accepted by MNRAS; 19 pgs, 10 figures. Supplementary figures of the power spectrum fits for each band can be viewed here:
Submitted: 2019-12-19, last modified: 2020-01-10
We analyze the 1D spatial power spectra of dust surface density and mid to far-infrared emission at $24-500\,\mu$m in the LMC, SMC, M31, and M33. By forward-modelling the point-spread-function (PSF) on the power spectrum, we find that nearly all power spectra have a single power-law and point source component. A broken power-law model is only favoured for the LMC 24 $\mu$m MIPS power spectrum and is due to intense dust heating in 30 Doradus. We also test for local power spectrum variations by splitting the LMC and SMC maps into $820$ pc boxes. We find significant variations in the power-law index with no strong evidence for breaks. The lack of a ubiquitous break suggests that the spatial power spectrum does not constrain the disc scale height. This contradicts claims of a break where the turbulent motion changes from 3D to 2D. The power spectrum indices in the LMC, SMC, and M31 are similar ($2.0-2.5$). M33 has a flatter power spectrum ($1.3$), similar to more distant spiral galaxies with a centrally-concentrated H$_{2}$ distribution. We compare the power spectra of HI, CO, and dust in M31 and M33 and find that HI power spectra are consistently flatter than CO power spectra. These results cast doubt on the idea that the spatial power spectrum traces large scale turbulent motion in nearby galaxies. Instead, we find that the spatial power spectrum is influenced by (1) the PSF on scales below $\sim3$ times the FWHM, (2) bright compact regions (30 Doradus), and (3) the global morphology of the tracer (an exponential CO disc).
[5]  [pdf] - 2042157
The Karl G. Jansky Very Large Array Sky Survey (VLASS). Science case and survey design
Comments: 34 pages, accepted by PASP (modified from prior version to address referee's and coauthor comments). (v2) Minor fixes to author list
Submitted: 2019-07-03, last modified: 2019-12-30
The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.
[6]  [pdf] - 2020603
Rotation of molecular clouds in M~51
Comments: Accepted by Astronomy and Astrophysics
Submitted: 2019-11-20
The grand-design spiral galaxy M~51 was observed at 40pc resolution in CO(1--0) by the PAWS project. A large number of molecular clouds were identified and we search for velocity gradients in two high signal-to-noise subsamples, containing 682 and 376 clouds. The velocity gradients are found to be systematically prograde oriented, as was previously found for the rather flocculent spiral M~33. This strongly supports the idea that the velocity gradients reflect cloud rotation, rather than more random dynamical forces, such as turbulence. Not only are the gradients prograde, but their $\frac{\partial v}{\partial x}$ and $\frac{\partial v}{\partial y}$ coefficients follow galactic shear in sign, although with a lower amplitude. No link is found between the orientation of the gradient and the orientation of the cloud. The values of the cloud angular momenta appear to be an extension of the values noted for galactic clouds despite the orders of magnitude difference in cloud mass. Roughly 30\% of the clouds show retrograde velocity gradients. For a strictly rising rotation curve, as in M~51, gravitational contraction would be expected to yield strictly prograde rotators within an axisymmetric potential. In M~51, the fraction of retrograde rotators is found to be higher in the spiral arms than in the disk as a whole. Along the leading edge of the spiral arms, a majority of the clouds are retrograde rotators. While this work should be continued on other nearby galaxies, the M~33 and M~51 studies have shown that clouds rotate and that they rotate mostly prograde, although the amplitudes are not such that rotational energy is a significant support mechanism against gravitation. In this work, we show that retrograde rotation is linked to the presence of a spiral gravitational potential.
[7]  [pdf] - 1996168
The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks
Comments: 11 pages, a contributed white paper prepared for Canada's 2020 Long Range Plan decadal process
Submitted: 2019-11-05, last modified: 2019-11-09
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.
[8]  [pdf] - 2026285
The lifecycle of molecular clouds in nearby star-forming disc galaxies
Comments: 39 pages, 14 figures; resubmitted to MNRAS after a favourable referee report (November 6, 2019)
Submitted: 2019-11-08
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$\alpha$ 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 $\Sigma_{\rm H_2}\geqslant8$M$_{\odot}$pc$^{-2}$, the GMC lifetime correlates with time-scales for galactic dynamical processes, whereas at $\Sigma_{\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$\alpha$ (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.
[9]  [pdf] - 1995601
Canada and the SKA from 2020-2030
Comments: 14 pages, 4 figures, 2020 Canadian Long-Range Plan (LRP2020) white paper
Submitted: 2019-11-08
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.
[10]  [pdf] - 1994371
The Next Generation Very Large Array
Comments: 11 pages; a contributed white paper for Canada's 2020 Long Range Plan decadal process
Submitted: 2019-11-04
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.
[11]  [pdf] - 1990407
The headlight cloud in NGC 628: An extreme giant molecular cloud in a typical galaxy disk
Comments: A&A accepted
Submitted: 2019-10-31, last modified: 2019-11-04
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.
[12]  [pdf] - 2026110
The Gas-Star Formation Cycle in Nearby Star-Forming Galaxies I. Assessment of Multi-scale Variations
Comments: Accepted to ApJ, 29 pages with 9 figures and 7 tables
Submitted: 2019-10-23
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.
[13]  [pdf] - 1984055
LRP2020: The cosmic origin and evolution of the elements
Comments: White paper submitted to the Canadian Long Range Plan 2020. Minor formatting changes relative to submitted version
Submitted: 2019-10-21
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.
[14]  [pdf] - 2025690
Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures
Comments: accepted to ApJ in Oct, 2019; in press
Submitted: 2019-08-12, last modified: 2019-10-18
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, $\beta$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $\beta \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.
[15]  [pdf] - 2026025
Mapping metallicity variations across nearby galaxy disks
Comments: 14 pages, 9 figures, 4 Appendices including image atlases, ApJ accepted
Submitted: 2019-10-16
The distribution of metals within a galaxy traces the baryon cycle and the buildup of galactic disks, but the detailed gas phase metallicity distribution remains poorly sampled. We have determined the gas phase oxygen abundances for 7,138 HII regions across the disks of eight nearby galaxies using VLT/MUSE optical integral field spectroscopy as part of the PHANGS-MUSE survey. After removing the first order radial gradients present in each galaxy, we look at the statistics of the metallicity offset (Delta O/H) and explore azimuthal variations. Across each galaxy, we find low (sigma=0.03-0.05 dex) scatter at any given radius, indicative of efficient mixing. We compare physical parameters for those HII regions that are 1 sigma outliers towards both enhanced and reduced abundances. Regions with enhanced abundances have high ionization parameter, higher Halpha luminosity, lower Halpha velocity dispersion, younger star clusters and associated molecular gas clouds show higher molecular gas densities. This indicates recent star formation has locally enriched the material. Regions with reduced abundances show increased Halpha velocity dispersions, suggestive of mixing introducing more pristine material. We observe subtle azimuthal variations in half of the sample, but can not always cleanly associate this with the spiral pattern. Regions with enhanced and reduced abundances are found distributed throughout the disk, and in half of our galaxies we can identify subsections of spiral arms with clearly associated metallicity gradients. This suggests spiral arms play a role in organizing and mixing the ISM.
[16]  [pdf] - 1977460
Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)
Comments: White paper E004 submitted to the Canadian Long Range Plan 2020
Submitted: 2019-10-09
(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.
[17]  [pdf] - 1976903
Canadian Investigations of the Interstellar Medium
Comments: White paper for the Canadian astronomy Long Range Plan process. 10 pages.
Submitted: 2019-10-08
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.
[18]  [pdf] - 2025840
CLOVER: Convnet Line-fitting Of Velocities in Emission-line Regions
Comments: Accepted for publication in ApJ
Submitted: 2019-09-18
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.
[19]  [pdf] - 1956586
The CCAT-Prime Submillimeter Observatory
Comments: Astro2020 APC White Paper
Submitted: 2019-09-05
The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.
[20]  [pdf] - 1979698
KFPA Examinations of Young STellar Object Natal Environments (KEYSTONE): Hierarchical Ammonia Structures in Galactic Giant Molecular Clouds
Comments: Accepted for publication in ApJ
Submitted: 2019-08-27, last modified: 2019-08-29
We present initial results from the K-band focal plane array Examinations of Young STellar Object Natal Environments (KEYSTONE) survey, a large project on the 100-m Green Bank Telescope mapping ammonia emission across eleven giant molecular clouds at distances of $0.9-3.0$ kpc (Cygnus X North, Cygnus X South, M16, M17, MonR1, MonR2, NGC2264, NGC7538, Rosette, W3, and W48). This data release includes the NH$_3$ (1,1) and (2,2) maps for each cloud, which are modeled to produce maps of kinetic temperature, centroid velocity, velocity dispersion, and ammonia column density. Median cloud kinetic temperatures range from $11.4\pm2.2$ K in the coldest cloud (MonR1) to $23.0\pm6.5$ K in the warmest cloud (M17). Using dendrograms on the NH$_3$ (1,1) integrated intensity maps, we identify 856 dense gas clumps across the eleven clouds. Depending on the cloud observed, $40-100\%$ of the clumps are aligned spatially with filaments identified in H$_2$ column density maps derived from SED-fitting of dust continuum emission. A virial analysis reveals that 523 of the 835 clumps ($\sim63\%$) with mass estimates are bound by gravity alone. We find no significant difference between the virial parameter distributions for clumps aligned with the dust-continuum filaments and those unaligned with filaments. In some clouds, however, hubs or ridges of dense gas with unusually high mass and low virial parameters are located within a single filament or at the intersection of multiple filaments. These hubs and ridges tend to host water maser emission, multiple 70$\mu$m-detected protostars, and have masses and radii above an empirical threshold for forming massive stars.
[21]  [pdf] - 1953599
Northern Galactic Molecular Cloud Clumps in Hi-GAL: Clump and Star Formation within Clouds
Comments: Accepted to ApJ
Submitted: 2019-06-29
We investigate how the properties of Galactic giant molecular clouds (GMCs) and their denser substructures (clumps) correlate with the local star formation rate. We trace clouds using the $^{12}$CO(3-2) transition, as observed by the CO High Resolution Survey (COHRS). We identify their constituent clumps using thermal dust emission, as observed by the Herschel infrared GALactic plane survey (Hi-GAL). We estimate star formation rates in these clouds using 70 $\mu$m emission. In total, we match 3,674 clumps to 473 clouds in position-position-velocity space spanning the Galactic longitude range $10^\circ<\ell<56^\circ$. We find that more massive clouds produce more clumps and more massive clumps. These clumps have average number densities an order of magnitude greater than their host clouds. We find a mean clump mass fraction of $0.20^{+0.13}_{-0.10}$. This mass fraction weakly varies with mass and mass surface density of clouds, and shows no clear dependence on the virial parameter and line width of the clouds. The average clump mass fraction is only weakly dependent upon Galactocentric radius. Although the scatter in our measured properties is significant, the star formation rate for clouds is independent of clump mass fraction. However, there is a positive correlation between the depletion times for clouds and clump mass fraction. We find a star formation efficiency per free fall time of $\epsilon_{\mathrm{ff}}=0.15\%$ for GMCs but $\epsilon_{\mathrm{ff}}=0.37\%$ for clumps.
[22]  [pdf] - 1920084
The EDGE-CALIFA Survey: Exploring the Star Formation Law through Variable Selection
Comments: Accepted by MNRAS
Submitted: 2019-06-05, last modified: 2019-06-25
We present a multilinear analysis to determine the significant predictors of star formation in galaxies using the combined EDGE-CALIFA sample of galaxies. We analyze 1845 kpc-scale lines of sight across 39 galaxies with molecular line emission measurements from EDGE combined with optical IFU data drawn from CALIFA. We use the Least Absolute Shrinkage and Selection Operator (LASSO) to identify significant factors in predicting star formation rates. We find that the local star formation rate surface density is increased by higher molecular gas surface densities and stellar surface densities. In contrast, we see lower star formation rates in systems with older stellar populations, higher gas- and stellar-phase metallicities and larger galaxy masses. We also find a significant increase in star formation rate with galactocentric radius normalized by the disk scale length, which suggests additional parameters regulating star formation rate not explored in this study.
[23]  [pdf] - 1938308
Droplets I: Pressure-Dominated Sub-0.1 pc Coherent Structures in L1688 and B18
Comments: Accepted by ApJ in April, 2019
Submitted: 2018-09-26, last modified: 2019-05-15
We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.
[24]  [pdf] - 1898029
TurbuStat: Turbulence Statistics in Python
Comments: Accepted in AJ. 21 pages, 8 figures
Submitted: 2019-04-23, last modified: 2019-04-25
We present TurbuStat (v1.0): a Python package for computing turbulence statistics in spectral-line data cubes. TurbuStat includes implementations of fourteen methods for recovering turbulent properties from observational data. Additional features of the software include: distance metrics for comparing two data sets; a segmented linear model for fitting lines with a break-point; a two-dimensional elliptical power-law model; multi-core fast-fourier-transform support; a suite for producing simulated observations of fractional Brownian Motion fields, including two-dimensional images and optically-thin HI data cubes; and functions for creating realistic world coordinate system information for synthetic observations. This paper summarizes the TurbuStat package and provides representative examples using several different methods. TurbuStat is an open-source package and we welcome community feedback and contributions.
[25]  [pdf] - 1875340
Dense gas is not enough: environmental variations in the star formation efficiency of dense molecular gas at 100pc scales in M51
Comments: 23 pages, 9 figures, accepted for publication in A&A
Submitted: 2019-02-27, last modified: 2019-04-09
It remains unclear what sets the efficiency with which molecular gas transforms into stars. Here we present a new VLA map of the spiral galaxy M51 in 33GHz radio continuum, an extinction-free tracer of star formation, at 3" scales (~100pc). We combined this map with interferometric PdBI/NOEMA observations of CO(1-0) and HCN(1-0) at matched resolution for three regions in M51 (central molecular ring, northern and southern spiral arm segments). While our measurements roughly fall on the well-known correlation between total infrared and HCN luminosity, bridging the gap between Galactic and extragalactic observations, we find systematic offsets from that relation for different dynamical environments probed in M51, e.g. the southern arm segment is more quiescent due to low star formation efficiency (SFE) of the dense gas, despite having a high dense gas fraction. Combining our results with measurements from the literature at 100pc scales, we find that the SFE of the dense gas and the dense gas fraction anti-correlate and correlate, respectively, with the local stellar mass surface density. This is consistent with previous kpc-scale studies. In addition, we find a significant anti-correlation between the SFE and velocity dispersion of the dense gas. Finally, we confirm that a correlation also holds between star formation rate surface density and the dense gas fraction, but it is not stronger than the correlation with dense gas surface density. Our results are hard to reconcile with models relying on a universal gas density threshold for star formation and suggest that turbulence and galactic dynamics play a major role in setting how efficiently dense gas converts into stars.
[26]  [pdf] - 1861394
Astro2020 Science White Paper: Physical Conditions in the Cold Gas of Local Galaxies
Comments: Submitted in response to the call for Astro2020 Science White Papers, 8 pages, 2 figures
Submitted: 2019-04-02
We describe a next major frontier in observational studies of galaxy evolution and star formation: linking the physical conditions in the cold, star-forming interstellar medium to host galaxy and local environment. The integrated gas content of galaxies has been surveyed extensively over the last decades. The link between environment and cold gas density, turbulence, excitation, dynamical state, and chemical makeup remain far less well understood. We know that these properties do vary dramatically and theoretical work posits a strong connection between the state of the gas, its ability to form stars, and the impact of stellar feedback. A next major step in the field will be to use sensitive cm-, mm-, and submm-wave spectroscopy and high resolution spectroscopic imaging to survey the state of cold gas across the whole local galaxy population. Such observations have pushed the capabilities of the current generation of telescopes. We highlight three critical elements for progress in the next decade: (1) robust support and aggressive development of ALMA, (2) the deployment of very large heterodyne receiver arrays on single dish telescopes, and (3) development of a new interferometric array that dramatically improves on current capabilities at cm- and mm-wavelengths (~ 1-115 GHz).
[27]  [pdf] - 1857290
Efficient methanol desorption in shear instability
Comments: 16 pages, 6 figures, submitted to the Astrophysical Journal. Comments welcome!
Submitted: 2019-03-27
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). While the dense core is outlined by NH2D emission, 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 one side of the core. 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.
[28]  [pdf] - 1851767
Comparing the Properties of GMCs in M33 from Simulations and Observations
Comments: 15 pages, 12 figures, accepted by MNRAS
Submitted: 2019-03-11
We compare the properties of clouds in simulated M33 galaxies to those observed in the real M33. We apply a friends of friends algorithm and CPROPS to identify clouds, as well as a pixel by pixel analysis. We obtain very good agreement between the number of clouds, and maximum mass of clouds. Both are lower than occurs for a Milky Way-type galaxy and thus are a function of the surface density, size and galactic potential of M33. We reproduce the observed dependence of molecular cloud properties on radius in the simulations, and find this is due to the variation in gas surface density with radius. The cloud spectra also show good agreement between the simulations and observations, but the exact slope and shape of the spectra depends on the algorithm used to find clouds, and the range of cloud masses included when fitting the slope. Properties such as cloud angular momentum, velocity dispersions and virial relation are also in good agreement between the simulations and observations, but do not necessarily distinguish between simulations of M33 and other galaxy simulations. Our results are not strongly dependent on the level of feedback used here (10 and 20%) although they suggest that 15% feedback efficiency may be optimal. Overall our results suggest that the molecular cloud properties are primarily dependent on the gas and mass surface density, and less dependent on the localised physics such as the details of stellar feedback, or the numerical code used.
[29]  [pdf] - 1864037
The Green Bank Ammonia Survey: A Virial Analysis of Gould Belt Clouds in Data Release 1
Comments: 35 pages, 8 tables, and 14 figures consisting of 16 .pdf files. Accepted for publication in the Astrophysical Journal
Submitted: 2019-03-08
We perform a virial analysis of starless dense cores in three nearby star-forming regions : L1688 in Ophiuchus, NGC 1333 in Perseus, and B18 in Taurus. Our analysis takes advantage of comprehensive kinematic information for the dense gas in all of these regions made publicly available through the Green Bank Ammonia Survey Data Release 1, which used to estimate internal support against collapse. We combine this information with ancillary data used to estimate other important properties of the cores, including continuum data from the James Clerk Maxwell Telescope Gould Belt Survey for core identification, core masses, and core sizes. Additionally, we used \textit{Planck} and \textit{Herschel}-based column density maps for external cloud weight pressure, and Five College Radio Astronomy Observatory $^{13}$CO observations for external turbulent pressure. Our self-consistent analysis suggests that many dense cores in all three star-forming regions are not bound by gravity alone, but rather require additional pressure confinement to remain bound. Unlike a recent, similar study in Orion~A, we find that turbulent pressure represents a significant portion of the external pressure budget. Our broad conclusion emphasizing the importance of pressure confinement in dense core evolution, however, agrees with earlier work.
[30]  [pdf] - 1838411
Relationship between the Line Width of the Atomic and Molecular ISM in M33
Comments: Accepted for publication in MNRAS; 20 pages, 15 figures
Submitted: 2019-02-13
We investigate how the spectral properties of atomic (HI) and molecular (H$_2$) gas, traced by CO(2-1), are related in M33 on $80$ pc scales. We find the HI and CO(2-1) velocity at peak intensity to be highly correlated, consistent with previous studies. By stacking spectra aligned to the velocity of HI peak intensity, we find that the CO line width ($\sigma_{\rm HWHM}=4.6\pm0.9$ km s$^{-1}$; $\sigma_{\rm HWHM}$ is the effective Gaussian width) is consistently smaller than the HI line width ($\sigma_{\rm HWHM}=6.6\pm0.1$ km s$^{-1}$), with a ratio of ${\sim}0.7$, in agreement with Druard et al. (2014). The ratio of the line widths remains less than unity when the data are smoothed to a coarser spatial resolution. In other nearby galaxies, this line width ratio is close to unity which has been used as evidence for a thick, diffuse molecular disk that is distinct from the thin molecular disk dominated by molecular clouds. The smaller line width ratio found here suggests that M33 has a marginal thick molecular disk. From modelling individual lines-of-sight, we recover a strong correlation between HI and CO line widths when only the HI located closest to the CO component is considered. The median line width ratio of the line-of-sight line widths is $0.56\pm0.01$. There is substantial scatter in the HI--CO(2-1) line width relation, larger than the uncertainties, that results from regional variations on $<500$ pc scales, and there is no significant trend in the line widths, or their ratios, with galactocentric radius. These regional line width variations may be a useful probe of changes in the local cloud environment or the evolutionary state of molecular clouds.
[31]  [pdf] - 1800052
The integrated properties of the molecular clouds from the JCMT CO(3-2) High Resolution Survey
Comments: 58 pages, 46 figure. Accepted for publication in MNRAS. A new, fully documented, version of SCIMES is public available on Github ( The molecular cloud integrated property catalog is public available on Zenodo (
Submitted: 2018-12-11
We define the molecular cloud properties of the Milky Way first quadrant using data from the JCMT CO(3-2) High Resolution Survey. We apply the Spectral Clustering for Interstellar Molecular Emission Segmentation (SCIMES) algorithm to extract objects from the full-resolution dataset, creating the first catalog of molecular clouds with a large dynamic range in spatial scale. We identify $>85\,000$ clouds with two clear sub-samples: $\sim35\,500$ well-resolved objects and $\sim540$ clouds with well-defined distance estimations. Only 35% of the cataloged clouds (as well as the total flux encompassed by them) appear enclosed within the Milky Way spiral arms. The scaling relationships between clouds with known distances are comparable to the characteristics of the clouds identified in previous surveys. However, these relations between integrated properties, especially from the full catalog, show a large intrinsic scatter ($\sim0.5$ dex), comparable to other cloud catalogs of the Milky Way and nearby galaxies. The mass distribution of molecular clouds follows a truncated-power law relationship over three orders of magnitude in mass with a form $dN/dM \propto M^{-1.7}$ with a clearly defined truncation at an upper mass of $M_0 \sim 3 \times 10^6~M_\odot$, consistent with theoretical models of cloud formation controlled by stellar feedback and shear. Similarly, the cloud population shows a power-law distribution of size with $dN/dR \propto R^{-2.8}$ with a truncation at $R_0 = 70$ pc.
[32]  [pdf] - 1811219
Do Spectroscopic Dense Gas Fractions Track Molecular Cloud Surface Densities?
Comments: 11 pages, 4 figures, accepted for publication in The Astrophysical Journal Letters
Submitted: 2018-11-19
We use ALMA and IRAM 30-m telescope data to investigate the relationship between the spectroscopically-traced dense gas fraction and the cloud-scale (120 pc) molecular gas surface density in five nearby, star-forming galaxies. We estimate the dense gas mass fraction at 650 pc and 2800 pc scales using the ratio of HCN (1-0) to CO (1-0) emission. We then use high resolution (120 pc) CO (2-1) maps to calculate the mass-weighted average molecular gas surface density within 650 pc or 2770 pc beam where the dense gas fraction is estimated. On average, the dense gas fraction correlates with the mass-weighted average molecular gas surface density. Thus, parts of a galaxy with higher mean cloud-scale gas surface density also appear to have a larger fraction of dense gas. The normalization and slope of the correlation do vary from galaxy to galaxy and with the size of the regions studied. This correlation is consistent with a scenario where the large-scale environment sets the gas volume density distribution, and this distribution manifests in both the cloud-scale surface density and the dense gas mass fraction.
[33]  [pdf] - 1805999
Forming Super Star Clusters in the Central Starburst of NGC 253
Comments: 20 pages, 5 figures, accepted for publication in the Astrophysical Journal
Submitted: 2018-04-05, last modified: 2018-10-24
NGC 253 hosts the nearest nuclear starburst. Previous observations show a region rich in molecular gas, with dense clouds associated with recent star formation. We used ALMA to image the 350 GHz dust continuum and molecular line emission from this region at 2 pc resolution. Our observations reveal ~14 bright, compact (~2-3 pc FWHM) knots of dust emission. Most of these sources are likely to be forming super star clusters (SSCs) based on their inferred dynamical and gas masses, association with 36 GHz radio continuum emission, and coincidence with line emission tracing dense, excited gas. One source coincides with a known SSC, but the rest remain invisible in Hubble near-infrared (IR) imaging. Our observations imply that gas still constitutes a large fraction of the overall mass in these sources. Their high brightness temperature at 350 GHz also implies a large optical depth near the peak of the IR spectral energy distribution. As a result, these sources may have large IR photospheres and the IR radiation force likely exceeds L/c. Still, their moderate observed velocity dispersions suggest that feedback from radiation, winds, and supernovae are not yet disrupting most sources. This mode of star formation appears to produce a large fraction of stars in the burst. We argue for a scenario in which this phase lasts ~1 Myr, after which the clusters shed their natal cocoons but continue to produce ionizing photons. The strong feedback that drives the observed cold gas and X-ray outflows likely occurs after the clusters emerge from this early phase.
[34]  [pdf] - 1767827
Science with an ngVLA: Molecular Clouds in Galaxies
Comments: To be published in the ASP Monograph Series, "Science With a Next-Generation VLA," ed. E. J. Murphy (ASP, San Francisco, CA)
Submitted: 2018-10-15
Stars form in cold clouds of predominantly molecular (H2) gas. We are just beginning to understand how the formation, properties, and destruction of these clouds varies across the universe. In this chapter, we describe how the thermal line imaging capabilities of the proposed next generation Very Large Array (ngVLA) could make major contributions to this field. Looking at CO emission, the proposed ngVLA would be able to quickly survey the bulk properties of molecular clouds across the whole nearby galaxy population. This includes many unique very nearby northern targets (e.g., Andromeda) inaccessible to ALMA. Such surveys offer a main observational constraint on the formation, destruction, lifetime, and star formation properties of clouds. Targeting specific regions, the ngVLA will also be able to heavily resolve clouds in the nearest galaxies. This will allow detailed studies of the substructure and kinematics --- and so the internal physics --- of clouds across different chemical and dynamical environments.
[35]  [pdf] - 1756142
Northern Galactic Molecular Cloud Clumps in Hi-GAL: Dense Gas Map and Environmental Trends
Comments: 13 pages, 12 figures, 2 tables, published in Monthly Notices of the Royal Astronomical Society
Submitted: 2018-09-24
In the quest to understand high-mass star formation, it is necessary to understand that from which the high-mass stars will form -- the dense molecular gas clouds and clumps. The Herschel infrared GALactic plane survey (Hi-GAL) is a comprehensive survey of thermal dust emission that can be used to characterize the properties and Galactic distribution of molecular gas. We have analyzed the survey maps within the Galactic longitude range $10^\circ < \ell < 56^\circ$ and have transformed these 2D maps into a 3D dense gas map of the Galactic plane using distance probability density functions. This range corresponds to the extent of the Galactic Ring Survey (GRS), which provided the majority of our line-of-sight velocities, and thus kinematic distances. In this section of Hi-GAL, we identify 19,886 clumps, out of which 10,124 have, at minimum, a kinematic distance, and 5,405 have their distances well constrained through Bayesian techniques. Of those with well-constrained distances, clump masses tend to decrease with Galactocentric radius, whereas clump radii are independent of Galactocentric radius.
[36]  [pdf] - 1748045
A Case Study of Triggered Star Formation in Cygnus X
Submitted: 2018-09-05
Radiative feedback from massive stars can potentially trigger star formation in the surrounding molecular gas. Inspired by the case of radiatively driven implosion in M16 or Eagle Nebula, we analyze a similar case of star formation observed in the Cygnus X region. We present new JCMT observations of $^{13}$CO(3-2) and C$^{18}$O(3-2) molecular lines of a cometary feature located at 50 pc north of the Cyg OB2 complex that was previously identified in $^{12}$CO(3-2) mapping. These data are combined with archival H$\alpha$, infrared, and radio continuum emission data, from which we measure the mass to be 110 M$_\odot$. We identify Cyg OB2 as the ionizing source. We measure the properties of two highly energetic molecular outflows and the photoionized rim. From this analysis, we argue the external gas pressure and gravitational energy dominate the internal pressure. The force balance along with previous simulation results and a close comparison with the case of Eagle Nebula favours a triggering scenario.
[37]  [pdf] - 1763078
The JCMT Gould Belt Survey: SCUBA-2 Data-Reduction Methods and Gaussian Source Recovery Analysis
Comments: Accepted for publication in ApJS
Submitted: 2018-08-23
The JCMT Gould Belt Survey was one of the first Legacy Surveys with the James Clerk Maxwell Telescope in Hawaii, mapping 47 square degrees of nearby (< 500 pc) molecular clouds in both dust continuum emission at 850 $\mu$m and 450 $\mu$m, as well as a more-limited area in lines of various CO isotopologues. While molecular clouds and the material that forms stars have structures on many size scales, their larger-scale structures are difficult to observe reliably in the submillimetre regime using ground-based facilities. In this paper, we quantify the extent to which three subsequent data-reduction methods employed by the JCMT GBS accurately recover emission structures of various size scales, in particular, dense cores which are the focus of many GBS science goals. With our current best data-reduction procedure, we expect to recover 100% of structures with Gaussian sigma sizes of $\le$30" and intensity peaks of at least five times the local noise for isolated peaks of emission. The measured sizes and peak fluxes of these compact structures are reliable (within 15% of the input values), but source recovery and reliability both decrease significantly for larger emission structures and for fainter peaks. Additional factors such as source crowding have not been tested in our analysis. The most recent JCMT GBS data release includes pointing corrections, and we demonstrate that these tend to decrease the sizes and increase the peak intensities of compact sources in our dataset, mostly at a low level (several percent), but occasionally with notable improvement.
[38]  [pdf] - 1730653
Kinematics of the Atomic ISM in M33 on 80 pc scales
Comments: MNRAS, in press. 30 pages, 24 figures
Submitted: 2018-06-28, last modified: 2018-08-08
We present new L-band ($1\mbox{--}2$ GHz) observations of the nearby spiral galaxy M33 with 80 pc resolution obtained with the Karl G. Jansky Very Large Array. The HI observations, combined with HI measurements from the Green Bank Telescope, improve the spectral resolution and sensitivity ($2.8$ K rms noise in a $0.2$ km s$^{-1}$ channel) compared to previous observations. We find individual profiles are usually non-Gaussian, harbouring line wings, multiple components, and asymmetries. Given this spectral complexity, we quantify the motions in the atomic ISM through moment analysis of the spectra and fits to aligned, stacked profiles. The measured value of the HI line width depends strongly on the method used, with the velocity stacked profiles aligned to the peak velocity giving the minimum value of $\sigma = 7$ km s$^{-1}$ and all other methods giving higher values ($\sigma\sim10$ km s$^{-1}$). All measurements of the line width show a shallow radial trend, with $\sigma$ decreasing by $\sim2$ km s$^{-1}$ from $R_{\rm gal}=0$ to $R_{\rm gal}=8$ kpc. We consider a number of energy sources that might maintain the line width against turbulent dissipation, but no single source is adequate. We find excess emission relative to a Gaussian in the stacked profile line wings, ranging from 9% to 26% depending on how the spectra are aligned. By splitting the line wings into symmetric and asymmetric components, we find that the lagging rotational disk accounts for one-third of the line wing flux. We also find emission far from the rotation axis of the galaxy in multiple discrete HI clouds, including a filament with a projected length of $\sim8$ kpc.
[39]  [pdf] - 1736466
A 50 pc scale view of star formation efficiency across NGC 628
Comments: 10 pages, 5 figures, accepted for publication in ApJ Letters
Submitted: 2018-07-30
Star formation is a multi-scale process that requires tracing cloud formation and stellar feedback within the local (<kpc) and global galaxy environment. We present first results from two large observing programs on ALMA and VLT/MUSE, mapping cloud scales (1arcsec = 47pc) in both molecular gas and star forming tracers across 90 kpc^2 of the central disk of NGC 628 to probe the physics of star formation. Systematic spatial offsets between molecular clouds and HII regions illustrate the time evolution of star-forming regions. Using uniform sampling of both maps on 50-500 pc scales, we infer molecular gas depletion times of 1-3 Gyr, but also find the increase of scatter in the star formation relation on small scales is consistent with gas and HII regions being only weakly correlated at the cloud (50 pc) scale. This implies a short overlap phase for molecular clouds and HII regions, which we test by directly matching our catalog of 1502 HII regions and 738 GMCs. We uncover only 74 objects in the overlap phase, and we find depletion times >1 Gyr, significantly longer than previously reported for individual star-forming clouds in the Milky Way. Finally, we find no clear trends that relate variations in the depletion time observed on 500 pc scales to physical drivers (metallicity, molecular and stellar mass surface density, molecular gas boundedness) on 50 pc scales.
[40]  [pdf] - 1714352
The Star Formation Efficiency per Free Fall Time in Nearby Galaxies
Comments: 10 pages, 3 figures, and 2 tables. Accepted for publication in ApJ Letters
Submitted: 2018-06-28
We estimate the star formation efficiency per gravitational free fall time, $\epsilon_{\rm ff}$, from observations of nearby galaxies with resolution matched to the typical size of a Giant Molecular Cloud. This quantity, $\epsilon_{\rm ff}$, is theoretically important but so far has only been measured for Milky Way clouds or inferred indirectly in a few other galaxies. Using new, high resolution CO imaging from the PHANGS-ALMA survey, we estimate the gravitational free-fall time at 60 to 120 pc resolution, and contrast this with the local molecular gas depletion time to estimate $\epsilon_{\rm ff}$. Assuming a constant thickness of the molecular gas layer ($H = 100$ pc) across the whole sample, the median value of $\epsilon_{\rm ff}$ in our sample is $0.7\%$. We find a mild scale-dependence, with higher $\epsilon_{\rm ff}$ measured at coarser resolution. Individual galaxies show different values of $\epsilon_{\rm ff}$, with the median $\epsilon_{\rm ff}$ ranging from $0.3\%$ to $2.6\%$. We find the highest $\epsilon_{\rm ff}$ in our lowest mass targets, reflecting both long free-fall times and short depletion times, though we caution that both measurements are subject to biases in low mass galaxies. We estimate the key systematic uncertainties, and show the dominant uncertainty to be the estimated line-of-sight depth through the molecular gas layer and the choice of star formation tracers.
[41]  [pdf] - 1712986
Assessing the Impact of Astrochemistry on Molecular Cloud Turbulence Statistics
Comments: 27 pages, 21 figures, accepted to ApJ
Submitted: 2018-05-24
We analyze hydrodynamic simulations of turbulent, star-forming molecular clouds that are post-processed with the photo-dissociation region astrochemistry code 3D-PDR. We investigate the sensitivity of 15 commonly applied turbulence statistics to post-processing assumptions, namely variations in gas temperature, abundance and external radiation field. We produce synthetic $^{12}$CO(1-0) and CI($^{3}$P$_{1}$-$^{3}$P$_{0}$) observations and examine how the variations influence the resulting emission distributions. To characterize differences between the datasets, we perform statistical measurements, identify diagnostics sensitive to our chemistry parameters, and quantify the statistic responses by using a variety of distance metrics. We find that multiple turbulent statistics are sensitive not only to the chemical complexity but also to the strength of the background radiation field. The statistics with meaningful responses include principal component analysis, spatial power spectrum and bicoherence. A few of the statistics, such as the velocity coordinate spectrum, are primarily sensitive to the type of tracer being utilized, while others, like the delta-variance, strongly respond to the background radiation field. Collectively, these findings indicate that more realistic chemistry impacts the responses of turbulent statistics and is necessary for accurate statistical comparisons between models and observed molecular clouds.
[42]  [pdf] - 1712960
The MALATANG Survey: the L_gas-L_IR correlation on sub-kiloparsec scale in six nearby star-forming galaxies as traced by HCN J=4-3 and HCO^+ J=4-3
Comments: 29 pages, 12 figures, accepted for publication in ApJ
Submitted: 2018-05-15
We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 \mu m/100 \mu m flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies.
[43]  [pdf] - 1712947
Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies
Comments: Accepted for publication in ApJ. 45 pages, 11 figures, 8 tables, 4 Appendices; key results summarized in Figure 10. Machine-readable table can be downloaded at prior to publication. For a brief video describing the main results of this paper, please see
Submitted: 2018-05-02
We measure the velocity dispersion, $\sigma$, and surface density, $\Sigma$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $\Sigma$ and $\sigma$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $\alpha_\mathrm{vir}\propto\sigma^2/\Sigma$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\propto\Sigma\,\sigma^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $\Sigma$, $\sigma$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $\sigma$-$\Sigma$ relation flattens at $\Sigma\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $\sigma$ for a given $\Sigma$ and high apparent $\alpha_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.
[44]  [pdf] - 1686808
Dense Gas, Dynamical Equilibrium Pressure, and Star Formation in Nearby Star-Forming Galaxies
Comments: 31 pages, 13 figures, accepted for publication in The Astrophysical Journal, email for access to data table before publication
Submitted: 2018-03-28
We use new ALMA observations to investigate the connection between dense gas fraction, star formation rate, and local environment across the inner region of four local galaxies showing a wide range of molecular gas depletion times. We map HCN (1-0), HCO$^+$ (1-0), CS (2-1), $^{13}$CO (1-0), and C$^{18}$O (1-0) across the inner few kpc of each target. We combine these data with short spacing information from the IRAM large program EMPIRE, archival CO maps, tracers of stellar structure and recent star formation, and recent HCN surveys by Bigiel et al. and Usero et al. We test the degree to which changes in the dense gas fraction drive changes in the SFR. $I_{HCN}/I_{CO}$ (tracing the dense gas fraction) correlates strongly with $I_{CO}$ (tracing molecular gas surface density), stellar surface density, and dynamical equilibrium pressure, $P_{DE}$. Therefore, $I_{HCN}/I_{CO}$ becomes very low and HCN becomes very faint at large galactocentric radii, where ratios as low as $I_{HCN}/I_{CO} \sim 0.01$ become common. The apparent ability of dense gas to form stars, $\Sigma_{SFR}/\Sigma_{dense}$ (where $\Sigma_{dense}$ is traced by the HCN intensity and the star formation rate is traced by a combination of H$\alpha$ and 24$\mu$m emission), also depends on environment. $\Sigma_{SFR}/\Sigma_{dense}$ decreases in regions of high gas surface density, high stellar surface density, and high $P_{DE}$. Statistically, these correlations between environment and both $\Sigma_{SFR}/\Sigma_{dense}$ and $I_{HCN}/I_{CO}$ are stronger than that between apparent dense gas fraction ($I_{HCN}/I_{CO}$) and the apparent molecular gas star formation efficiency $\Sigma_{SFR}/\Sigma_{mol}$. We show that these results are not specific to HCN.
[45]  [pdf] - 1627863
Detection of an OH 1665-MHz Maser in M33
Comments: 3 pages, 1 figure. Accepted in RNAAS
Submitted: 2018-01-31
We report one OH 1665-MHz detection in new observations with the NSF's Karl G. Jansky Very Large Array. To the authors' knowledge, this is the first OH maser detection in M33. The source occurs at a position of $\alpha=23.50069(9)^{\circ}$ and $\delta = 30.67984(9)^{\circ}$ (J2000) and is unresolved at the beam size ($9."7$). Two H$\alpha$ sources from the Moody et al. (2017) catalogue --- C1-1 and C1-2 --- are consistent within the beam size. We find a peak flux density of 24.1 mJy in a 1.5 km s$^{-1}$ channel. At the distance of M33, the source's specific luminosity is consistent with the brightest Galactic OH 1665-MHz masers.
[46]  [pdf] - 1626530
The JCMT Gould Belt Survey: A First Look at the Auriga-California Molecular Cloud with SCUBA-2
Comments: 39 pages (54 including the Appendix), 9 figures, 3 tables, in press in ApJ
Submitted: 2018-01-24
We present 850 and 450 micron observations of the dense regions within the Auriga-California molecular cloud using SCUBA-2 as part of the JCMT Gould Belt Legacy Survey to identify candidate protostellar objects, measure the masses of their circumstellar material (disk and envelope), and compare the star formation to that in the Orion A molecular cloud. We identify 59 candidate protostars based on the presence of compact submillimeter emission, complementing these observations with existing Herschel/SPIRE maps. Of our candidate protostars, 24 are associated with young stellar objects (YSOs) in the Spitzer and Herschel/PACS catalogs of 166 and 60 YSOs, respectively (177 unique), confirming their protostellar nature. The remaining 35 candidate protostars are in regions, particularly around LkHalpha 101, where the background cloud emission is too bright to verify or rule out the presence of the compact 70 micron emission that is expected for a protostellar source. We keep these candidate protostars in our sample but note that they may indeed be prestellar in nature. Our observations are sensitive to the high end of the mass distribution in Auriga-Cal. We find that the disparity between the richness of infrared star forming objects in Orion A and the sparsity in Auriga-Cal extends to the submillimeter, suggesting that the relative star formation rates have not varied over the Class II lifetime and that Auriga-Cal will maintain a lower star formation efficiency.
[47]  [pdf] - 1674832
Properties and rotation of molecular clouds in M 33
Comments: Accepted for publication in Astronomy and Astrophysics, 11 pages
Submitted: 2018-01-12
The sample of 566 molecular clouds identified in the CO(2--1) IRAM survey covering the disk of M~33 is explored in detail.The clouds were found using CPROPS and were subsequently catalogued in terms of their star-forming properties as non-star-forming (A), with embedded star formation (B), or with exposed star formation C.We find that the size-linewidth relation among the M~33 clouds is quite weak but, when comparing with clouds in other nearby galaxies, the linewidth scales with average metallicity.The linewidth and particularly the line brightness decrease with galactocentric distance.The large number of clouds makes it possible to calculate well-sampled cloud mass spectra and mass spectra of subsamples.As noted earlier, but considerably better defined here, the mass spectrum steepens (i.e. higher fraction of small clouds) with galactocentric distance.A new finding is that the mass spectrum of A clouds is much steeper than that of the star-forming clouds.Further dividing the sample, this difference is strong at both large and small galactocentric distances and the A vs C difference is a stronger effect than the inner/outer disk difference in mass spectra.Velocity gradients are identified in the clouds using standard techniques.The gradients are weak and are dominated by prograde rotation; the effect is stronger for the high signal-to-noise clouds.A discussion of the uncertainties is presented.The angular momenta are low but compatible with at least some simulations.The cloud and galactic gradients are similar; the cloud rotation periods are much longer than cloud lifetimes and comparable to the galactic rotation period.The rotational kinetic energy is 1-2\% of the gravitational potential energy and the cloud edge velocity is well below the escape velocity, such that cloud-scale rotation probably has little influence on the evolution of molecular clouds.
[48]  [pdf] - 1622438
Full-disc $^{13}$CO(1-0) mapping across nearby galaxies of the EMPIRE survey and the CO-to-H$_2$ conversion factor
Comments: accepted for publication in MNRAS
Submitted: 2018-01-09
Carbon monoxide (CO) provides crucial information about the molecular gas properties of galaxies. While $^{12}$CO has been targeted extensively, isotopologues such as $^{13}$CO have the advantage of being less optically thick and observations have recently become accessible across full galaxy discs. We present a comprehensive new dataset of $^{13}$CO(1-0) observations with the IRAM 30-m telescope of the full discs of 9 nearby spiral galaxies from the EMPIRE survey at a spatial resolution of $\sim$1.5kpc. $^{13}$CO(1-0) is mapped out to $0.7-1r_{25}$ and detected at high signal-to-noise throughout our maps. We analyse the $^{12}$CO(1-0)-to-$^{13}$CO(1-0) ratio ($\Re$) as a function of galactocentric radius and other parameters such as the $^{12}$CO(2-1)-to-$^{12}$CO(1-0) intensity ratio, the 70-to-160$\mu$m flux density ratio, the star-formation rate surface density, the star-formation efficiency, and the CO-to-H$_2$ conversion factor. We find that $\Re$ varies by a factor of 2 at most within and amongst galaxies, with a median value of 11 and larger variations in the galaxy centres than in the discs. We argue that optical depth effects, most likely due to changes in the mixture of diffuse/dense gas, are favored explanations for the observed $\Re$ variations, while abundance changes may also be at play. We calculate a spatially-resolved $^{13}$CO(1-0)-to-H$_2$ conversion factor and find an average value of $1.0\times10^{21}$ cm$^{-2}$ ($^{-1}$ over our sample with a standard deviation of a factor of 2. We find that $^{13}$CO(1-0) does not appear to be a good predictor of the bulk molecular gas mass in normal galaxy discs due to the presence of a large diffuse phase, but it may be a better tracer of the mass than $^{12}$CO(1-0) in the galaxy centres where the fraction of dense gas is larger.
[49]  [pdf] - 1641377
A Model for the Onset of Self-gravitation and Star Formation in Molecular Gas Governed by Galactic Forces: I. Cloud-scale Gas Motions
Comments: Accepted for publication in ApJ, 26 pages, 11 figures
Submitted: 2017-12-18
Modern extragalactic molecular gas surveys now reach the scales of star-forming giant molecular clouds (GMCs, 20-50 pc). Systematic variations in GMC properties with galaxy environment imply that clouds are not universally self-gravitating objects, decoupled from their surroundings. Here we reexamine the coupling of clouds to their environment and develop a model for 3D gas motions generated by forces arising with the galaxy gravitational potential defined by the background disk of stars and dark matter. We show that these motions can resemble or even exceed the motions needed to support gas against its own self-gravity throughout typical galaxy disks. The importance of the galactic potential in spiral arms and galaxy centers suggests that the response to self-gravity does not always dominate the motions of gas at GMC scales, with implications for observed gas kinematics, virial equilibrium and cloud morphology. We describe how a uniform treatment of gas motions in the plane and in the vertical direction synthesizes the two main mechanisms proposed to regulate star formation: vertical pressure equilibrium and shear/Coriolis forces as parameterized by Toomre Q~1. As the modeled motions are coherent and continually driven by the external potential, they represent support for the gas that is distinct from that conventionally attributed to turbulence, which decays rapidly and requires thus maintenance, e.g. via feedback from star formation. Thus our model suggests the galaxy itself can impose an important limit to star formation, as we explore in a second paper in this series.
[50]  [pdf] - 1630284
The EDGE-CALIFA survey: the influence of galactic rotation on the molecular depletion time across the Hubble sequence
Comments: 18 pages, 11 figures, accepted for publication in MNRAS
Submitted: 2017-12-10
We present a kpc-scale analysis of the relationship between the molecular depletion time ($\tau_\mathrm{dep}^\mathrm{mol}$) and the orbital time ($\tau_\mathrm{orb}$) across the field of 39 face-on local galaxies, selected from the EDGE-CALIFA sample. We find that, on average, 5% of the available molecular gas is converted into stars per orbital time, or $\tau_\mathrm{dep}^\mathrm{mol}\sim20\tau_\mathrm{orb}$. The resolved relation shows a scatter of $\sim0.5$ dex. The scatter is ascribable to galaxies of different morphologies that follow different $\tau_\mathrm{dep}^\mathrm{mol}-\tau_\mathrm{orb}$ relations which decrease in steepness from early- to late-types. The morphologies appear to be linked with the star formation rate surface density, the molecular depletion time, and the orbital time, but they do not correlate with the molecular gas content of the galaxies in our sample. We speculate that in our molecular gas rich, early-type galaxies, the morphological quenching (in particular the disc stabilization via shear), rather than the absence of molecular gas, is the main factor responsible for their current inefficient star formation.
[51]  [pdf] - 1626415
Mapping Jet-ISM Interactions in X-ray Binaries with ALMA: A GRS 1915$+$105 Case Study
Comments: 23 pages, 16 Figures, Accepted to MNRAS
Submitted: 2017-12-01
We present Atacama Large Millimetre/Sub-Millimetre Array (ALMA) observations of IRAS 19132+1035, a candidate jet-ISM interaction zone near the black hole X-ray binary (BHXB) GRS 1915+105. With these ALMA observations (combining data from the 12 m array and the Atacama Compact Array), we map the molecular line emission across the IRAS 19132+1035 region. We detect emission from the $^{12}$CO [$J=2-1$], $^{13}$CO [$\nu=0$, $J=2-1$], C$^{18}$O [$J=2-1$], ${\rm H}_{2}{\rm CO}$ [$J=3_{0,3}-2_{0,2}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,2}-2_{2,1}$], ${\rm H}_{2}{\rm CO}$ [$J=3_{2,1}-2_{2,0}$], SiO [$\nu=0$, $J=5-4$], CH$_3$OH [$J=4_{2,2}-3_{1,2}$], and CS [$\nu=0$, $J=5-4$] transitions. Given the morphological, spectral, and kinematic properties of this molecular emission, we present several lines of evidence that support the presence of a jet-ISM interaction at this site, including a jet-blown cavity in the molecular gas. This compelling new evidence identifies this site as a jet-ISM interaction zone, making GRS 1915$+$105 the third Galactic BHXB with at least one conclusive jet-ISM interaction zone. However, we find that this interaction occurs on much smaller scales than was postulated by previous work, where the BHXB jet does not appear to be dominantly powering the entire IRAS 19132+1035 region. Using estimates of the ISM conditions in the region, we utilize the detected cavity as a calorimeter to estimate the time-averaged power carried in the GRS 1915+105 jets of $(8.4^{+7.7}_{-8.1})\times10^{32}\,{\rm erg\,s}^{-1}$. Overall, our analysis demonstrates that molecular lines are excellent diagnostic tools to identify and probe jet-ISM interaction zones near Galactic BHXBs.
[52]  [pdf] - 1605074
The Green Bank Ammonia Survey: Unveiling the Dynamics of the Barnard 59 star-forming Clump
Comments: 18 pages, 18 figures
Submitted: 2017-11-06
Understanding the early stages of star formation is a research field of ongoing development, both theoretically and observationally. In this context, molecular data have been continuously providing observational constraints on the gas dynamics at different excitation conditions and depths in the sources. We have investigated the Barnard 59 core, the only active site of star formation in the Pipe Nebula, to achieve a comprehensive view of the kinematic properties of the source. These information were derived by simultaneously fitting ammonia inversion transition lines (1,1) and (2,2). Our analysis unveils the imprint of protostellar feedback, such as increasing line widths, temperature and turbulent motions in our molecular data. Combined with complementary observations of dust thermal emission, we estimate that the core is gravitationally bound following a virial analysis. If the core is not contracting, another source of internal pressure, most likely the magnetic field, is supporting it against gravitational collapse and limits its star formation efficiency.
[53]  [pdf] - 1598236
The Survey of Water and Ammonia in the Galactic Center (SWAG): Molecular Cloud Evolution in the Central Molecular Zone
Comments: 68 pages, 44 figures, accepted for publication in ApJS; slightly compressed images due to file size limitations, high resolution available at
Submitted: 2017-10-18
The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.
[54]  [pdf] - 1593709
The Green Bank Ammonia Survey: Observations of Hierarchical Dense Gas Structures in Cepheus-L1251
Comments: Accepted for publication in ApJ
Submitted: 2017-10-12
We use Green Bank Ammonia Survey observations of NH$_3$ (1,1) and (2,2) emission with 32'' FWHM resolution from a ~ 10 pc$^{2}$ portion of the Cepheus-L1251 molecular cloud to identify hierarchical dense gas structures. Our dendrogram analysis of the NH$_3$ data results in 22 top-level structures, which reside within 13 lower-level, parent structures. The structures are compact (0.01 pc $\lesssim R_{eff} \lesssim$ 0.1 pc) and are spatially correlated with the highest H$_2$ column density portions of the cloud. We also compare the ammonia data to a catalog of dense cores identified by higher-resolution (18.2'' FWHM) Herschel Space Observatory observations of dust continuum emission from Cepheus-L1251. Maps of kinetic gas temperature, velocity dispersion, and NH$_3$ column density, derived from detailed modeling of the NH$_3$ data, are used to investigate the stability and chemistry of the ammonia-identified and Herschel-identified structures. We show that the dust and dense gas in the structures have similar temperatures, with median $T_{dust}$ and $T_K$ measurements of 11.7 $\pm$ 1.1 K and 10.3 $\pm$ 2.0 K, respectively. Based on a virial analysis, we find that the ammonia-identified structures are gravitationally dominated, yet may be in or near a state of virial equilibrium. Meanwhile, the majority of the Herschel-identified dense cores appear to be not bound by their own gravity and instead confined by external pressure. CCS $(2_0-1_0)$ and HC$_5$N $(9-8)$ emission from the region reveal broader line widths and centroid velocity offsets when compared to the NH$_3$ (1,1) emission in some cases, likely due to these carbon-based molecules tracing the turbulent outer layers of the dense cores.
[55]  [pdf] - 1602607
ALMA Observations of a Quiescent Molecular Cloud in the Large Magellanic Cloud
Comments: Accepted by ApJ; 21 pages in AASTeX two-column style
Submitted: 2017-08-29, last modified: 2017-10-12
We present high-resolution (sub-parsec) observations of a giant molecular cloud in the nearest star-forming galaxy, the Large Magellanic Cloud. ALMA Band 6 observations trace the bulk of the molecular gas in $^{12}$CO(2-1) and high column density regions in $^{13}$CO(2-1). Our target is a quiescent cloud (PGCC G282.98-32.40, which we refer to as the "Planck cold cloud" or PCC) in the southern outskirts of the galaxy where star-formation activity is very low and largely confined to one location. We decompose the cloud into structures using a dendrogram and apply an identical analysis to matched-resolution cubes of the 30 Doradus molecular cloud (located near intense star formation) for comparison. Structures in the PCC exhibit roughly 10 times lower surface density and 5 times lower velocity dispersion than comparably sized structures in 30 Dor, underscoring the non-universality of molecular cloud properties. In both clouds, structures with relatively higher surface density lie closer to simple virial equilibrium, whereas lower surface density structures tend to exhibit super-virial line widths. In the PCC, relatively high line widths are found in the vicinity of an infrared source whose properties are consistent with a luminous young stellar object. More generally, we find that the smallest resolved structures ("leaves") of the dendrogram span close to the full range of line widths observed across all scales. As a result, while the bulk of the kinetic energy is found on the largest scales, the small-scale energetics tend to be dominated by only a few structures, leading to substantial scatter in observed size-linewidth relationships.
[56]  [pdf] - 1582499
The EDGE-CALIFA survey: Variations in the Molecular Gas Depletion Time in Local Galaxies
Comments: 17 pages, accepted in ApJ on Aug 24, 2017
Submitted: 2017-04-11, last modified: 2017-09-18
We present results from the EDGE survey, a spatially resolved CO(1-0) follow-up to CALIFA, an optical Integral Field Unit (IFU) survey of local galaxies. By combining the data products of EDGE and CALIFA, we study the variation in molecular gas depletion time ($\tau_{\rm dep}$) on kiloparsec scales in 52 galaxies. We divide each galaxy into two parts: the center, defined as the region within $0.1 \ R_{25}$, and the disk, defined as the region between $0.1$ and $0.7 \ R_{25}$. We find that 14 galaxies show a shorter $\tau_{\rm dep}$ ($\sim 1$ Gyr) in the center relative to that in the disk ($\tau_{\rm dep} \sim 2.4$ Gyrs), which means the central region in those galaxies is more efficient at forming stars per unit molecular gas mass. This finding implies that the centers with shorter $\tau_{\rm dep}$ resemble the intermediate regime between galactic disks and starburst galaxies. Furthermore, the central drop in $\tau_{\rm dep}$ is correlated with a central increase in the stellar surface density, suggesting that a shorter $\tau_{\rm dep}$ is associated with molecular gas compression by the stellar gravitational potential. We argue that varying the CO-to-H$_2$ conversion factor only exaggerates the central drop of $\tau_{\rm dep}$.
[57]  [pdf] - 1587772
CARMA Survey toward Infrared-bright Nearby Galaxies (STING). IV. Spatially Resolved 13CO in Spiral Galaxies
Comments: Accepted for publication in the Astrophysical Journal, 26 pages, 17 figures
Submitted: 2017-08-31
We present a $^{13}\mathrm{CO} (J = 1 \rightarrow 0)$ mapping survey of 12 nearby galaxies from the CARMA STING sample. The line intensity ratio $\mathcal{R} \equiv I[^{12}\mathrm{CO} (J = 1 \rightarrow 0)]/I[^{13}\mathrm{CO} (J = 1 \rightarrow 0)]$ is derived to study the variations in molecular gas properties. For 11 galaxies where it can be measured with high significance, the spatially resolved $\mathcal{R}$ on (sub-)kiloparsec scales varies by up to a factor of 3--5 within a galaxy. Lower $\mathcal{R}$ values are usually found in regions with weaker $^{12}\rm CO$. We attribute this apparent trend to a bias against measuring large $\mathcal{R}$ values when $^{12}\rm CO$ is weak. Limiting our analysis to the $^{12}\rm CO$ bright regions that are less biased, we do not find $\mathcal{R}$ on (sub)kpc scales correlate with galactocentric distance, velocity dispersion or the star formation rate. The lack of correlation between SFR and $\mathcal{R}$ indicates that the CO optical depth is not sensitive to stellar energy input, or that any such sensitivity is easily masked by other factors. Extending the analysis to all regions with $\rm ^{12}CO$ emission by spectral stacking, we find that 5 out of 11 galaxies show higher stacked $\mathcal{R}$ for galactocentric radii of $\gtrsim 1$ kpc and $\Sigma_{\mathrm{SFR}} \lesssim 0.1 \rm \ M_{sun} yr^{-1} kpc^{-2}$, which could result from a greater contribution from diffuse gas. Moreover, significant galaxy-to-galaxy variations are found in $\mathcal{R}$, but the global $\mathcal{R}$ does not strongly depend on dust temperature, inclination, or metallicity of the galaxy.
[58]  [pdf] - 1587231
The Green Bank Ammonia Survey: Dense Cores Under Pressure in Orion A
Comments: Accepted for publication in ApJ Aug 10, 2017
Submitted: 2017-08-17
We use gas temperature and velocity dispersion data from the Green Bank Ammonia Survey and core masses and sizes from the James Clerk Maxwell Telescope Gould Belt Survey to estimate the virial states of dense cores within the Orion A molecular cloud. Surprisingly, we find that almost none of the dense cores are sufficiently massive to be bound when considering only the balance between self-gravity and the thermal and non-thermal motions present in the dense gas. Including the additional pressure binding imposed by the weight of the ambient molecular cloud material and additional smaller pressure terms, however, suggests that most of the dense cores are pressure confined.
[59]  [pdf] - 1582436
The EDGE-CALIFA Survey: Interferometric Observations of 126 Galaxies with CARMA
Comments: Accepted by ApJ, 50 pp including figures
Submitted: 2017-04-08, last modified: 2017-08-16
We present interferometric CO observations made with the Combined Array for Millimeter-wave Astronomy (CARMA) of galaxies from the Extragalactic Database for Galaxy Evolution survey (EDGE). These galaxies are selected from the Calar Alto Legacy Integral Field Area (CALIFA) sample, mapped with optical integral field spectroscopy. EDGE provides good quality CO data (3$\sigma$ sensitivity $\Sigma_{\rm mol}$ $\sim$ 11 M$_\odot$ pc$^{-2}$ before inclination correction, resolution $\sim1.4$ kpc) for 126 galaxies, constituting the largest interferometric CO survey of galaxies in the nearby universe. We describe the survey, the data characteristics, the data products, and present initial science results. We find that the exponential scale-lengths of the molecular, stellar, and star-forming disks are approximately equal, and galaxies that are more compact in molecular gas than in stars tend to show signs of interaction. We characterize the molecular to stellar ratio as a function of Hubble type and stellar mass, present preliminary results on the resolved relations between the molecular gas, stars, and star formation rate, and discuss the dependence of the resolved molecular depletion time on stellar surface density, nebular extinction, and gas metallicity. EDGE provides a key dataset to address outstanding topics regarding gas and its role in star formation and galaxy evolution, which will be publicly available on completion of the quality assessment.
[60]  [pdf] - 1586047
Identifying Tools for Comparing Simulations and Observations of Spectral-line Data Cubes
Comments: 26 pages, 12 figures, accepted to MNRAS. If the data DOI link is not yet active, go here for the data:
Submitted: 2017-07-17
We present a statistical framework to compare spectral-line data cubes of molecular clouds and use the framework to perform an analysis of various statistical tools developed from methods proposed in the literature. We test whether our methods are sensitive to changes in the underlying physical properties of the clouds or whether their behaviour is governed by random fluctuations. We perform a set of 32 self-gravitating magnetohydrodynamic simulations that test all combinations of five physical parameters -- Mach number, plasma parameter, virial parameter, driving scales, and solenoidal driving fraction -- each of which can be set to a low or high value. We create mock observational data sets of ${\rm ^{13}CO}$(1-0) emission from each simulation. We compare these mock data to a those generated from a set of baseline simulations using pseudo-distance metrics based on 18 different statistical techniques that have previously been used to study molecular clouds. We analyze these results using methods from the statistical field of experimental design and find that several of the statistics can reliably track changes in the underlying physics. Our analysis shows that the interactions between parameters are often among the most significant effects. A small fraction of statistics are also sensitive to changes in magnetic field properties. We use this framework to compare the set of simulations to observations of three nearby star-forming regions: NGC 1333, Oph A, and IC 348. We find that no one simulation agrees significantly better with the observations, although there is evidence that the high Mach number simulations are more consistent with the observations.
[61]  [pdf] - 1579707
Towards a new classification of galaxies: principal component analysis of CALIFA circular velocity curves
Comments: Accepted for publication in MNRAS (Minor changes), 123 pages, 19 figures, 87 Tables (containing the basic properties of the 238 E1-Sdm galaxies; the five main Principal Component Eigenvectors; the five main Principal Components - PC_i; the Multi-Gaussian Expansion models - MGEs; the circular velocity curve models and their uncertainties)
Submitted: 2015-09-10, last modified: 2017-06-02
We present a galaxy classification system for 238 (E1-Sdm) CALIFA (Calar Alto Legacy Integral Field Area) galaxies based on the shapes and amplitudes of their circular velocity curves (CVCs). We infer the CVCs from the de-projected surface brightness of the galaxies, after scaling by a constant mass-to-light ratio based on stellar dynamics - solving axisymmetric Jeans equations via fitting the second velocity moment $V_{\mathrm{rms}}=\sqrt{V^2+\sigma^2}$ of the stellar kinematics. We use principal component analysis (PCA) applied to the CVC shapes to find characteristic features and use a $k$-means classifier to separate circular curves into classes. This objective classification method identifies four different classes, which we name slow-rising (SR), flat (FL), round-peaked (RP) and sharp-peaked (SP) circular curves. SR are typical for low-mass, late-type (Sb-Sdm), young, faint, metal-poor and disc-dominated galaxies. SP are typical for high-mass, early-type (E1-E7), old, bright, metal-rich and bulge-dominated galaxies. FL and RP appear presented by galaxies with intermediate mass, age, luminosity, metallicity, bulge-to-disk ratio and morphologies (E4-S0a, Sa-Sbc). The discrepancy mass factor, $f_d=1-M_{*}/M_{dyn}$, have the largest value for SR and SP classes ($\sim$ 74 per cent and $\sim$ 71 per cent, respectively) in contrast to the FL and RP classes (with $\sim$ 59 per cent and $\sim$ 61 per cent, respectively). Circular curve classification presents an alternative to typical morphological classification and appears more tightly linked to galaxy evolution.
[62]  [pdf] - 1582361
Thermal Feedback in the high-mass star and cluster forming region W51
Comments: Accepted in ApJ. Data are publicly available at . Repository of analysis scripts is at
Submitted: 2017-04-05, last modified: 2017-05-17
High-mass stars have generally been assumed to accrete most of their mass while already contracted onto the main sequence, but this hypothesis has not been observationally tested. We present ALMA observations of a 3 x 1.5 pc area in the W51 high-mass star-forming complex. We identify dust continuum sources and measure the gas and dust temperature through both rotational diagram modeling of CH3OH and brightness-temperature-based limits. The observed region contains three high-mass YSOs that appear to be at the earliest stages of their formation, with no signs of ionizing radiation from their central sources. The data reveal high gas and dust temperatures (T > 100 K) extending out to about 5000 AU from each of these sources. There are no clear signs of disks or rotating structures down to our 1000 AU resolution. The extended warm gas provides evidence that, during the process of forming, these high-mass stars heat a large volume and correspondingly large mass of gas in their surroundings, inhibiting fragmentation and therefore keeping a large reservoir available to feed from. By contrast, the more mature massive stars that illuminate compact H ii regions have little effect on their surrounding dense gas, suggesting that these main sequence stars have completed most or all of their accretion. The high luminosity of the massive protostars (L > 10^4 solar luminosities), combined with a lack of centimeter continuum emission from these sources, implies that they are not on the main sequence while they accrete the majority of their mass; instead, they may be bloated and cool.
[63]  [pdf] - 1573968
A Galactic Molecular Cloud Clump Catalog from Hi-GAL Data: Method and Initial Results Comparison to BGPS
Comments: 26 pages, 12 figures, 3 tables; published in the Astrophysical Journal
Submitted: 2017-05-05
As the precursors to stellar clusters, it is imperative that we understand the distribution and physical properties of dense molecular gas clouds and clumps. Such a study has been done with the ground-based Bolocam Galactic Plane Survey (BGPS). Now the Herschel infrared GALactic plane survey (Hi-GAL) allows us to do the same with higher quality data and complete coverage of the Galactic plane. We have made a pilot study comparing dense molecular gas clumps identified in the Hi-GAL and BGPS surveys, using six $2^\circ \times 2^\circ$ regions centered at Galactic longitudes of $\ell = 11^\circ$, 30$^\circ$, 41$^\circ$, 50$^\circ$, 202$^\circ$, and $217^\circ$. We adopted the BGPS methodology for identifying clumps and estimating distances, leading to 6198 clumps being identified in our substudy, with 995 of those having well-constrained distances. These objects were evenly distributed with Galactic longitude, a consequence of Hi-GAL being source confusion limited. These clumps range in mass from 10$^{-2}$ M$_\odot$ to 10$^{5}$ M$_\odot$, and have heliocentric distances of up to 16 kpc. When clumps found in both surveys are compared, we see that distances agree within 1 kpc and ratios of masses are of order unity. This serves as an external validation for BGPS and instills confidence as we move forward to cataloging the clumps from the entirety of Hi-GAL. In addition to the sources that were in common with BGPS, Hi-GAL found many additional sources, primarily due to the lack of atmospheric noise. We expect Hi-GAL to yield $2\times10^5$ clumps, with 20% having well-constrained distances, an order of magnitude above what was found in BGPS.
[64]  [pdf] - 1582936
Extreme Jet Ejections from the Black Hole X-ray Binary V404 Cygni
Comments: 25 pages, 11 figures, Accepted to MNRAS
Submitted: 2017-04-27
We present simultaneous radio through sub-mm observations of the black hole X-ray binary (BHXB) V404 Cygni during the most active phase of its June 2015 outburst. Our $4$ hour long set of overlapping observations with the Very Large Array, the Sub-millimeter Array, and the James Clerk Maxwell Telescope (SCUBA-2), covers 8 different frequency bands (including the first detection of a BHXB jet at $666 \,{\rm GHz}/450\mu m$), providing an unprecedented multi-frequency view of the extraordinary flaring activity seen during this period of the outburst. In particular, we detect multiple rapidly evolving flares, which reach Jy-level fluxes across all of our frequency bands. With this rich data set we performed detailed MCMC modeling of the repeated flaring events. Our custom model adapts the van der Laan synchrotron bubble model to include twin bi-polar ejections, propagating away from the black hole at bulk relativistic velocities, along a jet axis that is inclined to the line of sight. The emission predicted by our model accounts for projection effects, relativistic beaming, and the geometric time delay between the approaching and receding ejecta in each ejection event. We find that a total of 8 bi-polar, discrete jet ejection events can reproduce the emission that we observe in all of our frequency bands remarkably well. With our best fit model, we provide detailed probes of jet speed, structure, energetics, and geometry. Our analysis demonstrates the paramount importance of the mm/sub-mm bands, which offer a unique, more detailed view of the jet than can be provided by radio frequencies alone.
[65]  [pdf] - 1582744
The Green Bank Ammonia Survey (GAS): First Results of NH3 mapping the Gould Belt
Comments: 33 pages, 27 figures, accepted to ApJS. Datasets are publicly available:
Submitted: 2017-04-20
We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt star-forming regions with $A_V \gtrsim 7$ mag visible from the northern hemisphere in emission from NH$_3$ and other key molecular tracers. This first release includes the data for four regions in Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH$_3$ emission to dust continuum emission from Herschel, and find that the two tracers correspond closely. NH$_3$ is present in over 60\% of lines-of-sight with $A_V \gtrsim 7$ mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH$_3$ is detected toward ~ 40\% of lines-of-sight with $A_V \gtrsim 7$ mag. Moreover, we find that the NH$_3$ emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH$_3$ column densities through forward modeling of the hyperfine structure of the NH$_3$ (1,1) and (2,2) lines. We show that the NH$_3$ velocity dispersion, ${\sigma}_v$, and gas kinetic temperature, $T_K$, vary systematically between the regions included in this release, with an increase in both the mean value and spread of ${\sigma}_v$ and $T_K$ with increasing star formation activity. The data presented in this paper are publicly available.
[66]  [pdf] - 1564077
The Varying Mass Distribution of Molecular Clouds Across M83
Comments: 13 pages, accepted to MNRAS
Submitted: 2017-02-24
The work of Adamo et al. (2015) showed that the mass distributions of young massive stellar clusters were truncated above a maximum-mass scale in the nearby galaxy M83 and that this truncation mass varies with galactocentric radius. Here, we present a cloud-based analysis of ALMA CO($1\to 0$) observations of M83 to search for such a truncation mass in the molecular cloud population. We identify a population of 873 molecular clouds in M83 that is largely similar to those found in the Milky Way and Local Group galaxies, though clouds in the centre of the galaxy show high surface densities and enhanced turbulence, as is common for clouds in high-density nuclear environments. Like the young massive clusters, we find a maximum-mass scale for the molecular clouds that decreases radially in the galaxy. We find the most massive young massive cluster tracks the most massive molecular cloud with the cluster mass being $10^{-2}$ times that of the most massive molecular cloud. Outside the nuclear region of M83 ($R_{g}>0.5$ kpc), there is no evidence for changing internal conditions in the population of molecular clouds, with the average internal pressures, densities, and free-fall times remaining constant for the cloud population over the galaxy. This result is consistent with the bound cluster formation efficiency depending only on the large-scale properties of the ISM, rather than the internal conditions of individual clouds.
[67]  [pdf] - 1554012
The JCMT Nearby Galaxies Legacy Survey -- XI. -- Environmental Variations in the Atomic and Molecular Gas Radial Profiles of Nearby Spiral Galaxies
Comments: Accepted for publication in MNRAS, 12 pages, 6 figures
Submitted: 2017-02-07
We present an analysis of the radial profiles of a sample of 43 HI-flux selected spiral galaxies from the Nearby Galaxies Legacy Survey (NGLS) with resolved James Clerk Maxwell Telescope (JCMT) CO $J=3-2$ and/or Very Large Array (VLA) HI maps. Comparing the Virgo and non-Virgo populations, we confirm that the HI disks are truncated in the Virgo sample, even for these relatively HI-rich galaxies. On the other hand, the H$_{2}$ distribution is enhanced for Virgo galaxies near their centres, resulting in higher H$_{2}$ to HI ratios and steeper H$_{2}$ and total gas radial profiles. This is likely due to the effects of moderate ram pressure stripping in the cluster environment, which would preferentially remove low density gas in the outskirts while enhancing higher density gas near the centre. Combined with H$\alpha$ star formation rate data, we find that the star formation efficiency (SFR/H$_{2}$) is relatively constant with radius for both samples, but Virgo galaxies have a $\sim40\%$ lower star formation efficiency than non-Virgo galaxies.
[68]  [pdf] - 1534707
Dense Molecular Gas Tracers in the Outflow of the Starburst Galaxy NGC 253
Comments: Accepted for publication in the AAS Journals (ApJ)
Submitted: 2017-01-18
We present a detailed study of a molecular outflow feature in the nearby starburst galaxy NGC 253 using ALMA. We find that this feature is clearly associated with the edge of NGC 253's prominent ionized outflow, has a projected length of ~300 pc, with a width of ~50 pc and a velocity dispersion of ~40 km s^-1, consistent with an ejection from the disk about 1 Myr ago. The kinematics of the molecular gas in this feature can be interpreted (albeit not uniquely) as accelerating at a rate of 1 km s^-1 pc^-1. In this scenario, the gas is approaching escape velocity at the last measured point. Strikingly, bright tracers of dense molecular gas (HCN, CN, HCO+, CS) are also detected in the molecular outflow: We measure an HCN(1-0)/CO(1-0) line ratio of ~1/10 in the outflow, similar to that in the central starburst region of NGC 253 and other starburst galaxies. By contrast, the HCN/CO line ratio in the NGC 253 disk is significantly lower (~1/30), similar to other nearby galaxy disks. This strongly suggests that the streamer gas originates from the starburst, and that its physical state does not change significantly over timescales of ~1 Myr during its entrainment in the outflow. Simple calculations indicate that radiation pressure is not the main mechanism for driving the outflow. The presence of such dense material in molecular outflows needs to be accounted for in simulations of galactic outflows.
[69]  [pdf] - 1534694
The JCMT Gould Belt Survey: A First Look at IC 5146
Comments: 30 pages, 11 figures, 8 tables, accepted by ApJ
Submitted: 2017-01-17
We present 450 and 850 micron submillimetre continuum observations of the IC5146 star-forming region taken as part of the JCMT Gould Belt Survey. We investigate the location of bright submillimetre (clumped) emission with the larger-scale molecular cloud through comparison with extinction maps, and find that these denser structures correlate with higher cloud column density. Ninety-six individual submillimetre clumps are identified using FellWalker and their physical properties are examined. These clumps are found to be relatively massive, ranging from 0.5to 116 MSun with a mean mass of 8 MSun and a median mass of 3.7 MSun. A stability analysis for the clumps suggest that the majority are (thermally) Jeans stable, with M/M_J < 1. We further compare the locations of known protostars with the observed submillimetre emission, finding that younger protostars, i.e., Class 0 and I sources, are strongly correlated with submillimetre peaks and that the clumps with protostars are among the most Jeans unstable. Finally, we contrast the evolutionary conditions in the two major star-forming regions within IC5146: the young cluster associated with the Cocoon Nebula and the more distributed star formation associated with the Northern Streamer filaments. The Cocoon Nebula appears to have converted a higher fraction of its mass into dense clumps and protostars, the clumps are more likely to be Jeans unstable, and a larger fraction of these remaining clumps contain embedded protostars. The Northern Streamer, however, has a larger number of clumps in total and a larger fraction of the known protostars are still embedded within these clumps.
[70]  [pdf] - 1557673
Deuteration of ammonia in the starless core Ophiuchus/H-MM1
Comments: to appear in Astronomy & Astrophysics
Submitted: 2016-04-19, last modified: 2016-12-22
Ammonia and its deuterated isotopologues probe physical conditions in dense molecular cloud cores. With the aim of testing the current understanding of the spin-state chemistry of these molecules, we observed spectral lines of NH3, NH2D, NHD2, ND3, and N2D+ towards a dense, starless core in Ophiuchus with the APEX, GBT, and IRAM 30-m telescopes. The observations were interpreted using a gas-grain chemistry model combined with radiative transfer calculations. The chemistry model distinguishes between the different nuclear spin states of light hydrogen molecules, ammonia, and their deuterated forms. High deuterium fractionation ratios with NH2D/NH3=0.4, NHD2/NH2D=0.2, and ND3/NHD2=0.06 were found in the core. The observed ortho/para ratios of NH2D and NHD2 are close to the corresponding nuclear spin statistical weights. The chemistry model can approximately reproduce the observed abundances, but predicts uniformly too low ortho/para-NH2D, and too large ortho/para-NHD2 ratios. The longevity of N2H+ and NH3 in dense gas, which is prerequisite to their strong deuteration, can be attributed to the chemical inertia of N2 on grain surfaces. The discrepancies between the chemistry model and the observations are likely to be caused by the fact that the model assumes complete scrambling in principal gas-phase deuteration reactions of ammonia, which means that all the nuclei are mixed in reactive collisions. If, instead, these reactions occur through proton hop/hydrogen abstraction processes, statistical spin ratios are to be expected. The present results suggest that while the deuteration of ammonia changes with physical conditions and time, the nuclear spin ratios of ammonia isotopologues do not probe the evolutionary stage of a cloud.
[71]  [pdf] - 1533475
Optical Depth Estimates and Effective Critical Densities of Dense Gas Tracers in the Inner Parts of Nearby Galaxy Discs
Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2016-12-13
High critical density molecular lines like HCN(1-0) or HCO+(1-0) represent our best tool to study currently star-forming, dense molecular gas at extragalactic distances. The optical depth of these lines is a key ingredient to estimate the effective density required to excite emission. However, constraints on this quantity are even scarcer in the literature than measurements of the high density tracers themselves. Here, we combine new observations of HCN, HCO+ and HNC(1-0) and their optically thin isotopologues H13CN, H13CO+ and HN13C(1-0) to measure isotopologue line ratios. We use IRAM 30-m observations from the large program EMPIRE and new ALMA observations, which together target 6 nearby star-forming galaxies. Using spectral stacking techniques, we calculate or place strong upper limits on the HCN/H13CN, HCO+/H13CO+ and HNC/HN13C line ratios in the inner parts of these galaxies. Under simple assumptions, we use these to estimate the optical depths of HCN(1-0) and HCO+(1-0) to be \tau ~2-11 in the active, inner regions of our targets. The critical densities are consequently lowered to values between 5-20$\times 10^5$, 1-3$\times 10^5$ and 9$\times 10^4$ cm-3 for HCN, HCO+ and HNC, respectively. We study the impact of having different beam-filling factors, $\eta$, on these estimates and find that the effective critical densities decrease by a factor of $\frac{\eta_{12}}{\eta_{13}}\,\tau_{12}$. A comparison to existing work in NGC 5194 and NGC 253 shows HCN/H13CN and HCO+/H13CO+ ratios in agreement with our measurements within the uncertainties. The same is true for studies in other environments such as the Galactic Centre or nuclear regions of AGN-dominated nearby galaxies.
[72]  [pdf] - 1532894
Millimeter-Wave Line Ratios and Sub-beam Volume Density Distributions
Comments: 24 pages, 16 figure, Accepted for publication in the Astrophysical Journal, two online tables temporarily available at and
Submitted: 2016-11-29, last modified: 2016-12-01
We explore the use of mm-wave emission line ratios to trace molecular gas density when observations integrate over a wide range of volume densities within a single telescope beam. For observations targeting external galaxies, this case is unavoidable. Using a framework similar to that of Krumholz and Thompson (2007), we model emission for a set of common extragalactic lines from lognormal and power law density distributions. We consider the median density of gas producing emission and the ability to predict density variations from observed line ratios. We emphasize line ratio variations, because these do not require knowing the absolute abundance of our tracers. Patterns of line ratio variations have the prospect to illuminate the high-end shape of the density distribution, and to capture changes in the dense gas fraction and median volume density. Our results with and without a high density power law tail differ appreciably; we highlight better knowledge of the PDF shape as an important area. We also show the implications of sub-beam density distributions for isotopologue studies targeting dense gas tracers. Differential excitation often implies a significant correction to the naive case. We provide tabulated versions of many of our results, which can be used to interpret changes in mm-wave line ratios in terms of changes in the underlying density distributions.
[73]  [pdf] - 1532410
An Exploration of the Statistical Signatures of Stellar Feedback
Comments: 22 pages, 23 figures, accepted to ApJ
Submitted: 2016-11-02
All molecular clouds are observed to be turbulent, but the origin, means of sustenance, and evolution of the turbulence remain debated. One possibility is that stellar feedback injects enough energy into the cloud to drive observed motions on parsec scales. Recent numerical studies of molecular clouds have found that feedback from stars, such as protostellar outflows and winds, injects energy and impacts turbulence. We expand upon these studies by analyzing magnetohydrodynamic simulations of molecular clouds, including stellar winds, with a range of stellar mass-loss rates and magnetic field strengths. We generate synthetic $^{12}$CO(1-0) maps assuming that the simulations are at the distance of the nearby Perseus molecular cloud. By comparing the outputs from different initial conditions and evolutionary times, we identify differences in the synthetic observations and characterize these using common astrostatistics. We quantify the different statistical responses using a variety of metrics proposed in the literature. We find that multiple astrostatistics, including principal component analysis, the spectral correlation function and the velocity coordinate spectrum, are sensitive to changes in stellar mass-loss rates and/or time evolution. A few statistics, including the Cramer statistic and velocity coordinate spectrum, are sensitive to the magnetic field strength. These findings demonstrate that stellar feedback influences molecular cloud turbulence and can be identified and quantified observationally using such statistics.
[74]  [pdf] - 1528300
The JCMT Gould Belt Survey: First results from SCUBA-2 observations of the Cepheus Flare Region
Comments: Accepted for publication by MNRAS, 29 pages, 16 figures, 8 tables
Submitted: 2016-10-12
We present observations of the Cepheus Flare obtained as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Legacy Survey (GBLS) with the SCUBA-2 instrument. We produce a catalogue of sources found by SCUBA-2, and separate these into starless cores and protostars. We determine masses and densities for each of our sources, using source temperatures determined by the Herschel Gould Belt Survey. We compare the properties of starless cores in four different molecular clouds: L1147/58, L1172/74, L1251 and L1228. We find that the core mass functions for each region typically show shallower-than-Salpeter behaviour. We find that L1147/58 and L1228 have a high ratio of starless cores to Class II protostars, while L1251 and L1174 have a low ratio, consistent with the latter regions being more active sites of current star formation, while the former are forming stars less actively. We determine that, if modelled as thermally-supported Bonnor-Ebert spheres, most of our cores have stable configurations accessible to them. We estimate the external pressures on our cores using archival $^{13}$CO velocity dispersion measurements and find that our cores are typically pressure-confined, rather than gravitationally bound. We perform a virial analysis on our cores, and find that they typically cannot be supported against collapse by internal thermal energy alone, due primarily to the measured external pressures. This suggests that the dominant mode of internal support in starless cores in the Cepheus Flare is either non-thermal motions or internal magnetic fields.
[75]  [pdf] - 1513719
The inner mass distribution of late-type spiral galaxies from SAURON stellar kinematic maps
Comments: 44 pages, 60 figures, MNRAS accepted. The ADC-MCMC and JAM-MCMC python codes are available at: The Multi-Gaussian Expansion (MGE) results are also available in the Appendix
Submitted: 2016-09-27
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_{\phi}}/\sigma_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_{\phi}}/\sigma_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.
[76]  [pdf] - 1553879
The molecular gas mass of M33
Comments: Accepted for publication in Astronomy and Astrophysics
Submitted: 2016-09-13
[Abridged] Do some environments favor efficient conversion of molecular gas into stars? To answer this, we need to be able to estimate the H2 mass. Traditionally, this is done using CO and a few assumptions but the Herschel observations in the FIR make it possible to estimate the molecular gas mass independently of CO. Previous attempts to derive gas masses from dust emission suffered from biases. Generally, dust surface densities, HI column densities, and CO intensities are used to derive a gas-to-dust ratio (GDR) and the local CO intensity to H2 column density ratio (XCO), sometimes allowing for an additional CO-dark gas component (Kdark). We tested earlier methods, revealing degeneracies among the parameters, and then used a Bayesian formalism to derive the most likely values for each of the parameters mentioned above as a function of position in the nearby low metallicity spiral galaxy M33. The data are from the IRAM 30m CO(2-1) line, high-resolution HI and Herschel dust continuum observations. Solving for GDR, XCO, and Kdark in macro pixels 500 pc in size, we find that (i) allowing for CO-dark gas significantly improves fits; (ii) Kdark decreases with galactocentric distance; (iii) GDR is slightly higher than initially expected and increases with galactocentric distance; (iv) the total amount of dark gas closely follows the radially decreasing CO emission, as might be expected if the dark gas is H2 where CO is photodissociated. The total amount of H2, including dark gas, yields an average XCO of twice the galactic value of 2e20 cm^-2/(K km/s), 55% of this traced directly through CO. The rather constant fraction of dark gas suggests that there is no large population of diffuse H2 clouds (unrelated to GMCs) without CO emission. Unlike in large spirals, we detect no systematic radial trend in XCO, possibly linked to the absence of a radial decrease in CO line ratios.
[77]  [pdf] - 1501974
CARMA Large Area Star Formation Survey: Dense Gas in the Young L1451 Region of Perseus
Comments: Accepted to The Astrophysical Journal (ApJ), 45 pages, 24 figures (some with reduced resolution in this preprint); Project website is at
Submitted: 2016-06-28
We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.
[78]  [pdf] - 1475392
The JCMT Gould Belt Survey: A First Look at Southern Orion A with SCUBA-2
Comments: 31 Pages, 19 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society (MNRAS)
Submitted: 2016-06-28
We present the JCMT Gould Belt Survey's first look results of the southern extent of the Orion A Molecular Cloud ($\delta \leq -5\mathrm{:}31\mathrm{:}27.5$). Employing a two-step structure identification process, we construct individual catalogues for large-scale regions of significant emission labelled as islands and smaller-scale subregions called fragments using the 850 $\mu$m continuum maps obtained using SCUBA-2. We calculate object masses, sizes, column densities, and concentrations. We discuss fragmentation in terms of a Jeans instability analysis and highlight interesting structures as candidates for follow-up studies. Furthermore, we associate the detected emission with young stellar objects (YSOs) identified by Spitzer and Herschel. We find that although the population of active star-forming regions contains a wide variety of sizes and morphologies, there is a strong positive correlation between the concentration of an emission region and its calculated Jeans instability. There are, however, a number of highly unstable subregions in dense areas of the map that show no evidence of star formation. We find that only $\sim$72\% of the YSOs defined as Class 0+I and flat-spectrum protostars coincide with dense 850 $\mu$m emission structures (column densities $>3.7\times10^{21}\mathrm{\:cm}^{-2}$). The remaining 28\% of these objects, which are expected to be embedded in dust and gas, may be misclassified. Finally, we suggest that there is an evolution in the velocity dispersion of young stellar objects such that sources which are more evolved are associated with higher velocities.
[79]  [pdf] - 1470714
Cloud Structure of Galactic OB Cluster Forming Regions from Combining Ground and Space Based Bolometric Observations
Comments: ApJ accepted, Some figures in appendix are omitted
Submitted: 2016-06-24, last modified: 2016-06-26
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).
[80]  [pdf] - 1510269
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
Comments: Accepted for publication in the Astrophysical Journal, 36 pages (24+appendix), 21 figures (12+appendix), until publication high resolution version at
Submitted: 2016-06-22
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.
[81]  [pdf] - 1447767
The JCMT Gould Belt Survey: Evidence for Dust Grain Evolution in Perseus Star-forming Clumps
Comments: Accepted to ApJ, 15 figures, 3 tables
Submitted: 2016-05-19
The dust emissivity spectral index, $\beta$, is a critical parameter for deriving the mass and temperature of star-forming structures, and consequently their gravitational stability. The $\beta$ value is dependent on various dust grain properties, such as size, porosity, and surface composition, and is expected to vary as dust grains evolve. Here we present $\beta$, dust temperature, and optical depth maps of the star-forming clumps in the Perseus Molecular Cloud determined from fitting SEDs to combined Herschel and JCMT observations in the 160 $\mu$m, 250 $\mu$m, 350 $\mu$m, 500 $\mu$m, and 850 $\mu$m bands. Most of the derived $\beta$, and dust temperature values fall within the ranges of 1.0 - 2.7 and 8 - 20 K, respectively. In Perseus, we find the $\beta$ distribution differs significantly from clump to clump, indicative of grain growth. Furthermore, we also see significant, localized $\beta$ variations within individual clumps and find low $\beta$ regions correlate with local temperature peaks, hinting at the possible origins of low $\beta$ grains. Throughout Perseus, we also see indications of heating from B stars and embedded protostars, as well evidence of outflows shaping the local landscape.
[82]  [pdf] - 1443972
The JCMT Gould Belt Survey: Evidence for radiative heating and contamination in the W40 complex
Comments: 27 pages, 25 figures, 7 tables, 3 online catalogues
Submitted: 2016-05-16
We present SCUBA-2 450{\mu}m and 850{\mu}m observations of the W40 complex in the Serpens-Aquila region as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) of nearby star-forming regions. We investigate radiative heating by constructing temperature maps from the ratio of SCUBA-2 fluxes using a fixed dust opacity spectral index, {\beta} = 1.8, and a beam convolution kernel to achieve a common 14.8" resolution. We identify 82 clumps ranging between 10 and 36K with a mean temperature of 20{\pm}3K. Clump temperature is strongly correlated with proximity to the external OB association and there is no evidence that the embedded protostars significantly heat the dust. We identify 31 clumps that have cores with densities greater than 105cm{^{-3}}. Thirteen of these cores contain embedded Class 0/I protostars. Many cores are associated with bright-rimmed clouds seen in Herschel 70 {\mu}m images. From JCMT HARP observations of the 12CO 3-2 line, we find contamination of the 850{\mu}m band of up to 20 per cent. We investigate the free-free contribution to SCUBA-2 bands from large-scale and ultracompact H ii regions using archival VLA data and find the contribution is limited to individual stars, accounting for 9 per cent of flux per beam at 450 {\mu}m or 12 per cent at 850 {\mu}m in these cases. We conclude that radiative heating has potentially influenced the formation of stars in the Dust Arc sub-region, favouring Jeans stable clouds in the warm east and fragmentation in the cool west.
[83]  [pdf] - 1374382
Effects of Galactic Disc Inclination and Resolution on Observed GMC Properties and Larson's Scaling Relations
Comments: 12 pages, 5 figures, accepted for publication in MNRAS
Submitted: 2016-02-27
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.
[84]  [pdf] - 1354823
Dense Gas Fraction and Star Formation Efficiency Variations in the Antennae Galaxies
Comments: 9 pages, 4 figures, ApJ accepted
Submitted: 2015-10-19, last modified: 2016-02-09
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.
[85]  [pdf] - 1343007
The JCMT Nearby Galaxies Legacy Survey X. Environmental Effects on the Molecular Gas and Star Formation Properties of Spiral Galaxies
Comments: Accepted for publication in MNRAS, 25 pages, 14 figures
Submitted: 2015-12-17
We present a study of the molecular gas properties in a sample of 98 HI - flux selected spiral galaxies within $\sim25$ Mpc, using the CO $J=3-2$ line observed with the James Clerk Maxwell Telescope. We use the technique of survival analysis to incorporate galaxies with CO upper limits into our results. Comparing the group and Virgo samples, we find a larger mean H$_{2}$ mass in the Virgo galaxies, despite their lower mean HI mass. This leads to a significantly higher H$_{2}$ to HI ratio for Virgo galaxies. Combining our data with complementary H$\alpha$ star formation rate measurements, Virgo galaxies have longer molecular gas depletion times compared to group galaxies, due to their higher H$_{2}$ masses and lower star formation rates. We suggest that the longer depletion times may be a result of heating processes in the cluster environment or differences in the turbulent pressure. From the full sample, we find that the molecular gas depletion time has a positive correlation with the stellar mass, indicative of differences in the star formation process between low and high mass galaxies, and a negative correlation between the molecular gas depletion time and the specific star formation rate.
[86]  [pdf] - 1359019
The JCMT Gould Belt Survey: A First Look at Dense Cores in Orion B
Comments: Accepted for publication in ApJ; 28 pages, 19 figures. Data associated with the paper, including the full table 4, can be found at
Submitted: 2015-12-02
We present a first look at the SCUBA-2 observations of three sub-regions of the Orion B molecular cloud: LDN 1622, NGC 2023/2024, and NGC 2068/2071, from the JCMT Gould Belt Legacy Survey. We identify 29, 564, and 322 dense cores in L1622, NGC 2023/2024, and NGC 2068/2071 respectively, using the SCUBA-2 850 micron map, and present their basic properties, including their peak fluxes, total fluxes, and sizes, and an estimate of the corresponding 450 micron peak fluxes and total fluxes, using the FellWalker source extraction algorithm. Assuming a constant temperature of 20 K, the starless dense cores have a mass function similar to that found in previous dense core analyses, with a Salpeter-like slope at the high-mass end. The majority of cores appear stable to gravitational collapse when considering only thermal pressure; indeed, most of the cores which have masses above the thermal Jeans mass are already associated with at least one protostar. At higher cloud column densities, above 1-2 x 10^23 cm^-2, most of the mass is found within dense cores, while at lower cloud column densities, below 1 x 10^23 cm^-2, this fraction drops to 10% or lower. Overall, the fraction of dense cores associated with a protostar is quite small (<8%), but becomes larger for the densest and most centrally concentrated cores. NGC 2023 / 2024 and NGC 2068/2071 appear to be on the path to forming a significant number of stars in the future, while L1622 has little additional mass in dense cores to form many new stars.
[87]  [pdf] - 1411222
The Bolocam Galactic Plane Survey. XIV. Physical Properties of Massive Starless and Star Forming Clumps
Comments: Accepted for publication in ApJ; 33 pages; 22 figures; 7 tables
Submitted: 2015-11-27
We sort $4683$ molecular clouds between $10^\circ< \ell <65^\circ$ from the Bolocam Galactic Plane Survey based on observational diagnostics of star formation activity: compact $70$ $\mu{\rm m}$ sources, mid-IR color-selected YSOs, ${\rm H_2O}$ and ${\rm CH_3OH}$ masers, and UCHII regions. We also present a combined ${\rm NH_3}$-derived gas kinetic temperature and ${\rm H_2O}$ maser catalog for $1788$ clumps from our own GBT 100m observations and from the literature. We identify a subsample of $2223$ ($47.5\%$) starless clump candidates, the largest and most robust sample identified from a blind survey to date. Distributions of flux density, flux concentration, solid angle, kinetic temperature, column density, radius, and mass show strong ($>1$ dex) progressions when sorted by star formation indicator. The median starless clump candidate is marginally sub-virial ($\alpha \sim 0.7$) with $>75\%$ of clumps with known distance being gravitationally bound ($\alpha < 2$). These samples show a statistically significant increase in the median clump mass of $\Delta M \sim 170-370$ M$_\odot$ from the starless candidates to clumps associated with protostars. This trend could be due to (i) mass growth of the clumps at $\dot{M}\sim200-440$ Msun Myr$^{-1}$ for an average free-fall $0.8$ Myr time-scale, (ii) a systematic factor of two increase in dust opacity from starless to protostellar phases, (iii) and/or a variation in the ratio of starless to protostellar clump lifetime that scales as $\sim M^{-0.4}$. By comparing to the observed number of ${\rm CH_3OH}$ maser containing clumps we estimate the phase-lifetime of massive ($M>10^3$ M$_\odot$) starless clumps to be $0.37 \pm 0.08 \ {\rm Myr} \ (M/10^3 \ {\rm M}_\odot)^{-1}$; the majority ($M<450$ M$_\odot$) have phase-lifetimes longer than their average free-fall time.
[88]  [pdf] - 1530332
Bridging Galaxy Dynamics and Baryon Efficiency of 40 EDGE-CALIFA galaxies
Comments: 2 pages, 1 figure, Proceedings IAU Symposium No. 315, 2015
Submitted: 2015-10-15, last modified: 2015-10-25
We apply the Jeans Axisymmetric Multi-Gaussian Expansion method to the stellar kinematic maps of 40 Sa-Sd EDGE-CALIFA galaxies and derive their circular velocity curves (CVCs). The CVCs are classified using the Dynamical Classification method developed in Kalinova et al. (2015) . We also calculate the observational baryon efficiency, OBE, where $M_*/M_b=M_*/(M_*+M_{HI}+M_{H_2})$ of the galaxies using their stellar mass, total neutral hydrogen mass and total molecular gas from CO luminosities. Slow-rising, Flat and Round-peaked CVC types correspond to specific OBEs, stellar and dark matter (DM) halo mass values, while the Sharp-peaked CVCs span in the whole DM halo mass range of $10^{11}-10^{14} M_{\odot}$.
[89]  [pdf] - 1298215
Next Generation Very Large Array Memo No. 7 Science Working Group 2: "Galaxy Ecosystems": The Matter Cycle in and Around Galaxies
Comments: 40 pages. For more information or to get involved see
Submitted: 2015-10-21
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.
[90]  [pdf] - 1530338
The road toward a full, high resolution Molecular Cloud catalog of the Galaxy
Comments: 2 pages, 1 figure, Proceedings IAU Symposium No. 315, 2015
Submitted: 2015-10-19
The statistical description of Giant Molecular Cloud (GMC) properties relies heavily on the performance of automatic identification algorithms, which are often seriously affected by the survey design. The algorithm we designed, SCIMES (Spectral Clustering for Interstellar Molecular Emission Segmentation), is able to overcome some of these limitations by considering the cloud segmentation problem in the broad framework of the graph theory. The application of the code on the CO(3-2) High Resolution Survey (COHRS) data allowed for a robust decomposition of more than 12,000 objects in the Galactic Plane. Together with the wealth of Galactic Plane surveys of the recent years, this approach will help to open the door to a future, systematic cataloging of all discrete molecular features of our own Galaxy.
[91]  [pdf] - 1296363
Graph-based interpretation of the Molecular Interstellar Medium Segmentation
Comments: 25 pages, 18 figures, 3 tables, accepted for publication in MNRAS. The SCIMES package is available at and the documentation at
Submitted: 2015-10-14
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.
[92]  [pdf] - 1277169
What is a GMC? Are Observers and Simulators Discussing the Same Star-forming Clouds?
Comments: Accepted for publication in MNRAS
Submitted: 2015-08-07
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.
[93]  [pdf] - 1243880
Filament Identification through Mathematical Morphology
Comments: 18 pages, 14 figures, Accepted for publication in MNRAS. The FilFinder algorithm is available at and documentation at Maps of the identified filaments are available at and as online figures
Submitted: 2015-07-08
We present a new algorithm for detecting filamentary structure FilFinder. The algorithm uses the techniques of mathematical morphology for filament identification, presenting a complementary approach to current algorithms which use matched filtering or critical manifolds. Unlike other methods, FilFinder identifies filaments over a wide dynamic range in brightness. We apply the new algorithm to far infrared imaging data of dust emission released by the Herschel Gould Belt Survey team. Our preliminary analysis characterizes both filaments and fainter striations. We find a typical filament width of 0.09 pc across the sample, but the brightness varies from cloud to cloud. Several regions show a bimodal filament brightness distribution, with the bright mode (filaments) being an order of magnitude brighter than the faint mode (striations). Using the Rolling Hough Transform, we characterize the orientations of the striations in the data, finding preferred directions that agree with magnetic field direction where data are available. There is a suggestive but noisy correlation between typical filament brightness and literature values of the star formation rates for clouds in the Gould Belt.
[94]  [pdf] - 1286944
Variations in the Star Formation Efficiency of the Dense Molecular Gas across the Disks of Star-Forming Galaxies
Comments: Accepted for publication in the Astronomical Journal. 43 pages, 27 Figures
Submitted: 2015-06-01
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).
[95]  [pdf] - 1095663
The Bolocam Galactic Plane Survey. XIII. Physical Properties and Mass Functions of Dense Molecular Cloud Structures
Comments: Accepted for publication in ApJ; 32 pages, 21 figures
Submitted: 2015-04-06
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 $\alpha = 2.0\pm0.1$ for sources nearer than 6.5 kpc and $\alpha = 1.9\pm0.1$ for objects between 2 kpc and 10 kpc). The high-mass power-law indices are generally $1.85 \leq \alpha \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{\alpha} = 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%.
[96]  [pdf] - 974511
The Bolocam Galactic Plane Survey. XI. Temperatures and Substructure of Galactic Clumps Based on 350 micron Observations
Comments: 76 pages, 26 figures, 6 tables. Fig. 25 shows the complete sample of 107 maps. Accepted for publication in ApJS. Full tables 3 and 5 will be available online through ApJS. Typos corrected and other small changes on text after journal proof-review
Submitted: 2015-01-23, last modified: 2015-04-02
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.
[97]  [pdf] - 1231897
The JCMT Gould Belt Survey: SCUBA-2 observations of circumstellar disks in L 1495
Comments: 11 pages, 10 figures plus Appendix. MNRAS accepted
Submitted: 2015-02-27
We present 850$\mu$m and 450$\mu$m data from the JCMT Gould Belt Survey obtained with SCUBA-2 and characterise the dust attributes of Class I, Class II and Class III disk sources in L1495. We detect 23% of the sample at both wavelengths, with the detection rate decreasing through the Classes from I--III. The median disk mask is 1.6$\times 10^{-3}$M$_{\odot}$, and only 7% of Class II sources have disk masses larger than 20 Jupiter masses. We detect a higher proportion of disks towards sources with stellar hosts of spectral type K than spectral type M. Class II disks with single stellar hosts of spectral type K have higher masses than those of spectral type M, supporting the hypothesis that higher mass stars have more massive disks. Variations in disk masses calculated at the two wavelengths suggests there may be differences in dust opacity and/or dust temperature between disks with hosts of spectral types K to those with spectral type M.
[98]  [pdf] - 1449778
The JCMT Gould Belt Survey: First results from the SCUBA-2 observations of the Ophiuchus molecular cloud and a virial analysis of its prestellar core population
Comments: 34 pages, 20 figures, 5 tables. Accepted for publication in MNRAS
Submitted: 2015-02-20
In this paper we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N$_{2}$H$^{+}$ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialised. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N$_{2}$H$^{+}$ is a good tracer of the bound material of prestellar cores, although we find some evidence for N$_{2}$H$^{+}$ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C$^{18}$O and that traced by N$_{2}$H$^{+}$, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor-Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from southwest to northeast, indicating sequential star formation across the region.
[99]  [pdf] - 1095582
Measuring star formation with resolved observations: the test case of M33
Comments: 15 pages, 9 figures, 1 table. Accepted for publication in A&A
Submitted: 2015-02-04
Context. Measuring star formation at a local scale is important to constrain star formation laws. Yet, it is not clear whether and how the measure of star formation is affected by the spatial scale at which a galaxy is observed. Aims. We want to understand the impact of the resolution on the determination of the spatially resolved star formation rate (SFR) and other directly associated physical parameters such as the attenuation. Methods. We have carried out a multi-scale, pixel-by-pixel study of the nearby galaxy M33. Assembling FUV, Halpha, 8, 24, 70, and 100 micron maps, we have systematically compared the emission in individual bands with various SFR estimators from a resolution of 33 pc to 2084 pc. Results. We have found that there are strong, scale-dependent, discrepancies up to a factor 3 between monochromatic SFR estimators and Halpha+24 micron. The scaling factors between individual IR bands and the SFR show a strong dependence on the spatial scale and on the intensity of star formation. Finally, strong variations of the differential reddening between the nebular emission and the stellar continuum are seen, depending on the specific SFR (sSFR) and on the resolution. At the finest spatial scales, there is little differential reddening at high sSFR. The differential reddening increases with decreasing sSFR. At the coarsest spatial scales the differential reddening is compatible with the canonical value found for starburst galaxies. Conclusions. Our results confirm that monochromatic estimators of the SFR are unreliable at scales smaller than 1 kpc. Furthermore, the extension of local calibrations to high redshift galaxies presents non-trivial challenges as the properties of these systems may be poorly known.
[100]  [pdf] - 977887
Giant Molecular Clouds in the Early-Type Galaxy NGC4526
Comments: 21 pages, 16 figures. Accepted for publication in the Astrophysical Journal
Submitted: 2015-02-03
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 $\sigma_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 $\alpha \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 $\Sigma$ 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. $\sigma_v R^{-1/2} \propto \Sigma^{1/2}$.
[101]  [pdf] - 1449762
ALMA Multi-line Imaging of the Nearby Starburst Galaxy NGC 253
Comments: 20 pages, 10 figures, accepted to the Astrophysical Journal
Submitted: 2015-01-22
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 $\beta$ recombination lines that, like their $\alpha$ 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.
[102]  [pdf] - 918376
The dense gas mass fraction in the W51 cloud and its protoclusters
Comments: Accepted to A&A. 31 pages, 59 figures. Data available from and source code available from
Submitted: 2014-11-06, last modified: 2015-01-07
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 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.
[103]  [pdf] - 911134
The JCMT Gould Belt Survey: Evidence for radiative heating in Serpens MWC 297 and its influence on local star formation
Comments: 24 pages, 13 figures, 7 tables
Submitted: 2014-12-18, last modified: 2014-12-19
We present SCUBA-2 450micron and 850micron observations of the Serpens MWC 297 region, part of the JCMT Gould Belt Survey of nearby star-forming regions. Simulations suggest that radiative feedback influences the star-formation process and we investigate observational evidence for this by constructing temperature maps. Maps are derived from the ratio of SCUBA-2 fluxes and a two component model of the JCMT beam for a fixed dust opacity spectral index of beta = 1.8. Within 40 of the B1.5Ve Herbig star MWC 297, the submillimetre fluxes are contaminated by free-free emission with a spectral index of 1.03+-0.02, consistent with an ultra-compact HII region and polar winds/jets. Contamination accounts for 73+-5 per cent and 82+-4 per cent of peak flux at 450micron and 850micron respectively. The residual thermal disk of the star is almost undetectable at these wavelengths. Young Stellar Objects are confirmed where SCUBA-2 850micron clumps identified by the fellwalker algorithm coincide with Spitzer Gould Belt Survey detections. We identify 23 objects and use Tbol to classify nine YSOs with masses 0.09 to 5.1 Msun. We find two Class 0, one Class 0/I, three Class I and three Class II sources. The mean temperature is 15+-2K for the nine YSOs and 32+-4K for the 14 starless clumps. We observe a starless clump with an abnormally high mean temperature of 46+-2K and conclude that it is radiatively heated by the star MWC 297. Jeans stability provides evidence that radiative heating by the star MWC 297 may be suppressing clump collapse.
[104]  [pdf] - 1223541
Astrochemical Correlations in Molecular Clouds
Comments: Accepted by ApJ. 19 pages, 13 figures
Submitted: 2014-12-08
We investigate the spectral correlations between different species used to observe molecular clouds. We use hydrodynamic simulations and a full chemical network to study the abundances of over 150 species in typical Milky Way molecular clouds. We perform synthetic observations in order to produce emission maps of a subset of these tracers. We study the effects of different lines of sight and spatial resolution on the emission distribution and perform a robust quantitative comparison of the species to each other. We use the Spectral Correlation Function (SCF), which quantifies the root mean squared difference between spectra separated by some length scale, to characterize the structure of the simulated cloud in position-position-velocity (PPV) space. We predict the observed SCF for a broad range of observational tracers, and thus, identify homologous species. In particular, we show that the pairs C and CO, C$^{+}$ and CN, NH$_3$ and H$_2$CS have very similar SCFs. We measure the SCF slope variation as a function of beam size for all species and demonstrate that the beam size has a distinct effect on different species emission. However, for beams of up to 10", placing the cloud at 1 kpc, the change is not large enough to move the SCF slopes into different regions of parameter space. The results from this study provide observational guidance for choosing the best tracer to probe various cloud length scales.
[105]  [pdf] - 1223116
ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst
Comments: Accepted for publication in the Astrophysical Journal. 31 pages, 16 Figures
Submitted: 2014-11-11
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.
[106]  [pdf] - 921504
The Bolocam Galactic Plane Survey. XII. Distance Catalog Expansion Using Kinematic Isolation of Dense Molecular Cloud Structures With 13CO(1-0)
Comments: 28 pages, 19 figures. Accepted for publication in ApJ. Distance-Omnibus code available at
Submitted: 2014-11-10
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.
[107]  [pdf] - 1516791
The JCMT Legacy Survey of the Gould Belt: a molecular line study of the Ophiuchus molecular cloud
Comments: MNRAS in press 1 November 2014 - 29 pages
Submitted: 2014-11-05
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.
[108]  [pdf] - 1222555
CARMA Large Area Star Formation Survey: Structure and Kinematics of Dense Gas in Serpens Main
Comments: Accepted to ApJ. 38 pages, 16 figures, 5 tables
Submitted: 2014-10-13
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.
[109]  [pdf] - 1216683
CARMA Large Area Star Formation Survey: Project Overview with Analysis of Dense Gas Structure and Kinematics in Barnard 1
Comments: Accepted to The Astrophysical Journal (ApJ), 51 pages, 27 figures (some with reduced resolution in this preprint); Project website is at
Submitted: 2014-09-03
We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.
[110]  [pdf] - 1216503
The Interstellar Medium and Star Formation in Edge-On Galaxies. II. NGC 4157, 4565, and 5907
Comments: 26 pages, 26 figures, Accepted for publication in AJ
Submitted: 2014-08-25
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.
[111]  [pdf] - 853246
CARMA Large Area Star Formation Survey: Observational Analysis of Filaments in the Serpens South Molecular Cloud
Comments: 12 pages, 4 figures, published in ApJL (July 2014)
Submitted: 2014-07-02, last modified: 2014-07-17
We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.
[112]  [pdf] - 1208924
Power Spectrum Analysis of Polarized Emission from the Canadian Galactic Plane Survey
Comments: Accepted to ApJ; 17 pages
Submitted: 2014-04-10
Angular power spectra are calculated and presented for the entirety of the Canadian Galactic Plane Survey polarization dataset at 1.4 GHz covering an area of 1060 deg$^2$. The data analyzed are a combination of data from the 100-m Effelsberg Telescope, the 26-m Telescope at the Dominion Radio Astrophysical Observatory, and the Synthesis Telescope at the Dominion Radio Astrophysical Observatory, allowing all scales to be sampled down to arcminute resolution. The resulting power spectra cover multipoles from $\ell \approx 60$ to $\ell \approx 10^4$ and display both a power-law component at low multipoles and a flattening at high multipoles from point sources. We fit the power spectrum with a model that accounts for these components and instrumental effects. The resulting power-law indices are found to have a mode of 2.3, similar to previous results. However, there are significant regional variations in the index, defying attempts to characterize the emission with a single value. The power-law index is found to increase away from the Galactic plane. A transition from small-scale to large-scale structure is evident at $b= 9^{\circ}$, associated with the disk-halo transition in a 15$^{\circ}$ region around $l=108^{\circ}$. Localized variations in the index are found toward HII regions and supernova remnants, but the interpretation of these variations is inconclusive. The power in the polarized emission is anticorrelated with bright thermal emission (traced by H$\alpha$ emission) indicating that the thermal emission depolarizes background synchrotron emission.
[113]  [pdf] - 1202974
Comparing Simulated Emission from Molecular Clouds Using Experimental Design
Comments: Accepted by ApJ
Submitted: 2014-01-23
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.
[114]  [pdf] - 1180759
Shock Excited Molecules in NGC 1266: ULIRG conditions at the center of a Bulge Dominated Galaxy
Comments: Accepted for publication in ApJL 11/15/2013
Submitted: 2013-11-15
We investigate the far infrared spectrum of NGC 1266, a S0 galaxy that contains a massive reservoir of highly excited molecular gas. Using the SPIRE-FTS, we detect the $^{12}$CO ladder up to J=(13-12), [C I] and [N II] lines, and also strong water lines more characteristic of UltraLuminous IR Galaxies (ULIRGs). The 12CO line emission is modeled with a combination of a low-velocity C-shock and a PDR. Shocks are required to produce the H2O and most of the high-J 12CO emission. Despite having an infrared luminosity thirty times less than a typical ULIRG, the spectral characteristics and physical conditions of the ISM of NGC 1266 closely resemble those of ULIRGs, which often harbor strong shocks and large-scale outflows.
[115]  [pdf] - 1179482
CARMA Survey Toward Infrared-bright Nearby Galaxies (STING). III. The Dependence of Atomic and Molecular Gas Surface Densities on Galaxy Properties
Comments: 6 pages, to appear in ApJ Letters
Submitted: 2013-09-25
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.
[116]  [pdf] - 1173488
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
Comments: Accepted for publication in ApJS. The full velocity catalog and data tables (complete Tables 2,3, and 5) may be downloaded from Kinematic finder charts (0.5 deg x 0.5 deg) are also available at the same address
Submitted: 2013-08-19
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%).
[117]  [pdf] - 1251867
The Starburst-Driven Molecular Wind in NGC 253 and the Suppression of Star Formation
Comments: Appearing in the July 25 2013 issue of Nature
Submitted: 2013-07-23, last modified: 2013-07-29
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{\alpha} 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 {\eta}~1-3, establishing the importance of the starburst-driven wind in limiting the star formation activity and the final stellar content.
[118]  [pdf] - 1171639
The Bolocam Galactic Plane Survey IX: Data Release 2 and Outer Galaxy Extension
Comments: 36 pages, 16 figures, accepted to ApJS. Data available from
Submitted: 2013-05-24, last modified: 2013-07-26
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-\sigma 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.
[119]  [pdf] - 666734
The Bolocam Galactic Plane Survey. VIII. A Mid-Infrared Kinematic Distance Discrimination Method
Comments: Accepted for publication in ApJ; 29 pages, 16 figures. Data available prior to publication by contacting the authors
Submitted: 2013-04-16
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.
[120]  [pdf] - 1124270
Hierarchical Structure of Magnetohydrodynamic Turbulence In Position-Position-Velocity Space
Comments: submitted to ApJ
Submitted: 2012-06-20, last modified: 2013-04-10
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.
[121]  [pdf] - 1152051
Dense Gas in M33 (HerM33es)
Comments: 21 pages, 9 figures, published in A&A
Submitted: 2012-10-11, last modified: 2013-01-31
We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios.
[122]  [pdf] - 1158952
Molecular Gas and Star Formation in Nearby Disk Galaxies
Comments: 39 pages (incl. 10 page appendix), 19 figures, figures degraded for astro-ph. Until publication a higher resolution version available at - we suggest that one
Submitted: 2013-01-10
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.
[123]  [pdf] - 1152306
A Water Maser and Ammonia Survey of GLIMPSE Extended Green Objects (EGOs)
Comments: Astrophysical Journal, accepted. Emulateapj, 24 pages including 24 figures, plus 9 tables (including full content of online-only tables)
Submitted: 2012-10-19
We present the results of a Nobeyama 45-m water maser and ammonia survey of all 94 northern GLIMPSE Extended Green Objects (EGOs), a sample of massive young stellar objects (MYSOs) identified based on their extended 4.5 micron emission. We observed the ammonia (1,1), (2,2), and (3,3) inversion lines, and detect emission towards 97%, 63%, and 46% of our sample, respectively (median rms ~50 mK). The water maser detection rate is 68% (median rms ~0.11 Jy). The derived water maser and clump-scale gas properties are consistent with the identification of EGOs as young MYSOs. To explore the degree of variation among EGOs, we analyze subsamples defined based on MIR properties or maser associations. Water masers and warm dense gas, as indicated by emission in the higher-excitation ammonia transitions, are most frequently detected towards EGOs also associated with both Class I and II methanol masers. 95% (81%) of such EGOs are detected in water (ammonia(3,3)), compared to only 33% (7%) of EGOs without either methanol maser type. As populations, EGOs associated with Class I and/or II methanol masers have significantly higher ammonia linewidths, column densities, and kinetic temperatures than EGOs undetected in methanol maser surveys. However, we find no evidence for statistically significant differences in water maser properties (such as maser luminosity) among any EGO subsamples. Combining our data with the 1.1 mm continuum Bolocam Galactic Plane Survey, we find no correlation between isotropic water maser luminosity and clump number density. Water maser luminosity is weakly correlated with clump (gas) temperature and clump mass.
[124]  [pdf] - 1123969
The JCMT Nearby Galaxies Legacy Survey VIII. CO data and the L(CO3-2)-L(FIR) correlation in the SINGS sample
Comments: 37 pages, 107 figures; accepted for publication in MNRAS. Data and additional information are available at
Submitted: 2012-06-07
The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase.
[125]  [pdf] - 1123192
Cool and warm dust emission from M33 (HerM33es)
Comments: 12 pages, 14 figures, accepted for publication in A&A
Submitted: 2012-05-07
We study the far-infrared emission from the nearby spiral galaxy M33 in order to investigate the dust physical properties such as the temperature and the luminosity density across the galaxy. Taking advantage of the unique wavelength coverage (100, 160, 250, 350 and 500 micron) of the Herschel Space Observatory and complementing our dataset with Spitzer-IRAC 5.8 and 8 micron and Spitzer-MIPS 24 and 70 micron data, we construct temperature and luminosity density maps by fitting two modified blackbodies of a fixed emissivity index of 1.5. We find that the 'cool' dust grains are heated at temperatures between 11 and 28 K with the lowest temperatures found in the outskirts of the galaxy and the highest ones in the center and in the bright HII regions. The infrared/submillimeter total luminosity (5 - 1000 micron) is estimated to be 1.9x10^9 Lsun. 59% of the total luminosity of the galaxy is produced by the 'cool' dust grains (~15 K) while the rest 41% is produced by 'warm' dust grains (~55 K). The ratio of the cool-to-warm dust luminosity is close to unity (within the computed uncertainties), throughout the galaxy, with the luminosity of the cool dust being slightly enhanced in the center of the galaxy. Decomposing the emission of the dust into two components (one emitted by the diffuse disk of the galaxy and one emitted by the spiral arms) we find that the fraction of the emission in the disk in the mid-infrared (24 micron) is 21%, while it gradually rises up to 57% in the submillimeter (500 micron). We find that the bulk of the luminosity comes from the spiral arm network that produces 70% of the total luminosity of the galaxy with the rest coming from the diffuse dust disk. The 'cool' dust inside the disk is heated at a narrow range of temperatures between 18 and 15 K (going from the center to the outer parts of the galaxy).
[126]  [pdf] - 1117397
Low CO Luminosities in Dwarf Galaxies
Comments: Accepted for publication in the Astronomical Journal, 19 pages, 7 figures
Submitted: 2012-03-19
[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.
[127]  [pdf] - 1084714
CARMA Survey Toward Infrared-bright Nearby Galaxies (STING) II: Molecular Gas Star Formation Law and Depletion Time Across the Blue Sequence
Comments: To appear in ApJ, December 2011; 17 pages; 8 figures
Submitted: 2011-10-07, last modified: 2011-11-10
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.
[128]  [pdf] - 1083293
The Bolocam Galactic Plane Survey. VII. Characterizing the Properties of Massive Star-Forming Regions
Comments: 30 pages, 25 figures, 9 tables. Accepted for publication in ApJ. Full tables will be available via ApJ
Submitted: 2011-08-08
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.
[129]  [pdf] - 1076815
A Molecular Star Formation Law in the Atomic Gas Dominated Regime in Nearby Galaxies
Comments: Accepted for publication in the Astronomical Journal, 28 pages, 25 figures; footnote changed in version 2
Submitted: 2011-05-23, last modified: 2011-06-03
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.
[130]  [pdf] - 1076772
The Bolocam Galactic Plane Survey V: HCO+ and N2H+ Spectroscopy of 1.1 mm Dust Continuum Sources
Comments: 48 pages, 25 figures, Accepted for publication
Submitted: 2011-05-20
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$\sigma$. 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.
[131]  [pdf] - 1053489
Minimal HCN emission from Molecular Clouds in M33
Comments: Accepted for publication in MNRAS
Submitted: 2011-04-20
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.
[132]  [pdf] - 1052279
The CO-to-H2 Conversion Factor From Infrared Dust Emission Across the Local Group
Comments: Accepted for publication in the Astrophysical Journal, 15 pages, 7 figures
Submitted: 2011-02-22
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 \alpha_NGC6822 \approx 30 M_sun/pc^-2 (K km s^-1)^-1 and \alpha_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.
[133]  [pdf] - 1040636
CARMA Survey Toward Infrared-bright Nearby Galaxies (STING): Molecular Gas Star Formation Law in NGC4254
Comments: accepted for publication in ApJ, vol 729, March 10 2011 issue; 30 pages; 14 figures; revised version includes referee's comments; results unchanged
Submitted: 2010-09-16, last modified: 2011-01-25
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\alpha with azimuthally averaged extinction correction, mid-infrared 24 micron, combined H\alpha 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\alpha 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.
[134]  [pdf] - 270153
The Bolocam Galactic Plane Survey IV: 1.1 and 0.35 mm Dust Continuum Emission in the Galactic Center Region
Comments: 76 pages, 22 figures, published in ApJ:
Submitted: 2010-11-03
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 {\mu}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 {\mu}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 {\mu}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]
[135]  [pdf] - 1041699
The Bolocam Galactic Plane Survey: Survey Description and Data Reduction
Comments: Accepted for publication in Astrophysical Journal Supplement
Submitted: 2010-11-02
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.
[136]  [pdf] - 1041278
Unusually Luminous Giant Molecular Clouds in the Outer Disk of M33
Comments: Accepted for Publication in ApJ; 7 pages, 4 figures
Submitted: 2010-10-13
We use high spatial resolution (~7pc) CARMA observations to derive detailed properties for 8 giant molecular clouds (GMCs) at a galactocentric radius corresponding to approximately two CO scale lengths, or ~0.5 optical radii (r25), in the Local Group spiral galaxy M33. At this radius, molecular gas fraction, dust-to-gas ratio and metallicity are much lower than in the inner part of M33 or in a typical spiral galaxy. This allows us to probe the impact of environment on GMC properties by comparing our measurements to previous data from the inner disk of M33, the Milky Way and other nearby galaxies. The outer disk clouds roughly fall on the size-linewidth relation defined by extragalactic GMCs, but are slightly displaced from the luminosity-virial mass relation in the sense of having high CO luminosity compared to the inferred virial mass. This implies a different CO-to-H2 conversion factor, which is on average a factor of two lower than the inner disk and the extragalactic average. We attribute this to significantly higher measured brightness temperatures of the outer disk clouds compared to the ancillary sample of GMCs, which is likely an effect of enhanced radiation levels due to massive star formation in the vicinity of our target field. Apart from brightness temperature, the properties we determine for the outer disk GMCs in M33 do not differ significantly from those of our comparison sample. In particular, the combined sample of inner and outer disk M33 clouds covers roughly the same range in size, linewidth, virial mass and CO luminosity than the sample of Milky Way GMCs. When compared to the inner disk clouds in M33, however, we find even the brightest outer disk clouds to be smaller than most of their inner disk counterparts. This may be due to incomplete sampling or a potentially steeper cloud mass function at larger radii.
[137]  [pdf] - 1034821
The Scale Dependence of the Molecular Gas Depletion Time in M33
Comments: Accepted for publication in the Astrophysical Journal, 9 pages, 5 figures
Submitted: 2010-09-08
We study the Local Group spiral galaxy M33 to investigate how the observed scaling between the (kpc-averaged) surface density of molecular gas (\Sigma_H2) and recent star formation rate (\Sigma_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 (\tau_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 \tau_dep \sim 1 Gyr, with no dependence on type of region targeted. Below these scales, \tau_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 \tau_dep (i.e., high CO-to-Halpha flux ratio) and small apertures targeted toward Halpha peaks have very short \tau_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 \tau_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 \tau_dep (and thus the CO-to-Halpha ratio).
[138]  [pdf] - 1032775
The Bolocam Galactic Plane Survey -- III. Characterizing Physical Properties of Massive Star-Forming Regions in the Gemini OB1 Molecular Cloud
Comments: 26 pages, 19 figures, 8 tables. Accepted for publication by ApJ
Submitted: 2010-05-26
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.
[139]  [pdf] - 1032597
Cool gas and dust in M33: Results from the Herschel M33 extended survey (HERM33ES)
Comments: 5 pages, 5 figures Accepted for publication in Astronomy and Astrophysics
Submitted: 2010-05-19
We present an analysis of the first space-based far-IR-submm observations of M 33, which measure the emission from the cool dust and resolve the giant molecular cloud complexes. With roughly half-solar abundances, M33 is a first step towards young low-metallicity galaxies where the submm may be able to provide an alternative to CO mapping to measure their H$_2$ content. In this Letter, we measure the dust emission cross-section $\sigma$ using SPIRE and recent CO and \HI\ observations; a variation in $\sigma$ is present from a near-solar neighborhood cross-section to about half-solar with the maximum being south of the nucleus. Calculating the total H column density from the measured dust temperature and cross-section, and then subtracting the \HI\ column, yields a morphology similar to that observed in CO. The H$_2$/\HI\ mass ratio decreases from about unity to well below 10% and is about 15% averaged over the optical disk. The single most important observation to reduce the potentially large systematic errors is to complete the CO mapping of M 33.
[140]  [pdf] - 1032481
PACS and SPIRE photometer maps of M33: First results of the Herschel M33 extended survey (HERM33ES)
Comments: 5 pages, 3 figures, accepted for publication in the A&A Herschel Special Issue
Submitted: 2010-05-14
Within the framework of the HERM33ES key project, we are studying the star forming interstellar medium in the nearby, metal-poor spiral galaxy M33, exploiting the high resolution and sensitivity of Herschel. We use PACS and SPIRE maps at 100, 160, 250, 350, and 500 micron wavelength, to study the variation of the spectral energy distributions (SEDs) with galacto-centric distance. Detailed SED modeling is performed using azimuthally averaged fluxes in elliptical rings of 2 kpc width, out to 8 kpc galacto-centric distance. Simple isothermal and two-component grey body models, with fixed dust emissivity index, are fitted to the SEDs between 24 and 500 micron using also MIPS/Spitzer data, to derive first estimates of the dust physical conditions. The far-infrared and submillimeter maps reveal the branched, knotted spiral structure of M33. An underlying diffuse disk is seen in all SPIRE maps (250-500 micron). Two component fits to the SEDs agree better than isothermal models with the observed, total and radially averaged flux densities. The two component model, with beta fixed at 1.5, best fits the global and the radial SEDs. The cold dust component clearly dominates; the relative mass of the warm component is less than 0.3% for all the fits. The temperature of the warm component is not well constrained and is found to be about 60K plus/minus 10K. The temperature of the cold component drops significantly from about 24K in the inner 2 kpc radius to 13K beyond 6 kpc radial distance, for the best fitting model. The gas-to-dust ratio for beta=1.5, averaged over the galaxy, is higher than the solar value by a factor of 1.5 and is roughly in agreement with the subsolar metallicity of M33.
[141]  [pdf] - 1025194
Direct observation of a sharp transition to coherence in Dense Cores
Comments: Accepted for publication in ApJ Letters. 6 pages, 5 figures
Submitted: 2010-02-15
We present NH3 observations of the B5 region in Perseus obtained with the Green Bank Telescope (GBT). The map covers a region large enough (~11'x14') that it contains the entire dense core observed in previous dust continuum surveys. The dense gas traced by NH3(1,1) covers a much larger area than the dust continuum features found in bolometer observations. The velocity dispersion in the central region of the core is small, presenting subsonic non-thermal motions which are independent of scale. However, it is thanks to the coverage and high sensitivity of the observations that we present the detection, for the first time, of the transition between the coherent core and the dense but more turbulent gas surrounding it. This transition is sharp, increasing the velocity dispersion by a factor of 2 in less than 0.04 pc (the 31" beam size at the distance of Perseus, ~250 pc). The change in velocity dispersion at the transition is ~3 km/s/pc. The existence of the transition provides a natural definition of dense core: the region with nearly-constant subsonic non-thermal velocity dispersion. From the analysis presented here we can not confirm nor rule out a corresponding sharp density transition.
[142]  [pdf] - 1024808
The Effect of Projection on Derived Mass-Size and Linewidth-Size Relationships
Comments: 25 pages, including 7 Figures; Accepted for publication in ApJ
Submitted: 2010-01-25, last modified: 2010-02-10
Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound.
[143]  [pdf] - 1024774
Planetary nebulae in M33: probes of AGB nucleosynthesis and ISM abundances
Comments: 21 pages, 20 figures. Accepted for publication in MNRAS
Submitted: 2010-01-22
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.
[144]  [pdf] - 1024739
Dense Gas Tracers in Perseus: Relating the N2H+, NH3, and Dust Continuum Properties of Pre- and Proto-Stellar Cores
Comments: ApJ accepted. 45 pages, 8 figures
Submitted: 2010-01-20
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.
[145]  [pdf] - 901918
Ionized Gas Towards Molecular Clumps: Physical Properties of Massive Star Forming Regions
Comments: 67 pages, 16 figures, 8 tables, accepted for publication in ApJ
Submitted: 2009-10-01
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.
[146]  [pdf] - 901738
The Bolocam Galactic Plane Survey -- II. Catalog of the Image Data
Comments: Accepted by ApJS. Data release and full resolution manuscript at
Submitted: 2009-09-15, last modified: 2009-09-17
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.
[147]  [pdf] - 26070
Variable Sodium Absorption in a Low-Extinction Type Ia Supernova
Comments: 16 pages, 10 figures (8 in color), 5 tables. Accepted for publication in ApJ
Submitted: 2009-07-07
Recent observations have revealed that some Type Ia supernovae exhibit narrow, time-variable Na I D absorption features. The origin of the absorbing material is controversial, but it may suggest the presence of circumstellar gas in the progenitor system prior to the explosion, with significant implications for the nature of the supernova progenitors. We present the third detection of such variable absorption, based on six epochs of high-resolution spectroscopy of the Type Ia supernova SN 2007le from Keck and the HET. The data span ~3 months, from 5 days before maximum light to 90 days after maximum. We find that one component of the Na D absorption lines strengthened significantly with time, indicating a total column density increase of ~2.5 x 10^12 cm^-2. The changes are most prominent after maximum light rather than at earlier times when the UV flux from the SN peaks. As with SN 2006X, we detect no change in the Ca II H&K lines over the same time period, rendering line-of-sight effects improbable and suggesting a circumstellar origin for the absorbing material. Unlike the previous two SNe exhibiting variable absorption, SN 2007le is not highly reddened (E_B-V = 0.27 mag), also pointing toward circumstellar rather than interstellar absorption. Photoionization models show that the data are consistent with a dense (10^7 cm^-3) cloud or clouds of gas located ~0.1 pc from the explosion. These results broadly support the single-degenerate scenario previously proposed to explain the variable absorption, with mass loss from a nondegenerate companion star responsible for providing the circumstellar gas. We also present tentative evidence for narrow Halpha emission associated with the SN, which will require followup observations at late times to confirm. [abridged]
[148]  [pdf] - 24820
The Perils of Clumpfind: The Mass Spectrum of Sub-structures in Molecular Clouds
Comments: 5 pages, 3 figures. Accepted for publication in ApJ Letters
Submitted: 2009-06-01
We study the mass spectrum of sub-structures in the Perseus Molecular Cloud Complex traced by 13CO (1-0), finding that $dN/dM\propto M^{-2.4}$ for the standard Clumpfind parameters. This result does not agree with the classical $dN/dM\propto M^{-1.6}$. To understand this discrepancy we study the robustness of the mass spectrum derived using the Clumpfind algorithm. Both 2D and 3D Clumpfind versions are tested, using 850 $\mu$m dust emission and 13CO spectral-line observations of Perseus, respectively. The effect of varying threshold is not important, but varying stepsize produces a different effect for 2D and 3D cases. In the 2D case, where emission is relatively isolated (associated with only the densest peaks in the cloud), the mass spectrum variability is negligible compared to the mass function fit uncertainties. In the 3D case, however, where the 13CO emission traces the bulk of the molecular cloud, the number of clumps and the derived mass spectrum are highly correlated with the stepsize used. The distinction between "2D" and "3D" here is more importantly also a distinction between "sparse" and "crowded" emission. In any "crowded" case, Clumpfind should not be used blindly to derive mass functions. Clumpfind's output in the "crowded" case can still offer a statistical description of emission useful in inter-comparisons, but the clump-list should not be treated as a robust region decomposition suitable to generate a physically-meaningful mass function. We conclude that the 13CO mass spectrum depends on the observations resolution, due to the hierarchical structure of MC.
[149]  [pdf] - 315575
Dense Cores in Perseus: The Influence of Stellar Content and Cluster Environment
Comments: 28 pages, 17 figures. Accepted to ApJ
Submitted: 2009-02-03
We present the chemistry, temperature, and dynamical state of a sample of 193 dense cores or core candidates in the Perseus Molecular cloud and compare the properties of cores associated with young stars and clusters with those which are not. The combination of our NH3 and CCS observations with previous millimeter, sub-millimeter, and Spitzer data available for this cloud enable us both to determine core properties precisely and to accurately classify cores as starless or protostellar. The properties of cores in different cluster environments and before-and-after star formation provide important constraints on simulations of star-formation, particularly under the paradigm that the essence of star formation is set by the turbulent formation of prestellar cores. We separate the influence of stellar content from that of cluster environment and find that cores within clusters have (1) higher kinetic temperatures and (2) lower fractional abundances of CCS and NH3. Cores associated with protostars have (1) slightly higher kinetic temperatures (2) higher NH3 excitation temperatures), (3) are at higher column density, have (4) slightly more non-thermal/turbulent NH3 linewidths, have (5) higher masses and have (6) lower fractional abundance of CCS. We find that neither cluster environment nor protostellar content makes a significant difference to the dynamical state of cores as estimated by the virial parameter -- most cores in each category are gravitationally bound. Overall, cluster environment and protostellar content have a smaller influence on the properties of the cores than is typically assumed, and the variation within categories is larger than the differences between categories.
[150]  [pdf] - 17478
Resolved Properties of Extragalactic Giant Molecular Clouds
Comments: 4 pages, 3 figures, to appear in proceedings for IAU 255, "Low-Metallicity Star Formation: From the First Stars to Dwarf Galaxies"
Submitted: 2008-10-15
Giant molecular clouds (GMCs) are the major reservoirs of molecular gas in galaxies, and the starting point for star formation. As such, their properties play a key role in setting the initial conditions for the formation of stars. We present a comprehensive combined inteferometric/single-dish study of the resolved GMC properties in a number of extragalactic systems, including both normal and dwarf galaxies. We find that the extragalactic GMC properties measured across a wide range of environments, characterized by the Larson relations, are to first order remarkably compatible with those in the Milky Way. Using these data to investigate trends due to galaxy metallicity, we find that: 1) these measurements do not accord with simple expectations from photoionization-regulated star formation theory, 2) there is no trend in the virial CO-to-H2 conversion factor on the spatial scales studied, and 3) there are measurable departures from the Galactic Larson relations in the Small Magellanic Cloud -- the object with the lowest metallicity in the sample -- where GMCs have velocity dispersions that are too small for their sizes. I will discuss the stability of these clouds in the light of our recent far-infrared analysis of this galaxy, and I will contrast the results of the virial and far-infrared studies on the issue of the CO-to-H2 conversion factor and what they tell us about the structure of molecular clouds in primitive galaxies.
[151]  [pdf] - 315137
The Gas Temperature of Starless Cores in Perseus
Comments: Accepted to ApJ; 13 pages, including 3 tables and three figures
Submitted: 2008-09-30
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.
[152]  [pdf] - 14051
The Resolved Properties of Extragalactic Giant Molecular Clouds
Comments: 21 pages, 7 figures, 4 tables (one of them electronic). The Astrophysical Journal, accepted. Revised to reflect changes made to proofs
Submitted: 2008-06-30, last modified: 2008-08-14
We use high spatial resolution observations of CO to systematically measure the resolved size-line width, luminosity-line width, luminosity-size, and the mass-luminosity relations of Giant Molecular Clouds (GMCs) in a variety of extragalactic systems. Although the data are heterogeneous we analyze them in a consistent manner to remove the biases introduced by limited sensitivity and resolution, thus obtaining reliable sizes, velocity dispersions, and luminosities. We compare the results obtained in dwarf galaxies with those from the Local Group spiral galaxies. We find that extragalactic GMC properties measured across a wide range of environments are very much compatible with those in the Galaxy. We use these results to investigate metallicity trends in the cloud average column density and virial CO-to-H2 factor. We find that these measurements do not accord with simple predictions from photoionization-regulated star formation theory, although this could be due to the fact that we do not sample small enough spatial scales or the full gravitational potential of the molecular cloud. We also find that the virial CO-to-H2 conversion factor in CO-bright GMCs is very similar to Galactic, and that the excursions do not show a measurable metallicity trend. We contrast these results with estimates of molecular mass based on far-infrared measurements obtained for the Small Magellanic Cloud, which systematically yield larger masses, and interpret this discrepancy as arising from large H2 envelopes that surround the CO-bright cores. We conclude that GMCs identified on the basis of their CO emission are a unique class of object that exhibit a remarkably uniform set of properties from galaxy to galaxy (abridged).
[153]  [pdf] - 12413
The Mass Distribution and Lifetime of Prestellar Cores in Perseus, Serpens, and Ophiuchus
Comments: 27 pages, 15 figures, 5 tables, accepted for publication in ApJ. Version with full resolution figures available at
Submitted: 2008-05-07
We present an unbiased census of starless cores in Perseus, Serpens, and Ophiuchus, assembled by comparing large-scale Bolocam 1.1 mm continuum emission maps with Spitzer c2d surveys. We use the c2d catalogs to separate 108 starless from 92 protostellar cores in the 1.1 mm core samples from Enoch et al. (2006), Young et al. (2006), and Enoch et al. (2007). A comparison of these populations reveals the initial conditions of the starless cores. Starless cores in Perseus have similar masses but larger sizes and lower densities on average than protostellar cores, with sizes that suggest density profiles substantially flatter than r^-2. By contrast, starless cores in Serpens are compact and have lower masses than protostellar cores; future star formation will likely result in lower mass objects than the currently forming protostars. Comparison to dynamical masses estimated from the NH_3 survey of Perseus cores by Rosolowsky et al. (2007) suggests that most of the starless cores are likely to be gravitationally bound, and thus prestellar. The combined prestellar core mass distribution includes 108 cores and has a slope of -2.3+/-0.4 for M>0.8 Msun. This slope is consistent with recent measurements of the stellar initial mass function, providing further evidence that stellar masses are directly linked to the core formation process. We place a lower limit on the core-to-star efficiency of 25%. There are approximately equal numbers of prestellar and protostellar cores in each cloud, thus the dense prestellar core lifetime must be similar to the lifetime of embedded protostars, or 4.5x10^5 years, with a total uncertainty of a factor of two. Such a short lifetime suggests a dynamic, rather than quasi-static, core evolution scenario, at least at the relatively high mean densities (n>2x10^4 cm^-3) to which we are sensitive.
[154]  [pdf] - 10272
Structural Analysis of Molecular Clouds: Dendrograms
Comments: 15 pages, 16 figures. Accepted to ApJ. Paper will full resolution figures available at
Submitted: 2008-02-20
We demonstrate the utility of dendrograms at representing the essential features of the hierarchical structure of the isosurfaces for molecular line data cubes. The dendrogram of a data cube is an abstraction of the changing topology of the isosurfaces as a function of contour level. The ability to track hierarchical structure over a range of scales makes this analysis philosophically different from local segmentation algorithms like CLUMPFIND. Points in the dendrogram structure correspond to specific volumes in data cubes defined by their bounding isosurfaces. We further refine the technique by measuring the properties associated with each isosurface in the analysis allowing for a multiscale calculation of molecular gas properties. Using COMPLETE 13CO(1-0) data from the L1448 region in Perseus and mock observations of a simulated data cube, we identify regions that have a significant contribution by self-gravity to their energetics on a range of scales. We find evidence for self-gravitation on all spatial scales in L1448 though not in all regions. In the simulated observations, nearly all of the emission is found in objects that would be self-gravitating if gravity were included in the simulation. We reconstruct the size-line width relationship within the data cube using the dendrogram-derived properties and find it follows the standard relation: s_v ~ R^0.58. Finally, we show that constructing the dendrogram of CO J=1-0 emission from the Orion-Monoceros region allows for the identification of giant molecular clouds in a blended molecular line data set using only a physically motivated definition (self-gravitating clouds with masses 5x10^4 Msun.
[155]  [pdf] - 8156
A Metallicity Map of M33
Comments: 4 pages, 2 figures, uses asp2006.sty. To appear in the proceedings of the conference "Formation and Evolution of Galaxy Disks", Rome, Italy, 1-5 October 2007, eds. J. G. Funes and E. M. Corsini
Submitted: 2007-12-14
We present initial results from the M33 Metallicity Project. Out of the thousands of cataloged HII regions in M33, only ~30 have electron-temperature based abundances in the literature. We have obtained Keck spectroscopy of a sample of ~200 HII regions in M33, with 61 detections of the [O III] 4363 A line that can be used for determining electron temperatures, including measurements at small galactocentric radii where auroral lines are generally difficult to detect. We find an oxygen abundance gradient of -0.027 +/- 0.012 dex/kpc, in agreement with infrared measurements of the neon abundance gradient but much shallower than most previous oxygen gradient measurements. There is substantial intrinsic scatter of 0.11 dex in the metallicity at any given radius in M33, which imposes a fundamental limit on the accuracy of gradient measurements that rely on small samples of objects. Finally, we present a two-dimensional map of oxygen abundances across the southern half of M33 and discuss the evidence for deviations from axisymmetry.
[156]  [pdf] - 7449
The M33 Metallicity Project: Resolving the Abundance Gradient Discrepancies in M33
Comments: 10 pages, Accepted to ApJ
Submitted: 2007-11-27
We present a new determination of the metallicity gradient in M33, based on Keck/LRIS measurements of oxygen abundances using the temperature-sensitive emission line [OIII] 4363 A in 61 HII regions. These data approximately triple the sample of direct oxygen abundances in M33. We find a central abundance of 12 + log(O/H) = 8.36+/-0.04 and a slope of -0.027+/-0.012 dex/kpc, in agreement with infrared measurements of the neon abundance gradient but much shallower than most previous oxygen gradient measurements. There is substantial intrinsic scatter of 0.11 dex in the metallicity at any given radius in M33, which imposes a fundamental limit on the accuracy of gradient measurements that rely on small samples of objects. We also show that the ionization state of neon does not follow the ionization state of oxygen as is commonly assumed, suggesting that neon abundance measurements from optical emission lines require careful treatment of the ionization corrections.
[157]  [pdf] - 6632
An Ammonia Spectral Atlas of Dense Cores in Perseus
Comments: 19 pages; Accepted to ApJS; version with high resolution figures available at ; online data at
Submitted: 2007-11-02
We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the (sub)millimeter sources as well as at several positions with no associated continuum emission. For each detection, we have measured physical properties by fitting a simple model to every spectral line simultaneously. Where appropriate, we have refined the model by accounting for low optical depths, multiple components along the line of sight and imperfect coupling to the GBT beam. For the cores in Perseus, we find a typical kinetic temperature of T=11 K, a typical column density of N(NH3)~ 10^14.5 /cm^2 and velocity dispersions ranging from sigma_v = 0.07 km/s to 0.7 km/s. However, many cores with velocity dispersions > 0.2 km/s show evidence for multiple velocity components along the line of sight.
[158]  [pdf] - 1944079
A Catalog of HI Clouds in the Large Magellanic Cloud
Comments: 24 pages, 15 figures, ApJS, in press
Submitted: 2007-06-11, last modified: 2007-06-28
A 21 cm neutral hydrogen interferometric survey of the Large Magellanic Cloud (LMC) combined with the Parkes multi-beam HI single-dish survey clearly shows that the HI gas is distributed in the form of clumps or clouds. The HI clouds and clumps have been identified using a thresholding method with three separate brightness temperature thresholds ($T_b$). Each catalog of HI cloud candidates shows a power law relationship between the sizes and the velocity dispersions of the clouds roughly following the Larson Law scaling $\sigma_v \propto R^{0.5}$, with steeper indices associated with dynamically hot regions. The clouds in each catalog have roughly constant virial parameters as a function mass suggesting that that the clouds are all in roughly the same dynamical state, but the values of the virial parameter are significantly larger than unity showing that turbulent motions dominate gravity in these clouds. The mass distribution of the clouds is a power law with differential indices between -1.6 and -2.0 for the three catalogs. In contrast, the distribution of mean surface densities is a log-normal distribution.
[159]  [pdf] - 1868
Hierarchical Star-Formation in M33: Fundamental properties of the star-forming regions
Comments: 12 pages, 16 figures, accepted MNRAS
Submitted: 2007-06-04
Star-formation within galaxies appears on multiple scales, from spiral structure, to OB associations, to individual star clusters, and often sub-structure within these clusters. This multitude of scales calls for objective methods to find and classify star-forming regions, regardless of spatial size. To this end, we present an analysis of star-forming groups in the local group spiral galaxy M33, based on a new implementation of the Minimum Spanning Tree (MST) method. Unlike previous studies which limited themselves to a single spatial scale, we study star-forming structures from the effective resolution limit (~20pc) to kpc scales. We find evidence for a continuum of star-forming group sizes, from pc to kpc scales. We do not find a characteristic scale for OB associations, unlike that found in previous studies, and we suggest that the appearance of such a scale was caused by spatial resolution and selection effects. The luminosity function of the groups is found to be well represented by a power-law with an index, -2, similar to that found for clusters and GMCs. Additionally, the groups follow a similar mass-radius relation as GMCs. The size distribution of the groups is best described by a log-normal distribution and we show that within a hierarchical distribution, if a scale is selected to find structure, the resulting size distribution will have a log-normal distribution. We find an abrupt drop of the number of groups outside a galactic radius of ~4kpc, suggesting a change in the structure of the star-forming ISM, possibly reflected in the lack of GMCs beyond this radius. (abridged)
[160]  [pdf] - 316893
High Resolution Molecular Gas Maps of M33
Comments: 14 pages, 14 figures, accepted in ApJ. Some tables poorly typeset in emulateapj; see source files for raw data
Submitted: 2007-02-28
New observations of CO (J=1->0) line emission from M33, using the 25 element BEARS focal plane array at the Nobeyama Radio Observatory 45-m telescope, in conjunction with existing maps from the BIMA interferometer and the FCRAO 14-m telescope, give the highest resolution (13'') and most sensitive (RMS ~ 60 mK) maps to date of the distribution of molecular gas in the central 5.5 kpc of the galaxy. A new catalog of giant molecular clouds (GMCs) has a completeness limit of 1.3 X 10^5 M_sun. The fraction of molecular gas found in GMCs is a strong function of radius in the galaxy, declining from 60% in the center to 20% at galactocentric radius R_gal ~ 4 kpc. Beyond that radius, GMCs are nearly absent, although molecular gas exists. Most (90%) of the emission from low mass clouds is found within 100 pc projected separation of a GMC. In an annulus 2.1< R_gal <4.1 kpc, GMC masses follow a power law distribution with index -2.1. Inside that radius, the mass distribution is truncated, and clouds more massive than 8 X 10^5 M_sun are absent. The cloud mass distribution shows no significant difference in the grand design spiral arms versus the interarm region. The CO surface brightness ratio for the arm to interarm regions is 1.5, typical of other flocculent galaxies.
[161]  [pdf] - 85013
Giant Molecular Clouds in M31: I - Molecular Cloud Properties
Comments: Accepted for publication in the Astrophysical Journal. 13 pages, 7 figures
Submitted: 2006-09-14
We present Berkeley Illinois Maryland Association (BIMA) millimeter interferometer observations of giant molecular clouds (GMCs) along a spiral arm in M31. The observations consist of a survey using the compact configuration of the interferometer and follow-up, higher-resolution observations on a subset of the detections in the survey. The data are processed using an analysis algorithm designed to extract GMCs and correct their derived properties for observational biases thereby facilitating comparison with Milky Way data. The algorithm identifies 67 GMCs of which 19 have sufficient signal-to-noise to accurately measure their properties. The GMCs in this portion of M31 are indistinguishable from those found in the Milky Way, having a similar size-line width relationship and distribution of virial parameters, confirming the results of previous, smaller studies. The velocity gradients and angular momenta of the GMCs are comparable to the values measured in M33 and the Milky Way; and, in all cases, are below expected values based on the local galactic shear. The studied region of M31 has a similar interstellar radiation field, metallicity, Toomre Q parameter, and midplane volume density as the inner Milky Way, so the similarity of GMC populations between the two systems is not surprising.
[162]  [pdf] - 82850
The Evolution of the ISM in the Mildly Disturbed Spiral Galaxy NGC 4647
Comments: ApJ, accepted
Submitted: 2006-06-16
We present matched-resolution maps of HI and CO emission in the Virgo Cluster spiral NGC 4647. The galaxy shows a mild kinematic disturbance in which one side of the rotation curve flattens but the other side continues to rise. This kinematic asymmetry is coupled with a dramatic asymmetry in the molecular gas distribution but not in the atomic gas. An analysis of the gas column densities and the interstellar pressure suggests that the H2/HI surface density ratio on the east side of the galaxy is three times higher than expected from the hydrostatic pressure contributed by the mass of the stellar disk. We discuss the probable effects of ram pressure, gravitational interactions, and asymmetric potentials on the interstellar medium and suggest it is likely that a m=1 perturbation in the gravitational potential could be responsible for all of the galaxy's features. Kinematic disturbances of the type seen here are common, but the curious thing about NGC 4647 is that the molecular distribution appears more disturbed than the HI distribution. Thus it is the combination of the two gas phases that provides such interesting insight into the galaxy's history and into models of the interstellar medium.
[163]  [pdf] - 81723
The Role of Pressure in GMC Formation II: The H_2 - Pressure Relation
Comments: 14 pages, 7 figures, Accepted to the Astrophysical Journal
Submitted: 2006-05-01
We show that the ratio of molecular to atomic gas in galaxies is determined by hydrostatic pressure and that the relation between the two is nearly linear. The pressure relation is shown to be good over three orders of magnitude for 14 galaxies including dwarfs, HI-rich, and H_2-rich galaxies as well as the Milky Way. The sample spans a factor of five in mean metallicity. The rms scatter of individual points of the relation is only about a factor of two for all the galaxies, though some show much more scatter than others. Using these results, we propose a modified star formation prescription based on pressure determining the degree to which the ISM is molecular. The formulation is different in high and low pressure regimes defined by whether the gas is primarily atomic or primarily molecular. This formulation can be implemented in simulations and provides a more appropriate treatment of the outer regions of spiral galaxies and molecule-poor systems such as dwarf irregulars and damped Lyman-alpha systems.
[164]  [pdf] - 80176
Giant Molecular Clouds in Local Group Galaxies
Comments: 17 pages, 13 figures, to appear in Protostars and Planets V. A full color version with higher resolution figures is available at
Submitted: 2006-02-27
We present the first comparative study of extragalactic GMCs using complete data sets for entire galaxies and a uniform set of reduction and analysis techniques. We present results based on CO observations for the LMC, SMC, M33, M31, IC10 and the nucleus of M64, and make comparisons with archival Milky Way observations. Our sample includes large spirals and dwarf irregulars with metallicities that vary by an order of magnitude. GMCs in HI rich galaxies are seen to be well-correlated with HI filaments that pervade the galactic disks, suggesting that they form from pre-existing HI structures. Virial estimates of the ratio of CO line strength to H2 column density, XCO, suggests that a value of 4 x 10^20 cm^-2 (K km/s)^-1 is a good value to use in most galaxies (except the SMC) if the GMCs are virialized. However, if the clouds are only marginally self-gravitating, as appears to be the case judging from their appearance, half the virial value may be more appropriate. There is no clear trend of XCO with metallicity. The clouds within a galaxy are shown to have the about the same H2 surface density and differences between galaxies seem to be no more than a factor of ~2. We show that hydrostatic pressure appears to be the main factor in determining what fraction of atomic gas is turned into molecules. In the high-pressure regions often found in galactic centers, the observed properties of GMCs appear to be different from those in the found in the Local Group. From the association of tracers of star formation with GMCs in the LMC, we find that about 1/4 of the GMCs exhibit no evidence of star formation and we estimate that the lifetime of a typical GMC in these galaxies is 20--30 Myr.
[165]  [pdf] - 79564
Bias-free Measurement of Giant Molecular Cloud Properties
Comments: Accepted to PASP (19 pgs., 12 figures). The submission describes an IDL software package available from
Submitted: 2006-01-31
(abridged) We review methods for measuring the sizes, line widths, and luminosities of giant molecular clouds (GMCs) in molecular-line data cubes with low resolution and sensitivity. We find that moment methods are robust and sensitive -- making full use of both position and intensity information -- and we recommend a standard method to measure the position angle, major and minor axis sizes, line width, and luminosity using moment methods. Without corrections for the effects of beam convolution and sensitivity to GMC properties, the resulting properties may be severely biased. This is particularly true for extragalactic observations, where resolution and sensitivity effects often bias measured values by 40% or more. We correct for finite spatial and spectral resolutions with a simple deconvolution and we correct for sensitivity biases by extrapolating properties of a GMC to those we would expect to measure with perfect sensitivity. The resulting method recovers the properties of a GMC to within 10% over a large range of resolutions and sensitivities, provided the clouds are marginally resolved with a peak signal-to-noise ratio greater than 10. We note that interferometers systematically underestimate cloud properties, particularly the flux from a cloud. The degree of bias depends on the sensitivity of the observations and the (u,v) coverage of the observations. In the Appendix to the paper we present a conservative, new decomposition algorithm for identifying GMCs in molecular-line observations. This algorithm treats the data in physical rather than observational units, does not produce spurious clouds in the presence of noise, and is sensitive to a range of morphologies. As a result, the output of this decomposition should be directly comparable among disparate data sets.
[166]  [pdf] - 75568
The Mass Spectra of Giant Molecular Clouds in the Local Group
Comments: 8 pages, 3 figures, accepted to PASP
Submitted: 2005-08-31
We reanalyze the catalogs of molecular clouds in the Local Group to determine the parameters of their mass distributions in a uniform manner. The analysis uses the error-in-variables method of parameter estimation which accounts not only for the variance of the sample when drawn from a parent distribution but also for errors in the mass measurements. Testing the method shows that it recovers the underlying properties of cumulative mass distribution without bias while accurately reflecting uncertainties in the parameters. Clouds in the inner disk of the Milky Way follow a truncated power-law distribution with index \gamma=-1.5 +/- 0.1 and maximum mass of 10^6.5 M_sun. The distributions of cloud mass for the outer Milky Way and M33 show significantly steeper indices (\gamma_OMW= -2.1 +/- 0.2 and \gamma_M33= -2.9 +/- 0.4) with no evidence of a cutoff. The mass distribution of clouds in the Large Magellanic Cloud has a marginally steeper distribution than the inner disk of the Milky Way (\gamma=-1.7 +/- 0.2) and also shows evidence of a truncation with a maximum mass of 10^6.5 M_sun. The mass distributions of molecular clouds vary dramatically across the Local Group, even after accounting for the systematic errors that arise in comparing heterogeneous data and catalogs. These differences should be accounted for in studies that aim to reproduce the molecular cloud mass distributions or in studies that use the mass spectrum as a parameter in a model.
[167]  [pdf] - 1468735
Giant Molecular Clouds in M64
Comments: Accepted for publication in ApJ, 22 pages, 21 figures
Submitted: 2005-01-18
(abridged) We investigate the properties of Giant Molecular Clouds (GMCs) in the molecule-rich galaxy M64 (NGC 4826). In M64, the mean surface density of molecular gas is 2N(H_2) ~ 10^22 over a 2 kpc region, equal to the surface densities of individual GMCs in the Milky Way. We observed the J=1->0 transitions of CO, 13CO, and HCN. The line ratio W_CO/W_13CO for 200 pc < R_gal < 800 pc is comparable to that found in the Milky Way and increases significantly outside this region, in part due to a large contribution to the CO emission from diffuse gas, which composes 25% of the molecular mass in the galaxy. We developed a modified CLUMPFIND algorithm to decompose the 13CO emission into 25 resolved clouds. The clouds have a luminosity--linewidth relationship L ~ DV^{2.2+/-0.4}, substantially different from the Milky Way trend reported by Solomon et al. (1987): L ~ DV^5. Similarly, the clouds have a linewidth--size relationship of DV ~ R_e^{1.0 +/- 0.3} compared to DV ~ R_e^0.5 in the Milky Way. Estimates of the kinetic and binding energies of the clouds suggest that the clouds are self-gravitating and significantly overpressured with respect to the remainder of the ISM in M64. The M64 clouds have a mean surface density at least 2.5 times larger than observed in Local Group GMCs, and the surface density is not independent of mass as it is in the Local Group: \Sigma_H2 ~ M^{0.7+/-0.2}. The clouds are correlated with the recombination emission from the galaxy, implying that they are star forming; the rate is comparable to that in other galaxies despite the increased densities of the clouds. We note that the internal pressures of clouds in several galaxies scales with the external pressure exerted on the clouds by the ambient ISM: P_int ~ P_ext^{0.75+/-0.05}.
[168]  [pdf] - 69119
The Initial Mass Function of GMCs
Comments: 9 pages, 5 figures, to appear the proceedings of "IMF@50: The Initial Mass Function 50 Years Later" (eds. E. Corbelli, F. Palla, and H. Zinnecker)
Submitted: 2004-11-17
The properties of GMCs in several Local Group galaxies are quantified and compared. It is found that the mass spectrum of GMCs varies from galaxy to galaxy. The variations are significant and do not appear to be the result of systematic uncertainties. Nevertheless, it appears that all of the GMCs follow the same size-linewidth and mass-linewidth relations with little scatter. The power law indices of these relations imply that the GMCs are self-gravitating, and that the mean surface density of Local Group GMCs is approximately constant. This, in turn, implies that the mean internal pressure of GMCs is also constant. If the IMF of stars is determined by a Jeans instability, this constant internal pressure suggests that the distribution of stellar masses does not vary significantly in galactic disks when averaged over suitably large areas. Thus, although the distribution of GMC masses produced by various Local Group galaxies is quite variable, the large-scale properties of the GMCs is not.
[169]  [pdf] - 66754
GMC Formation in M33: the Role of Pressure
Comments: 8 Pages, to appear in Star Formation in the Interstellar Medium Workshop, ASP Conference Series (F. C. Adams, D. Johnstone, D. N. C. Lin, & E. Ostriker, eds.)
Submitted: 2004-08-13
We present a complete map of the Giant Molecular Clouds in M33 and show that they lie exclusively on filaments of HI. The GMCs are localized on the extended HI filaments, and HI filaments extend well beyond the radius where the GMCs are found, suggesting that the H_2 forms from the HI and not vice-versa. We propose that the fraction of neutral gas that is molecular at a particular radius is determined primarily by the ambient hydrostatic pressure at that radius. This leads to a prediction that the transition radius, the radius where the mean surface density of HI and H_2 are equal, occurs at a constant value of the stellar surface density. We show that for a sample of 26 galaxies, the stellar surface density at the transition radius is constant to 40%.
[170]  [pdf] - 66750
A Statistical Investigation of HI in the Magellanic Bridge
Comments: 15 pages, Latex. using emulateapj.sty. 11 Figures. Accepted for publication in ApJ
Submitted: 2004-08-13
We present results from two statistical analyses applied to an neutral hydrogen (HI) dataset of the nearby tidal bridge in the Magellanic System. Primarily, analyses of the Spatial Power Spectrum suggest that the Magellanic Bridge, historically considered to be a single contiguous feature, may in fact be a projection of two kinematically and morphologically distinct structures. The southern and more obviously turbulent parts appear to show structure organized similarly to the adjacent Small Magellanic Cloud (SMC), while the northern regions are shown to be relatively deficient in large scale power. The extent of modification to the spatial power index by the velocity fluctuations is also highly variant across these parts of the Bridge. We find again that the northern part appears distinct from the southern part and from the SMC, in that the power spectrum is significantly more affected by slower velocity perturbations. We also probe the rate of spectral variation of the HI by measuring the Spectral Correlation Function over selected regions. The results from this analysis highlight a tendency for the HI spectra within the bright parts of the Bridge to have a more persistent correlation in the E-W direction than in the N-S direction. These results are considered to be quantitative evidence for the tidal processes which are thought to have been active throughout the evolution of the Magellanic Bridge
[171]  [pdf] - 66334
The Role of Pressure in GMC Formation
Comments: 4 pages, 1 Figure, Accepted to ApJ Letters
Submitted: 2004-07-22
We examine the hypothesis that hydrostatic pressure alone determines the ratio of atomic to molecular gas averaged over a particular radius in disk galaxies. The hypothesis implies that the transition radius, the location where the ratio is unity, should always occur at the same value of stellar surface density in all galaxies. We examine data for 28 galaxies and find that the stellar surface density at the transition radius is indeed constant to 40% at a value of 120 M_sun/pc^2. If the hypothesis can be confirmed at all radii within a large range of galaxy types and metallicities, combining it with the observed constancy of the star formation rate with H_2 surface density may enable a physically motivated star formation prescription with wide applicability.
[172]  [pdf] - 58035
Giant Molecular Clouds in M33 -- II. High Resolution Observations
Comments: 16 Pages, 16 Figures, Accepted for Publication in the Astrophysical Journal
Submitted: 2003-07-16, last modified: 2003-10-15
We present 12CO(1->0) observations of 45 giant molecular clouds in M33 made with the BIMA array. The observations have a linear resolution of 20 pc, sufficient to measure the sizes of most GMCs in the sample. We place upper limits on the specific angular momentum of the GMCs and find the observed values to be nearly an order of magnitude below the values predicted from simple formation mechanisms. The velocity gradients across neighboring, high-mass GMCs appear preferentially aligned on scales less than 500 pc. If the clouds are rotating, 40% are counter-rotating with respect to the galaxy. GMCs require a braking mechanism if they form from the large scale radial accumulation of gas. These observations suggest that molecular clouds form locally out of atomic gas with significant braking by magnetic fields to dissipate the angular momentum imparted by galactic shear. The observed GMCs share basic properties with those found in the Galaxy such as similar masses, sizes, and linewidths as well as a constant surface density of 120 M_sun pc^{-2}. The size-linewidth relationship follows Delta V ~ r^{0.45\pm 0.02}, consistent with that found in the Galaxy. The cloud virial masses imply that the CO-to-H_2 conversion factor has a value of 2 X 10^20 H2 cm^-2/(K km/s) and does not change significantly over the disk of M33 despite a change of 0.8 dex in the metallicity.
[173]  [pdf] - 58667
Giant Molecular Clouds in M33 - I. BIMA All Disk Survey
Comments: 19 pages, Accepted for Publication in the Astrophysical Journal Supplement
Submitted: 2003-08-21
We present the first interferometric CO(J=1->0) map of the entire H-alpha disk of M33. The 13" diameter synthesized beam corresponds to a linear resolution of 50 pc, sufficient to distinguish individual giant molecular clouds (GMCs). From these data we generated a catalog of 148 GMCs with an expectation that no more than 15 of the sources are spurious. The catalog is complete down to GMC masses of 1.5 X 10^5 M_sun and contains a total mass of 2.3 X 10^7 M_sun. Single dish observations of CO in selected fields imply that our survey detects ~50% of the CO flux, hence that the total molecular mass of M33 is 4.5 X 10^7 M_sun, approximately 2% of the HI mass. The GMCs in our catalog are confined largely to the central region (R < 4 kpc). They show a remarkable spatial and kinematic correlation with overdense HI filaments; the geometry suggests that the formation of GMCs follows that of the filaments. The GMCs exhibit a mass spectrum dN/dM ~ M^(-2.6 +/- 0.3), considerably steeper than that found in the Milky Way and in the LMC. Combined with the total mass, this steep function implies that the GMCs in M33 form with a characteristic mass of 7 X 10^4 M_sun. More than 2/3 of the GMCs have associated HII regions, implying that the GMCs have a short quiescent period. Our results suggest the rapid assembly of molecular clouds from atomic gas, with prompt onset of massive star formation.
[174]  [pdf] - 38691
The Effects of Noise and Sampling on the Spectral Correlation Function
Comments: 28 pages, 12 figures included. To be published in The Astrophysical Journal
Submitted: 2000-10-17
The effects of noise and sampling on the ``Spectral Correlation Function'' (SCF) introduced by Rosolowsky et al. 1999 are studied using observational data, numerical simulations of magneto-hydrodynamic turbulence, and simple models of Gaussian spectral line profiles. The most significant innovations of this paper are: i) the normalization of the SCF based on an analytic model for the effect of noise; ii) the computation of the SCF as a function of the spatial lag between spectra within a map. By computing of the dependence of the SCF on the spatial lag, S_o(dr), we have been able to conclude that: 1) S_o(dr) is a power law, with slope in the range of scales l_i-l_o. 2) The correlation outer scale, l_o, is determined by the size of the map, and no evidence for a true departure from self-similarity on large scales has been found. 3) The correlation inner scale, l_i, is a true estimate of the smallest self-similar scale in a map. 4) The spectral slope in a given region, is independent of velocity resolution (above a minimum resolution threshold), spatial resolution, and average spectrum quality. 5) Molecular transitions which trace higher gas density yield larger values of the spectral slope than transitions tracing lower gas density. 6) Nyquist sampling, bad pixels in detector arrays, and reference sharing data acquisition need to be taken into account for a correct determination of the SCF at lag dr=1. The value of the spectral slope, however, can be computed correctly without a detailed knowledge of observational procedures.
[175]  [pdf] - 105858
The Spectral Correlation Function -- A New Tool for Analyzing Spectral-Line Maps
Comments: 29 pages, including 4 figures (tar file submitted as source) See also:
Submitted: 1999-03-30, last modified: 1999-04-20
The "spectral correlation function" analysis we introduce in this paper is a new tool for analyzing spectral-line data cubes. Our initial tests, carried out on a suite of observed and simulated data cubes, indicate that the spectral correlation function [SCF] is likely to be a more discriminating statistic than other statistical methods normally applied. The SCF is a measure of similarity between neighboring spectra in the data cube. When the SCF is used to compare a data cube consisting of spectral-line observations of the ISM with a data cube derived from MHD simulations of molecular clouds, it can find differences that are not found by other analyses. The initial results presented here suggest that the inclusion of self-gravity in numerical simulations is critical for reproducing the correlation behavior of spectra in star-forming molecular clouds.