Showing 1–50 of 110 results for author: Israel, F P

The Physical Drivers and Observational Tracers of CO-to-H2 Conversion Factor Variations in Nearby Barred Galaxy Centers

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Munan Gong, Alberto D. Bolatto, I-Da Chiang, Adam K. Leroy, Antonio Usero, Simon C. O. Glover, Ralf S. Klessen, Daizhong Liu, Miguel Querejeta, Eva Schinnerer, Frank Bigiel, Yixian Cao, Melanie Chevance, Cosima Eibensteiner, Kathryn Grasha, Frank P. Israel, Eric J. Murphy, Lukas Neumann, Hsi-An Pan, Francesca Pinna, Mattia C. Sormani, J. D. T. Smith , et al. (2 additional authors not shown)

Abstract: The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $α_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NG… ▽ More The CO-to-H$_2$ conversion factor ($α_\rm{CO}$) is central to measuring the amount and properties of molecular gas. It is known to vary with environmental conditions, and previous studies have revealed lower $α_\rm{CO}$ in the centers of some barred galaxies on kpc scales. To unveil the physical drivers of such variations, we obtained ALMA Band 3, 6, and 7 observations toward the inner 2 kpc of NGC 3627 and NGC 4321 tracing $^{12}$CO, $^{13}$CO, and C$^{18}$O lines on 100 pc scales. Our multi-line modeling and Bayesian likelihood analysis of these datasets reveal variations of molecular gas density, temperature, optical depth, and velocity dispersion, which are among the key drivers of $α_\rm{CO}$. The central 300 pc nuclei in both galaxies show strong enhancement of temperature $T_\rm{k}>100$ K and density $n_\rm{H_2}>10^3$ cm$^{-3}$. Assuming a CO-to-H$_2$ abundance of $3\times10^{-4}$, we derive 4-15 times lower $α_\rm{CO}$ than the Galactic value across our maps, which agrees well with previous kpc-scale measurements. Combining the results with our previous work on NGC 3351, we find a strong correlation of $α_\rm{CO}$ with low-J $^{12}$CO optical depths ($τ_\rm{CO}$), as well as an anti-correlation with $T_\rm{k}$. The $τ_\rm{CO}$ correlation explains most of the $α_\rm{CO}$ variation in the three galaxy centers, whereas changes in $T_\rm{k}$ influence $α_\rm{CO}$ to second order. Overall, the observed line width and $^{12}$CO/$^{13}$CO 2-1 line ratio correlate with $τ_\rm{CO}$ variation in these centers, and thus they are useful observational indicators for $α_\rm{CO}$ variation. We also test current simulation-based $α_\rm{CO}$ prescriptions and find a systematic overprediction, which likely originates from the mismatch of gas conditions between our data and the simulations. △ Less

Submitted 10 April, 2023; originally announced April 2023.

Comments: Accepted for publication in ApJ; 30 pages of main text + 3 appendices

Central molecular zones in galaxies: multitransition survey of dense gas tracers HCN, HNC, and HCO+

Authors: F. P. Israel

Abstract: New HCN, HNC, and HCO+ measurements of 46 normal galaxies in transitions up to J=4-3 are included in a multitransition database covering HCN and HCO+ (130 galaxies) and HNC (94 galaxies). The near-linear luminosity relations are dominated by distance effects and do not reflect galaxy physical properties. Individual galaxies show significant dispersion in both their luminosity and line ratio. Line… ▽ More New HCN, HNC, and HCO+ measurements of 46 normal galaxies in transitions up to J=4-3 are included in a multitransition database covering HCN and HCO+ (130 galaxies) and HNC (94 galaxies). The near-linear luminosity relations are dominated by distance effects and do not reflect galaxy physical properties. Individual galaxies show significant dispersion in both their luminosity and line ratio. Line ratios do not correlate well with either luminosities or other line ratios. Only the normalized J-transition ladders of HCN and HCO+ and the J=1-0 12CO/13CO isotopologue ratio are positively correlated with CO and far infrared (FIR) luminosity. The HCN and HCO+ molecules have very similar intensities and trace the same gas. In galaxies dominated by an active nucleus, HCO+ intensities are depressed relative to HCN intensities. Only a small fraction of the CO emission is directly associated with gas emitting in HCN and HCO, yet a significant fraction of even that gas appears to be translucent molecular gas. In the observed galaxy centers, the HCN/CO line intensity ratio is not a proxy for the dense gas fraction. Likewise, the FIR/HCN and FIR/CO ratios are not proxies for the star formation efficiency. The observed molecular line emission is fully consistent with UV-photon heating boosted by significant mechanical heating. The molecular gas sampled by HCN and HCO+ has low kinetic temperatures T(kin)=10-50 K, low densities n(H)=10^4-10^5 cm^(-3), and low optical depths in the ground-state lines. Most of the gas sampled by CO has densities lower by one to two orders of magnitude. For a mechanical heating fraction of 0.5, a modest energy input of only G=300 Go is required. A proper understanding of star formation requires a more appropriate determination of the gas mass than provided by the intensities of individual HCN or CO transitions. △ Less

Submitted 1 December, 2022; originally announced December 2022.

Comments: 19 pages, 11 figures, 8 tables, accepted by Astronomy & Astrophysics

Journal ref: A&A 671, A59 (2023)

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

Authors: Yu-Hsuan Teng, Karin M. Sandstrom, Jiayi Sun, Adam K. Leroy, L. Clifton Johnson, Alberto D. Bolatto, J. M. Diederik Kruijssen, Andreas Schruba, Antonio Usero, Ashley T. Barnes, Frank Bigiel, Guillermo A. Blanc, Brent Groves, Frank P. Israel, Daizhong Liu, Erik Rosolowsky, Eva Schinnerer, J. D. Smith, Fabian Walter

Abstract: 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 sp… ▽ More 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. △ Less

Submitted 28 December, 2021; v1 submitted 10 November, 2021; originally announced November 2021.

Comments: 31 pages, 17 figures; accepted for publication in ApJ; font size of text and figures adjusted

The peculiar nebula Simeis 57: II. Distance, nature and excitation

Authors: L. H. T. Oudshoorn, F. P. Israel, J. Brinchmann, M. B. C. Kloppenburg, A. G. A. Brown, J. Bally, T. R. Gull, P. T. Boyd

Abstract: Simeis 57 (HS 191) is an optically bright nebula in the Cygnus X region with a peculiar appearance that suggests an outflow from a rotating source. Newly obtained observations and archival data reveal Simeis 57 as a low-density ($n_{e}\,\sim\,100$ cm$^{-3}$) nebula with an east-to-west excitation gradient. The extinction of the nebula is $A_{V}\,\leq$ 2 mag. The nebula is recognizable but not prom… ▽ More Simeis 57 (HS 191) is an optically bright nebula in the Cygnus X region with a peculiar appearance that suggests an outflow from a rotating source. Newly obtained observations and archival data reveal Simeis 57 as a low-density ($n_{e}\,\sim\,100$ cm$^{-3}$) nebula with an east-to-west excitation gradient. The extinction of the nebula is $A_{V}\,\leq$ 2 mag. The nebula is recognizable but not prominent in mid- and far-infrared images. In its direction, half a dozen small CO clouds have been identified at $V_{LSR}$ = + 5 km s$^{-1}$. One of these coincides with both the optical nebula and a second CO cloud at the nebular velocity $V_{LSR}\,\approx$ -10 km $^{-1}$. No luminous stars are embedded in these molecular clouds, nor are any obscured by them and no sufficiently luminous stars are found in the immediate vicinity of the nebula. Instead, all available data points to the evolved star HD 193793 = WR 140 (an O4-5 supergiant and WC7 Wolf-Rayet binary) as the source of excitation, notwithstanding its large separation of $50'$, about 25 pc at the stellar distance of 1.7 kpc. Simeis 57 appears to be a part of a larger structure surrounding the HI void centered on HD 193793. △ Less

Submitted 29 April, 2021; originally announced April 2021.

A Giant Molecular Cloud Catalog in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

Authors: E. R. Miura, D. Espada, A. Hirota, C. Henkel, S. Verley, M. Kobayashi, S. Matsushita, F. P. Israel, B. Vila-Vilaro, K. Morokuma-Matsui, J. Ott, C. Vlahakis, A. B. Peck, S. Aalto, M. Hogerheijde, N. Neumayer, D. Iono, K. Kohno, H. Takemura, S. Komugi

Abstract: We present the first census of giant molecular clouds (GMCs) complete down to 10$^6 M_{\odot}$ and within the inner 4 kpc of the nearest giant elliptical and powerful radio galaxy, Centaurus A. We identified 689 GMCs using CO(1--0) data with 1" spatial resolution ($\sim 20$ pc) and 2 km/s velocity resolution obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). The $I$(CO)-… ▽ More We present the first census of giant molecular clouds (GMCs) complete down to 10$^6 M_{\odot}$ and within the inner 4 kpc of the nearest giant elliptical and powerful radio galaxy, Centaurus A. We identified 689 GMCs using CO(1--0) data with 1" spatial resolution ($\sim 20$ pc) and 2 km/s velocity resolution obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). The $I$(CO)-$N$(H$_2$) conversion factor based on the virial method is $X_{\rm CO}$ = $(2 \pm 1 )\times10^{20}$ cm$^{-2}$(K km/s)$^{-1}$ for the entire molecular disk, consistent with that of the disks of spiral galaxies including the Milky Way, and $X_{\rm CO}$ = $(5 \pm 2)\times10^{20}$ cm$^{-2}$(K km/s)$^{-1}$ for the circumnuclear disk (CND, within a galactocentric radius of 200 pc). We obtained the GMC mass spectrum distribution and find that the best-truncated power-law fit for the whole molecular disk, with index $γ\simeq -2.41 \pm 0.02$ and upper cutoff mass $\sim 1.3 \times 10^{7} M_{\odot}$, is also in agreement with that of nearby disk galaxies. A trend is found in the mass spectrum index from steep to shallow as we move to inner radii. Although the GMCs are in an elliptical galaxy, the general GMC properties in the molecular disk are as in spiral galaxies. However, in the CND, large offsets in the line-width-size scaling relations ($\sim$ 0.3 dex higher than those in the GMCs in the molecular disk), a different $X_{\rm CO}$ factor, and the shallowest GMC mass distribution shape ($γ= -1.1 \pm 0.2$) all suggest that there the GMCs are most strongly affected by the presence of the AGN and/or shear motions. △ Less

Submitted 26 April, 2021; originally announced April 2021.

Comments: 19 pages, 11 figures, 3 tables, accepted for publication in MNRAS

Central molecular zones in galaxies: 12CO-to-13CO ratios, carbon budget, and X factors

Authors: F. P. Israel

Abstract: This paper presents ground-based 12CO and 13CO measurements of 126 nearby galaxy centers in various J transitions. More than 60 galaxies were measured in at least four lines. The average relative intensities of the first four 12CO transitions are 1.00 : 0.92 : 0.70 : 0.57. The average 12CO-to-13CO ratios are 13.0, 11.6, and 12.8 for the first three transitiions. The sizes of central CO concentrati… ▽ More This paper presents ground-based 12CO and 13CO measurements of 126 nearby galaxy centers in various J transitions. More than 60 galaxies were measured in at least four lines. The average relative intensities of the first four 12CO transitions are 1.00 : 0.92 : 0.70 : 0.57. The average 12CO-to-13CO ratios are 13.0, 11.6, and 12.8 for the first three transitiions. The sizes of central CO concentrations are well defined in maps, but poorly determined by multi-aperture photometry. Using radiative transfer models (RADEX), we derived model gas parameters, where the assumed carbon elemental abundances and carbon gas depletion onto dust are the main causes of uncertainty. The new CO data and published [CI] and [CII] data imply that CO, Co, and C+ each represent about one-third of the gas-phase carbon in the molecular interstellar medium in galaxy centers Their mean beam-averaged central molecular hydrogen column density is N(H2) = (1.5+/-0.2) x 10^(21) cm^(2) and CO-to-H2 conversion factors are typically ten times lower than the `standard' Milky Way disk value, with a mean X(CO) = (1.9+/-0.2) x 10^(19) cm^(2)/K km/s. The corresponding [CI]-to-H2 factor is five times higher than X(CO), with X[CI] = (9+/-2) x 10^(19) cm^(2)/K km/s. No unique conversion factor can be determined for [CII]. The low molecular gas content of galaxy centers relative to their CO intensities is explained in roughly equal parts by higher central gas-phase carbon abundances, elevated gas temperatures, and larger gas velocity dispersions than found in galaxy disks. △ Less

Submitted 14 January, 2020; originally announced January 2020.

Comments: Main paper 23 pages including 17 figures, appendices 18 pages

Journal ref: A&A 635, A131 (2020)

Star Formation Efficiencies at Giant Molecular Cloud Scales in the Molecular Disk of the Elliptical Galaxy NGC 5128 (Centaurus A)

Authors: D. Espada, S. Verley, R. E. Miura, F. P. Israel, C. Henkel, S. Matsushita, B. Vila-Vilaro, J. Ott, K. Morokuma-Matsui, A. B. Peck, A. Hirota, S. Aalto, A. C. Quillen, M. R. Hogerheijde, N. Neumayer, C. Vlahakis, D. Iono, K. Kohno

Abstract: We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previ… ▽ More We present ALMA CO(1-0) observations toward the dust lane of the nearest elliptical and radio galaxy, NGC 5128 (Centaurus A), with high angular resolution ($\sim$ 1 arcsec, or 18 pc), including information from large to small spatial scales and total flux. We find a total molecular gas mass of 1.6$\times$10$^9$ $M_\odot$ and we reveal the presence of filamentary components more extended than previously seen, up to a radius of 4 kpc. We find that the global star formation rate is $\sim$1 \Msol yr$^{-1}$, which yields a star formation efficiency (SFE) of 0.6 Gyr$^{-1}$ (depletion time $τ=$1.5 Gyr), similar to those in disk galaxies. We show the most detailed view to date (40\,pc resolution) of the relation between molecular gas and star formation within the stellar component of an elliptical galaxy, from several kpc scale to the circumnuclear region close to the powerful radio jet. Although on average the SFEs are similar to those of spiral galaxies, the circumnuclear disk (CND) presents SFEs of 0.3 Gyr$^{-1}$, lower by a factor of 4 than the outer disk. The low SFE in the CND is in contrast to the high SFEs found in the literature for the circumnuclear regions of some nearby disk galaxies with nuclear activity, probably as a result of larger shear motions and longer AGN feedback. The higher SFEs in the outer disk suggests that only central molecular gas or filaments with sufficient density and strong shear motions will remain in $\sim$1 Gyr, which will later result in the compact molecular distributions and low SFEs usually seen in other giant ellipticals with cold gas. △ Less

Submitted 16 December, 2019; v1 submitted 4 June, 2019; originally announced June 2019.

Comments: Minor update. 20 pages, 11 figures, accepted for publication in ApJ

Herschel water maps towards the vicinity of the black hole Sgr A*

Authors: J. Armijos-Abendaño, J. Martín-Pintado, M. A. Requena-Torres, E. González-Alfonso, R. Güsten, A. Weiß, A. I. Harris, F. P. Israel, C. Kramer, J. Stutzki, P. van der Werf

Abstract: Aims: We study the spatial distribution and kinematics of water emission in a ~64 pc$^2$ region of the Galactic Center (GC) around Sgr A*. We also analyze the water excitation to derive the physical conditions and water abundances in the CND and the `quiescent clouds'. Methods: We presented the integrated intensity maps of the ortho 1$_{10}-1_{01}$, and para 2$_{02}-1_{11}$ and 1$_{11}-0_{00}$ wat… ▽ More Aims: We study the spatial distribution and kinematics of water emission in a ~64 pc$^2$ region of the Galactic Center (GC) around Sgr A*. We also analyze the water excitation to derive the physical conditions and water abundances in the CND and the `quiescent clouds'. Methods: We presented the integrated intensity maps of the ortho 1$_{10}-1_{01}$, and para 2$_{02}-1_{11}$ and 1$_{11}-0_{00}$ water transitions observed with the HIFI instrument on board Herschel. To study the water excitation we used ground state ortho and para H$_2^{18}$O transitions. In our study, we also used SPIRE continuum measurements of the CND. Using a non-LTE radiative transfer code, the water line profiles and dust continuum were modeled. We also used a rotating ring model to reproduce the CND kinematics represented by the PV diagram. Results: We identify the water emission arising from the CND, the Western Streamer, and the 20 and 50 km s$^{-1}$ clouds. The ortho water maps show absorption structures in the range of [-220,10] km s$^{-1}$. The PV diagram shows that the 2$_{02}-1_{11}$ H$_2$O emission traces the CND. We derive high X$_{H_2O}$ of $\sim$(0.1-1.3)$\times$10$^{-5}$, V$_t$ of 14-23 km s$^{-1}$ and T$_d$ of 15-45 K for the CND, and the lower X$_{\rm H_2O}$ of 4$\times$10$^{-8}$ and V$_t$ of 9 km s$^{-1}$ for the 20 km s$^{-1}$ cloud. Collisional excitation and dust effects are responsible for the water excitation in the southwest lobe of the CND and the 20 km s$^{-1}$ cloud, whereas only collisions can account for the water excitation in the northeast lobe of the CND. We propose that the water vapor in the CND is caused by grain sputtering by shocks of 10-20 km s$^{-1}$, with some contribution of high temperature and cosmic-ray chemistries plus a PDR chemistry. The low X$_{\rm H_2O}$ derived for the 20 km s$^{-1}$ cloud could be partially a consequence of the water freeze-out on grains. △ Less

Submitted 13 February, 2019; originally announced February 2019.

Comments: 15 pages, 11 figures

Journal ref: A&A 624, A112 (2019)

Planck Observations of M33

Authors: C. T. Tibbs, F. P. Israel, R. J. Laureijs, J. A. Tauber, B. Partridge, M. W. Peel, L. Fauvet

Abstract: We have performed a comprehensive investigation of the global integrated flux density of M33 from radio to ultraviolet wavelengths, finding that the data between $\sim$100 GHz and 3 THz are accurately described by a single modified blackbody curve with a dust temperature of $T_\mathrm{dust}$ = 21.67$\pm$0.30 K and an effective dust emissivity index of $β_\mathrm{eff}$ = 1.35$\pm$0.10, with no indi… ▽ More We have performed a comprehensive investigation of the global integrated flux density of M33 from radio to ultraviolet wavelengths, finding that the data between $\sim$100 GHz and 3 THz are accurately described by a single modified blackbody curve with a dust temperature of $T_\mathrm{dust}$ = 21.67$\pm$0.30 K and an effective dust emissivity index of $β_\mathrm{eff}$ = 1.35$\pm$0.10, with no indication of an excess of emission at millimeter/sub-millimeter wavelengths. However, sub-dividing M33 into three radial annuli, we found that the global emission curve is highly degenerate with the constituent curves representing the sub-regions of M33. We also found gradients in $T_\mathrm{dust}$ and $β_\mathrm{eff}$ across the disk of M33, with both quantities decreasing with increasing radius. Comparing the M33 dust emissivity with that of other Local Group members, we find that M33 resembles the Magellanic Clouds rather than the larger galaxies, i.e., the Milky Way and M31. In the Local Group sample, we find a clear correlation between global dust emissivity and metallicity, with dust emissivity increasing with metallicity. A major aspect of this analysis is the investigation into the impact of fluctuations in the Cosmic Microwave Background (CMB) on the integrated flux density spectrum of M33. We found that failing to account for these CMB fluctuations would result in a significant over-estimate of $T_\mathrm{dust}$ by $\sim$5 K and an under-estimate of $β_\mathrm{eff}$ by $\sim$0.4. △ Less

Submitted 27 March, 2018; originally announced March 2018.

Comments: Accepted for publication in MNRAS

The State-of-Play of Anomalous Microwave Emission (AME) Research

Authors: Clive Dickinson, Y. Ali-Haïmoud, A. Barr, E. S. Battistelli, A. Bell, L. Bernstein, S. Casassus, K. Cleary, B. T. Draine, R. Génova-Santos, S. E. Harper, B. Hensley, J. Hill-Valler, Thiem Hoang, F. P. Israel, L. Jew, A. Lazarian, J. P. Leahy, J. Leech, C. H. López-Caraballo, I. McDonald, E. J. Murphy, T. Onaka, R. Paladini, M. W. Peel , et al. (8 additional authors not shown)

Abstract: Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot b… ▽ More Anomalous Microwave Emission (AME) is a component of diffuse Galactic radiation observed at frequencies in the range $\approx 10$-60 GHz. AME was first detected in 1996 and recognised as an additional component of emission in 1997. Since then, AME has been observed by a range of experiments and in a variety of environments. AME is spatially correlated with far-IR thermal dust emission but cannot be explained by synchrotron or free-free emission mechanisms, and is far in excess of the emission contributed by thermal dust emission with the power-law opacity consistent with the observed emission at sub-mm wavelengths. Polarization observations have shown that AME is very weakly polarized ($\lesssim 1$%). The most natural explanation for AME is rotational emission from ultra-small dust grains ("spinning dust"), first postulated in 1957. Magnetic dipole radiation from thermal fluctuations in the magnetization of magnetic grain materials may also be contributing to the AME, particularly at higher frequencies ($\gtrsim 50$ GHz). AME is also an important foreground for Cosmic Microwave Background analyses. This paper presents a review and the current state-of-play in AME research, which was discussed in an AME workshop held at ESTEC, The Netherlands, June 2016. △ Less

Submitted 26 February, 2018; v1 submitted 22 February, 2018; originally announced February 2018.

Comments: Accepted for publication in New Astronomy Reviews. Summary of AME workshop held at ESTEC, The Netherlands, June 2016, 40 pages, 18 figures. Updated to approximately match published version

First Results from the $Herschel$ and ALMA Spectroscopic Surveys of the SMC: The Relationship Between [CII]-bright Gas and CO-bright Gas at Low Metallicity

Authors: Katherine E. Jameson, Alberto D. Bolatto, Mark Wolfire, Steven R. Warren, Rodrigo Herrera-Camus, Kevin Croxall, Eric Pellegrini, John-David Smith, Monica Rubio, Remy Indebetouw, Frank P. Israel, Margaret Meixner, Julia Roman-Duval, Jacco Th. van Loon, Erik Muller, Celia Verdugo, Hans Zinnecker, Yoko Okada

Abstract: The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverag… ▽ More The Small Magellanic Cloud (SMC) provides the only laboratory to study the structure of molecular gas at high resolution and low metallicity. We present results from the Herschel Spectroscopic Survey of the SMC (HS$^{3}$), which mapped the key far-IR cooling lines [CII], [OI], [NII], and [OIII] in five star-forming regions, and new ALMA 7m-array maps of $^{12}$CO and $^{13}$CO $(2-1)$ with coverage overlapping four of the five HS$^{3}$ regions. We detect [CII] and [OI] throughout all of the regions mapped. The data allow us to compare the structure of the molecular clouds and surrounding photodissociation regions using $^{13}$CO, CO, [CII], and [OI] emission at $<10$" ($<3$ pc) scales. We estimate Av using far-IR thermal continuum emission from dust and find the CO/[CII] ratios reach the Milky Way value at high A$_{V}$ in the centers of the clouds and fall to $\sim{1/5-1/10}\times$ the Milky Way value in the outskirts, indicating the presence of translucent molecular gas not traced by bright CO emission. We estimate the amount of molecular gas traced by bright [CII] emission at low A$_{V}$ and bright CO emission at high A$_{V}$. We find that most of the molecular gas is at low A$_{V}$ and traced by bright [CII] emission, but that faint CO emission appears to extend to where we estimate the H$_{2}$-to-HI transition occurs. By converting our H$_{2}$ gas estimates to a CO-to-H$_{2}$ conversion factor ($X_{CO}$), we show that $X_{CO}$ is primarily a function of A$_{V}$, consistent with simulations and models of low metallicity molecular clouds. △ Less

Submitted 10 January, 2018; originally announced January 2018.

Comments: Accepted for publication in ApJ

ALMA Observations of the Physical and Chemical Conditions in Centaurus A

Authors: Mark McCoy, Jürgen Ott, David S. Meier, Sébastien Muller, Daniel Espada, Sergio Martín, Frank P. Israel, Christian Henkel, Violette Impellizzeri, Susanne Aalto, Philip G. Edwards, Andreas Brunthaler, Nadine Neumayer, Alison B. Peck, Paul van der Werf, Ilana Feain

Abstract: Centaurus A, with its gas-rich elliptical host galaxy, NGC 5128, is the nearest radio galaxy at a distance of 3.8 Mpc. Its proximity allows us to study the interaction between an active galactic nucleus, radio jets, and molecular gas in great detail. We present ALMA observations of low J transitions of three CO isotopologues, HCN, HCO$^{+}$, HNC, CN, and CCH toward the inner projected 500 pc of NG… ▽ More Centaurus A, with its gas-rich elliptical host galaxy, NGC 5128, is the nearest radio galaxy at a distance of 3.8 Mpc. Its proximity allows us to study the interaction between an active galactic nucleus, radio jets, and molecular gas in great detail. We present ALMA observations of low J transitions of three CO isotopologues, HCN, HCO$^{+}$, HNC, CN, and CCH toward the inner projected 500 pc of NGC 5128. Our observations resolve physical sizes down to 40 pc. By observing multiple chemical probes, we determine the physical and chemical conditions of the nuclear interstellar medium of NGC 5128. This region contains molecular arms associated with the dust lanes and a circumnuclear disk (CND) interior to the molecular arms. The CND is approximately 400 pc by 200 pc and appears to be chemically distinct from the molecular arms. It is dominated by dense gas tracers while the molecular arms are dominated by $^{12}$CO and its rare isotopologues. The CND has a higher temperature, elevated CN/HCN and HCN/HNC intensity ratios, and much weaker $^{13}$CO and C$^{18}$O emission than the molecular arms. This suggests an influence from the AGN on the CND molecular gas. There is also absorption against the AGN with a low velocity complex near the systemic velocity and a high velocity complex shifted by about 60 km s$^{-1}$. We find similar chemical properties between the CND in emission and both the low and high velocity absorption complexes implying that both likely originate from the CND. If the HV complex does originate in the CND, then that gas would correspond to gas falling toward the supermassive black hole. △ Less

Submitted 1 December, 2017; originally announced December 2017.

HIFI Spectroscopy of ${\rm H_2O}$ submm Lines in Nuclei of Actively Star Forming Galaxies

Authors: L. Liu, A. Weiss, J. P. Perez-Beaupuits, R. Güsten, D. Liu, Y. Gao, K. M. Menten, P. van der Werf, F. P. Israel, A. Harris, J. Martin-Pintado, M. A. Requena-Torres, J. Stutzki

Abstract: We present a systematic survey of multiple velocity-resolved H$_2$O spectra using Herschel/HIFI towards nine nearby actively star forming galaxies. The ground-state and low-excitation lines (E$_{\rm up}\,\le 130\,{\rm K}$) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines ($130\,{\rm K}\, \le\, E_{\rm up}\,\le\,350\,{\rm K}$) typically displ… ▽ More We present a systematic survey of multiple velocity-resolved H$_2$O spectra using Herschel/HIFI towards nine nearby actively star forming galaxies. The ground-state and low-excitation lines (E$_{\rm up}\,\le 130\,{\rm K}$) show profiles with emission and absorption blended together, while absorption-free medium-excitation lines ($130\,{\rm K}\, \le\, E_{\rm up}\,\le\,350\,{\rm K}$) typically display line shapes similar to CO. We analyze the HIFI observation together with archival SPIRE/PACS H$_2$O data using a state-of-the-art 3D radiative transfer code which includes the interaction between continuum and line emission. The water excitation models are combined with information on the dust- and CO spectral line energy distribution to determine the physical structure of the interstellar medium (ISM). We identify two ISM components that are common to all galaxies: A warm ($T_{\rm dust}\,\sim\,40-70\,{\rm K}$), dense ($n({\rm H})\,\sim\,10^5-10^6\,{\rm cm^{-3}}$) phase which dominates the emission of medium-excitation H$_2$O lines. This gas phase also dominates the FIR emission and the CO intensities for $J_{\rm up} > 8$. In addition a cold ($T_{\rm dust}\,\sim\,20-30\,{\rm K}$), dense ($n({\rm H})\sim\,10^4- 10^5\,{\rm cm^{-3}}$) more extended phase is present. It outputs the emission in the low-excitation H$_2$O lines and typically also produces the prominent line absorption features. For the two ULIRGs in our sample (Arp 220 and Mrk 231) an even hotter and more compact (R$_s\,\le\,100$ pc) region is present which is possibly linked to AGN activity. We find that collisions dominate the water excitation in the cold gas and for lines with $E_{\rm up}\le300\,{\rm K}$ and $E_{\rm up}\le800\,{\rm K}$ in the warm and hot component, respectively. Higher energy levels are mainly excited by IR pumping. △ Less

Submitted 16 July, 2017; originally announced July 2017.

Comments: Accepted by ApJ, in press

Disentangling the Circumnuclear Environs of Centaurus A: III. An Inner Molecular Ring, Nuclear Shocks and the CO to warm H2 interface

Authors: D. Espada, S. Matsushita, R. E. Miura, F. P. Israel, N. Neumayer, S. Martin, C. Henkel, T. Izumi, D. Iono, S. Aalto, J. Ott, A. B. Peck, A. C. Quillen, K. Kohno

Abstract: We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist with… ▽ More We present the distribution and kinematics of the molecular gas in the circumnuclear disk (CND, 400 pc x 200 pc) of Centaurus A with resolutions of ~5 pc (0.3 arcsec) and shed light onto the mechanism feeding the Active Galactic Nucleus (AGN) using CO(3-2), HCO+(4-3), HCN(4-3), and CO(6-5) observations obtained with ALMA. Multiple filaments or streamers of tens to a hundred parsec scale exist within the CND, which form a ring-like structure with an unprojected diameter of 9 x 6 arcsec (162pc x 108pc) and a position angle PA = 155deg. Inside the nuclear ring, there are two leading and straight filamentary structures with lengths of about 30-60pc at PA = 120deg on opposite sides of the AGN, with a rotational symmetry of 180deg and steeper position-velocity diagrams, which are interpreted as nuclear shocks due to non-circular motions. Along the filaments, and unlike other nearby AGNs, several dense molecular clumps present low HCN/HCO+(4-3) ratios (~0.5). The filaments abruptly end in the probed transitions at r = 20pc from the AGN, but previous near-IR H2 (J=1-0) S(1) maps show that they continue in an even ~1000 K), winding up in the form of nuclear spirals, and forming an inner ring structure with another set of symmetric filaments along the N-S direction and within r = 10pc. The molecular gas is governed primarily by non-circular motions, being the successive shock fronts at different scales where loss of angular momentum occurs, a mechanism which may feed efficiently powerful radio galaxies down to parsec scales. △ Less

Submitted 18 June, 2017; originally announced June 2017.

Comments: 46 pages. Accepted for publication in ApJ

The outflow of gas from the Centaurus A circumnuclear disk: atomic spectral line maps from Herschel-PACS and APEX

Authors: F. P. Israel, R. Guesten, R. Meijerink, M. A. Requena-Torres, J. Stutzki

Abstract: The physical state of the gas in the central 500 pc of NGC~5128 (the radio galaxy Centaurus A - Cen A), was investigated using the far-infrared fine-structure lines of carbon, oxygen, and nitrogen, as well as the CO(4-3) molecular line. The circumnuclear disk (CND) is traced by emission from dust and the neutral gas ([CI] and CO). A gas outflow with a line-of-sight velocity of 60 km/s is evident i… ▽ More The physical state of the gas in the central 500 pc of NGC~5128 (the radio galaxy Centaurus A - Cen A), was investigated using the far-infrared fine-structure lines of carbon, oxygen, and nitrogen, as well as the CO(4-3) molecular line. The circumnuclear disk (CND) is traced by emission from dust and the neutral gas ([CI] and CO). A gas outflow with a line-of-sight velocity of 60 km/s is evident in both species. The center of the CND is bright in [OI], [OIII], and [CII]; [OI]63mu emission dominates that of [CII] even though it is absorbed with optical depths of 1.0-1.5. The outflow is well-traced by the [NII] and [NIII] lines and also seen in the [CII] and [OIII] lines that peak in the center. Ionized gas densities are moderate in the CND and low everywhere else. Neutral gas densities range from 4000 per cm3 (outflow, extended thin disk ETD) to 20 000 per cm3 (CND). The CND radiation field is weak compared to the ETD starburst field. The outflow has a much stronger radiation field. The total mass of all the CND gas is 9 x 10^(7) M(o) and the mass of the outflowing gas is only 15%-30% of that. The outflow most likely originates from the shock-dominated CND cavity surrounding the central black hole. With a factor of three uncertainty, the mass outflow rate is about 2 M(o)/yr, a thousand times higher than the accretion rate of the black hole. Without replenishment, the CND will be depleted in 15-120 million years. However, the outflow velocity is well below the escape velocity. △ Less

Submitted 17 November, 2016; originally announced November 2016.

Comments: 14 pages, 9 figures. Accepted by Astronomy & Astrophysics

Journal ref: A&A 599, A53 (2017)

The LOFAR Two-metre Sky Survey - I. Survey Description and Preliminary Data Release

Authors: T. W. Shimwell, H. J. A. Röttgering, P. N. Best, W. L. Williams, T. J. Dijkema, F. de Gasperin, M. J. Hardcastle, G. H. Heald, D. N. Hoang, A. Horneffer, H. Intema, E. K. Mahony, S. Mandal, A. P. Mechev, L. Morabito, J. B. R. Oonk, D. Rafferty, E. Retana-Montenegro, J. Sabater, C. Tasse, R. J. van Weeren, M. Brüggen, G. Brunetti, K. T. Chyży, J. E. Conway , et al. (47 additional authors not shown)

Abstract: The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolutio… ▽ More The LOFAR Two-metre Sky Survey (LoTSS) is a deep 120-168 MHz imaging survey that will eventually cover the entire Northern sky. Each of the 3170 pointings will be observed for 8 hrs, which, at most declinations, is sufficient to produce ~5arcsec resolution images with a sensitivity of ~0.1mJy/beam and accomplish the main scientific aims of the survey which are to explore the formation and evolution of massive black holes, galaxies, clusters of galaxies and large-scale structure. Due to the compact core and long baselines of LOFAR, the images provide excellent sensitivity to both highly extended and compact emission. For legacy value, the data are archived at high spectral and time resolution to facilitate subarcsecond imaging and spectral line studies. In this paper we provide an overview of the LoTSS. We outline the survey strategy, the observational status, the current calibration techniques, a preliminary data release, and the anticipated scientific impact. The preliminary images that we have released were created using a fully-automated but direction-independent calibration strategy and are significantly more sensitive than those produced by any existing large-area low-frequency survey. In excess of 44,000 sources are detected in the images that have a resolution of 25arcsec, typical noise levels of less than 0.5 mJy/beam, and cover an area of over 350 square degrees in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45d00m00s to 57d00m00s). △ Less

Submitted 8 November, 2016; originally announced November 2016.

Comments: 23 pages, 15 figures, 3 tables, accepted for publication in A&A

Journal ref: A&A 598, A104 (2017)

Constraining cloud parameters using high density gas tracers in galaxies

Authors: M. V. Kazandjian, I. Pelupessy, R. Meijerink, F. P. Israel, C. M. Coppola, M. J. F. Rosenberg, M. Spaans

Abstract: Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic… ▽ More Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion σ= 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3. We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77. △ Less

Submitted 8 September, 2016; originally announced September 2016.

Comments: 13 pages, 11 figures

Modelling Mechanical Heating in Star-Forming Galaxies: CO and 13CO Line Ratios as Sensitive Probes

Authors: M. V. Kazandjian, I. Pelupessy, R. Meijerink, F. P. Israel, M. Spaans

Abstract: We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider t… ▽ More We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models. We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$. △ Less

Submitted 8 September, 2016; originally announced September 2016.

Comments: 17 pages, 12 figures

Planck intermediate results. XLV. Radio spectra of northern extragalactic radio sources

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, R. Battye, K. Benabed, G. J. Bendo, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, C. Burigana, R. C. Butler , et al. (180 additional authors not shown)

Abstract: Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the fla… ▽ More Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts. △ Less

Submitted 16 June, 2016; originally announced June 2016.

Comments: 39 pages, 113 figures. Accepted by A&A

Continuum Observations of M51 and M83 at 1.1 mm with AzTEC

Authors: W. F. Wall, I. Puerari, R. Tilanus, F. P. Israel, J. E. Austermann, I. Aretxaga, G. Wilson, M. Yun, K. S. Scott, T. A. Perera, C. M. Roberts, D. H. Hughes

Abstract: We observed the spiral galaxies M51 and M83 at 20" spatial resolution with the bolometer array AzTEC on the JCMT in the 1.1$\,$mm continuum, recovering the extended emission out to galactocentric radii of more than 12 kpc in both galaxies. The 1.1 mm-continuum fluxes are 5.6+/-0.7 and 9.9+/-1.4 Jy, with associated gas masses estimated at 9.4 X 10^9 Mo and 7.2 X 10^9 Mo for M51 and M83, respectivel… ▽ More We observed the spiral galaxies M51 and M83 at 20" spatial resolution with the bolometer array AzTEC on the JCMT in the 1.1$\,$mm continuum, recovering the extended emission out to galactocentric radii of more than 12 kpc in both galaxies. The 1.1 mm-continuum fluxes are 5.6+/-0.7 and 9.9+/-1.4 Jy, with associated gas masses estimated at 9.4 X 10^9 Mo and 7.2 X 10^9 Mo for M51 and M83, respectively. In the interarm regions of both galaxies the N(H2)/I(CO) (or X-factor) ratios exceed those in the arms by factors of ~1.5-2. In the inner disks of both galaxies, the X-factor is about 1 X 10^20 cm^-2 / (K km s^-1). In the outer parts, the CO-dark molecular gas becomes more important. While the spiral density wave in M51 appears to influence the interstellar medium and stars in a similar way, the bar potential in M83 influences the interstellar medium and the stars differently. We confirm the result of Foyle et al. (2010) that the arms merely heighten the star formation rate and the gas surface density in the same proportion. Our maps reveal a threshold gas surface density for an SFR increase by two or more orders of magnitude. In both galaxy centers, the molecular gas depletion time is about 1 Gyr climbing to 10-20 Gyr at radii of 6-8 kpc. This is consistent with an inside-out depletion of the molecular gas in the disks of spiral galaxies. △ Less

Submitted 24 March, 2016; originally announced March 2016.

Comments: 34 pages, 22 figures. Accepted by MNRAS

Carbon gas in SMC low-metallicity star-forming regions

Authors: M. A. Requena-Torres, F. P. Israel, Y. Okada, R. Guesten, J. Stutzki, C. Risacher, R. Simon, H. Zinnecker

Abstract: This paper presents [CII], [CI] and CO emission line maps of the star-forming regions N66, N25+N26, and N88 in the metal-poor Local Group dwarf galaxy SMC. The spatial and velocity structure of the large HII region N66 reveals an expanding ring of shocked molecular gas centered on the exciting star cluster NGC346, whereas a more distant dense molecular cloud is being eroded by UV radiation from th… ▽ More This paper presents [CII], [CI] and CO emission line maps of the star-forming regions N66, N25+N26, and N88 in the metal-poor Local Group dwarf galaxy SMC. The spatial and velocity structure of the large HII region N66 reveals an expanding ring of shocked molecular gas centered on the exciting star cluster NGC346, whereas a more distant dense molecular cloud is being eroded by UV radiation from the same cluster. In the N25+N26 and N88 maps, diffuse [CII] emission at a relatively low surface brightness extends well beyond the compact boundaries of the bright emission associated with the [CII] regions. In all regions, the distribution of this bright [CII] emission and the less prominent [CI] emission closely follows the outline of the CO complexes, but the intensity of the [CII] and [CI] emission is generally anticorrelated, which can be understood by the action of photodissociation and photoionization processes. Notwithstanding the overall similarity of CO and [CII] maps, the intensity ratio of these lines varies significantly, mostly due to changes in CO brightness. [CII] emission line profiles are up to 50% wider in velocity than corresponding CO profiles. A radiative transfer analysis shows that the [CII] line is the dominant tracer of (CO-dark) molecular hydrogen in the SMC. CO emission traces only a minor fraction of the total amount of gas. The similarity of the spatial distribution and line profile shape, and the dominance of molecular gas associated with [[CII] rather than CO emission imply that in the low-metallicity environment of the SMC the small amount of dense molecular gas traced by CO is embedded in the much more extended molecular gas traced only by [CII] emission. The contribution from neutral atomic and ionized hydrogen zones is negligible in the star-forming regions observed. △ Less

Submitted 24 February, 2016; originally announced February 2016.

Comments: 17 pages, 16 figures, Accepted in A&A

Journal ref: A&A 589, A28 (2016)

The JCMT Nearby Galaxies Legacy Survey X. Environmental Effects on the Molecular Gas and Star Formation Properties of Spiral Galaxies

Authors: Angus Mok, C. D. Wilson, J. Golding, B. E. Warren, F. P. Israel, S. Serjeant, J. H. Knapen, J. R. Sanchez-Gallego, P. Barmby, G. J. Bendo, E. Rosolowsky, P. van der Werf

Abstract: We present a study of the molecular gas properties in a sample of 98 HI - flux selected spiral galaxies within $\sim25$ Mpc, using the CO $J=3-2$ line observed with the James Clerk Maxwell Telescope. We use the technique of survival analysis to incorporate galaxies with CO upper limits into our results. Comparing the group and Virgo samples, we find a larger mean H$_{2}$ mass in the Virgo galaxies… ▽ More We present a study of the molecular gas properties in a sample of 98 HI - flux selected spiral galaxies within $\sim25$ Mpc, using the CO $J=3-2$ line observed with the James Clerk Maxwell Telescope. We use the technique of survival analysis to incorporate galaxies with CO upper limits into our results. Comparing the group and Virgo samples, we find a larger mean H$_{2}$ mass in the Virgo galaxies, despite their lower mean HI mass. This leads to a significantly higher H$_{2}$ to HI ratio for Virgo galaxies. Combining our data with complementary H$α$ star formation rate measurements, Virgo galaxies have longer molecular gas depletion times compared to group galaxies, due to their higher H$_{2}$ masses and lower star formation rates. We suggest that the longer depletion times may be a result of heating processes in the cluster environment or differences in the turbulent pressure. From the full sample, we find that the molecular gas depletion time has a positive correlation with the stellar mass, indicative of differences in the star formation process between low and high mass galaxies, and a negative correlation between the molecular gas depletion time and the specific star formation rate. △ Less

Submitted 17 December, 2015; originally announced December 2015.

Comments: Accepted for publication in MNRAS, 25 pages, 14 figures

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

Authors: Katherine E. Jameson, Alberto D. Bolatto, Adam K. Leroy, Margaret Meixner, Julia Roman-Duval, Karl Gordon, Annie Hughes, Frank P. Israel, Monica Rubio, Remy Indebetouw, Suzanne C. Madden, Caroline Bot, Sacha Hony, Diane Cormier, Eric W. Pellegrini, Maud Galametz, George Sonneborn

Abstract: 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 g… ▽ More 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. △ Less

Submitted 25 April, 2016; v1 submitted 27 October, 2015; originally announced October 2015.

Comments: 24 pages, 14 figures, accepted for publication in ApJ. FITS files of the dust-based estimates of the H2 column densities for the LMC and SMC (shown in Figures 2 and 3) will be available online through ApJ

Neutral carbon and CO in 76 (U)LIRGs and starburst galaxy centers A method to determine molecular gas properties in luminous galaxies

Authors: F. P. Israel, M. J. F. Rosenberg, P. van der Werf

Abstract: We present fluxes in both neutral carbon [CI] lines at the centers of 76 galaxies with FIR luminosities between 10^{9} and 10^{12} L(o) obtained with Herschel-SPIRE and with ground-based facilities, along with the J=7-6, J=4-3, J=2-1 12CO and J=2-1 13CO line fluxes. We investigate whether these lines can be used to characterize the molecular ISM of the parent galaxies in simple ways and how the mo… ▽ More We present fluxes in both neutral carbon [CI] lines at the centers of 76 galaxies with FIR luminosities between 10^{9} and 10^{12} L(o) obtained with Herschel-SPIRE and with ground-based facilities, along with the J=7-6, J=4-3, J=2-1 12CO and J=2-1 13CO line fluxes. We investigate whether these lines can be used to characterize the molecular ISM of the parent galaxies in simple ways and how the molecular gas properties define the model results. In most starburst galaxies, the [CI]/13CO flux ratio is much higher than in Galactic star-forming regions, and it is correlated to the total FIR luminosity. The [CI](1-0)/CO(4-3), the [CI](2-1) (2-1)/CO(7-6), and the [CI] (2-1)/(1-0) flux ratios are also correlated, and trace the excitation of the molecular gas. In the most luminous infrared galaxies (LIRGs), the ISM is fully dominated by dense and moderately warm gas clouds that appear to have low [C]/[CO] and [13CO]/[12CO] abundances. In less luminous galaxies, emission from gas clouds at lower densities becomes progressively more important, and a multiple-phase analysis is required to determine consistent physical characteristics. Neither the CO nor the [CI] velocity-integrated line fluxes are good predictors of H2 column densities in individual galaxies, and X(CI) conversion factors are not superior to X(CO) factors. The methods and diagnostic diagrams outlined in this paper also provide a new and relatively straightforward means of deriving the physical characteristics of molecular gas in high-redshift galaxies up to z=5, which are otherwise hard to determine. △ Less

Submitted 29 April, 2015; originally announced April 2015.

Journal ref: A&A 578, A95 (2015)

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

Authors: M. J. F. Rosenberg, P. P. van der Werf, S. Aalto, L. Armus, V. Charmandaris, T. Díaz-Santos, A. S. Evans, J. Fischer, Y. Gao, E. González-Alfonso, T. R. Greve, A. I. Harris, C. Henkel, F. P. Israel, K. G. Isaak, C. Kramer, R. Meijerink, D. A. Naylor, D. B. Sanders, H. A. Smith, M. Spaans, L. Spinoglio, G. J. Stacey, I. Veenendaal, S. Veilleux , et al. (5 additional authors not shown)

Abstract: (Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10… ▽ More (Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important. △ Less

Submitted 13 January, 2015; originally announced January 2015.

Planck intermediate results. XXV. The Andromeda Galaxy as seen by Planck

Authors: Planck Collaboration, P. A. R. Ade, N. Aghanim, M. Arnaud, M. Ashdown, J. Aumont, C. Baccigalupi, A. J. Banday, R. B. Barreiro, N. Bartolo, E. Battaner, R. Battye, K. Benabed, G. J. Bendo, A. Benoit-Lévy, J. -P. Bernard, M. Bersanelli, P. Bielewicz, A. Bonaldi, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, C. Burigana, R. C. Butler , et al. (180 additional authors not shown)

Abstract: The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Plan… ▽ More The Andromeda Galaxy (M31) is one of a few galaxies that has sufficient angular size on the sky to be resolved by the Planck satellite. Planck has detected M31 in all of its frequency bands, and has mapped out the dust emission with the High Frequency Instrument, clearly resolving multiple spiral arms and sub-features. We examine the morphology of this long-wavelength dust emission as seen by Planck, including a study of its outermost spiral arms, and investigate the dust heating mechanism across M31. We find that dust dominating the longer wavelength emission ($\gtrsim 0.3\,$mm) is heated by the diffuse stellar population (as traced by 3.6$\,μ$m emission), with the dust dominating the shorter wavelength emission heated by a mix of the old stellar population and star-forming regions (as traced by 24$\,μ$m emission). We also fit spectral energy distributions (SEDs) for individual 5' pixels and quantify the dust properties across the galaxy, taking into account these different heating mechanisms, finding that there is a linear decrease in temperature with galactocentric distance for dust heated by the old stellar population, as would be expected, with temperatures ranging from around 22$\,$K in the nucleus to 14$\,$K outside of the 10$\,$kpc ring. Finally, we measure the integrated spectrum of the whole galaxy, which we find to be well-fitted with a global dust temperature of ($18.2\pm1.0$)$\,$K with a spectral index of $1.62\pm0.11$ (assuming a single modified blackbody), and a significant amount of free-free emission at intermediate frequencies of 20-60$\,$GHz, which corresponds to a star formation rate of around $0.12$M$_\odot\,$yr$^{-1}$. We find a $2.3\,σ$ detection of the presence of spinning dust emission, with a 30$\,$GHz amplitude of $0.7\pm0.3\,$Jy, which is in line with expectations from our Galaxy. △ Less

Submitted 8 September, 2015; v1 submitted 21 July, 2014; originally announced July 2014.

Comments: 26 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 582, A28 (2015)

Molecular gas heating in Arp 299

Authors: M. J. F. Rosenberg, R. Meijerink, F. P. Israel, P. P. van der Werf, E. M. Xilouris, A. Weiß

Abstract: Understanding the heating and cooling mechanisms in nearby (Ultra) luminous infrared galaxies can give us insight into the driving mechanisms in their more distant counterparts. Molecular emission lines play a crucial role in cooling excited gas, and recently, with Herschel Space Observatory we have been able to observe the rich molecular spectrum. CO is the most abundant and one of the brightest… ▽ More Understanding the heating and cooling mechanisms in nearby (Ultra) luminous infrared galaxies can give us insight into the driving mechanisms in their more distant counterparts. Molecular emission lines play a crucial role in cooling excited gas, and recently, with Herschel Space Observatory we have been able to observe the rich molecular spectrum. CO is the most abundant and one of the brightest molecules in the Herschel wavelength range. CO transitions are observed with Herschel, and together, these lines trace the excitation of CO. We study Arp 299, a colliding galaxy group, with one component harboring an AGN and two more undergoing intense star formation. For Arp 299 A, we present PACS spectrometer observations of high-J CO lines up to J=20-19 and JCMT observations of $^{13}$CO and HCN to discern between UV heating and alternative heating mechanisms. There is an immediately noticeable difference in the spectra of Arp 299 A and Arp 299 B+C, with source A having brighter high-J CO transitions. This is reflected in their respective spectral energy line distributions. We find that photon-dominated regions (PDRs) are unlikely to heat all the gas since a very extreme PDR is necessary to fit the high-J CO lines. In addition, this extreme PDR does not fit the HCN observations, and the dust spectral energy distribution shows that there is not enough hot dust to match the amount expected from such an extreme PDR. Therefore, we determine that the high-J CO and HCN transitions are heated by an additional mechanism, namely cosmic ray heating, mechanical heating, or X-ray heating. We find that mechanical heating, in combination with UV heating, is the only mechanism that fits all molecular transitions. We also constrain the molecular gas mass of Arp 299 A to 3e9 Msun and find that we need 4% of the total heating to be mechanical heating, with the rest UV heating. △ Less

Submitted 8 July, 2014; originally announced July 2014.

Journal ref: A&A 568, A90 (2014)

The nature of the low-frequency emission of M51: First observations of a nearby galaxy with LOFAR

Authors: D. D. Mulcahy, A. Horneffer, R. Beck, G. Heald, A. Fletcher, A. Scaife, B. Adebahr, J. M. Anderson, A. Bonafede, M. Brüggen, G. Brunetti, K. T. Chyży, J. Conway, R-J. Dettmar, T. Enßlin, M. Haverkorn, C. Horellou, M. Iacobelli, F. P. Israel, H. Junklewitz, W. Jurusik, J. Köhler, M. Kuniyoshi, E. Orrú, R. Paladino , et al. (3 additional authors not shown)

Abstract: The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies… ▽ More The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates. △ Less

Submitted 22 July, 2014; v1 submitted 4 July, 2014; originally announced July 2014.

Comments: 20 pages, 18 figures, accepted for publication in A&A

Journal ref: A&A 568, A74 (2014)

Dust and Gas in the Magellanic Clouds from the HERITAGE Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm Excess Emission

Authors: Karl D. Gordon, Julia Roman-Duval, Caroline Bot, Margaret Meixner, Brian Babler, Jean-Philippe Bernard, Alberto Bolatto, Martha L. Boyer, Geoffrey C. Clayton, Charles Engelbracht, Yasuo Fukui, Maud Galametz, Frederic Galliano, Sacha Hony, Annie Hughes, Remy Indebetouw, Frank P. Israel, Katie Jameson, Akiko Kawamura, Vianney Lebouteiller, Aigen Li, Suzanne C. Madden, Mikako Matsuura, Karl Misselt, Edward Montiel , et al. (11 additional authors not shown)

Abstract: The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two… ▽ More The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micron. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power- law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submm excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micron. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micron submm excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our the fitting results shows that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submm excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 +/- 1.7) x 10^5 and (8.3 +/- 2.1) times 10^4 M(sun) for the LMC and SMC, respectively. We find significant correlations between the submm excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations. △ Less

Submitted 23 June, 2014; originally announced June 2014.

Comments: 20 pages, 9 figures, accepted for publication in the ApJ

Molecular gas heating mechanisms, and star formation feedback in merger/s tarbursts: NGC 6240 and Arp 193 as case studies

Authors: Padelis P. Papadopoulos, Zhi-Yu Zhang, E. M. Xilouris, Axel Weiss, Paul van der Werf, Frank. P. Israel, Thomas. R. Greve, Kate G. Isaak, Yu Gao

Abstract: We used the SPIRE/FTS instrument aboard the Herschel Space Observatory (HSO) to obtain the Spectral Line Energy Distributions (SLEDs) of CO from J=4-3 to J=13-12 of Arp 193 and NGC 6240, two classical merger/starbursts selected from our molecular line survey of local Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11} L_{sol}). The high-J CO SLEDs are then combined with ground-based low-J CO, {13}C… ▽ More We used the SPIRE/FTS instrument aboard the Herschel Space Observatory (HSO) to obtain the Spectral Line Energy Distributions (SLEDs) of CO from J=4-3 to J=13-12 of Arp 193 and NGC 6240, two classical merger/starbursts selected from our molecular line survey of local Luminous Infrared Galaxies (LIRGs: L_{IR}>=10^{11} L_{sol}). The high-J CO SLEDs are then combined with ground-based low-J CO, {13}CO, HCN, HCO+, CS line data and used to probe the thermal and dynamical states of their large molecular gas reservoirs. We find the two CO SLEDs strongly diverging from J=4-3 onwards, with NGC6240 having a much higher CO line excitation than Arp193, despite their similar low-J CO SLEDs and L_{FIR}/L_{CO,1-0}, L_{HCN}/L_{CO} (J=1-0) ratios (proxies of star formation efficiency and dense gas mass fraction). In Arp193, one of the three most extreme starbursts in the local Universe, the molecular SLEDs indicate a small amount ~(5-15)% of dense gas (n>=10^{4}cm^{-3}) unlike NGC6240 where most of the molecular gas (~(60-70)%) is dense n~(10^4-10^5)cm^{-3}. Strong star-formation feedback can drive this disparity in their dense gas mass fractions, and also induce extreme thermal and dynamical states for the molecular gas.In NGC6240, and to a lesser degree in Arp193, we find large molecular gas masses whose thermal states cannot be maintained by FUV photons from Photon Dominated Regions (PDRs). We argue that this may happen often in metal-rich merger/starbursts, strongly altering the initial conditions of star formation. ALMA can now directly probe these conditions across cosmic epoch, and even probe their deeply dust-enshrouded outcome, the stellar IMF averaged over galactic evolution. △ Less

Submitted 24 April, 2014; originally announced April 2014.

Comments: 54 pages, 19 figures, accepted for publication in The Astrophysical Journal

Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios

Authors: M. V. Kazandjian, R. Meijerink, I. Pelupessy, F. P. Israel, M. Spaans

Abstract: CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an… ▽ More CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates. △ Less

Submitted 27 March, 2014; originally announced March 2014.

Comments: 38 pages, to appear in A&A

The molecular circumnuclear disk (CND) in Centaurus A: A multi-transition CO and [CI] survey with Herschel, APEX, JCMT, and SEST

Authors: F. P. Israel, R. Guesten, R. Meijerink, A. F. Loenen, M. A. Requena-Torres, J. Stutzki, P. van der Werf, A. Harris, C. Kramer, J. Martin-Pintado, A. Weiss

Abstract: We present new CO and C^o line measurements of the compact circumnuclear disk in the center of NGC 128 (Centaurus~A) obtained with the Herschel Space Observatory, as well as SEST, JCMT, and APEX. The Cen A center CO ladder is quite different from those of either star-burst galaxies or AGNs. In addition, the relative intensity of the central Cen A [CI] emission lines is much greate than that in any… ▽ More We present new CO and C^o line measurements of the compact circumnuclear disk in the center of NGC 128 (Centaurus~A) obtained with the Herschel Space Observatory, as well as SEST, JCMT, and APEX. The Cen A center CO ladder is quite different from those of either star-burst galaxies or AGNs. In addition, the relative intensity of the central Cen A [CI] emission lines is much greate than that in any other galaxy. The CO surface brightness of the compact circumnuclear disk (CND) is significantly higher than that of the much more extended thin disk (ETD) in the same line of sight. Our LVG and PDR/XDR models suggest that much of the CND gas is relatively cool (25 - 80 K) and not very dense (~ 300 cm^{-3}) if the heating is by UV photons, although there is some gas in both the CND and the ETD with a much higher density of ~30 000 cm^{-3}. Finally, there is also high-excitation, high-density phase in the CND (but not in the ETD), either in the form of an extreme PDR but more likely in the form of an XDR. The total gas mass of the Cen A CND is 8.4 x 10^{7} M(sun), uncertain by a factor of two. The CO-H2 conversion factor is 4 x 10^{20} K km/s, also within a factor of two. △ Less

Submitted 5 February, 2014; originally announced February 2014.

Comments: 14 pages, 5 figures; A&A in press

Radiative and mechanical feedback into the molecular gas of NGC 253

Authors: M. J. F. Rosenberg, M. V. Kazandjian, P. P. van der Werf, F. P. Israel, R. Meijerink, A. Weiß, M. A. Requena-Torres, R. Güsten

Abstract: Starburst galaxies are undergoing intense periods of star formation. Understanding the heating and cooling mechanisms in these galaxies can give us insight to the driving mechanisms that fuel the starburst. Molecular emission lines play a crucial role in the cooling of the excited gas. With SPIRE on the Herschel Space Observatory we have observed the rich molecular spectrum towards the central reg… ▽ More Starburst galaxies are undergoing intense periods of star formation. Understanding the heating and cooling mechanisms in these galaxies can give us insight to the driving mechanisms that fuel the starburst. Molecular emission lines play a crucial role in the cooling of the excited gas. With SPIRE on the Herschel Space Observatory we have observed the rich molecular spectrum towards the central region of NGC 253. CO transitions from J=4-3 to 13-12 are observed and together with low-J line fluxes from ground based observations, these lines trace the excitation of CO. By studying the CO excitation ladder and comparing the intensities to models, we investigate whether the gas is excited by UV radiation, X-rays, cosmic rays, or turbulent heating. Comparing the $^{12}$CO and $^{13}$CO observations to large velocity gradient models and PDR models we find three main ISM phases. We estimate the density, temperature,and masses of these ISM phases. By adding $^{13}$CO, HCN, and HNC line intensities, we are able to constrain these degeneracies and determine the heating sources. The first ISM phase responsible for the low-J CO lines is excited by PDRs, but the second and third phases, responsible for the mid to high-J CO transitions, require an additional heating source. We find three possible combinations of models that can reproduce our observed molecular emission. Although we cannot determine which of these are preferable, we can conclude that mechanical heating is necessary to reproduce the observed molecular emission and cosmic ray heating is a negligible heating source. We then estimate the mass of each ISM phase; $6\times 10^7$ M$_\odot$ for phase 1 (low-J CO lines), $3\times 10^7$ M$_\odot$ for phase 2 (mid-J CO lines), and $9\times 10^6$ M$_\odot$ for phase 3 (high-J CO lines) for a total system mass of $1\times10^{8}$ M$_\odot$. △ Less

Submitted 20 January, 2014; originally announced January 2014.

A CO J=3-2 map of M51 with HARP-B: Radial properties of the spiral structure

Authors: C. Vlahakis, P. van der Werf, F. P. Israel, R. P. J. Tilanus

Abstract: We present the first complete CO J=3-2 map of the nearby grand-design spiral galaxy M51 (NGC 5194), at a spatial resolution of ~600 pc, obtained with the HARP-B instrument on the James Clerk Maxwell Telescope. The map covers the entire optical galaxy disk and out to the companion NGC 5195, with CO J=3-2 emission detected over an area of ~9'x6' (~21x14 kpc). We describe the CO J=3-2 integrated inte… ▽ More We present the first complete CO J=3-2 map of the nearby grand-design spiral galaxy M51 (NGC 5194), at a spatial resolution of ~600 pc, obtained with the HARP-B instrument on the James Clerk Maxwell Telescope. The map covers the entire optical galaxy disk and out to the companion NGC 5195, with CO J=3-2 emission detected over an area of ~9'x6' (~21x14 kpc). We describe the CO J=3-2 integrated intensity map and combine our results with maps of CO J=2-1, CO J=1-0 and other data from the literature to investigate the variation of the molecular gas, atomic gas and polycyclic aromatic hydrocarbon (PAH) properties of M51 as a function of distance along the spiral structure on sub-kpc scales. We find that for the CO J=3-2 and CO J=2-1 transitions there is a clear difference between the variation of arm and inter-arm emission with galactocentric radius, with the inter-arm emission relatively constant with radius and the contrast between arm and inter-arm emission decreasing with radius. For CO J=1-0 and HI the variation with radius shows a similar trend for the arm and inter-arm regions, and the arm-inter-arm contrast appears relatively constant with radius. We investigate the variation of CO line ratios (J=3-2/2-1, J=2-1/1-0 and J=3-2/1-0) as a function of distance along the spiral structure. Line ratios are consistent with the range of typical values for other nearby galaxies in the literature. The highest CO J=3-2/2-1 line ratios are found in the central ~1 kpc and in the spiral arms and the lowest line ratios in the inter-arm regions.We find no clear evidence of a trend with radius for the spiral arms but for the inter-arm regions there appears to be a trend for all CO line ratios to increase with radius. We find a strong relationship between the ratio of CO J=3-2 intensity to stellar continuum-subtracted 8mu PAH surface brightness and the CO J=3-2 intensity that appears to vary with radius. △ Less

Submitted 27 April, 2013; originally announced April 2013.

Comments: 26 pages, 18 figures, 7 tables

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

Authors: Cinthya N. Herrera, Monica Rubio, Alberto D. Bolatto, Francois Boulanger, Frank P. Israel, Fredrik T. Rantakyro

Abstract: 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… ▽ More 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. △ Less

Submitted 1 April, 2013; originally announced April 2013.

Comments: A&A accepted, 13 pages, 10 figures

Gas and dust cooling along the major axis of M33 (HerM33es): ISO/LWS CII observations

Authors: C. Kramer, J. Abreu-Vicente, S. Garcia-Burillo, M. Relano, S. Aalto, M. Boquien, J. Braine, C. Buchbender, P. Gratier, F. P. Israel, T. Nikola, M. Roellig, S. Verley, P. van der Werf, E. M. Xilouris

Abstract: We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis o… ▽ More We aim to better understand the heating of the gas by observing the prominent gas cooling line [CII] at 158um in the low-metallicity environment of the Local Group spiral galaxy M33 at scales of 280pc. In particular, we aim at describing the variation of the photoelectric heating efficiency with galactic environment. In this unbiased study, we used ISO/LWS [CII] observations along the major axis of M33, in combination with Herschel PACS and SPIRE continuum maps, IRAM 30m CO 2-1 and VLA HI data to study the variation of velocity integrated intensities. The ratio of [CII] emission over the far-infrared continuum is used as a proxy for the heating efficiency, and models of photon-dominated regions are used to study the local physical densities, FUV radiation fields, and average column densities of the molecular clouds. The heating efficiency stays constant at 0.8% in the inner 4.5kpc radius of the galaxy where it starts to increase to reach values of ~3% in the outskirts at about 6kpc radial distance. The rise of efficiency is explained in the framework of PDR models by lowered volume densities and FUV fields, for optical extinctions of only a few magnitudes at constant metallicity. In view of the significant fraction of HI emission stemming from PDRs, and for typical pressures found in the Galactic cold neutral medium (CNM) traced by HI emission, the CNM contributes ~15% to the observed [CII] emission in the inner 2kpc radius of M33. The CNM contribution remains largely undetermined in the south, while positions between 2 and 7.3kpc radial distance in the north of M33 show a contribution of ~40%+-20%. △ Less

Submitted 31 March, 2013; originally announced April 2013.

Comments: (16 pages, a version with higher resolution figures may be downloaded from http://www.iram.es/IRAMES/hermesWiki/HermesPublications)

The thermal dust emission in the N158-N159-N160 (LMC) star forming complex mapped by Spitzer, Herschel and LABOCA

Authors: M. Galametz, S. Hony, F. Galliano, S. C. Madden, M. Albrecht, C. Bot, D. Cormier, C. Engelbracht, Y. Fukui, F. P. Israel, A. Kawamura, V. Lebouteiller, A. Li, M. Meixner, K. Misselt, E. Montiel, K. Okumura, P. Panuzzo, J. Roman- Duval, M. Rubio, M. Sauvage, J. P. Seale, M. Sewilo, J. Th. van Loon

Abstract: We present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and r… ▽ More We present a study of the infrared/submm emission of the LMC star forming complex N158-N159-N160. Combining observations from the Spitzer Space Telescope (3.6-70um), the Herschel Space Observatory (100-500um) and LABOCA (870um) allows us to work at the best angular resolution available now for an extragalactic source. We observe a remarkably good correlation between SPIRE and LABOCA emission and resolve the low surface brightnesses emission. We use the Spitzer and Herschel data to perform a resolved Spectral Energy Distribution (SED) modelling of the complex. Using MBB, we derive a global emissivity index beta_c of 1.47. If beta cold is fixed to 1.5, we find an average temperature of 27K. We also apply the Galliano et al. (2011) modelling technique (and amorphous carbon to model carbon dust) to derive maps of the star formation rate, the mean starlight intensity, the fraction of PAHs or the dust mass surface density of the region. We observe that the PAH fraction strongly decreases in the HII regions. This decrease coincides with peaks in the mean radiation field intensity map. The dust surface densities follow the FIR distribution, with a total dust mass of 2.1x10^4 Msolar (2.8 times less than when using graphite grains) in the resolved elements we model. We find a non-negligible amount of dust in the molecular cloud N159 South (showing no massive SF). We also investigate the drivers of the Herschel/PACS and SPIRE submm colours as well as the variations in the gas-to-dust mass ratio (G/D) and the XCO conversion factor in the region N159. We finally model individual regions to analyse variations in the SED shape across the complex and the 870um emission in more details. No measurable submm excess emission at 870um seems to be detected in these regions. △ Less

Submitted 12 February, 2013; originally announced February 2013.

Comments: Accepted for publication in MNRAS; 23 pages, 11 figures, 2 tables

The excitation of near-infrared H2 emission in NGC 253

Authors: M. J. F. Rosenberg, P. P. van der Werf, F. P. Israel

Abstract: Because of its large angular size and proximity to the Milky Way, NGC 253, an archetypal starburst galaxy, provides an excellent laboratory to study the intricacies of this intense episode of star formation. We aim to characterize the excitation mechanisms driving the emission in NGC 253. Specifically we aim to distinguish between shock excitation and UV excitation as the dominant driving mechanis… ▽ More Because of its large angular size and proximity to the Milky Way, NGC 253, an archetypal starburst galaxy, provides an excellent laboratory to study the intricacies of this intense episode of star formation. We aim to characterize the excitation mechanisms driving the emission in NGC 253. Specifically we aim to distinguish between shock excitation and UV excitation as the dominant driving mechanism, using Brγ, H_2 and [FeII] as diagnostic emission line tracers. Using SINFONI observations, we create linemaps of Brγ, [FeII]_{1.64}, and all detected H_2 transitions. By using symmetry arguments of the gas and stellar gas velocity field, we find a kinematic center in agreement with previous determinations. The ratio of the 2-1 S(1) to 1-0 S(1) H_2 transitions can be used as a diagnostic to discriminate between shock and fluorescent excitation. Using the 1-0 S(1)/2-1 S(1) line ratio as well as several other H_2 line ratios and the morphological comparison between H_2 and Brγand [FeII], we find that excitation from UV photons is the dominant excitation mechanisms throughout NGC 253. We employ a diagnostic energy level diagram to quantitatively differentiate between mechanisms. We compare the observed energy level diagrams to PDR and shock models and find that in most regions and over the galaxy as a whole, fluorescent excitation is the dominant mechanism exciting the H_2 gas. We also place an upper limit of the percentage of shock excited H_2 at 29%. We find that UV radiation is the dominant excitation mechanism for the H_2 emission. The H_2 emission does not correlate well with Brγbut closely traces the PAH emission, showing that not only is H_2 fluorescently excited, but it is predominately excited by slightly lower mass stars than O stars which excite Brγ, such as B stars. △ Less

Submitted 7 December, 2012; originally announced December 2012.

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

Authors: R. Meijerink, L. E. Kristensen, A. Weiss, P. P. van der Werf, F. Walter, M. Spaans, A. F. Loenen, J. Fischer, F. P. Israel, K. Isaak, P. P. Papadopoulos, S. Aalto, L. Armus, V. Charmandaris, K. M. Dasyra, T. Diaz-Santos, A. Evans, Y. Gao, E. Gonzalez-Alfonso, R. Guesten, C. Henkel, C. Kramer, S. Lord, J. Martin-Pintado, D. Naylor , et al. (6 additional authors not shown)

Abstract: We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation… ▽ More We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks. △ Less

Submitted 28 November, 2012; originally announced November 2012.

Comments: Accepted by Astrophysical Journal Letters

Dense Gas in M33 (HerM33es)

Authors: C. Buchbender, C. Kramer, M. Gonzalez-Garcia, F. P. Israel, S. García-Burillo, P. van der Werf, J. Braine, E. Rosolowsky, B. Mookerjea, S. Aalto, M. Boquien, P. Gratier, C. Henkel, G. Quintana-Lacaci, S. Verley, F. van der Tak

Abstract: We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial di… ▽ More We aim to better understand the emission of molecular tracers of the diffuse and dense gas in giant molecular clouds and the influence that metallicity, optical extinction, density, far-UV field, and star formation rate have on these tracers. Using the IRAM 30m telescope, we detected HCN, HCO+, 12CO, and 13CO in six GMCs along the major axis of M33 at a resolution of ~ 114pc and out to a radial distance of 3.4kpc. Optical, far-infrared, and submillimeter data from Herschel and other observatories complement these observations. To interpret the observed molecular line emission, we created two grids of models of photon-dominated regions, one for solar and one for M33-type subsolar metallicity. The observed HCO+/HCN line ratios range between 1.1 and 2.5. Similarly high ratios have been observed in the Large Magellanic Cloud. The HCN/CO ratio varies between 0.4% and 2.9% in the disk of M33. The 12CO/13CO line ratio varies between 9 and 15 similar to variations found in the diffuse gas and the centers of GMCs of the Milky Way. Stacking of all spectra allowed HNC and C2H to be detected. The resulting HCO+/HNC and HCN/HNC ratios of ~ 8 and 6, respectively, lie at the high end of ratios observed in a large set of (ultra-)luminous infrared galaxies. HCN abundances are lower in the subsolar metallicity PDR models, while HCO+ abundances are enhanced. For HCN this effect is more pronounced at low optical extinctions. The observed HCO+/HCN and HCN/CO line ratios are naturally explained by subsolar PDR models of low optical extinctions between 4 and 10 mag and of moderate densities of n = 3x10^3 - 3x10^4 cm^-3, while the FUV field strength only has a small effect on the modeled line ratios. The line ratios are almost equally well reproduced by the solar-metallicity models, indicating that variations in metallicity only play a minor role in influencing these line ratios. △ Less

Submitted 31 January, 2013; v1 submitted 11 October, 2012; originally announced October 2012.

Comments: 21 pages, 9 figures, published in A&A

Spectrally resolved CII emission in M~33 ({\tt HerM33es}): Physical conditions and kinematics around BCLMP 691

Authors: J. Braine, P. Gratier, C. Kramer, F. P. Israel, F. van der Tak, B. Mookerjea, M. Boquien, F. Tabatabaei, P. van der Werf, C. Henkel

Abstract: This work presents high spectral resolution observations of the \CII\ line at 158 \micron, one of the major cooling lines of the interstellar medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an \HII\ region far north (3.3 kpc) in the disk of M 33, the \CII\ and CO line profiles show similar velocities within $0.5 \kms$, while the \HI\ line velocities are… ▽ More This work presents high spectral resolution observations of the \CII\ line at 158 \micron, one of the major cooling lines of the interstellar medium, taken with the HIFI heterodyne spectrometer on the Herschel satellite. In BCLMP 691, an \HII\ region far north (3.3 kpc) in the disk of M 33, the \CII\ and CO line profiles show similar velocities within $0.5 \kms$, while the \HI\ line velocities are systematically shifted towards lower rotation velocities by $\sim 5\kms$. Observed at the same $12"$ angular resolution, the \CII\ lines are broader than those of CO by about 50% but narrower than the \HI\ lines. The \CII\ line to far-infrared continuum ratio suggests a photoelectric heating efficiency of 1.1%. The data, together with published models indicate a UV field $G_0 \sim 100$ in units of the solar neighborhood value, a gas density $n_H \sim 1000 \cc$, and a gas temperature $T\sim 200$ K. Adopting these values, we estimate the C$^+$ column density to be $N_{C^+} \approx 1.3 \times 10^{17} \cmt$. The \CII\ emission comes predominantly from the warm neutral region between the \HII\ region and the cool molecular cloud behind it. From published abundances, the inferred C$^+$ column corresponds to a hydrogen column density of $N_H \sim 2 \times 10^{21} \cmt$. The CO observations suggest that $N_H = 2 N_{H_2} \sim 3.2 \times 10^{21} \cmt$ and 21cm measurements, also at $12"$ resolution, yield $N_\HI \approx 1.2 \times 10^{21} \cmt$ within the \CII\ velocity range. Thus, some H$_2$ not detected in CO must be present, in agreement with earlier findings based on the SPIRE 250 -- 500 $μ$m emission. △ Less

Submitted 13 July, 2012; originally announced July 2012.

Comments: accepted in A&A

The JCMT Nearby Galaxies Legacy Survey VIII. CO data and the L(CO3-2)-L(FIR) correlation in the SINGS sample

Authors: C. D. Wilson, B. E. Warren, F. P. Israel, S. Serjeant, D. Attewell, G. J. Bendo, H. M. Butner, P. Chanial, D. L. Clements, J. Golding, V. Heesen, J. Irwin, J. Leech, H. E. Matthews, S. Muhle, A. M. J. Mortier, G. Petitpas, J. R. Sanchez-Gallego, E. Sinukoff, K. Shorten, B. K. Tan, R. P. J. Tilanus, A. Usero, M. Vaccari, T. Wiegert , et al. (35 additional authors not shown)

Abstract: The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the… ▽ More The James Clerk Maxwell Telescope Nearby Galaxies Legacy Survey (NGLS) comprises an HI-selected sample of 155 galaxies spanning all morphological types with distances less than 25 Mpc. We describe the scientific goals of the survey, the sample selection, and the observing strategy. We also present an atlas and analysis of the CO J=3-2 maps for the 47 galaxies in the NGLS which are also part of the Spitzer Infrared Nearby Galaxies Survey. We find a wide range of molecular gas mass fractions in the galaxies in this sample and explore the correlation of the far-infrared luminosity, which traces star formation, with the CO luminosity, which traces the molecular gas mass. By comparing the NGLS data with merging galaxies at low and high redshift which have also been observed in the CO J=3-2 line, we show that the correlation of far-infrared and CO luminosity shows a significant trend with luminosity. This trend is consistent with a molecular gas depletion time which is more than an order of magnitude faster in the merger galaxies than in nearby normal galaxies. We also find a strong correlation of the L(FIR)/L(CO3-2) ratio with the atomic to molecular gas mass ratio. This correlation suggests that some of the far-infrared emission originates from dust associated with atomic gas and that its contribution is particularly important in galaxies where most of the gas is in the atomic phase. △ Less

Submitted 7 June, 2012; originally announced June 2012.

Comments: 37 pages, 107 figures; accepted for publication in MNRAS. Data and additional information are available at http://www.physics.mcmaster.ca/~wilson/www_xfer/NGLS/

Diagnostics of the Molecular Component of PDRs with Mechanical Heating

Authors: M. V. Kazandjian, R. Meijerink, I. Pelupessy, F. P. Israel, M. Spaans

Abstract: Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the… ▽ More Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR. △ Less

Submitted 22 March, 2012; originally announced March 2012.

Comments: 26 pages, 14 figures in the text and 13 figures as supplementary material. Accepted for publication in A&A

[FeII] as a tracer of Supernova Rate in Near-by Starburst Galaxies

Authors: M. J. F. Rosenberg, P. P. van der Werf, F. P. Israel

Abstract: Supernovae play an integral role in the feedback of processed material into the ISM of galaxies and are responsible for much of the chemical enrichment of the universe. The rate of supernovae can also reveal the star formation histories. Supernova rates are usually measured through the non-thermal radio continuum luminosity; however, a correlation between near-infrared [FeII] emission and supernov… ▽ More Supernovae play an integral role in the feedback of processed material into the ISM of galaxies and are responsible for much of the chemical enrichment of the universe. The rate of supernovae can also reveal the star formation histories. Supernova rates are usually measured through the non-thermal radio continuum luminosity; however, a correlation between near-infrared [FeII] emission and supernova remnants has also been noted. We aim to find a quantitative relationship between the [FeII] at 1.26 um ([FeII]$_{1.26}$) luminosity and supernova rate in a sample of 11 near-by starburst galaxy centers. We perform a pixel-pixel analysis of this correlation on SINFONI data cubes. Using Br$γ$ equivalent width and luminosity as the only observational inputs into the Starburst 99 model, we derive the supernova rate at each pixel and thus create maps of supernova rates. We then compare these morphologically and quantitatively to the [FeII]$_{1.26}$ luminosity. We have found that a strong linear and morphological correlation exists between supernova rate and [FeII]$_{1.26}$ on a pixel-pixel basis: \[ log\frac{ν_{SNrate}}{yr^{-1}pc^{-2}} = 1.01 \pm 0.2\ast log\frac{[FeII]_{1.26}}{erg s^{-1}pc^{-2}} - 41.17 \pm 0.9\] This relation is valid for normal star forming galaxies but breaks down for extreme ultra luminous galaxies. The supernova rates derived from the Starburst 99 model are in good agreement with the radio-derived supernova rates, which underlines the strength of using [FeII] emission as a tracer of supernova rate. With the strong correlation found in this sample of galaxies, we conclude that [FeII]$_{1.26}$ emission can be generally used to derive accurate supernova rates on either a pixel-pixel or integrated galactic basis. △ Less

Submitted 13 February, 2012; originally announced February 2012.

[FeII] as a tracer supernova rate

Authors: M. J. F. Rosenberg, P. P. van der Werf, F. P. Israel

Abstract: Supernovae play an integral role in the feedback of processed material into the ISMof galaxies and are responsible for most of the chemical enrichment of the universe. The rate of supernovae can also reveal the star formation histories. Supernova rate is usually measured through the non-thermal radio continuum luminosity, but in this paper we establish a quantitative relationship between the [FeII… ▽ More Supernovae play an integral role in the feedback of processed material into the ISMof galaxies and are responsible for most of the chemical enrichment of the universe. The rate of supernovae can also reveal the star formation histories. Supernova rate is usually measured through the non-thermal radio continuum luminosity, but in this paper we establish a quantitative relationship between the [FeII]1.26 luminosity and supernova rate in a sample of 11 near-by starburst galaxies. SINFONI data cubes are used to perform a pixel pixel analysis of this correlation. Using Br equivalent width and luminosity as the only observational inputs into Starburst 99, the supernova rate is derived at each pixel and a map of supernova rate is created. This is then compared morphologically and quantitatively to [FeII]1.26 luminosity map. We find a strong linear and morphological correlation between supernova rate and [FeII]1.26 on a pixel-pixel basis: log SNrate yr-1 pc-2 = (1.01 \pm 0.2) \ast log[FeII]1.26 ergs-1 pc-2 - 41.17 \pm 0.9 The Starburst 99 derived supernova rates are also in good agreement with the radio derived supernova rates, which further demonstrates the strength of [FeII] as a tracer of supernova rate. With the strong correlation found in this sample of galaxies, we now qualitatively use [FeII]1.26 to derive supernova rate on either a pixel-pixel or integrated galactic basis. △ Less

Submitted 16 January, 2012; originally announced January 2012.

Comments: Conference proceedings for Astronomy at High Angular Resolution 2011: The Central Kiloparsec in Galactic Nuclei

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

Authors: Karin M. Sandstrom, Alberto D. Bolatto, Caroline Bot, B. T. Draine, James G. Ingalls, Frank P. Israel, James M. Jackson, Adam K. Leroy, Aigen Li, Mónica Rubio, Joshua D. Simon, J. D. T. Smith, Snežana Stanimirović, A. G. G. M. Tielens, Jacco Th. van Loon

Abstract: 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 in… ▽ More 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. △ Less

Submitted 5 September, 2011; originally announced September 2011.

Comments: Submitted to ApJ

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

Authors: Alberto D. Bolatto, Adam K. Leroy, Katherine Jameson, Eve Ostriker, Karl Gordon, Brandon Lawton, Snezana Stanimirovic, Frank P. Israel, Suzanne C. Madden, Sacha Hony, Karin M. Sandstrom, Caroline Bot, Monica Rubio, P. Frank Winkler, Julia Roman-Duval, Jacco Th. van Loon, Joana M. Oliveira, Remy Indebetouw

Abstract: 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 effic… ▽ More 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. △ Less

Submitted 8 July, 2011; originally announced July 2011.

Comments: Accepted for publication in the Astrophysical Journal. 21 pages and 11 figures in apjemulate style

C+ Emission from the Magellanic Clouds II. [CII] maps of star-forming regions LMC-N 11, SMC-N 66, and several others

Authors: F. P. Israel, P. R. Maloney

Abstract: We study the 158 micron [CII] fine-structure line emission from star-forming regions as a function of metallicity. We have measured and mapped the [CII] emission from the very bright HII region complexes N 11 in the LMC and N 66 in the SMC, as well as the SMC HII regions N 25, N 27, N 83/N 84, and N 88, with the FIFI instrument on the Kuiper Airborne Observatory. In both the LMC and SMC, the ratio… ▽ More We study the 158 micron [CII] fine-structure line emission from star-forming regions as a function of metallicity. We have measured and mapped the [CII] emission from the very bright HII region complexes N 11 in the LMC and N 66 in the SMC, as well as the SMC HII regions N 25, N 27, N 83/N 84, and N 88, with the FIFI instrument on the Kuiper Airborne Observatory. In both the LMC and SMC, the ratio of the [CII] line to the CO line and to the far-infrared continuum emission is much higher than seen almost anywhere else, including Milky Way star-forming regions and whole galaxies. In the low metallicity, low dust-abundance environment of the LMC and the SMC, UV mean free path lengths are much greater than those in the higher-metallicity Milky Way. The increased photoelectric heating efficiencies cause significantly greater relative [CII] line emission strengths. At the same time, similar decreases in PAH abundances have the opposite effect, by diminishing photoelectric heating rates. Consequently, in low-metallicity environments the relative [CII] strengths are high but exhibit little further dependence on actual metallicity. Relative [CII] strengths are slightly higher in the LMC than in the SMC, which has both lower dust and lower PAH abundances. △ Less

Submitted 11 April, 2011; originally announced April 2011.

Comments: To appear in Astronomy & Astrophysics

The JCMT Nearby Galaxies Legacy Survey VI: The distribution of gas and star formation in M81

Authors: J. R. Sánchez-Gallego, J. H. Knapen, J. S. Heiner, C. D. Wilson, B. E. Warren, R. J. Allen, M. Azimlu, P. Barmby, G. J. Bendo, S. Comerón, F. P. Israel, S. Serjeant, R. P. J. Tilanus, C. Vlahakis, P. van der Werf

Abstract: We present a first complete 12 CO J=3-2 map of M81, observed as part of the Nearby Galaxies Legacy Survey being carried out at the James Clerk Maxwell Telescope. We detect 9 regions of significant CO emission located at different positions within the spiral arms, and confirm that the global CO emission in the galaxy is low. We combine these data with a new H-alpha map obtained using the Isaac Newt… ▽ More We present a first complete 12 CO J=3-2 map of M81, observed as part of the Nearby Galaxies Legacy Survey being carried out at the James Clerk Maxwell Telescope. We detect 9 regions of significant CO emission located at different positions within the spiral arms, and confirm that the global CO emission in the galaxy is low. We combine these data with a new H-alpha map obtained using the Isaac Newton Telescope and archival HI, 24 microns and FUV images to uncover a correlation between the molecular gas and star forming regions in M81. For the nine regions detected in CO J=3-2, we combine our CO J=3-2 data with existing CO J=1-0 data to calculate line ratios. We find that the ratio J=(3-2)/(1-0) is in agreement with the range of typical values found in the literature (0.2-0.8). Making reasonable assumptions, this allows us to constrain the hydrogen density to the range (10^3-10^4) cm^{-3}. We also estimated the amount of hydrogen produced in photo-dissociation regions near the locations where CO J=3-2 was detected. △ Less

Submitted 19 December, 2010; v1 submitted 23 November, 2010; originally announced November 2010.

Comments: Accepted for publication in Astronomy & Astrophysics. 11 pages, 4 figures

The JCMT Nearby Galaxies Legacy Survey IV. Velocity Dispersions in the Molecular Interstellar Medium in Spiral Galaxies

Authors: C. D. Wilson, B. E. Warren, J. Irwin, J. H. Knapen, F. P. Israel, S. Serjeant, D. Attewell, G. J. Bendo, E. Brinks, H. M. Butner, D. L. Clements, J. Leech, H. E. Matthews, S. Muehle, A. M. J. Mortier, T. J. Parkin, G. Petitpas, B. K. Tan, R. P. J. Tilanus, A. Usero, M. Vaccari, P. van der Werf, T. Wiegert, M. Zhu

Abstract: An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 12 nearby spiral galaxies reveals low velocity dispersions in the molecular component of the interstellar medium. The three lowest luminosity galaxies show a relatively flat velocity dispersion as a function of radius while the remaining nine galaxies show a central peak with a radial fall-off within 0.2-0.4 r(25).… ▽ More An analysis of large-area CO J=3-2 maps from the James Clerk Maxwell Telescope for 12 nearby spiral galaxies reveals low velocity dispersions in the molecular component of the interstellar medium. The three lowest luminosity galaxies show a relatively flat velocity dispersion as a function of radius while the remaining nine galaxies show a central peak with a radial fall-off within 0.2-0.4 r(25). Correcting for the average contribution due to the internal velocitydispersions of a population of giant molecular clouds, the average cloud-cloud velocity dispersion across the galactic disks is 6.1 +/- 1.0 km/s (standard deviation 2.9 km/s), in reasonable agreement with previous measurements for the Galaxy andM33. The cloud-cloud velocity dispersion derived from the CO data is on average two times smaller than the HI velocity dispersion measured in the same galaxies. The low cloud-cloudvelocity dispersion implies that the molecular gas is the critical component determining the stability of the galactic disk against gravitational collapse, especially in those regions of the disk which are H2 dominated. The cloud-cloud velocity dispersion shows a significant positivecorrelation with both the far-infrared luminosity, which traces the star formation activity, and the K-band absolute magnitude, which traces the total stellar mass. For three galaxies in the Virgo cluster, smoothing the data to a resolution of 4.5 kpc (to match the typical resolution of high redshift CO observations) increases the measured velocity dispersion by roughly a factor of two, comparable to the dispersion measured recently in a normal galaxy at z=1. This comparison suggests that the mass and star formation rate surface densities may be similar in galaxies from z=0-1 and that the high star formation rates seen at z=1 may be partly due to the presence of physically larger molecular gas disks. △ Less

Submitted 16 September, 2010; originally announced September 2010.

Comments: 16 pages, 23 figures. Accepted to MNRAS