Molecular depletion times and the CO-to-H2 conversion factor in metal-poor galaxies

Tracing molecular hydrogen content with carbon monoxide in low-metallicity galaxies has been exceedingly difficult. Here we present a new effort, with IRAM 30-m observations of 12CO(1-0) of a sample of 8 dwarf galaxies having oxygen abundances ranging from 12+logO/H=7.7 to 8.4. CO emission is detected in all galaxies, including the most metal-poor galaxy of our sample (0.1 Zsun); to our knowledge this is the largest number of 12CO(1-0) detections ever reported for galaxies with 12+logO/H<=8 (0.2 Zsun) outside the Local Group. We calculate stellar masses (Mstar) and star-formation rates (SFRs), and analyze our results by combining our observations with galaxy samples from the literature. Extending previous results for a correlation of the molecular gas depletion time, tau(dep), with Mstar and specific SFR (sSFR), we find a variation in tau(dep) of a factor of 200 or more (from <50 Myr to 10 Gyr) over a spread of 1000 in sSFR and Mstar. We exploit the variation of tau(dep) to constrain the CO-to-H2 mass conversion factor alpha(CO) at low metallicity, and assuming a power-law variation find alpha(CO) \propto (Z/Zsun)^1.9, similar to results based on dust continuum measurements compared with gas mass. By including HI measurements, we show that the fraction of total gas mass relative to the baryonic mass is higher in galaxies that are metal poor, of low mass, and of high sSFR. Finally, comparisons of the data with star-formation models of the molecular gas phases suggest that, at metallicities Z/Zsun<=0.2, there are some discrepancies with model predictions.Comment: 18 pages, 15 figures, accepted for publication in A&

ALMA observations of cool dust in a low-metallicity starburst, SBS0335-052

We present Atacama Large Millimeter/submillimeter Array (ALMA) Cycle 0 Band 7 observations of an extremely metal-poor dwarf starburst galaxy in the Local Universe, SBS0335-052 (12+log(O/H)~7.2). With these observations, dust is detected at 870micron (ALMA Band 7), but 87% of the flux in this band is due to free-free emission from the starburst. We have compiled a spectral energy distribution (SED) of SBS0335-052 that spans almost 6 orders of magnitude in wavelength and fit it with a spherical dust shell heated by a single-age stellar population; the best-fit model gives a dust mass of (3.8+/-0.6)x10^4 Msun. We have also constructed a SED including Herschel archival data for IZw18, another low-metallicity dwarf starburst (12+log(O/H)=7.17), and fit it with a similar model to obtain a dust mass of (3.4+/-1.0)x10^2 Msun. Compared with their atomic gas mass, the dust mass of SBS0335-052 far exceeds the prediction of a linear trend of dust-to-gas mass ratio with metallicity, while IZw18 falls far below. We use gas scaling relations to assess a putative missing gas component in both galaxies and find that the missing, possibly molecular, gas in SBS0335-052 is a factor of 6 times higher than the value inferred from the observed HI column density; in IZw18 the missing component is 4 times smaller. Ultimately, despite their similarly low metallicity, the differences in gas and dust column densities in SBS0335-052 and IZw18 suggest that metal abundance does not uniquely define star-formation processes. At some level, self-shielding and the survival of molecules may depend just as much on gas and dust column density as on metallicity. The effects of low metallicity may at least be partially compensated for by large column densities in the interstellar medium.Comment: 15 pages, 11 figures, accepted for publication in A&

ALMA reveals a warm and compact starburst around a heavily obscured supermassive black hole at z=4.75

We report ALMA Cycle 0 observations at 1.3mm of LESS J033229.4-275619 (XID403), an Ultraluminous Infrared Galaxy at $z=4.75$ in the Chandra Deep Field South hosting a Compton-thick QSO. The source is not resolved in our data at a resolution of $\sim$0.75 arcsec, placing an upper-limit of 2.5 kpc to the half-light radius of the continuum emission from heated-dust. After deconvolving for the beam size, however, we found a $\sim3\sigma$ indication of an intrinsic source size of $0.27\pm0.08$ arcsec (Gaussian FWHM), which would correspond to $r_{half}\sim0.9\pm0.3$ kpc. We build the far-IR SED of XID403 by combining datapoints from both ALMA and Herschel and fit it with a modified blackbody spectrum. For the first time, we measure the dust temperature $T_d=58.5\pm5.3$ K in this system, which is comparable to what has been observed in other high-z submillimeter galaxies. The measured star formation rate is SFR=$1020\pm150$ $M_{\odot}$ yr$^{-1}$, in agreement with previous estimates at lower S/N. Based on the measured SFR and source size, we constrain the SFR surface density to be $\Sigma_{SFR}>26\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ ($\sim200\;M_{\odot}$yr$^{-1}$kpc$^{-2}$ for $r_{half}\sim0.9$ kpc). The compactness of this starburst is comparable to what has been observed in other local and high-z starburst galaxies. If the gas mass measured from previous [CII] and CO(2-1) observations at low resolution is confined within the same dust region, assuming $r_{half}\sim0.9\pm0.3$ kpc, this would produce a column density of $N_H\sim0.3-1.1\times10^{24}$cm$^{-2}$ towards the central SMBH, similar to the column density of $\approx1.4\times10^{24}$cm$^{-2}$ measured from the X-rays. Then, in principle, if both gas and dust were confined on sub-kpc scales, this would be sufficient to produce the observed X-ray column density without any need of a pc-scale absorber [abridged].Comment: 11 pages, 5 figures. Accepted for publication in A&

AGN impact on the molecular gas in galactic centres as probed by CO lines

We present a detailed analysis of the X-ray, infrared, and carbon monoxide (CO) emission for a sample of 35 local (z ≤ 0.15), active (LX ≥ 1042 erg s-1) galaxies. Our goal is to infer the contribution of far-ultraviolet (FUV) radiation from star formation (SF), and X-ray radiation from the active galactic nuclei (AGNs), respectively, producing photodissociation regions (PDRs) and X-ray-dominated regions (XDRs), to the molecular gas heating. To this aim, we exploit the CO spectral line energy distribution (CO SLED) as traced by Herschel, complemented with data from single-dish telescopes for the low-J lines, and high-resolution ALMA images of the mid-J CO emitting region. By comparing our results to the Schmidt-Kennicutt relation, we find no evidence for AGN influence on the cold and low-density gas on kpc-scales. On nuclear (r = 250 pc) scales, we find weak correlations between the CO line ratios and either the FUV or X-ray fluxes: this may indicate that neither SF nor AGN radiation dominates the gas excitation, at least at r = 250 pc. From a comparison of the CO line ratios with PDR and XDR models, we find that PDRs can reproduce observations only in presence of extremely high gas densities (n &gt; 105 cm-3). In the XDR case, instead, the models suggest moderate densities (n ≈ 102-4 cm-3). We conclude that a mix of the two mechanisms (PDR for the mid-J, XDR, or possibly shocks for the high-J) is necessary to explain the observed CO excitation in active galaxies

ALMA resolves the torus of NGC 1068: continuum and molecular line emission

We have used the Atacama Large Millimeter Array (ALMA) to map the emission of the CO(6-5) molecular line and the 432 {\mu}m continuum emission from the 300 pc-sized circumnuclear disk (CND) of the nearby Seyfert 2 galaxy NGC 1068 with a spatial resolution of ~4 pc. These observations spatially resolve the CND and, for the first time, image the dust emission, the molecular gas distribution, and the kinematics from a 7-10 pc-diameter disk that represents the submillimeter counterpart of the putative torus of NGC 1068. We fitted the nuclear spectral energy distribution of the torus using ALMA and near and mid-infrared (NIR/MIR) data with CLUMPY models. The mass and radius of the best-fit solution for the torus are both consistent with the values derived from the ALMA data alone: Mgas_torus=(1+-0.3)x10^5 Msun and Rtorus=3.5+-0.5 pc. The dynamics of the molecular gas in the torus show non-circular motions and enhanced turbulence superposed on the rotating pattern of the disk. The kinematic major axis of the CO torus is tilted relative to its morphological major axis. By contrast with the nearly edge-on orientation of the H2O megamaser disk, we have found evidence suggesting that the molecular torus is less inclined (i=34deg-66deg) at larger radii. The lopsided morphology and complex kinematics of the torus could be the signature of the Papaloizou-Pringle instability, long predicted to likely drive the dynamical evolution of active galactic nuclei (AGN) tori.Comment: Final version accepted by the Astrophysical Journal Letters (ApJLetters) on April 27th 2016, 6 pages, 5 figure

Italian Science Case for ALMA Band 2+3

The Premiale Project "Science and Technology in Italy for the upgraded ALMA Observatory - iALMA" has the goal of strengthening the scientific, technological and industrial Italian contribution to the Atacama Large Millimeter/submillimeter Array (ALMA), the largest ground based international infrastructure for the study of the Universe in the microwave. One of the main objectives of the Science Working Group (SWG) inside iALMA, the Work Package 1, is to develop the Italian contribution to the Science Case for the ALMA Band 2 or Band 2+3 receiver. ALMA Band 2 receiver spans from ~67 GHz (bounded by an opaque line complex of ozone lines) up to 90 GHz which overlaps with the lower frequency end of ALMA Band 3. Receiver technology has advanced since the original definition of the ALMA frequency bands. It is now feasible to produce a single receiver which could cover the whole frequency range from 67 GHz to 116 GHz, encompassing Band 2 and Band 3 in a single receiver cartridge, a so called Band 2+3 system. In addition, upgrades of the ALMA system are now foreseen that should double the bandwidth to 16 GHz. The science drivers discussed below therefore also discuss the advantages of these two enhancements over the originally foreseen Band 2 system.Comment: 43 pages, 21 figure

Molecular line emission in NGC1068 imaged with ALMA. I An AGN-driven outflow in the dense molecular gas

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.Comment: Version accepted for publication in A&A (June 4th). Accepted version. References (3) added and minor typos corrected. 24 pages, 20 figure

The volume densities of giant molecular clouds in M83

Using observed GALEX far-ultraviolet (FUV) fluxes and VLA images of the 21-cm HI column densities, along with estimates of the local dust abundances, we measure the volume densities of a sample of actively star-forming giant molecular clouds (GMCs) in the nearby spiral galaxy M83 on a typical resolution scale of 170 pc. Our approach is based on an equilibrium model for the cycle of molecular hydrogen formation on dust grains and photodissociation under the influence of the FUV radiation on the cloud surfaces of GMCs. We find a range of total volume densities on the surface of GMCs in M83, namely 0.1 - 400 cm-3 inside R25, 0.5 - 50 cm-3 outside R25 . Our data include a number of GMCs in the HI ring surrounding this galaxy. Finally, we discuss the effects of observational selection, which may bias our results.Comment: 9 pages, 11 figure

Molecular line emission in NGC1068 imaged with ALMA: II. The chemistry of the dense molecular gas

We present a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.Comment: Accepted by A&A; please contact the author for a better version of the pdf where the resolution and positions of figures are as they will appear in the Journa