Source author record

Alberto D. Bolatto

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

62works
9topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

62 published item(s)

preprint2023arXiv

X-ray Properties of NGC 253's Starburst-Driven Outflow

We analyze image and spectral data from $\approx$365~ks of observations from the {\it Chandra} X-ray Observatory of the nearby, edge-on starburst galaxy NGC 253 to constrain properties of the hot phase of the outflow. We focus our analysis on the $-$1.1 to $+$0.63 kpc region of the outflow and define several regions for spectral extraction where we determine best-fit temperatures and metal abundances. We find that the temperatures and electron densities peak in the central $\sim$250 pc region of the outflow and decrease with distance. These temperature and density profiles are in disagreement with an adiabatic spherically expanding starburst wind model and suggest the presence of additional physics such as mass loading and non-spherical outflow geometry. Our derived temperatures and densities yield few-Myr cooling times in the nuclear region, which may imply that the hot gas can undergo bulk radiative cooling as it escapes along the minor axis. Our metal abundances of O, Ne, Mg, Si, S, and Fe all peak in the central region and decrease with distance along the outflow, with the exception of Ne which maintains a flat distribution. The metal abundances indicate significant dilution outside of the starburst region. We also find estimates on the mass outflow rates which are $2.8\:M_{\odot}/\rm{yr}$ in the northern outflow and $3.2\:M_{\odot}/\rm{yr}$ in the southern outflow. Additionally, we detect emission from charge exchange and find it has a significant contribution ($20-42$\%) to the total broad-band ($0.5-7$ keV) X-ray emission in the central and southern regions of the outflow.

preprint2022arXiv

Cuspy dark matter density profiles in massive dwarf galaxies

Rotation curves of galaxies probe their total mass distributions, including dark matter. Dwarf galaxies are excellent systems to investigate the dark matter density distribution, as they tend to have larger fractions of dark matter compared to higher mass systems. The core-cusp problem describes the discrepancy found in the slope of the dark matter density profile in the centres of galaxies ($β^*$) between observations of dwarf galaxies (shallower cores) and dark matter-only simulations (steeper cusps). We investigate $β^*$ in six nearby spiral dwarf galaxies for which high-resolution CO $J=1-0$ data were obtained with ALMA. We derive rotation curves and decompose the mass profile of the dark matter using our CO rotation curves as a tracer of the total potential and 4.5$μ$m photometry to define the stellar mass distribution. We find $\langleβ^*\rangle = 0.6$ with a standard deviation of $\pm0.1$ among the galaxies in this sample, in agreement with previous measurements in this mass range. The galaxies studied are on the high stellar mass end of dwarf galaxies and have cuspier profiles than lower mass dwarfs, in agreement with other observations. When the same definition of the slope is used, we observe steeper slopes than predicted by the FIRE and NIHAO simulations. This may signal that these relatively massive dwarfs underwent stronger gas inflows toward their centres than predicted by these simulations, that these simulations over-predict the frequency of accretion or feedback events, or that a combination of these or other effects are at work.

preprint2022arXiv

Developing a New Generation of Integrated Micro-Spec Far Infrared Spectrometers for the EXperiment for Cryogenic Large-Aperture Intensity Mapping (EXCLAIM)

The current state of far-infrared astronomy drives the need to develop compact, sensitive spectrometers for future space and ground-based instruments. Here we present details of the $\rm μ$-Spec spectrometers currently in development for the far-infrared balloon mission EXCLAIM. The spectrometers are designed to cover the $\rm 555 - 714\ μ$m range with a resolution of $\rm R\ =\ λ/ Δλ =\ 512$ at the $\rm 638\ μ$m band center. The spectrometer design incorporates a Rowland grating spectrometer implemented in a parallel plate waveguide on a low-loss single-crystal Si chip, employing Nb microstrip planar transmission lines and thin-film Al kinetic inductance detectors (KIDs). The EXCLAIM $\rm μ$-Spec design is an advancement upon a successful $\rm R = 64\ μ$-Spec prototype, and can be considered a sub-mm superconducting photonic integrated circuit (PIC) that combines spectral dispersion and detection. The design operates in a single $M{=}2$ grating order, allowing one spectrometer to cover the full EXCLAIM band without requiring a multi-order focal plane. The EXCLAIM instrument will fly six spectrometers, which are fabricated on a single 150 mm diameter Si wafer. Fabrication involves a flip-wafer-bonding process with patterning of the superconducting layers on both sides of the Si dielectric. The spectrometers are designed to operate at 100 mK, and will include 355 Al KID detectors targeting a goal of NEP ${\sim}8\times10^{-19}$ $\rm W/\sqrt{Hz}$. We summarize the design, fabrication, and ongoing development of these $\rm μ$-Spec spectrometers for EXCLAIM.

preprint2022arXiv

Effects of CO-dark Gas on Measurements of Molecular Cloud Stability and the Size-Linewidth Relationship

Stars form within molecular clouds, so characterizing the physical states of molecular clouds is key in understanding the process of star formation. Cloud structure and stability is frequently assessed using metrics including the virial parameter and Larson (1981) scaling relationships between cloud radius, velocity dispersion, and surface density. Departures from the typical Galactic relationships between these quantities have been observed in low-metallicity environments. The amount of H$_2$ gas in cloud envelopes without corresponding CO emission is expected to be high under these conditions; therefore, this "CO-dark" gas could plausibly be responsible for the observed variations in cloud properties. We derive simple corrections that can be applied to empirical clump properties (mass, radius, velocity dispersion, surface density, and virial parameter) to account for CO-dark gas in clumps following power-law and Plummer mass density profiles. We find that CO-dark gas is not likely to be the cause of departures from Larson's relationships in low-metallicity regions, but that virial parameters may be systematically overestimated. We demonstrate that correcting for CO-dark gas is critical for accurately comparing the dynamical state and evolution of molecular clouds across diverse environments.

preprint2022arXiv

Snowmass2021 Cosmic Frontier White Paper: Observational Facilities to Study Dark Matter

We present an overview of future observational facilities that will significantly enhance our understanding of the fundamental nature of dark matter. These facilities span a range of observational techniques including optical/near-infrared imaging and spectroscopy, measurements of the cosmic microwave background, pulsar timing, 21-cm observations of neutral hydrogen at high redshift, and the measurement of gravitational waves. Such facilities are a critical component of a multi-pronged experimental program to uncover the nature of dark matter, while often providing complementary measurements of dark energy, neutrino physics, and inflation.

preprint2022arXiv

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

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

preprint2021arXiv

HI mapping of the Leo Triplet: Morphologies and kinematics of tails and bridges

A fully-sampled and hitherto highest resolution and sensitivity observation of neutral hydrogen (HI) in the Leo Triplet (NGC 3628, M 65/NGC 3623, and M 66/NGC 3627) reveals six HI structures beyond the three galaxies. We present detailed results of the morphologies and kinematics of these structures, which can be used for future simulations. In particular, we detect a two-arm structure in the plume of NGC 3628 for the first time, which can be explained by a tidal interaction model. The optical counterpart of the plume is mainly associated with the southern arm. The connecting part (base) of the plume (directed eastwards) with NGC 3628 is located at the blueshifted (western) side of NGC 3628. Two bases appear to be associated with the two arms of the plume. A clump with reversed velocity gradient (relative to the velocity gradient of M 66) and a newly detected tail, i.e. M 66SE, is found in the southeast of M 66. We suspect that M 66SE represents gas from NGC 3628 which was captured by M 66 in the recent interaction between the two galaxies. Meanwhile gas is falling toward M 66, resulting in features already previously observed in the southeastern part of M 66, e.g. large line widths and double peaks. An upside-down `Y'-shaped HI gas component (M 65S) is detected in the south of M 65 which suggests that M 65 may also have been involved in the interaction. We strongly encourage modern hydrodynamical simulations of this interacting group of galaxies to reveal the origin of the gaseous debris surrounding all three galaxies.

preprint2021arXiv

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

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

preprint2021arXiv

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

The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is critical to studying molecular gas and star formation in galaxies. The value of $α_\rm{CO}$ has been found to vary within and between galaxies, but the specific environmental conditions that cause these variations are not fully understood. Previous observations on $\sim$kpc scales revealed low values of $α_\rm{CO}$ in the centers of some barred spiral galaxies, including NGC 3351. We present new ALMA Band 3, 6, and 7 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales in the inner $\sim$2 kpc of NGC 3351. Using multi-line radiative transfer modeling and a Bayesian likelihood analysis, we infer the H$_2$ density, kinetic temperature, CO column density per line width, and CO isotopologue abundances on a pixel-by-pixel basis. Our modeling implies the existence of a dominant gas component with a density of $2{-}3\times10^3$ $\rm{cm^{-3}}$ in the central ${\sim}$1 kpc and a high temperature of 30$-$60 K near the nucleus and near the contact points that connect to the bar-driven inflows. Assuming a CO/H$_2$ abundance of $3\times10^{-4}$, our analysis yields $α_\rm{CO}{\sim}0.5{-}2.0$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$ with a decreasing trend with galactocentric radius in the central $\sim$1 kpc. The inflows show a substantially lower $α_\rm{CO} < 0.1$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, likely due to lower optical depths caused by turbulence or shear in the inflows. Over the whole region, this gives an intensity-weighted $α_\rm{CO}$ of ${\sim}1.5$ $\rm{M_\odot\,(K~km~s^{-1}~pc^2)^{-1}}$, which is similar to previous dust modeling based results at kpc scales. This suggests that low $α_\rm{CO}$ on kpc scales in the centers of some barred galaxies may be due to the contribution of low optical depth CO emission in bar-driven inflows.

preprint2020arXiv

Cool Outflows in Galaxies and their Implications

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.

preprint2020arXiv

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

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

preprint2020arXiv

Evolution of Stellar Feedback in HII Regions

Stellar feedback is needed to produce realistic giant molecular clouds (GMCs) and galaxies in simulations, but due to limited numerical resolution, feedback must be implemented using subgrid models. Observational work is an important means to test and anchor these models, but limited studies have assessed the relative dynamical role of multiple feedback modes, particularly at the earliest stages of expansion when HII regions are still deeply embedded. In this paper, we use multiwavelength (radio, infrared, and X-ray) data to measure the pressures associated with direct radiation ($P_{\rm dir}$), dust-processed radiation ($P_{\rm IR}$), photoionization heating ($P_{\rm HII}$), and shock-heating from stellar winds ($P_{\rm X}$) in a sample of 106 young, resolved HII regions with radii $\lesssim$0.5 pc to determine how stellar feedback drives their expansion. We find that the $P_{\rm IR}$ dominates in 84% of the regions and that the median $P_{\rm dir}$ and $P_{\rm HII}$ are smaller than the median $P_{\rm IR}$ by factors of $\approx 6$ and $\approx 9$, respectively. Based on the radial dependences of the pressure terms, we show that HII regions transition from $P_{\rm IR}$-dominated to $P_{\rm HII}$-dominated at radii of $\sim$3 pc. We find a median trapping factor of $f_{\rm trap} \sim$ 8 without any radial dependence for the sample, suggesting this value can be adopted in sub-grid feedback models. Moreover, we show that the total pressure is greater than the gravitational pressure in the majority of our sample, indicating that the feedback is sufficient to expel gas from the regions.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2020arXiv

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

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

preprint2019arXiv

The EDGE-CALIFA Survey: Using Optical Extinction to Probe the Spatially-Resolved Distribution of Gas in Nearby Galaxies

We present an empirical relation between the cold gas surface density ($Σ_{\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 $Σ_{\rm gas}$-${\rm A_V}$ relation on global, radial, and kpc (spaxel) scales. We determine ${\rm A_V}$ from the Balmer decrement (H$α$/H$β$). We find that the best fit for this relation is $Σ_{\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 $Σ_{\rm gas}$/ ${\rm A_V}$ ratio on radial and spaxel scales as a function of ${\rm EW(Hα)}$. We find that at larger values of ${\rm EW(Hα)}$ (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 $Σ_{\rm gas}$/${\rm A_V}$ ratio and the ionized gas metallicity. We contrast our estimates of $Σ_{\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 $Σ_{\rm gas}$-$Σ_*$ relation. These tests show that optical extinction is a reliable proxy for estimating $Σ_{\rm gas}$ in the absence of direct sub/millimeter observations of the cold gas.

preprint2019arXiv

Using [CII] 158 micron Emission From Isolated ISM Phases as a Star-Formation Rate Indicator

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.

preprint2018arXiv

Cloud-Scale Molecular Gas Properties in 15 Nearby Galaxies

We measure the velocity dispersion, $σ$, and surface density, $Σ$, of the molecular gas in nearby galaxies from CO spectral line cubes with spatial resolution $45$-$120$ pc, matched to the size of individual giant molecular clouds. Combining $11$ galaxies from the PHANGS-ALMA survey with $4$ targets from the literature, we characterize ${\sim}30,000$ independent sightlines where CO is detected at good significance. $Σ$ and $σ$ show a strong positive correlation, with the best-fit power law slope close to the expected value for resolved, self-gravitating clouds. This indicates only weak variation in the virial parameter $α_\mathrm{vir}\proptoσ^2/Σ$, which is ${\sim}1.5$-$3.0$ for most galaxies. We do, however, observe enormous variation in the internal turbulent pressure $P_\mathrm{turb}\proptoΣ\,σ^2$, which spans ${\sim}5\rm\;dex$ across our sample. We find $Σ$, $σ$, and $P_\mathrm{turb}$ to be systematically larger in more massive galaxies. The same quantities appear enhanced in the central kpc of strongly barred galaxies relative to their disks. Based on sensitive maps of M31 and M33, the slope of the $σ$-$Σ$ relation flattens at $Σ\lesssim10\rm\;M_\odot\,pc^{-2}$, leading to high $σ$ for a given $Σ$ and high apparent $α_\mathrm{vir}$. This echoes results found in the Milky Way, and likely originates from a combination of lower beam filling factors and a stronger influence of local environment on the dynamical state of molecular gas in the low density regime.

preprint2018arXiv

The Interstellar Medium and Star Formation in Dwarf Galaxies

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.

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2016arXiv

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

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

preprint2015arXiv

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

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

preprint2015arXiv

CARMA CO Observations of Three Extremely Metal-Poor, Star-Forming Galaxies

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.

preprint2015arXiv

Spatially extended and high-velocity dispersion molecular component in spiral galaxies: single-dish vs. interferometric observations

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.

preprint2015arXiv

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

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

preprint2015arXiv

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

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

preprint2014arXiv

ALFALFA Discovery of the Nearby Gas-rich Dwarf Galaxy Leo P. V. Neutral Gas Dynamics and Kinematics

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.

preprint2014arXiv

ALMA Reveals the Molecular Medium Fueling the Nearest Nuclear Starburst

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

preprint2014arXiv

Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. II. Gas-to-Dust Ratio Variations across ISM Phases

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.

preprint2014arXiv

Extended HCN and HCO$^{+}$ emission in the starburst galaxy M82

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.

preprint2014arXiv

Search for [CII] emission in z=6.5-11 star-forming galaxies

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(μ/15)^{-1}{\rm L}_{\odot}$ and an upper limit in the FIR luminosity of $<6.1\times10^{10}\times(μ/15)^{-1}{\rm L}_{\odot}$ (where $μ$ 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.

preprint2014arXiv

The CARMA Paired Antenna Calibration System: Atmospheric Phase Correction for Millimeter Wave Interferometry and its Application to Mapping the Ultraluminous Galaxy Arp 193

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.

preprint2014arXiv

The Role of Stellar Feedback in the Dynamics of HII Regions

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.

preprint2013arXiv

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

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

preprint2013arXiv

Dust May Be More Rare Than Expected in Metal Poor Galaxies

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

preprint2013arXiv

Dusty OB Stars in the Small Magellanic Cloud - I: Optical Spectroscopy Reveals Predominantly Main-Sequence OB Stars

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.

preprint2013arXiv

Dusty OB stars in the Small Magellanic Cloud - II: Extragalactic Disks or Examples of the Pleiades Phenomenon?

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)

preprint2013arXiv

Formation of Molecular Clouds and Global Conditions for Star Formation

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.

preprint2013arXiv

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

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

preprint2013arXiv

Millimeter dust emission compared with other mass estimates in N11 molecular clouds in the LMC

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.

preprint2013arXiv

Misalignment of Magnetic Fields and Outflows in Protostellar Cores

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.

preprint2013arXiv

The CO-to-H2 Conversion Factor

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

preprint2013arXiv

The EGNoG Survey: Gas Excitation in Normal Galaxies at z~0.3

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.

preprint2013arXiv

The EGNoG Survey: Molecular Gas in Intermediate-Redshift Star-Forming Galaxies

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.

preprint2013arXiv

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

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

preprint2012arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

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

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

preprint2011arXiv

What Drives the Expansion of Giant HII Regions?: A Study of Stellar Feedback in 30 Doradus

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.

preprint2010arXiv

Atmospheric Phase Correction using CARMA-PACS: High Angular Resolution Observations of the FU-Orionis star PP 13S*

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.

preprint2010arXiv

LABOCA observations of giant molecular clouds in the south west region of the Small Magellanic Cloud

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.

preprint2010arXiv

Spitzer Analysis of HII Region Complexes in the Magellanic Clouds: Determining a Suitable Monochromatic Obscured Star Formation Indicator

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.

preprint2010arXiv

The Spitzer Survey of the Small Magellanic Cloud (S3MC): Insights into the Life-Cycle of Polycyclic Aromatic Hydrocarbons

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.

preprint2009arXiv

Infrared Dark Clouds in the Small Magellanic Cloud?

We have applied the unsharp-masking technique to the 24 $μ$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 $μ$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 $μ$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.

preprint2008arXiv

Measuring Dust Production in the Small Magellanic Cloud Core-Collapse Supernova Remnant 1E 0102.2-7219

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.

preprint2006arXiv

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

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