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

Bolatto, Alberto

Normalized to: Bolatto, A.

210 article(s) in total. 865 co-authors, from 1 to 94 common article(s). Median position in authors list is 4,0.

[1]  [pdf] - 2068910
Plateau de Bure High-z Blue-Sequence Survey 2 (PHIBSS2): Search for Secondary Sources, CO Luminosity Functions in the Field, and the Evolution of Molecular Gas Density through Cosmic Time
Comments: 29 pages (including 9 appendix pages), 14 figures and 6 tables, submitted to ApJ, replaced with version accepted by ApJ
Submitted: 2019-08-05, last modified: 2020-03-23
We report on the results of a search for serendipitous sources in CO emission in 110 cubes targeting CO(2-1), CO(3-2), and CO(6-5) at z ~ 1-2 from the second Plateau de Bure High-z Blue-Sequence Survey (PHIBSS2). The PHIBSS2 observations were part of a 4-year legacy program at the IRAM Plateau de Bure Interferometer aimed at studying early galaxy evolution from the perspective of molecular gas reservoirs. We present a catalog of 67 candidate secondary sources from this search, with 45 out of the 110 data cubes showing sources in addition to the primary target that appear to be field detections, unrelated to the central sources. This catalog includes the redshifts, line widths, fluxes, as well as an estimation of their reliability based on their false positive probability. We perform a search in the 3D-HST/CANDELS catalogs for the secondary CO detections and tentatively find that ~64% of these have optical counterparts, which we use to constrain their redshifts. Finally, we use our catalog of candidate CO detections to derive the CO(2-1), CO(3-2), CO(4-3), CO(5-4), and CO(6-5) luminosity functions over a range of redshifts, as well as the molecular gas mass density evolution. Despite the different methodology, these results are in very good agreement with previous observational constraints derived from blind searches in deep fields. They provide an example of the type of "deep field" science that can be carried out with targeted observations.
[2]  [pdf] - 2052672
Dynamical Equilibrium in the Molecular ISM in 28 Nearby Star-Forming Galaxies
Comments: 28 pages + 3 appendices, ApJ in press. See for a short video describing the main results. Data tables available at prior to publication
Submitted: 2020-02-20
We compare the observed turbulent pressure in molecular gas, $P_\mathrm{turb}$, to the required pressure for the interstellar gas to stay in equilibrium in the gravitational potential of a galaxy, $P_\mathrm{DE}$. To do this, we combine arcsecond resolution CO data from PHANGS-ALMA with multi-wavelength data that traces the atomic gas, stellar structure, and star formation rate (SFR) for 28 nearby star-forming galaxies. We find that $P_\mathrm{turb}$ correlates with, but almost always exceeds the estimated $P_\mathrm{DE}$ on kiloparsec scales. This indicates that the molecular gas is over-pressurized relative to the large-scale environment. We show that this over-pressurization can be explained by the clumpy nature of molecular gas; a revised estimate of $P_\mathrm{DE}$ on cloud scales, which accounts for molecular gas self-gravity, external gravity, and ambient pressure, agrees well with the observed $P_\mathrm{turb}$ in galaxy disks. We also find that molecular gas with cloud-scale ${P_\mathrm{turb}}\approx{P_\mathrm{DE}}\gtrsim{10^5\,k_\mathrm{B}\,\mathrm{K\,cm^{-3}}}$ in our sample is more likely to be self-gravitating, whereas gas at lower pressure appears more influenced by ambient pressure and/or external gravity. Furthermore, we show that the ratio between $P_\mathrm{turb}$ and the observed SFR surface density, $\Sigma_\mathrm{SFR}$, is compatible with stellar feedback-driven momentum injection in most cases, while a subset of the regions may show evidence of turbulence driven by additional sources. The correlation between $\Sigma_\mathrm{SFR}$ and kpc-scale $P_\mathrm{DE}$ in galaxy disks is consistent with the expectation from self-regulated star formation models. Finally, we confirm the empirical correlation between molecular-to-atomic gas ratio and kpc-scale $P_\mathrm{DE}$ reported in previous works.
[3]  [pdf] - 2050535
Cool Outflows in Galaxies and their Implications
Comments: Review accepted for publication in Astronomy & Astrophysics Reviews, 168 pages, 34 figures
Submitted: 2020-02-18
Neutral-atomic and molecular outflows are a common occurrence in galaxies, near and far. They operate over the full extent of their galaxy hosts, from the innermost regions of galactic nuclei to the outermost reaches of galaxy halos. They carry a substantial amount of material that would otherwise have been used to form new stars. These cool outflows may have a profound impact on the evolution of their host galaxies and environments. This article provides an overview of the basic physics of cool outflows, a comprehensive assessment of the observational techniques and diagnostic tools used to characterize them, a detailed description of the best-studied cases, and a more general discussion of the statistical properties of these outflows in the local and distant universe. The remaining outstanding issues that have not yet been resolved are summarized at the end of the review to inspire new research directions.
[4]  [pdf] - 2050303
The Experiment for Cryogenic Large-aperture Intensity Mapping (EXCLAIM)
Comments: 10 pages, 4 figures
Submitted: 2019-12-15
The EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM) is a cryogenic balloon-borne instrument that will survey galaxy and star formation history over cosmological time scales. Rather than identifying individual objects, EXCLAIM will be a pathfinder to demonstrate an intensity mapping approach, which measures the cumulative redshifted line emission. EXCLAIM will operate at 420-540 GHz with a spectral resolution R=512 to measure the integrated CO and [CII] in redshift windows spanning 0 < z < 3.5. CO and [CII] line emissions are key tracers of the gas phases in the interstellar medium involved in star-formation processes. EXCLAIM will shed light on questions such as why the star formation rate declines at z < 2, despite continued clustering of the dark matter. The instrument will employ an array of six superconducting integrated grating-analog spectrometers (micro-spec) coupled to microwave kinetic inductance detectors (MKIDs). Here we present an overview of the EXCLAIM instrument design and status.
[5]  [pdf] - 2046390
Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: The KINGFISH Sample
Comments: ApJ, accepted
Submitted: 2019-12-10
Dust and starlight are modeled for the KINGFISH project galaxies. With data from 3.6 micron to 500 micron, models are strongly constrained. For each pixel in each galaxy we estimate (1) dust surface density; (2) q_PAH, the dust mass fraction in PAHs; (3) distribution of starlight intensities heating the dust; (4) luminosity emitted by the dust; and (5) dust luminosity from regions with high starlight intensity. The models successfully reproduce both global and resolved spectral energy distributions. We provide well-resolved maps for the dust properties. As in previous studies, we find q_PAH to be an increasing function of metallicity, above a threshold Z/Z_sol approx 0.15. Dust masses are obtained by summing the dust mass over the map pixels; these "resolved" dust masses are consistent with the masses inferred from model fits to the global photometry. The global dust-to-gas ratios obtained from this study correlate with galaxy metallicities. Systems with Z/Z_sol > 0.5 have most of their refractory elements locked up in dust, whereas when Z/Z_sol < 0.3 most of these elements tend to remain in the gas phase. Within galaxies, we find that q_PAH is suppressed in regions with unusually warm dust with nu L_nu(70 um) > 0.4L_dust. With knowledge of one long-wavelength flux density ratio (e.g., f_{160}/f_{500}), the minimum starlight intensity heating the dust (U_min) can be estimated to within ~50%. For the adopted dust model, dust masses can be estimated to within ~0.07 dex accuracy using the 500 micron luminosity nu L_nu(500) alone. There are additional systematic errors arising from the choice of dust model, but these are hard to estimate. These calibrated prescriptions may be useful for studies of high-redshift galaxies.
[6]  [pdf] - 2026414
The EDGE-CALIFA Survey: Using Optical Extinction to Probe the Spatially-Resolved Distribution of Gas in Nearby Galaxies
Comments: 12 pages, 10 figures. Submitted to MNRAS
Submitted: 2019-11-21
We present an empirical relation between the cold gas surface density ($\Sigma_{\rm gas}$) and the optical extinction (${\rm A_V}$) in a sample of 103 galaxies from the Extragalactic Database for Galaxy Evolution (EDGE) survey. This survey provides CARMA interferometric CO observations for 126 galaxies included in the Calar Alto Legacy Integral Field Area (CALIFA) survey. The matched, spatially resolved nature of these data sets allows us to derive the $\Sigma_{\rm gas}$-${\rm A_V}$ relation on global, radial, and kpc (spaxel) scales. We determine ${\rm A_V}$ from the Balmer decrement (H$\alpha$/H$\beta$). We find that the best fit for this relation is $\Sigma_{\rm gas} ({\rm M_\odot pc^{-2}})\sim~26~\times~ {\rm A_V}({\rm mag})$, and that it does not depend on the spatial scale used for the fit. However, the scatter in the fits increases as we probe smaller spatial scales, reflecting the complex relative spatial distributions of stars, gas, and dust. We investigate the $\Sigma_{\rm gas}$/ ${\rm A_V}$ ratio on radial and spaxel scales as a function of ${\rm EW(H\alpha)}$. We find that at larger values of ${\rm EW(H\alpha)}$ (i.e., actively star-forming regions) this ratio tend to converge to the value expected for dust-star mixed geometries ($\sim$ 30 $\mathrm{M_{\odot} \,pc^{-2}\,mag^{-1}}$). On radial scales, we do not find a significant relation between the $\Sigma_{\rm gas}$/${\rm A_V}$ ratio and the ionized gas metallicity. We contrast our estimates of $\Sigma_{\rm gas}$ using ${\rm A_V}$ with compilations in the literature of the gas fraction on global and radial scales as well as with well known scaling relations such as the radial star-formation law and the $\Sigma_{\rm gas}$-$\Sigma_*$ relation. These tests show that optical extinction is a reliable proxy for estimating $\Sigma_{\rm gas}$ in the absence of direct sub/millimeter observations of the cold gas.
[7]  [pdf] - 2026349
The Dark Matter Distributions in Low-Mass Disk Galaxies. II. The Inner Density Profiles
Submitted: 2019-11-13
Dark matter-only simulations predict that dark matter halos have steep, cuspy inner density profiles, while observations of dwarf galaxies find a range of inner slopes that are often much shallower. There is debate whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. In Paper 1 of this series, we obtained high-resolution integral field H$\alpha$ observations for 26 dwarf galaxies with $M_*=10^{8.1}-10^{9.7}\textrm{M}_\odot$. We derived rotation curves from our observations, which we use here to construct mass models. We model the total mass distribution as the sum of a generalized Navarro-Frenk-White (NFW) dark matter halo and the stellar and gaseous components. Our analysis of the slope of the dark matter density profile focuses on the inner 300-800 pc, chosen based on the resolution of our data and the region resolved by modern hydrodynamical simulations. The inner slope measured using ionized and molecular gas tracers is consistent, and it is additionally robust to the choice of stellar mass-to-light ratio. We find a range of dark matter profiles, including both cored and cuspy slopes, with an average of $\rho_{\rm DM}\sim r^{-0.74\pm 0.07}$, shallower than the NFW profile, but steeper than those typically observed for lower-mass galaxies with $M_*\sim 10^{7.5}\textrm{M}_\odot$. Simulations that reproduce the observed slopes in those lower-mass galaxies also produce slopes that are too shallow for galaxies in our mass range. We therefore conclude that supernova feedback models do not yet provide a fully satisfactory explanation for the observed trend in dark matter slopes.
[8]  [pdf] - 2026307
Mapping Spatial Variations of HI Turbulent Properties in the Small and Large Magellanic Cloud
Comments: 17 pages, 14 figures, accepted to the ApJ
Submitted: 2019-11-11
We developed methods for mapping spatial variations of the spatial power spectrum (SPS) and structure function (SF) slopes, with a goal of connecting neutral hydrogen (HI) statistical properties with the turbulent drivers. The new methods were applied on the HI observations of the Small and Large Magellanic Clouds (SMC and LMC). In the case of the SMC, we find highly uniform turbulent properties of HI, with no evidence for local enhancements of turbulence due to stellar feedback. Such properties could be caused by a significant turbulent driving on large-scales. Alternatively, a significant line-of-sight depth of the SMC could be masking out localized regions with a steeper SPS slope caused by stellar feedback. In contrast to the SMC, the LMC HI shows a large diversity in terms of its turbulent properties. Across most of the LMC, the small-scale SPS slope is steeper than the large-scale slope due to the presence of the HI disk. On small spatial scales, we find several areas of localized steepening of the SPS slope around major HII regions, with the 30 Doradus region being the most prominent. This is in agreement with predictions from numerical simulations which suggest steepening of the SPS slope due to stellar feedback eroding and destroying interstellar clouds. We also find localized steepening of the large-scale SPS slope in the outskirts of the LMC. This is likely caused by the flaring of the HI disk, or alternatively ram-pressure stripping of the LMC disk due to the interactions with the surrounding halo gas.
[9]  [pdf] - 2025997
Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator
Comments: 19 pages, 8 figures, 8 tables
Submitted: 2019-10-11
The brightest observed emission line in many star-forming galaxies is the [CII] 158 micron line, making it detectable up to z~7. In order to better understand and quantify the [CII] emission as a tracer of star-formation, the theoretical ratio between the [NII] 205 micron emission and the [CII] 158 micron emission has been employed to empirically determine the fraction of [CII] emission that originates from the ionized and neutral phases of the ISM. Sub-kiloparsec measurements of the [CII] 158 micron and [NII] 205 micron line in nearby galaxies have recently become available as part of the Key Insights in Nearby Galaxies: a Far Infrared Survey with Herschel (KINGFISH) and Beyond the Peak (BtP) programs. With the information from these two far-infrared lines along with the multi-wavelength suite of KINGFISH data, a calibration of the [CII] emission line as a star formation rate indicator and a better understanding of the [CII] deficit are pursued. [CII] emission is also compared to PAH emission in these regions to compare photoelectric heating from PAH molecules to cooling by [CII] in the neutral and ionized phases of the ISM. We find that the [CII] emission originating in the neutral phase of the ISM does not exhibit a deficit with respect to the infrared luminosity and is therefore preferred over the [CII] emission originating in the ionized phase of the ISM as a star formation rate indicator for the normal star-forming galaxies included in this sample.
[10]  [pdf] - 2005508
A Search for Intrinsic HI 21 cm and OH 18 cm Absorption Toward Compact Radio Sources
Comments: 28 pages, 15 figures. Accepted for publication in ApJS
Submitted: 2019-09-26
We present the results of a large search for intrinsic HI 21 cm and OH 18 cm absorption in 145 compact radio sources in the redshift range 0.02< z <3.8 with the Green Bank Telescope. We re-detect HI 21 cm absorption toward six known absorption systems but detect no new HI or OH absorption in 102 interference-free sources. 79 sources have not previously been observed for HI 21 cm absorption. We recover a mean optical depth limit of $\tau_{3\sigma}<0.023$ for all the non-detections in the survey. Our results do not support the high intrinsic absorption rates found by previous studies in compact radio sources at low redshift. Our results do, however, support the hypothesis proposed by Curran et al. (2008) that high ultraviolet (UV) luminosity active galactic nuclei (AGN) do not show intrinsic HI 21 cm absorption, confirming a threshold of $L_{\rm UV} = 10^{23}$ W Hz$^{-1}$, above which our intrinsic absorption fraction is zero (54 sources). The exact nature of the UV luminosity effect on HI absorption systems remains ambiguous. We additionally find no statistical correlation between the 1.4 GHz radio luminosity or the source size and the 21~cm absorption detection rate. We attribute the lack of intrinsic absorption in our survey to the UV luminosity effect caused by an optical selection bias and a decreased column density sensitivity with increasing redshift due to lower radio continuum flux densities, high radio frequency interference, and higher telescope system temperatures at low frequencies.
[11]  [pdf] - 1958257
An ATCA Survey of HI Absorption in the Magellanic Clouds I: HI Gas Temperature Measurements in the Small Magellanic Cloud
Comments: 22 pages, 16 figures, accepted for publication in ApJS
Submitted: 2019-08-13
We present the first results from the Small Magellanic Cloud portion of a new Australia Telescope Compact Array (ATCA) HI absorption survey of both of the Magellanic Clouds, comprising over 800 hours of observations. Our new HI absorption line data allow us to measure the temperature and fraction of cold neutral gas in a low metallicity environment. We observed 22 separate fields, targeting a total of 55 continuum sources against 37 of which we detected HI absorption; from this we measure a column density weighted mean average spin temperature of $<T_{s}>=150$ K. Splitting the spectra into individual absorption line features, we estimate the temperatures of different gas components and find an average cold gas temperature of $\sim{30}$ K for this sample, lower than the average of $\sim{40}$ K in the Milky Way. The HI appears to be evenly distributed throughout the SMC and we detect absorption in $67\%$ of the lines of sight in our sample, including some outside the main body of the galaxy ($N_{\text{HI}}>2\times{10^{21}}$ cm$^{-2}$). The optical depth and temperature of the cold neutral atomic gas shows no strong trend with location spatially or in velocity. Despite the low metallicity environment, we find an average cold gas fraction of $\sim{20\%}$, not dissimilar from that of the Milky Way.
[12]  [pdf] - 1945845
The molecular outflow in NGC253 at a resolution of two parsecs
Comments: 27 pages, 12 figures, accepted for publication in ApJ
Submitted: 2019-07-01
We present 0.15'' (~2.5pc) resolution ALMA CO(3-2) observations of the starbursting center in NGC253. Together with archival ALMA CO(1-0) and CO(2-1) data we decompose the emission into a disk and non-disk component. We find ~7-16% of the CO luminosity to be associated with the non-disk component ($1.2-4.2 \times 10^7$ K km s$^{-1}$ pc$^2$). The total molecular gas mass in the center of NGC253 is $\sim 3.6 \times 10^8$ M$_\odot$ with $\sim 0.5 \times 10^8$ M$_\odot$ (~15%) in the non-disk component. These measurements are consistent across independent mass estimates through three CO transitions. The high-resolution CO(3-2) observations allow us to identify the molecular outflow within the non-disk gas. Using a starburst conversion factor, we estimate the deprojected molecular mass outflow rate, kinetic energy and momentum in the starburst of NGC253. The deprojected molecular mass outflow rate is in the range ~14-39 M$_\odot$ yr$^{-1}$ with an uncertainty of 0.4dex. The large spread arises due to different interpretations of the kinematics of the observed gas while the errors are due to unknown geometry. The majority of this outflow rate is contributed by distinct outflows perpendicular to the disk, with a significant contribution by diffuse molecular gas. This results in a mass loading factor $\eta = \dot{M}_\mathrm{out} / \dot{M}_\mathrm{SFR}$ in the range $\eta \sim 8-20$ for gas ejected out to ~300pc. We find the kinetic energy of the outflow to be $\sim 2.5-4.5 \times 10^{54}$ erg and ~0.8dex typical error which is ~0.1% of the total or ~8% of the kinetic energy supplied by the starburst. The outflow momentum is $4.8-8.7 \times 10^8$ M$_\odot$ km s$^{-1}$ (~0.5dex error) or ~2.5-4% of the kinetic momentum released into the ISM by feedback. The unknown outflow geometry and launching sites are the primary source of uncertainty in this study.
[13]  [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.
[14]  [pdf] - 1938490
EMPIRE: The IRAM 30-m Dense Gas Survey of Nearby Galaxies
Comments: Accepted for publication in ApJ
Submitted: 2019-06-20
We present EMPIRE, an IRAM 30-m large program that mapped $\lambda = 3{-}4$ mm dense gas tracers at $\sim 1{-}2\,$kpc resolution across the whole star-forming disk of nine nearby, massive, spiral galaxies. We describe the EMPIRE observing and reduction strategies and show new whole-galaxy maps of HCN(1-0), HCO$^+$(1-0), HNC(1-0) and CO(1-0). We explore how the HCN-to-CO and IR-to-HCN ratios, observational proxies for the dense gas fraction and dense gas star formation efficiency, depend on host galaxy and local environment. We find that the fraction of dense gas correlates with stellar surface density, gas surface density, molecular-to-atomic gas ratio, and dynamical equilibrium pressure. In EMPIRE, the star formation rate per unit dense gas anti-correlates with these same environmental parameters. Thus, although dense gas appears abundant the central regions of many spiral galaxies, this gas appears relatively inefficient at forming stars. These results qualitatively agree with previous work on nearby galaxies and the Milky Way's Central Molecular Zone. To first order, EMPIRE demonstrates that the conditions in a galaxy disk set the gas density distribution and that the dense gas traced by HCN shows an environment-dependent relation to star formation. However, our results also show significant ($\pm 0.2$ dex) galaxy-to-galaxy variations. We suggest that gas structure below the scale of our observations and dynamical effects likely also play an important role.
[15]  [pdf] - 1958132
The EDGE-CALIFA Survey: Evidence for Pervasive Extraplanar Diffuse Ionized Gas in Nearby Edge-On Galaxies
Comments: 36 pages (including 13 appendix pages), 17 figures, resubmitted to ApJ after one round of refereeing
Submitted: 2019-05-13
We investigate the prevalence, properties, and kinematics of extraplanar diffuse ionized gas (eDIG) in a sample of 25 edge-on galaxies selected from the CALIFA survey. We measure ionized gas scale heights from ${\rm H\alpha}$ and find that 90% have measurable scale heights with a median of $0.8^{+0.7}_{-0.4}$ kpc. From the ${\rm H\alpha}$ kinematics, we find that 60% of galaxies show a decrease in the rotation velocity as a function of height above the midplane. This lag is characteristic of eDIG, and we measure a median lag of 21 km s$^{-1}$ kpc$^{-1}$ which is comparable to lags measured in the literature. We also investigate variations in the lag with radius. $\rm H{\small I}$ lags have been reported to systematically decrease with galactocentric radius. We find both increasing and decreasing ionized gas lags with radius, as well as a large number of galaxies consistent with no radial lag variation, and investigate these results in the context of internal and external origins for the lagging ionized gas. We confirm that the ${\rm [S{\small II}]}$/${\rm H\alpha}$ and ${\rm [N{\small II}]}$/${\rm H\alpha}$ line ratios increase with height above the midplane as is characteristic of eDIG. The ionization of the eDIG is dominated by star-forming complexes (leaky ${\rm H{\small II}}$ regions). We conclude that the lagging ionized gas is turbulent ejected gas likely resulting from star formation activity in the disk as opposed to gas in the stellar thick disk or bulge. This is further evidence for the eDIG being a product of stellar feedback and for the pervasiveness of this WIM-like phase in many local star-forming galaxies.
[16]  [pdf] - 1861434
Astro2020 Science White Paper: Cold Gas Outflows, Feedback, and the Shaping of Galaxies
Comments: Submitted as a Science White Paper to the Astro2020 Decadal Survey. arXiv admin note: substantial text overlap with arXiv:1810.06737
Submitted: 2019-04-03
There is wide consensus that galaxy outflows are one of the most important processes determining the evolution of galaxies through cosmic time, for example playing a key role in shaping the galaxy mass function. Our understanding of outflows and their drivers, however, is in its infancy --- this is particularly true for the cold (neutral atomic and molecular) phases of outflows, which present observational and modeling challenges. Here we outline several key open questions, briefly discussing the requirements of the observations necessary to make progress, and the relevance of several existing and planned facilities. It is clear that galaxy outflows, and particularly cold outflows, will remain a topic of active research for the next decade and beyond.
[17]  [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).
[18]  [pdf] - 1853644
The evolution of the cosmic molecular gas density
Comments: 7 pages, 2 figures, Science White paper submitted to Astro2020 Decadal Survey
Submitted: 2019-03-20
One of the last missing pieces in the puzzle of galaxy formation and evolution through cosmic history is a detailed picture of the role of the cold gas supply in the star-formation process. Cold gas is the fuel for star formation, and thus regulates the buildup of stellar mass, both through the amount of material present through a galaxy's gas mass fraction, and through the efficiency at which it is converted to stars. Over the last decade, important progress has been made in understanding the relative importance of these two factors along with the role of feedback, and the first measurements of the volume density of cold gas out to redshift 4, (the "cold gas history of the Universe") has been obtained. To match the precision of measurements of the star formation and black-hole accretion histories over the coming decades, a two orders of magnitude improvement in molecular line survey speeds is required compared to what is possible with current facilities. Possible pathways towards such large gains include significant upgrades to current facilities like ALMA by 2030 (and beyond), and eventually the construction of a new generation of radio-to-millimeter wavelength facilities, such as the next generation Very Large Array (ngVLA) concept.
[19]  [pdf] - 1842572
The Dark Matter Distributions in Low-Mass Disk Galaxies. I. H$\alpha$ Observations Using the Palomar Cosmic Web Imager
Submitted: 2019-02-25
Dark-matter-only simulations predict that dark matter halos have cusp-like inner density profiles, while observations of low-mass galaxies have found a range of inner slopes that are typically much shallower. It is still not well established whether this discrepancy can be explained by baryonic feedback or if it may require modified dark matter models. To better understand the diversity of dark matter profiles in dwarf galaxies, we undertook a survey of 26 low-mass galaxies ($\log M_*/\textrm{M}_\odot = 8.4-9.8$, $v_{\rm max} = 50-140$ km s$^{-1}$) within 30 Mpc using the Palomar Cosmic Web Imager, which is among the largest integral field spectroscopic surveys of its type. In this paper, we derive H$\alpha$ velocity fields for the full sample with a typical spatial resolution of $\sim$160 pc. We extract rotation curves and verify their robustness to several choices in the analysis. We present a method for improving the velocity precision obtained from image slicing spectrographs using narrowband H$\alpha$ images. For 11 galaxies, we compare the H$\alpha$ velocity fields to CO kinematics measured using CARMA, finding the maps to be in good agreement. The standard deviation of the difference is typically $\sim$7 km s$^{-1}$, comparable to the level of turbulence in the interstellar medium, showing that the two tracers have substantially the same bulk kinematics. In a companion paper, we will use the rotation curves produced here to construct mass models of the galaxies and determine their dark matter density profiles.
[20]  [pdf] - 1882526
A diversity of starburst-triggering mechanisms in interacting galaxies and their signatures in CO emission
Comments: submitted to A&A
Submitted: 2019-02-06
The physical origin of enhanced star formation activity in interacting galaxies remains an open question. Knowing whether starbursts are triggered by an increase of the quantity of dense gas or an increase of the star formation efficiency would improve our understanding of galaxy evolution and allow to transpose the results obtained in the local Universe to high redshift galaxies. In this paper, we analyze a parsec-resolution simulation of an Antennae-like model of interacting galaxies. We find that the interplay of physical processes has complex and important variations in time and space, through different combinations of mechanisms like tides, shear and turbulence. These can have similar imprints on observables like depletion time and CO emission. The densest gas only constitutes the tail of the density distribution of some clouds, but exists in large excess in others. The super-linearity of the star formation rate dependence on gas density implies that this excess translates into a reduction of depletion times, and thus leads to a deviation from the classical star formation regime, visible up to galactic scales. These clouds are found in all parts of the galaxies, but their number density varies from one region to the next, due to different cloud assembly mechanisms. Therefore, the dependence of cloud and star formation-related quantities (like CO flux and depletion time) on the scale at which they are measured also varies across the galaxies. We find that the $\alpha_{\rm CO}$ conversion factor between the CO luminosity and molecular gas mass has even stronger spatial than temporal variations in a system like the Antennae. These results raise a number of cautionary notes for the interpretation of observations of unresolved star-forming regions, but also predict that the diversity of environments for star formation will be better captured by the future generations of instruments.
[21]  [pdf] - 1826878
PHIBSS2: survey design and z=0.5-0.8 results. Molecular gas reservoirs during the winding-down of star formation
Comments: Accepted for publication in A&A (28 pages, 9 figures, 4 tables)
Submitted: 2018-12-19, last modified: 2019-01-17
Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.
[22]  [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.
[23]  [pdf] - 1815030
Testing Feedback Regulated Star Formation in Gas Rich, Turbulent Disk Galaxies
Comments: Accepted to ApJ
Submitted: 2018-11-07
In this paper we compare the molecular gas depletion times and mid-plane hydrostatic pressure in turbulent, star forming disk galaxies to internal properties of these galaxies. For this analysis we use 17 galaxies from the DYNAMO sample of nearby ($z\sim0.1$) turbulent disks. We find a strong correlation, such that galaxies with lower molecular gas depletion time (t_{dep}) have higher gas velocity dispersion ($\sigma$). Within the scatter of our data, our observations are consistent with the prediction that depletion time is inversely proportionate to velocity dispersion made in theories of feedback-regulated star formation. We also show a strong, single power-law correlation between mid-plane pressure (P) and star formation rate surface density (Sigma_{SFR}), which extends for 6 orders of magnitude in pressure. Disk galaxies with lower pressure are found to be roughly in agreement with theoretical predictions. However, in galaxies with high pressure we find P/Sigma_{SFR} values that are significantly larger than theoretical predictions. Our observations could be explained with any of the following: (1) the correlation of \Sigma_{SFR}-P is significantly sub-linear; (2) the momentum injected from star formation feedback (so-called p*/m*) is not a single, universal value; or (3) alternate sources of pressure support are important in gas-rich disk galaxies. Finally using published survey results, we find that our results are consistent with the cosmic evolution of t_{dep}(z) and velocity dispersion. Our interpretation of these results is that the cosmic evolution of t_{dep} may be regulated not just by the supply of gas, but also the internal regulation of star formation via feedback.
[24]  [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.
[25]  [pdf] - 1769124
Science with an ngVLA: The ngVLA Science Case and Associated Science Requirements
Comments: 12 pages, 5 figures; To be published in the ASP Monograph Series, "Science with a Next-Generation VLA", ed. E. J. Murphy (ASP, San Francisco, CA). arXiv admin note: text overlap with arXiv:1711.09960
Submitted: 2018-10-15
The science case and associated science requirements for a next-generation Very Large Array (ngVLA) are described, highlighting the five key science goals developed out of a community-driven vision of the highest scientific priorities in the next decade. Building on the superb cm observing conditions and existing infrastructure of the VLA site in the U.S. Southwest, the ngVLA is envisaged to be an interferometric array with more than 10 times the sensitivity and spatial resolution of the current VLA and ALMA, operating at frequencies spanning $\sim1.2 - 116$\,GHz with extended baselines reaching across North America. The ngVLA will be optimized for observations at wavelengths between the exquisite performance of ALMA at submm wavelengths, and the future SKA-1 at decimeter to meter wavelengths, thus lending itself to be highly complementary with these facilities. The ngVLA will be the only facility in the world that can tackle a broad range of outstanding scientific questions in modern astronomy by simultaneously delivering the capability to: (1) unveil the formation of Solar System analogues; (2) probe the initial conditions for planetary systems and life with astrochemistry; (3) characterize the assembly, structure, and evolution of galaxies from the first billion years to the present; (4) use pulsars in the Galactic center as fundamental tests of gravity; and (5) understand the formation and evolution of stellar and supermassive blackholes in the era of multi-messenger astronomy.
[26]  [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.
[27]  [pdf] - 1767834
Science with an ngVLA: How Do Cold Gas Outflows Shape 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
The ngVLA will obtain breakthrough observations of the cold phases (molecular and atomic) of galactic outflows. These observations will be key to study the driving mechanisms of cold winds, the fate of the gas, the mass-loss rates, the conditions that trigger cold outflows, the fraction of gas that escapes galaxies, the effects of winds on suppressing star formation and slowing black hole growth, the statistical properties of outflows, and the relation between large scale outflows and UFOs and BAL outflows. Finally, we highlight the ngVLA complementarity with existing and planned facilities, such as ALMA/SKA, OST/SPICA, and Lynx/ATHENA.
[28]  [pdf] - 1767828
Science with an ngVLA: Gas Density Across the Universe
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
Gas density is widely believed to play a governing role in star formation. However, the exact role of density in setting the star formation rate remains debated. We also lack a general theory that explains how the gas density distribution in galaxies is set. The primary factor preventing the resolution of these issues is the limited number of observations of the gas density distribution across diverse environments. Centimeter- and millimeter-wave spectroscopy offer the most promising way forward in this field, but the key density-sensitive transitions are faint compared to the capabilities of current telescopes. In this chapter, we describe how a next-generation Very Large Array (ngVLA) represents the natural next step forward in this sensitivity-limited field. Such a facility would provide a crucial link between the `Milky Way' and `Extragalactic' views of star formation and dramatically advance our understanding of the drive and role of gas density in galaxies, building on current results from ALMA, NOEMA, the Green Bank Telescope, and other current facilities working in this area.
[29]  [pdf] - 2025327
The Interstellar Medium and Star Formation in Dwarf Galaxies
Comments: Invited review at IAU Symposium 344 "Dwarf Galaxies: From the Deep Universe to the Present"
Submitted: 2018-10-12
This is a brief review of our understanding of the properties of the interstellar medium (ISM) in dwarf galaxies in connection to their star formation activity. What are the dominant phases of the ISM in these objects? How do the properties of these phases depend on the galaxy properties? What do we know about their cold gas content and its link to star formation activity? Does star formation proceed differently in these galaxies? How does star formation feedback operate in dwarf galaxies? The availability of observations from space-based facilities such as FUSE, Spitzer, Herschel, and Fermi, as well as observatories such as SOFIA and ALMA, is allowing us to make significant strides in our understanding of these questions.
[30]  [pdf] - 1791151
ALMA observations of HCN and HCO+ outflows in the merging galaxy NGC 3256
Comments: 22 pages, 12 figures, 5 tables, accepted for publication in APJ
Submitted: 2018-10-10, last modified: 2018-10-12
We report ~2" resolution Atacama Large Millimeter/submillimeter Array observations of the HCN(1-0), HCO+(1-0), CO(1-0), CO(2-1), and CO(3-2) lines towards the nearby merging double-nucleus galaxy NGC 3256. We find that the high density gas outflow traced in HCN(1-0) and HCO+(1-0) emission is co-located with the diffuse molecular outflow emanating from the southern nucleus, where a low-luminosity active galactic nucleus (AGN) is believed to be the dominant source of the far-infrared luminosity. On the other hand, the same lines were undetected in the outflow region associated with the northern nucleus, whose primary heating source is likely related to starburst activity without obvious signs of AGN. Both HCO+(1-0)/CO(1-0) line ratio (i.e. dense gas fraction) and the CO(3-2)/CO(1-0) line ratio are larger in the southern outflow (0.20$\pm$0.04 and 1.3$\pm$0.2, respectively) than in the southern nucleus (0.08$\pm$0.01, 0.7$\pm$0.1, respectively). By investigating these line ratios for each velocity component in the southern outflow, we find that the dense gas fraction increases and the CO(3-2)/CO(1-0) line ratio decreases towards the largest velocity offset. This suggests the existence of a two-phase (diffuse and clumpy) outflow. One possible scenario to produce such a two-phase outflow is an interaction between the jet and the interstellar medium, which possibly triggers shocks and/or star formation associated with the outflow.
[31]  [pdf] - 1783882
Spatially Resoved 12CO(2-1)/12CO(1-0) in the Starburst Galaxy NGC 253: Assessing Optical Depth to Constrain the Molecular Mass Outflow Rate
Comments: Accepted to ApJ on 31 August 2018, 14 pages, 8 figures
Submitted: 2018-09-04
We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO(1-0) and 12CO(2-1) in the central 40" (680 pc) of the nuclear starburst galaxy NGC 253, including its molecular outflow. We measure the ratio of brightness temperature for CO(2-1)/CO(1-0), r_21, in the central starburst and outflow-related features. We discuss how r_21 can be used to constrain the optical depth of the CO emission, which impacts the inferred mass of the outflow and consequently the molecular mass outflow rate. We find r_21 less than or equal to 1 throughout, consistent with a majority of the CO emission being optically-thick in the outflow, as it is in the starburst. This suggests that the molecular outflow mass is 3-6 times larger than the lower limit reported for optically thin CO emission from warm molecular gas. The implied molecular mass outflow rate is 25-50 solar masses per year, assuming that conversion factor for the outflowing gas is similar to our best estimates for the bulk of the starburst. This is a factor of 9-19 times larger than the star formation rate in NGC 253. We see tentative evidence for an extended, diffuse CO(2-1) component.
[32]  [pdf] - 1822832
The Molecular and Ionized Gas Phases of an AGN-driven Outflow in a Typical Massive Galaxy at z=2
Comments: 20 pages, 14 figures. Submitted to ApJ
Submitted: 2018-07-18, last modified: 2018-07-21
Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming, disk galaxies and quiescent spheroidal systems. Yet, our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work we use new deep, high-spatial resolution ALMA CO(3-2) and archival VLT/SINFONI H$\alpha$ observations to study the molecular and ionized components of the AGN-driven outflow in zC400528 ---a massive, main sequence galaxy at z=2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient and extends for at least ~10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale at least twice shorter than that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star-formation activity in zC400528.
[33]  [pdf] - 1779582
Evidence for Cosmic-Ray Escape in the Small Magellanic Cloud using Fermi Gamma-rays
Comments: 12 pages, 8 figures, submitted to ApJ
Submitted: 2018-07-17
Galaxy formation simulations demonstrate that cosmic-ray (CR) feedback may be important in the launching of galactic-scale winds. CR protons dominate the bulk of the CR population, yet most observational constraints of CR feedback come from synchrotron emission of CR electrons. In this paper, we present an analysis of 105 months of Fermi Gamma-ray Space Telescope observations of the Small Magellanic Cloud (SMC), with the aim of exploring CR feedback and transport in an external galaxy. We produce maps of the 2-300 GeV emission and detect statistically significant, extended emission along the Bar and the Wing, where active star formation is occurring. Gamma-ray emission is not detected above 13 GeV, and we set stringent upper-limits on the flux above this energy. We find the best fit to the gamma-ray spectrum is a single-component model with a power-law of index $\Gamma=-2.11\pm0.06\pm0.06$ and an exponential cutoff energy of $E_{\rm c} =13.1\pm5.1\pm1.6$ GeV. We assess the relative contribution of pulsars and CRs to the emission, and we find that pulsars may produce up to 14$^{+4}_{-2}$% of the flux above 100 MeV. Thus, we attribute most of the gamma-ray emission (based on its spectrum and morphology) to CR interactions with the ISM. We show that the gamma-ray emissivity of the SMC is five times smaller than that of the Milky Way and that the SMC is far below the calorimetric limit, where all CR protons experience pion losses. We interpret these findings as evidence that CRs are escaping the SMC via advection and diffusion.
[34]  [pdf] - 1705260
The EDGE-CALIFA Survey: Molecular and Ionized Gas Kinematics in Nearby Galaxies
Comments: 37 pages, 18 figures, accepted for publication in ApJ
Submitted: 2018-04-16, last modified: 2018-05-03
We present a comparative study of molecular and ionized gas kinematics in nearby galaxies. These results are based on observations from the EDGE survey, which measured spatially resolved $^{12}$CO(J=1-0) in 126 nearby galaxies. Every galaxy in EDGE has corresponding resolved ionized gas measurements from CALIFA. Using a sub-sample of 17 rotation dominated, star-forming galaxies where precise molecular gas rotation curves could be extracted, we derive CO and H$\alpha$ rotation curves using the same geometric parameters out to $\gtrsim$1 $R_e$. We find that $\sim$75% of our sample galaxies have smaller ionized gas rotation velocities than the molecular gas in the outer part of the rotation curve. In no case is the molecular gas rotation velocity measurably lower than that of the ionized gas. We suggest that the lower ionized gas rotation velocity can be attributed to a significant contribution from extraplanar diffuse ionized gas in a thick, turbulence supported disk. Using observations of the H$\gamma$ transition also available from CALIFA, we measure ionized gas velocity dispersions and find that these galaxies have sufficiently large velocity dispersions to support a thick ionized gas disk. Kinematic simulations show that a thick disk with a vertical rotation velocity gradient can reproduce the observed differences between the CO and H$\alpha$ rotation velocities. Observed line ratios tracing diffuse ionized gas are elevated compared to typical values in the midplane of the Milky Way. In galaxies affected by this phenomenon, dynamical masses measured using ionized gas rotation curves will be systematically underestimated.
[35]  [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.
[36]  [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.
[37]  [pdf] - 1648953
The EDGE-CALIFA survey: validating stellar dynamical mass models with CO kinematics
Comments: 22 pages, 18 figures, Accepted for publication in MNRAS
Submitted: 2018-03-06
Deriving circular velocities of galaxies from stellar kinematics can provide an estimate of their total dynamical mass, provided a contribution from the velocity dispersion of the stars is taken into account. Molecular gas (e.g., CO) on the other hand, is a dynamically cold tracer and hence acts as an independent circular velocity estimate without needing such a correction. In this paper we test the underlying assumptions of three commonly used dynamical models, deriving circular velocities from stellar kinematics of 54 galaxies (S0-Sd) that have observations of both stellar kinematics from the CALIFA survey, and CO kinematics from the EDGE survey. We test the Asymmetric Drift Correction (ADC) method, as well as Jeans, and Schwarzschild models. The three methods each reproduce the CO circular velocity at 1Re to within 10%. All three methods show larger scatter (up to 20%) in the inner regions (R < 0.4Re) which may be due to an increasingly spherical mass distribution (which is not captured by the thin disk assumption in ADC), or non-constant stellar M/L ratios (for both the JAM and Schwarzschild models). This homogeneous analysis of stellar and gaseous kinematics validates that all three models can recover Mdyn at 1Re to better than 20%, but users should be mindful of scatter in the inner regions where some assumptions may break down.
[38]  [pdf] - 1663318
HI Kinematics Along The Minor Axis of M82
Comments: ApJ accepted, 12 pages, 8 figures
Submitted: 2018-02-12, last modified: 2018-02-13
M82 is one of the best studied starburst galaxies in the local universe, and is consequently a benchmark for studying star formation feedback at both low and high redshift. We present new VLA HI observations that reveal the cold gas kinematics along the minor axis in unprecedented detail. This includes the detection of HI up to 10 kpc along the minor axis toward the South and beyond 5 kpc to the North. A surprising aspect of these observations is that the line-of-sight HI velocity decreases substantially from about 120 km/s to 50 km/s from 1.5 to 10 kpc off the midplane. The velocity profile is not consistent with the HI gas cooling from the hot wind. We demonstrate that the velocity decrease is substantially greater than the deceleration expected from gravitational forces alone. If the HI consists of a continuous population of cold clouds, some additional drag force must be present, and the magnitude of the drag force places a joint constraint on the ratio of the ambient medium to the typical cloud size and density. We also show that the HI kinematics are inconsistent with a simple conical outflow centered on the nucleus, but instead require the more widespread launch of the HI over the ~1 kpc extent of the starburst region. Regardless of the launch mechanism for the HI gas, the observed velocity decrease along the minor axis is sufficiently great that the HI may not escape the halo of M82. We estimate the HI outflow rate is much less than 1 M$_{\odot}$ per year at 10 kpc off the midplane.
[39]  [pdf] - 1637458
PHIBSS: Unified Scaling Relations of Gas Depletion Time and Molecular Gas Fractions
Comments: Now accepted by ApJ; Figures and tables have been updated, and an appendix has been added
Submitted: 2017-02-03, last modified: 2018-01-11
This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.
[40]  [pdf] - 1634314
First Results from the $Herschel$ and ALMA Spectroscopic Surveys of the SMC: The Relationship Between [CII]-bright Gas and CO-bright Gas at Low Metallicity
Comments: Accepted for publication in ApJ
Submitted: 2018-01-10
The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverage overlapping four of the five HS$^{3}$ regions. We detect [CII] and [OI] throughout all of the regions mapped. The data allow us to compare the structure of the molecular clouds and surrounding photodissociation regions using $^{13}$CO, CO, [CII], and [OI] emission at $<10$" ($<3$ pc) scales. We estimate Av using far-IR thermal continuum emission from dust and find the CO/[CII] ratios reach the Milky Way value at high A$_{V}$ in the centers of the clouds and fall to $\sim{1/5-1/10}\times$ the Milky Way value in the outskirts, indicating the presence of translucent molecular gas not traced by bright CO emission. We estimate the amount of molecular gas traced by bright [CII] emission at low A$_{V}$ and bright CO emission at high A$_{V}$. We find that most of the molecular gas is at low A$_{V}$ and traced by bright [CII] emission, but that faint CO emission appears to extend to where we estimate the H$_{2}$-to-HI transition occurs. By converting our H$_{2}$ gas estimates to a CO-to-H$_{2}$ conversion factor ($X_{CO}$), we show that $X_{CO}$ is primarily a function of A$_{V}$, consistent with simulations and models of low metallicity molecular clouds.
[41]  [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.
[42]  [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.
[43]  [pdf] - 1626410
Probing the baryon cycle of galaxies with SPICA mid- and far-infrared observations
Comments: Accepted for publication in PASA; 19 pages, 8 figures, 1 table
Submitted: 2017-11-30
The SPICA mid and far-infrared telescope will address fundamental issues in our understanding of star formation and ISM physics in galaxies. A particular hallmark of SPICA is the outstanding sensitivity enabled by the cold telescope, optimized detectors, and wide instantaneous bandwidth throughout the mid- and far-infrared. The spectroscopic, imaging and polarimetric observations that SPICA will be able to collect will help in clarifying the complex physical mechanisms which underlie the baryon cycle of galaxies. In particular: (i) The access to a large suite of atomic and ionic fine-structure lines for large samples of galaxies will shed light on the origin of the observed spread in star formation rates within and between galaxies. (ii) Observations of HD rotational lines (out to $\sim$10 Mpc) and fine structure lines such as [CII] 158 $\mu$m (out to $\sim$100 Mpc) will clarify the main reservoirs of interstellar matter in galaxies, including phases where CO does not emit. (iii) Far-infrared spectroscopy of dust and ice features will address uncertainties in the mass and composition of dust in galaxies, and the contributions of supernovae to the interstellar dust budget will be quantified by photometry and monitoring of supernova remnants in nearby galaxies. (iv) Observations of far-infrared cooling lines such as [OI] 63 $\mu$m from star-forming molecular clouds in our Galaxy will evaluate the importance of shocks to dissipate turbulent energy. The paper concludes with requirements for the telescope and instruments, and recommendations for the observing strategy.
[44]  [pdf] - 1593816
Key Science Goals for the Next Generation Very Large Array (ngVLA): Report from the ngVLA Science Advisory Council
Comments: ngVLA memo 19
Submitted: 2017-11-27
This document describes some of the fundamental astrophysical problems that require observing capabilities at millimeter- and centimeter wavelengths well beyond those of existing, or already planned, telescopes. The results summarized in this report follow a solicitation from the National Radio Astronomy Observatory to develop key science cases for a future U. S.-led radio telescope, the "next generation Very Large Array" (ngVLA). The ngVLA will have roughly 10 times the collecting area of the Jansky VLA, operate at frequencies from 1 GHz to 116 GHz with up to 20 GHz of bandwidth, possess a compact core for high surface-brightness sensitivity, and extended baselines of at least hundreds of kilometers and ultimately across the continent to provide high-resolution imaging. The ngVLA builds on the scientific and technical legacy of the Jansky VLA and ALMA, and will be designed to provide the next leap forward in our understanding of planets, galaxies, and black holes.
[45]  [pdf] - 1614859
The Parsec-Scale Relationship Between I_CO and A_V in Local Molecular Clouds
Comments: Accepted for publication in MNRAS, 16 pages, 9 figures in main text, manuscript also includes a 24 page atlas showing each cloud. Full resolution available here:
Submitted: 2017-10-23
We measure the parsec-scale relationship between integrated CO intensity (I_CO) and visual extinction (A_V) in 24 local molecular clouds using maps of CO emission and dust optical depth from Planck. This relationship informs our understanding of CO emission across environments, but clean Milky Way measurements remain scarce. We find uniform I_CO for a given A_V, with the results bracketed by previous studies of the Pipe and Perseus clouds. Our measured I_CO-A_V relation broadly agrees with the standard Galactic CO-to-H2 conversion factor, the relation found for the Magellanic clouds at coarser resolution, and numerical simulations by Glover & Clark (2016). This supports the idea that CO emission primarily depends on shielding, which protects molecules from dissociating radiation. Evidence for CO saturation at high A_V and a threshold for CO emission at low A_V varies remains uncertain due to insufficient resolution and ambiguities in background subtraction. Resolution of order 0.1 pc may be required to measure these features. We use this I_CO-AV relation to predict how the CO-to-H2 conversion factor (X_CO) would change if the Solar Neighborhood clouds had different dust-to-gas ratio (metallicity). The calculations highlight the need for improved observations of the CO emission threshold and HI shielding layer depth. They are also sensitive to the shape of the column density distribution. Because local clouds collectively show a self-similar distribution, we predict a shallow metallicity dependence for X_CO down to a few tenths of solar metallicity. However, our calculations also imply dramatic variations in cloud-to-cloud X_CO at subsolar metallicity.
[46]  [pdf] - 1593682
Galaxy evolution studies with the SPace IR telescope for Cosmology and Astrophysics (SPICA): the power of IR spectroscopy
Comments: Paper accepted for publication on PASA on 21st September 2017, as part of the SPICA Special Issue
Submitted: 2017-10-05
IR spectroscopy in the range 12-230 micron with the SPace IR telescope for Cosmology and Astrophysics (SPICA) will reveal the physical processes that govern the formation and evolution of galaxies and black holes through cosmic time, bridging the gap between the James Webb Space Telescope (JWST) and the new generation of Extremely Large Telescopes (ELTs) at shorter wavelengths and the Atacama Large Millimeter Array (ALMA) at longer wavelengths. SPICA, with its 2.5-m telescope actively-cooled to below 8K, will obtain the first spectroscopic determination, in the mid-IR rest-frame, of both the star-formation rate and black hole accretion rate histories of galaxies, reaching lookback times of 12 Gyr, for large statistically significant samples. Densities, temperatures, radiation fields and gas-phase metallicities will be measured in dust-obscured galaxies and active galactic nuclei (AGN), sampling a large range in mass and luminosity, from faint local dwarf galaxies to luminous quasars in the distant Universe. AGN and starburst feedback and feeding mechanisms in distant galaxies will be uncovered through detailed measurements of molecular and atomic line profiles. SPICA's large-area deep spectrophotometric surveys will provide mid-IR spectra and continuum fluxes for unbiased samples of tens of thousands of galaxies, out to redshifts of z~6. Furthermore, SPICA spectroscopy will uncover the most luminous galaxies in the first few hundred million years of the Universe, through their characteristic dust and molecular hydrogen features.
[47]  [pdf] - 1589349
SPICA and the Chemical Evolution of Galaxies: The Rise of Metals and Dust
Comments: In press. This paper belongs to the "SPICA Special Issue" on PASA
Submitted: 2017-10-05
The physical processes driving the chemical evolution of galaxies in the last $\sim 11\, \rm{Gyr}$ cannot be understood without directly probing the dust-obscured phase of star-forming galaxies and active galactic nuclei. This phase, hidden to optical tracers, represents the bulk of star formation and black hole accretion activity in galaxies at $1 < z < 3$. Spectroscopic observations with a cryogenic infrared (IR) observatory like SPICA will be sensitive enough to peer through the dust-obscured regions of galaxies and access the rest-frame mid- to far-IR range in galaxies at high-$z$. This wavelength range contains a unique suite of spectral lines and dust features that serve as proxies for the abundances of heavy elements and the dust composition, providing tracers with a feeble response to both extinction and temperature. In this work, we investigate how SPICA observations could be exploited to understand key aspects in the chemical evolution of galaxies: the assembly of nearby galaxies based on the spatial distribution of heavy element abundances, the global content of metals in galaxies reaching the knee of the luminosity function up to $z \sim 3$, and the dust composition of galaxies at high-$z$. Possible synergies with facilities available in the late 2020s are also discussed.
[48]  [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}$.
[49]  [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.
[50]  [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.
[51]  [pdf] - 1586237
Gas Content and Kinematics in Clumpy, Turbulent Star-forming Disks
Comments: 10 pages, 6 figures. Accepted to ApJ
Submitted: 2017-07-21
We present molecular gas mass estimates for a sample of 13 local galaxies whose kinematic and star forming properties closely resemble those observed in $z\approx 1.5$ main-sequence galaxies. Plateau de Bure observations of the CO[1-0] emission line and Herschel Space Observatory observations of the dust emission both suggest molecular gas mass fractions of ~20%. Moreover, dust emission modeling finds $T_{dust}<$30K, suggesting a cold dust distribution compared to their high infrared luminosity. The gas mass estimates argue that $z\sim$0.1 DYNAMO galaxies not only share similar kinematic properties with high-z disks, but they are also similarly rich in molecular material. Pairing the gas mass fractions with existing kinematics reveals a linear relationship between $f_{gas}$ and $\sigma$/$v_{c}$, consistent with predictions from stability theory of a self-gravitating disk. It thus follows that high gas velocity dispersions are a natural consequence of large gas fractions. We also find that the systems with lowest depletion times ($\sim$0.5 Gyr) have the highest ratios of $\sigma$/$v_{c}$ and more pronounced clumps, even at the same high molecular gas fraction.
[52]  [pdf] - 1584737
ALMA observations of N83C in the early stage of star formation in the Small Magellanic Cloud
Comments: 12 pages, 4 figures, 2 tables, accepted for publication in ApJ
Submitted: 2017-06-15, last modified: 2017-07-14
We have performed Atacama Large Millimeter/submillimeter Array (ALMA) observations in $^{12}$CO($J=2-1$), $^{13}$CO($J=2-1$), C$^{18}$O($J=2-1$), $^{12}$CO($J=3-2$), $^{13}$CO($J=3-2$), and CS($J=7-6$) lines toward the active star-forming region N83C in the Small Magellanic Cloud (SMC), whose metallicity is $\sim$ 1/5 of the Milky Way (MW). The ALMA observations first reveal sub-pc scale molecular structures in $^{12}$CO($J=2-1$) and $^{13}$CO($J=2-1$) emission. We found strong CO peaks associated with young stellar objects (YSOs) identified by the $Spitzer$ Space Telescope, and also found that overall molecular gas is distributed along the edge of the neighboring HII region. We derived a gas density of $\sim 10^4$ cm$^{-3}$ in molecular clouds associated with YSOs based on the virial mass estimated from $^{12}$CO($J=2-1$) emission. This high gas density is presumably due to the effect of the HII region under the low-metallicity (accordingly small-dust content) environment in the SMC; far-UV radiation from the HII region can easily penetrate and photo-dissociate the outer layer of $^{12}$CO molecules in the molecular clouds, and thus only the innermost parts of the molecular clouds are observed even in $^{12}$CO emission. We obtained the CO-to-H$_2$ conversion factor $X_{\rm CO}$ of $7.5 \times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$ in N83C based on virial masses and CO luminosities, which is four times larger than that in the MW, 2 $\times 10^{20}$ cm$^{-2}$ (K km s$^{-1}$)$^{-1}$. We also discuss the difference in the nature between two high-mass YSOs, each of which is associated with a molecular clump with a mass of about a few $\times 10^3 M_{\odot}$.
[53]  [pdf] - 1585938
The Origins of [CII] Emission in Local Star-forming Galaxies
Comments: 8 pages, accepted for publication in ApJ
Submitted: 2017-07-14
The [CII] 158um fine-structure line is the brightest emission line observed in local star-forming galaxies. As a major coolant of the gas-phase interstellar medium, [CII] balances the heating, including that due to far-ultraviolet photons, which heat the gas via the photoelectric effect. However, the origin of [CII] emission remains unclear, because C+ can be found in multiple phases of the interstellar medium. Here we measure the fractions of [CII] emission originating in the ionized and neutral gas phases of a sample of nearby galaxies. We use the [NII] 205um fine-structure line to trace the ionized medium, thereby eliminating the strong density dependence that exists in the ratio of [CII]/[NII] 122um. Using the FIR [CII] and [NII] emission detected by the KINGFISH and Beyond the Peak Herschel programs, we show that 60-80% of [CII] emission originates from neutral gas. We find that the fraction of [CII] originating in the neutral medium has a weak dependence on dust temperature and the surface density of star formation, and a stronger dependence on the gas-phase metallicity. In metal-rich environments, the relatively cooler ionized gas makes substantially larger contributions to total [CII] emission than at low abundance, contrary to prior expectations. Approximate calibrations of this metallicity trend are provided.
[54]  [pdf] - 1585765
Spatial variations of turbulent properties of neutral hydrogen gas in the Small Magellanic Cloud using structure function analysis
Comments: 10 pages, 7 figures, accepted for publication in ApJ
Submitted: 2017-07-10
We investigate spatial variations of turbulent properties in the Small Magellanic Cloud (SMC) by using neutral hydrogen HI observations. With the goal of testing the importance of stellar feedback on HI turbulence, we define central and outer SMC regions based on the star formation rate (SFR) surface density, as well as the HI integrated intensity. We use the structure function and the Velocity Channel Analysis (VCA) to calculate the power-law index (gamma) for both underlying density and velocity fields in these regions. In all cases, our results show essentially no difference in gamma between the central and outer regions. This suggests that HI turbulent properties are surprisingly homogeneous across the SMC when probed at a resolution of 30 pc. Contrary to recent suggestions from numerical simulations, we do not find a significant change in gamma due to stellar feedback as traced by the SFR surface density. This could be due to the stellar feedback being widespread over the whole of the SMC, but more likely due to a large-scale gravitational driving of turbulence. We show that the lack of difference between central and outer SMC regions can not be explained by the high optical depth HI.
[55]  [pdf] - 1585161
Cloud Scale ISM Structure and Star Formation in M51
Comments: Accepted for publication in the Astrophysical Journal, 22 pages, 11 figures, email for access to data table before publication, for high resolution version before publication:
Submitted: 2017-06-26
We compare the structure of molecular gas at $40$ pc resolution to the ability of gas to form stars across the disk of the spiral galaxy M51. We break the PAWS survey into $370$ pc and $1.1$ kpc resolution elements, and within each we estimate the molecular gas depletion time ($\tau_{\rm Dep}^{\rm mol}$), the star formation efficiency per free fall time ($\epsilon_{\rm ff}$), and the mass-weighted cloud-scale (40 pc) properties of the molecular gas: surface density, $\Sigma$, line width, $\sigma$, and $b\equiv\Sigma/\sigma^2\propto\alpha_{\rm vir}^{-1}$, a parameter that traces the boundedness of the gas. We show that the cloud-scale surface density appears to be a reasonable proxy for mean volume density. Applying this, we find a typical star formation efficiency per free-fall time, $\epsilon_{ff} \left( \left< \Sigma_{40pc} \right> \right) \sim 0.3{-}0.36\%$, lower than adopted in many models and found for local clouds. More, the efficiency per free fall time anti-correlates with both $\Sigma$ and $\sigma$, in some tension with turbulent star formation models. The best predictor of the rate of star formation per unit gas mass in our analysis is $b \equiv \Sigma / \sigma^2$, tracing the strength of self gravity, with $\tau_{\rm Dep}^{\rm mol} \propto b^{-0.9}$. The sense of the correlation is that gas with stronger self-gravity (higher $b$) forms stars at a higher rate (low $\tau_{\rm Dep}^{\rm mol}$). The different regions of the galaxy mostly overlap in $\tau_{\rm Dep}^{\rm mol}$ as a function of $b$, so that low $b$ explains the surprisingly high $\tau_{\rm Dep}^{\rm mol}$ found towards the inner spiral arms found by by Meidt et al. (2013).
[56]  [pdf] - 1584160
Quasar Feedback in the Ultraluminous Infrared Galaxy F11119+3257: Connecting the Accretion Disk Wind with the Large-Scale Molecular Outflow
Comments: 15 pages, 6 figures, 4 tables, accepted for publication in ApJ
Submitted: 2017-06-01
In Tombesi et al. (2015), we reported the first direct evidence for a quasar accretion disk wind driving a massive molecular outflow. The target was F11119+3257, an ultraluminous infrared galaxy (ULIRG) with unambiguous type-1 quasar optical broad emission lines. The energetics of the accretion disk wind and molecular outflow were found to be consistent with the predictions of quasar feedback models where the molecular outflow is driven by a hot energy-conserving bubble inflated by the inner quasar accretion disk wind. However, this conclusion was uncertain because the energetics were estimated from the optically thick OH 119 um transition profile observed with Herschel. Here, we independently confirm the presence of the molecular outflow in F11119+3257, based on the detection of broad wings in the CO(1-0) profile derived from ALMA observations. The broad CO(1-0) line emission appears to be spatially extended on a scale of at least ~7 kpc from the center. Mass outflow rate, momentum flux, and mechanical power of (80-200) R_7^{-1} M_sun/yr, (1.5-3.0) R_7^{-1} L_AGN/c, and (0.15-0.40)% R_7^{-1} L_AGN are inferred from these data, assuming a CO-to-H_2 conversion factor appropriate for a ULIRG (R_7 is the radius of the outflow normalized to 7 kpc and L_AGN is the AGN luminosity). These rates are time-averaged over a flow time scale of 7x10^6 yrs. They are similar to the OH-based rates time-averaged over a flow time scale of 4x10^5 yrs, but about a factor 4 smaller than the local ("instantaneous"; <10^5 yrs) OH-based estimates cited in Tombesi et al. The implications of these new results are discussed in the context of time-variable quasar-mode feedback and galaxy evolution. The need for an energy-conserving bubble to explain the molecular outflow is also re-examined.
[57]  [pdf] - 1561195
Connecting Clump Sizes in Turbulent Disk Galaxies to Instability Theory
Comments: Accepted to Apj Letters
Submitted: 2017-03-01
In this letter we study the mean sizes of Halpha clumps in turbulent disk galaxies relative to kinematics, gas fractions, and Toomre Q. We use 100~pc resolution HST images, IFU kinematics, and gas fractions of a sample of rare, nearby turbulent disks with properties closely matched to z~1.5-2 main-sequence galaxies (the DYNAMO sample). We find linear correlations of normalized mean clump sizes with both the gas fraction and the velocity dispersion-to-rotation velocity ratio of the host galaxy. We show that these correlations are consistent with predictions derived from a model of instabilities in a self-gravitating disk (the so-called "violent disk instability model"). We also observe, using a two-fluid model for Q, a correlation between the size of clumps and self-gravity driven unstable regions. These results are most consistent with the hypothesis that massive star forming clumps in turbulent disks are the result of instabilities in self-gravitating gas-rich disks, and therefore provide a direct connection between resolved clump sizes and this in situ mechanism.
[58]  [pdf] - 1546390
Updated 34-Band Photometry for the SINGS/KINGFISH Samples of Nearby Galaxies
Comments: Accepted for publication in the Astrophysical Journal; object, sky, and contaminant apertures are available from the first author
Submitted: 2017-02-16
We present an update to the ultraviolet-to-radio database of global broadband photometry for the 79 nearby galaxies that comprise the union of the KINGFISH (Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel) and SINGS (Spitzer Infrared Nearby Galaxies Survey) samples. The 34-band dataset presented here includes contributions from observational work carried out with a variety of facilities including GALEX, SDSS, PS, NOAO, 2MASS, WISE, Spitzer, Herschel, Planck, JCMT, and the VLA. Improvements of note include recalibrations of previously-published SINGS BVRcIc and KINGFISH far-infrared/submillimeter photometry. Similar to previous results in the literature, an excess of submillimeter emission above model predictions is seen primarily for low-metallicity dwarf/irregular galaxies. This 34-band photometric dataset for the combined KINGFISH$+$SINGS sample serves as an important multi-wavelength reference for the variety of galaxies observed at low redshift. A thorough analysis of the observed spectral energy distributions is carried out in a companion paper.
[59]  [pdf] - 1636597
Molecular gas at supernova local environments unveiled by EDGE
Comments: MNRAS accepted, 17 pages, 8 Figures, 4 Tables
Submitted: 2017-02-09
CO observations allow estimations of the gas content of molecular clouds, which trace the reservoir of cold gas fuelling star formation, as well as to determine extinction via H$_2$ column density, N(H$_2$). Here, we studied millimetric and optical properties at 26 supernovae (SNe) locations of different types in a sample of 23 nearby galaxies by combining molecular $^{12}$C$^{16}$O (J = 1 $\rightarrow$ 0) resolved maps from the EDGE survey and optical Integral Field Spectroscopy from the CALIFA survey. We found an even clearer separation between type II and type Ibc SNe in terms of molecular gas than what we found in the optical using H$\alpha$ emission as a proxy for current SF rate, which reinforces the fact that SNe Ibc are more associated with SF-environments. While A$_V$ at SN locations is similar for SNe II and SNe Ibc, and higher compared to SNe Ia, N(H$_2$) is significantly higher for SNe Ibc than for SNe II and SNe Ia. When compared to alternative extinction estimations directly from SN photometry and spectroscopy, we find that our SNe Ibc have also redder color excess but showed standard Na I D absorption pseudo-equivalent widths ($\sim$1 \AA). In some cases we find no extinction when estimated from the environment, but high amounts of extinction when measured from SN observations, which suggests that circumstellar material or dust sublimation may be playing a role. This work serves as a benchmark for future studies combining last generation millimeter and optical IFS instruments to reveal the local environmental properties of extragalactic SNe.
[60]  [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.
[61]  [pdf] - 1534388
Physical Properties of Molecular Clouds at 2 parsec Resolution in the Low-Metallicity Dwarf Galaxy NGC 6822 and the Milky Way
Comments: Accepted for publication in The Astrophysical Journal; 22 pages, 10 figures, 7 tables
Submitted: 2017-01-10
We present the ALMA survey of CO(2-1) emission from the 1/5 solar metallicity, Local Group dwarf galaxy NGC 6822. We achieve high (0.9 arcsec ~ 2 pc) spatial resolution while covering large area: four 250 pc x 250 pc regions that encompass ~2/3 of NGC 6822's star formation. In these regions, we resolve ~150 compact CO clumps that have small radii (~2-3 pc), narrow line width (~1 km/s), and low filling factor across the galaxy. This is consistent with other recent studies of low metallicity galaxies, but here shown with a 15 times larger sample. At parsec scales, CO emission correlates with 8 micron emission better than with 24 micron emission and anti-correlates with Halpha, so that PAH emission may be an effective tracer of molecular gas at low metallicity. The properties of the CO clumps resemble those of similar-size structures in Galactic clouds except of slightly lower surface brightness and CO-to-H2 ratio ~1-2 times the Galactic value. The clumps exist inside larger atomic-molecular complexes with masses typical for giant molecular cloud. Using dust to trace H2 for the entire complex, we find CO-to-H2 to be ~20-25 times the Galactic value, but with strong dependence on spatial scale and variations between complexes that may track their evolutionary state. The H2-to-HI ratio is low globally and only mildly above unity within the complexes. The SFR-to-H2 ratio is ~3-5 times higher in the complexes than in massive disk galaxies, but after accounting for the bias from targeting star-forming regions, we conclude that the global molecular gas depletion time may be as long as in massive disk galaxies.
[62]  [pdf] - 1534188
Thermal Pressure in the Cold Neutral Medium of Nearby Galaxies
Comments: Accepted for publication in ApJ on Dec. 20th, 2016; 17 pages; 7 figures and 1 table
Submitted: 2017-01-05
Dynamic and thermal processes regulate the structure of the multi-phase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure ($P_{\rm th}$) of the diffuse, neutral gas. By combining [C II] 158 $\mu$m, HI, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of $\sim$1 kiloparsec. We find a distribution of thermal pressures in the $P_{\rm th}/k\sim10^3-10^5$ K cm$^{-3}$ range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of $P_{\rm th}/k\approx3600$ K cm$^{-3}$. These results are consistent with thermal pressure measurements using other observational methods. We find that $P_{\rm th}$ increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral medium could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about $\sim30$% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.
[63]  [pdf] - 1542784
The Radio Spectral Energy Distribution and Star Formation Rate Calibration in Galaxies
Comments: 29 pages, 15 figures, accepted for publication in Astrophysical Journal
Submitted: 2016-11-05, last modified: 2016-12-28
We study the spectral energy distribution (SED) of the radio continuum emission from the KINGFISH sample of nearby galaxies to understand the energetics and origin of this emission. Effelsberg multi-wavelength observations at 1.4GHz, 4.8GHz, 8.5GHz, and 10.5GHz combined with archive data allow us, for the first time, to determine the mid-radio continuum (1-10 GHz, MRC) bolometric luminosities and further present calibration relations vs. the monochromatic radio luminosities. The 1-10 GHz radio SED is fitted using a Bayesian Markov Chain Monte Carlo (MCMC) technique leading to measurements for the nonthermal spectral index and the thermal fraction f_th with mean values of alpha_nt=0.97+-0.16 (0.79+-0.15 for the total spectral index) and f_th= 10% +- 9% at 1.4 GHz. The MRC luminosity changes over ~3 orders of magnitude in the sample. The thermal emission is responsible for ~23% of the MRC on average. We also compare the extinction-corrected diagnostics of star formation rate with the thermal and nonthermal radio tracers and derive the first star formation calibration relations using the MRC radio luminosity. The nonthermal spectral index flattens with increasing star formation rate surface density, indicating the effect of the star formation feedback on the cosmic ray electron population in galaxies. Comparing the radio and IR SEDs, we find that the FIR-to-MRC ratio could decrease with star formation rate, due to the amplification of the magnetic fields in star forming regions. This particularly implies a decrease in the ratio at high redshifts, where mostly luminous/star forming galaxies are detected.
[64]  [pdf] - 1561046
PHIBSS: Exploring the Dependence of the CO-H$_2$ Conversion Factor on Total Mass Surface Density at ${\it z} < 1.5$
Comments: Accepted to MNRAS. 16 pages, 8 figures
Submitted: 2016-11-14, last modified: 2016-12-15
We present an analysis of the relationship between the CO-H$_{2}$ conversion factor ($\alpha_{\rm CO}$) and total mass surface density ($\Sigma_{\rm tot}$) in star-forming galaxies at $z < 1.5$. Our sample, which is drawn from the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) and the CO Legacy Database for GASS (COLD GASS), includes 'normal,' massive star-forming galaxies that dominate the evolution of the cosmic star formation rate (SFR) at this epoch and probe the $\Sigma_{\rm tot}$ regime where the strongest variation in $\alpha_{\rm CO}$ is observed. We constrain $\alpha_{\rm CO}$ via existing CO observations, measurements of the star formation rate, and an assumed molecular gas depletion time ($t_{\rm dep}$=$M_{\rm gas}$/SFR) --- the latter two of which establish the total molecular gas mass independent of the observed CO luminosity. For a broad range of adopted depletion times, we find that $\alpha_{\rm CO}$ is independent of total mass surface density, with little deviation from the canonical Milky Way value. This runs contrary to a scenario in which $\alpha_{\rm CO}$ decreases as surface density increases within the extended clouds of molecular gas that potentially fuel clumps of star formation in $z\sim1$ galaxies, similar to those observed in local ULIRGs. Instead, our results suggest that molecular gas, both at $z\sim0$ and $z\sim1$, is primarily in the form of self-gravitating molecular clouds. While CO observations suggest a factor of $\sim3$ reduction in the average molecular gas depletion time between $z \sim 0$ and $z\sim1$, we find that, for typical galaxies, the structure of molecular gas and the process of star formation at $z \sim 1$ is otherwise remarkably similar to that observed in local star-forming systems.
[65]  [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.
[66]  [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.
[67]  [pdf] - 1532450
The Spatially Resolved [CII] Cooling Line Deficit in Galaxies
Comments: 17 pages, 5 figures, accepted for publication in ApJ
Submitted: 2016-11-04
We present [CII] 158um measurements from over 15,000 resolved regions within 54 nearby galaxies of the KINGFISH program to investigate the so-called [CII] "line cooling deficit" long known to occur in galaxies with different luminosities. The [CII]/TIR ratio ranges from above 1% to below 0.1% in the sample, with a mean value of 0.48+-0.21%. We find that the surface density of 24um emission dominates this trend, with [CII]/TIR dropping as nuInu{24um} increases. Deviations from this overall decline are correlated with changes in the gas phase metal abundance, with higher metallicity associated with deeper deficits at a fixed surface brightness. We supplement the local sample with resolved [CII] measurements from nearby luminous infrared galaxies and high redshift sources from z=1.8-6.4, and find that star formation rate density drives a continuous trend of deepening [CII] deficit across six orders of magnitude in SFRD. The tightness of this correlation suggests that an approximate star formation rate density can be estimated directly from global measurements of [CII]/TIR, and a relation is provided to do so. Several low-luminosity AGN hosts in the sample show additional and significant central suppression of [CII]/TIR, but these deficit enhancements occur not in those AGN with the highest X-ray luminosities, but instead those with the highest central starlight intensities. Taken together, these results demonstrate that the [CII] cooling line deficit in galaxies likely arises from local physical phenomena in interstellar gas.
[68]  [pdf] - 1505088
The Case for a Publicly Available, Well-Instrumented GBT Operating at 20-115 GHz
Submitted: 2016-10-27, last modified: 2016-10-31
A well-instrumented Green Bank Telescope (GBT) operating at high frequency represents a unique scientific resource for the US community. As a filled-aperture, 100m-diameter telescope, the GBT is ideally suited to fast mapping of extended, low surface brightness emission with excellent instantaneous frequency coverage. This capability makes the GBT a key facility for a range of cutting edge science described in this document, only possible at these frequencies. We note that the ability to perform the necessary observations is unique and highly complementary to the capabilities offered by interferometers, and should be preserved. We argue that rather than divesting from this exceptional resource, it makes sense for the US community to invest moderately to maintain GBT operations and to instrument it in an optimal manner, enabling it to become an extraordinary complement to existing and future radio interferometers. Adequately instrumented, the GBT would be a pillar for 20-115 GHz science in the US and the world.
[69]  [pdf] - 1528221
The Molecular Wind in the Nearest Seyfert Galaxy Circinus Revealed by ALMA
Comments: Accepted to ApJ, 18 pages, 15 figures
Submitted: 2016-09-20
We present ALMA observations of the inner 1' (1.2 kpc) of the Circinus galaxy, the nearest Seyfert. We target CO (1-0) in the region associated with a well-known multiphase outflow driven by the central active galactic nucleus (AGN). While the geometry of Circinus and its outflow make disentangling the latter difficult, we see indications of outflowing molecular gas at velocities consistent with the ionized outflow. We constrain the mass of the outflowing molecular gas to be 1.5e5 -5.1e6 solar masses, yielding a molecular outflow rate of 0.35-12.3 solar masses per year. The values within this range are comparable to the star formation rate in Circinus, indicating that the outflow indeed regulates star formation to some degree. The molecular outflow in Circinus is considerably lower in mass and energetics than previously-studied AGN-driven outflows, especially given its high ratio of AGN luminosity to bolometric luminosity. The molecular outflow in Circinus is, however, consistent with some trends put forth in Cicone et al. (2014), including a linear relation between kinetic power and AGN luminosity, as well as its momentum rate vs. bolometric luminosity (although the latter places Circinus among the starburst galaxies in that sample). We detect additional molecular species including CN and C17O.
[70]  [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.
[71]  [pdf] - 1542668
The Ionized Gas in Nearby Galaxies as Traced by the [NII] 122 and 205 \mu m Transitions
Comments: 18 pages, 9 figures, accepted for publication in The Astrophysical Journal
Submitted: 2016-05-10
The [NII] 122 and 205 \mu m transitions are powerful tracers of the ionized gas in the interstellar medium. By combining data from 21 galaxies selected from the Herschel KINGFISH and Beyond the Peak surveys, we have compiled 141 spatially resolved regions with a typical size of ~1 kiloparsec, with observations of both [NII] far-infrared lines. We measure [NII] 122/205 line ratios in the ~0.6-6 range, which corresponds to electron gas densities $n_e$~1-300 cm$^{-3}$, with a median value of $n_e$=30 cm$^{-3}$. Variations in the electron density within individual galaxies can be as a high as a factor of ~50, frequently with strong radial gradients. We find that $n_e$ increases as a function of infrared color, dust-weighted mean starlight intensity, and star formation rate surface density ($\Sigma_{SFR}$). As the intensity of the [NII] transitions is related to the ionizing photon flux, we investigate their reliability as tracers of the star formation rate (SFR). We derive relations between the [NII] emission and SFR in the low-density limit and in the case of a log-normal distribution of densities. The scatter in the correlation between [NII] surface brightness and $\Sigma_{SFR}$ can be understood as a property of the $n_e$ distribution. For regions with $n_e$ close to or higher than the [NII] line critical densities, the low-density limit [NII]-based SFR calibration systematically underestimates the SFR since [NII] emission is collisionally quenched. Finally, we investigate the relation between [NII] emission, SFR, and $n_e$ by comparing our observations to predictions from the MAPPINGS-III code.
[72]  [pdf] - 1457258
Comparing [CII], HI, and CO dynamics of nearby galaxies
Comments: Accepted for publication in the Astronomical Journal
Submitted: 2016-04-29
The HI and CO components of the interstellar medium (ISM) are usually used to derive the dynamical mass M_dyn of nearby galaxies. Both components become too faint to be used as a tracer in observations of high-redshift galaxies. In those cases, the 158 $\mu$m line of atomic carbon [CII] may be the only way to derive M_dyn. As the distribution and kinematics of the ISM tracer affects the determination of M_dyn, it is important to quantify the relative distributions of HI, CO and [CII]. HI and CO are well-characterised observationally, however, for [CII] only very few measurements exist. Here we compare observations of CO, HI, and [CII] emission of a sample of nearby galaxies, drawn from the HERACLES, THINGS and KINGFISH surveys. We find that within R_25, the average [CII] exponential radial profile is slightly shallower than that of the CO, but much steeper than the HI distribution. This is also reflected in the integrated spectrum ("global profile"), where the [CII] spectrum looks more like that of the CO than that of the HI. For one galaxy, a spectrally resolved comparison of integrated spectra was possible; other comparisons were limited by the intrinsic line-widths of the galaxies and the coarse velocity resolution of the [CII] data. Using high-spectral-resolution SOFIA [CII] data of a number of star forming regions in two nearby galaxies, we find that their [CII] linewidths agree better with those of the CO than the HI. As the radial extent of a given ISM tracer is a key input in deriving M_dyn from spatially unresolved data, we conclude that the relevant length-scale to use in determining M_dyn based on [CII] data, is that of the well-characterised CO distribution. This length scale is similar to that of the optical disk.
[73]  [pdf] - 1432909
The Relationship Between Molecular Gas, HI, and Star Formation in the Low-Mass, Low-Metallicity Magellanic Clouds
Comments: 24 pages, 14 figures, accepted for publication in ApJ. FITS files of the dust-based estimates of the H2 column densities for the LMC and SMC (shown in Figures 2 and 3) will be available online through ApJ
Submitted: 2015-10-27, last modified: 2016-04-25
The Magellanic Clouds provide the only laboratory to study the effect of metallicity and galaxy mass on molecular gas and star formation at high (~20 pc) resolution. We use the dust emission from HERITAGE Herschel data to map the molecular gas in the Magellanic Clouds, avoiding the known biases of CO emission as a tracer of H$_{2}$. Using our dust-based molecular gas estimates, we find molecular gas depletion times of ~0.4 Gyr in the LMC and ~0.6 SMC at 1 kpc scales. These depletion times fall within the range found for normal disk galaxies, but are shorter than the average value, which could be due to recent bursts in star formation. We find no evidence for a strong intrinsic dependence of the molecular gas depletion time on metallicity. We study the relationship between gas and star formation rate across a range in size scales from 20 pc to ~1 kpc, including how the scatter in molecular gas depletion time changes with size scale, and discuss the physical mechanisms driving the relationships. We compare the metallicity-dependent star formation models of Ostriker, McKee, and Leroy (2010) and Krumholz (2013) to our observations and find that they both predict the trend in the data, suggesting that the inclusion of a diffuse neutral medium is important at lower metallicity.
[74]  [pdf] - 1411511
The EMPIRE Survey: Systematic Variations in the Dense Gas Fraction and Star Formation Efficiency from Full-Disk Mapping of M51
Comments: 7 pages, 5 figures, ApJL accepted
Submitted: 2016-04-18
We present the first results from the EMPIRE survey, an IRAM large program that is mapping tracers of high density molecular gas across the disks of nine nearby star-forming galaxies. Here, we present new maps of the 3-mm transitions of HCN, HCO+, and HNC across the whole disk of our pilot target, M51. As expected, dense gas correlates with tracers of recent star formation, filling the "luminosity gap" between Galactic cores and whole galaxies. In detail, we show that both the fraction of gas that is dense, f_dense traced by HCN/CO, and the rate at which dense gas forms stars, SFE_dense traced by IR/HCN, depend on environment in the galaxy. The sense of the dependence is that high surface density, high molecular gas fraction regions of the galaxy show high dense gas fractions and low dense gas star formation efficiencies. This agrees with recent results for individual pointings by Usero et al. 2015 but using unbiased whole-galaxy maps. It also agrees qualitatively with the behavior observed contrasting our own Solar Neighborhood with the central regions of the Milky Way. The sense of the trends can be explained if the dense gas fraction tracks interstellar pressure but star formation occurs only in regions of high density contrast.
[75]  [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.
[76]  [pdf] - 1415023
Using CO to Measure Molecular Masses
Comments: Proceedings for invited presentation to the 6th Zermatt Symposium, "Conditions and Impact of Star Formation: from Lab to Space", Zermatt Switzerland September 7-11, 2015
Submitted: 2015-12-07
With an increased appreciation for the role of gas in galaxy evolution, there is renewed interest in measuring gas masses for galaxies. I review some of the basic concepts in using CO to determine molecular masses, and discuss some of the recent work.
[77]  [pdf] - 1347678
The dust properties and physical conditions of the interstellar medium in the LMC massive star forming complex N11
Comments: 24 pages, 16 figures, accepted for publication in MNRAS
Submitted: 2015-11-23, last modified: 2015-11-25
We combine Spitzer and Herschel data of the star-forming region N11 in the Large Magellanic Cloud to produce detailed maps of the dust properties in the complex and study their variations with the ISM conditions. We also compare APEX/LABOCA 870um observations with our model predictions in order to decompose the 870um emission into dust and non-dust (free-free emission and CO(3-2) line) contributions. We find that in N11, the 870um can be fully accounted for by these 3 components. The dust surface density map of N11 is combined with HI and CO observations to study local variations in the gas-to-dust mass ratios. Our analysis leads to values lower than those expected from the LMC low-metallicity as well as to a decrease of the gas-to-dust mass ratio with the dust surface density. We explore potential hypotheses that could explain the low observed gas-to-dust mass ratios (variations in the XCO factor, presence of CO-dark gas or of optically thick HI or variations in the dust abundance in the dense regions). We finally decompose the local SEDs using a Principal Component Analysis (i.e. with no a priori assumption on the dust composition in the complex). Our results lead to a promising decomposition of the local SEDs in various dust components (hot, warm, cold) coherent with that expected for the region. Further analysis on a larger sample of galaxies will follow in order to understand how unique this decomposition is or how it evolves from one environment to another.
[78]  [pdf] - 1530359
On the Variation of Gas Depletion Time
Comments: 4 pages, 2 figures, Proceedings of IAU Symposium No. 315, 2015
Submitted: 2015-11-09
We present an updated status of the EDGE project, which is a survey of 125 local galaxies in the $^{12}$CO($1-0$) and $^{13}$CO($1-0$) lines. We combine the molecular data of the EDGE survey with the stellar and ionized gas maps of the CALIFA survey to give a comprehensive view of the dependence of the star formation efficiency, or equivalently, the molecular gas depletion time, on various local environments, such as the stellar surface density, metallicity, and radius from the galaxy center. This study will provide insight into the parameters that drive the star formation efficiency in galaxies at $z \sim 0$.
[79]  [pdf] - 1292426
CARMA CO Observations of Three Extremely Metal-Poor, Star-Forming Galaxies
Comments: Accepted for publication in ApJ
Submitted: 2015-10-12
We present sensitive CO (J = 1 - 0) emission line observations of three metal-poor dwarf irregular galaxies Leo P (Z ~ 3% Z_Solar), Sextans A (Z ~ 7.5% Z_Solar), and Sextans B (Z ~ 7.5% Z_Solar), all obtained with the Combined Array for Millimeter-wave Astronomy (CARMA) interferometer. While no CO emission was detected, the proximity of the three systems allows us to place very stringent (4 sigma) upper limits on the CO luminosity (L_CO) in these metal-poor galaxies. We find the CO luminosities to be L_CO < 2900 K km/s pc^2 for Leo P, L_CO < 12400 K km/s pc^2 for Sextans A, and L_CO < 9700 K km/s pc^2 for Sextans B. Comparison of our results with recent observational estimates of the factor for converting between L_CO and the mass of molecular hydrogen, as well as theoretical models, provides further evidence that either the CO-to-H_2 conversion factor increases sharply as metallicity decreases, or that stars are forming in these three galaxies very efficiently, requiring little molecular hydrogen.
[80]  [pdf] - 1275020
The Multi-Phase Cold Fountain in M82 Revealed by a Wide, Sensitive Map of the Molecular ISM
Comments: 36 pages including 25 figures and 4 appendices. Accepted for publication in the Astrophysical Journal. For a high resolution version please visit
Submitted: 2015-09-09
We present a wide area (~ 8 x 8 kpc), sensitive map of CO (2-1) emission around the nearby starburst galaxy M82. Molecular gas extends far beyond the stellar disk, including emission associated with the well-known outflow as far as 3 kpc from M82's midplane. Kinematic signatures of the outflow are visible in both the CO and HI emission: both tracers show a minor axis velocity gradient and together they show double peaked profiles, consistent with a hot outflow bounded by a cone made of a mix of atomic and molecular gas. Combining our CO and HI data with observations of the dust continuum, we study the changing properties of the cold outflow as it leaves the disk. While H_2 dominates the ISM near the disk, the dominant phase of the cool medium changes as it leaves the galaxy and becomes mostly atomic after about a kpc. Several arguments suggest that regardless of phase, the mass in the cold outflow does not make it far from the disk; the mass flux through surfaces above the disk appears to decline with a projected scale length of ~ 1-2 kpc. The cool material must also end up distributed over a much wider angle than the hot outflow based on the nearly circular isophotes of dust and CO at low intensity and the declining rotation velocities as a function of height from the plane. The minor axis of M82 appears so striking at many wavelengths because the interface between the hot wind cavity and the cool gas produces Halpha, hot dust, PAH emission, and scattered UV light. We also show the level at which a face-on version of M82 would be detectable as an outflow based on unresolved spectroscopy. Finally, we consider multiple constraints on the CO-to-H$_2$ conversion factor, which must change across the galaxy but appears to be only a factor of ~ 2 lower than the Galactic value in the outflow.
[81]  [pdf] - 1280998
Faint CO Line Wings in Four Star-Forming (Ultra)luminous Infrared Galaxies
Comments: 9 pages, 5 figures, 1 appendix; accepted for publication in the Astrophysical Journal
Submitted: 2015-08-13
We report the results of a search for large velocity width, low-intensity line wings - a commonly used signature of molecular outflows - in four low redshift (ultra)luminous infrared galaxies (U/LIRGs) that appear to be dominated by star formation. The targets were drawn from a sample of fourteen such galaxies presented in Chung et al. (2011), who showed the stacked CO spectrum of the sample to exhibit 1000 km/s-wide line wings. We obtained sensitive, wide bandwidth imaging of our targets using the IRAM Plateau de Bure Interferometer. We detect each target at very high significance but do not find the claimed line wings in these four targets. Instead, we constrain the flux in the line wings to be only a few percent. Casting our results as mass outflow rates following Cicone et al. (2014) we show them to be consistent with a picture in which very high mass loading factors preferentially occur in systems with high AGN contributions to their bolometric luminosity. We identify one of our targets, IRAS05083 (VII Zw 31), as a candidate molecular outflow.
[82]  [pdf] - 1273232
High Resolution Imaging of PHIBSS z~2 Main Sequence Galaxies in CO J=1-0
Comments: To appear in ApJ
Submitted: 2015-07-20, last modified: 2015-07-22
We present Karl G. Jansky Very Large Array observations of the CO J=1-0 transition in a sample of four $z\sim2$ main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies at their redshift, and are part of the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) which imaged them in CO J=3-2. Two galaxies are imaged here at high signal-to-noise, allowing determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using these and published measurements, we show that the CO and optical disks have similar sizes in main-sequence galaxies, and in the galaxy where we can compare CO J=1-0 and J=3-2 sizes we find these are also very similar. Assuming a Galactic CO-to-H$_2$ conversion, we measure surface densities of $\Sigma_{mol}\sim1200$ M$_\odot$pc$^{-2}$ in projection and estimate $\Sigma_{mol}\sim500-900$ M$_\odot$pc$^{-2}$ deprojected. Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios $r_{31}$ that are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J=1-0 and J=3-2 indicate that both transitions sample the same kinematics, implying that their emission is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active star-formation. We suggest that $r_{31}$ in very actively star-forming galaxies is likely an indicator of how well mixed the star formation activity and the molecular reservoir are.
[83]  [pdf] - 1045078
A Cryogenic Space Telescope for Far-Infrared Astrophysics: A Vision for NASA in the 2020 Decade
Comments: This document contains a large paper on CALISTO, as well as 6 small scientific white papers, all prepared in response to NASA's Cosmic Origins Program Analysis Group (PAG) request for input on future mission concepts
Submitted: 2015-05-20
Many of the transformative processes in the Universe have taken place in regions obscured by dust, and are best studied with far-IR spectroscopy. We present the Cryogenic-Aperture Large Infrared-Submillimeter Telescope Observatory (CALISTO), a 5-meter class, space-borne telescope actively cooled to 4 K, emphasizing moderate-resolution spectroscopy in the crucial 35 to 600 micron band. CALISTO will enable NASA and the world to study the rise of heavy elements in the Universe's first billion years, chart star formation and black hole growth in dust-obscured galaxies through cosmic time, and conduct a census of forming planetary systems in our region of the Galaxy. CALISTO will capitalize on rapid progress in both format and sensitivity of far-IR detectors. Arrays with a total count of a few 100,000 detector pixels will form the heart of a suite of imaging spectrometers in which each detector reaches the photon background limit. This document contains a large overview paper on CALISTO, as well as six 2-3 page scientific white papers, all prepared in response to NASA's Cosmic Origins Program Analysis Group (COPAG's) request for input on future mission concepts. The Far-IR Science Interest Group will meet from 3-5 June 2015 with the intention of reaching consensus on the architecture for the Far-IR Surveyor mission. This white paper describes one of the architectures to be considered by the community. One or more companion papers will describe alternative architectures.
[84]  [pdf] - 1043211
The Relationship between CO Emission and Visual Extinction Traced by Dust Emission in the Magellanic Clouds
Comments: 20 pages, 11 figures, submitted to MNRAS
Submitted: 2015-04-10
To test the theoretical understanding that finding bright CO emission depends primarily on dust shielding, we investigate the relationship between CO emission ($I_{\rm CO}$) and the amount of dust (estimated from IR emission and expressed as "$A_V$") across the Large Magellanic Cloud, the Small Magellanic Cloud, and the Milky Way. We show that at our common resolution of 10 pc scales, $I_{\rm CO}$ given a fixed line-of-sight $A_V$ is similar across all three systems despite the difference in metallicity. We find some evidence for a secondary dependence of $I_{\rm CO}$ on radiation field; in the LMC, $I_{\rm CO}$ at a given $A_V$ is smaller in regions of high $T_{\rm dust}$, perhaps because of an increased photodissociating radiation field. We suggest a simple but useful picture in which the CO-to-H$_2$ conversion factor (\xco) depends on two separable factors: (1) the distribution of gas column densities, which maps to an extinction distribution via a dust-to-gas ratio; and (2) the dependence of $I_{\rm CO}$ on $A_V$. Assuming that the probability distribution function (PDF) of local Milky Way clouds is universal, this approach predicts a dependence of \xco\ on $Z$ between $Z^{-1}$ and $Z^{-2}$ above about a third solar metallicity. Below this metallicity, CO emerges from only the high column density parts of the cloud and so depends very sensitively on the adopted PDF and the H$_2$/{\sc Hi} prescription. The PDF of low metallicity clouds is thus of considerable interest and the uncertainty associated with even an ideal prescription for \xco\ at very low metallicity will be large.
[85]  [pdf] - 1222697
The influence of supernova remnants on the interstellar medium in the Large Magellanic Cloud seen at 20--600 $\mu$m wavelengths
Comments: Accepted to Astrophysical Journal
Submitted: 2014-10-21, last modified: 2015-03-06
We present the analysis of supernova remnants (SNRs) in the Large Magellanic Cloud (LMC) and their influence on the environment at far-infrared (FIR) and submillimeter wavelengths. We use new observations obtained with the {\it Herschel} Space Observatory and archival data obtained with the {\it Spitzer} Space Telescope, to make the first FIR atlas of these objects. The SNRs are not clearly discernible at FIR wavelengths, however their influence becomes apparent in maps of dust mass and dust temperature, which we constructed by fitting a modified black-body to the observed spectral energy distribution in each sightline. Most of the dust that is seen is pre-existing interstellar dust in which SNRs leave imprints. The temperature maps clearly reveal SNRs heating surrounding dust, while the mass maps indicate the removal of 3.7$^{+7.5}_{-2.5}$ M$_{\odot}$ of dust per SNR. This agrees with the calculations by others that significant amounts of dust are sputtered by SNRs. Under the assumption that dust is sputtered and not merely pushed away, we estimate a dust destruction rate in the LMC of $0.037^{+0.075}_{-0.025}$ M$_\odot$ yr$^{-1}$ due to SNRs, yielding an average lifetime for interstellar dust of $2^{+4.0}_{-1.3}\times10^7$ yr. We conclude that sputtering of dust by SNRs may be an important ingredient in models of galactic evolution, that supernovae may destroy more dust than they produce, and that they therefore may not be net producers of long lived dust in galaxies.
[86]  [pdf] - 1224139
Spatially extended and high-velocity dispersion molecular component in spiral galaxies: single-dish vs. interferometric observations
Comments: 12 pages, 5 figures, Accepted to AJ
Submitted: 2015-01-22
Recent studies of the molecular medium in nearby galaxies have provided mounting evidence that the molecular gas can exist in two phases: one that is clumpy and organized as molecular clouds and another one that is more diffuse. This last component has a higher velocity dispersion than the clumpy one. In order to investigate these two molecular components further, we compare the fluxes and line widths of CO in NGC 4736 and NGC 5055, two nearby spiral galaxies for which high-quality interferometric as well as single-dish data sets are available. Our analysis leads to two main results: 1) Employing three different methods, we determine the flux recovery of the interferometer as compared to the single-dish to be within a range of 35-74% for NGC4736 and 81-92% for NGC5055, and 2) when focusing on high (SNR>5) lines of sight, the single-dish line widths are larger by ~(40+-20)% than the ones derived from interferometric data; which is in agreement with stacking all lines of sight. These results point to a molecular gas component that is distributed over spatial scales larger than 30"(~1kpc), and is therefore filtered out by the interferometer. The available observations do not allow us to distinguish between a truly diffuse gas morphology and a uniform distribution of small clouds that are separated by less than the synthesized beam size (~3" or ~100pc), as they would both be invisible for the interferometer. This high velocity dispersion component has a dispersion similar to what is found in the atomic medium, as traced through observations of the HI line.
[87]  [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.
[88]  [pdf] - 1216669
Combined CO & Dust Scaling Relations of Depletion Time and Molecular Gas Fractions with Cosmic Time, Specific Star Formation Rate and Stellar Mass
Comments: ApJ accepted
Submitted: 2014-09-03, last modified: 2014-12-05
We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.
[89]  [pdf] - 1223213
Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. II. Gas-to-Dust Ratio Variations across ISM Phases
Submitted: 2014-11-17
The spatial variations of the gas-to-dust ratio (GDR) provide constraints on the chemical evolution and lifecycle of dust in galaxies. We examine the relation between dust and gas at 10-50 pc resolution in the Large and Small Magellanic Clouds (LMC and SMC) based on Herschel far-infrared (FIR), H I 21 cm, CO, and Halpha observations. In the diffuse atomic ISM, we derive the gas-to-dust ratio as the slope of the dust-gas relation and find gas-to-dust ratios of 380+250-130 in the LMC, and 1200+1600-420 in the SMC, not including helium. The atomic-to-molecular transition is located at dust surface densities of 0.05 Mo pc-2 in the LMC and 0.03 Mo pc-2 in the SMC, corresponding to AV ~ 0.4 and 0.2, respectively. We investigate the range of CO-to-H2 conversion factor to best account for all the molecular gas in the beam of the observations, and find upper limits on XCO to be 6x1020 cm-2 K-1 km-1 s in the LMC (Z=0.5Zo) at 15 pc resolution, and 4x 1021 cm-2 K-1 km-1 s in the SMC (Z=0.2Zo) at 45 pc resolution. In the LMC, the slope of the dust-gas relation in the dense ISM is lower than in the diffuse ISM by a factor ~2, even after accounting for the effects of CO-dark H2 in the translucent envelopes of molecular clouds. Coagulation of dust grains and the subsequent dust emissivity increase in molecular clouds, and/or accretion of gas-phase metals onto dust grains, and the subsequent dust abundance (dust-to-gas ratio) increase in molecular clouds could explain the observations. In the SMC, variations in the dust-gas slope caused by coagulation or accretion are degenerate with the effects of CO-dark H2. Within the expected 5--20 times Galactic XCO range, the dust-gas slope can be either constant or decrease by a factor of several across ISM phases. Further modeling and observations are required to break the degeneracy between dust grain coagulation, accretion, and CO-dark H2.
[90]  [pdf] - 1223124
Dust continuum emission as a tracer of gas mass in galaxies
Comments: Accepted for publication in ApJ
Submitted: 2014-11-11
We use a sample of 36 galaxies from the KINGFISH (Herschel IR), HERACLES (IRAM CO), and THINGS (VLA HI) surveys to study empirical relations between Herschel infrared (IR) luminosities and the total mass of the interstellar gas (H2+HI). Such a comparison provides a simple empirical relationship without introducing the uncertainty of dust model fitting. We find tight correlations, and provide fits to these relations, between Herschel luminosities and the total gas mass integrated over entire galaxies, with the tightest, almost linear, correlation found for the longest wavelength data (SPIRE500). However, we find that accounting for the gas-phase metallicity (affecting the dust-to-gas ratio) is crucial when applying these relations to low-mass, and presumably high-redshift, galaxies. The molecular (H2) gas mass is found to be better correlated with the peak of the IR emission (e.g. PACS160), driven mostly by the correlation of stellar mass and mean dust temperature. When examining these relations as a function of galactocentric radius we find the same correlations, albeit with a larger scatter, up to a radius of 0.7 r_25 (within which most of the galaxy's baryonic mass resides). However, beyond this radius the same correlations no longer hold, with the gas mass (predominantly HI) increasing relative to the infrared emission. The tight relations found for the bulk of the galaxy's baryonic content suggest that the total gas masses of disk-like (non-merging) galaxies can be inferred from far-infrared continuum measurements in situations where only the latter are available, e.g. in ALMA continuum observations of high-redshift galaxies.
[91]  [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.
[92]  [pdf] - 906411
Extended HCN and HCO$^{+}$ emission in the starburst galaxy M82
Comments: 10 pages, 6 figures, 5 tables. Accepted for publication in the ApJ
Submitted: 2014-10-23
We mapped 3 mm continuum and line emission from the starburst galaxy M82 using the Combined Array for Research in Millimeter-wave Astronomy. We targeted the HCN, HCO$^{+}$, HNC, CS and HC$_{3}$N lines, but here we focus on the HCN and HCO$^{+}$ emission. The map covers a field of 1.2' with a ~5" resolution. The HCN and HCO$^{+}$ observations are combined with single dish images. The molecular gas in M82 had been previously found to be distributed in a molecular disk, coincident with the central starburst, and a galactic scale outflow which originates in the central starburst. With the new short spacings-corrected maps we derive some of the properties of the dense molecular gas in the base of the outflow. From the HCN and HCO$^{+}$ J=(1-0) line emission, and under the assumptions of the gas being optically thin and in local thermodynamic equilibrium, we place lower limits to the amount of dense molecular gas in the base of the outflow. The lower limits are $7\times10^{6}$ $M_{\odot}$ and $21\times10^{6}$ $M_{\odot}$, or $\gtrsim2\%$ of the total molecular mass in the outflow. The kinematics and spatial distribution of the dense gas outside the central starburst suggests that it is being expelled through chimneys. Assuming a constant outflow velocity, the derived outflow rate of dense molecular gas is $\geq0.3$ $M_{\odot}$ yr$^{-1}$, which would lower the starburst lifetime by $\geq5\%$. The energy required to expel this mass of dense gas is $(1-10)\times10^{52}$ erg.
[93]  [pdf] - 2038171
The CARMA Paired Antenna Calibration System: Atmospheric Phase Correction for Millimeter Wave Interferometry and its Application to Mapping the Ultraluminous Galaxy Arp 193
Comments: 19 pages, 18 figures
Submitted: 2014-10-21
Phase fluctuations introduced by the atmosphere are the main limiting factor in attaining diffraction limited performance in extended interferometric arrays at millimeter and submillimeter wavelengths. We report the results of C-PACS, the Combined Array for Research in Millimeter-Wave Astronomy Paired Antenna Calibration System. We present a systematic study of several hundred test observations taken during the 2009-2010 winter observing season where we utilize CARMA's eight 3.5-m antennas to monitor an atmospheric calibrator while simultaneously acquiring science observations with 6.1-m and 10.4-m antennas on baselines ranging from a few hundred meters to ~2 km. We find that C-PACS is systematically successful at improving coherence on long baselines under a variety of atmospheric conditions. We find that the angular separation between the atmospheric calibrator and target source is the most important consideration, with consistently successful phase correction at CARMA requiring a suitable calibrator located $\lesssim$6$^\circ$ away from the science target. We show that cloud cover does not affect the success of C-PACS. We demonstrate C-PACS in typical use by applying it to the observations of the nearby very luminous infrared galaxy Arp 193 in $^{12}$CO(2-1) at a linear resolution of ~70 pc (0.12" x 0.18"), 3 times better than previously published molecular maps of this galaxy. We resolve the molecular disk rotation kinematics and the molecular gas distribution and measure the gas surface densities and masses on 90 pc scales. We find that molecular gas constitutes $\sim30\%$ of the dynamical mass in the inner 700 pc of this object with a surface density $\sim10^4 M_\odot$ pc$^{-2}$; we compare these properties to those of the starburst region of NGC 253.
[94]  [pdf] - 1222236
[CII] 158 $\mu$m Emission as a Star Formation Tracer
Comments: 25 pages, 14 figures, accepted for publication in ApJ
Submitted: 2014-09-24
The [CII] 157.74 $\mu$m transition is the dominant coolant of the neutral interstellar gas, and has great potential as a star formation rate (SFR) tracer. Using the Herschel KINGFISH sample of 46 nearby galaxies, we investigate the relation of [CII] surface brightness and luminosity with SFR. We conclude that [CII] can be used for measurements of SFR on both global and kiloparsec scales in normal star-forming galaxies in the absence of strong active galactic nuclei (AGN). The uncertainty of the $\Sigma_{\rm [CII]}-\Sigma_{\rm SFR}$ calibration is $\pm$0.21 dex. The main source of scatter in the correlation is associated with regions that exhibit warm IR colors, and we provide an adjustment based on IR color that reduces the scatter. We show that the color-adjusted $\Sigma_{\rm[CII]}-\Sigma_{\rm SFR}$ correlation is valid over almost 5 orders of magnitude in $\Sigma_{\rm SFR}$, holding for both normal star-forming galaxies and non-AGN luminous infrared galaxies. Using [CII] luminosity instead of surface brightness to estimate SFR suffers from worse systematics, frequently underpredicting SFR in luminous infrared galaxies even after IR color adjustment (although this depends on the SFR measure employed). We suspect that surface brightness relations are better behaved than the luminosity relations because the former are more closely related to the local far-UV field strength, most likely the main parameter controlling the efficiency of the conversion of far-UV radiation into gas heating. A simple model based on Starburst99 population-synthesis code to connect SFR to [CII] finds that heating efficiencies are $1\%-3\%$ in normal galaxies.
[95]  [pdf] - 969777
Cool dust heating and temperature mixing in nearby star-forming galaxies
Comments: 28 pages, 20 figures, accepted for publication in A&A
Submitted: 2014-09-20
Physical conditions of the interstellar medium in galaxies are closely linked to the ambient radiation field and the heating of dust grains. In order to characterize dust properties in galaxies over a wide range of physical conditions, we present here the radial surface brightness profiles of the entire sample of 61 galaxies from Key Insights into Nearby Galaxies: Far-Infrared Survey with Herschel (KINGFISH). The main goal of our work is the characterization of the grain emissivities, dust temperatures, and interstellar radiation fields responsible for heating the dust. After fitting the dust and stellar radial profiles with exponential functions, we fit the far-infrared spectral energy distribution (SED) in each annular region with single-temperature modified black bodies using both variable (MBBV) and fixed (MBBF) emissivity indices beta, as well as with physically motivated dust models. Results show that while most SED parameters decrease with radius, the emissivity index beta also decreases with radius in some galaxies, but in others is increasing, or rising in the inner regions and falling in the outer ones. Despite the fixed grain emissivity (average beta~ 2.1) of the physically-motivated models, they are well able to accommodate flat spectral slopes with beta<= 1. We find that flatter slopes (beta<= 1.5) are associated with cooler temperatures, contrary to what would be expected from the usual Tdust-beta degeneracy. This trend is related to variations in Umin since beta and Umin are very closely linked over the entire range in Umin sampled by the KINGFISH galaxies: low Umin is associated with flat beta<=1. Both these results strongly suggest that the low apparent \beta values (flat slopes) in MBBV fits are caused by temperature mixing along the line-of-sight, rather than by intrinsic variations in grain properties. Abstract truncated for arXiv.
[96]  [pdf] - 1179398
The Role of Stellar Feedback in the Dynamics of HII Regions
Comments: 19 pages, 8 figures; accepted by ApJ
Submitted: 2013-09-20, last modified: 2014-09-15
Stellar feedback is often cited as the biggest uncertainty in galaxy formation models today. This uncertainty stems from a dearth of observational constraints as well as the great dynamic range between the small scales (<1 pc) where the feedback occurs and the large scales of galaxies (>1 kpc) that are shaped by this feedback. To bridge this divide, in this paper we aim to assess observationally the role of stellar feedback at the intermediate scales of HII regions. In particular, we employ multiwavelength data to examine several stellar feedback mechanisms in a sample of 32 HII regions in the Large and Small Magellanic Clouds (LMC and SMC, respectively). Using optical, infrared, radio, and X-ray images, we measure the pressures exerted on the shells from the direct stellar radiation, the dust-processed radiation, the warm ionized gas, and the hot X-ray emitting gas. We find that the warm ionized gas dominates over the other terms in all of the sources, although two have comparable dust-processed radiation pressures to their warm gas pressures. The hot gas pressures are comparatively weak, while the direct radiation pressures are 1-2 orders of magnitude below the other terms. We discuss the implications of these results, particularly highlighting evidence for hot gas leakage from the HII shells and regarding the momentum deposition from the dust-processed radiation to the warm gas. Furthermore, we emphasize that similar observational work should be done on very young HII regions to test whether direct radiation pressure and hot gas can drive the dynamics at early times.
[97]  [pdf] - 1215169
Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm Excess Emission
Comments: 20 pages, 9 figures, accepted for publication in the ApJ
Submitted: 2014-06-23
The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.
[98]  [pdf] - 1209803
Extreme gas fractions in clumpy, turbulent disk galaxies at z~0.1
Comments: Accepted to ApJ Letters
Submitted: 2014-05-28
In this letter we report the discovery of CO fluxes, suggesting very high gas fractions in three disk galaxies seen in the nearby Universe (z ~ 0.1). These galaxies were investigated as part of the DYnamics of Newly Assembled Massive Objects (DYNAMO) survey. High-resolution Hubble Space Telescope imaging of these objects reveals the presence of large star forming clumps in the bodies of the galaxies, while spatially resolved spectroscopy of redshifted Halpha reveals the presence of high dispersion rotating disks. The internal dynamical state of these galaxies resembles that of disk systems seen at much higher redshifts (1 < z < 3). Using CO(1-0) observations made with the Plateau de Bure Interferometer, we find gas fractions of 20-30% and depletion times of tdep ~ 0.5 Gyr (assuming a Milky Way-like CO conversion factor). These properties are unlike those expected for low- redshift galaxies of comparable specific star formation rate, but they are normal for their high-z counterparts. DYNAMO galaxies break the degeneracy between gas fraction and redshift, and we show that the depletion time per specific star formation rate for galaxies is closely tied to gas fraction, independent of redshift. We also show that the gas dynamics of two of our local targets corresponds to those expected from unstable disks, again resembling the dynamics of high-z disks. These results provide evidence that DYNAMO galaxies are local analogues to the clumpy, turbulent disks, which are often found at high redshift.
[99]  [pdf] - 1208288
A far-IR view of the starburst driven superwind in NGC 2146
Comments: 16 pages, 13 figures, accepted for publication in ApJ
Submitted: 2014-03-10, last modified: 2014-05-05
NGC 2146, a nearby luminous infrared galaxy (LIRG), presents evidence for outflows along the disk minor axis in all gas phases (ionized, neutral atomic and molecular). We present an analysis of the multi-phase starburst driven superwind in the central 5 kpc as traced in spatially resolved spectral line observations, using far-IR Herschel PACS spectroscopy, to probe the effects on the atomic and ionized gas, and optical integral field spectroscopy to examine the ionized gas through diagnostic line ratios. We observe an increased ~250 km/s velocity dispersion in the [OI] 63 micron, [OIII] 88 micron, [NII] 122 micron and [CII] 158 micron fine-structure lines that is spatially coincident with high excitation gas above and below the disk. We model this with a slow ~200 km/s shock and trace the superwind to the edge of our field of view 2.5 kpc above the disk. We present new SOFIA 37 micron observations to explore the warm dust distribution, and detect no clear dust entrainment in the outflow. The stellar kinematics appear decoupled from the regular disk rotation seen in all gas phases, consistent with a recent merger event disrupting the system. We consider the role of the superwind in the evolution of NGC 2146 and speculate on the evolutionary future of the system. Our observations of NGC 2146 in the far-IR allow an unobscured view of the wind, crucial for tracing the superwind to the launching region at the disk center, and provide a local analog for future ALMA observations of outflows in high redshift systems.
[100]  [pdf] - 1209074
ALFALFA Discovery of the Nearby Gas-rich Dwarf Galaxy Leo P. V. Neutral Gas Dynamics and Kinematics
Comments: 38 pages, 11 figures, Accepted for publication in the Astronomical Journal
Submitted: 2014-04-21
We present new HI spectral line imaging of the extremely metal-poor, star-forming dwarf irregular galaxy Leo P. Our HI images probe the global neutral gas properties and the local conditions of the interstellar medium (ISM). The HI morphology is slightly elongated along the optical major-axis. We do not find obvious signatures of interaction or infalling gas at large spatial scales. The neutral gas disk shows obvious rotation, although the velocity dispersion is comparable to the rotation velocity. The rotation amplitude is estimated to be V_c = 15 +/- 5 km/s. Within the HI radius probed by these observations, the mass ratio of gas to stars is roughly 2:1, while the ratio of the total mass to the baryonic mass is ~15:1. We use this information to place Leo P on the baryonic Tully-Fisher relation, testing the baryonic content of cosmic structures in a sparsely populated portion of parameter space that has hitherto been occupied primarily by dwarf spheroidal galaxies. We detect the signature of two temperature components in the neutral ISM of Leo P; the cold and warm components have characteristic velocity widths of 4.2 +/- 0.9 km/s and 10.1 +/- 1.2 km/s, corresponding to kinetic temperatures of ~1100 K and ~6200 K, respectively. The cold HI component is unresolved at a physical resolution of 200 pc. The highest HI surface densities are observed in close physical proximity to the single HII region. A comparison of the neutral gas properties of Leo P with other extremely metal-deficient (XMD) galaxies reveals that Leo P has the lowest neutral gas mass of any known XMD, and that the dynamical mass of Leo P is more than two orders of magnitude smaller than any known XMD with comparable metallicity.
[101]  [pdf] - 1202744
Search for [CII] emission in z=6.5-11 star-forming galaxies
Comments: 12 pages, 9 figures, accepted for publication in the Astrophysical Journal
Submitted: 2014-01-14
We present the search for the [CII] emission line in three $z>6.5$ Lyman-alpha emitters (LAEs) and one J-Dropout galaxy using the Combined Array for Research in Millimeter-wave Astronomy (CARMA) and the Plateau de Bure Interferometer (PdBI). We observed three bright $z\sim6.5-7$ LAEs discovered in the SUBARU deep field (SDF) and the Multiple Imaged lensed $z\sim 11$ galaxy candidate found behind the galaxy cluster MACSJ0647.7+7015. For the LAEs IOK-1 ($z=6.965$), SDF J132415.7+273058 ($z=6.541$) and SDF J132408.3+271543 ($z=6.554$) we find upper limits for the [CII] line luminosity of $<2.05$, $<4.52$ and $<10.56\times10^{8}{\rm L}_{\odot}$ respectively. We find upper limits to the FIR luminosity of the galaxies using a spectral energy distribution template of the local galaxy NGC 6946 and taking into account the effects of the Cosmic Microwave Background on the mm observations. For IOK-1, SDF J132415.7+273058 and SDF J132408.3+271543 we find upper limits for the FIR luminosity of $<2.33$, $3.79$ and $7.72\times10^{11}{\rm L}_{\odot}$ respectively. For the lensed galaxy MACS0647-JD, one of the highest redshift galaxy candidate to date with $z_{\rm ph}=10.7^{+0.6}_{-0.4}$ we put an upper limit in the [CII] emission of $<1.36\times10^{8}\times(\mu/15)^{-1}{\rm L}_{\odot}$ and an upper limit in the FIR luminosity of $<6.1\times10^{10}\times(\mu/15)^{-1}{\rm L}_{\odot}$ (where $\mu$ is the magnification factor). We explore the different conditions relevant for the search for [CII] emission in high redshift galaxies as well as the difficulties for future observations with ALMA and CCAT.
[102]  [pdf] - 1202178
Formation of Molecular Clouds and Global Conditions for Star Formation
Comments: 25 pages, 7 figures, Accepted for publication as a chapter in Protostars and Planets VI, University of Arizona Press (2014), eds. H. Beuther, R. Klessen, C. Dullemond, Th. Henning
Submitted: 2013-12-11
Giant molecular clouds (GMCs) are the primary reservoirs of cold, star-forming molecular gas in the Milky Way and similar galaxies, and thus any understanding of star formation must encompass a model for GMC formation, evolution, and destruction. These models are necessarily constrained by measurements of interstellar molecular and atomic gas, and the emergent, newborn stars. Both observations and theory have undergone great advances in recent years, the latter driven largely by improved numerical simulations, and the former by the advent of large-scale surveys with new telescopes and instruments. This chapter offers a thorough review of the current state of the field.
[103]  [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.
[104]  [pdf] - 734730
Dust May Be More Rare Than Expected in Metal Poor Galaxies
Comments: 13 pages, 4 figures, accepted for publication in Nature
Submitted: 2013-10-17
'Normal' galaxies observed at z>6, when the Universe was <1 billion years old, thus far show no evidence of the cold dust that accompanies star formation in the local Universe, where the dust-to-gas mass ratio is 1%. A prototypical example is 'Himiko' (z=6.6), which a mere 840 Myr after the Big Bang is forming stars at a rate of 30-100 Msun/yr, yielding a mass assembly time M^{star}/SFR 150x10^6 yr. Himiko is estimated to have a low fraction (2-3% of the Solar value) of elements heavier than helium (metallicity), and although its gas mass cannot be asserted at this time its dust-to-stellar mass ratio is constrained to be <0.05%. The local galaxy I Zw 18, with a metallicity 4% solar and forming stars less rapidly than Himiko but still vigorously for its mass (M^{star}/SFR 1.6x10^9 yr), is also very dust deficient and perhaps one of the best analogues of primitive galaxies accessible to detailed study. Here we report observations of dust emission from I Zw 18 from which we determine its dust mass to be 450-1800 Msun, yielding a dust-to-stellar mass ratio \approx 10^{-6}-10^{-5} and a dust-to-gas mass ratio 3.2-13x10^{-6}. If I Zw 18 is a reasonable analog of Himiko, then Himiko's dust mass is \approx 50,000 Msun, a factor of 100 below the current upper limit. These numbers are considerably uncertain, but if most high-z galaxies are more like Himiko than like the quasar host SDSS J114816.64+525150.3, then the prospects for detecting the gas and dust in them are much poorer than hitherto anticipated.
[105]  [pdf] - 1171992
Indications of a sub-linear and non-universal Kennicutt-Schmidt relationship
Comments: 7 pages, 2 figures, 1 table. Updated to match MNRAS accepted version
Submitted: 2013-06-12, last modified: 2013-10-10
We estimate the parameters of the Kennicutt-Schmidt (KS) relationship, linking the star formation rate (Sigma_SFR) to the molecular gas surface density (Sigma_mol), in the STING sample of nearby disk galaxies using a hierarchical Bayesian method. This method rigorously treats measurement uncertainties, and provides accurate parameter estimates for both individual galaxies and the entire population. Assuming standard conversion factors to estimate Sigma_SFR and Sigma_mol from the observations, we find that the KS parameters vary between galaxies, indicating that no universal relationship holds for all galaxies. The KS slope of the whole population is 0.76, with the 2sigma range extending from 0.58 to 0.94. These results imply that the molecular gas depletion time is not constant, but varies from galaxy to galaxy, and increases with the molecular gas surface density. Therefore, other galactic properties besides just Sigma_mol affect Sigma_SFR, such as the gas fraction or stellar mass. The non-universality of the KS relationship indicates that a comprehensive theory of star formation must take into account additional physical processes that may vary from galaxy to galaxy.
[106]  [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.
[107]  [pdf] - 1158309
The CO-to-H2 Conversion Factor and Dust-to-Gas Ratio on Kiloparsec Scales in Nearby Galaxies
Comments: accepted to ApJ July 9, 2013, 63 pages, 2-column ApJ format, p38-63 are an image atlas
Submitted: 2012-12-05, last modified: 2013-08-18
We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby galaxies.
[108]  [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.
[109]  [pdf] - 1172957
Carbon Monoxide in the Cold Debris of Supernova 1987A
Comments: Accepted to the Astrophysical Journal Letters, 6 pages, 3 figures
Submitted: 2013-07-24
We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.
[110]  [pdf] - 1159071
The CO-to-H2 Conversion Factor
Comments: To appear in ARAA, Annual Reviews of Astronomy and Astrophysics, ARAA issue 51. This instance is compiled in ApJ format, 35 pages, 12 figures (most in color)
Submitted: 2013-01-15, last modified: 2013-06-08
CO line emission represents the most accessible and widely used tracer of the molecular interstellar medium. This renders the translation of observed CO intensity into total H2 gas mass critical to understand star formation and the interstellar medium in our Galaxy and beyond. We review the theoretical underpinning, techniques, and results of efforts to estimate this CO-to-H2 "conversion factor," Xco, in different environments. In the Milky Way disk, we recommend a conversion factor Xco = 2x10^{20} cm^-2/(K km/s)^-1 with +/-30% uncertainty. Studies of other "normal galaxies" return similar values in Milky Way-like disks, but with greater scatter and systematic uncertainty. Departures from this Galactic conversion factor are both observed and expected. Dust-based determinations, theoretical arguments, and scaling relations all suggest that Xco increases with decreasing metallicity, turning up sharply below metallicity ~1/3-1/2 solar in a manner consistent with model predictions that identify shielding as a key parameter. Based on spectral line modeling and dust observations, Xco appears to drop in the central, bright regions of some but not all galaxies, often coincident with regions of bright CO emission and high stellar surface density. This lower Xco is also present in the overwhelmingly molecular interstellar medium of starburst galaxies, where several lines of evidence point to a lower CO-to-H2 conversion factor. At high redshift, direct evidence regarding the conversion factor remains scarce; we review what is known based on dynamical modeling and other arguments.
[111]  [pdf] - 1171499
Dusty OB stars in the Small Magellanic Cloud - II: Extragalactic Disks or Examples of the Pleiades Phenomenon?
Comments: Accepted for publication in ApJ, 23 pages, 11 figures, 8 tables
Submitted: 2013-05-21
We use mid-infrared Spitzer spectroscopy and far-infrared Herschel photometry for a sample of twenty main sequence O9--B2 stars in the Small Magellanic Cloud (SMC) with strong 24 micron excesses to investigate the origin of the mid-IR emission. Either debris disks around the stars or illuminated patches of dense interstellar medium (ISM) can cause such mid-IR emission. In a companion paper, Paper I, we use optical spectroscopy to show that it is unlikely for any of these sources to be classical Be stars or Herbig Ae/Be stars. We focus our analysis on debris disks and cirrus hot spots. We find three out of twenty stars to be significantly extended in the mid-IR, establishing them as cirrus hot spots. We then fit the IR spectral energy distributions to determine dust temperatures and masses. We find the dust masses in the SMC stars to be larger than for any known debris disks, although this evidence against the debris disk hypothesis is circumstantial. Finally, we created a local comparison sample of bright mid-IR OB stars in the Milky Way (MW) by cross-matching the WISE and Hipparcos catalogs. All such local stars in the appropriate luminosity range that can be unambiguously classified are young stars with optical emission lines or are spatially resolved by WISE with sizes too large to be plausible debris disk candidates. We conclude that the very strong mid-IR flux excesses are most likely explained as cirrus hot spots, although we cannot rigorously rule out that a small fraction of the sample is made up of debris disks or transition disks. We present suggestive evidence that bow-shock heating around runaway stars may be a contributing mechanism to the interstellar emission. These sources, interpreted as cirrus hot spots, offer a new localised probe of diffuse interstellar dust in a low metallicity environment. (Abridged)
[112]  [pdf] - 1171498
Dusty OB Stars in the Small Magellanic Cloud - I: Optical Spectroscopy Reveals Predominantly Main-Sequence OB Stars
Comments: Accepted for publication in ApJ; 38 pages, 18 figures, 7 tables
Submitted: 2013-05-21
We present the results of optical spectroscopic follow-up of 125 candidate main sequence OB stars in the Small Magellanic Cloud (SMC) that were originally identified in the S3MC infrared imaging survey as showing an excess of emission at 24 um indicative of warm dust, such as that associated with a transitional or debris disks. We use these long-slit spectra to investigate the origin of the 24 um emission and the nature of these stars. A possible explanation for the observed 24 um excess, that these are emission line stars with dusty excretion disks, is disproven for the majority of our sources. We find that 88 of these objects are normal stars without line emission, with spectral types mostly ranging from late-O to early-B; luminosity classes a sub-set of our sample indicate that most are main-sequence stars. We further identify 17 emission-line stars, 7 possible emission-line stars, and 5 other stars with forbidden-line emission. We discover a new O6 Iaf star; it exhibits strong HeII 4686 A emission but relatively weak NIII 4640 A emission that we attribute to the lower nitrogen abundance in the SMC. Two other objects are identified with planetary nebulae, one with a young stellar object, and two with X-ray binaries. To shed additional light on the nature of the observed 24 um excess we use optical and infrared photometry to estimate the dust properties of the objects with normal O and B star spectra and compare these properties to those of a sample of hot spots in the Galactic interstellar medium (ISM). We find that the dust properties of the dusty OB star sample resemble the properties of the Galactic sample of hot spots. Some may be runaway systems with bow-shocks resulting from a large velocity difference between star and ISM. We further investigate the nature of these dusty OB stars in a companion paper presenting mid-infrared spectroscopy and additional imaging.
[113]  [pdf] - 657717
Misalignment of Magnetic Fields and Outflows in Protostellar Cores
Comments: 6 pages, 2 figures, 1 table
Submitted: 2012-12-03, last modified: 2013-04-30
We present results of 1.3 mm dust polarization observations toward 16 nearby, low-mass protostars, mapped with ~2.5" resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of ~1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.
[114]  [pdf] - 780622
Towards a resolved Kennicutt-Schmidt law at high redshift
Comments: 6 pages, 4 figures, 2 tables, accepted by Astronomy and Astrophysics
Submitted: 2013-01-03, last modified: 2013-04-29
Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.
[115]  [pdf] - 1165761
Clumping and the Interpretation of kpc-Scale Maps of the Interstellar Medium: Smooth HI and Clumpy, Variable H2 Surface Density
Comments: 6 pages, 2 figures, accepted for publication in ApJ Letters
Submitted: 2013-04-04
Many recent models consider the structure of individual interstellar medium (ISM) clouds as a way to explain observations of large parts of galaxies. To compare such models to observations, one must understand how to translate between surface densities observed averaging over large (~kpc) scales and surface densities on the scale of individual clouds (~pc scale), which are treated by models. We define a "clumping factor" that captures this translation as the ratio of the mass-weighted surface density, which is often the quantity of physical interest, to the area-weighted surface density, which is observed. We use high spatial resolution (sub-kpc) maps of CO and HI emission from nearby galaxies to measure the clumping factor of both atomic and molecular gas. The molecular and atomic ISM exhibit dramatically different degrees of clumping. As a result, the ratio H2/HI measured at ~kpc resolution cannot be trivially interpreted as a cloud-scale ratio of surface densities. HI emission appears very smooth, with a clumping factor of only ~1.3. Based on the scarce and heterogeneous high resolution data available, CO emission is far more clumped with a widely variable clumping factor, median ~7 for our heterogeneous data. Our measurements do not provide evidence for a universal mass-weighted surface density of molecular gas, but also cannot conclusively rule out such a scenario. We suggest that a more sophisticated treatment of molecular ISM structure, one informed by high spatial resolution CO maps, is needed to link cloud-scale models to kpc-scale observations of galaxies.
[116]  [pdf] - 1165679
PHIBSS: molecular gas, extinction, star formation and kinematics in the z=1.5 star forming galaxy EGS13011166
Comments: submitted to the Astrophysical Journal
Submitted: 2013-04-02
We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.
[117]  [pdf] - 1165654
Millimeter dust emission compared with other mass estimates in N11 molecular clouds in the LMC
Comments: A&A accepted, 13 pages, 10 figures
Submitted: 2013-04-01
CO and dust emission at millimeter wavelengths are independent tracers of cold interstellar matter, which have seldom been compared on the scale of GMCs in other galaxies. In this study, and for the first time in the Large Magellanic Cloud, we compute the molecular cloud masses from the mm emission of the dust and compare them with the masses derived from their CO luminosity and virial theorem. We present CO (J=1-0,2-1) and 1.2 mm continuum observations of the N11 star forming region in the LMC obtained with the SEST telescope and the SIMBA bolometer, respectively. We use the CO data to identify individual molecular clouds and measure their physical properties. The correlations between the properties of the N11 clouds are in agreement with those found in earlier studies in the LMC that sample a larger set of clouds and a larger range of cloud masses. For the N11 molecular clouds, we compare the masses estimated from the CO luminosity (Xco\Lco), the virial theorem (Mvir) and the millimeter dust luminosity (L_d). The measured ratios Lco/Mvir and L_d/Mvir constrain the Xco and K_d (dust emissivity at 1.2 mm per unit gas mass) parameters as a function of the virial parameter a_vir. The comparison between the different mass estimates yields a Xco-factor of 8.8x10^20 cm-2/(K km s-1) x a_vir and a K_d parameter of 1.5x10^-3 cm2/g x a_vir. We compare our N11 results with a similar analysis for molecular clouds in the Gould's Belt. We do not find in N11 a large discrepancy between the dust mm and virial masses as reported in earlier studies of molecular clouds in the SMC. The ratio between L_d and Mvir in N11 is half of that measured for Gould's Belt clouds, which can be accounted for by a factor of two lower gas-to-dust mass ratio, as the difference in gas metallicities. If the two samples have similar a_vir values, this result implies that their dust far-IR properties are also similar.
[118]  [pdf] - 1165625
The EGNoG Survey: Molecular Gas in Intermediate-Redshift Star-Forming Galaxies
Comments: Accepted for publication in the Astrophysical Journal; 29 pages, 20 figures, 6 tables
Submitted: 2013-03-29
We present the Evolution of molecular Gas in Normal Galaxies (EGNoG) survey, an observational study of molecular gas in 31 star-forming galaxies from z=0.05 to z=0.5, with stellar masses of (4-30)x10^10 M_Sun and star formation rates of 4-100 M_Sun yr^-1. This survey probes a relatively un-observed redshift range in which the molecular gas content of galaxies is expected to have evolved significantly. To trace the molecular gas in the EGNoG galaxies, we observe the CO(1-0) and CO(3-2) rotational lines using the Combined Array for Research in Millimeter-wave Astronomy (CARMA). We detect 24 of 31 galaxies and present resolved maps of 10 galaxies in the lower redshift portion of the survey. We use a bimodal prescription for the CO to molecular gas conversion factor, based on specific star formation rate, and compare the EGNoG galaxies to a large sample of galaxies assembled from the literature. We find an average molecular gas depletion time of 0.76 \pm 0.54 Gyr for normal galaxies and 0.06 \pm 0.04 Gyr for starburst galaxies. We calculate an average molecular gas fraction of 7-20% at the intermediate redshifts probed by the EGNoG survey. By expressing the molecular gas fraction in terms of the specific star formation rate and molecular gas depletion time (using typical values), we also calculate the expected evolution of the molecular gas fraction with redshift. The predicted behavior agrees well with the significant evolution observed from z~2.5 to today.
[119]  [pdf] - 1157432
The Ionized Circumstellar Envelopes of Orion Source I and the Becklin Neugebauer Object
Comments: 8 pages, 5 figures; accepted for publication in ApJ
Submitted: 2012-10-31, last modified: 2013-01-27
The 229 GHz (lambda 1.3mm) radio emission from Orion-KL was mapped with up to 0.14" angular resolution with CARMA, allowing measurements of the flux densities of Source I ('SrcI') and the Becklin-Neugebauer Object (BN), the 2 most massive stars in this region. We find integrated flux densities of 310 +/- 45 mJy for SrcI and 240 +/- 35 mJy for BN. SrcI is optically thick even at 229 GHz. No trace of the H30alpha recombination line is seen in its spectrum, although the v_2=1, 5(5,0)-6(4,3) transition of H2O, 3450 K above the ground state, is prominent. SrcI is elongated at position angle 140 degrees, as in 43 GHz images. These results are most easily reconciled with models in which the radio emission from SrcI arises via the H- free-free opacity in a T < 4500 K disk, as considered by Reid et al. (2007). By contrast, the radio spectrum of BN is consistent with p+/e- free-free emission from a dense (n_e ~ 5x10^7 cm^{-3}), but otherwise conventional, hypercompact HII region. The source is becoming optically thin at 229 GHz, and the H30alpha recombination line, at VLSR = 23.2 +/- 0.5 km/sec, is prominent in its spectrum. A Lyman continuum flux of 5x10^{45} photons/sec, consistent with that expected from a B star, is required to maintain the ionization. Supplementary 90 GHz observations were made to measure the H41alpha and H42alpha recombination lines toward BN. Published 43 and 86 GHz data suggest that SrcI brightened with respect to BN over the 15 year period from 1994 to 2009.
[120]  [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.
[121]  [pdf] - 1158814
The EGNoG Survey: Gas Excitation in Normal Galaxies at z~0.3
Comments: Accepted for publication in the Astrophysical Journal, to appear January 2013; 18 pages, 10 figures, 6 tables
Submitted: 2013-01-03
As observations of molecular gas in galaxies are pushed to lower star formation rate galaxies at higher redshifts, it is becoming increasingly important to understand the conditions of the gas in these systems to properly infer their molecular gas content. The rotational transitions of the carbon monoxide (CO) molecule provide an excellent probe of the gas excitation conditions in these galaxies. In this paper we present the results from the gas excitation sample of the Evolution of molecular Gas in Normal Galaxies (EGNoG) survey at the Combined Array for Research in Millimeter-wave Astronomy (CARMA). This subset of the full EGNoG sample consists of four galaxies at z~0.3 with star formation rates of 40-65 M_Sun yr^-1 and stellar masses of ~2x10^11 M_Sun. Using the 3 mm and 1 mm bands at CARMA, we observe both the CO(1-0) and CO(3-2) transitions in these four galaxies in order to probe the excitation of the molecular gas. We report robust detections of both lines in three galaxies (and an upper limit on the fourth), with an average line ratio, r_31 = L'_CO(3-2) / L'_CO(1-0), of 0.46 \pm 0.07 (with systematic errors \lesssim 40%), which implies sub-thermal excitation of the CO(3-2) line. We conclude that the excitation of the gas in these massive, highly star-forming galaxies is consistent with normal star-forming galaxies such as local spirals, not starbursting systems like local ultra-luminous infrared galaxies. Since the EGNoG gas excitation sample galaxies are selected from the main sequence of star-forming galaxies, we suggest that this result is applicable to studies of main sequence galaxies at intermediate and high redshifts, supporting the assumptions made in studies that find molecular gas fractions in star forming galaxies at z~1-2 to be an order of magnitude larger than what is observed locally.
[122]  [pdf] - 1158564
The VIRUS-P Exploration of Nearby Galaxies (VENGA): The Xco Gradient in NGC 628
Comments: 10 pages, 7 figures, accepted for publication in The Astrophysical Journal
Submitted: 2012-12-17
We measure the radial profile of the 12CO(1-0) to H_2 conversion factor (Xco) in NGC 628. The H\alpha emission from the VENGA integral field spectroscopy is used to map the star formation rate surface density (\Sigma_{SFR}). We estimate the molecular gas surface density (\Sigma_{H2}) from \Sigma_{SFR} by inverting the molecular star formation law (SFL), and compare it to the CO intensity to measure Xco. We study the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers (H\alpha vs. far-UV plus 24\mu m), and CO maps from different telescopes (single-dish and interferometers). The observed Xco profile is robust against these systematics, drops by a factor of 2 from R~7 kpc to the center of the galaxy, and is well fit by a gradient \Delta log(Xco)=0.06\pm0.02 dex kpc^-1. We study how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical models show that the gradient in Xco can be explained by a combination of decreasing metallicity, and decreasing \Sigma_{H2} with radius. Photoelectric heating from the local UV radiation field appears to contribute to the decrease of Xco in higher density regions. Our results show that galactic environment plays an important role at setting the physical conditions in star forming regions, in particular the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We caution against adopting a single Xco value when large changes in gas surface density or metallicity are present.
[123]  [pdf] - 1158058
The Molecular Gas Density in Galaxy Centers and How It Connects to Bulges
Comments: Accepted for publication in ApJ
Submitted: 2012-11-26
In this paper we present gas density, star formation rate, stellar masses, and bulge disk decompositions for a sample of 60 galaxies. Our sample is the combined sample of BIMA SONG, CARMA STING, and PdBI NUGA surveys. We study the effect of using CO-to-H_2 conversion factors that depend on the CO surface brightness, and also that of correcting star formation rates for diffuse emission from old stellar populations. We estimate that star formation rates in bulges are typically lower by 20% when correcting for diffuse emission. We find that over half of the galaxies in our sample have molecular gas surface density >100 M_sun pc^-2. We find a trend between gas density of bulges and bulge Sersic index; bulges with lower Sersic index have higher gas density. Those bulges with low Sersic index (pseudobulges) have gas fractions that are similar to that of disks. We also find that there is a strong correlation between bulges with the highest gas surface density and the galaxy being barred. However, we also find that classical bulges with low gas surface density can be barred as well. Our results suggest that understanding the connection between the central surface density of gas in disk galaxies and the presence of bars should also take into account the total gas content of the galaxy and/or bulge Sersic index. Indeed, we find that high bulge Sersic index is the best predictor of low gas density inside the bulge (not barredness of the disk). Finally, we show that when using the corrected star formation rates and gas densities, the correlation between star formation rate surface density and gas surface density of bulges is similar to that of disks.
[124]  [pdf] - 1158017
PHIBSS: molecular gas content and scaling relations in z~1-3 normal star forming galaxies
Comments: Submitted to the Astrophysical Journal; 14 figures
Submitted: 2012-11-25
We present PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.
[125]  [pdf] - 1124904
Modeling Dust and Starlight in Galaxies Observed by Spitzer and Herschel: NGC 628 and NGC 6946
Comments: To be published in Apj, September 2012. See the full version at
Submitted: 2012-07-17
We characterize the dust in NGC628 and NGC6946, two nearby spiral galaxies in the KINGFISH sample. With data from 3.6um to 500um, dust models are strongly constrained. Using the Draine & Li (2007) dust model, (amorphous silicate and carbonaceous grains), for each pixel in each galaxy we estimate (1) dust mass surface density, (2) dust mass fraction contributed by polycyclic aromatic hydrocarbons (PAH)s, (3) distribution of starlight intensities heating the dust, (4) total infrared (IR) luminosity emitted by the dust, and (5) IR luminosity originating in regions with high starlight intensity. We obtain maps for the dust properties, which trace the spiral structure of the galaxies. The dust models successfully reproduce the observed global and resolved spectral energy distributions (SEDs). The overall dust/H mass ratio is estimated to be 0.0082+/-0.0017 for NGC628, and 0.0063+/-0.0009 for NGC6946, consistent with what is expected for galaxies of near-solar metallicity. Our derived dust masses are larger (by up to a factor 3) than estimates based on single-temperature modified blackbody fits. We show that the SED fits are significantly improved if the starlight intensity distribution includes a (single intensity) "delta function" component. We find no evidence for significant masses of cold dust T<12K. Discrepancies between PACS and MIPS photometry in both low and high surface brightness areas result in large uncertainties when the modeling is done at PACS resolutions, in which case SPIRE, MIPS70 and MIPS160 data cannot be used. We recommend against attempting to model dust at the angular resolution of PACS.
[126]  [pdf] - 1118174
Dust-to-Gas Ratio in the Extremely Metal Poor Galaxy I ZW 18
Comments: Accepted for publication in ApJ: 8 pages, 4 figures, 1 table
Submitted: 2012-04-20
The blue compact dwarf galaxy I Zw 18 is one of the most metal poor systems known in the local Universe (12 + log(O/H) $=$ 7.17). In this work we study I Zw 18 using data from {\it Spitzer}, {\it Herschel Space Telescope} and IRAM Plateau de Bure Interferometer. Our data set includes the most sensitive maps of I Zw 18, to date, in both, the far infrared and the CO $J=1\rightarrow0$ transition. We use dust emission models to derive a dust mass upper limit of only M$_{dust}\leq1.1\times10^4$ M$_{\odot}$ ($3\sigma$ limit). This upper limit is driven by the non-detection at 160 $\mu$m, and it is a factor of 4-10 times smaller than previous estimates (depending upon the model used). We also estimate an upper limit to the total dust-to-gas mass ratio of M$_{Dust}$/M$_{gas}\leq5.0\times10^{-5}$. If a linear correlation between the dust-to-gas mass ratio and metallicity (measure as O/H) were to hold, we would expect a ratio of 3.9$\times10^{-4}$. We also show that the infrared SED is similar to that of starbursting systems.
[127]  [pdf] - 1117726
A Study of Heating and Cooling of the ISM in NGC 1097 with Herschel-PACS and Spitzer-IRS
Comments: 17 pages, 14 figures, 5 tables; accepted for publication in ApJ
Submitted: 2012-04-02
NGC 1097 is a nearby Seyfert 1 galaxy with a bright circumnuclear starburst ring, a strong large-scale bar and an active nucleus. We present a detailed study of the spatial variation of the far infrared (FIR) [CII]158um and [OI]63um lines and mid-infrared H2 emission lines as tracers of gas cooling, and of the polycyclic aromatic hydrocarbon (PAH) bands as tracers of the photoelectric heating, using Herschel-PACS, and Spitzer-IRS infrared spectral maps. We focus on the nucleus and the ring, and two star forming regions (Enuc N and Enuc S). We estimated a photoelectric gas heating efficiency ([CII]158um+[OI]63um)/PAH in the ring about 50% lower than in Enuc N and S. The average 11.3/7.7um PAH ratio is also lower in the ring, which may suggest a larger fraction of ionized PAHs, but no clear correlation with [CII]158{\mu}m/PAH(5.5 - 14um) is found. PAHs in the ring are responsible for a factor of two more [CII]158um and [OI]63um emission per unit mass than PAHs in the Enuc S. SED modeling indicates that at most 25% of the FIR power in the ring and Enuc S can come from high intensity photodissociation regions (PDRs), in which case G0 ~ 10^2.3 and nH ~ 10^3.5 cm^-3 in the ring. For these values of G0 and nH PDR models cannot reproduce the observed H2 emission. Much of the the H2 emission in the starburst ring could come from warm regions in the diffuse ISM that are heated by turbulent dissipation or shocks.
[128]  [pdf] - 1116609
Estimating the Star Formation Rate at 1 kpc Scales in Nearby Galaxies
Comments: 34 pages, 24 figures, accepted for publication in the Astronomical Journal
Submitted: 2012-02-13
Using combinations of H\alpha, ultraviolet (UV), and infrared (IR) emission, we estimate the star formation rate (SFR) surface density, \Sigma_SFR, at 1 kpc resolution for 30 disk galaxies that are targets of the IRAM HERACLES CO survey. We present a new physically-motivated IR spectral energy distribution-based approach to account for possible contributions to 24\mum emission not associated with recent star formation. Considering a variety of "reference" SFRs from the literature, we revisit the calibration of the 24\mum term in hybrid (UV+IR or H\alpha+IR) tracers. We show that the overall calibration of this term remains uncertain at the factor of two level because of the lack of wide-field, robust reference SFR estimates. Within this uncertainty, published calibrations represent a reasonable starting point for 1 kpc-wide areas of star-forming disk galaxies but we re-derive and refine the calibration of the IR term in these tracers to match our resolution and approach to 24\mum emission. We compare a large suite of \Sigma_SFR estimates and find that above \Sigma_SFR \sim 10^-3 M_\odot yr^-1 kpc^-2 the systematic differences among tracers are less than a factor of two across two orders of magnitude dynamic range. We caution that methodology and data both become serious issues below this level. We note from simple model considerations that focusing on a part of a galaxy dominated by a single stellar population the intrinsic uncertainty in H\alpha and FUV-based SFRs are \sim 0.3 and \sim 0.5 dex.
[129]  [pdf] - 1092738
Resolving the Far-IR Line Deficit: Photoelectric Heating and Far-IR Line Cooling in NGC 1097 and NGC 4559
Comments: 19 pages, 14 pages, accepted for publication in ApJ
Submitted: 2012-01-04
The physical state of interstellar gas and dust is dependent on the processes which heat and cool this medium. To probe heating and cooling of the ISM over a large range of infrared surface brightness, on sub-kiloparsec scales, we employ line maps of [C \ii] 158 $\mu$m, [O \one] 63 $\mu$m, and [N \ii] 122 $\mu$m in NGC 1097 and NGC 4559, obtained with the PACS spectrometer onboard {\it Herschel}. We matched new observations to existing Spitzer-IRS data that trace the total emission of polycyclic aromatic hydrocarbons (PAHs). We confirm at small scales in these galaxies that the canonical measure of photoelectric heating efficiency, ([C \ii] + [O \one])/TIR, decreases as the far-infrared color, $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m), increases. In contrast, the ratio of far-infrared (far-IR) cooling to total PAH emission, ([C \ii] + [O \one])/PAH, is a near constant $\sim$6% over a wide range of far-infrared color, 0.5 \textless\ $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\lesssim$ 0.95. In the warmest regions, where $\nu f_\nu$(70 $\mu$m)/$\nu f_\nu$(100 $\mu$m) $\gtrsim$ 0.95, the ratio ([C \ii] + [O \one])/PAH drops rapidly to 4%. We derived representative values of the local UV radiation density, $G_0$, and the gas density, $n_H$, by comparing our observations to models of photodissociation regions. The ratio $G_0/n_H$, derived from fine-structure lines, is found to correlate with the mean dust-weighted starlight intensity, $<U>$ derived from models of the IR SED. Emission from regions that exhibit a line deficit is characterized by an intense radiation field, indicating that small grains are susceptible to ionization effects. We note that there is a shift in the 7.7 / 11.3 $\mu$m PAH ratio in regions that exhibit a deficit in ([C \ii] + [O \one])/PAH, suggesting that small grains are ionized in these environments.
[130]  [pdf] - 1092174
Herschel Far-Infrared and Sub-millimeter Photometry for the KINGFISH Sample of Nearby Galaxies
Comments: Accepted for publication in ApJ
Submitted: 2011-12-05
New far-infrared and sub-millimeter photometry from the Herschel Space Observatory is presented for 61 nearby galaxies from the Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel (KINGFISH) sample. The spatially-integrated fluxes are largely consistent with expectations based on Spitzer far-infrared photometry and extrapolations to longer wavelengths using popular dust emission models. Dwarf irregular galaxies are notable exceptions, as already noted by other authors, as their 500um emission shows evidence for a sub-millimeter excess. In addition, the fraction of dust heating attributed to intense radiation fields associated with photo-dissociation regions is found to be (21+/-4)% larger when Herschel data are included in the analysis. Dust masses obtained from the dust emission models of Draine & Li are found to be on average nearly a factor of two higher than those based on single-temperature modified blackbodies, as single blackbody curves do not capture the full range of dust temperatures inherent to any galaxy. The discrepancy is largest for galaxies exhibiting the coolest far-infrared colors.
[131]  [pdf] - 1091782
KINGFISH -- Key Insights on Nearby Galaxies: A Far-Infrared Survey with Herschel: Survey Description and Image Atlas
Comments: 52 pages, 10 figures, 2 tables. Accepted for publication in PASP. A pdf copy of the image atlas can be downloaded separately at the KINGFISH public website,
Submitted: 2011-11-18
The KINGFISH project (Key Insights on Nearby Galaxies: a Far-Infrared Survey with Herschel) is an imaging and spectroscopic survey of 61 nearby (d < 30 Mpc) galaxies, chosen to cover a wide range of galaxy properties and local interstellar medium (ISM) environments found in the nearby Universe. Its broad goals are to characterize the ISM of present-day galaxies, the heating and cooling of their gaseous and dust components, and to better understand the physical processes linking star formation and the ISM. KINGFISH is a direct descendant of the Spitzer Infrared Nearby Galaxies Survey (SINGS), which produced complete Spitzer imaging and spectroscopic mapping and a comprehensive set of multi-wavelength ancillary observations for the sample. The Herschel imaging consists of complete maps for the galaxies at 70, 100, 160, 250, 350, and 500 microns. The spectal line imaging of the principal atomic ISM cooling lines ([OI]63um, [OIII]88um, [NII]122,205um, and [CII]158um) covers the subregions in the centers and disks that already have been mapped in the mid-infrared with Spitzer. The KINGFISH and SINGS multi-wavelength datasets combined provide panchromatic mapping of the galaxies sufficient to resolve individual star-forming regions, and tracing the important heating and cooling channels of the ISM, across a wide range of local extragalactic ISM environments. This paper summarizes the scientific strategy for KINGFISH, the properties of the galaxy sample, the observing strategy, and data processing and products. It also presents a combined Spitzer and Herschel image atlas for the KINGFISH galaxies, covering the wavelength range 3.6 -- 500 microns. All imaging and spectroscopy data products will be released to the Herschel user generated product archives.
[132]  [pdf] - 1077204
The metallicity dependence of the CO {\rightarrow} H_2 conversion factor in z>1 star forming galaxies
Comments: ApJ in press
Submitted: 2011-06-10, last modified: 2011-11-11
We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor {\alpha}CO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log({\alpha}CO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.
[133]  [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.
[134]  [pdf] - 449148
The Spitzer Spectroscopic Survey of the Small Magellanic Cloud (S4MC): Probing the Physical State of Polycyclic Aromatic Hydrocarbons in a Low-Metallicity Environment
Comments: Submitted to ApJ
Submitted: 2011-09-05
We present results of mid-infrared spectroscopic mapping observations of six star-forming regions in the Small Magellanic Cloud from the Spitzer Spectroscopic Survey of the SMC (S4MC). We detect the mid-IR emission from polycyclic aromatic hydrocarbons (PAHs) in all of the mapped regions, greatly increasing the range of environments where PAHs have been spectroscopically detected in the SMC. We investigate the variations of the mid-IR bands in each region and compare our results to studies of the PAH bands in the SINGS sample and in a sample of low-metallicity starburst galaxies. PAH emission in the SMC is characterized by low ratios of the 6-9 micron features relative to the 11.3 micron feature and weak 8.6 and 17.0 micron features. Interpreting these band ratios in the light of laboratory and theoretical studies, we find that PAHs in the SMC tend to be smaller and less ionized than those in higher metallicity galaxies. Based on studies of PAH destruction, we argue that a size distribution shifted towards smaller PAHs cannot be the result of processing in the interstellar medium, but instead reflects differences in the formation of PAHs at low metallicity. Finally, we discuss the implications of our observations for our understanding of the PAH life-cycle in low-metallicity galaxies---namely that the observed deficit of PAHs may be a consequence of PAHs forming with smaller average sizes and therefore being more susceptible to destruction under typical interstellar medium conditions.
[135]  [pdf] - 1077423
The Emission by Dust and Stars of Nearby Galaxies in the Herschel KINGFISH Survey
Comments: 19 pages, 12 figures; published in ApJ
Submitted: 2011-06-20, last modified: 2011-08-16
Using new far-infrared imaging from the Herschel Space Observatory with ancillary data from ultraviolet to submillimeter wavelengths, we estimate the total emission from dust and stars of 62 nearby galaxies in the KINGFISH survey in a way that is as empirical and model-independent as possible. We collect and exploit these data in order to measure from the spectral energy distributions (SEDs) precisely how much stellar radiation is intercepted and re-radiated by dust, and how this quantity varies with galaxy properties. By including SPIRE data, we are more sensitive to emission from cold dust grains than previous analyses at shorter wavelengths, allowing for more accurate estimates of dust temperatures and masses. The dust/stellar flux ratio, which we measure by integrating the SEDs, has a range of nearly three decades. The inclusion of SPIRE data shows that estimates based on data not reaching these far-IR wavelengths are biased low. We find that the dust/stellar flux ratio varies with morphology and total IR luminosity. We also find that dust/stellar flux ratios are related to gas-phase metallicity, while the dust/stellar mass ratios are less so. The substantial scatter between dust/stellar flux and dust/stellar mass indicates that the former is a poor proxy of the latter. Comparing the dust/stellar flux ratios and dust temperatures, we show that early-types tend to have slightly warmer temperatures than spiral galaxies, which may be due to more intense interstellar radiation fields, or to different dust grain compositions. Finally, we show that early-types and early-type spirals have a strong correlation between the dust/stellar flux ratio and specific star formation rate, which suggests that the relatively bright far-IR emission of some of these galaxies is due to ongoing star formation and the radiation field from older stars.
[136]  [pdf] - 1077188
Disk and Envelope Structure in Class 0 Protostars: II. High Resolution Millimeter Mapping of the Serpens Sample
Comments: Accepted to the Astrophysical Journal Supplements
Submitted: 2011-06-09, last modified: 2011-08-01
We present high-resolution CARMA 230 GHz continuum imaging of nine deeply embedded protostars in the Serpens Molecular Cloud, including six of the nine known Class 0 protostars in Serpens. This work is part of a program to characterize disk and envelope properties for a complete sample of Class 0 protostars in nearby low-mass star forming regions. Here we present CARMA maps and visibility amplitudes as a function of uv-distance for the Serpens sample. Observations are made in the B, C, D, and E antenna configurations, with B configuration observations utilizing the CARMA Paired Antenna Calibration System. Combining data from multiple configurations provides excellent uv-coverage (4-500 klam), allowing us to trace spatial scales from 1e2 to 1e4 AU. We find evidence for compact disk components in all of the observed Class 0 protostars, suggesting that disks form at very early times (t<0.2 Myr) in Serpens. We make a first estimate of disk masses using the flux at 50 klam, where the contribution from the envelope should be negligible, assuming an unresolved disk. The resulting disk masses range from 0.04 Msun to 1.7 Msun, with a mean of approximately 0.2 Msun. Our high resolution maps are also sensitive to binary or multiple sources with separations > 250 AU, but significant evidence of multiplicity on scales <2000 AU is seen in only one source.
[137]  [pdf] - 389420
Physical conditions of the interstellar medium of high-redshift, strongly lensed submillimetre galaxies from the Herschel-ATLAS
Comments: 13 pages, 6 figures, published in MNRAS
Submitted: 2011-05-19, last modified: 2011-07-22
We present Herschel-SPIRE Fourier Transform Spectrometer (FTS) and radio follow-up observations of two Herschel-ATLAS (H-ATLAS) detected strongly lensed distant galaxies. In one of the targeted galaxies H-ATLAS J090311.6+003906 (SDP.81) we detect [OIII] 88\mum and [CII] 158\mum lines at a signal-to-noise ratio of ~5. We do not have any positive line identification in the other fainter target H-ATLAS J091305.0-005343 (SDP.130). Currently SDP.81 is the faintest sub-mm galaxy with positive line detections with the FTS, with continuum flux just below 200 mJy in the 200-600 \mum wavelength range. The derived redshift of SDP.81 from the two detections is z=3.043 +/-0.012, in agreement with ground-based CO measurements. This is the first detection by Herschel of the [OIII] 88\mum line in a galaxy at redshift higher than 0.05. Comparing the observed lines and line ratios with a grid of photo-dissociation region (PDR) models with different physical conditions, we derive the PDR cloud density n ~ 2000 cm^{-3} and the far-UV ionizing radiation field G_0 ~ 200 (in units of the Habing field -- the local Galactic interstellar radiation field of 1.6x10^{-6} W/m^2). Using the CO derived molecular mass and the PDR properties we estimate the effective radius of the emitting region to be 500-700 pc. These characteristics are typical for star-forming, high redshift galaxies. The radio observations indicate that SDP.81 deviates significantly from the local FIR/radio correlation, which hints that some fraction of the radio emission is coming from an AGN. The constraints on the source size from millimiter-wave observations put a very conservative upper limit of the possible AGN contribution to less than 33%. These indications, together with the high [OIII]/FIR ratio and the upper limit of [OI] 63\mum/[CII] 158\mum suggest that some fraction of the ionizing radiation is likely to originate from an AGN.
[138]  [pdf] - 1078137
Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud (SAGE-SMC). I. Overview
Comments: 16 pages, 11 figures, AJ, in press
Submitted: 2011-07-21
The Small Magellanic Cloud (SMC) provides a unique laboratory for the study of the lifecycle of dust given its low metallicity (~1/5 solar) and relative proximity (~60 kpc). This motivated the SAGE-SMC (Surveying the Agents of Galaxy Evolution in the Tidally-Stripped, Low Metallicity Small Magellanic Cloud) Spitzer Legacy program with the specific goals of studying the amount and type of dust in the present interstellar medium, the sources of dust in the winds of evolved stars, and how much dust is consumed in star formation. This program mapped the full SMC (30 sq. deg.) including the Body, Wing, and Tail in 7 bands from 3.6 to 160 micron using the IRAC and MIPS instruments on the Spitzer Space Telescope. The data were reduced, mosaicked, and the point sources measured using customized routines specific for large surveys. We have made the resulting mosaics and point source catalogs available to the community. The infrared colors of the SMC are compared to those of other nearby galaxies and the 8 micron/24 micron ratio is somewhat lower and the 70 micron/160 micron ratio is somewhat higher than the average. The global infrared spectral energy distribution shows that the SMC has ~3X lower aromatic emission/PAH (polycyclic aromatic hydrocarbon) abundances compared to most nearby galaxies. Infrared color-magnitude diagrams are given illustrating the distribution of different asymptotic giant branch stars and the locations of young stellar objects. Finally, the average spectral energy distribution (SED) of HII/star formation regions is compared to the equivalent Large Magellanic Cloud average HII/star formation region SED. These preliminary results are expanded in detail in companion papers.
[139]  [pdf] - 1077453
The Water Vapor Spectrum of APM 08279+5255: X-Ray Heating and Infrared Pumping over Hundreds of Parsecs
Comments: Astrophysical Journal, Letters, Accepted
Submitted: 2011-06-21, last modified: 2011-07-18
We present the rest-frame 200--320 \mm\ spectrum of the z=3.91 quasar \apm, obtained with Z-Spec at the Caltech Submillimeter Observatory. In addition to the \jeight\ to \jthirteen\ CO rotational transitions which dominate the CO cooling, we find six transitions of water originating at energy levels ranging up to 643 K. Most are first detections at high redshift, and we have confirmed one transition with CARMA. The CO cooling is well-described by our XDR model, assuming L$_{\rm 1-100\,keV}\sim1\times10^{46}\rm\,erg\,s^{-1}$, and that the gas is distributed over a 550-pc sizescale, per the now-favored $\mu$=4 lensing model. The total observed cooling in water corresponds to 6.5$\times10^{9}$ \ls, comparable to that of CO. We compare the water spectrum with that of Mrk 231, finding that the intensity ratios among the high-lying lines are similar, but with a total luminosity scaled up by a factor of $\sim$50. Using this scaling, we estimate an average water abundance relative to \hh\ of 1.4$\times10^{-7}$, a good match to the prediction of the chemical network in the XDR model. As with Mrk 231, the high-lying water transitions are excited radiatively via absorption in the rest-frame far-infrared, and we show that the powerful dust continuum in \apm\ is more than sufficient to pump this massive reservoir of warm water vapor.
[140]  [pdf] - 1078018
Radio Interferometric Planet Search II: Constraints on sub-Jupiter-Mass Companions to GJ 896A
Comments: Accepted for publication in ApJ
Submitted: 2011-07-15
We present results from the Radio Interferometric Planet (RIPL) search for compan- ions to the nearby star GJ 896A. We present 11 observations over 4.9 years. Fitting astrometric parameters to the data reveals a residual with peak-to-peak amplitude of ~ 3 mas in right ascension. This residual is well-fit by an acceleration term of 0.458 \pm 0.032 mas/y^2. The parallax is fit to an accuracy of 0.2 mas and the proper motion terms are fit to accuracies of 0.01 mas/y. After fitting astrometric and acceleration terms residuals are 0.26 mas in each coordinate, demonstrating that stellar jitter does not limit the ability to carry out radio astrometric planet detection and characterization. The acceleration term originates in part from the companion GJ 896B but the amplitude of the acceleration in declination is not accurately predicted by the orbital model. The acceleration sets a mass upper limit of 0.15 MJ at a semi-major axis of 2 AU for a planetary companion to GJ 896A. For semi-major axes between 0.3 and 2 AU upper limits are determined by the maximum angular separation; the upper limits scale from the minimum value in proportion to the inverse of the radius. Upper limits at larger radii are set by the acceleration and scale as the radius squared. An improved solution for the stellar binary system could improve the exoplanet mass sensitivity by an order of magnitude.
[141]  [pdf] - 1077858
The State of the Gas and the Relation Between Gas and Star Formation at Low Metallicity: the Small Magellanic Cloud
Comments: Accepted for publication in the Astrophysical Journal. 21 pages and 11 figures in apjemulate style
Submitted: 2011-07-08
We compare atomic gas, molecular gas, and the recent star formation rate (SFR) inferred from H-alpha in the Small Magellanic Cloud (SMC). By using infrared dust emission and local dust-to-gas ratios, we construct a map of molecular gas that is independent of CO emission. This allows us to disentangle conversion factor effects from the impact of metallicity on the formation and star formation efficiency of molecular gas. On scales of 200 pc to 1 kpc we find a characteristic molecular gas depletion time of ~1.6 Gyr, similar to that observed in the molecule-rich parts of large spiral galaxies on similar spatial scales. This depletion time shortens on much larger scales to ~0.6 Gyr because of the presence of a diffuse H-alpha component, and lengthens on much smaller scales to ~7.5 Gyr because the H-alpha and H2 distributions differ in detail. We estimate the systematic uncertainties in our measurement to be a factor of 2-3. We suggest that the impact of metallicity on the physics of star formation in molecular gas has at most this magnitude. The relation between SFR and neutral (H2+HI) gas surface density is steep, with a power-law index ~2.2+/-0.1, similar to that observed in the outer disks of large spiral galaxies. At a fixed total gas surface density the SMC has a 5-10 times lower molecular gas fraction (and star formation rate) than large spiral galaxies. We explore the ability of the recent models by Krumholz et al. (2009) and Ostriker et al. (2010) to reproduce our observations. We find that to explain our data at all spatial scales requires a low fraction of cold, gravitationally-bound gas in the SMC. We explore a combined model that incorporates both large scale thermal and dynamical equilibrium and cloud-scale photodissociation region structure and find that it reproduces our data well, as well as predicting a fraction of cold atomic gas very similar to that observed in the SMC.
[142]  [pdf] - 1076853
Calibrating Extinction-Free Star Formation Rate Diagnostics with 33GHz Free-Free Emission in NGC6946
Comments: 18 pages, 7 Figures, Accepted for publication in ApJ
Submitted: 2011-05-24
Abridged: Using free-free emission measured in the Ka-band (26-40GHz) for 10 star-forming regions in the nearby galaxy NGC6946, including its starbursting nucleus, we compare a number of SFR diagnostics that are typically considered to be unaffected by interstellar extinction: i.e., non-thermal radio (i.e., 1.4GHz), total infrared (IR; 8-1000um), and warm dust (i.e., 24um) emission, along with the hybrid (obscured + unobscured) indicators of H\alpha+24um and UV+IR. The 33GHz free-free emission is assumed to provide the most accurate measure of the current SFR. Among the extranuclear star-forming regions, the 24um, H\alpha+24um and UV+IR SFR calibrations are in good agreement with the 33GHz free-free SFRs. However, each of the SFR calibrations relying on some form of dust emission overestimate the nuclear SFR by a factor of ~2. This is more likely the result of excess dust heating through an accumulation of non-ionizing stars associated with an extended episode of star formation in the nucleus rather than increased competition for ionizing photons by dust. SFR calibrations using the non-thermal radio continuum yield values which only agree with the free-free SFRs for the nucleus, and underestimate the SFRs from the extranuclear star-forming regions by a factor of ~2. This result likely arises from the CR electrons decaying within the starburst region with negligible escape compared to the young extranuclear star-forming regions. Finally, we find that the SFRs estimated using the total 33GHz emission agree well with the free-free SFRs due to the large thermal fractions present at these frequencies even when local diffuse backgrounds are not removed. Thus, rest-frame 33GHz observations may act as a reliable method to measure the SFRs of galaxies at increasingly high redshift without the need of ancillary radio data to account for the non-thermal emission.
[143]  [pdf] - 1034273
What Drives the Expansion of Giant HII Regions?: A Study of Stellar Feedback in 30 Doradus
Comments: 14 pages, 17 figures; ApJ in press
Submitted: 2010-08-13, last modified: 2011-03-03
Observations show that star formation is an inefficient and slow process. This result can be attributed to the injection of energy and momentum by stars that prevents free-fall collapse of molecular clouds. The mechanism of this stellar feedback is debated theoretically: possible sources of pressure include the classical warm HII gas, the hot gas generated by shock-heating from stellar winds and supernovae, direct radiation of stars, and the dust-processed radiation field trapped inside the HII shell. In this paper, we measure observationally the pressures associated with each component listed above across the giant HII region 30 Doradus in the Large Magellanic Cloud. We exploit high-resolution, multi-wavelengh images (radio, infrared, optical, and X-ray) to map these pressures as a function of position. We find that radiation pressure dominates within 75 pc of the central star cluster, R136, while the HII gas pressure dominates at larger radii. By contrast, the dust-processed radiation pressure and hot gas pressure are generally weak and not dynamically important, although the hot gas pressure may have played a more significant role at early times. Based on the low X-ray gas pressures, we demonstrate that the hot gas is only partially confined and must be leaking out the HII shell. Additionally, we consider the implications of a dominant radiation pressure on the early dynamics of 30 Doradus.
[144]  [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.
[145]  [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.
[146]  [pdf] - 1051578
The Displaced Dusty ISM of NGC 3077: Tidal Stripping in the M 81 Triplet
Comments: ApJ Letters, in print
Submitted: 2011-01-20
We present the detection of extended (~30 kpc^2) dust emission in the tidal \hi\ arm near NGC 3077 (member of the M\,81 triplet) using SPIRE on board Herschel. Dust emission in the tidal arm is typically detected where the \hi\ column densities are >10^21 cm^-2. The SPIRE band ratios show that the dust in the tidal arm is significantly colder (~13 K) than in NGC 3077 itself (~31 K), consistent with the lower radiation field in the tidal arm. The total dust mass in the tidal arm is ~1.8 x 10^6 M_sun (assuming beta=2), i.e. substantially larger than the dust mass associated with NGC 3077 (~2 x 10^5 M_sun). Where dust is detected, the dust-to-gas ratio is 6+/-3 x 10^-3, consistent within the uncertainties with what is found in NGC 3077 and nearby spiral galaxies with Galactic metallicities. The faint HII regions in the tidal arm can not be responsible for the detected enriched material and are not the main source of the dust heating in the tidal arm. We conclude that the interstellar medium (atomic HI, molecules and dust) in this tidal feature was pre-enriched and stripped off NGC 3077 during its recent interaction (~3 x 10^8 yr ago) with M 82 and M 81. This implies that interaction can efficiently remove heavy elements and enriched material (dust, molecular gas) from galaxies. As interactions were more frequent at large lookback times, it is conceivable that they could substantially contribute (along with galactic outflows) to the enrichment of the intergalactic medium.
[147]  [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.
[148]  [pdf] - 1040809
Atmospheric Phase Correction using CARMA-PACS: High Angular Resolution Observations of the FU-Orionis star PP 13S*
Comments: 12 pages, 9 figures, accepted for publication in ApJ
Submitted: 2010-09-23
We present 0.15" resolution observations of the 227 GHz continuum emission from the circumstellar disk around the FU-Orionis star PP 13S*. The data were obtained with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) Paired Antenna Calibration System (C-PACS), which measures and corrects the atmospheric delay fluctuations on the longest baselines of the array in order to improve the sensitivity and angular resolution of the observations. A description of the C-PACS technique and the data reduction procedures are presented. C-PACS was applied to CARMA observations of PP 13S*, which led to a factor of 1.6 increase in the observed peak flux of the source, a 36% reduction in the noise of the image, and a 52% decrease in the measured size of the source major axis. The calibrated complex visibilities were fitted with a theoretical disk model to constrain the disk surface density. The total disk mass from the best fit model corresponds to 0.06 \msun, which is larger than the median mass of a disk around a classical T Tauri star. The disk is optically thick at a wavelength of 1.3 mm for orbital radii less than 48 AU. At larger radii, the inferred surface density of the PP 13S* disk is an order of magnitude lower than that needed to develop a gravitational instability.
[149]  [pdf] - 1034606
LABOCA observations of giant molecular clouds in the south west region of the Small Magellanic Cloud
Comments: 12 pages, 4 figures, Accepted in A&A
Submitted: 2010-09-01
The amount of molecular gas is a key for understanding the future star formation in a galaxy. Because H2 is difficult to observe directly in dense and cold clouds, tracers like CO are used. However, at low metallicities especially, CO only traces the shielded interiors of the clouds. mm dust emission can be used as a tracer to unveil the total dense gas masses. The comparison of masses deduced from the continuum SIMBA 1.2 mm emission and virial masses in a sample of giant molecular clouds (GMCs), in the SW region of the Small Magellanic Cloud (SMC), showed a discrepancy that is in need of an explanation. This study aims at better assessing possible uncertainties on the dust emission observed in the sample of GMCs from the SMC and focuses on the densest parts of the GMCs where CO is detected. New observations were obtained with the LABOCA camera on the APEX telescope. All GMCs previously observed in CO are detected and their emission at 870microns is compared to ancillary data. The different contributions to the sub-mmm emission are estimated, as well as dust properties, in order to deduce molecular cloud masses precisely. The (sub-)mm emission observed in the GMCs in the SW region of the SMC is dominated by dust emission and masses are deduced for the part of each cloud where CO is detected and compared to the virial masses. The mass discrepancy between both methods is confirmed at 870microns with the LABOCA observations: the virial masses are on average 4 times smaller than the masses of dense gas deduced from dust emission, contrary to what is observed for equivalent clouds in our Galaxy. The origin of this mass discrepancy in the SMC remains unkown. The direct interpretation of this effect is that the CO linewidth used to compute virial masses do not measure the full velocity distribution of the gas. Geometrical effects and uncertainties on the dust properties are also discussed.
[150]  [pdf] - 1032937
Herschel Inventory of The Agents of Galaxy Evolution (HERITAGE): the Large Magellanic Cloud dust
Comments: To be published in Astronomy & Astrophysics, Herschel First Results Issue
Submitted: 2010-06-04
The HERschel Inventory of The Agents of Galaxy Evolution (HERITAGE) of the Magellanic Clouds will use dust emission to investigate the life cycle of matter in both the Large and Small Magellanic Clouds (LMC and SMC). Using the Herschel Space Observatory's PACS and SPIRE photometry cameras, we imaged a 2x8 square degree strip through the LMC, at a position angle of ~22.5 degrees as part of the science demonstration phase of the Herschel mission. We present the data in all 5 Herschel bands: PACS 100 and 160 {\mu}m and SPIRE 250, 350 and 500 {\mu}m. We present two dust models that both adequately fit the spectral energy distribution for the entire strip and both reveal that the SPIRE 500 {\mu}m emission is in excess of the models by 6 to 17%. The SPIRE emission follows the distribution of the dust mass, which is derived from the model. The PAH-to-dust mass (f_PAH) image of the strip reveals a possible enhancement in the LMC bar in agreement with previous work. We compare the gas mass distribution derived from the HI 21 cm and CO J=1-0 line emission maps to the dust mass map from the models and derive gas-to-dust mass ratios (GDRs). The dust model, which uses the standard graphite and silicate optical properties for Galactic dust, has a very low GDR = 65(+15,-18) making it an unrealistic dust model for the LMC. Our second dust model, which uses amorphous carbon instead of graphite, has a flatter emissivity index in the submillimeter and results in a GDR = 287(+25,-42) that is more consistent with a GDR inferred from extinction.
[151]  [pdf] - 1032562
Cold dust clumps in dynamically hot gas
Comments: 4 pages, 4 figures, to appear in A&A special issue
Submitted: 2010-05-18
We present clumps of dust emission from Herschel observations of the Large Magellanic Cloud (LMC) and their physical and statistical properties. We catalog cloud features seen in the dust emission from Herschel observations of the LMC, the Magellanic type irregular galaxy closest to the Milky Way, and compare these features with HI catalogs from the ATCA+Parkes HI survey. Using an automated cloud-finding algorithm, we identify clouds and clumps of dust emission and examine the cumulative mass distribution of the detected dust clouds. The mass of cold dust is determined from physical parameters that we derive by performing spectral energy distribution fits to 250, 350, and 500 micronm emission from SPIRE observations using DUSTY and GRASIL radiative transfer calculation with dust grain size distributions for graphite/silicate in low-metallicity extragalactic environments. The dust cloud mass spectrum follows a power law distribution with an exponent of gamma=-1.8 for clumps larger than 400 solar mass and is similar to the HI mass distribution. This is expected from the theory of ISM structure in the vicinity of star formation.
[152]  [pdf] - 1032476
Dust/gas correlations from Herschel Observations
Submitted: 2010-05-14
Spitzer and IRAS observations of the LMC suggest an excess of FIR emission with respect to the gas surface density traced by 12CO and HI 21 cm emission lines. This "FIR excess" is noticeable near molecular clouds in the LMC, and has usually been interpreted as the presence of a self-shielded H2 component not traced by CO molecular clouds' envelopes. Based on Herschel observations, we examine the correlation between gas and dust at higher resolution than previously achieved. We consider three additional causes for the FIR excess: X factor, FIR dust emissivity, and gas-to-dust ratio variations between the diffuse and dense phases of the ISM. We examine the structure of NT80 and NT71, two molecular clouds detected in the NANTEN 12CO survey of the LMC. Dust surface density maps are derived from the HERITAGE data. The gas phase is traced by MAGMA 12CO and ATCA HI 21 cm observations of the LMC. The dust emissivity, gas-to-dust ratio, and X factor required to match the dust and gas surface densities are derived, and their correlations with the dust surface density are examined. The dust surface density is spatially correlated with the atomic and molecular gas phases. The dust temperature is consistently lower in the dense phase of the ISM than in the diffuse phase. We confirm variations in the ratio of FIR emission to gas surface density derived from HI and CO observations. There is an excess of FIR emission, spatially correlated with regions of intermediate HI and dust surface densities (Av = 1-2), and little or no CO. While there is no significant trend in the dust emissivity or gas-to-dust ratio with dust surface density, the X factor is enhanced at Av = 1-2. We conclude that H2 envelopes not traced by CO and X factor variations close to the CO boundary are more likely to cause these deviations between FIR emission and gas surface density than gas-to-dust ratio or emissivity variations.
[153]  [pdf] - 1032471
Determining dust temperatures and masses in the Herschel era: the importance of observations longward of 200 micron
Comments: 5 pages, 3 figurers, A&A, in press
Submitted: 2010-05-14
The properties of the dust grains (e.g., temperature and mass) can be derived from fitting far-IR SEDs (>100 micron). Only with SPIRE on Herschel has it been possible to get high spatial resolution at 200 to 500 micron that is beyond the peak (~160 micron) of dust emission in most galaxies. We investigate the differences in the fitted dust temperatures and masses determined using only <200 micron data and then also including >200 micron data (new SPIRE observations) to determine how important having >200 micron data is for deriving these dust properties. We fit the 100 to 350 micron observations of the Large Magellanic Cloud (LMC) point-by-point with a model that consists of a single temperature and fixed emissivity law. The data used are existing observations at 100 and 160 micron (from IRAS and Spitzer) and new SPIRE observations of 1/4 of the LMC observed for the HERITAGE Key Project as part of the Herschel Science Demonstration phase. The dust temperatures and masses computed using only 100 and 160 micron data can differ by up to 10% and 36%, respectively, from those that also include the SPIRE 250 & 350 micron data. We find that an emissivity law proportional to lambda^-1.5 minimizes the 100-350 micron fractional residuals. We find that the emission at 500 micron is ~10% higher than expected from extrapolating the fits made at shorter wavelengths. We find the fractional 500 micron excess is weakly anti-correlated with MIPS 24 micron flux and the total gas surface density. This argues against a flux calibration error as the origin of the 500 micron excess. Our results do not allow us to distinguish between a systematic variation in the wavelength dependent emissivity law or a population of very cold dust only detectable at lambda > 500 micron for the origin of the 500 micron excess.
[154]  [pdf] - 1026817
Mapping far-IR emission from the central kiloparsec of NGC 1097
Comments: Accepted for publication in the A&A Herschel Special Edition
Submitted: 2010-05-11
Using photometry of NGC 1097 from the Herschel PACS (Photodetector Array Camera and Spectrometer) instrument, we study the resolved properties of thermal dust continuum emission from a circumnuclear starburst ring with a radius ~ 900 pc. These observations are the first to resolve the structure of a circumnuclear ring at wavelengths that probe the peak (i.e. lambda ~ 100 micron) of the dust spectral energy distribution. The ring dominates the far-infrared (far-IR) emission from the galaxy - the high angular resolution of PACS allows us to isolate the ring's contribution and we find it is responsible for 75, 60 and 55% of the total flux of NGC 1097 at 70, 100 and 160 micron, respectively. We compare the far-IR structure of the ring to what is seen at other wavelengths and identify a sequence of far-IR bright knots that correspond to those seen in radio and mid-IR images. The mid- and far-IR band ratios in the ring vary by less than +/- 20% azimuthally, indicating modest variation in the radiation field heating the dust on ~ 600 pc scales. We explore various explanations for the azimuthal uniformity in the far-IR colors of the ring including a lack of well-defined age gradients in the young stellar cluster population, a dominant contribution to the far-IR emission from dust heated by older (> 10 Myr) stars and/or a quick smoothing of local enhancements in dust temperature due to the short orbital period of the ring. Finally, we improve previous limits on the far-IR flux from the inner ~ 600 pc of NGC 1097 by an order of magnitude, providing a better estimate of the total bolometric emission arising from the active galactic nucleus and its associated central starburst.
[155]  [pdf] - 1026799
Enhanced dust heating in the bulges of early-type spiral galaxies
Comments: Accepted for publication on the A&A Herschel Special Issue
Submitted: 2010-05-11
Stellar density and bar strength should affect the temperatures of the cool (T ~ 20-30 K) dust component in the inner regions of galaxies, which implies that the ratio of temperatures in the circumnuclear regions to the disk should depend on Hubble type. We investigate the differences between cool dust temperatures in the central 3 kpc and disk of 13 nearby galaxies by fitting models to measurements between 70 and 500 microns. We attempt to quantify temperature trends in nearby disk galaxies, with archival data from Spitzer/MIPS and new observations with Herschel/SPIRE, which were acquired during the first phases of the Herschel observations for the KINGFISH (key insights in nearby galaxies: a far-infrared survey with Herschel) sample. We fit single-temperature modified blackbodies to far-infrared and submillimeter measurements of the central and disk regions of galaxies to determine the temperature of the component(s) emitting at those wavelengths. We present the ratio of central-region-to-disk-temperatures of the cool dust component of 13 nearby galaxies as a function of morphological type. We find a significant temperature gradient in the cool dust component in all galaxies, with a mean center-to-disk temperature ratio of 1.15 +/- 0.03. The cool dust temperatures in the central ~3 kpc of nearby galaxies are 23(+/-3)% hotter for morphological types earlier than Sc, and only 9(+/-3)% hotter for later types. The temperature ratio is also correlated with bar strength, with only strongly barred galaxies having a ratio over 1.2. The strong radiation field in the high stellar density of a galactic bulge tends to heat the cool dust component to higher temperatures, at least in early-type spirals with relatively large bulges, especially when paired with a strong bar.
[156]  [pdf] - 1026794
The Herschel revolution: unveiling the morphology of the high mass star formation sites N44 and N63 in the LMC
Comments: Accepted for publication in A&A letters (Herschel special issue)
Submitted: 2010-05-11
We study the structure of the medium surrounding sites of high-mass star formation to determine the interrelation between the HII regions and the environment from which they were formed. The density distribution of the surroundings is key in determining how the radiation of the newly formed stars interacts with the surrounds in a way that allows it to be used as a star formation tracer. We present new Herschel/SPIRE 250, 350 and 500 mum data of LHA 120-N44 and LHA 120-N63 in the LMC. We construct average spectral energy distributions (SEDs) for annuli centered on the IR bright part of the star formation sites. The annuli cover ~10-~100 pc. We use a phenomenological dust model to fit these SEDs to derive the dust column densities, characterise the incident radiation field and the abundance of polycyclic aromatic hydrocarbon molecules. We see a factor 5 decrease in the radiation field energy density as a function of radial distance around N63. N44 does not show a systematic trend. We construct a simple geometrical model to derive the 3-D density profile of the surroundings of these two regions. Herschel/SPIRE data have proven very efficient in deriving the dust mass distribution. We find that the radiation field in the two sources behaves very differently. N63 is more or less spherically symmetric and the average radiation field drops with distance. N44 shows no systematic decrease of the radiation intensity which is probably due to the inhomogeneity of the surrounding molecular material and to the complex distribution of several star forming clusters in the region.
[157]  [pdf] - 1026755
Far-Infrared Line Imaging of the Starburst Ring in NGC 1097 with the Herschel/PACS Spectrometer
Comments: Accepted for publication on the A&A Herschel Special Issue
Submitted: 2010-05-10
NGC 1097 is a nearby SBb galaxy with a Seyfert nucleus and a bright starburst ring. We study the physical properties of the interstellar medium (ISM) in the ring using spatially resolved far-infrared spectral maps of the circumnuclear starburst ring of NGC 1097, obtained with the PACS spectrometer on board the Herschel Space Telescope. In particular, we map the important ISM cooling and diagnostic emission lines of [OI] 63 $\mu$m, [OIII] 88 $\mu$m, [NII] 122 $\mu$m, [CII] 158 $\mu$m and [NII] 205 $\mu$m. We observe that in the [OI] 63 $\mu$m, [OIII] 88 $\mu$m, and [NII] 122 $\mu$m line maps, the emission is enhanced in clumps along the NE part of the ring. We observe evidence of rapid rotation in the circumnuclear ring, with a rotation velocity of ~220$ km s$^{-1}$ (inclination uncorrected) measured in all lines. The [OI] 63 $\mu$m/[CII] 158 $\mu$m ratio varies smoothly throughout the central region, and is enhanced on the northeastern part of the ring, which may indicate a stronger radiation field. This enhancement coincides with peaks in the [OI] 63 $\mu$m and [OIII] 88 $\mu$m maps. Variations of the [NII] 122 $\mu$m/[NII] 205 $\mu$m ratio correspond to a range in the ionized gas density between 150 and 400 cm$^{-3}$.
[158]  [pdf] - 1025964
A Study of the Gas-Star Formation Relation over Cosmic Time
Comments: accepted for publication in MNRAS
Submitted: 2010-03-26, last modified: 2010-05-07
We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.
[159]  [pdf] - 169061
Spitzer Analysis of HII Region Complexes in the Magellanic Clouds: Determining a Suitable Monochromatic Obscured Star Formation Indicator
Comments: 21 pages, 12 figures, 4 tables. Will be published in ApJ.
Submitted: 2010-04-30
HII regions are the birth places of stars, and as such they provide the best measure of current star formation rates (SFRs) in galaxies. The close proximity of the Magellanic Clouds allows us to probe the nature of these star forming regions at small spatial scales. We aim to determine the monochromatic IR band that most accurately traces the bolometric IR flux (TIR), which can then be used to estimate an obscured SFR. We present the spatial analysis, via aperture/annulus photometry, of 16 LMC and 16 SMC HII region complexes using the Spitzer IRAC and MIPS bands. UV rocket data and SHASSA H-alpha data are also included. We find that nearly all of the LMC and SMC HII region SEDs peak around 70um, from ~10 to ~400 pc from the central sources. As a result, the sizes of HII regions as probed by 70um is approximately equal to the sizes as probed by TIR (about 70 pc in radius); the radial profile of the 70um flux, normalized by TIR, is constant at all radii (70um ~ 0.45 TIR); the 1-sigma standard deviation of the 70um fluxes, normalized by TIR, is a lower fraction of the mean (0.05 to 0.12 out to ~220 pc) than the normalized 8, 24, and 160um normalized fluxes (0.12 to 0.52); and these results are invariant between the LMC and SMC. From these results, we argue that 70um is the most suitable IR band to use as a monochromatic obscured star formation indicator because it most accurately reproduces the TIR of HII regions in the LMC and SMC and over large spatial scales. We also explore the general trends of the 8, 24, 70, and 160um bands in the LMC and SMC HII region SEDs, radial surface brightness profiles, sizes, and normalized (by TIR) radial flux profiles. We derive an obscured SFR equation that is modified from the literature to use 70um luminosity, SFR(Mo/yr) = 9.7(0.7)x10^{-44} L(70)(ergs/s), which is applicable from 10 to 300 pc distance from the center of an HII region.
[160]  [pdf] - 1025893
The Spitzer Survey of the Small Magellanic Cloud (S3MC): Insights into the Life-Cycle of Polycyclic Aromatic Hydrocarbons
Comments: accepted for publication in ApJ, abstract abridged
Submitted: 2010-03-23
We present the results of modeling dust SEDs across the SMC with the aim of mapping the distribution of PAHs in a low-metallicity environment. Using Spitzer Survey of the SMC (S3MC) photometry from 3.6-160 um over the main star-forming regions of the Wing and Bar along with spectral mapping from 5-38 um from the Spitzer Spectroscopic Survey of the SMC (S4MC) in selected regions, we model the dust SED and emission spectrum to determine the fraction of dust in PAHs across the SMC. We use the regions of overlapping photometry and spectroscopy to test the reliability of the PAH fraction as determined from SED fits alone. The PAH fraction in the SMC is low compared to the Milky Way and variable--with relatively high fractions (q_PAH~1-2%) in molecular clouds and low fractions in the diffuse ISM (<q_PAH>=0.6%). We use the map of PAH fraction across the SMC to test a number of ideas regarding the production, destruction and processing of PAHs in the ISM. We find weak or no correlation between the PAH fraction and the distribution of carbon AGB stars, the location of supergiant H I shells and young SN remnants, or the turbulent Mach number. We find that the PAH fraction is correlated with CO intensity, peaks in the dust surface density and the molecular gas surface density as determined from 160 um emission. The PAH fraction is high in regions of active star-formation, as predicted by its correlation with molecular gas, but is suppressed in H II regions. Because the PAH fraction in the diffuse ISM is generally very low--in accordance with previous work on modeling the SED of the SMC--and the PAH fraction is relatively high in molecular regions, we suggest that PAHs are destroyed in the diffuse ISM of the SMC and/or PAHs are forming in molecular clouds. We discuss the implications of these observations for our understanding of the PAH life cycle, particularly in low-metallicity galaxies.
[161]  [pdf] - 1025112
High molecular gas fractions in normal massive star forming galaxies in the young Universe
Comments: To Appear in Nature Feb 11, 2010; Supplementary Information included here
Submitted: 2010-02-10
Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.
[162]  [pdf] - 1018584
The SAURON project - XV. Modes of star formation in early-type galaxies and the evolution of the red sequence
Comments: Accepted for publication in MNRAS. 48 pages, 36 figures. Revision fixes typos, including column headers in Table 2. The revised paper with full resolution figures is available at
Submitted: 2009-12-01, last modified: 2009-12-21
We combine SAURON integral field data of a representative sample of local early-type, red sequence galaxies with Spitzer/IRAC imaging in order to investigate the presence of trace star formation in these systems. With the Spitzer data, we identify galaxies hosting low-level star formation, as traced by PAH emission, with measured star formation rates that compare well to those estimated from other tracers. This star formation proceeds according to established scaling relations with molecular gas content, in surface density regimes characteristic of disk galaxies and circumnuclear starbursts. We find that star formation in early-type galaxies happens exclusively in fast-rotating systems and occurs in two distinct modes. In the first, star formation is a diffuse process, corresponding to widespread young stellar populations and high molecular gas content. The equal presence of co- and counter-rotating components in these systems strongly implies an external origin for the star-forming gas, and we argue that these star formation events may be the final stages of (mostly minor) mergers that build up the bulges of red sequence lenticulars. In the second mode of star formation, the process is concentrated into well-defined disk or ring morphologies, outside of which the host galaxies exhibit uniformly evolved stellar populations. This implies that these star formation events represent rejuvenations within previously quiescent stellar systems. Evidence for earlier star formation events similar to these in all fast rotating early-type galaxies suggests that this mode of star formation may be common to all such galaxies, with a duty cycle of roughly 1/10, and likely contributes to the embedded, co-rotating inner stellar disks ubiquitous in this population.
[163]  [pdf] - 901743
Tracing the Bipolar Outflow from Orion Source I
Comments: 4 pages, 4 figures; accepted for publication in Ap J Letters
Submitted: 2009-09-15
Using CARMA, we imaged the 87 GHz SiO v=0 J=2-1 line toward Orion-KL with 0.45 arcsec angular resolution. The maps indicate that radio source I drives a bipolar outflow into the surrounding molecular cloud along a NE--SW axis, in agreement with the model of Greenhill et al. (2004). The extended high velocity outflow from Orion-KL appears to be a continuation of this compact outflow. High velocity gas extends farthest along a NW--SE axis, suggesting that the outflow direction changes on time scales of a few hundred years.
[164]  [pdf] - 1003481
Infrared Dark Clouds in the Small Magellanic Cloud?
Comments: accepted for publication in the Astronomical Journal
Submitted: 2009-08-16
We have applied the unsharp-masking technique to the 24 $\mu$m image of the Small Magellanic Cloud (SMC), obtained with the Spitzer Space Telescope, to search for high-extinction regions. This technique has been used to locate very dense and cold interstellar clouds in the Galaxy, particularly infrared dark clouds (IRDCs). Fifty five candidate regions of high-extinction, namely high-contrast regions (HCRs), have been identified from the generated decremental contrast image of the SMC. Most HCRs are located in the southern bar region and mainly distributed in the outskirts of CO clouds, but most likely contain a significant amount of H2. HCRs have a peak-contrast at 24 $\mu$m of 2 - 2.5 % and a size of 8 - 14 pc. This corresponds to the size of typical and large Galactic IRDCs, but Galactic IRDCs are 2 - 3 times darker at 24 $\mu$m than our HCRs. To constrain the physical properties of the HCRs, we have performed NH3, N2H+, HNC, HCO+, and HCN observations toward one of the HCRs, HCR LIRS36-EAST, using the Australia Telescope Compact Array and the Mopra single-dish radio telescope. We did not detect any molecular line emission, however, our upper limits to the column densities of molecular species suggest that HCRs are most likely moderately dense with n ~ 10^{3} cm-3. This volume density is in agreement with predictions for the cool atomic phase in low metallicity environments. We suggest that HCRs may be tracing clouds at the transition from atomic to molecule-dominated medium, and could be a powerful way to study early stages of gas condensation in low metallicity galaxies. Alternatively, if made up of dense molecular clumps < 0.5 pc in size, HCRs could be counterparts of Galactic IRDCs, and/or regions with highly unusual abundance of very small dust grains.
[165]  [pdf] - 26271
The Structure of a Low-Metallicity Giant Molecular Cloud Complex
Comments: 17 pages, 10 figures (including appendix), accepted for publication in the Astrophysical Journal
Submitted: 2009-07-13
To understand the impact of low metallicities on giant molecular cloud (GMC) structure, we compare far infrared dust emission, CO emission, and dynamics in the star-forming complex N83 in the Wing of the Small Magellanic Cloud. Dust emission (measured by Spitzer as part of the S3MC and SAGE-SMC surveys) probes the total gas column independent of molecular line emission and traces shielding from photodissociating radiation. We calibrate a method to estimate the dust column using only the high-resolution Spitzer data and verify that dust traces the ISM in the HI-dominated region around N83. This allows us to resolve the relative structures of H2, dust, and CO within a giant molecular cloud complex, one of the first times such a measurement has been made in a low-metallicity galaxy. Our results support the hypothesis that CO is photodissociated while H2 self-shields in the outer parts of low-metallicity GMCs, so that dust/self shielding is the primary factor determining the distribution of CO emission. Four pieces of evidence support this view. First, the CO-to-H2 conversion factor averaged over the whole cloud is very high 4-11 \times 10^21 cm^-2/(K km/s), or 20-55 times the Galactic value. Second, the CO-to-H2 conversion factor varies across the complex, with its lowest (most nearly Galactic) values near the CO peaks. Third, bright CO emission is largely confined to regions of relatively high line-of-sight extinction, A_V >~ 2 mag, in agreement with PDR models and Galactic observations. Fourth, a simple model in which CO emerges from a smaller sphere nested inside a larger cloud can roughly relate the H2 masses measured from CO kinematics and dust.
[166]  [pdf] - 26182
Radio Interferometric Planet Search I: First Constraints on Planetary Companions for Nearby, Low-Mass Stars from Radio Astrometry
Comments: accepted for publication in ApJ
Submitted: 2009-07-09
Radio astrometry of nearby, low-mass stars has the potential to be a powerful tool for the discovery and characterization of planetary companions. We present a Very Large Array survey of 172 active M dwarfs at distances of less than 10 pc. Twenty nine stars were detected with flux densities greater than 100 microJy. We observed 7 of these stars with the Very Long Baseline Array at milliarcsecond resolution in three separate epochs. With a detection threshold of 500 microJy in images of sensitivity 1 sigma ~ 100 microJy, we detected three stars three times (GJ 65B, GJ896A, GJ 4247), one star twice (GJ 285), and one star once (GJ 803). Two stars were undetected (GJ 412B and GJ 1224). For the four stars detected in multiple epochs, residuals from the optically-determined proper motions have an rms deviation of ~0.2 milliarcseconds, consistent with statistical noise limits. Combined with previous optical astrometry, these residuals provide acceleration upper limits that allow us to exclude planetary companions more massive than 3-6 M_Jup at a distance of ~1 AU with a 99% confidence level.
[167]  [pdf] - 26133
An HST/NICMOS view of the prototypical giant HII region NGC604 in M33
Comments: Paper presented in the Workshop "Young massive star clusters: initial conditions and environments" (Granada, Spain - Sept 2007). Astrophysics & Space Science in press, 7 pages, 4 figures
Submitted: 2009-07-09
We present the first high-spatial resolution near-infrared (NIR) imaging of NGC 604, obtained with the NICMOS camera aboard the Hubble Space Telescope (HST). These NICMOS broadband images reveal new NIR point sources, clusters, and diffuse structures. We found an excellent spatial correlation between the 8.4 GHz radio continuum and the 2.2mu-m nebular emission. Moreover, massive young stellar object candidates appear aligned with these radio peaks, reinforcing the idea that those areas are star-forming regions. Three different scaled OB associations are recognized in the NICMOS images. The brightest NIR sources in our images have properties that suggest that they are red supergiant stars, of which one of them was previously known. This preliminary analysis of the NICMOS images shows the complexity of the stellar content of the NGC 604 nebula.
[168]  [pdf] - 21851
How do Galaxies Accrete Gas and Form Stars?
Comments: 8 pages, Science white paper for the Astro2010 Decadal Survey, 2 color figures
Submitted: 2009-02-26
Great strides have been made in the last two decades in determining how galaxies evolve from their initial dark matter seeds to the complex structures we observe at z=0. The role of mergers has been documented through both observations and simulations, numerous satellites that may represent these initial dark matter seeds have been discovered in the Local Group, high redshift galaxies have been revealed with monstrous star formation rates, and the gaseous cosmic web has been mapped through absorption line experiments. Despite these efforts, the dark matter simulations that include baryons are still unable to accurately reproduce galaxies. One of the major problems is our incomplete understanding of how a galaxy accretes its baryons and subsequently forms stars. Galaxy formation simulations have been unable to accurately represent the required gas physics on cosmological timescales, and observations have only just begun to detect the star formation fuel over a range of redshifts and environments. How galaxies obtain gas and subsequently form stars is a major unsolved, yet tractable problem in contemporary extragalactic astrophysics. In this paper we outline how progress can be made in this area in the next decade.
[169]  [pdf] - 18625
The Dust-to-Gas Ratio in the Small Magellanic Cloud Tail
Comments: 6 pages, 3 figures, ApJ Letters, in press, (version with full resolution figures at
Submitted: 2008-11-17
The Tail region of the Small Magellanic Cloud (SMC) was imaged using the MIPS instrument on the Spitzer Space Telescope as part of the SAGE-SMC Spitzer Legacy. Diffuse infrared emission from dust was detected in all the MIPS bands. The Tail gas-to-dust ratio was measured to be 1200 +/- 350 using the MIPS observations combined with existing IRAS and HI observations. This gas-to-dust ratio is higher than the expected 500-800 from the known Tail metallicity indicating possible destruction of dust grains. Two cluster regions in the Tail were resolved into multiple sources in the MIPS observations and local gas-to-dust ratios were measured to be ~440 and ~250 suggests dust formation and/or significant amounts of ionized gas in these regions. These results support the interpretation that the SMC Tail is a tidal tail recently stripped from the SMC that includes gas, dust, and young stars.
[170]  [pdf] - 315195
Measuring Dust Production in the Small Magellanic Cloud Core-Collapse Supernova Remnant 1E 0102.2-7219
Comments: submitted to ApJ
Submitted: 2008-10-15, last modified: 2008-10-16
We present mid-infrared spectral mapping observations of the core-collapse supernova remnant 1E 0102.2-7219 in the Small Magellanic Cloud using the InfraRed Spectrograph (IRS) on the Spitzer Space Telescope. The remnant shows emission from fine structure transitions of neon and oxygen as well as continuum emission from dust. Comparison of the mid-IR dust emission with observations at x-ray, radio and optical wavelengths shows that the dust is associated with the supernova ejecta and is thus newly formed in the remnant. The spectrum of the newly formed dust is well reproduced by a model that includes 3x10^-3 solar masses of amorphous carbon dust at 70 K and 2x10^-5 solar masses of Mg2SiO4 (forsterite) at 145 K. Our observations place a lower limit on the amount of dust in the remnant since we are not sensitive to the cold dust in the unshocked ejecta. We compare our results to observations of other core-collapse supernovae and remnants, particularly Cas A where very similar spectral mapping observations have been carried out. We observe a a factor of ~10 less dust in E 0102 than seen in Cas A, although the amounts of amorphous carbon and forsterite are comparable.
[171]  [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.
[172]  [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).
[173]  [pdf] - 3933
HI and OH absorption in the lensing galaxy of MG J0414+0534
Comments: Accepted by MNRAS Letters, 5 (and a bit) pages, 5 figures
Submitted: 2007-08-12, last modified: 2007-09-12
We report the detection of \HI 21-cm absorption in the $z=0.96$ early-type lensing galaxy towards MG J0414+0534 with the Green Bank Telescope. The absorption, with total $N_{\rm HI}=1.6 \times 10^{18} (T_{\rm s}/f) {\rm cm}^{-2}$, is resolved into two strong components, probably due to the two strongest lens components, which are separated by 0.4\arcsec. Unlike the other three lenses which have been detected in \HI, J0414+0534 does not exhibit strong OH absorption, giving a OH/\HI column density ratio of $N_{\rm OH}/N_{\rm HI}\lapp10^{-6}$ (for $T_{\rm s}=100$ K, $T_{\rm x}=10$ K and $f_{\rm HI}=f_{\rm OH}=1$). This underabundance of molecular gas may indicate that the extreme optical--near-IR colour ($V-K=10.26$) along the line-of-sight is not due to the lens. We therefore suggest that despite the strong upper limits on molecular absorption at the quasar redshift, as traced by millimetre lines, the extinction occurs primarily in the quasar host galaxy.
[174]  [pdf] - 3448
The Spitzer Survey of the Small Magellanic Cloud: Discovery of Embedded Protostars in the HII Region NGC 346
Comments: 13 pages, 5 figures (3 in color), 2 tables. Accepted for publication in ApJ
Submitted: 2007-07-26
We use Spitzer Space Telescope observations from the Spitzer Survey of the Small Magellanic Cloud (S3MC) to study the young stellar content of N66, the largest and brightest HII region in the SMC. In addition to large numbers of normal stars, we detect a significant population of bright, red infrared sources that we identify as likely to be young stellar objects (YSOs). We use spectral energy distribution (SED) fits to classify objects as ordinary (main sequence or red giant) stars, asymptotic giant branch stars, background galaxies, and YSOs. This represents the first large-scale attempt at blind source classification based on Spitzer SEDs in another galaxy. We firmly identify at least 61 YSOs, with another 50 probable YSOs; only one embedded protostar in the SMC was reported in the literature prior to the S3MC. We present color selection criteria that can be used to identify a relatively clean sample of YSOs with IRAC photometry. Our fitted SEDs indicate that the infrared-bright YSOs in N66 have stellar masses ranging from 2 Msun to 17 Msun, and that approximately half of the objects are Stage II protostars, with the remaining YSOs roughly evenly divided between Stage I and Stage III sources. We find evidence for primordial mass segregation in the HII region, with the most massive YSOs being preferentially closer to the center than lower-mass objects. Despite the low metallicity and dust content of the SMC, the observable properties of the YSOs appear consistent with those in the Milky Way. Although the YSOs are heavily concentrated within the optically bright central region of N66, there is ongoing star formation throughout the complex and we place a lower limit on the star formation rate of 3.2 x 10^-3 Msun/yr over the last ~1 Myr.
[175]  [pdf] - 3120
Detection of $1.6\times 10^{10}$ M$_\odot$ of molecular gas in the host galaxy of the $z=5.77$ SDSS quasar J0927+2001
Comments: 12 pages, 2 figures, to appear in ApJ Letters
Submitted: 2007-07-16
We have detected emission by the CO 5-4 and 6-5 rotational transitions at $z = 5.7722\pm 0.0006$ from the host galaxy of the SDSS quasar J0927+2001 using the Plateau de Bure interferometer. The peak line flux density for the CO 5-4 line is $0.72 \pm 0.09$ mJy, with a line FWHM = $610 \pm 110$ km s$^{-1}$. The implied molecular gas mass is $(1.6 \pm 0.3) \times 10^{10}$ M$_\odot$. We also detect the 90 GHz continuum at $0.12 \pm 0.03$ mJy, consistent with a 47K dust spectrum extrapolated from higher frequencies. J0927+2001 is the second example of a huge molecular gas reservoir within the host galaxy of a quasar within 1 Gyr of the big bang. Observations of J0927+2001 are consistent with a massive starburst coeval with a bright quasar phase in the galaxy, suggesting the rapid formation of both a super-massive black hole through accretion, and the stellar host spheroid, at a time close to the end of cosmic reionization.
[176]  [pdf] - 3023
New insights on the dense molecular gas in NGC253 as traced by HCN and HCO+
Comments: 9 pages, 12 figures, ApJ in press (volume 666 September)
Submitted: 2007-07-12
We have imaged the central ~1kpc of the circumnuclear starburst disk in the galaxy NGC253 in the HCN(1-0), HCO+(1-0), and CO(1-0) transitions at 60pc resolution using the Owens Valley Radio Observatory Millimeter-Wavelength Array (OVRO). We have also obtained Atacama Pathfinder Experiment (APEX) observations of the HCN(4-3) and the HCO+(4-3) lines of the starburst disk. We find that the emission from the HCN(1-0) and HCO+(1-0) transitions, both indicators of dense molecular gas, trace regions which are non-distinguishable within the uncertainties of our observations. Even though the continuum flux varies by more than a factor 10 across the starburst disk, the HCN/HCO+ ratio is constant throughout the disk, and we derive an average ratio of 1.1+/-0.2. From an excitation analysis we find that all lines from both molecules are subthermally excited and that they are optically thick. This subthermal excitation implies that the observed HCN/HCO+ line ratio is sensitive to the underlying chemistry. The constant line ratio thus implies that there are no strong abundance gradients across the starburst disk of NGC253. This finding may also explain the variations in L'(HCN)/L'(HCO+) between different star forming galaxies both nearby and at high redshifts.
[177]  [pdf] - 2273
A Parallactic Distance of 389 +24/-21 parsecs to the Orion Nebula Cluster from Very Long Baseline Array Observations
Comments: 10 pages, 4 figures, emulateapj, accepted to ApJ
Submitted: 2007-06-15
We determine the parallax and proper motion of the flaring, non-thermal radio star GMR A, a member of the Orion Nebula Cluster, using Very Long Baseline Array observations. Based on the parallax, we measure a distance of 389 +24/-21 parsecs to the source. Our measurement places the Orion Nebula Cluster considerably closer than the canonical distance of 480 +/- 80 parsecs determined by Genzel et al. (1981). A change of this magnitude in distance lowers the luminosities of the stars in the cluster by a factor of ~ 1.5. We briefly discuss two effects of this change--an increase in the age spread of the pre-main sequence stars and better agreement between the zero-age main-sequence and the temperatures and luminosities of massive stars.
[178]  [pdf] - 2067
Self-consistent massive disks in triaxial dark matter halos
Comments: ApJ, in press
Submitted: 2007-06-10
Galactic disks in triaxial dark matter halos become deformed by the elliptical potential in the plane of the disk in such a way as to counteract the halo ellipticity. We develop a technique to calculate the equilibrium configuration of such a disk in the combined disk-halo potential, which is based on the method of Jog (2000) but accounts for the radial variation in both the halo potential and the disk ellipticity. This crucial ingredient results in qualitatively different behavior of the disk: the disk circularizes the potential at small radii, even for a reasonably low disk mass. This effect has important implications for proposals to reconcile cuspy halo density profiles with low surface brightness galaxy rotation curves using halo triaxiality. The disk ellipticities in our models are consistent with observational estimates based on two-dimensional velocity fields and isophotal axis ratios.
[179]  [pdf] - 1494
SubmilliJansky Transients in Archival Radio Observations
Comments: Accepted for ApJ; full quality figures available at
Submitted: 2007-05-21
[ABRIDGED] We report the results of a 944-epoch survey for transient sources with archival data from the Very Large Array spanning 22 years with a typical epoch separation of 7 days. Observations were obtained at 5 or 8.4 GHz for a single field of view with a full-width at half-maximum of 8.6' and 5.1', respectively, and achieved a typical point-source detection threshold at the beam center of ~300 microJy per epoch. Ten transient sources were detected with a significance threshold such that only one false positive would be expected. Of these transients, eight were detected in only a single epoch. Two transients were too faint to be detected in individual epochs but were detected in two-month averages. None of the ten transients was detected in longer-term averages or associated with persistent emission in the deep image produced from the combination of all epochs. The cumulative rate for the short timescale radio transients above 370 microJy at 5 and 8.4 GHz is 0.07 < R < 40 deg^-2 yr^-1, where the uncertainty is due to the unknown duration of the transients, 20 min < t_char < 7 days. A two-epoch survey for transients will detect 1.5 +/- 0.4 transient per square degrees above a flux density of 370 microJy. Two transients are associated with galaxies at z=0.040 and z=0.249. These may be similar to the peculiar Type Ib/c radio supernova SN 1998bw associated with GRB 980428. Six transients have no counterparts in the optical or infrared (R=27, Ks=18). The hosts and progenitors of these transients are unknown.
[180]  [pdf] - 171
The Low CO Content of the Extremely Metal Poor Galaxy I Zw 18
Comments: 5 pages in emulateapj style, accepted by the Astrophysical Journal
Submitted: 2007-04-06
We present sensitive molecular line observations of the metal-poor blue compact dwarf I Zw 18 obtained with the IRAM Plateau de Bure interferometer. These data constrain the CO J=1-0 luminosity within our 300 pc (FWHM) beam to be L_CO < 1 \times 10^5 K km s^-1 pc^2 (I_CO < 1 K km s^-1), an order of magnitude lower than previous limits. Although I Zw 18 is starbursting, it has a CO luminosity similar to or less than nearby low-mass irregulars (e.g. NGC 1569, the SMC, and NGC 6822). There is less CO in I Zw 18 relative to its B-band luminosity, HI mass, or star formation rate than in spiral or dwarf starburst galaxies (including the nearby dwarf starburst IC 10). Comparing the star formation rate to our CO upper limit reveals that unless molecular gas forms stars much more efficiently in I Zw 18 than in our own galaxy, it must have a very low CO-to-H_2 ratio, \sim 10^-2 times the Galactic value. We detect 3mm continuum emission, presumably due to thermal dust and free-free emission, towards the radio peak.
[181]  [pdf] - 38
Radio Astrometric Detection and Characterization of Extra-Solar Planets: A White Paper Submitted to the NSF ExoPlanet Task Force
Submitted: 2007-04-02
The extraordinary astrometric accuracy of radio interferometry creates an important and unique opportunity for the discovery and characterization of exo-planets. Currently, the Very Long Baseline Array can routinely achieve better than 100 microarcsecond accuracy, and can approach 10 microarcsecond with careful calibration. We describe here RIPL, the Radio Interferometric PLanet search, a new program with the VLBA and the Green Bank 100 m telescope that will survey 29 low-mass, active stars over 3 years with sub-Jovian planet mass sensitivity at 1 AU. An upgrade of the VLBA bandwidth will increase astrometric accuracy by an order of magnitude. Ultimately, the colossal collecting area of the Square Kilometer Array could push astrometric accuracy to 1 microarcsecond, making detection and characterizaiton of Earth mass planets possible. RIPL and other future radio astrometric planet searches occupy a unique volume in planet discovery and characterization parameter space. The parameter space of astrometric searches gives greater sensitivity to planets at large radii than radial velocity searches. For the VLBA and the expanded VLBA, the targets of radio astrometric surveys are by necessity nearby, low-mass, active stars, which cannot be studied efficiently through the radial velocity method, coronagraphy, or optical interferometry. For the SKA, detection sensitivity will extend to solar-type stars. Planets discovered through radio astrometric methods will be suitable for characterization through extreme adaptive optics. The complementarity of radio astrometric techniques with other methods demonstrates that radio astrometry can play an important role in the roadmap for exoplanet discovery and characterization.
[182]  [pdf] - 87073
The Spitzer Survey of the Small Magellanic Cloud: FIR Emission and Cold Gas in the SMC
Comments: 24 pages, 16 figures, accepted by ApJ. A version with nice figures is available at
Submitted: 2006-11-22
We present new far infrared maps of the Small Magellanic Cloud (SMC) at 24, 70, and 160 microns obtained as part of the Spitzer Survey of the Small Magellanic Cloud (S3MC,Bolatto et al. 2006). These maps cover most of the active star formation in the SMC Bar and the more quiescent Wing. We combine our maps with literature data to derive the dust surface density across the SMC. We find a total dust mass of Mdust = 3 10^5 Msun, implying a dust-to-hydrogen ratio over the region studied of log D/H = -2.86, or 1-to-700, which includes H_2. Assuming the dust to trace the total gas column, we derive H_2 surface densities across the SMC. We find a total H_2 mass M_H2 = 3.2 10^7 Msun in a distribution similar to that of the CO, but more extended. We compare profiles of CO and H_2 around six molecular peaks and find that on average H_2 is more extended than CO by a factor of \sim 1.3. The implied CO-to-H_2 conversion factor over the whole SMC is XCO = 13 +/- 1 10^21 cm^-2 (K km/s)^-1. Over the volume occupied by CO we find a lower conversion factor, XCO = 6 +/- 1 10^21 cm^-2 (K km/s)^-1, which is still a few times larger than that found using virial mass methods. The molecular peaks have H_2 surface densities \Sigma_H2 \approx 180 +/- 30 Msun pc^-2, similar to those in Milky Way GMCs, and correspondingly low extinctions, A_V \sim 1 - 2 mag. To reconcile these measurements with predictions by the theory of photoionization-regulated star formation, which requires A_V \sim 6, the GMCs must be \sim 3 times smaller than our 46 pc resolution element. We find that for a given hydrostatic gas pressure, the SMC has a 2 - 3 times lower ratio of molecular to atomic gas than spiral galaxies. Combined with the lower mean densities in the SMC this may explain why this galaxy has only 10% of its gas in the molecular phase.
[183]  [pdf] - 316495
The Spitzer Survey of the Small Magellanic Cloud: S3MC Imaging and Photometry in the Mid- and Far-Infrared Wavebands
Comments: Accepted for publication in The Astrophysical Journal. Given the draconian figure file-size limits implemented in astro-ph, readers are encouraged to download the manuscript with full quality images from
Submitted: 2006-08-25
We present the initial results from the Spitzer Survey of the Small Magellanic Cloud (S3MC), which imaged the star-forming body of the Small Magellanic Cloud (SMC) in all seven MIPS and IRAC wavebands. We find that the F_8/F_24 ratio (an estimate of PAH abundance) has large spatial variations and takes a wide range of values that are unrelated to metallicity but anticorrelated with 24 um brightness and F_24/F_70 ratio. This suggests that photodestruction is primarily responsible for the low abundance of PAHs observed in star-forming low-metallicity galaxies. We use the S3MC images to compile a photometric catalog of ~400,000 mid- and far-infrared point sources in the SMC. The sources detected at the longest wavelengths fall into four main categories: 1) bright 5.8 um sources with very faint optical counterparts and very red mid-infrared colors ([5.8]-[8.0]>1.2), which we identify as YSOs. 2) Bright mid-infrared sources with mildly red colors (0.16<[5.8]-[8.0]<0.6), identified as carbon stars. 3) Bright mid-infrared sources with neutral colors and bright optical counterparts, corresponding to oxygen-rich evolved stars. And, 4) unreddened early B stars (B3 to O9) with a large 24 um excess. This excess is reminiscent of debris disks, and is detected in only a small fraction of these stars (<5%). The majority of the brightest infrared point sources in the SMC fall into groups one to three. We use this photometric information to produce a catalog of 282 bright YSOs in the SMC with a very low level of contamination (~7%).
[184]  [pdf] - 81043
Infrared Sources in the Small Magellanic Cloud: First Results
Comments: 4 pages, 2 figures, uses asp2004.sty. To appear in the proceedings of the conference "Infrared Diagnostics of Galaxy Evolution", Pasadena, CA, 14-16 November 2005, ed. R. Chary
Submitted: 2006-03-30
We have imaged the entire Small Magellanic Cloud (SMC), one of the two nearest star-forming dwarf galaxies, in all seven IRAC and MIPS bands. The low mass and low metallicity (1/6 solar) of the SMC make it the best local analog for primitive galaxies at high redshift. By studying the properties of dust and star formation in the SMC at high resolution, we can gain understanding of similar distant galaxies that can only be observed in much less detail. In this contribution, we present a preliminary analysis of the properties of point sources detected in the Spitzer Survey of the Small Magellanic Cloud (S^3MC). We find ~400,000 unresolved or marginally resolved sources in our IRAC images, and our MIPS 24 micron mosaic contains ~17,000 point sources. Source counts decline rapidly at the longer MIPS wavelengths. We use color-color and color-magnitude diagrams to investigate the nature of these objects, cross-correlate their positions with those of known sources at other wavelengths, and show examples of how these data can be used to identify interesting classes of objects such as carbon stars and young stellar objects. For additional examples of some of the questions that can be studied with these data, please see the accompanying contributions by the other members of our team. The mosaic images and point source catalogs we have made have been released to the public on our website (
[185]  [pdf] - 79632
Molecular Gas in the Low Metallicity, Star Forming Dwarf IC 10
Comments: 21 pages, 14 figures, Accepted to ApJ
Submitted: 2006-02-02
We present a complete survey of CO 1->0 emission in the Local Group dwarf irregular IC 10. The survey, conducted with the BIMA interferometer, covers the stellar disk and a large fraction of the extended HI envelope with the sensitivity and resolution necessary to detect individual giant molecular clouds (GMCs) at the distance of IC 10 (950 kpc). We find 16 clouds with a total CO luminosity of 1 x 10^6 K km s^-1 pc^2, equivalent to 4 x 10^6 Msun of molecular gas using the Galactic CO-to-H2 conversion factor. Observations with the ARO 12m find that BIMA may resolve out as much as 50% of the CO emission, and we estimate the total CO luminosity as 2.2 x 10^6 K km s^-1 pc^2. We measure the properties of 14 GMCs from high resolution OVRO data. These clouds are very similar to Galactic GMCs in their sizes, line widths, luminosities, and CO-to-H2 conversion factors despite the low metallicity of IC 10 (Z ~ 1/5 Zsun). Comparing the BIMA survey to the atomic gas and stellar content of IC 10 we find that most of the CO emission is coincident with high surface density HI. IC 10 displays a much higher star formation rate per unit molecular (H2) or total (HI+H2) gas than most galaxies. This could be a real difference or may be an evolutionary effect - the star formation rate may have been higher in the recent past.
[186]  [pdf] - 76367
Spitzer Space Telescope detection of the young supernova remnant 1E 0102.2-7219
Comments: 5 pages, accepted by ApJL
Submitted: 2005-09-26
We present infrared observations of the young, oxygen-rich supernova remnant 1E 0102.2-7219 (E0102) in the Small Magellanic Cloud, obtained with the Spitzer Space Telescope. The remnant is detected at 24 um but not at 8 or 70 um and has a filled morphology with two prominent filaments. We find evidence for the existence of up to 8x10^-4 Msolar of hot dust (T~120 K) associated with the remnant. Most of the hot dust is located in the central region of E0102 which appears significantly enhanced in infrared and radio continuum emission relative to the X--ray emission. Even if all of the hot dust was formed in the explosion of E0102, the estimated mass of dust is at least 100 times lower that what is predicted by some recent theoretical models.
[187]  [pdf] - 73972
High Excitation Molecular Gas in the Magellanic Clouds
Comments: 18 pages, 3 figures, 5 tables. To appear in ApJ
Submitted: 2005-06-23
We present the first survey of submillimeter CO 4-3 emission in the Magellanic Clouds. The survey is comprised of 15 6'x6' maps obtained using the AST/RO telescope toward the molecular peaks of the Large and Small Magellanic Clouds. We have used these data to constrain the physical conditions in these objects, in particular their molecular gas density and temperature. We find that there are significant amounts of molecular gas associated with most of these molecular peaks, and that high molecular gas temperatures are pervasive throughout our sample. We discuss whether this may be due to the low metallicities and the associated dearth of gas coolants in the Clouds, and conclude that the present sample is insufficient to assert this effect.
[188]  [pdf] - 71093
The Molecular ISM of Dwarf Galaxies on Kiloparsec Scales: A New Survey for CO in Northern, IRAS-detected Dwarf Galaxies
Comments: 28 pages, 14 figures, 5 tables, ApJ accepted
Submitted: 2005-02-15
We present a new survey for CO in dwarf galaxies using the Kitt Peak 12m telescope. We observed the central regions of 121 northern dwarfs with IRAS detections and no known CO emission. We detect CO in 28 of these galaxies and marginally detect another 16, increasing by about 50% the number of such galaxies known to have significant CO emission. The galaxies we detect are comparable in mass to the LMC, although somewhat brighter in CO and fainter in the FIR. Within dwarfs, we find that the CO luminosity, L_CO, is most strongly correlated with the K-band and the far infrared luminosities. There are also strong correlations with the radio continuum and B-band luminosities, and linear diameter. We suggest that L_CO and L_K correlate well because the stellar component of a galaxy dominates the midplane gravitational field and thus sets the pressure of the atomic gas, which controls the formation of H_2 from HI. We compare our sample with more massive galaxies and find that dwarfs and large galaxies obey the same relationship between CO and the 1.4 GHz radio continuum (RC) surface brightness. This relationship is well described by a Schmidt Law with Sigma_RC proportional to Sigma_CO^1.3. Therefore, dwarf galaxies and large spirals exhibit the same relationship between molecular gas and star formation rate (SFR). We find that this result is robust to moderate changes in the RC-to-SFR and CO-to-H_2 conversion factors. Our data appear to be inconsistent with large (order of magnitude) variations in the CO-to-H_2 conversion factor in the star forming molecular gas. [abridged]
[189]  [pdf] - 70120
Warm Molecular Gas Traced with CO J=7->6 in the Galaxy's Central 2 Parsecs: Dynamical Heating of the Circumnuclear Disk
Comments: 31 pages, including 5 figures, accepted for publication in ApJ
Submitted: 2005-01-02
We present an 11 arcsec resolution map of the central two parsecs of the Galaxy in the CO J =7->6 rotational transition. The CO emission shows rotation about Sgr A*, but also evidence for non-circular turbulent motion and a clumpy morphology. We combine our dataset with available CO measurements to model the physical conditions in the disk. We find that the molecular gas in the region is both warm and dense, with T~200-300 K, n_H2~50,000-70,000 cm^-3. The mass of warm molecular gas we measure in the central two parsecs is at least 2000 M_solar, about 20 times the UV-excited atomic gas mass, ruling out an UV heating scenario for the molecular material. We compare the available spectral tracers with theoretical models and conclude that molecular gas is heated with magneto-hydrodynamic shocks with v~10-20 kms and B~0.3-0.5 mG. Using the conditions derived with the CO analysis, we include the other important coolants--neutral oxygen and molecular hydrogen--to estimate the total cooling budget of the molecular material. We derive a mass to luminosity ratio of 2-3 M_solar/ L_solar, which is consistent with the total power dissipated via turbulent decay in 0.1 pc cells with v_rms~15 kms. These size and velocity scales are comparable to the observed clumping scale and the velocity dispersion. At this rate, the material near Sgr A* its dissipating its orbital energy on an orbital timescale, and cannot last for more than a few orbits. Our conclusions support a scenario in which the features near Sgr A* such as the CND and northern arm are generated by infalling clouds with low specific angular momentum.
[190]  [pdf] - 69439
High-Resolution Measurements of the Halos of Four Dark Matter-Dominated Galaxies: Deviations from a Universal Density Profile
Comments: 22 pages, 7 tables, 9 figures (7 in color; Figure 1 resolution degraded). Accepted for publication in ApJ
Submitted: 2004-12-01
We derive rotation curves for four nearby, low-mass spiral galaxies and use them to constrain the shapes of their dark matter density profiles. This analysis is based on high-resolution two-dimensional Halpha velocity fields of NGC 4605, NGC 5949, NGC 5963, and NGC 6689 and CO velocity fields of NGC 4605 and NGC 5963. In combination with our previous study of NGC 2976, the full sample of five galaxies contains density profiles that span the range from alpha_dm = 0 to alpha_dm = 1.20, where alpha_dm is the power law index describing the central density profile. The scatter in alpha_dm from galaxy to galaxy is 0.44, three times as large as in Cold Dark Matter (CDM) simulations, and the mean density profile slope is alpha_dm = 0.73, shallower than that predicted by the simulations. These results call into question the hypothesis that all galaxies share a universal dark matter density profile. We show that one of the galaxies in our sample, NGC 5963, has a cuspy density profile that closely resembles those seen in CDM simulations, demonstrating that while galaxies with the steep central density cusps predicted by CDM do exist, they are in the minority. In spite of these differences between observations and simulations, the relatively cuspy density profiles we find do not suggest that this problem represents a crisis for CDM. Improving the resolution of the simulations and incorporating additional physics may resolve the remaining discrepancies. We also find that four of the galaxies contain detectable radial motions in the plane of the galaxy. We investigate the hypothesis that these motions are caused by a triaxial dark matter halo, and place lower limits on the ellipticity of the orbits in the plane of the disk of 0.043 - 0.175.
[191]  [pdf] - 67022
Detection of 21 Centimeter HI Absorption at z = 0.78 in a Survey of Radio Continuum Sources
Comments: AASTeX, 6 pages, 2 figures, accepted by ApJ Letters
Submitted: 2004-08-27
We report the detection of a deep broad HI 21 cm absorption system at z = 0.78 toward the radio source [HB89] 2351+456 (4C+45.51) at z = 1.992. The HI absorption was identified in a blind spectral line survey conducted at the Green Bank Telescope spanning 0.63 < z < 1.10 toward a large sample of radio continuum sources. The HI column density is N(HI) = 2.35 x 10^19 (T_s/f) cm^-2, where T_s is the spin temperature and f is the continuum covering factor of the absorbing gas. For T_s/f > 8.5 K, this system is by definition a damped Ly alpha absorption system (N(HI) >= 2 x 10^20 cm^-2). The line is unusually broad, with a FWHM of 53 km/s and a full span of 163 km/s, suggesting a physically extended HI gas structure. Radio surveys identify damped Ly alpha systems in a manner that bypasses many of the selection effects present in optical/UV surveys, including dust extinction and the atmospheric cutoff for z < 1.65. The smooth broad profile of this HI 21 cm absorption system is similar to the z = 0.89 HI absorption toward PKS 1830-211, which suggests that the absorber toward [HB89] 2351+456 is also a gravitational lens and a molecular absorption system. However, very long baseline interferometry and Hubble Space Telescope observations show little evidence for gravitational lensing, and BIMA millimeter observations show no HCO+ (1-2) or HCN (1-2) absorption down to tau = 0.15 (3 sigma) in 5 km/s channels. Although this radio damped Ly alpha selection technique would include dusty, molecule-rich systems, [HB89] 2351+456 appears to be a ``vanilla'' HI 21 cm absorber.
[192]  [pdf] - 64424
Molecular fraction limits in damped Lyman-alpha absorption systems
Comments: Accepted by MNRAS Supplement available from
Submitted: 2004-04-26
We have used the Green Bank Telescope (GBT) and Berkeley-Illinois-Maryland Association (BIMA) array to search for redshifted millimetre absorption in a sample of damped Lyman-alpha absorption systems (DLAs). This brings the number of published systems searched from 18 to 30. In 17 cases we reach $3\sigma$ limits of $\tau\leq0.1$, which is a significant improvement over the previous searches and more than sufficient to detect the 4 known redshifted millimetre absorbers ($\tau\gapp1$). While the CO rotational (millimetre) column density limits obtained are weaker than the electronic (optical) limits, they may provide useful limits below the atmospheric cut-off for the Lyman and Werner \MOLH-bands in the UV ($z_{\rm abs}\lapp1.8$). Using a model for the DLA metallicity evolution combined with assumed HCO$^+$/\MOLH ~and CO/\MOLH \~conversion ratios, we use the molecular column density limits to calculate plausible \MOLH ~molecular fraction limits. Finally, we use these results to discuss the feasibility of detecting rotational CO transitions in DLAs with the next generation of large radio telescopes.
[193]  [pdf] - 63750
The Faintness of the 158 um [CII] Transition in the z=6.42 Quasar SDSS J1148+5251
Comments: Five pages, one figure, accepted by the Astrophysical Journal Letters
Submitted: 2004-03-23
We report the non-detection of the [CII] 157.74 um transition in the z=6.42 quasar SDSS J1148+5251 after 37.5 hours of integration with the James Clerk Maxwell Telescope. This transition is the main cooling line of the star-forming interstellar medium, and usually the brightest FIR line in galaxies. Our observed RMS of 1.3 mK in the Ta* scale translates to L([CII])<2.6 x 10^9 Lsun. Using a recent estimate of the far-infrared continuum of this quasar, we derive for SDSS J1148+5251 L([CII])/L(FIR)<5 x 10^-4, a ratio similar to that observed in local ultra-luminous infrared galaxies but considerably smaller than what is typical in nearby normal and starburst galaxies. This indicates that the small L([CII])/L(FIR) ratio observed locally in luminous far-infrared objects also persists at the highest redshifts.
[194]  [pdf] - 1536172
Molecular Gas in Nearby Dwarf Galaxies: Single Dish and Interferometric Results
Comments: 4 pages, 2 figures. To appear in the proceedings of the 4th Cologne-Bonn Zermatt Symposium "The Dense Interstellar Medium in Galaxies", which took place in September 2003
Submitted: 2003-11-11
We present results from an extensive survey of molecular gas in nearby dwarf galaxies. This survey consists of a single-dish search for CO emission at the UASO Kitt Peak 12m telescope and interferometric follow-up using the Berkeley-Illinois-Maryland Array. We have detected 41 new sources, approximately doubling the number of detections of such objects in the literature. In this paper we discuss two results from the survey that lead us to agree with recent suggestions that, above a metallicity of 12+log[O/H]~8.0, there appears to be little variation in the CO-to-H2 conversion factor for star-forming molecular gas. We find two major pieces of evidence for this result: 1) the virial mass to molecular mass ratio found from the interferometer observations in a variety of sources shows no clear trend with metallicity; and 2) dwarf galaxies appear to exhibit a CO-to-radio continuum relationship that is similar to that observed in larger, presumably metal-rich, galaxies leaving little room for a large systematic variation of the CO-to-H2 conversion factor with metallicity in the star-forming gas.
[195]  [pdf] - 60839
Dark Matter in Dwarf Galaxies: High Resolution Observations
Comments: 6 pages, 4 figures. To appear in the proceedings of IAU Symposium 220 "Dark Matter in Galaxies", which took place in July 2003
Submitted: 2003-11-11
We present observations and analysis of rotation curves and dark matter halo density profiles in the central regions of four nearby dwarf galaxies. This observing program has been designed to overcome some of the limitations of other rotation curve studies that rely mostly on longslit spectra. We find that these objects exhibit the full range of central density profiles between constant density and NFW halos. This result suggests that there is a distribution of central density slopes rather than a unique halo density profile.
[196]  [pdf] - 59856
Dark Matter in Dwarf Galaxies: Latest Density Profile Results
Comments: 10 pages, 4 figures (1 in color), uses newpasp.sty. To appear in the proceedings of the conference "Satellites and Tidal Streams", La Palma, Spain, 26-30 May 2003, eds. F. Prada, D. Martinez-Delgado, T. Mahoney
Submitted: 2003-10-07
We present high-resolution two-dimensional velocity fields in Halpha and CO of the nearby dwarf galaxy NGC 2976. Our observations were made at both higher spatial resolution (~75 pc) and higher velocity resolution (13 km/s in Halpha and 2 km/s in CO) than most previous studies. We show that NGC 2976 has a very shallow dark matter density profile, with rho(r) lying between rho ~ r^-0.3 and rho ~ r^0. We carefully test the effects of systematic uncertainties on our results, and demonstrate that well-resolved, two-dimensional velocity data can eliminate many of the systematic problems that beset longslit observations. We also present a preliminary analysis of the velocity field of NGC 5963, which appears to have a nearly NFW density profile.
[197]  [pdf] - 58556
A Giant Outburst at Millimeter Wavelengths in the Orion Nebula
Comments: 29 pages, 7 figures, accepted for publication in ApJ
Submitted: 2003-08-15
BIMA observations of the Orion nebula discovered a giant flare from a young star previously undetected at millimeter wavelengths. The star briefly became the brightest compact object in the nebula at 86 GHz. Its flux density increased by more than a factor of 5 on a timescale of hours, to a peak of 160 mJy. This is one of the most luminous stellar radio flares ever observed. Remarkably, the Chandra X-ray observatory was in the midst of a deep integration of the Orion nebula at the time of the BIMA discovery; the source's X-ray flux increased by a factor of 10 approximately 2 days before the radio detection. Follow-up radio observations with the VLA and BIMA showed that the source decayed on a timescale of days, then flared again several times over the next 70 days, although never as brightly as during the discovery. Circular polarization was detected at 15, 22, and 43 GHz, indicating that the emission mechanism was cyclotron. VLBA observations 9 days after the initial flare yield a brightness temperature Tb > 5 x 10^7 K at 15 GHz. Infrared spectroscopy indicates the source is a K5V star with faint Br gamma emission, suggesting that it is a weak-line T Tauri object. Zeeman splitting measurements in the infrared spectrum find B ~ 2.6 +/- 1.0 kG. The flare is an extreme example of magnetic activity associated with a young stellar object. These data suggest that short observations obtained with ALMA will uncover hundreds of flaring young stellar objects in the Orion region.
[198]  [pdf] - 57867
High-Resolution Measurements of the Dark Matter Halo of NGC 2976: Evidence for a Shallow Density Profile
Comments: 30 pages, 4 tables, 13 figures (7 in color; Figures 1 and 3 are low-resolution to save space). Accepted for publication in ApJ. Version with full-resolution figures available at (46 MB)
Submitted: 2003-07-08
We have obtained two-dimensional velocity fields of the dwarf spiral galaxy NGC 2976 in Halpha and CO. The high spatial (~75 pc) and spectral (13 km/s and 2 km/s, respectively) resolution of these observations, along with our multicolor optical and near-infrared imaging, allow us to measure the shape of the density profile of the dark matter halo with good precision. We find that the total (baryonic plus dark matter) mass distribution of NGC 2976 follows a rho_tot ~ r^(-0.27 +/- 0.09) power law out to a radius of 1.8 kpc, assuming that the observed radial motions provide no support. The density profile attributed to the dark halo is even shallower, consistent with a nearly constant density of dark matter over the entire observed region. A maximal disk fit yields an upper limit to the K-band stellar mass-to-light ratio (M*/L_K) of 0.09^{+0.15}_{-0.08} M_sun/L_sun,K (including systematic uncertainties), with the caveat that for M*/L_K > 0.19 M_sun/L_sun,K the dark matter density increases with radius, which is unphysical. Assuming 0.10 M_sun/L_sun,K < M*/L_K < 0.19 M_sun/L_sun,K, the dark matter density profile lies between rho_dm ~ r^-0.17 and rho_dm ~ r^-0.01. Therefore, independent of any assumptions about the stellar disk or the functional form of the density profile, NGC 2976 does not contain a cuspy dark matter halo. We also investigate some of the systematic effects that can hamper rotation curve studies, and show that 1) longslit rotation curves are far more vulnerable to systematic errors than two-dimensional velocity fields, 2) NGC 2976 contains large radial motions at small radii, and 3) the Halpha and CO velocity fields of NGC 2976 agree within their uncertainties. [slightly abridged]
[199]  [pdf] - 57225
Unusual CO Line Ratios and kinematics in the N83/N84 Region of the Small Magellanic Cloud
Comments: 29 pages, aastex preprint style, 7 B&W figures (some degraded to fit the file size constraints), accepted for publication in ApJ, retrieve paper with full resolution figures from
Submitted: 2003-06-06
We present new CO 1-0 and 2-1 observations of the N83/N84 molecular cloud complex in the south-east Wing of the Small Magellanic Cloud (SMC). While the 2-1/1-0 integrated line brightness ratio (in temperature units) is uniformly 0.9 throughout most of the complex, we find two distinct regions with unusually high ratios 2-1/1-0 > 2. These regions are associated with the N84D nebula and with the inside of the 50 pc expanding molecular shell N83. This shell is spatially coincident with the NGC 456 stellar association and the HFPK2000-448 radio continuum/X-ray source tentatively classified as a supernova remnant. We explore possible causes for the high ratios observed and conclude that the CO emission probably arises from an ensemble of small (R~0.1 pc), warm (T~40 K) clumps. Analysis of the CO shell parameters suggests that it is wind--driven and has an age of slightly more than 2 million years. We have also used this dataset to determine the CO-to-H2 conversion factor in the SMC, an especially interesting measurement because of the low metallicity of this source (~1/9 solar). Surprisingly, after comparing the CO luminosities of clouds in N83/N84 with their virial masses, we find a CO-to-H2 conversion factor Xco only 1.9 times larger than what we obtain when applying the same algorithm to solar metallicity clouds in the Milky Way and M 33. This result fits into the emerging pattern that CO observations with high linear resolution suggest nearly Galactic values of Xco in a wide range of environments.
[200]  [pdf] - 53661
CO (J=7->6) Observations of NGC 253: Cosmic Ray Heated Warm Molecular Gas
Comments: 12+ pages, 4 figures, to appear in The Astrophysical Journal (now includes tables)
Submitted: 2002-12-11, last modified: 2003-01-13
We report observations of the CO J=7->6 transition toward the starburst nucleus of NGC 253. This is the highest-excitation CO measurement in this source to date, and allows an estimate of the molecular gas excitation conditions. Comparison of the CO line intensities with a large velocity gradient, escape probability model indicates that the bulk of the 2-5 x 10^7 solar masses of molecular gas in the central 180 pc is highly excited. A model with T ~ 120 K, n_H_2 ~ 4.5 x 10^4 cm^-3 is consistent with the observed CO intensities as well as the rotational H2 lines observed with ISO. The inferred mass of warm, dense molecular gas is 10--30 times the atomic gas mass as traced through its [CII] and [OI] line emission. This large mass ratio is inconsistent with photodissociation region models where the gas is heated by far-UV starlight. It is also not likely that the gas is heated by shocks in outflows or cloud-cloud collisions. We conclude that the best mechanism for heating the gas is cosmic rays, which provide a natural means of uniformly heating the full volume of molecular clouds. With the tremendous supernova rate in the nucleus of NGC 253, the CR heating rate is at least ~800 times greater than in the Galaxy, more than sufficient to match the cooling observed in the CO lines.
[201]  [pdf] - 49233
SPIFI: a direct-detection imaging spectrometer for submillimeter wavelengths
Comments: 36 pages with 10 figures
Submitted: 2002-05-10
The South Pole Imaging Fabry-Perot Interferometer (SPIFI) is the first instrument of its kind -- a direct-detection imaging spectrometer for astronomy in the submillimeter band. SPIFI's focal plane is a square array of 25 silicon bolometers cooled to 60 mK; the spectrometer consists of two cryogenic, scanning Fabry-Perot interferometers in series with a 60 mK bandpass filter. The instrument operates in the short submillimeter windows (350 microns, 450 microns) available from the ground, with spectral resolving power selectable between 500 and 10000. At present, SPIFI's sensitivity is within a factor 1.5-3 of the photon background limit, comparable to the best heterodyne spectrometers. The instrument's large bandwidth and mapping capability provide substantial advantages for specific astrophysical projects, including deep extragalactic observations. In this article we present the motivation for and design of SPIFI and its operational characteristics on the telescope.
[202]  [pdf] - 44026
The Density Profile of the Dark Matter Halo of NGC 4605
Comments: 8 pages, 5 figures. Revised version (minor changes), accepted for publication in The Astrophysical Journal. Also available at
Submitted: 2001-08-02, last modified: 2001-08-30
We have obtained ~100 pc resolution CO and ~60 pc resolution Halpha observations of the dwarf spiral galaxy NGC 4605. We use them to derive a high resolution rotation curve and study the central density profile of NGC 4605's dark matter halo. We find that these observations do not agree with the predictions of most high resolution Cold Dark Matter calculations. We investigate two extreme cases: 1) NGC 4605 has a maximal exponential disk, that we model using K-band observations and remove to study the structure of its dark matter halo, and 2) NGC 4605 is dark matter dominated and its disk is dynamically negligible. Because the mass-to-light ratio of the maximal disk is already very low we favor the first solution, which indicates the halo has one component with a density profile rho ~ r^{-0.6} out to R~2.8 kpc. In the second case, the rotation curve requires the presence of two components: a small (~600 pc) core surrounded by a much steeper rho ~ r^{-1.1} halo. Removal of intermediate (submaximal) disks does not ameliorate the discrepancy between the predictions and the observations.
[203]  [pdf] - 1468202
CO (J = 4-3) and [C I] Observations of the Carina Molecular Cloud Complex
Comments: 27 text pages, 15 figures
Submitted: 2001-01-16
We present large area, fully-sampled maps of the Carina molecular cloud complex in the CO (J = 4 - 3) and neutral carbon [CI] 3P1 - 3P0 transitions. These data were obtained using the 1.7 meter diameter Antarctic Submillimeter Telescope and Remote Observatory (AST/RO). The maps cover an area of approximately 3 square degrees with a uniform 1' spatial sampling. Analysis of these data, in conjunction with CO (J = 1 - 0) data from the Columbia CO survey and the IRAS HIRES continuum maps for the same region, suggests that the spiral density wave shock associated with the Carina spiral arm may be playing an important role in the formation and dissociation of the cloud complex, as well as in maintaining the internal energy balance of the clouds in this region. Massive stars form at the densest regions of the molecular cloud complex. The winds and outflows associated with these stars have a disrupting effect on the complex and inject mechanical energy into the parent clouds, while the UV radiation from the young stars also heat the parent clouds. The present set of data suggests, however, that massive stars alone may not account for the energetics of the clouds in the Carina region. The details of the data and the correlation among the various data sets hint at the possible role that the spiral density wave shock plays in feeding interstellar turbulence and in heating molecular clouds.
[204]  [pdf] - 37569
The Antarctic Submillimeter Telescope and Remote Observatory (AST/RO)
Comments: 57 pages, 15 figures. Submitted to PASP
Submitted: 2000-08-16, last modified: 2000-12-15
AST/RO, a 1.7 m diameter telescope for astronomy and aeronomy studies at wavelengths between 200 and 2000 microns, was installed at the South Pole during the 1994-1995 Austral summer. The telescope operates continuously through the Austral winter, and is being used primarily for spectroscopic studies of neutral atomic carbon and carbon monoxide in the interstellar medium of the Milky Way and the Magellanic Clouds. The South Pole environment is unique among observatory sites for unusually low wind speeds, low absolute humidity, and the consistent clarity of the submillimeter sky. Four heterodyne receivers, an array receiver, three acousto-optical spectrometers, and an array spectrometer are installed. A Fabry-Perot spectrometer using a bolometric array and a Terahertz receiver are in development. Telescope pointing, focus, and calibration methods as well as the unique working environment and logistical requirements of the South Pole are described.
[205]  [pdf] - 37755
AST/RO Observations of Atomic Carbon near the Galactic Center
Comments: 11 pages, 6 figures
Submitted: 2000-08-28
We present a coarsely-sampled map of the region |l|<=2 degrees, |b|<=0.1 degrees in the 492 GHz 3P1-->3P0 fine structure transition of neutral carbon, observed with the Antarctic Submillimeter Telescope and Remote Observatory (AST/RO). The distribution of C[I] emission is similar on the large scale to that of CO J=1-->0. On average, the ratio of the integrated intensities, I_{C I}/I_{12CO}, is higher in the Galactic disk than in the Galactic Center region. This result is accounted for by the absorption of 12CO within the clouds located in the outer Galactic disk. I_{C I}/I_{12CO} is surprisingly uniform over the variety of environments near the Galactic Center. On average, C[I] is optically thin (or as optically thin as 13CO(J=1-->0)) even in the dense molecular clouds of the Galactic Center region.
[206]  [pdf] - 37371
First Detection of Submillimeter [C I] Emission in the Small Magellanic Cloud
Comments: 12 pages including 1 table and 2 figures, using AASTEX macros v4.0, to appear in The Astrophysical Journal Letters
Submitted: 2000-08-02
We report the first detection of [CI] (3P1-3P0) emission at 609 um in a region of the Small Magellanic Cloud (N27). Environments poor in heavy elements and dust such as the SMC are thought to be dominated by photodissociation regions. This is the lowest metallicity source where submillimeter neutral carbon emission has been detected. Studying the [CI]/CO intensity ratio in several sources spanning more than an order of magnitude in metallicity, Z, we find that the [CI]/CO ratio increases for decreasing Z. The existence of such a trend points to a photodissociation origin for most of neutral carbon in molecular clouds, in agreement with standard PDR models. We also report ISO FIR spectroscopic observations of N27, and use them to derive its physical properties. Comparison between the density and radiation field revealed by FIR diagnostics (n~300-1000 cm-3, Xuv~30-100 Xo) and those derived from millimeter and submillimeter data (n~10^5 cm-3, Xuv<30 Xo) suggests that the FIR lines originate in more diffuse gas, and are perhaps dominated by the interclump medium. Regardless of the cause, analysis of the FIR and mm-submm data produces a discrepancy of two orders of magnitude for the density of this source.
[207]  [pdf] - 36987
Carbon in the N159/N160 Complex of the Large Magellanic Cloud
Comments: 41 pages including 5 tables and 9 figures, using AASTEX macros v4.0, because of space restrictions some figures are reduced quality bitmaps, download a preprint with publication quality figures from, to appear in The Astrophysical Journal
Submitted: 2000-07-11
We present a study of carbon in N159/N160, an HII region complex in the low metallicity Large Magellanic Cloud. We have mapped this region, which comprises four distinct molecular clouds spanning a wide range of star-formation activity, in four transitions: 13CO (1-0) CO (2-1) and (4-3), and [CI] (1-0). Combining these data with existing [CII] observations provides a complete picture of the predominant forms of carbon in the gas phase of the ISM. The new CO (2-1) data show that the complex is immersed in an envelope of extended, low level emission, undetected by previous (1-0) mapping efforts. The CO (2-1)/(1-0) ratio in this envelope is >~3, a value consistent with optically thin CO emission. The envelope is also relatively bright in [CI] and [CII], and calculations show that it is mostly photodissociated: it appears to be translucent (Av<1). Neutral carbon emission in the complex unexpectedly peaks at the quiescent southern cloud (N159S). In the northern portion of the map (the N160 nebula), the HII regions prominent in [CII] correspond to holes in the [CI] distribution. Overall we find that, while the [CII]/CO ratio is enhanced with respect to similar complexes in the Milky Way, the [CI]/CO ratio appears to be similar or reduced.
[208]  [pdf] - 109118
Submillimeter observations of IC 10: the dust properties and neutral carbon content of a low metallicity starburst
Comments: 31 pages, 9 figures, 5 tables Prepared using AASLaTeX macros v4.0 Accepted for publication in the Astrophysical Journal
Submitted: 1999-11-01
We present submillimeter observations of the Local Group, metal-poor, irregular dwarf galaxy IC 10, directly relevant to the interaction between interstellar medium and star formation activity in primeval galaxies. Using the JCMT, we have observed the fine structure neutral carbon transition at 492 GHz and the rotational J=3-2 transition of 12CO and 13CO in the most massive giant molecular cloud complex in this galaxy, IC10-SE. We find that, although the I([CII])/I(CO) ratio for this object is a factor of 4 larger than the typical Milky Way value, its [C~I] to CO intensity ratio is I([CI])/I(CO)~18 is similar (only about 50% larger) to that of the Milky Way. Modelling of the [CII]/CO and [CI]/CO intensity ratios with metallicity indicates that, if C+ and C are chiefly produced by UV photodissociation, both ratios should increase sharply with decreasing metallicity (and consequently diminished UV shielding; Bolatto, Jackson, and Ingalls 1999). These data then suggest a different origin for an important fraction of C in these clouds, unrelated to photodissociation. We have also mapped the 850 um continuum in this region using SCUBA. Employing these data in conjunction with KAO and IRAM measurements we find an extremely low emissivity exponent, b~0.5. We conclude that this low exponent is most likely due to the destruction of small dust grains, brought about by the increased penetration of UV radiation in the low metallicity ISM. If a low b in the submillimeter is a general property of metal-poor systems then the interpretation of millimeter and submillimeter surveys of high-z galaxies should be revised.
[209]  [pdf] - 104288
A semianalytical model for the observational properties of the dominant carbon species at different metallicities
Comments: 28 pages, 1 table, 7 Postscript figures, to appear in ApJ vol. 513, March 1999
Submitted: 1998-12-09
We present a simple, semianalytical approach to compute the line emission of C+, C and CO in Photodissociation Regions of varying metallicity that depends on very few parameters and naturally incorporates the clumpiness of the interstellar medium into the calculations. We use it to study the dependence of the I_[C II]/I_CO and I_[C I]/I_CO line ratios on the metallicity. We show how to modify it to include the effects of density and radiation field, and we have compared it with observational data. We find that the model explains the observed trend of enhanced [C II]/CO line ratio with decreasing metallicity as the natural result of the augmented fraction of photodissociated gas in a clump. We show that enhanced [C II]/CO ratios can be produced by lowering the upper limit of the clump size distribution, as may happen in 30 Doradus. We also find that the available data favors a [C I]/CO intensity ratio essentially independent of metallicity, albeit the paucity of observations does not exclude other possibilities. This is difficult to understand if most of the C is produced by UV photons as is the case with C+. Finally, we study the prediction of the model for the trend of the X_CO factor with metallicity. Comparison with previous observational studies yields good agreement.
[210]  [pdf] - 91616
The detectability of planetary companions of compact galactic objects from their effects on microlensed lightcurves of distant stars
Comments: 11 pages, Format: compressed and uuencoded Poscript file (for recovering the original Poscript file run first "uudecode filename" where 'filename' is the name you gave to the file where you stored this message, then "uncompress". A file named '' will appear. In case of problem please contact, Harvard-Smithsonian Center for Astrophysics preprint# 3891 (to appear in The Astrophysical Journal of November 20, 1994)
Submitted: 1994-07-11
We discuss a possible method for detection of dark companions of galactic objects of stellar mass. Such binary systems are likely to occur in the galactic disk and possibly also in the halo. The high incidence of binary and higher-multiplicity systems in the solar neighborhood, if indicative of the galactic disk at large, implies that current searches for the gravitational microlensing signature of massive compact objects in our galaxy would yield a significant fraction of binary systems. Our calculations suggest that 40% of the lightcurves that will be obtained in such searches may be sufficiently perturbed to reveal, if sufficiently well-sampled, the presence of a compact dark companion of Jovian mass orbiting a primary. The likelihood of occurrence of perturbed lightcurves must also be taken into account by systematic search programs, to improve their event detection efficiency. The statistics of perturbed lensing events, if carefully interpreted, could yield estimates of the incidence of binary systems with low mass ratios, including that of systems with planets.